帶有彈簧-質量-圓柱殼耦合結構自由振動分析

胡　浩， 李正良， 于　偉

(重慶大學 土木工程學院，重慶　400045)

?

帶有彈簧-質量-圓柱殼耦合結構自由振動分析

胡浩， 李正良， 于偉

(重慶大學 土木工程學院，重慶400045)

摘要：采用子結構導納法研究了簡支邊界條件下帶有多根彈簧-集中質量-圓柱殼耦合結構的自由振動。根據(jù)已有結果，通過求解多根彈簧等效剛度的推導思路，采用柔度法得到各子結構矩陣元素的柔度系數(shù)，進而求得集中質量帶有多根彈簧的圓柱殼耦合結構自振頻率及模態(tài)理論公式，并與已有文獻結果作對比，證明了理論推導的正確性。應用理論方法進一步求解了集中質量帶有三根彈簧圓柱殼耦合結構的自振頻率及振型，并與建立的ANSYS有限元模型分析結果對比，二者結果的誤差可忽略不計，再次驗證所提出理論的合理性及正確性。

關鍵詞：彈簧-集中質量-圓柱殼；導納法；自由振動；等效剛度

板殼耦合結構廣泛應用于航天器、導彈、壓力容器、建筑和潛艇等結構中。對于復雜結構，通常分解為多個子結構進行分析研究，現(xiàn)階段對板殼耦合結構振動的研究大多采用子結構導納法[1-4]，即從柔度概念出發(fā)研究耦合結構振動。趙芝梅等[5-6]在傳統(tǒng)導納法的基礎上提出了改進的子結構線導納法，同時將其用于鋪板與圓柱殼耦合結構的振動響應分析，并進一步研究了多點激勵下加筋板殼耦合結構的振動特性。

另一類常見的耦合結構即為彈簧-板殼結構和彈簧-附加質量-板殼結構，該類結構也得到了廣泛的研究。Allaei等[7]采用點導納法研究了帶有徑向彈簧的圓環(huán)(輪胎)的振動特性，但與圓環(huán)相連的彈簧的另一端都假設成固定。Azimi[8]應用模態(tài)擴展法及導納法分析了任意徑向位置處具有彈性支撐和剛性支撐的圓板的振動特征。Amabili[9]應用導納法研究了帶有集中質量空腹及充液圓柱殼的振動，并與試驗結果對比，吻合良好，證明了該理論的正確性。Achong[10]利用導納法分析了帶有附加質量的圓板及扁殼，把該理論方法應用于聲道鋼盤鼓，并與其他方法及試驗所得結果對比，吻合較好。Amabili等[11]基于Flügge薄殼理論利用Rayleigh-Ritz法研究了具有附加質量及彈性約束的鉚接圓柱殼振動特征。在此基礎上，Amabili等[12]進一步分析了集中質量旋轉慣性對空腹及充液圓柱殼體頻率的影響，理論及試驗結果表明：集中質量旋轉慣性除了對殼體的局部模態(tài)起重要影響外，一般情況下可以忽略。

對附帶質量彈簧的殼體的研究主要集中于質量附著在殼體上，或者質量通過一根或者兩根共線彈簧連接在殼體上[13]，該種連接比較簡單。本文在已有研究成果的基礎上進一步分析了集中質量與圓柱殼體通過多根彈簧連接的耦合結構的振動特性，推導了該彈簧-質量-殼體耦合結構的自振頻率及位移響應公式，并與已有結果及ANSYS有限元結果對比，證明了本文理論推導的有效性及正確性。

1理論分析

1.1子結構柔度系數(shù)推導

對附加質量-彈簧-圓柱殼耦合結構建立坐標系見圖1，A為圓柱殼;L為殼體長度;R為半徑;h為厚度;ρ為材料密度;ν為泊松比;E為彈性模量；集中質量B懸浮在圓柱殼中心，M為質量;不考慮質量自重的影響，其通過L根彈簧與殼體連接，集中質量及各彈簧均處于同一高度x0處，每根彈簧與水平方向夾角為θi，其剛度系數(shù)為ki。

基于導納法，即一般的柔度法理論[13-14]，把圓柱殼和集中質量當作兩個子系統(tǒng)A和B,各個系統(tǒng)受力情況見圖2，假設殼體與彈簧連接處殼體受到的外力為FAi，殼體在該處位移為ZAi，對應彈簧的另一端與集中質量連接處質量受到的外力為FBi，質量在該處位移為ZBi。易得如下關系式：

圖1　圓柱殼結構模型Fig.1 Geometry of a cylindrical shell

圖2　受力模型Fig.2 Model of forces

(1)

(2)

由式(1)易得式(3)和式(4)：

(3)

(4)

由式(2)、式(3)和式(4)易得式(5)：

(5)

(6)

式中：αij及βij為柔度系數(shù);具體表達式分別見式(10)及式(17);i,j取為1,2,…,L。

令系數(shù)矩陣行列式為0，可求得耦合結構頻率值：

