王新杰 岳洪浩 鄧宗全

（哈爾濱工業(yè)大學(xué) 機(jī)電工程學(xué)院，哈爾濱 150001）

自由邊界拋物薄殼的模態(tài)振型函數(shù)

王新杰 岳洪浩 鄧宗全

（哈爾濱工業(yè)大學(xué) 機(jī)電工程學(xué)院，哈爾濱 150001）

基于membrane理論，推導(dǎo)出一組完全滿足拋物薄殼自由邊界條件的模態(tài)振型函數(shù).通過模態(tài)分析實(shí)驗(yàn)得到自由邊界拋物薄殼的低階模態(tài)實(shí)驗(yàn)振型，將實(shí)驗(yàn)振型與有限元分析振型及依據(jù)模態(tài)振型函數(shù)所得到的理論振型進(jìn)行對比，三者能夠較好的吻合，表明所推導(dǎo)的模態(tài)振型函數(shù)可正確描述拋物薄殼的模態(tài)形狀.

拋物薄殼;模態(tài)振型函數(shù);membrane理論;模態(tài)分析

拋物殼結(jié)構(gòu)具有聚焦特性和良好的空氣動力學(xué)性能，在航空航天領(lǐng)域得到廣泛的應(yīng)用，如大型薄壁通信天線、太空光學(xué)反射鏡和火箭整流罩等.而拋物殼結(jié)構(gòu)的振動抑制問題，成為工程應(yīng)用和設(shè)計(jì)中的難點(diǎn)和重點(diǎn).近年來眾多學(xué)者通過智能材料對拋物殼結(jié)構(gòu)進(jìn)行主動振動控制從而達(dá)到抑制振動的目的［1－2］.而拋物殼結(jié)構(gòu)的模態(tài)振型反映了結(jié)構(gòu)對外載荷響應(yīng)的相對幅度、相位和振動形狀，是研究板殼結(jié)構(gòu)智能監(jiān)測、傳感/作動器設(shè)計(jì)與布局的重要依據(jù)和基礎(chǔ).

針對拋物殼的動力學(xué)問題的研究起步較早，眾多學(xué)者作了大量的研究工作［3－9］，但針對拋物殼模態(tài)振型所開展的研究則相對較少.1960年，Lin和Lee推導(dǎo)了基于彎矩假設(shè)的自由邊界條件下的模態(tài)形狀方程［10］，但實(shí)際上該模態(tài)形狀方程只滿足彎矩假設(shè)條件而并不滿足自由邊界條件.1964年，Hoppmann等人采用實(shí)驗(yàn)的方法描繪出不同邊界條件下拋物殼的振動節(jié)線［11］.1973年，Glockner和Tawadros對拋物殼在簡支和固支條件下的自由振動進(jìn)行了實(shí)驗(yàn)，并給出拋物殼模態(tài)形狀的離散實(shí)驗(yàn)值［12］.2005 年，Kang和 Leissa應(yīng)用三維Ritz方法分析了不同厚度拋物殼彈性結(jié)構(gòu)的自由振動特性，利用包含有厚度參數(shù)在內(nèi)的多項(xiàng)式來表示拋物殼的模態(tài)形狀函數(shù)，但并沒有給出拋物殼在相應(yīng)邊界條件下模態(tài)形狀函數(shù)的具體形式［13］.綜上所述，針對自由邊界拋物殼的模態(tài)振型函數(shù)仍需進(jìn)一步的理論研究，開展對應(yīng)的模態(tài)實(shí)驗(yàn)分析，進(jìn)行驗(yàn)證非常有必要.

本文以自由邊界拋物薄殼的模態(tài)振型為研究對象.首先基于membrane簡化，對自由邊界條件下拋物薄殼的模態(tài)振型函數(shù)進(jìn)行理論推導(dǎo);然后通過模態(tài)分析實(shí)驗(yàn)，對所推導(dǎo)出的模態(tài)振型函數(shù)進(jìn)行驗(yàn)證.

1 拋物殼的薄膜力與薄膜剪力

