曹雄濤，華宏星

（上海交通大學(xué) 機(jī)械系統(tǒng)與振動(dòng)國家重點(diǎn)實(shí)驗(yàn)室，上海 200240）

圓錐殼體在很多工程應(yīng)用中被廣泛使用。在航空航天飛行器、潛艇中，圓錐殼體是一個(gè)常用的結(jié)構(gòu)單元。Rayleigh-Ritz法被用來研究圓錐殼的自由振動(dòng)［1，2］。Galerkin 法 被 用 來 研 究 圓 錐 殼 的 振 動(dòng) 特性［3－7］。其圓錐殼體的模型是基于 Love-Timoshenko的薄殼理論。Shu［8］運(yùn)用廣義微分積分法研究了層合圓錐殼的自由振動(dòng)。Liew［9］用kp－Ritz法分析了薄錐殼的自由振動(dòng)，討論了不同邊界條件下的頻率特性。Crenwelge和 Muster［10］研究了正交加筋圓錐殼體的振動(dòng)，這是早期的加筋錐殼的研究工作。能量法和假設(shè)模態(tài)法被用來求解縱橫加筋殼體的振動(dòng)。一些研究人員［1，11］也用這種方法分析了正交加筋復(fù)合殼體的自由振動(dòng)。這種方法的不便之處在于控制方程不能顯式獲得。Mecitoglu［12，13］研究了加筋各向同性錐殼的自由振動(dòng)。通過平均能量法，環(huán)筋和縱骨的慣性和剛度融合到錐殼里面。一個(gè)簡(jiǎn)單的正交各向異性圓錐殼得到了。由于這種理論的局限性，這種模型僅在低頻范圍內(nèi)是有效的。融合筋法（smeared theory）也被Goldfeld［14］用來研究螺旋筋加強(qiáng)錐殼的屈曲穩(wěn)定性分析。Tong［15］研究了復(fù)合錐殼的自由振動(dòng)，解是通過冪級(jí)數(shù)法獲得的。但是其結(jié)果并不是精確的解析解，盡管有較好的精度。Caresta和Kessissoglou［16］分析了各向同性錐 － 柱組合殼的自由振動(dòng)特性。錐殼通過冪級(jí)數(shù)法求解。在周向波數(shù)為1，自由－自由邊界條件下，Donnell-Mushtari殼體方程產(chǎn)生錯(cuò)誤的頻率。傳遞矩陣法［17.18］被用來研究錐殼的固有頻率。有限元法［19－21］被用來分析錐殼的自由振動(dòng)。Sofiyev［22］研究了功能梯度（FGM）錐殼的穩(wěn)定性和非線性振動(dòng)。

圓錐殼體聲輻射的研究工作比較少。Guo［23］利用多尺度技術(shù)研究了彈性波在流體載荷下錐殼中的傳播。將流體載荷作用下的圓柱殼的控制方程通過攝動(dòng)來研究錐殼的動(dòng)力特性。這種研究?jī)H僅是一種漸近分析。Caresta和Kessissoglou［24］用冪級(jí)數(shù)序列求解了流體作用下錐殼的動(dòng)力響應(yīng)。流體通過將圓錐分段成小錐殼，這些小錐殼被當(dāng)作圓柱殼，然后確定每一小段圓柱的流體載荷。然而，對(duì)于母線方向以冪級(jí)數(shù)序列表示的位移場(chǎng)，他們的模型中，沒有給出母線方向的波數(shù)km是如何確定的。冪級(jí)數(shù)序列似乎不適合于聲輻射，因?yàn)閮缂?jí)數(shù)序列不能直接同波數(shù)相關(guān)。而且對(duì)于流體載荷，km是至關(guān)重要的，這一點(diǎn)能從本文中看到。

在本文研究中，波傳播法被用來分析流體加載下的正交各向異性圓錐殼體的自由振動(dòng)。波傳播法被一些研究人員用來分析圓柱殼的振動(dòng)特性［25－27］。根據(jù)邊界條件，母線方向的波數(shù)km是已知的。這種性質(zhì)對(duì)于處理流體載荷是方便的。圓錐的流體載荷用多個(gè)小段圓柱的流體載荷來逼近。應(yīng)用殼體和流體間的界面邊界條件和求解圓柱坐標(biāo)系下的Helmholtz方程，就分段獲得了圓錐殼體的流體載荷。給出了確定流體載荷的詳細(xì)過程。圓錐殼體的運(yùn)動(dòng)方程是基于Reissner－Naghdi薄殼理論，并且利用Galerkin法求解。研究了圓錐的自由振動(dòng)和流體載荷對(duì)圓錐的影響。對(duì)于穩(wěn)態(tài)振動(dòng)，在本文中時(shí)間因子e－iωt略去不寫。

