在动静载荷作用下S型波纹管膜片非线性稳定性分析

期刊菜单

在动静载荷作用下S型波纹管膜片非线性稳定性分析
Nonlinear Stability Analysis of S-Type Bellows Diaphragm under Dynamic and Static Loading

DOI:10.12677/IJM.2021.102012,PDF,HTML,XML,下载: 513浏览: 779国家自然科学基金支持
作者:韩明君,李鸿瑞^*,周朝逾,席晓峰,马连生：兰州理工大学理学院，甘肃兰州
关键词:S型波纹管；膜片；Floquet指数；稳定性；S-Shaped Bellows；Diaphragm；Floquet Exponent；Stability

摘要:研究了S型波纹管膜片的非线性稳定性问题，运用拟壳法将S型波纹管膜片当作有初挠度的圆环薄板的复合结构，利用薄壳的非线性弯曲理论，得到S型波纹管膜片在动静载荷作用下的非线性动力学方程组。在边界条件、连续条件下用Galerkin得到膜片非线性系统的受迫振动方程。用Floquet指数判断了该系统发生分叉的条件，讨论了系统在平衡点领域的稳定性问题。

Abstract:he S-shaped bellows diaphragm is regarded as the combined structure of the circular plate with initial deflection by method of simulated shell. Using the nonlinear bending theory of thin shell, the nonlinear dynamic equations of the S-shaped bellows diaphragm under static-dynamic load are obtained. According to the boundary condition and continuous condition, the forced vibration equation of diaphragm nonlinear system is obtained by Galerkin. The conditions for bifurcation of the system were determined by using Floquet index, and the stability at the equilibrium point of the system was discussed.

文章引用：韩明君, 李鸿瑞, 周朝逾, 席晓峰, 马连生. 在动静载荷作用下S型波纹管膜片非线性稳定性分析[J]. 力学研究, 2021, 10(2): 116-126. https://doi.org/10.12677/IJM.2021.102012

1. 引言

波纹管是密封系统中重要的密封元件，有着可以弥补动磨损等因素导致的轴向位移等特点，所以，波纹管干气密封也广泛应用于各种工业领域 [1]。但由于波纹管结构复杂，计算难度较大，刘静等 [2] 通过建立波纹管有限元模型，对比分析了单双层U型波纹管在液压胀形过程中的变形情况。高洪彪 [3] 通过数值模拟的方法对横向载荷作用下波纹管的变形情况进行了分析研究。阿斯耶姆·肖开提等 [4] 利用有限元分析了波纹管在机械密封中温度和应力变化。Ghader等 [5] 采用有限元法对某金属波纹管进行了数值模拟，并最终得到了最佳的设计参数。Krovvidi等 [6] 根据EJMA标准完成了高温下波纹管的初步设计，解决了波纹管的疲劳、面内、柱内失稳等破坏模式。随后，根据RCC-MR给出的准则，使用有限元分析对初步设计进行了分析验证。Hu [7] 在钱伟长研究基础上给出了，并研究了制作方法等原因引起的波纹管的厚度分布变化对波纹管刚度的影响。Chand等 [8] 通过板理论和梁理论对不同的波纹管在轴向载荷下的弯曲变形情况进行了研究分析。李树勋等 [9] 通过运用Ansys建立了U型波纹管模型，研究了波纹管在受到外压载荷和拉压位移载荷的作用下的波纹管的疲劳特性，以及在外载荷作用对波纹管疲劳寿命的影响进行了分析。宋林红等 [10] 通过疲劳模型的建立，运用Smith-Watson Topper法研究了了U型波纹管在不同工况下的疲劳寿命情况。王占彬 [11] 等通过模仿工作中的边界条件，运用数值分析的方法研究了波纹管的稳定性等问题。韩明君 [12] 等通过数值解析的方法运用拟板法及Floquet指数法，对波纹管的稳定性问题进行了研究分析。P. Dinesh Babu [13] 等通过有限元的方法对金属波纹管在其静态下的性能进行了分析研究，同时考虑了温度因素对性能的影像。Aleksander K [14] 等通过得出控制方程，采用Lagrange基本原理以及有限差分法对波纹管的应力应变的问题进行了分析研究。

以上运用多种方法对波纹管的不同问题进行了分析研究，但对S型波纹管膜片的稳定性问题的分析研究的相对比较少。文章通过研究波纹管在动静载荷下的状态，分析了S型波纹管的稳定性，并对在平衡点处的稳定性做出了分析。

2. 基本方程和边界条件

2.1. 基本方程的建立

