基于波動法的靜水壓力下環肋圓柱殼耦合振動特性研究

梁 斌，李 戎，劉小宛，NODA Nao－Aki，徐紅玉

(1.河南科技大學 土木工程學院，河南 洛陽 471023;2.九州工業大學 機械系，日本 北九州 804－8550)

由于肋骨加強的圓柱殼體結構能夠在減小自身重量的前提下達到增強結構剛度和強度的作用，加肋圓柱殼已被廣泛應用至潛艇耐壓船體結構、航空航天等領域，研究環肋圓柱殼在靜水壓力下的振動特性對于研究水下圓柱殼體結構的力學行為有著非常重要的意義。Junger［1－2］最早對浸沒在聲場中的各向同性圓柱殼進行研究。Zhang等［3－4］利用波動法研究了考慮流體影響的圓柱殼的振動特性。Loy等［5］利用Rayleigh－Ritz法研究了僅含有一條環肋時圓柱殼的振動特性。Liu等［6］考慮了靜水壓力對水下環肋圓柱殼輸入功率的影響。Gan等［7］以環肋圓柱殼為研究對象討論了波動法在對此類問題進行分析時的計算精度問題。Tatianna等［8］采用 Flügge 殼體理論研究了側向壓力下圓柱殼體的穩定性。Zhu等［9］基于波動法討論了通過殼體的固有頻率和靜水壓力來得到不加肋水下圓柱殼臨界壓力的方法的可行性。

本文在Flügge理論的基礎上，考慮流體的影響，采用波動法得到耦合系統的振動方程。再通過變換軸向波數，利用牛頓迭代法得到連續變化的靜水壓力下、不同邊界條件時水下圓柱殼的固有頻率。通過分別退化計算不考慮液體影響的環肋圓柱殼和水下不加肋圓柱殼的固有頻率，驗證了本文方法的有效性和正確性。作為算例，計算并分析了肋條尺寸和數目、靜水壓力、邊界條件等因素對殼體固有頻率的影響。

1 力學模型

假定環肋圓柱殼(外加肋)如圖1所示，半徑為R，長度為L，壁厚為h，肋條截面長度為hr，寬度為 br，肋條間距為d2。將正交坐標系(x，θ，z)建立在殼體中面上，其中x，θ和z分別為殼體的軸向、環向和徑向坐標。

圖1 環肋圓柱殼模型Fig.1 Geometry of a ring-stiffened cylindrical shell

2 理論推導

由Flügge［10］理論可以得到考慮靜水壓力時圓柱殼的運動方程

圓柱殼的振動位移函數可以表示為

式中，Um，Vm，Wm分別表示 x，θ，z方向的波幅，ω 為固有角頻率。

考慮液體影響的圓柱殼滿足聲波方程，在柱坐標系(x，θ，r)中，流體的運動方程為

式中，t為時間，c為流體的聲速，ψ為聲壓，r坐標沿殼體徑向選取。而在流體中，滿足聲波方程(3)的聲壓場在柱坐標系中的形式，可以表示為

式(5)中，Ω，CL和CF分別為無量綱頻率參數，圓柱殼體的聲速及流體的聲速。由于在殼體外壁與流體的接觸界面上，流體徑向位移與殼體徑向位移必須相等。該耦合條件表達式為

針對不同情況，代入具體的邊界條件至km［11］，求解式(8)，可以得到

式中，P1(ω)和P2(ω)分別為關于未知數ω的多項式。在一定外壓下，采用牛頓迭代法即可得到給定靜水壓力時環肋圓柱殼的固有頻率值。當FL=0，Q=0時，即可還原為不考慮流體影響情況下環肋圓柱殼的固有頻率的計算。

3 算例與討論

表1計算了靜水壓力下，不考慮肋條影響時三種邊界條件下圓柱殼的固有頻率。通過與Zhu［9］的結果進行對比，驗證了本文在計算靜水壓力下殼體固有頻率時的有效性和正確性。殼體的幾何參數為:E=2.1×1011N/m2，μ =0.3，ρ=7 850 kg/m3，L/R=20，h/R=0.01，R=1 m，m=1，n=2。

