葛仁余,　張金輪,　韓有民,　牛忠榮,　程長征,　楊智勇

(1.安徽工程大學(xué) 力學(xué)重點(diǎn)實(shí)驗(yàn)室,安徽 蕪湖　241000; 2.合肥工業(yè)大學(xué) 土木與水利工程學(xué)院,安徽 合肥　230009)

?

非均質(zhì)變截面桿軸向振動的插值矩陣法分析

葛仁余1,張金輪1,韓有民1,牛忠榮2,程長征2,楊智勇2

(1.安徽工程大學(xué) 力學(xué)重點(diǎn)實(shí)驗(yàn)室,安徽 蕪湖241000; 2.合肥工業(yè)大學(xué) 土木與水利工程學(xué)院,安徽 合肥230009)

文章分析了非均質(zhì)變截面彈性桿軸向自由振動問題,該類桿件的彈性模量、橫截面面積和單位桿長質(zhì)量是截面位置的函數(shù)，對于函數(shù)的具體形式無任何限制條件。運(yùn)用插值矩陣法可一次性地計(jì)算出各階振動固有頻率,它們具有同階精度,并可同時獲取相應(yīng)的振型函數(shù)。數(shù)值計(jì)算結(jié)果與現(xiàn)有結(jié)果對比表明了文中方法的有效性和良好的計(jì)算精度。

變截面桿;軸向振動;固有頻率;插值矩陣法

變截面桿是工程結(jié)構(gòu)中應(yīng)用廣泛的承力構(gòu)件,對于變截面桿振動問題,已有研究成果尚缺乏一種計(jì)算簡單且精度較高的方法。分析結(jié)構(gòu)振動的固有頻率問題,李茲法是常用的近似計(jì)算方法,但計(jì)算精度取決于事先假設(shè)的振型函數(shù),對于高階頻率的計(jì)算往往不能達(dá)到滿意的效果。文獻(xiàn)[1-3]用不同的方法對均質(zhì)變截面桿進(jìn)行了大量的研究，這些解法不同程度地存在計(jì)算方法復(fù)雜、求解精度不好的缺點(diǎn)；文獻(xiàn)[4]建立了均質(zhì)變截面直桿縱振動的頻率方程,將小波分析理論與微分求積法相結(jié)合求解微分方程,對其振動頻率進(jìn)行了分析；文獻(xiàn)[5]根據(jù)均質(zhì)變截面直桿的縱振與等截面直梁的彎曲問題的一些相似性質(zhì),把求變截面直桿縱振函數(shù)變?yōu)榍罅簭澢鷷r的截面上的彎矩問題,對應(yīng)于桿的橫截面面積的某些變化規(guī)律,得到了求桿的縱振函數(shù)及固有頻率的計(jì)算公式。

隨著航空航天工業(yè)的飛速發(fā)展,非均質(zhì)變截面直桿振動問題也引起了人們的普遍重視。在工程結(jié)構(gòu)的優(yōu)化設(shè)計(jì)中,為了達(dá)到強(qiáng)度和質(zhì)量的最佳分布或滿足結(jié)構(gòu)功能的要求,常常使用非均質(zhì)變截面直桿類構(gòu)件。但非均質(zhì)變截面彈性直桿的軸向自由振動的控制方程是變系數(shù)二階偏微分方程,現(xiàn)有的解析方法只能得到一些特殊類型的非均質(zhì)變截面直桿的振動方程的精確解,例如楔形體和錐體、截去尖端的楔形體和錐體、彎曲剛度隨長度坐標(biāo)以任意冪函數(shù)變化的懸臂直桿等。文獻(xiàn)[6]獲得了計(jì)算微弱楔形直桿固有頻率的簡便公式,但該公式對截面變化很大的楔形直桿不能使用；文獻(xiàn)[7]獲得了彈性模量、密度和截面面積以指數(shù)形式變化時的直桿縱振和扭振的一個通解,但求解固有頻率時計(jì)算量特別大且解的形式比較復(fù)雜；文獻(xiàn)[8]采用攝動法得到了該類直桿縱振的固有頻率和振型,但也只適用于彈性模量、密度和截面面積變化微弱的直桿；文獻(xiàn)[9]借助積分方程理論得出一種計(jì)算非均質(zhì)變截面直桿的振動頻率和振型的方法。

