部分入水懸臂梁的自由彎曲振動固有特性研究

劉少文，林澤峰，楊翊仁

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

懸臂梁的振動是一類重要的力學(xué)模型。由于實際需要，許多真實物體是放置在水中的，例如橋墩，它們可以被認為是懸臂梁置于水中，對其振動分析時，水的作用是不能忽視的。懸臂梁在水中的固有特性研究是一個具有理論和實際意義的課題，它屬于流固耦合動力學(xué)分析領(lǐng)域。在現(xiàn)有研究中，許多情況下考慮水深與梁高相等，然后用無水情況下的梁自由振動的振型函數(shù)作為基函數(shù)，將有水振型函數(shù)按無水振型函數(shù)作級數(shù)展開［1－2］，或者是考慮部分有水情況時直接將兩段方程分別列出來，求出兩段的帶系數(shù)的精確解，然后根據(jù)梁的邊界條件和兩段對接條件列出以待定系數(shù)為未知量的線性方程組進行求解［3］。本文對不同水深情況下梁振動的固有頻率的變化情況進行了探討，并使用Galerkin離散法和以ANSYS軟件為平臺的有限元法兩種方法進行計算和結(jié)果對比，給出了圓截面梁的前兩階固有頻率隨水深的變化規(guī)律［4－15］。

1 振動方程的建立

懸臂梁坐標(biāo)示意圖如圖1所示?？紤]懸臂梁豎直放置，其長度為H，入水部分為0～h，懸臂梁橫截面為半徑為a的圓，則由振動力學(xué)理論可知梁在空氣中的橫向自由振動的振動方程為:

圖1 懸臂梁坐標(biāo)示意圖

梁作簡諧振動，則令

代入式(1)可得

其中:EI為梁的彎曲剛度;ρ為梁的密度;A為橫截面積。

對于入水懸臂梁部分，本文進行如下假設(shè):水為不可壓縮、無黏性的理想流體，并且忽略水表面自由波的影響。則流體運動的速度勢ψ滿足Laplace方程，其在柱坐標(biāo)中的形式為

使用分離變量法，令

代入方程(4)，并根據(jù)相應(yīng)邊界條件可得［3］:

其中:h為水的深度;ρ0為水的密度。因此入水部分梁的彎曲自由振動方程為

自由振動時，梁振型隨時間變化為諧振規(guī)律，消去時間項則可將方程(8)變?yōu)槿缦滦问?

綜合方程(1)和方程(9)，可以得到該懸臂梁的橫向振動方程為:

2 振動方程的求解

采用Galerkin法離散方程(10)。根據(jù)Galerkin法的原理，若系統(tǒng)振動方程為L(X)－λM(X)=0，L，則 M為2p次微分算子。選取比較函數(shù)X1，X2，…，Xn，令，則應(yīng)該滿足

本文選取:

其中:

式(12)存在非零解的條件為

解此廣義特征值問題，變化水深參數(shù)h就可以得到部分入水懸臂梁的前兩階固有頻率隨水深的變化規(guī)律。

3 算例

取a=0.02 m，H=1 m，ρ=2 400 kg/m3，ρ0=1 000 kg/m3，E=2.44 GPa，水深變化可以分別取h=0，0.2，0.4，0.6，0.8，1 m，計算得到部分入水懸臂梁的前兩階固有頻率見表1。

為了對比計算結(jié)果的精度，本文還使用了有限元方法分析該系統(tǒng)，分析過程在ANSYS軟件平臺上進行。

本算例中，懸臂梁使用solid45單元模擬，水使用fluid30單元模擬，水和梁的幾何模型均為圓柱體。為了使分析結(jié)果更加精確，水模型的半徑設(shè)置為梁半徑的10倍。將梁單元和水單元接觸的地方使用FSI方法進行耦合，施加相應(yīng)的約束條件后變化水深參數(shù)，分別進行模態(tài)分析就能得到梁的前兩階固有頻率隨水深的變化規(guī)律，如表2所示。