1 圓錐殼體的運(yùn)動(dòng)方程

一個(gè)薄壁正交各向異性圓錐殼顯示在圖1中。l1和l分別是圓錐殼體小端母線方向的長度和大端母線方向的長度。u1，u2和u3是殼體中面在母線方向α1，圓周方向α2和法向方向α3的位移。曲線坐標(biāo)和位移場(chǎng)顯示在圖1中。

圖1 圓錐殼體曲線坐標(biāo)和中面位移Fig.1 The curvilinear coordinates and displacements of middle surface for conical shell

其中，q1，q2，q3和pa分別是表面力載荷作用在α1，α2，α3方向和流體載荷。Lij是微分算子，它們是：

其中，I1是正交各向異性圓錐殼體的面密度，其值為I1=ρch。ρc和h是圓錐殼的體積密度和厚度。α是圓錐殼體的半錐角。正交各向異性圓錐殼體的材料參數(shù)Aij和Dij定義如下：

其中，E1，E2，ν12，ν21和G12是正交各向異性圓錐殼體材料的彈性常數(shù)。正交各向異性材料的Poisson比有如下關(guān)系：

正交各向異性圓錐殼體的位移場(chǎng)用波傳播的形式可以寫為：

其中，km是母線方向波傳播因子，其值由邊界條件決定。通常km能選為與錐殼邊界條件相似的Euler梁的波傳播因子［27］。通過把方程（4）代入（1），方程（1）是解耦的對(duì)于參數(shù)α2，但對(duì)于參數(shù)α1是耦合的。通過運(yùn)用Galerkin法到方程（1），加權(quán)積分的Galerkin過程能寫為：

像Lam［7］那樣處理：在方程（5）中的耦合項(xiàng)被忽略，例如：cos［km（α1－l1）］cos［kd（α1－l1）］和sin［km（α1－l1）］sin［kd（α1－l1）］，m≠d。這些耦合項(xiàng)的貢獻(xiàn)非常小。

2 流體載荷和積分過程

Moon和Spencer［28］用分離變量法研究過Helmholtz方程的解，指出Helmholtz方程能在40種正交曲線坐標(biāo)系下用分離變量法求解。Helmholtz方程在錐殼的正交曲線坐標(biāo)系中是不能用分離變量法來處理的。因此，一個(gè)分段流體加載方法用來近似圓錐殼體的流體載荷，如圖2所示。

2.1 流體載荷

流體介質(zhì)中的聲壓滿足柱坐標(biāo)系下的Helmholtz方程：

面向分布式光伏虛擬集群的有源配電網(wǎng)多級(jí)調(diào)控//竇曉波,常莉敏,倪春花,段向梅,葛浦東,吳在軍//(3):21

其中，流體域中的波數(shù)k0為ω/c，c是流體中的聲速。Laplace算子為：

其中，α1，α2和α3分別是軸向、周向和徑向坐標(biāo)。第j個(gè)圓錐段（被近似為圓柱殼段）與流體間的邊界條件為：

將方程（6），方程（8）和方程組（4）的第三個(gè)方程聯(lián)立起來，運(yùn)用Sommerfeld遠(yuǎn)場(chǎng)輻射條件，就獲得了第j個(gè)圓錐段的流體載荷：

圖2 圓錐殼段的劃分和一個(gè)流體加載的圓柱殼Fig.2 Strips of conical shells and a strip of fluid loaded circular cylindrical shell

2.2 殼體表面壓力在Galerkin法中的處理

如果圓錐殼體被劃分為s段圓柱殼，錐殼的流體載荷pa由s段以方程（9）的形式表示的流體載荷組成。運(yùn)用方程（9），流體載荷在方程（5）中的貢獻(xiàn)pmnε可以寫為：

式中的積分很容易獲得。其中，一些耦合項(xiàng)像前面的處理那樣被忽略。這些耦合項(xiàng)的貢獻(xiàn)很小。

3 流體加載圓錐殼體的自由振動(dòng)

設(shè)外力q1，q2和q3為零，將方程（4），方程（9）代入方程（5），得到：