本文從非均質(zhì)變截面桿件一般形式的振動微分方程出發(fā),將非均質(zhì)變截面桿件振動頻率的計(jì)算轉(zhuǎn)化為常微分方程組的特征值求解,運(yùn)用插值矩陣法求解常微分方程組特征值問題,可獲得桿件軸向振動所有頻率,同時相應(yīng)的振型一并解出。插值矩陣法原理簡單、易于計(jì)算、實(shí)用性好、計(jì)算精度高,能靈活適用于任意變截面彈性直桿的軸向自由振動問題。

1　非均質(zhì)變截面桿軸向振動微分方程

1.1振動微分方程的建立

一變截面桿件如圖1所示。

圖1　變截面桿

設(shè)桿件固有軸向振動位移為w(x,t),則非均質(zhì)變截面桿軸向振動的平衡微分方程[3]為：

(1)

將(1)式展開,則有：

(2)

當(dāng)非均質(zhì)變截面桿做固有振動時,設(shè)其解為:

(3)

其中,W(x)為桿件軸向振動的振型函數(shù);ω為桿件軸向振動的固有頻率。將(3)式代入(2)式可得非均質(zhì)變截面桿軸向振動的微分方程為:

(4)

其中,E(x)、A(x)、ρ(x)分別為變截面桿彈性模量、截面面積和密度,皆為截面位置x的函數(shù)。

設(shè)ξ=x/l,則有0≤x≤l,ξ∈[0,1],單位化后(4)式轉(zhuǎn)化為:

(5)

由于(5)式中含有ω2項(xiàng),引入新變量g(ξ),消去(5)式中含ω2的非線性項(xiàng),組成新的方程組如下:

(6)

當(dāng)材料為均質(zhì)材料時,E(ξ)=E,ρ(ξ)=ρ,由(6)式可得均質(zhì)材料變截面桿軸向振動的方程為:

(7)

從而,非均質(zhì)和均質(zhì)變截面桿軸向振動固有頻率分析轉(zhuǎn)化為求解常微分方程組(6)式、(7)式的常微分方程組特征值問題。

1.2軸向振動微分方程的邊界條件

將(3)式代入下列3種情形,于是將桿端邊界條件轉(zhuǎn)化為以下表達(dá)形式。

情形1一端固定、一端自由的桿件軸向振動邊界條件為:

(8)

情形2兩端自由的變截面桿件,其軸向振動邊界條件為:

(9)

情形3大頭固定、小頭有彈性支承(設(shè)其彈性剛度為k)變截面桿件,其軸向振動邊界條件為:

(10)

至此,非均質(zhì)變截面桿件軸向振動頻率ω的計(jì)算歸結(jié)為求解在邊界條件(8)～(10) 式下線性常微分方程組(6) 式的特征值問題。一般需要數(shù)值方法求解這類問題,本文采用插值矩陣法計(jì)算[10-12]。

2　常微分方程組特征值問題求解

針對一般的常微分方程組兩點(diǎn)邊值問題,A.Φ.斯米爾諾夫用拉格朗日多項(xiàng)式插值創(chuàng)立了積分矩陣法,求解區(qū)間最多分為7段。文獻(xiàn)[11]在積分矩陣法基礎(chǔ)上建立了插值矩陣法,采用了分段多項(xiàng)式插值,使得求解區(qū)間可以任意等分,計(jì)算值可達(dá)到充分高的精度。假設(shè)一個含r個方程的線性常微分方程組特征值問題為(11)式,相應(yīng)的邊值條件為(12)式,即

(11)

(12)

插值矩陣法將求解區(qū)間[xa,xb]剖分為n+1個節(jié)點(diǎn)、n個子區(qū)間,在每個子區(qū)間上使用低階分段多項(xiàng)式函數(shù)插值逼近待求函數(shù),以常微分方程組中出現(xiàn)的最高階導(dǎo)數(shù)在離散節(jié)點(diǎn)上的值作為離散系統(tǒng)的未知參數(shù),形成代數(shù)特征值方程組,解之可以獲得特征值與相應(yīng)的特征函數(shù)。該方法可以同時求解出常微分方程組特征值問題中出現(xiàn)的所有函數(shù)和其各階導(dǎo)數(shù),且具有同階精度。本文采用插值矩陣法來求解非均質(zhì)變截面桿件軸向振動頻率的常微分方程組特征值問題。

3　算例分析

算例1設(shè)一直桿的彈性模量、橫截面面積、密度均以指數(shù)函數(shù)變化,即

