組合L型變截面梁的自由振動特性分析

牛國華， 王剛鋒， 王 劍

(1. 長安大學(xué) 工程機(jī)械學(xué)院， 西安 710064； 2. 長安大學(xué) 道路施工技術(shù)與裝備教育部重點實驗室， 西安 710064)

作為工程結(jié)構(gòu)基本構(gòu)件，梁結(jié)構(gòu)被廣泛應(yīng)用于飛機(jī)、車輛等交通運輸產(chǎn)品[1-2]。組合梁結(jié)構(gòu)一般由多根梁組成，且在各梁節(jié)點處具有連續(xù)條件。組合梁的振動問題是工程結(jié)構(gòu)中的基本問題，得到組合梁振動的解析解對于產(chǎn)品設(shè)計及動力學(xué)分析具有重要的意義。

目前，有很多學(xué)者提出了一些組合梁結(jié)構(gòu)振動特性的計算方法。王民等[3]提出了一種組合階梯梁振動特性的計算方法，并將其方法應(yīng)用于刀柄裝夾銑刀系統(tǒng)的振動特性分析中。張文福等[4]以Timoshenko連續(xù)梁模擬索-桁架組合梁結(jié)構(gòu)，利用能量變分法分析了索-桁架組合梁結(jié)構(gòu)的固有振動。羅延忠[5]采用折算系數(shù)法建立了組合梁振動的運動微分方程，并分析了組合梁的振動特性。Wang[6]基于歐拉-伯努利梁理論，開發(fā)了一種用于求解微分方程得到固有頻率的算法，并將其運用與非對稱截面組合梁自由振動中。Georgiades 等[7-9]推導(dǎo)了考慮轉(zhuǎn)動慣量影響的結(jié)構(gòu)的線性運動方程，并將其應(yīng)用于L型梁結(jié)構(gòu)進(jìn)行非平面運動分析。以上方法在計算效率上已有很大提高，但大多數(shù)應(yīng)用于結(jié)構(gòu)簡單、截面形狀單一的組合梁中。

此外，也有很多學(xué)者提出了一些變截面梁振動特性計算方法。鮑四元等[10]基于譜幾何法和Hamilton原理，利用矩陣特征值問題研究了任意邊界條件下連續(xù)多段梁的振動特性。葛仁余等[11]基于Timoshenko梁理論，利用插值矩陣法計算軸向功能梯度變截面梁的固有頻率。王劍等[12]利用導(dǎo)納綜合法與改進(jìn)的傳遞矩陣法，計算了存在質(zhì)量偏心的階梯梁模型的彎縱耦合振動特性。馬一江等[13]基于傳遞矩陣法，提出了一種計算含多條裂紋變截面簡支梁固有頻率的方法。劉鵬等[14]基于歐拉-伯努利梁理論，利用樣條有限點法計算了截面高寬度沿軸線性變化的變截面梁振動特性。Li[15]提出了一種新的精確方法，并將其應(yīng)用于具有任意數(shù)量裂紋和集中質(zhì)量的階梯梁進(jìn)行自由振動分析。Lee[16]利用基于歐拉-伯努利梁理論和Timoshenko梁理論的Chebyshev-tau方法，對階梯梁進(jìn)行自由振動分析。Failla 等[17]提出了一種求解非連續(xù)歐拉梁振動特性計算方法，并將其應(yīng)用于具有內(nèi)部平移和旋轉(zhuǎn)彈簧的階梯梁。Mao[18]提出使用ADM分解方法來研究階梯梁自由振動，并將其應(yīng)用到求解非均勻梁的振動問題。上述文獻(xiàn)針對變截面梁振動特性提出了多種計算方法，但還存在計算效率低、精度差及推導(dǎo)過程繁瑣等問題。另外，針對組合型變截面梁的研究還相對較少。

本文針對L型變截面梁的橫向自由振動問題，推導(dǎo)了其固有頻率及振型函數(shù)的求解過程并給出了相應(yīng)算例。根據(jù)相鄰梁段的力平衡條件和位移連續(xù)條件，建立整個組合梁的特征矩陣，得到解析形式的頻率方程。通過求解頻率方程即可得到L型變截面梁橫向振動的固有頻率，進(jìn)一步求得振型函數(shù)。通過給出具體算例并與有限元結(jié)果對比，驗證了本文方法的有效性。

1 L型梁自由振動的計算方法