(7)

對于兩端垂向自由簡支約束圓柱殼，軸向、切向及徑向陣型函數(shù)分別假設如下：

(8)

式中：系數(shù)Amn、Bmn、Cmn為振動幅值;m，n分別為軸向半波數(shù)及周向波數(shù);φ為模態(tài)相位角。

由文獻[13]可知單點諧波激勵下，圓柱殼位移：

(9)

圓柱殼柔度系數(shù)：

(10)

(11)

集中質量柔度系數(shù)：

(12)

式中：ωmn及ω分別為圓柱殼的自振頻率和耦合結構的自振頻率。

借鑒文獻[15]對煤氣柜中活塞彈簧等效剛度的推導思路，考慮彈簧的徑向變形，本文從柔度法角度出發(fā)，在已知某根彈簧變形量的前提下，推導求得任意一根彈簧壓縮量或者伸長量，即可知質量上柔度系數(shù)βij。

圖3為彈簧變形值計算簡圖。忽略集中質量點大小，把質量當作一個質點，質量中心點水平向右移動單位位移，同時水平彈簧壓縮量也是單位值，與該根水平彈簧夾角為φ值的彈簧變形量為dφ，斜彈簧的初始位置為C2坐標，變形后移動到C3坐標處，對應處的斜彈簧位置由C4移到C5坐標，易得如下幾何關系：

圖3　彈簧變形值計算簡圖Fig.3 Model of spring deflection

(13)

(14)

y3-y2=1

(15)

根據(jù)式(13)～式(15)關系，得彈簧變形量為：

(16)

所以，已知某根彈簧的變形量，根據(jù)式(16)即可求得與此彈簧成任意夾角的彈簧變形量。結合柔度的概念，考慮彈簧的拉壓方向性，易得集中質量柔度系數(shù)βij如下：

(17)

式中：θi、θj為第i，j根彈簧與水平方向的夾角，取值范圍是0～2π。

1.2 圓柱殼自振頻率求解

根據(jù)Flügge薄殼理論[16]，基于式(8)假設的振型函數(shù)，可得兩端簡支圓柱殼頻率方程：

(18)

(19)

(20)

由式(18)求得圓柱殼頻率ωmn，易得振型幅值比例關系：

(21)

故總位移表達式可由式(9)按多個單點諧波激勵下的位移響應線性疊加得到，如下：

(22)

其中Nmn參見式(11)。

1.3 理論推導的正確性驗證

由式(5)和式(6)，假設質量與圓柱殼體通過一根彈簧相連，則

(23)

當彈簧剛度無窮大，即k1→∞時，可認為質量塊直接附著在殼體上；當彈簧剛度很小，即k1→0時，可認為彈簧-質量塊對結構振動無影響，即ω=ωmn。式(23)與文獻[13]中推導公式相同，說明文獻[13]中為本文推導的特例，證明了本文方法的正確性。

2有限元模型建立及其數(shù)值分析

2.1建立有限元模型

為進一步應用本文方法，建立了集中質量帶有三根彈簧的兩端簡支圓柱殼耦合結構有限元模型，彈簧均勻分布，相互間夾角為120°，模型參數(shù)為：M=200 g，k1=k2=k3=4*105N/m，x0=262 mm，E=206 GPa，ρ=7 700 kg/m3，ν=0.3，L=664 mm，R=175 mm，h=1.0 mm。ANSYS模型中不考慮質量塊的重力，只考慮平面內慣性力，彈簧采用COMBIN14單元模擬一維軸向自由度，殼體采用SHELL63單元模擬。三維ANSYS有限元模型見圖4。耦合結構振型圖見圖5。

圖4　三維ANSYS有限元模型Fig.4 Finite element model

圖5　耦合結構振型圖Fig.5 Modal of coupled cylindrical shell

2.2計算結果與有限元結果對比

計算耦合結構頻率及模態(tài)，此時彈簧根數(shù)L=3，通過式(18)，利用Rayleigh-Ritz法計算得到單獨簡支圓柱殼32階模態(tài)頻率值ωmn(m=1,2,3,4;n=1,2,…,8)，再由軟件Mathematica求解式(7)，得到耦合結構頻率值ω，最后通過式(22)求解耦合結構的頻率及模態(tài)，見圖5。

本文主要列出了m=1時的前6階彎曲模態(tài)作為比較分析，圖5(a)～圖5(f)為耦合圓柱殼結構的振型圖，圖5(a2)～圖5(f2)為彈簧-質量高度處殼體徑向位移，圖中頻率值為本文計算結果，括號內頻率值為有限元分析結果。通過頻率值對比可見，本文計算結果與有限元結果很接近，誤差很小，對于工程而言，誤差可以忽略，再次為本文方法的合理性提供了依據(jù)。

3結論

本文基于已有文獻的研究成果，進一步推導了多根彈簧-質量-圓柱殼耦合結構自振模態(tài)及頻率值的理論公式，與已有文獻作對比，證明了本文理論推導的正確性，并建立有限元模型算例與本文理論值對比，二者所得結果誤差很小，再次表明本文方法的合理性。