定義拋物殼的曲線正交坐標(biāo)系為（φ，ψ，α3）.φ，ψ和α3分別表示拋物殼的子午線方向、圓周方向和橫截方向，如圖1所示.用c和a分別表示拋物殼的高度及開口半徑;并引入?yún)?shù)b，且b=a2/（2c）.φ*表示子午線方向的邊界弧度，且φ*=arctan（2c/a）.則子午線曲率半徑 Rφ=b/cos3φ，圓周曲率半徑Rψ=b/cosφ.

圖1 拋物殼的曲線坐標(biāo)系

設(shè)拋物殼中曲面在φ，ψ和α3方向的位移分量分別為uφ，uψ和u3，則拋物殼薄膜力、薄膜剪力與位移的關(guān)系可表示為［14］

2 自由邊界拋物薄殼模態(tài)振型函數(shù)

當(dāng)殼體的厚度與曲率半徑的比值小于0.05時(shí)，可采用 membrane理論［15］，忽略殼體的彎矩、扭矩及橫向剪力.本節(jié)將采用membrane理論對自由邊界拋物薄殼的模態(tài)振型函數(shù)進(jìn)行推導(dǎo).

對于自由邊界拋物殼，其邊界上的所有力與力矩都應(yīng)為零，由于本文采用了membrane理論，忽略了殼體的彎矩、扭矩及橫向剪力，因此邊界約束條件簡化為邊界上（即φ=φ*時(shí)）沿著子午線和圓周方向的薄膜力以及薄膜剪力為零，即

依據(jù)式（1）～式（3），如要滿足約束方程（4），可令

假設(shè)彈性拋物殼上任意點(diǎn)以一個(gè)固有頻率進(jìn)行簡諧振動，而其位移可表示為以下形式:

式中，k表示模態(tài)階數(shù);ωk表示第k階固有頻率;Uφk，Uψk和 U3k分別表示子午線方向、圓周方向和橫截方向的模態(tài)振型函數(shù).

將式（8）～式（10）代入式（5）～式（7）可得

首先將橫截方向的模態(tài)振型函數(shù)寫為

式中f3（φ）為變量φ的函數(shù).

將式（14）代入式（11）得:

將式（14）和式（15）代入式（12）得:

式（18）的通解為

式中c1，c2為任意常數(shù).

由于拋物殼的不同階次振型的最大幅值并不相同，因此這里引入模態(tài)的振動幅值系數(shù)Ak，則3個(gè)方向上的模態(tài)振型函數(shù)可以表示為

為了保證當(dāng) φ=φ*時(shí)，式（5）～式（7）成立，Uφk修改為

由式（19）可以看出，F(xiàn)3（φ）=sin φtan－kφ 和F3（φ）=sinφtankφ均為方程式（18）的特解.

當(dāng) F3（φ）=sin φtan－kφ 時(shí)，依據(jù)式（21）～ 式（23）可求得

針對式（26）的表達(dá)式來說，當(dāng)φ=0時(shí)，U3k的值為無窮大，當(dāng)φ由0→φ*時(shí)，tan－kφ越來越小，U3k的值逐漸減少.也就是說，殼的頂點(diǎn)（φ=0）處，其振動幅值為無窮大;而沿著子午線方向由頂點(diǎn)到邊緣（0→φ ）移動時(shí)，其振動幅值是逐漸減小的.這些均不符合自由邊界條件下拋物殼振動的實(shí)際情況，因此該組振型函數(shù)被排除掉.

當(dāng) F3（φ）=sin φtankφ 時(shí)，依據(jù)式（21）～ 式（23）可求得:

經(jīng)分析，式（29）所描述的振型函數(shù)符合自由邊界條件下拋物殼振動的實(shí)際情況.依據(jù)式（27）～式（29）所表述的振型函數(shù)可得到拋物薄殼各階理論振型的3D圖形.

3 拋物薄殼模態(tài)分析實(shí)驗(yàn)

3.1 拋物薄殼物理模型及模態(tài)分析測試平臺

本文以DSM Somos-14120樹脂為原料，采用激光快速成型法制作了拋物薄殼的物理模型.DSM Somos-14120樹脂的機(jī)械性能指標(biāo)為:彈性模量 Y=2.8×109N/m2，質(zhì)量密度 ρ=1 120 kg/m3（25℃），泊松比 μ=0.35.拋物薄殼的尺寸參數(shù)為:開口半徑為 a=0.2 m，高度 c=0.1 m，厚度h=1×10－3m.

