軸向運(yùn)動簡支梁振動響應(yīng)分析

齊亞峰，劉 寧，楊國來

(1.中國人民解放軍92941部隊，遼寧 葫蘆島 125000； 2.南京理工大學(xué)，南京 210094)

【機(jī)械制造與檢測技術(shù)】

軸向運(yùn)動簡支梁振動響應(yīng)分析

齊亞峰1，劉 寧2，楊國來2

(1.中國人民解放軍92941部隊，遼寧 葫蘆島 125000； 2.南京理工大學(xué)，南京 210094)

軸向運(yùn)動系統(tǒng)廣泛存在于諸如火炮發(fā)射系統(tǒng)、動力傳送帶、帶鋸、空中纜車索道等工程領(lǐng)域，理論上可將其抽象為軸向運(yùn)動弦線和梁的力學(xué)模型，具有重要理論意義和工程應(yīng)用價值?；贓uler-Bernoulli梁理論，建立了移動質(zhì)量作用下軸向運(yùn)動簡支梁振動方程，利用Galerkin法對方程進(jìn)行離散，采用Newmark-β法對方程進(jìn)行數(shù)值求解。計算了移動質(zhì)量對軸向運(yùn)動梁動態(tài)響應(yīng)的影響，結(jié)果表明：相比軸向運(yùn)動梁的自由振動，移動質(zhì)量的作用使軸向運(yùn)動梁更多的模態(tài)被激起；移動質(zhì)量速度越小，頻域成分越豐富；在移動質(zhì)量作用前半階段，梁的振動頻率較低，在后半階段振動頻率較高。

軸向運(yùn)動梁；移動質(zhì)量；橫向振動；振動響應(yīng)

由于移動質(zhì)量作用下梁的振動和軸向運(yùn)動系統(tǒng)的動態(tài)響應(yīng)在工程上具有重要的研究價值，受到研究人員的廣泛關(guān)注。移動質(zhì)量-梁系統(tǒng)常見于重要的工程結(jié)構(gòu)中，如車-橋耦合系統(tǒng)、彈-炮耦合系統(tǒng)及貨物-龍門吊耦合系統(tǒng)等；而軸向運(yùn)動系統(tǒng)在工業(yè)上也有廣泛的應(yīng)用，如空中索道、傳送帶、紙張、帶鋸及高樓升降機(jī)纜繩等，計及抗彎剛度時這類系統(tǒng)可模型化為軸向運(yùn)動梁。Michaltsos研究了載荷不同運(yùn)動速度下橋梁的動力響應(yīng)[1]；李煒明等對移動質(zhì)量在連續(xù)速度變化下的簡支梁動力響應(yīng)進(jìn)行了研究[2]；王穎澤等研究了多個移動質(zhì)量作用下身管的振動[3]；劉寧等研究了移動質(zhì)量作用下軸向運(yùn)動懸臂梁振動特性[4]；Wu基于有限元法對移動質(zhì)量作用下軸向運(yùn)動梁的動態(tài)響應(yīng)進(jìn)行了數(shù)值模擬[5]；陳力群、楊曉東等對軸向運(yùn)動梁和弦線的振動進(jìn)行了分析[6]；丁虎、陳立群推導(dǎo)了軸向運(yùn)動梁控制方程的幾種不同形式，通過數(shù)值方法分析了各方程的差異[7]?，F(xiàn)有研究很少將移動質(zhì)量與軸向運(yùn)動系統(tǒng)結(jié)合起來分析，由于生產(chǎn)生活中大部分的軸向運(yùn)動系統(tǒng)都與移動載荷相互作用，因此研究軸向運(yùn)動系統(tǒng)需要考慮移動載荷的影響。對于移動質(zhì)量作用下兩端簡支軸向運(yùn)動梁的動態(tài)響應(yīng)研究還未見有報道。

針對上述問題，本研究建立了移動質(zhì)量作用下軸向運(yùn)動簡支梁的振動方程，采用Galerkin截斷方法分離時間和空間變量，得到以模態(tài)坐標(biāo)表示的二階常微分方程組，通過Newmark-β法進(jìn)行數(shù)值求解，研究了移動質(zhì)量對梁動態(tài)響應(yīng)的影響。

1 軸向運(yùn)動梁橫向振動方程

圖1所示為移動質(zhì)量(滑塊)作用下軸向運(yùn)動簡支梁的示意圖，其中u為梁的軸向運(yùn)動速度，L為梁的長度，P為軸向力，x為軸向空間坐標(biāo)，x軸取靜止梁變形前的軸線，y(x,t)為梁的橫向位移，m為滑塊質(zhì)量，v為滑塊運(yùn)動速度，ξ為滑塊位移。本文僅考慮運(yùn)動梁的橫向振動，略去高速運(yùn)動梁軸向位移對橫向振動的影響。根據(jù)Euler-Bernoulli梁模型和Newton第二定律，移動質(zhì)量作用下軸向運(yùn)動梁的控制方程為

ρA(y,tt+2uy,xt+u2y,xx)-Py,xx+EIy,xxxx=