表1 不同邊界條件下水下圓柱殼固有頻率的對比分析Tab.1 Comparison offundamental natural frequencies under different hydrostatic pressures in various boundary conditions

表2計算了兩端簡支邊界條件下，不考慮液體影響時環肋圓柱殼的固有頻率。通過與參考文獻［7，12］的結果進行對比，驗證了本文在計算環肋圓柱殼的固有頻率時的有效性和正確性。

由于當靜水壓力超過臨界壓力時，結構會失穩，再進行固有頻率的計算已無意義，所以對圓柱殼容器進行結構分析時，需要考慮其臨界壓力大小。而對于靜水壓力下的環肋圓柱殼，當殼體的固有頻率為0時，則可認為此時所對應的靜水壓力即為環肋圓柱殼臨界荷載的彈性理論解。因此，本文在進行水下環肋圓柱殼的振動特性分析時，靜水壓力數值均控制在其殼體分別對應的臨界壓力以下。作為算例，本文分別計算了連續變化的靜水壓力下，變換h/R、軸向半波數、肋條尺寸和數目、邊界條件等多種影響因素時圓柱殼的固有頻率變化曲線，見圖2至圖10。殼體的幾何參數為:E=2.077 88 × 1011N/m2，μ =0.317 756，ρ=8 166 kg/m3，L/R=20，R=1 m，m=1，hr=0.005 m，br=0.003 m，肋條數 Nr=20。

表2 環肋圓柱殼固有頻率的對比分析Tab.2 Comparison of natural frequencies of stiffened cylindrical shell

為了更好的表現肋條對于圓柱殼的固有頻率的影響，下面首先以兩端簡支邊界條件下的環肋圓柱殼作為研究對象，分析不考慮靜水壓力時肋條對殼體固有頻率的影響，h/R=0.002。圖2分別以不加肋和環肋圓柱殼作為研究對象，給出了不同軸向半波數時，兩端簡支邊界條件下圓柱殼的固有頻率的變化規律曲線。從圖中可以看出，加肋圓柱殼的固有頻率明顯大于不加肋圓柱殼的固有頻率，且隨著周向波數的逐漸增大，肋條對固有頻率的影響逐漸顯現;與不加肋圓柱殼相比，肋條對殼體的固有頻率的影響主要表現在周向波數較大的情況下;軸向半波數對固有頻率的影響主要表現在周向波數較小的情況下。

圖3給出了不同肋條數目(當殼體長度一定時，增減肋條數目等同于改變肋條間距)時，殼體固有頻率的變化曲線。可以看出，隨著肋條數目的不斷增多(肋條間距的逐漸減小)，殼體的固有頻率逐漸增大，肋條數目對固有頻率的影響主要表現在周向波數較大的情況下;隨著周向波數的逐漸增大，固有頻率先減小后增大。

圖4以肋條截面長度和寬度比作為出發點，研究了肋條尺寸對殼體固有頻率的影響，br=0.003 m。由圖可知，殼體的固有頻率隨著肋條截面長度和寬度比的增大而增大，且其比例值越大，頻率增長幅度約明顯;長度和寬度比對固有頻率的影響主要表現在周向波數較大的情況下;與肋條數目(肋條間距)相比，肋條的截面長度和寬度比對固有頻率的影響更為顯著。

圖5給出了不同h/R時環肋圓柱殼的固有頻率曲線。如圖所示，隨著周向波數的逐漸增大，殼體的固有頻率先減小后增大;固有頻率隨著h/R的增大而增大，且h/R對固有頻率的影響主要表現在周向波數較大的情況下。

圖2 不同軸向半波數m時圓柱殼的固有頻率的變化規律曲線Fig.2 Variation of natural frequencies of cylindrical shells for some different axial half wave numbers