S型波纹管膜片为圆环轴对称结构，取一半如图1所示，右边虚线为对称轴，a为波纹管膜片的外径，c为内径，b为对称中心到S型弧外侧的距离，d为对称中心到S型弧内侧的距离，h为对称中心到S型弧中间的距离，f为矢高，膜片的内外边界分别自由和固定。在膜片上施加一个均布载荷q，在内径边缘施加集中力P，如图2所示。因波纹管结构的膜片是内外焊接而成，相对变形可简化为一端固定一端自由。

Figure 1. Simplified diagram diaphragm’s model

图1. 膜片模型简化图

Figure 2. Mechanical model of diaphragm

图2. 膜片力学模型图

初挠度为：

$\begin{array}{l} w_{01} = - f \frac{{[r - (d + 3 R)]}^{2}}{R^{2}}, \frac{d w_{01}}{d r} = - f \frac{2 [r - (d + 3 R)]}{R^{2}} \\ w_{02} = f \frac{{[r - (d + R)]}^{2}}{R^{2}}, \frac{d w_{02}}{d r} = f \frac{2 [r - (d + R)]}{R^{2}} \end{array}$

S型波纹管膜片在静态载荷下的基本方程组：

$D \frac{1}{r} \frac{d}{d r} r \frac{d}{d r} \frac{1}{r} \frac{d}{d r} r \frac{d w_{d}}{d r} = q_{j} + \frac{1}{r} \frac{d}{d r} [r N_{r j} (\frac{d w_{j}}{d r} + λ_{1} \frac{d w_{0}}{d r} + λ_{2} \frac{d w_{0}}{d r})]$ (1)

$\frac{1}{E h} r \frac{d}{d r} \frac{1}{r} \frac{d}{d r} (r^{2} N_{r j}) = - \frac{1}{2} (\frac{d w_{j}}{d r} + λ_{1} \frac{d w_{0}}{d r} + λ_{2} \frac{d w_{0}}{d r}) \frac{d w_{d}}{d r}$ (2)

将S型波纹管的膜片在动静载荷作用下发生的非线性大挠度看作波纹管膜片具有的初挠度，并将惯性项与阻尼项计算其中，计算出在动静载荷下轴对称非线性基本方程。