參 考 文 獻

[ 1 ] Wilken I D, Soedel W. The receptance method applied to ring-stiffened cylindricalshells: analysis of modal characteristics[J]. Journal of Sound and Vibration,1976,44(4):563-576.

[ 2 ] Huang D T, Soedel W. Study of the forced vibration of shell-plate combinations using the reacceptance method[J].Journal of Sound and Vibration, 1993, 162(2): 341-369.

[ 3 ] Huang D T, Soedel W. Natural frequencies and modes of a circular plate welded to a circular cylindrical shell at arbitrary axial positions[J]. Journal of Sound and Vibration,1993, 162(3):403-427.

[ 4 ] Huang D T. Influences of small curvatures on the modal characteristics of the joined hermetic shell structures[J]. Journal of Sound and Vibration,2000, 238(1): 85-111.

[ 5 ] 趙芝梅，王敏慶. 板殼結構振動功率流的子結構線導納法研究[J]. 工程力 學, 2012,29(8): 297-300.

ZHAO Zhi-mei, WANG Min-qing. A structure line receptance method for analysisof power flow in a clylindrical shell with a floor partition[J]. Journal of Engineering Mechanics, 2012,29(8):297-300.

[ 6 ] 趙芝梅，盛美萍，楊 陽. 多點激勵下加筋板殼耦合結構的振動特性[J]. 工程力學, 2013,30(11): 239-244.

ZHAO Zhi-mei, SHENG Mei-ping, YANG Yang. Vibtation of a stiffened shell with a ribbed plate under multi-point excitations[J]. Journal of Engineering Mechanics,2013,30(11):239-244.

[ 7 ] Allaei D, Soedel W,Yang T Y. Eigenvalues of rings with radial spring attachments[J].Journal of Sound andVibration,1987,121(3):547-561.

[ 8 ] Azimi S. Free vibration of circular plates with elastic or rigid interior support[J]. Journal of Sound andVibration,1988,120(1):37-52.

[ 9 ] Amabili M. Free vibration of a fluid-filled circular cylindrical shell with lumped masses attached, using the receptance method[J]. Journal of Shock and Vibration,1996,3(3):159-167.

[10] Achong A. Vibrational analysis of mass loaded plates and shallow shells by the receptance method with application to the steelpan[J]. Journal of Sound and Vibration, 1996,191(2): 207-217.

[11] Amabili M, Garziera R. Vibrations of circular cylindrical shells with nonuniform constraints, elastic bed and added mass; Part I: empty and fluid-filled shells[J]. Journal of Fluids and Structures,2000,14: 669-690.

[12] Amabili M, Garziera R, Carra S. The effect of rotary inertia of added masses on vibrations of empty and fluid-filled circular cylindrical shells[J]. Journal of Fluids and Structures,2005,21: 449-458.

[13] Soedel W. Vibrations of shells and plates[M]. Marcel Dekker, New York,1993.

[14] Shyh-Chin Huang, Kao-An Lin. A general receptance theory for vibrations of combined structures[C]//Twelfth International Congress on Sound and Vibration,11-14,July,2005,Lisbon.

[15] 賈冬云，曹平周，姜德進. 活塞對大型儲氣柜柜體抗震性能的影響分析[J]. 合肥工業(yè)大學學報,2011,34(1):109-113.

JIA Dong-yun,CAO Ping-zhou, JIANG De-jin.Study of damping effect of piston on seismic performance of large-scale gasholder[J]. Journal of HeFei University of Technology, 2011,34(1):109-113.

[16] Leissa A W. Vibrations of shells[M]. NASA SP-288,Washington, DC: Government Printing Office,1973.

Free vibration analysis for a spring-mass-cylindrical shell coupled structure

HUHao,LIZheng-liang,YUWei(College of Civil Engineering, Chongqing University, Chongqing 400045, China)

Abstract:The free vibration of a spring-mass-cylindrical shell coupled structure under simply supported boundary condition was investigated with the substructure receptance method. Through deriving the equivalent stiffness of several springs, the compliance coefficients of matrix elements of sub-structures were derived with the compliance method. Theoretical formulas for the natural frequencies and modes of the coupled structure via the receptance method, they were compared with those in other reference to verify the validity of the theoretical derivation here. The numerical results for natural frequencies and modal shapes of a three springs-lumped mass-sylinder shell coupled structure obtained with the proposed method were compared with those using the finite element software ANSYS to verify the reasonability and correctness of the theory presented here.

Key words:springs-lumped mass-cylindrical shell; receptance method; free vibration; equivalent stiffness

中圖分類號:TB123

文獻標志碼：A

DOI:10.13465/j.cnki.jvs.2016.07.032

收稿日期：2014-12-03修改稿收到日期：2015-04-20

基金項目：國家自然科學基金項目(51278511)

第一作者 胡浩 男，博士生，1986年生