圖3 不同肋條數目時圓柱殼的固有頻率的變化規律曲線Fig.3 Variation of natural frequencies ofring-stiffened cylindrical shells for the number of Nr

圖4 不同hr/br時圓柱殼的固有頻率的變化規律曲線Fig.4 Variation of natural frequencies for a ring-stiffened cylindrical shells with different hr/br ratios

圖5 不同h/R時環肋圓柱殼的固有頻率變化曲線Fig.5 Variation of natural frequency for a submerged ring-stiffened cylindrical shell with different h/R ratios

圖6 不同軸向半波數m時圓柱殼的耦合頻率和非耦合頻率的變化規律曲線Fig.6 Variation of coupled and uncoupled frequencies of cylindrical shells for some different axial half wave numbers

圖7 兩端簡支邊界條件下水下圓柱殼的固有頻率變化曲線Fig.7 Variation for a simply support submerged cylindrical shell

圖6 －10給出了考慮液體影響時環肋圓柱殼的固有頻率變化規律曲線，為了更好表現流體因素對殼體固有頻率的影響。圖6給出了僅考慮流體影響，靜水壓力Q=0時的固有頻率曲線。可以看出聲壓對殼體的固有頻率有顯著的影響;隨著周向波數的逐漸增大，耦合頻率逐漸由非耦合頻率的約1/5增大至非耦合頻率的1/3;軸向半波數對固有頻率的影響主要表現在周向波數較小的情況下。

選用h/R=0.01。圖7給出了靜水壓力下，兩端簡支環肋圓柱殼的固有頻率曲線。可以看出，隨著靜水壓力的逐漸增大，殼體的固有頻率逐漸減小;固有頻率隨著周向波數的增大而增大;當n=2，Q=80 kPa時，固有頻率接近0，此時的靜水壓力即為殼體的臨界壓力。

基于上圖結論，接下來將n=2時殼體固有頻率作為研究對象。圖8給出了靜水壓力下，不同肋條數目時殼體的固有頻率。從圖中可以看出，隨著靜水壓力的不斷增大，殼體的固有頻率逐漸減小;不同肋條數目時殼體的固有頻率非常接近，即肋條數目對殼體的固有頻率的影響并不明顯，當靜水壓力接近80 kPa，固有頻率下降速度明顯加快，不同固有頻率之間的差別開始顯現;當Nr=0時，殼體的固有頻率下降速度變快。

圖9給出了靜水壓力下，不同邊界條件時環肋圓柱殼的固有頻率變化曲線。可以看出，殼體的固有頻率隨著靜水壓力的逐漸增大而減小，接近臨界壓力時固有頻率下降速度變快;五種邊界條件中，兩端固定邊界條件下的固有頻率最大，一端自由一端滑移邊界條件下的固有頻率最小;一端自由一端滑移和一端固定一端自由這兩種邊界條件下的固有頻率非常接近;當靜水壓力為50 kPa時，一端自由一端滑移邊界條件下的固有頻率接近0，達到臨界靜水壓力。

圖10給出了靜水壓力下，不同h/R時殼體的固有頻率。由圖可知，殼體厚度對固有頻率的影響非常顯著;隨著靜水壓力的增大，殼體的固有頻率逐漸減小，當h/R=0.002時固有頻率下降速度最快;殼體的固有頻率隨著h/R的增大而增大，且固有頻率的下降速度隨著h/R的增大而迅速變緩。

圖8 靜水壓力下不同肋條數目時圓柱殼的固有頻率的變化曲線Fig.8 Variation of natural frequencies of ring-stiffened cylindrical shells for the number of Nr under hydrostatic pressure

圖9 不同邊界條件下環肋圓柱殼的固有頻率曲線Fig.9 Variation of natural frequency for a submerged ring-stiffened cylindrical shell with different boundary conditions

圖10 不同h/R時水下環肋圓柱殼的固有頻率變化曲線Fig.10 Variation of natural frequency for a submerged ring-stiffened cylindrical shell with different h/R ratios under hydrostatic pressure