-m(g+y,tt+2vy,xt+v2y,xx)·δ(x-ξ)

(1)

其中：ρ為材料密度；A為梁的橫截面積；EI為梁抗彎剛度；δ為Dirac函數(shù)，逗號表示對其后面的變量求偏導(dǎo)。

圖1 移動質(zhì)量-軸向運(yùn)動梁振動力學(xué)模型

簡支梁邊界條件和初始條件為

y(0,t)=y(L,t)=0

y,xx(0,t)=y,xx(L,t)=0

(2)

y(x,0)=α(x)，y,t(x,0)=β(x)

(3)

2 數(shù)值方法

為了使連續(xù)介質(zhì)系統(tǒng)運(yùn)動方程(1)化為有限維的形式，采用Galerkin截斷法離散控制方程。取兩端鉸支靜止梁的前N階模態(tài)φi(x)作為試函數(shù)，這樣邊界條件(2)自動滿足，方便了后面的Galerkin截斷處理。方程(1)的解可表示為如下形式

φi(x)ηi(t)

(4)

φi(x)=sin(iπx) (i=1,2,…,N)

(5)

其中，ηi(t)為待定系數(shù)，式(4)寫成矩陣形式

y(x,t)=φT·η

(6)

將式(6)代入方程(1)，兩邊左乘φ(x)，并對梁全長L積分，注意試函數(shù)的正交性，最后按η整理得

(7)

式中:M為質(zhì)量矩陣；C為阻尼矩陣；K為剛度矩陣；F為載荷向量。各系數(shù)矩陣表達(dá)式為：

φφTdx+mφ(ξ)φT(ξ)

(8)

(9)

(10)

F=-mgφ(ξ)

(11)

采用Newmark-β法直接積分求解時變系數(shù)方程(7)，利用Matlab語言編寫計算程序。

3 數(shù)值結(jié)果

為了驗(yàn)證模型和數(shù)值方法的可靠性，分別以文獻(xiàn)[5]中移動質(zhì)量作用下無軸向運(yùn)動簡支梁動態(tài)響應(yīng)和文獻(xiàn)[7]中軸向運(yùn)動簡支梁自由振動為算例，計算所需參數(shù)如表1所示。文獻(xiàn)[5]是考慮移動載荷慣性力影響，基于有限元法建立具有附加系數(shù)矩陣的時變動力學(xué)模型，求解得到的結(jié)果，在本文中作為梁無軸向運(yùn)動的對比算例。文獻(xiàn)[7]采用有限差分法計算軸向運(yùn)動梁的振動響應(yīng)，在本文中作為軸向運(yùn)動梁自由振動的對比算例。算例1和算例2中梁的初始橫向速度為零，初始橫向位移為：

α1(x)=0

(12)

α2(x)=0.064x3(1-x)3

(13)

在執(zhí)行算例1時只需在模型和算法中令梁的軸向速度為零，執(zhí)行算例2時只需令移動質(zhì)量為零。

表1 算例參數(shù)

圖2為采用算例1參數(shù)計算得到的移動質(zhì)量作用下簡支梁跨中撓度響應(yīng)曲線，撓度峰值為1.189 mm，文獻(xiàn)[5]中計算的撓度峰值為1.18 mm，相對誤差為0.7%。圖3(a)為根據(jù)算例2計算得到的軸向運(yùn)動梁跨中無量綱撓度響應(yīng)曲線，振幅為0.001，圖3(b)是撓度曲線的頻譜分析結(jié)果，一階主振頻為111.8 Hz，與文獻(xiàn)[7]計算結(jié)果一致。由上述兩個算例可以看出，本文建立的模型和數(shù)值方法是可靠的。

下面在算例2的基礎(chǔ)上研究移動質(zhì)量速度對軸向運(yùn)動簡支梁振動響應(yīng)的影響，梁的軸向運(yùn)動速度為10 m/s。設(shè)滑塊的質(zhì)量為1.2 kg，考察滑塊的3種運(yùn)動速度，即分別以小于、等于和大于梁軸向速度進(jìn)行滑動。3組速度分別取為5 m/s，10 m/s和20 m/s，計算移動質(zhì)量作用期間簡支梁動態(tài)響應(yīng)，即仿真總時間取為滑塊在梁上的滑動時間，跨中撓度曲線及其FFT頻譜分析如圖4～圖6所示。

由前文算例2結(jié)果可知，無移動質(zhì)量作用時梁的前兩階模態(tài)被激起，梁作周期性有規(guī)律的自由振動，而由圖4～圖6可知，加載移動質(zhì)量后梁的振動無明顯的規(guī)律性，振動幅值減小，振動頻率比梁自由振動豐富，說明移動質(zhì)量使更多的振頻被激起。由圖4～圖6的時域響應(yīng)還可發(fā)現(xiàn)，在移動質(zhì)量作用前半階段，梁的振動頻率較低，而后半階段振動頻率較高，同時由頻域響應(yīng)可看出移動質(zhì)量速度越小，頻域成分越豐富。