其中，系數(shù)Cij是非常復(fù)雜的表達(dá)式，涉及到材料參數(shù)、角頻率和殼體的幾何參數(shù)。僅僅系數(shù)C11在附錄中給出，其它系數(shù)可以用相似的方法得到。對(duì)于真空中的圓錐殼，特征方程（11）能夠用線性代數(shù)的方法求解；對(duì)于流體加載的圓錐殼，用迭代法求解。

4 數(shù)值結(jié)果

在本節(jié)中，利用前面的理論，研究了各向同性和正交各向異性圓錐殼體的自由振動(dòng)。對(duì)于方程（11），有三類不同的固有頻率。最低的頻率對(duì)應(yīng)錐殼的徑向運(yùn)動(dòng)，其它兩個(gè)固有頻率主要代表縱向運(yùn)動(dòng)或周向運(yùn)動(dòng)。無量綱頻率定義為：

比率γ12定義為：

固有頻率以ωc的形式給出，并且給出的頻率對(duì)應(yīng)錐殼徑向振動(dòng)模態(tài)。錐殼浸沒在水中。水中聲速為c=1 500 m/s，質(zhì)量密度為ρ=1 000 kg/m3。錐殼被劃分為s=20段圓柱殼來逼近錐殼的流體載荷。

4.1 模型的有效性驗(yàn)證

將頻率參數(shù)同Lam［6］給出的結(jié)果作了比較。各向同性錐殼的參數(shù)為h/b=0.01，ν=0.3，（l－l1）sinα/b=0.25。表1給出了真空中簡(jiǎn)支錐殼的固有頻率。母線方向的波數(shù)為mπ/（l－l1）。

表1 簡(jiǎn)支錐殼的頻率 ωc，h/b=0.01，ν=0.3，（l－l1）sinα/b=0.25，m=1Tab.1 The nondimensional frequency ωcof a simply supported conical shell h/b=0.01，ν =0.3，（l－l1）sinα/b=0.25，m=1

可以看出當(dāng)前給出的結(jié)果同Lam給出的結(jié)果符合得很好。錐殼的頻率隨周向波數(shù)的變化規(guī)律為：初始時(shí)逐漸減小，然后逐漸增加。

下一個(gè)例子用來比較流體載荷下很小錐角的圓錐殼和相似結(jié)構(gòu)圓柱殼的固有頻率。張［26］利用波傳播法分析了流體載荷作用下簡(jiǎn)支圓柱殼的固有頻率。圓柱殼的參數(shù)如下：長度、半徑和厚度分別是20 m，1 m和0.01 m；質(zhì)量密度、Poisson比和彈性模量分別為：7 850 kg/m3，0.3 和2.1 ×1011Pa。圓錐小端半徑為 1 m，大端半徑為1.04 m。錐殼半錐角為0.114 591 6°。l1為500 m，l為520 m。圓錐在母線方向的長度為20 m。材料參數(shù)同圓柱一樣。錐殼和柱殼都浸沒在水中。圖3給出了流體載荷作用下的錐殼和柱殼的固有頻率。對(duì)比顯示小錐角錐殼與柱殼的流固耦合振動(dòng)固有頻率符合得很好。表明當(dāng)前方法處理流體載荷是有效的。圖3顯示：在低周向波數(shù)時(shí)，流體中的錐殼的固有頻率下降了對(duì)應(yīng)于真空中錐殼的頻率大約一半。

4.2 各向同性圓錐殼體的振動(dòng)

對(duì)于各向同性圓錐，彈性模量E為2.1×1011Pa，泊松比ν為0.3，密度ρc為7 800 kg/m3和厚度h為0.006 m。錐殼小端半徑a是0.2 m，大端半徑b是0.8 m。l1和l分別是0.66 m 和2.667 m。錐角為34.916°。錐殼兩端邊界條件為簡(jiǎn)支－簡(jiǎn)支（SS）。

圖4顯示了SS邊界條件下錐殼的耦合和未耦合固有頻率。未耦合和耦合模態(tài)分別對(duì)應(yīng)真空中的錐殼和浸沒在流體中的錐殼。從圖中能夠看出流體和錐殼是強(qiáng)烈耦合的。在低周向波數(shù)時(shí)耦合頻率幾乎下降到未耦合頻率的一半。在高周向波數(shù)時(shí)，耦合頻率大約下降了四分之一。在分析流體與錐殼相互作用時(shí)，流體的影響必須要考慮。對(duì)于兩端簡(jiǎn)支的錐殼，圖5顯示了錐殼分段數(shù)s對(duì)耦合頻率收斂特性的影響。s=5與s=20的頻率值是相同的，表明頻率已經(jīng)收斂。耦合頻率收斂速度很快，s的取值比較靈活。