4 結論

本文基于Flügge理論并考慮流體影響，采用波動法研究了靜水壓力下環肋圓柱殼的耦合振動特性。

(1)環肋圓柱殼的固有頻率明顯大于不加肋圓柱殼的固有頻率;肋條對固有頻率的影響主要表現在周向波數較大的情況下;軸向半波數對固有頻率的影響主要表現在周向波數較小的情況下。

(2)肋條數目對固有頻率的影響主要表現在周向波數較大的情況下;殼體的固有頻率隨著肋條截面長度和寬度比的增大而增大，且其比例值越大，頻率增長幅度約明顯;長度和寬度比對固有頻率的影響主要表現在周向波數較大的情況下;與肋條數目(肋條間距)相比，肋條的截面長度和寬度比對固有頻率的影響更為顯著。

(3)固有頻率隨著h/R的增大而增大，且h/R對固有頻率的影響主要表現在周向波數較大的情況下;當考慮流體影響時，隨著靜水壓力的增大，殼體的固有頻率逐漸減小，當h/R=0.002時固有頻率下降速度最快;殼體的固有頻率隨著h/R的增大而增大，且固有頻率的下降速度隨著h/R的增大而迅速變緩。

(4)隨著靜水壓力的逐漸增大，殼體的固有頻率逐漸減小;當靜水壓力接近臨近壓力時，固有頻率下降速度明顯加快。

(5)五種邊界條件中，兩端固定邊界條件下的固有頻率最大，一端自由一端滑移邊界條件下的固有頻率最小;一端自由一端滑移和一端固定一端自由這兩種邊界條件下的固有頻率非常接近;當靜水壓力為50 kPa時，一端自由一端滑移邊界條件下的固有頻率接近0，達到臨界靜水壓力。

［1］Junger M C.Radiation loading of cylindrical and spherical surfaces［J］.Acoust.Soc.Am.，1952，24(3):288 －289.

［2］Junger M C.Vibrations of elastic shells in a fluid medium and the associated radiation of sound ［J］.Appl.Mech.，1952，19(3):439 －445.

［3］Zhang X M，Liu G R，Lam K Y.Coupled vibration analysis of fluid-filled cylindrical shells using the wave propagation approach［J］.Applied Acoustics，2001，62(3):229 －243.

［4］Zhang X M.Frequency analysis of submerged cylindrical shells with the wave propagation approach［J］.International Journal of Mechanical Sciences，2002，44(7):1259－1273.

［5］Loy C T，Lam K Y.Vibration of cylindrical shells with ring support［J］.International Journal of Mechanical Sciences，1997，39(4):455 －471.

［6］ Liu Z Z，Li T Y，Zhu X，et al.Effect of hydrostatic pressure on input power flow in submerged ring-stiffened cylindrical shells［J］.Journal of Ship Mechanics.2011，15(3):301－312.

［7］ Gan L，Li X B，Zhang Z.Free vibration analysis of ringstiffened cylindrical shells using wave propagation approach［J］.Journal of Sound and Vibration，2009，326(3 －5):633－646.

［8］Vodenitcharova T， Ansourian P. Buckling of circular cylindrical shells subject to uniform lateral pressure［J］.Eng.Struct.，1996，18(8):604 －614.

［9］Zhu X，Ye W B，Li T Y，et al，The elastic critical pressure prediction of submerged cylindrical shell using wave propagation method［J］.Ocean Engineering，2013，58:22－26.

［10］ Flügge W.Stresses in shells［M］.Second Edition.Springer -Verlag，New York，1973.

［11］Zhang X M，Liu G R，Lam K Y.Coupled vibration analysis of fluid-filled cylindrical shells using the wave propagation approach［J］.Applied Acoustics，2001，62(3):229 －243.

［12］Jafari A A，Bagheri M.Free vibration of non-uniformly ring stiffened cylindrical shells using analytical，experimental and numerical methods ［J］.Thin-Walled Structures，2006，44(1):82－90.