圖2 算例1梁跨中撓度曲線

圖3 算例2梁跨中撓度曲線及頻譜分析

圖4 滑塊速度5 m/s時梁跨中撓度曲線及FFT分析

圖5 滑塊速度10 m/s時梁跨中撓度曲線及FFT分析

圖6 滑塊速度20 m/s時梁跨中撓度曲線及FFT分析

4 結(jié)論

目前國內(nèi)外對軸向運(yùn)動梁的動力學(xué)研究已得到了許多重要成果，但是這些成果大部分是基于梁自身的振動響應(yīng)，考慮移動質(zhì)量作用等實(shí)際工程問題較少。本文建立了移動質(zhì)量作用下軸向運(yùn)動梁的運(yùn)動方程，通過數(shù)值方法研究了移動質(zhì)量對軸向運(yùn)動梁動態(tài)響應(yīng)的影響，發(fā)現(xiàn)移動質(zhì)量的作用使梁的振幅減小，梁跨中撓度曲線頻譜分析表明，移動質(zhì)量使梁更多的頻率成分被激起，并且移動質(zhì)量運(yùn)動速度越小，梁的頻率成分越豐富。移動質(zhì)量對軸向運(yùn)動梁動態(tài)響應(yīng)的影響規(guī)律對實(shí)際工程問題中軸向運(yùn)動系統(tǒng)的振動控制、測試分析具有重要的指導(dǎo)意義。

[1] MICHALTSOS G T.Dynamic Behavior of a Single-Span Beam Subjected to Loads Moving with Variable Speeds[J].Journal of Sound and Vibration,2002,258(2):359-372.

[2] 李煒明，駱漢賓，朱宏平,等.移動質(zhì)量速度對簡支梁動力響應(yīng)的影響[J].華中科技大學(xué)學(xué)報:自然科學(xué)版，2008，36(9)：117-120.

[3] 王穎澤，張小兵.移動質(zhì)量作用下膨脹波發(fā)射裝置的振動響應(yīng)分析[J].航空動力學(xué)報，2009，24(8)：1714-1719.

[4] 劉寧，楊國來.移動質(zhì)量作用下軸向運(yùn)動懸臂梁振動特性分析[J].振動與沖擊學(xué)報，2012，31(3)：106-109.

[5] WU J J，WHITTAKER A R，CARTMELL M P.Dynamic Responses of Structures to Moving Bodies Using Combined Finite Element and Analytical Methods[J].International Journal of Mechanical Sciences 2001,43:2555-2579.

[6] CHEN Liqun,YANG Xiaodong.Nonlinear Free Transverse Vibration of an Axially Moving Beam:Comparison of Two Models[J].Journal of Sound and Vibration,2007，299(1):348-354.

[7] DING Hu，CHEN Liqun.On Two Transverse Nonlinear Models of Axially Moving Beams[J].Science in China Series E:Technological Science，2009，52(3)：743-751.

(責(zé)任編輯唐定國)

Vibration Analysis of Axially Moving Simply-Supported Beam

QI Ya-feng1, LIU Ning2, YANG Guo-lai2

(1.The No.92941stTroop of PLA, Huludao 125000, China; 2.Nanjing University of Science and Technology, Nanjing 210094, China)

Axially moving systems are widely used in engineering, such as gun launching system, power transmission band, belt saws, aerial cable tramways, etc. If the bending stiffness is concern, this kind of systems can be modeled as axially moving beams and has important theoretical significance and engineering application value. A partial-differential equation governing the transverse vibration of an axially moving simply supported beam with lumped mass was derived from Euler-Bernoulli beam model. The governing equation was resolved to a set of second order ordinary differential equations using the Galerkin method. Then the equations were calculated based on Newmark-βmethod. The effect of the moving mass on the dynamic response of the axially moving beam was studied. The results show that the lumped mass arouses more vibration modes of the simply supported moving beam than that of free vibration. Further, it is found that the decrease of the moving speed of the mass along the beam increases the frequency components of the bean vibration. In addition, the vibration frequencies of the beam increase when the lumped mass moves toward the end point. The beam vibration frequency is low, and the vibration frequency is higher in the second half of the beam.

axially moving beam; moving mass; transverse vibration; vibration response

2016-08-05；

江蘇省自然科學(xué)基金(BK20140789)；中央高?；究蒲袠I(yè)務(wù)費(fèi)專項資金(30915118826)

齊亞峰(1979—)，男，主要從事艦炮實(shí)驗(yàn)技術(shù)研究。

10.11809/scbgxb2016.12.029

齊亞峰，劉寧，楊國來.軸向運(yùn)動簡支梁振動響應(yīng)分析[J].兵器裝備工程學(xué)報，2016(12):126-129.

format:QI Ya-feng, LIU Ning, YANG Guo-lai.Vibration Analysis of Axially Moving Simply-Supported Beam[J].Journal of Ordnance Equipment Engineering,2016(12):126-129.

TJ301

A

2096-2304(2016)12-0126-04

修回日期：2016-09-10