為滿足自由邊界條件，采用橡皮筋將拋物薄殼懸吊起來，將殼體的頂點(diǎn)確定為懸吊點(diǎn)，以使懸吊方式對邊界條件的影響減到最小.此外，加速度傳感器的質(zhì)量和數(shù)量也會對精密柔性殼結(jié)構(gòu)的動態(tài)特性產(chǎn)生重要影響，所以實(shí)驗(yàn)中只設(shè)置一個(gè)響應(yīng)參考點(diǎn)，并選擇ICP加速度計(jì)（M353B18）對響應(yīng)參考點(diǎn)進(jìn)行測量.利用ICP力錘對拋物薄殼進(jìn)行激勵，通過NI動態(tài)信號采集卡（NI-PXI44 72）采集數(shù)據(jù).圖2所示為模態(tài)分析測試平臺.

圖2 自由邊界拋物薄殼模態(tài)分析測試平臺

3.2 模態(tài)振型測試及分析

在拋物薄殼模型上沿圓周方向和子午線方向等弧度交叉劃分48個(gè)節(jié)點(diǎn)，應(yīng)用移動激勵點(diǎn)法測試實(shí)驗(yàn)?zāi)Ｐ偷膭討B(tài)特性，并通過ModalVIEW模態(tài)分析軟件進(jìn)行頻響函數(shù)分析，最終提取自由邊界條件下拋物殼的低階振型.將拋物殼前三階實(shí)驗(yàn)振型與理論振型、ANSYS分析振型對比，如圖3～圖5所示.

圖3 自由邊界拋物薄殼第一階模態(tài)形狀對比

圖4 自由邊界拋物薄殼第二階模態(tài)形狀對比

圖5 自由邊界拋物薄殼第三階模態(tài)形狀對比

通過模態(tài)振型比較可以看出，理論振型、實(shí)驗(yàn)振型和ANSYS分析振型在各階模態(tài)均能準(zhǔn)確對應(yīng)，因此可以確認(rèn)基于membrane簡化推導(dǎo)的自由邊界拋物薄殼模態(tài)振型函數(shù)能夠正確描述各階模態(tài)振動的實(shí)際形狀.

4 結(jié)論

基于membrane簡化，推導(dǎo)出一組完全滿足拋物薄殼自由邊界條件的模態(tài)振型函數(shù)，依據(jù)所推導(dǎo)的模態(tài)振型函數(shù)可得到自由邊界拋物薄殼上的理論振型.利用模態(tài)分析軟件ModalVIEW得到自由邊界拋物薄殼的低階模態(tài)實(shí)驗(yàn)振型.將理論振型、實(shí)驗(yàn)振型和ANSYS分析振型進(jìn)行對比，三者能夠較好的吻合，驗(yàn)證了所推導(dǎo)模態(tài)振型函數(shù)的正確性.

（References）

［1］Tzou Hornsen，Ding Jianghong.Actuator placement and micro-actuation efficiency of adaptive paraboloidal shells［J］.Journal of Dynamic Systems，Measurement and Control，2003，125:577－584

［2］Yue Honghao，Deng Zongquan，Tzou Hornsen.Optimal actuator locations and precision micro-control actions on free paraboloidal membrane shells［J］.Communications in Nonlinear Science and Numerical Simulation，2008，13:2298－2307

［3］羅祖道，潘紀(jì)浩.拋物旋轉(zhuǎn)扁殼的軸對稱彎曲［J］.力學(xué)學(xué)報(bào)，1963（3）:200－218

Luo Zudao，Pan Jihao.Axisymmetrical bending of shall-ow paraboloidal shell［J］.Acta Mechanica Sinica，1963（3）:200－218（in Chinese）

［4］孫博華.拋物旋轉(zhuǎn)扁殼的一般彎曲問題［J］.力學(xué)學(xué)報(bào)，1989，21（2）:245－249

Sun Bohua.The generalbending problem ofa shallow paraboloidal shell［J］.Acta Mechanica Sinica，1989，21（2）:245－249（in Chinese）