$\begin{array}{l} D \frac{1}{r} \frac{d}{d r} r \frac{d}{d r} \frac{1}{r} \frac{d}{d r} r \frac{d w_{d}}{d r} \\ = q_{d} + \frac{1}{r} \frac{d}{d r} [r N_{r j} \frac{d w_{d}}{d r} + r N_{r d} (\frac{d w_{j}}{d r} + \frac{d w_{d}}{d r} + λ_{1} \frac{d w_{0}}{d r} + λ_{2} \frac{d w_{0}}{d r})] - c^{'} \frac{d w_{d}}{d t} - γ \frac{d^{2} w_{d}}{d t^{2}} \end{array}$ (3)

$\frac{1}{E h} r \frac{d}{d r} \frac{1}{r} \frac{d}{d r} (r^{2} N_{r d}) = - \frac{1}{2} (\frac{d w_{d}}{d r} + \frac{d w_{j}}{d r} + λ_{1} \frac{d w_{0}}{d r} + λ_{2} \frac{d w_{0}}{d r}) \frac{d w_{d}}{d r}$ (4)

“ $λ_{1}$ ”和“ $λ_{2}$ ”为引入的两个关于a的阶梯函数，分别为：

$λ_{1}, λ_{2}$ 取值：

$λ_{1} = {\begin{cases} 1, \frac{c}{a} \leq x \leq \frac{d}{a} \\ 0, 其他 \end{cases}$ ， $λ_{2} = {\begin{cases} - 1, \frac{h}{a} \leq x \leq \frac{b}{a} \\ 1, \frac{d}{a} \leq x \leq \frac{h}{a} \\ 0, 其他 \end{cases}$

E代表杨氏模量， $D = \frac{E h^{3}}{12 (1 - μ^{2})}$ ， $c^{'}$ 代表S型波纹管膜片的阻尼系数， $w_{0}$ 代表初挠度， $γ$ 代表S

型波纹管膜片的单位面积的体密度， $μ$ 代表泊松比。

对方程(3) (4)进行哈密尔顿原理和变分可得，非线性动力学变分方程和协调方程。

$\begin{array}{l} \int_{t_{1}}^{t_{2}} \int_{c}^{a} {D \frac{1}{r} \frac{\partial}{\partial r} r \frac{\partial}{\partial r} \frac{1}{r} \frac{\partial}{\partial r} r \frac{\partial}{\partial r} (w_{d}) - q_{d} - \frac{1}{r} \frac{\partial}{\partial r} [r N_{r j} \frac{\partial w_{d}}{\partial r} + r N_{r d} (\frac{\partial w_{j}}{\partial r} + \frac{\partial w_{d}}{\partial r} + λ_{1} \frac{\partial w_{0}}{\partial r} + λ_{2} \frac{\partial w_{0}}{\partial r})] \\ - \frac{P}{2 π} + c \frac{\partial w_{d}}{\partial t} + γ \frac{\partial^{2} w_{d}}{\partial t^{2}}} r δ w_{d} d r d t = 0 \end{array}$ (5)

$\frac{1}{E h} r \frac{\partial}{\partial r} \frac{1}{r} \frac{\partial}{\partial r} (r^{2} N_{r d}) = - [\frac{1}{2} {(\frac{\partial w_{d}}{\partial r})}^{2} + \frac{\partial w_{j}}{\partial r} \frac{\partial w_{d}}{\partial r} + λ_{1} \frac{\partial w_{d}}{\partial r} \frac{\partial w_{0}}{\partial r} + λ_{2} \frac{\partial w_{d}}{\partial r} \frac{\partial w_{0}}{\partial r}]$ (6)

边界条件为：

$r = a, w_{d (1)} = 0, \frac{d w_{d (1)}}{d r} = 0, \frac{d}{d r} (r N_{r (1)}) - μ N_{r (1)} = 0$ (7)

自由周边处：

$r = c, \frac{d^{2} w_{d (4)}}{d r^{2}} + \frac{μ}{r} \frac{d w_{d (4)}}{d r} = 0, N_{r (4)} = 0$ (8)

连续处：

$r = b, w_{d (1)} = w_{d (2)}, \frac{d w_{d (1)}}{d r} = \frac{d w_{d (2)}}{d r}, M_{r (1)} = M_{r (2)}, N_{r (1)} = N_{r (2)}, N_{θ (1)} = N_{θ (2)}$ (9)

$r = h, w_{d (2)} = w_{d (3)}, \frac{d w_{d (2)}}{d r} = \frac{d w_{d (3)}}{d r}, M_{r (2)} = M_{r (3)}, N_{r (2)} = N_{r (3)}, N_{θ (2)} = N_{θ (3)}$ (10)

$r = d, w_{d (3)} = w_{d (4)}, \frac{d w_{d (3)}}{d r} = \frac{d w_{d (4)}}{d r}, M_{r (3)} = M_{r (4)}, N_{r (3)} = N_{r (4)}, N_{θ (3)} = N_{θ (4)}$ (11)

2.2. 基本方程和边界条件的无量纲化

$\begin{array}{l} x = \frac{r}{a}, W = \frac{w}{h}, S = \frac{12 a (1 - μ^{2}) r N_{r} a}{E h^{3}}, β Q = \frac{12 a^{4} (1 - μ^{2}) q}{E h^{4}}, k = f \frac{2 a^{2}}{R^{2} h}, ϕ_{1} = f \frac{2 a (d + 3 R)}{R^{2} h}, \\ ϕ_{2} = - f \frac{2 a (d + R)}{R^{2} h}, c^{″} = \frac{12 a^{4} (1 - μ^{2})}{E h^{3}} c^{'}, γ_{1} = \frac{12 a^{4} (1 - μ^{2})}{E h^{3}} γ, α = 6 (1 - μ^{2}) \end{array}$

把以上量纲带入(5) (6)，计算求出无量纲变分方程(12)：

$\begin{array}{l} \int_{0}^{2 π} \int_{\frac{c}{a}}^{1} {L_{1} (W_{d}) - Q - \frac{1}{x} \frac{\partial}{\partial x} [S_{j} \frac{\partial W_{d}}{\partial x} + S_{d} (\frac{\partial W_{j}}{\partial x} + \frac{\partial W_{d}}{\partial x} + λ_{1} (ϕ_{1} - k x) + λ_{2} (k x - ϕ_{2}))] \\ + c^{″} \frac{\partial W_{d}}{\partial τ} + γ_{1} \frac{\partial^{2} W_{d}}{\partial τ^{2}}} x W_{d} d x d τ = 0 \end{array}$ (12)

无量纲变形协调方程：

$L_{2} (x S_{d}) = - [α {(\frac{\partial W_{d}}{\partial x})}^{2} + 2 α \frac{\partial W_{j}}{\partial x} \frac{\partial W_{d}}{\partial x} + 2 α \frac{\partial W_{d}}{\partial x} (λ_{1} (ϕ_{1} - k x) + λ_{2} (k x - ϕ_{2}))]$ (13)

其中 $L_{1} = \frac{1}{x} \frac{\partial}{\partial x} x \frac{\partial}{\partial x} \frac{1}{x} \frac{\partial}{\partial x} x \frac{\partial}{\partial x}$ ， $L_{2} = x \frac{\partial}{\partial x} \frac{1}{x} \frac{\partial}{\partial x}$ 。

由无量纲表示的边界条件为：

固定周边处：

$x = 1, W_{d(1)} = 0, \frac{d W_{d (1)}}{d x} = 0, \frac{d S_{d (1)}}{d x} - μ S_{d (1)} = 0$ (14)

自由周边处：

$x = \frac{c}{a}, \frac{d^{2} W_{d (4)}}{d x^{2}} + \frac{μ}{x} \frac{d W_{d (4)}}{d x} = 0, S_{d (4)} = 0$ (15)

连续处：

$\begin{array}{l} x = \frac{b}{a}, W_{d (1)} = W_{d (2)}, \frac{d^{2} W_{d (1)}}{d x^{2}} + \frac{μ}{x} \frac{d W_{d (1)}}{d x} = \frac{d^{2} W_{d (2)}}{d x^{2}} + \frac{μ}{x} \frac{d W_{d (2)}}{d x}, \\ \frac{d W_{d (1)}}{d x} = \frac{d W_{d (2)}}{d x}, S_{d (1)} = S_{d (2)}, \frac{d}{d x} (x S_{d (1)}) = \frac{d}{d x} (x S_{d (2)}) \end{array}$ (16)

$\begin{array}{l} x = \frac{h}{a}, W_{2 (2)} = W_{d (3)}, \frac{d^{2} W_{d (2)}}{d x^{2}} + \frac{μ}{x} \frac{d W_{d (2)}}{d x} = \frac{d^{2} W_{d (3)}}{d x^{2}} + \frac{μ}{x} \frac{d W_{d (3)}}{d x}, \\ \frac{d W_{d (2)}}{d x} = \frac{d W_{d (2)}}{d x}, S_{d (2)} = S_{d (3)}, \frac{d}{d x} (x S_{d (2)}) = \frac{d}{d x} (x S_{d (3)}) \end{array}$ (17)

$\begin{array}{l} x = \frac{d}{a}, W_{d (3)} = W_{d (4)}, \frac{d^{2} W_{d (3)}}{d x^{2}} + \frac{μ}{x} \frac{d W_{d (3)}}{d x} = \frac{d^{2} W_{d (4)}}{d x^{2}} + \frac{μ}{x} \frac{d W_{d (4)}}{d x}, \\ \frac{d W_{d (3)}}{d x} = \frac{d W_{d (4)}}{d x}, S_{d (3)} = S_{d (4)}, \frac{d}{d x} (x S_{d (3)}) = \frac{d}{d x} (x S_{d (4)}) \end{array}$ (18)

3. 受迫振动方程

由S型波纹管的非线性弯曲可知W_j和S_j。 $W_{j (1)}, W_{j (2)}, W_{j (3)}, W_{j (4)}$ 比较繁长，这里不再进行展开。

$\begin{matrix} S_{j (1)} = \frac{1}{x} (- \frac{1}{12288} a Q_{1}^{2} x^{8} - \frac{1}{384} a Q_{1} c_{14} x^{6} - \frac{1}{64} a x^{4} Q_{1} c_{15} - \frac{1}{32} a x^{4} c_{14}^{2} - \frac{1}{2} a c_{14} c_{15} \ln (x) x^{2} \\ + \frac{1}{4} a c_{14} c_{15} x^{2} + \frac{1}{2} a c_{15}^{2} \ln (x)) + \frac{1}{2} c_{16} x + \frac{c_{17}}{x} \end{matrix}$ (19)

$\begin{matrix} S_{j (2)} = \frac{1}{x} (- \frac{1}{12288} a Q_{1}^{2} x^{8} + \frac{1}{192} a x^{6} k Q_{1} - \frac{1}{384} a Q_{1} n_{14} x^{6} - \frac{1}{120} a Q_{1} ϕ_{1} x^{5} - \frac{1}{64} a x^{4} Q_{1} n_{15} \\ + \frac{1}{8} a x^{4} k n_{14} - \frac{1}{32} a x^{4} n_{14}^{2} + a n_{15} \ln (x) k x^{2} - \frac{1}{2} a n_{15} k x^{2} - \frac{1}{2} a n_{14} n_{15} \ln (x) x^{2} \\ + \frac{1}{4} a n_{14} n_{15} x^{2} - \frac{1}{3} a n_{14} ϕ_{1} x^{3} + 2 a n_{15} ϕ_{1} x + \frac{1}{2} a n_{15}^{2} \ln (x)) + \frac{1}{2} n_{16} x + \frac{n_{17}}{x} \end{matrix}$ (20)

$\begin{matrix} S_{j (3)} = \frac{1}{x} (- \frac{1}{12288} a Q_{1}^{2} x^{8} - \frac{1}{192} a x^{6} k Q_{1} - \frac{1}{384} a Q_{1} m_{14} x^{6} + \frac{1}{120} a Q_{1} ϕ_{1} x^{5} - \frac{1}{64} a x^{4} Q_{1} m_{15} \\ - \frac{1}{8} a x^{4} k m_{14} - \frac{1}{32} a x^{4} m_{14}^{2} - a m_{15} \ln (x) k x^{2} + \frac{1}{2} a m_{15} k x^{2} - \frac{1}{2} a m_{14} m_{15} \ln (x) x^{2} \\ + \frac{1}{4} a m_{14} m_{15} x^{2} + \frac{1}{3} a m_{14} ϕ_{2} x^{3} - 2 a m_{15} ϕ_{1} x + \frac{1}{2} a m_{15}^{2} \ln (x)) + \frac{1}{2} m_{16} x + \frac{m_{17}}{x} \end{matrix}$ (21)

$\begin{matrix} S_{j (4)} = \frac{1}{x} (- \frac{1}{12288} a Q_{1}^{2} x^{8} - \frac{1}{384} a Q_{1} h_{14} x^{6} - \frac{1}{64} a x^{4} Q_{1} h_{15} - \frac{1}{32} a x^{4} h_{14}^{2} - \frac{1}{2} a h_{14} h_{15} \ln (x) x^{2} \\ + \frac{1}{4} a h_{14} h_{15} x^{2} + \frac{1}{2} a h_{15}^{2} \ln (x)) + \frac{1}{2} h_{16} x + \frac{h_{17}}{x} \end{matrix}$ (22)

通过边界条件可求得S型波纹管膜片在动载荷下的挠度：

$W_{d} = (\frac{1}{2} x^{4} + h_{21} x^{2} + h_{22} \ln x + h_{23}) f (t)$ (23)

将W_j，W_d带入方程(13)得到S_d，再将W_j，W_d，S_d，S_j带入方程(12)通过Glalerkin法并化简后可得：

$\frac{d^{2} f (t)}{d t^{2}} + φ_{0} \frac{d f (t)}{d t} + φ_{1} f (t) + φ_{2} f^{2} (t) + φ_{3} f^{3} (t) = g \cos (Ω t)$ (24)

令 $φ_{1} = ω^{2}$ ，取 $t = ω τ$ ， $f (t) = \frac{ω}{\sqrt{φ_{3}}} η (τ)$ 。

由(24)可得到S型波纹管膜片线性系统的受迫振动方程：

$\frac{d^{2} η (τ)}{d τ^{2}} + α^{'} \frac{d η (τ)}{d τ} + η (τ) - χ^{'} η^{2} (τ) + η^{3} (τ) = g^{'} \cos (\frac{Ω}{ω} τ)$ (25)

其中 $α^{'} = \frac{φ_{0}}{ω}, χ^{'} = - \frac{φ_{2}}{ω \sqrt{φ_{3}}}, g^{'} = \frac{g \sqrt{φ_{3}}}{ω^{2}}$ 。

4. 对自由振动方程求解

方程(25)的等价方程为：

${\begin{cases} {η^{'}}_{1} (τ) = η_{2} (τ) \\ {η^{'}}_{2} (τ) = - α^{'} η_{1} (τ) - η^{3} (τ) + χ^{'} η_{1}^{2} (τ) - η_{1}^{3} (τ) + g^{'} \cos (\frac{Ω}{ω} τ) \end{cases}$ (26)

令方程(26) $α^{'} = 0$ 和 $g^{'} = 0$ ，则方程为：

${\begin{cases} {η^{'}}_{1} (τ) = η_{2} (τ) \\ {η^{'}}_{2} (τ) = - η^{3} (τ) + χ^{'} η_{1}^{2} (τ) - η_{1}^{3} (τ) \end{cases}$ (27)

得该系统的Hamilton初积分为：

$H (η_{1} (τ), η_{2} (τ)) = \frac{{({η^{'}}_{1} (τ))}^{2}}{2} + \frac{η_{1}^{2} (τ)}{2} - \frac{χ^{'}}{3} η_{1}^{3} (τ) + \frac{1}{4} η_{1}^{4} (τ)$ (28)

另H = 0，则(28)可以化简为：

$η_{2}^{2} = {({η^{'}}_{1})}^{2} = - η_{1}^{2} + \frac{2}{3} χ^{'} η_{1}^{3} - \frac{1}{2} η_{1}^{4}$ (29)

对(29)求解可得：

$\frac{d η_{1} (τ)}{d τ} = \pm η_{1} \sqrt{- \frac{η_{1}^{2}}{2} + \frac{2}{3} χ^{'} η_{1} - 1}$ ，则 $τ = \pm (\arcsin \frac{\frac{2}{3} χ^{'} η_{1} - 2}{| η_{1} | \sqrt{\frac{2}{3} χ^{'}} - 2})$

令 $u = \frac{2}{3} χ^{'}$ ， $v = \sqrt{{(\frac{2}{3} χ^{'})}^{2} - 2} = \sqrt{u^{2} - 2}$ ，则 $τ = \pm (\arcsin \frac{u η_{1} - 2}{v η_{1}} + C)$ 。

取C = 0，可以得到：

$η_{1} (τ) = \frac{2}{u \pm v \sin τ}$ ， $η_{2} (τ) = \frac{\mp 2 v \cos τ}{{(u \pm v \sin τ)}^{2}}$

经验证：

$(η_{1} (τ), η_{2} (τ)) = (\frac{2}{u + v \sin τ}, \frac{- 2 v \cos τ}{{(u + v \sin τ)}^{2}})$

$(η_{1} (τ), η_{2} (τ)) = (\frac{2}{u - v \sin τ}, \frac{2 v \cos τ}{{(u - v \sin τ)}^{2}})$

5. 研究平衡点处的稳定性

在不存在外激励情况下( $g^{'} = 0$ )，可将方程(25)看作是：

${\begin{cases} {η^{'}}_{1} (τ) = η_{2} (τ) \\ {η^{'}}_{2} (τ) = - α^{'} η_{2} (τ) - η_{1} (τ) + χ^{'} η_{1}^{2} (τ) - η_{1}^{3} (τ) \end{cases}$ (30)

令 ${\begin{cases} f_{1} (η_{1}, η_{2}) = η_{2} (τ) \\ f_{2} (η_{1}, η_{2}) = - α^{'} η_{2} (τ) - η_{1} (τ) + χ^{'} η_{1}^{2} (τ) - η_{1}^{3} (τ) \end{cases}$ ；

Jacobi矩阵是：

$[\begin{matrix} \frac{\partial f_{1}}{\partial η_{1}} & \frac{\partial f_{1}}{\partial η_{2}} \\ \frac{\partial f_{2}}{\partial η_{1}} & \frac{\partial f_{2}}{\partial η_{2}} \end{matrix}] = [\begin{matrix} 0 & 1 \\ - 1 + 2 χ^{'} η_{1} - 3 η_{1}^{2} & - α^{'} \end{matrix}]$

显然在 $2 \leq χ^{'}$ 有三平衡点， $(0, 0), (\frac{α^{″}}{2} \pm \sqrt{{(\frac{χ^{'}}{2})}^{2} - 1}, 0)$ 。

以下依据Floquet指数来分析平衡点处的稳定性

1) 在平衡点是 $(0, 0)$ 情况，Jacobi矩阵为 $[\begin{matrix} 0 & 1 \\ - 1 & - α^{'} \end{matrix}]$ 。

特征方程是 $| \begin{array}{l} 0 - δ 1 \\ - 1 - α^{'} - δ \end{array} | = 0$ ，解得 $δ_{1, 2} = \frac{- α^{'} \pm \sqrt{{(α^{'})}^{2} - 4}}{2}$ 。

在 $0 < α^{'} < 2$ 的情况下， $δ_{1, 2} = \frac{- α^{'}}{2} \pm i \sqrt{1 - {(\frac{α^{'}}{2})}^{2}}$ 是不等的复数，复平面上有稳定交点，图3。

$α^{'} = 0.5$ $χ^{'} = 2$

Figure 3. Phase diagram ( $0 < α^{'} < 2$ )

图3. 相图( $0 < α^{'} < 2$ )

在 $α^{'} = 2$ 的情况下， $δ_{1, 2} = - 1$ 是相等负数，平衡点是临界结点，图4。

$α^{'} = 2$ $χ^{'} = 2$

Figure 4. Phase diagram ( $α^{'} = 2$ )

图4. 相图( $α^{'} = 2$ )

在 $α^{'} > 2$ 的情况下， $δ_{1, 2} = - \frac{α^{'}}{2} \pm \sqrt{{(\frac{α^{'}}{2})}^{2} - 1}$ 是不等负数，平衡点是稳定结点，图5。

$α^{'} = 4$ $χ^{'} = 2$

Figure 5. Phase diagram ( $α^{'} > 2$ )

图5. 相图( $α^{'} > 2$ )

在 $α^{'} = 0$ 的情况下， $δ_{1, 2} = \pm i$ 是纯虚数，求得曲线是极限环，图6。

$χ^{'} = 2$

Figure 6. Phase diagram ( $α^{'} = 0$ )

图6. 相图( $α^{'} = 0$ )

2) 在平衡点 $(\frac{χ^{'}}{2} + \sqrt{{(\frac{χ^{'}}{2})}^{2} - 1}, 0)$ 处

Jacobi矩阵是 $[\begin{matrix} 0 & 1 \\ \frac{4 - {(χ^{'})}^{2} - χ^{'} \sqrt{{(χ^{'})}^{2} - 4}}{2} & - α^{'} \end{matrix}]$ ，特征方程是 $| \begin{array}{l} 0 - δ 1 \\ 2 - \frac{{(χ^{'})}^{2}}{2} - \frac{χ^{'}}{2} \sqrt{{(χ^{'})}^{2} - 4} - α^{'} - δ \end{array} | = 0$ ，解得 $δ_{1, 2} = \frac{- α^{'} \pm \sqrt{{(α^{'})}^{2} + (8 - 2 {(χ^{'})}^{2} - 2 χ^{'} \sqrt{{(χ^{'})}^{2} - 4})}}{2}$ 。

在 $0 < α^{'} < \sqrt{- 8 + 2 {(χ^{'})}^{2} - 2 χ^{'} \sqrt{{(χ^{'})}^{2} - 4}}$ 情况下， $δ_{1, 2} = \frac{- α^{'}}{2} \pm \frac{i}{2} \sqrt{- {(α^{'})}^{2} - 8 + 2 {(χ^{'})}^{2} + 2 χ^{'} \sqrt{{(χ^{'})}^{2} - 4}}$ 是不等复数，复平面有稳定交点；

在 $α^{'} = \sqrt{- 8 + 2 {(χ^{'})}^{2} + 2 χ^{'} \sqrt{{(χ^{'})}^{2} - 4}}$ 情况下， $δ_{1, 2} = - \frac{α^{'}}{2}$ 是相等负数， $(\frac{χ^{'}}{2} + \sqrt{{(\frac{χ^{'}}{2})}^{2} - 1}, 0)$ 是临界结点；

在 $α^{'} > \sqrt{- 8 + 2 {(χ^{'})}^{2} + 2 χ^{'} \sqrt{{(χ^{'})}^{2} - 4}}$ 情况下， $δ_{1, 2} = \frac{- α^{'} \pm \sqrt{{(α^{'})}^{2} - (8 - 2 {(χ^{'})}^{2} + 2 χ^{'} \sqrt{{(χ^{'})}^{2} - 4})}}{2}$ 是不等负数， $(\frac{χ^{'}}{2} + \sqrt{{(\frac{χ^{'}}{2})}^{2} - 1}, 0)$ 是稳定结点；

在 $α^{'} = 0$ 情况下，特征根 $δ_{1, 2} = \pm i \sqrt{- 8 + 2 {(χ^{'})}^{2} + 2 χ^{'} \sqrt{{(χ^{'})}^{2} - 4}}$ 是纯虚数，有Hopf分岔。

3) 在平衡点 $(\frac{χ^{'}}{2} - \sqrt{{(\frac{χ^{'}}{2})}^{2} - 1}, 0)$ 处

Jacobi矩阵 $[\begin{matrix} 0 & 1 \\ \frac{4 - {(χ^{'})}^{2} + χ^{'} \sqrt{{(χ^{'})}^{2} - 4}}{2} & - α^{'} \end{matrix}]$ ， $(\frac{χ^{'}}{2} - \sqrt{{(\frac{χ^{'}}{2})}^{2} - 1}, 0)$ 处的特征根是 $δ_{1, 2} = \frac{- α^{'} \pm \sqrt{{(α^{'})}^{2} - 8 + 2 {(χ^{'})}^{2} - 2 χ^{'} \sqrt{{(χ^{'})}^{2} - 4}}}{2}$ 。 $δ_{1, 2}$ 是不等负数，平衡点是稳定结点。

6. 小结

1) S型波纹管非线性力学自由振动方程，准确解

$(η_{1} (τ), η_{2} (τ)) = (\frac{2}{u + v \sin τ}, \frac{- 2 v \cos τ}{{(u + v \sin τ)}^{2}})$

$(η_{1} (τ), η_{2} (τ)) = (\frac{2}{u - v \sin τ}, \frac{2 v \cos τ}{{(u - v \sin τ)}^{2}})$

2) S型波纹管在动静载荷下，平衡点 $(0, 0)$ 处，在 $0 < α^{'} < 2$ 的情况下，根是不等复数，有稳定的交点；在 $α^{'} = 2$ 的情况下，根是相等负数，平衡点是临界点；在 $α^{'} > 2$ 的情况，根为不等负数，平衡点为稳定结点；在 $α^{'} = 0$ 的情况下，发生Hopf分岔。

基金项目

国家自然科学基金项目(11472123)，国家自然科学基金项目(11862012)。

NOTES

^*通讯作者。

参考文献

[1]	张静, 丁雪兴, 黄义仿, 等. 干气密封焊接金属波纹管的结构动态分析[J]. 化工机械, 2010, 37(2): 183-185.
[2]	刘静, 黄卫良, 李霄, 等. 单双层U形不锈钢波纹管液压胀形变形特征对比研究[J]. 锻压技术, 2016, 41(5): 24-28.
[3]	高洪彪. 横向载荷下弯曲波纹管的变形与吸能特性[D]: [硕士学位论文]. 太原: 太原理工大学, 2016.
[4]	阿斯耶姆•肖开提, 刘焕海. 焊接金属波纹管机械密封温度场及热变形有限元分析研究[J]. 科技视界, 2016(4): 25-26.
[5]	Faraji, G., Besharati, M.K., Mosavi, M., et al. (2008) Experimental and Finite Element Analysis of Parameters in Manufacturing of Metal Bellows. The International Journal of Advanced Manufacturing Technology, 38, 641-648. https://doi.org/10.1007/s00170-007-1122-9
[6]	Kumar Krovvidi, S.C.S.P., Padmakumar, G. and Bhaduri, A.K. (2016) Design and Analysis of Bellows for High Temperature Application as Per RCC-MR. Transactions of the Indian Institute of Metals, 69, 531-535. https://doi.org/10.1007/s12666-015-0813-4
[7]	Liang, H. (1993) On Problems of U-Shaped Bellows with Non-linear Deformation of Large Axisymmetrical Deflection (II)—Counting Variation of Thickness Distribution. Applied Mathematics and Mechanics, 14, 559-564. https://doi.org/10.1007/BF02451365
[8]	Chand, S. and Garg, S.B.L. (1981) Design Formulas for Expansion Bellows. Journal of Mechanical Design, 103, 881-891. https://doi.org/10.1115/1.3255001
[9]	李树勋, 张丽芳, 吕兴, 等. 二级波纹管截止阀用波纹管结构强度及疲劳可靠性研究[J]. 流体机械, 2018, 46(11): 39-44.
[10]	宋林红, 黄乃宁, 等. 金属波纹管疲劳寿命的有限元分析[J]. 管道技术与设备, 2008(11): 16-18.
[11]	王占彬, 周浩洋, 石朝锋, 等. 一种U型波纹管强度及稳定性分析方法研究[J]. 强度与环境, 2012, 39(4): 21-26.
[12]	韩明君, 王伟兵, 李鸿瑞, 等. 单圆弧波纹管膜片的非线性大变形分析[J]. 工程力学, 2020, 37(1): 26-33.
[13]	Dinesh Babu, P., Keerthi Prasath, S., Barani Dharan, M., Venkat Raman, C., Narayanan, R. and Ganesh, K.C. (2014) Analysis of Static Mechanical Behaviour of Metal Bellows Using Finite Element Modeling. Applied Mechanics and Materials, 592-594, 996-1000. https://doi.org/10.4028/www.scientific.net/AMM.592-594.996
[14]	Belyaev, A.K., Zinovieva, T.V. and Smirnov, K.K. (2017) Theoretical and Experimental Studies of the Stress-Strain State of Expansion Bellows as Elastic Shells. St. Petersburg Polytechnical University Journal: Physics and Mathematics, 3, 7-14. https://doi.org/10.1016/j.spjpm.2017.03.003

为你推荐

友情链接