(13)

其中,α、β、γ為常數(shù);E0、A0、μ0分別為ξ=0處的彈性模量、橫截面面積和單位桿長的質(zhì)量。將(13)式代入(6)式得:

(14)

(15)

令

(16)

將(16)式代入(15)式可得:

(17)

引入新函數(shù)g(ξ)=λW(ξ),則(17)式變換為：

(18)

由插值矩陣法求解常微分方程組(18)式,可以獲得非均質(zhì)變截面桿在一端固定另一端自由和兩端固定的邊界條件下微分方程的特征值λ,再由(16)式可以獲得軸向振動的固有頻率如下：

(19)

本文解與文獻(xiàn)[2]解、文獻(xiàn)[13]解、Bessel解、精確解的對比見表1、表2所列。由表1、表2可見,當(dāng)m=0,ε=1/2時,本文計(jì)算結(jié)果與精確解吻合良好,插值矩陣法計(jì)算精度比文獻(xiàn)[2]解和文獻(xiàn)[13]解要高。同時本文方法還給出了m、ε為其他任意值時的計(jì)算結(jié)果,而m、ε為其他任意值時獲得精確解是十分困難的,所以這也是本文方法的一個顯著優(yōu)點(diǎn)。

表1　算例1一端固定、另一端自由變截面桿軸向振動前8階固有頻率　Hz

表2　算例1兩端固定變截面桿軸向振動前8階固有頻率　Hz

算例2設(shè)非均質(zhì)變截面彈性直桿的參數(shù)為:

(20)

其中,E0、A0、μ0分別為ξ=0處的彈性模量、橫截面積和單位桿長質(zhì)量。將(20)式代入(6)式得：

(21)

(22)

引入新函數(shù)g(ξ)=λW(ξ),則(22)式變換為：

(23)

因此,非均質(zhì)變截面彈性直桿軸向振動固有頻率及其相應(yīng)振型的計(jì)算,轉(zhuǎn)化為由插值矩陣法求解常微分方程組(23)式特征值問題。由本文方法計(jì)算出一端固定、另一端自由的非均質(zhì)變截面桿軸向振動前8階固有頻率,并與文獻(xiàn)[2]和文獻(xiàn)[14]計(jì)算結(jié)果比較,見表3所列。由表3可見,本文方法計(jì)算值與已有結(jié)果吻合良好,再次證明了本文方法的有效性。由本文方法計(jì)算兩端固定非均質(zhì)變截面桿軸向振動前8階固有頻率,結(jié)果見表4所列。

表3　算例2一端固定、另一端自由變截面桿軸向振動前8階固有頻率　Hz

表4　算例2兩端固定變截面桿軸向振動前8階固有頻率　Hz

本文方法可計(jì)算出(3)式中非均質(zhì)變截面桿軸向振動的前若干階固有頻率ωn,同時相應(yīng)的振型W(ξ)也一并解出。表3中前4階固有頻率所對應(yīng)的振型函數(shù)W(ξ)的分布曲線如圖2所示。

圖2　變截面桿軸向振動前4階固有頻率對應(yīng)的振型

4　結(jié)　　論

本文從非均質(zhì)變截面桿件一般形式的振動方程出發(fā),采用變量代換法,將關(guān)于非均質(zhì)變截面直桿軸向振動固有頻率的一組非線性常微分方程的特征值問題轉(zhuǎn)化為線性特征問題,用插值矩陣法求解,獲得了軸向振動固有頻率及相應(yīng)的振型。本文方法可一次性地計(jì)算出非均質(zhì)變截面直桿軸向振動所有固有頻率和相應(yīng)的振型函數(shù),克服了其他方法的一些缺點(diǎn),如計(jì)算方法復(fù)雜、求解精度不好等；求出的同一階振型函數(shù)及其各階導(dǎo)數(shù)的計(jì)算精度是同階的,在利用振型的一階導(dǎo)數(shù)計(jì)算振動結(jié)構(gòu)應(yīng)力場時,這是一個顯著優(yōu)點(diǎn)；本文方法數(shù)值計(jì)算結(jié)果與現(xiàn)有結(jié)果吻合良好,表明了本文方法的有效性和良好的計(jì)算精度；本文方法適用于任意變截面彈性直桿的軸向自由振動問題,在計(jì)算復(fù)雜工程結(jié)構(gòu)振動的固有頻率和振型方面具有較高的計(jì)算精度和一定的工程應(yīng)用價值。