1.1 梁的橫向振動微分方程

基于歐拉-伯努利梁理論，各段梁長度為Li，水平梁段數(shù)為m，總段數(shù)為n的梁，自由振動的運動微分方程為

i=1,2,…，n

(1)

式中：y(x,t)為梁的橫向振動撓度函數(shù)；Ei為第i段梁的彈性模量；Ii為第i段梁的截面慣性矩；ρiAi為第i段梁的的線密度。L型變截面梁模型如圖1所示。

圖1 L型變截面梁模型

由于梁每段的EiIi和ρiAi為常數(shù)，式(1)可表示為

i=1,2,…，n

(2)

使用分離變量法，式(2)的通解為

y(x,t)=φ(x)T(t)

(3)

其中

T(t)=e-iωt

(4)

式中，ω為橫向振動的固有圓頻率。

φ(x)為梁橫向振動的模態(tài)振型函數(shù)，可表示為

φ(x)=C1sinβx+C2cosβx+C3sinhβx+

C4coshβx

(5)

將第i段梁的模態(tài)振型函數(shù)寫成

(6)

1.2 邊界條件

L型變截面梁兩端的邊界條件分別為一端固定一端自由。固定端位移和轉(zhuǎn)角為零，可得到方程組

(7)

自由端彎矩和剪力為零，可得到方程組

(8)

式中：y′,y″,y′″表示對位移求一階、二階、三階偏導(dǎo)。

1.3 平衡和連續(xù)條件

L型梁的變截面處存在力平衡和位移連續(xù)條件，對于橫梁及豎梁，相鄰兩段梁連續(xù)處左側(cè)截面與右側(cè)截面的位移、轉(zhuǎn)角、彎矩、剪力均相等，可以得到矩陣

(9)

對于連接橫梁與豎梁處的直角連接點，平衡及連續(xù)條件如圖2所示。

圖2 直角連接點處的平衡及連續(xù)條件

由圖2所示的直角連接點處的連續(xù)條件，可得位移及轉(zhuǎn)角關(guān)系如下

(10)

由圖2所示的直角連接點處的平衡條件，可得彎矩及剪力關(guān)系如下

(11)

綜上可得到直角連接點處的平衡及連續(xù)條件如下

(12)

1.4 固有頻率及振型函數(shù)的求解

結(jié)合邊界條件和連續(xù)條件，構(gòu)造特征矩陣H

(13)

矩陣H的各元素表達(dá)式如下推導(dǎo)。

由邊界條件式(7)、(8)可得

h11=h13=h22=h24=1

h3,4n-3=-cosβnLn，h3,4n-2=-sinβnLn

h3,4n-1=coshβnLn，h3,4n=sinhβnLn

由平衡及連續(xù)條件式(9)可得

h4i+1,4i-3=cosβiLi，h4i+1,4i-2=sinβiLi

h4i+1,4i-1=coshβiLi，h4i+1,4i=sinhβiLi

h4i+2,4i-3=-βisinβiLi，h4i+2,4i-2=βicosβiLi

h4i+2,4i-1=βisinhβiLi，h4i+2,4i=βicoshβiLi

h4i+1,4i+1=h4i+1,4i+3=-1，h4i+2,4i+2=h4i+2,4i+4=-βi+1

其中，i=1,2,…,m-1,m+1,m+2,…,n-1。

由直角連接處的平衡及連續(xù)條件式(12)可得

h4m+2,4m-3=-βmsinβmLm，h4m+2,4m-2=βmcosβmLm

h4m+2,4m-1=βmsinhβmLm，h4m+2,4m=βicoshβmLm

h4m+1,4m+1=h4m+1,4m+3=1，h4m+2,4m+2=h4m+2,4m+4=-βm+1

其余元素皆為0。

特征矩陣H滿足

Hψ=0

(14)

其中

(15)

為位移系數(shù)向量。為了求解系統(tǒng)的固有頻率，上式需有非零解，即矩陣H的行列式必須為零。通過令特征矩陣H的特征值為零，可以求解出固有圓頻率ω。將固有圓頻率ω代回特征矩陣并求解方程組，可以得到第1段～第n段的振型函數(shù)的待定系數(shù)，最終得到各段梁的模態(tài)振型函數(shù)。

2 算例分析與驗證

2.1 梁的主要參數(shù)

