求解變截面梁振動(dòng)特性的假設(shè)模態(tài)法

馬艷龍，李映輝

(西南交通大學(xué)力學(xué)與工程學(xué)院，成都 610031)

工程中很多結(jié)構(gòu)都可簡化為變截面梁進(jìn)行計(jì)算，如煙筒、風(fēng)力機(jī)塔架、風(fēng)機(jī)葉片等。目前針對(duì)變截面梁振動(dòng)特性的研究有很多成果［1－8］。周叮［4］對(duì)曲率半徑沿冪指數(shù)變化的一類變截面梁給出了其橫向振動(dòng)的精確解析解。王曉臣［5］用有限元法分析了截面沿梁軸向線性及拋物線性變化的矩形截面梁的振動(dòng)特性。崔燦等［6－7］采用將變截面梁分段等效的方法，給出了變截面梁振動(dòng)的半解析法和鐵木辛柯梁的快速計(jì)算方法。

本文針對(duì)自由端附有集中質(zhì)量的任意變截面梁，在使用假設(shè)模態(tài)法的基礎(chǔ)上得出新的質(zhì)量矩陣，由此提出一種該類變截面梁的頻率和振型計(jì)算方法。

1 模型介紹

考慮左端固支，右端附有一集中質(zhì)量m0的變截面梁，彎曲剛度為EI(x)，橫截面積為A(x)，梁的中心軸位移表示其撓度w(x)。變截面梁的結(jié)構(gòu)簡圖如圖1所示。

圖1 變截面梁結(jié)構(gòu)簡圖

在小變形條件下，梁的橫向振動(dòng)方程為

邊界條件為:

2 模態(tài)函數(shù)求解

由于懸臂梁為變截面，且右端有一集中質(zhì)量，不能使用均勻梁的模態(tài)函數(shù)，為此設(shè)w(x，t)=φ(x)eiωt。其中:φ(x)為待求模態(tài)函數(shù);ω 為振動(dòng)的固有頻率。代入式(1)得

設(shè)模態(tài)函數(shù)可表示為

其中:φi(x)為選取的滿足邊界條件的線性無關(guān)函數(shù);a1，a2，…，aN為待定系數(shù)。將式(4)代入式(3)可得

對(duì)式(5)兩邊同乘以 φj(x)，然后在［0，l］上積分可得

其中:aT=［a1，a2，…，aN］;M 為廣義質(zhì)量矩陣;K為廣義剛度矩陣，其元素定義為

方程(6)為關(guān)于a的線性方程組，若使其有非零解，則需滿足

由方程(9)可解得 ω =ω1，ω2，…，ωn為變截面梁的各階固有頻率。將ωi代入式(6)可求得相應(yīng)的aT=［a1，a2，…，aN］，再由式(4)求得變截面梁的模態(tài)函數(shù)。

3 實(shí)例及應(yīng)用

選取1.25MW SEC-WO2-1250型的風(fēng)力機(jī)高塔作為原型［9］。該塔高 98.12 m，塔底直徑為3.9 m，壁厚 67 mm，塔頂直徑為 2.55 mm，壁厚20 mm，塔底與塔頂之間的截面形式為均勻漸變，塔頂機(jī)艙、輪轂及槳葉總質(zhì)量為83155 kg。塔材料為鋼，密度 ρ為7900 kg/m3，彈性模量 E為210 GPa，泊松比 υ 為0.3。

計(jì)算中，選擇滿足其邊界條件的線性無關(guān)函數(shù):其中c1i和c2i為φi(x)的系數(shù)。根據(jù)邊界條件(2)可解得c1i和c2i的取值:

φi(x)的個(gè)數(shù)越多，求得高塔的模態(tài)函數(shù)φ(x)及固有頻率ωi越精確，相應(yīng)求解的復(fù)雜程度也越高，本文取式(4)中的N=8。

為檢驗(yàn)該方法的精確度，對(duì)高塔模型進(jìn)行有限元建模分析。采用有限元法、本文方法及文獻(xiàn)［8］方法所得的高塔橫向振動(dòng)的固有頻率及對(duì)比見表1。

表1 高塔模型前4階固有頻率

將固有頻率ω代入式(6)求得a，并將a代入式(4)求得高塔的各階振型。用有限元法計(jì)算該塔筒的前4階振型和同等材料、同等高度、直徑為3.2 m、壁厚為0.04 m的等截面梁的振型。前4階的振型對(duì)比如圖2所示。

圖2 2種模型的前4階振型

4 結(jié)束語

研究了自由端附有集中質(zhì)量的變截面梁的振動(dòng)特性，給出了該類變截面梁的固有頻率方程及振型的計(jì)算方法。通過工程實(shí)例的計(jì)算，說明了該方法的有效性。

［1］曹志遠(yuǎn).變截面塔式結(jié)構(gòu)的固有振動(dòng)分析［J］.振動(dòng)與沖擊，1990，36(1):40－46.

［2］Ece M C，Aydogdu M，Taskin V.Vibration of a variable cross-section beam［J］.Mechanics research communications，2007，34(2):78－84.

［3］Atmane H A，Tounsi A，Meftah S A.Free vibration behavior exponential functionally graded beams with varying cross-section［J］.Jounral of Vibration and Control，2011，17(2):311－318.

［4］周叮.一類變截面梁橫向自由振動(dòng)的精確解析解［J］.振動(dòng)與沖擊，1996，3(15):13－15.

［5］王曉臣.變截面梁有限元分析［J］.浙江工業(yè)大學(xué)學(xué)報(bào)，2012，36(3):311－315.

［6］崔燦，蔣晗，李映輝.變截面梁橫向振動(dòng)特性半解析法［J］.振動(dòng)與沖擊，2012，31(4):85－88.

［7］崔燦，李映輝.變截面鐵木辛柯梁振動(dòng)特性快速計(jì)算方法［J］.動(dòng)力學(xué)與控制學(xué)報(bào)，2012，10(2):258－262.

［8］李映輝，崔燦.含變軸力的變截面梁振動(dòng)特性計(jì)算方法［J］.重慶理工大學(xué)學(xué)報(bào):自然科學(xué)版，2012(5):28－32.

［9］張自立，陳建兵，李杰.圓球減震裝置對(duì)風(fēng)力發(fā)電高塔的振動(dòng)控制研究［J］.地震工程與工程振動(dòng)，2012，32(3):144－149.