4.3 正交各向異性圓錐殼體的自由振動(dòng)

對(duì)于正交各向異性圓錐殼，在計(jì)算中使用兩個(gè)殼體模型。兩個(gè)模型的材料參數(shù)是E1=E=2.1×1011Pa，E2=E1/γ12，ν12=0.3，G1=E1/［2（1+ν12）］和ρc=7 800 kg/m3。流體介質(zhì)是水。兩個(gè)模型中的邊界條件是SS。第一個(gè)錐殼模型的幾何參數(shù)像4.2節(jié)那樣保持不變。第二個(gè)錐殼模型的幾何參數(shù)與第一個(gè)模型相同，只是下面幾個(gè)參數(shù)改變了：錐殼大端半徑b為1.2 m，l1為0.4 m，l為2.4 m，錐角為60°。圖6－圖9顯示了比率γ12對(duì)有無流體加載下錐殼固有頻率的影響。通過比較圖6－圖7和圖8－圖9，可以看到流體載荷改變了錐殼的自由振動(dòng)特征。在低周向波數(shù)時(shí)，流體載荷使得錐殼的頻率下降了一半；在高周向波數(shù)時(shí)，錐殼的頻率下降了四分之一。低周向波數(shù)時(shí)，流體載荷強(qiáng)烈影響著錐殼的動(dòng)力特性。當(dāng)比率γ12從2增加到20時(shí)，圖6－圖9顯示ωc下降，并且周向波數(shù)的增加輕微地改變了ωc。這是因?yàn)楫?dāng)γ12增加時(shí)，與母線方向剛度相比，周向剛度變?nèi)趿?。這種變?nèi)醯闹芟騽偠确潘闪藢?duì)ωc的貢獻(xiàn)。對(duì)于大的γ12，ωc主要依賴于母線方向剛度。

圖6 γ12對(duì)ωc的影響，半錐角為α=17.458°，沒有流體載荷，m=1Fig.6 Effect of γ12on ωcfor conical shell with half vertex angle α =17.458°，without fluid load，m=1

圖7 γ12對(duì)ωc的影響，半錐角為 α =17.458°，有流體載荷，m=1Fig.7 Effect of γ12on ωcfor fluid loaded conical shell with half vertex angle α =17.458°，m=1

圖8 γ12對(duì)ωc的影響，半錐角為α=30°，沒有流體載荷，m=1Fig.8 Effect of γ12on ωcfor conical shell with half vertex angle α =30°，without fluid load，m=1

圖9 γ12對(duì) ωc的影響，半錐角為α=30°，有流體載荷，m=1Fig.9 Effect of γ12on ωcfor fluid loaded conical shell with half vertex angle α =30°，m=1

5 結(jié)論

給出一個(gè)包括流體載荷作用的圓錐殼自由振動(dòng)理論模型。一個(gè)技巧而又系統(tǒng)的錐殼自由振動(dòng)的解通過波傳播法和Galerkin法得到。將圓錐殼分段用圓柱殼近似，流體載荷易于處理。波傳播法適合于確定流體與錐殼的耦合效應(yīng)。因?yàn)殄F殼和聲場(chǎng)都體現(xiàn)了波動(dòng)。本文的數(shù)值結(jié)果同其他研究者的結(jié)果相符，表明了給出的模型是有效的。研究了不同邊界條件下錐殼自由振動(dòng)特征。分析了彈性模量比率γ12對(duì)正交各向異性圓錐殼體頻率ωc的影響。γ12越大，ωc就越小。數(shù)值結(jié)果顯示流體載荷強(qiáng)烈影響結(jié)構(gòu)的動(dòng)力特性。在低周向波數(shù)時(shí)，錐殼的耦合頻率下降了大約一半；高周向波數(shù)時(shí)，錐殼的耦合頻率下降了大約四分之一。本文的研究可以應(yīng)用到潛艇艉部錐殼的聲輻射上。

［1］ Talebitooti M，Ghayour M，Ziaei-Rad S，et al.Free vibrations of rotating composite conical shells with stringer and ring stiffeners［J］.Archive of Applied Mechanics，2010，8：1 －15.

［2］ Li F M，Kishimoto K，Huang W H.The calculations of natural frequencies and forced vibration responses of conical shell using the rayleigh-ritz method［J］.Mechanics Research Communications，2009，36（5）：595 －602.