［5］Kayran A，Vinson J R，Ardic E S.A method for the calculation of natural frequencies of orthotropic axisymmetri-cally loaded shells of revolution［J］.Journal of Vibration and Acoustics，1994，116:16－25

［6］張志良，程昌鈞.旋轉(zhuǎn)薄殼轉(zhuǎn)點(diǎn)頻段的軸對稱振動解［J］.固體力學(xué)學(xué)報(bào)，2006，27（2）:630－636

Zhang Zhiliang，Cheng Changjun.Turning-point solution for axisymmetric shell vibrations［J］.Acta Mechanica Solida Sinica，2006，27（2）:630－636（in Chinese）

［7］陳星文，張若京.旋轉(zhuǎn)薄殼自由振動中3類廣義相關(guān)函數(shù)的求解［J］.計(jì)算機(jī)輔助工程，2008，17（1）:36－38

Chen Xingwen，Zhang Ruojing.Solution to three kinds of generalized related functions in free vibration of revolution shell［J］.Computer Aided Engineering，2008，17（1）:36－38（in Chinese）

［8］張志良，程昌鈞.旋轉(zhuǎn)薄殼振動耦合轉(zhuǎn)點(diǎn)問題的一致有效解［J］.中國科學(xué) G 輯，2009，39（8）:1116－1125

Zhang Zhiliang，Cheng Changjun.Uniformly valid solutions of the coupling turning-point problem in revolution shell vibration［J］.Science in China Series G，2009，39（8）:1116－1125（in Chinese）

［9］Al-Khatib O J，Buchanan G R.Free vibration of a paraboloidal shell of revolution including shear deformation and rotary inertia［J］.Thin-Walled Structures，2010，48（3）:223－232

［10］Lin Y K，Lee F A.Vibrations of thin paraboloidal shells of revolution［J］.Journal of Applied Mechanics，1960，27（4）:743－746

［11］Hoppmann W H，Cohen M I，Kunukkasseril V X.Elastic vibrations of paraboloidal shells of revolution［J］.Journal of the A-coustic Society of America，1964，36（2）:349－353

［12］Glockner P G，Tawadros K Z.Experiments on free vibration of shells of revolution［J］.Experimental Mechanics，1973，13（10）:411－421

［13］Kang J H，Leissa A W.Free vibration analysis of complete paraboloidal shells of revolution with variable thickness and solid paraboloids from a three-dimensional theory［J］.Computers and Structures，2005，83:2594－2608

［14］Yue Honghao，Deng Zongquan，Tzou Hornsen.Distributed signal analysis of free-floating paraboloidal membrane shells［J］.Journal of Sound and Vibration，2007，304:625－639

［15］曹富新.簡明彈性力學(xué)［M］.沈陽:遼寧科學(xué)技術(shù)出版社，1984:269－270

Cao Fuxin.Concise elasticity［M］.Shenyang:Liaoning Science and Technology Publishing House，1984:269－270（in Chinese）

Mode shape functions of free paraboloidal thin shells

Wang Xinjie Yue Honghao Deng Zongquan

（School of Mechatronic Engineering，Harbin Institute of Technology，Harbin 150001，China）

Based on the membrane theory，a set of mode shape functions for paraboloidal thin shell with free boundary condition were formulated.Low mode shapes of free paraboloidal thin shell were obtained by the way of modal experiment analysis.The theoretical mode shapes based on the mode shape functions were compared with the mode shapes through finite element method and experimental method.The results show that the theoretical mode shapes agree well with the mode shapes through finite element method and experimental method.Then mode shapes of paraboloidal thin shell can be correctly described by the mode shape functions.

paraboloidal thin shell;mode shape function;membrane theory;modal analysis

V 214;TU 33

A

1001-5965（2012）03-0389-05

2010-11-15;< class="emphasis_bold">網(wǎng)絡(luò)出版時(shí)間:

時(shí)間:2012-03-09 10:37

www.cnki.net/kcms/detail/11.2625.V.20120309.1037.016.html

國家自然科學(xué)基金資助項(xiàng)目（50705017，50675043）;高等學(xué)校學(xué)科創(chuàng)新引智計(jì)劃資助項(xiàng)目（B07018）

王新杰（1982－），男，河南郟縣人，博士生，wxj_hit@163.com.

（編 輯:文麗芳）