表1 采用Galerkin法計算的不同水深下的前兩階固有頻率

表2 采用ANSYS計算的不同水深下的前兩階固有頻率

將表1和表2中的相應(yīng)數(shù)據(jù)進行擬合，分別得到兩種方法計算的1階固有頻率和2階固有頻率隨水深的變化情況曲線，見圖2和3。

圖2 不同水深1階固有頻率對比曲線

圖3 不同水深2階固有頻率對比曲線

從圖2和3中可以看出:懸臂梁的前兩階固有頻率隨水深的增加而下降。為了更加清晰地得到下降的具體情況，將不同水深下的彎曲梁固有頻率同無水情況下固有頻率作比值，令比值為c，可以得到比值c隨水深的變化情況，如圖4和5所示。

圖4 不同水深的一階固有頻率比值對比曲線

圖5 不同水深的二階固有頻率比值對比曲線

4 結(jié)論

通過計算結(jié)果可知:隨著水深的增加，圓截面懸臂梁彎曲振動的前兩階固有頻率逐漸減小。滿水時，采用Galerkin法計算的1階固有頻率減為無水頻率的85.91%，采用ANSYS軟件計算的則為無水頻率的84.40%。同樣，2階固有頻率減為無水頻率的 85.97%(Galerkin法計算結(jié)果)和84.71%(ANSYS軟件計算結(jié)果)。計算結(jié)果表明:水對于懸臂梁振動的影響相當(dāng)于分布附加質(zhì)量的作用，附加質(zhì)量隨著水深增加而增大，附加質(zhì)量的影響降低了梁的固有振動頻率。

［1］居初榮，曾心傳.彈性固體與液體的耦聯(lián)振動分析［M］.北京:地震出版社，1983.

［2］楊萬理.深水橋梁動水壓力分析方法研究［D］.成都:西南交通大學(xué)，2013.

［3］張悉德.部分埋入水中懸臂圓柱體的彎曲自由振動［J］.應(yīng)用數(shù)學(xué)和力學(xué)，1982(4):537－546.

［4］周叮.兩側(cè)受液時懸臂梁的自由振動分析［J］.工程力學(xué)，1991(3):107－115.

［5］張敏.橋墩與河水流固耦合振動分析［D］.大連:大連交通大學(xué)，2006.

［6］楊吉新，黨慧慧，雷凡.考慮水作用的橋墩自振特性計算方法對比分析［J］.世界地震工程，2010，26(3):91－95.

［7］陳貴清，楊翊仁.板狀梁結(jié)構(gòu)的流致振動研究［J］.學(xué)術(shù)動態(tài)，2003(2):31－33.

［8］劉振宇.深水橋梁的地震響應(yīng)分析［D］.成都:西南交通大學(xué)，2007.

［9］蘇海東.分析流固耦振的新方法及其應(yīng)用研究［D］.武漢:華中科技大學(xué)，2007.

［10］楊笑冬.基于ANSYS的懸臂梁模態(tài)分析［J］.機電一體化，2008(6):58－60.

［11］張阿漫，戴紹仕.流固耦合動力學(xué)［M］.北京:國防工業(yè)出版社，2011.

［12］Sharan S K.Hydrodynamic loadings due to the motion of large offshore structures［J］.Computers ＆ Structures，1989，32(6):132－137.

［13］戴大農(nóng)，王勖成，杜慶華.流固耦合系統(tǒng)動力響應(yīng)的模態(tài)分析理論［J］.固體力學(xué)學(xué)報，1990，11(4).

［14］劉雪萊.基于流固耦合的冷卻風(fēng)扇葉片模態(tài)及疲勞分析［J］.重慶理工大學(xué)學(xué)報:自然科學(xué)版，2014(9):25－30.

［15］田彥偉，崔曉娜.三維拉普拉斯方程的求解［J］.安陽師范學(xué)院學(xué)報，2007，10(5):3－5.