將表1中組合L型變截面梁的參數(shù)代入第一章所提出的公式及特征矩陣中，定義掃描初始頻率及終止頻率，使用MATLAB軟件掃描出組合梁的前10階固有頻率，再將固有頻率代回特征矩陣求出各段梁的振型系數(shù)，最后繪制出振型曲線，結(jié)果見下文。

表1 組合L型變截面梁的主要參數(shù)

2.2 收斂性分析

(1) 本文方法在建立特征矩陣時，隨著組合梁總段數(shù)n的增加，矩陣維數(shù)隨之增大且大多數(shù)元素為0，故針對此問題有必要做收斂性分析。

分別取n為10、16及20，即水平梁及豎直梁各5、8及10段時，對組合梁振動特性做收斂性分析。當(dāng)n取10時，組合梁結(jié)構(gòu)的物理參數(shù)如表1所示；當(dāng)n取16、20時，組合梁結(jié)構(gòu)的物理參數(shù)為水平梁為變截面(相鄰兩段相互交替，交替規(guī)律與n取10時相同)，豎直梁為等截面，各段梁長度、密度、及彈性模量均與表1相同。將3種組合梁的參數(shù)采用上述求解方法，得到前10階固有頻率，與有限元法相比各階頻率誤差如圖3所示。

圖3 本文方法與有限元法相比各階頻率誤差

由圖3可知，當(dāng)n取10和16時，各階頻率與有限元法相比誤差均小于1%，但當(dāng)n為20時，各階頻率與有限元法相比誤差有的已超出1%。因此為了確保計算精度，本文方法適用于n<16的組合梁模型中。

(2) 本文方法是基于局部坐標(biāo)下建立的平衡和連續(xù)方程，針對該方法下的高階模態(tài)有必要作收斂性分析，結(jié)果如圖4所示。

圖4 收斂性分析結(jié)果

由圖4可以看出，固有圓頻率在48階之后出現(xiàn)發(fā)散。由相關(guān)文獻(xiàn)[19-20]分析可知，在車輛航天等機(jī)械結(jié)構(gòu)中，低階模態(tài)對結(jié)構(gòu)影響較為重要。因此本文方法在計算組合L型梁結(jié)構(gòu)振動特性時仍具有較高的準(zhǔn)確性和收斂性，是滿足實際工程需要的。

2.3 與有限元法的對比

使用表1中組合梁的參數(shù)，利用ANSYS對L型變截面梁建模及有限元分析。有限元計算中，采用beam188單元模擬。有限元分析結(jié)果與本文方法計算結(jié)果如表2所示。

表2 本文方法與有限元法固有頻率對比

從表2中可以看出，本文方法與有限元法固有頻率相比吻合良好，各階誤差均在0.063%以內(nèi)，從而驗證了本文方法的有效性。

圖5給出了組合梁的前九階振型，可以看出本文結(jié)果與有限元結(jié)果振型曲線完全吻合，從而驗證了該方法的準(zhǔn)確性。

(a) 一階振型

3 結(jié) 論

本文研究了組合L型變截面梁的自由振動特性，主要工作及結(jié)論如下：

(1) 基于歐拉-伯努利梁理論，提出了一種考慮截面變化的L型梁自由振動特性計算方法，通過分析組合梁相鄰梁段之間的平衡及連續(xù)條件，構(gòu)建了整個組合梁的特征矩陣，通過求解特征矩陣得到了組合梁的固有頻率及振型曲線。

(2) 通過與有限元計算結(jié)果對比，驗證了本文所提出方法的有效性及準(zhǔn)確性；并對計算結(jié)果進(jìn)行了收斂性及有效性分析，結(jié)果表明本文方法適用于n<16的組合梁模型。

(3) 由于特征矩陣組裝方便的特點，使得本文所提出方法在處理幾何和物理參數(shù)變化的組合梁，如漸變截面L型梁的振動特性問題時具有較大的優(yōu)勢，進(jìn)而可以為組合梁結(jié)構(gòu)的設(shè)計、分析及優(yōu)化等問題提供理論基礎(chǔ)。

附錄A

(A.1)

(A.2)

(A.3)

(A.4)

其中，i=1,2,…,m-1,m+1,m+2,…,n-1

(A.5)

(A.6)

(A.7)

其余各元素均為0。

角標(biāo)含義: 例如H1:4,1:4中，前面的1:4表示第1行～第4行，后面的1:4表示第1列～第4列，其余同上。