[1]侯祥林,范煒,賈連光.變截面壓桿臨界荷載的迭代算法[J].哈爾濱工業(yè)大學(xué)學(xué)報(bào),2011,43(增刊1): 237-240.

[2]黎明安,王忠民.非均質(zhì)變截面彈性直桿縱向自由振動的差分解法[J].力學(xué)與實(shí)踐,1992,14(2): 32-35.

[3]張瑞平.可展為冪級數(shù)的變截面彈性直桿的縱向自由振動分析[J].機(jī)械科學(xué)與技術(shù),2000,25 (增刊1): 61-62.

[4]張菊梅,趙鳳群,黨曉敏.變截面彈性直桿縱振動分析的小波-DQ法[J].力學(xué)與實(shí)踐,2010,32(4): 71-73.

[5]袁鎰吾,李志堅(jiān).變截面直桿縱向自由振動的一種解[J].工程力學(xué),1996,13(3):105-113.

[6]劉先志.推算楔形直桿縱振與扭振自主頻率的一個方法[J].中國科學(xué),1977,6(2): 536-546.

[7]劉先志.關(guān)于非均質(zhì)變截面彈性直桿的縱向自主振動[J].應(yīng)用數(shù)學(xué)和力學(xué),1980,1(2):237-245.

[8]周叮.非均質(zhì)變截面彈性直桿縱向自由振動的漸近解法[J].力學(xué)與實(shí)踐,1990,12(4):16-19.

[9]許明田,周學(xué)圣,程德林.一種研究非均質(zhì)變截面彈性直桿振動的新方法[J].山東工業(yè)大學(xué)學(xué)報(bào),1992,22(2):7-12.

[10]葛仁余,程長征,牛忠榮,等.各向同性材料切口尖端熱流密度奇性特征分析[J].合肥工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,37(3):281-285.

[11]牛忠榮.兩點(diǎn)邊值問題的插值矩陣法及誤差分析[J].合肥工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),1993,16(1):25-32.

[12]楊智勇,牛忠榮,葛仁余,等.功能梯度疊層厚板彎曲半解析求解[J].合肥工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2016,39(3):355-359.

[13]吳曉,姚春梅,張龍庭.再談“非均質(zhì)變截面彈性直桿的縱向自主振動”[J].陜西工學(xué)院學(xué)報(bào),1995,11(1): 26-29.

[14]龐家駒.機(jī)械振動習(xí)題集[M].北京:清華大學(xué)出版社,1982:235-236.

(責(zé)任編輯張淑艷)

Longitudinal vibration analysis of nonuniform elastic straight bar with variable cross-section by interpolating matrix method

GE Renyu1,ZHANG Jinlun1,HAN Youmin1,NIU Zhongrong2,CHENG Changzheng2,YANG Zhiyong2

(1.Key Laboratory for Mechanics, Anhui Polytechnic University, Wuhu 241000, China; 2.School of Civil and Hydraulic Engineering, Hefei University of Technology, Hefei 230009, China)

This paper deals with the longitudinal vibration problems of functions of section location of elastic straight bar including elasticity modulus, section area and mass per unit length under the normal conditions. There is no restriction to the specific form of function. By applying the interpolating matrix method to solving the established equations, all the natural frequencies of longitudinal vibration companying with the corresponding vibration mode functions of bar with variable cross-section are calculated at a time. All the calculated natural frequencies of longitudinal vibration have the same high accuracy. The numerical results show that the computed results from the presented method have very high accuracy compared with the existing solutions.

variable cross-section bar; longitudinal vibration; natural frequency; interpolating matrix method

2015-03-30；

2016-06-18

國家自然科學(xué)基金資助項(xiàng)目(11272111;11372094);安徽省高校自然科學(xué)研究重點(diǎn)資助項(xiàng)目(KJ2016A055)

葛仁余(1969-),男,安徽合肥人,博士,安徽工程大學(xué)副教授;

牛忠榮(1957-),男,安徽合肥人,博士,合肥工業(yè)大學(xué)教授,博士生導(dǎo)師;

10.3969/j.issn.1003-5060.2016.08.016

O326

A

1003-5060(2016)08-1084-05

程長征(1979-),男,安徽太湖人,博士,合肥工業(yè)大學(xué)教授,博士生導(dǎo)師.