［3］Hua L.Influence of boundary conditions on the free vibrations of rotating truncated circular multi-layered conical shells［J］.Composites Part B：Engineering，2000，31（4）：265 －275.

［4］Hua L.Frequency characteristics of a rotating truncated circular layered conical shell［J］.Composite Structures，2000，50（1）：59－68.

［5］ Lam K Y，Hua L.Vibration analysis of a rotating truncated circular conical shell［J］.International Journal of Solids and Structures，1997，34（17）：2183 －2197.

［6］ Lam K Y，Hua L.Influence of boundary conditions on the frequency characteristics of a rotating truncated circular conical shell［J］.Journal of Sound and Vibration，1999，223（2）：171－195.

［7］ Lam K Y，Hua L.On free vibration of a rotating truncated circular orthotropic conical shell［J］.Composites Part B：Engineering，1999，30（2）：135 －144.

［8］Shu C.Free vibration analysis of composite laminated conical shells by generalized differential quadrature［J］.Journal of Sound and Vibration，1996，194（4）：587 －604.

［9］ Liew K M，Ng T Y，Zhao X.Free vibration analysis of conical shells via the element-free kp-ritz method［J］.Journal of Sound and Vibration，2005，281：3 －5，627 －645.

［10］ Crenwelge Jr O E，Muster D.Free vibrations of ring-andstringer-stiffened conical shells［J］.Journal of Acoustical Soc America，1969，46（1，2）：176 －185.

［11］Zhao X，Liew K M，Ng T Y.Vibrations of rotating cross-ply laminated circular cylindrical shells with stringer and ring stiffeners［J］.International Journal of Solids and Structures，2001，39（2）：529－545.

［12］ Mecitoglu Z.Vibration characteristics of a stiffened conical shell［J］.Journal of Sound and Vibration，1996，197（2）：191－206.

［13］ Mecitoglu Z.Governing equations of a stiffened laminated inhomogeneous conical shell［J］.AIAA Journal，1996，34（10）：2118－2125.

［14］ Goldfeld Y.Elastic buckling and imperfection sensitivity of generally stiffened conical shells［J］.AIAA Journal，2007，45（3）：721－729.

［15］ Tong L.Free vibration of composite laminated conical shells［J］.International Journal of Mechanical Sciences，1993，35（1）：47－61.

［16］ Caresta M，Kessissoglou N J.Free vibrational characteristics of isotropic coupled cylindrical-conical shells［J］.Journal of Sound and Vibration，2010，329（6）：733 －751.

［17］ Irie T，Yamada G，Muramoto Y.Free vibration of joined conical-cylindricalshells［J］. JournalofSound and Vibration，1984，95（1）：31 －39.

［18］ Irie T，Yamada G，Tanaka K.Natural frequencies of truncated conical shells［J］.Journal of Sound and Vibration，1984，92（3）：447－453.

［19］PatelB P， GanapathiM，KamatS. Free vibration characteristics of laminated composite joined conicalcylindrical shells［J］.Journal of Sound and Vibration，2000，237（5）：920－930.

［20］ Buchanan G R，Wong F I.Frequencies and mode shapes for thick truncated hollow cones［J］.International Journal of Mechanical Sciences，2001，43（12）：2815 －2832.

［21］ El Damatty A A，Saafan M S，Sweedan A M I.Dynamic characteristics of combined conical-cylindrical shells［J］.Thin-Walled Structures，2005，43（9）：1380 －1397.

［22］ Sofiyev A H，Kuruoglu N，Halilov H M.The vibration and stability of non-homogeneous orthotropic conical shells with clamped edges subjected to uniform external pressures［J］.Applied Mathematical Modelling，2010，34：1807 －1822.

［23］ Guo Y P.Fluid-loading effects on waves on conical shells［J］.Journal of Acoustical Society of Americ，1994，97（2）：1061－1066.

［24］ Caresta M，Kessissoglou N J.Vibration of fluid loaded conical shells［J］.Journal of the Acoustical Society of America，2008，124（4）：2068－2077.

［25］ Gan L，Li X，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.

［26］ 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.

［27］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.

［28］ Parry Moon D E S.Field theory handbook［M］.Berlin Springer-Verlag，1988.

附錄同其它系數(shù)Cij相比，系數(shù)C11有最簡(jiǎn)單的表達(dá)式。其它系數(shù)Cij能夠用相似的方法得到。

其中，sine積分函數(shù)si（x）定義為：

cosine積分函數(shù)ci（x）定義為：