軸向移動懸臂梁高效動力學建模及頻率響應(yīng)分析1)

華洪良 廖振強 張相炎

(南京理工大學機械工程學院，南京210094)

軸向移動懸臂梁高效動力學建模及頻率響應(yīng)分析1)

華洪良 廖振強2)張相炎

(南京理工大學機械工程學院，南京210094)

軸向移動梁動力學問題具有廣泛的工程應(yīng)用背景，如：機械手、機床主軸、武器身管等.計算軸向移動梁動力學響應(yīng)是評估結(jié)構(gòu)動力學性能以及最終指導(dǎo)結(jié)構(gòu)設(shè)計的一個重要手段.采用Rayleigh-Ritz法、拉格朗日方程推導(dǎo)了軸向移動懸臂梁時變動力學方程.選取冪級數(shù)函數(shù)構(gòu)造試函數(shù)對軸向移動系統(tǒng)動力問題進行求解.冪級數(shù)函數(shù)良好的積分與微分性能，使得推導(dǎo)容易以矩陣的形式快速進行，便于符號運算軟件直接生成MATLAB程序.由于MATLAB基本數(shù)據(jù)單位為矩陣，符號軟件生成的程序只需經(jīng)過簡單修改便可進行動力學計算.大大縮短了軸向移動梁從建模到動力學分析的時間，過程十分高效.通過四組算例，將本文方法計算得到的動力學響應(yīng)與文獻數(shù)據(jù)進行對比，對該方法準確性進行了驗證，并給出了冪級數(shù)函數(shù)擬合階數(shù)的選取原則.以此為基礎(chǔ)，研究了軸向移動梁的頻率響應(yīng)特性.分為考慮重力與忽略重力兩種情況，討論了軸向振動幅度對其頻率響應(yīng)特性的影響.

計算動力學，冪級數(shù)基函數(shù)，頻響分析,軸向移動梁

引言

軸向移動動力系統(tǒng)具有廣泛的工程應(yīng)用背景，常見的有：機械手、機床主軸、武器身管[1-3]、傳送帶[4-6]、衛(wèi)星結(jié)構(gòu)等.近年來，軸向移動系統(tǒng)動力學是一研究熱點.由于動力響應(yīng)在評估結(jié)構(gòu)動力學性能以及最終結(jié)構(gòu)設(shè)計中不可或缺，因此，一種快速、準確的建模與響應(yīng)分析方法在軸向移動系統(tǒng)動力學研究過程中至關(guān)重要.

在軸向移動系統(tǒng)動力學研究中，通常將軸向移動系統(tǒng)簡化為軸向移動弦或軸向移動梁.根據(jù)已有文獻，目前軸向移動系統(tǒng)動力學研究主要集中于如下幾類問題：軸向移動彈性梁與旋轉(zhuǎn)剛性支撐的耦合動力學[7-9]，軸向移動黏彈性梁 (或弦)[10-14]、外流體軸向移動梁[19-22]、含有摩擦邊界的軸向移動梁[23-24]、軸向移動梁的參數(shù)振動[25-26]、軸向移動梁橫向振動抑制技術(shù)研究[27-28]、軸向與橫向耦合振動研究[29-31]等.

在軸向移動系統(tǒng)動力學建模過程中，主要有兩種方法.第一種為分析力學方法[7-10,29-30,32-35]，并且應(yīng)用最為廣泛的分析力學方法為第二類拉格朗日方程.在該方法中，首先采用能量法對系統(tǒng)進行能量分析，得到相應(yīng)的拉格朗日函數(shù)，經(jīng)離散后采用拉格朗日方程推導(dǎo)得到其運動方程.第二種為彈性力學方法[12,14-24,26-28,31,34,37].采用該方法時，首先對研究對象取分離體進行受力分析，得到其控制方程，經(jīng)離散后直接得到其運動方程.

在離散過程中，兩種主要的離散技術(shù)為有限單元法[8-9,32-33,38]與假設(shè)模態(tài)法[34-37].在假設(shè)模態(tài)法中，需要一系列滿足所有(或者部分)邊界條件的基函數(shù)來構(gòu)造一個試函數(shù)對原函數(shù)進行逼近.在軸向移動梁控制方程離散過程中，通常選取其振型函數(shù)來構(gòu)造試函數(shù).梁的振型函數(shù)由一系列三角函數(shù)與雙曲函數(shù)構(gòu)成，形式復(fù)雜，其復(fù)雜程度取決于其具體邊界條件.在常見的邊界條件中，只有簡支邊界梁的振型函數(shù)相對簡單，而固支邊界、自由邊界、懸臂邊界等梁的振型函數(shù)較復(fù)雜.在建模過程中，通常需要采用振型函數(shù)之間的正交性對方程進行解耦.然而，由于梁的軸向運動，在其動力學建模過程中會產(chǎn)生諸多不同于定常結(jié)構(gòu)的陀螺項，利用目前已知的模態(tài)函數(shù)正交性條件尚無法使這些陀螺項完全解耦[39]，因此軸向移動梁動力學建模尚無統(tǒng)一方法.

基于上述因素，本文嘗試采用冪級數(shù)作為基函數(shù)構(gòu)造試函數(shù)對軸向移動系統(tǒng)動力問題進行求解.由于冪級數(shù)基函數(shù)具有較好的積分與微分性能，更重要的是對任何兩個基函數(shù)之間的正交性并不敏感.這一特點，使得積分演算變得非常迅速，并且對于不同問題均很容易以矩陣的形式進行推導(dǎo)，推導(dǎo)過程更為簡潔.采用矩陣形式的最主要優(yōu)勢在于：能夠采用符號運算軟件推導(dǎo)并直接生成MATLAB程序.由于MATLAB基本數(shù)據(jù)單位為矩陣，符號軟件生成的程序只需經(jīng)過簡單修改便可進行動力學計算.這樣大大縮短了軸向移動梁從建模到動力學分析的時間，提高了效率.

本文主要工作：

(1)基于歐拉梁理論，首先推導(dǎo)了軸向移動梁動能與勢能函數(shù)，考慮時變軸向移動速度.隨后，采用Rayleigh-Ritz法與冪級數(shù)基函數(shù)對系統(tǒng)動能與勢能函數(shù)進行離散，并通過拉格朗日方程推導(dǎo)得到系統(tǒng)運動方程.

(2)將本文計算得到的動力學響應(yīng)與相關(guān)文獻數(shù)據(jù)進行比較，對本文方法與模型的準確性進行驗證.并給出冪級數(shù)函數(shù)擬合階數(shù)的選取原則.

(3)研究軸向往復(fù)運動頻率對橫向振動的影響，即頻率響應(yīng)分析(身管武器應(yīng)用背景).分為考慮重力與忽略重力兩種情況進行討論.

1 動力學建模

考慮一軸向移動懸臂梁模型，如圖1.梁總長LB，懸臂長L(t).楊氏模量、慣性矩、密度、截面積分別為E,I,ρ,A.梁橫向位移采用w(x(t),t)進行描述.

梁上任意一點發(fā)生變形后，其位移矢量為

圖1 軸向移動懸臂梁Fig.1 Axially moving cantilever beam

將式(1)對時間求導(dǎo)得到速度矢量

式中，點和撇分別表示對時間、空間求偏導(dǎo).根據(jù)式(2)，系統(tǒng)動能為

梁彎曲變形產(chǎn)生的應(yīng)變能為

考慮梁軸向運動存在加速度，其軸向慣性力產(chǎn)生的拉伸應(yīng)變能為

根據(jù)式(3)～式(5)，得到系統(tǒng)拉格朗日函數(shù)

2 離散與運動方程

首先構(gòu)造試函數(shù)

對梁的橫向位移進行離散.式中

在式(8)中

為滿足懸臂梁邊界條件的冪級數(shù)函數(shù).式中，n1為冪級數(shù)函數(shù)階數(shù)；αji為冪級數(shù)函數(shù)系數(shù)，確定αji的一個簡易方法就是對懸臂梁振型函數(shù)進行擬合.在擬合前，首先將懸臂梁振型函數(shù)

映射至 1-1空間.式中，cosλi+chλi/2shλisinλi為一系列幅值映射系數(shù)；λi為特征方程cosλichλi+1=0的根.經(jīng)過該映射后，梁長度與自由端橫向位移均被無量綱化為1.將梁長度映射為1后，有利于后續(xù)縮放為梁實際懸臂長度進行動力學計算.將梁自由端橫向位移映射為1后，結(jié)果后處理將變得十分簡單，梁自由端真實的橫向位移即為各基函數(shù)坐標的和.

式(10)中冪級數(shù)系數(shù)可根據(jù)式(11)擬合得到.對于自適應(yīng)任意懸臂長度的基函數(shù)可以將式(10)沿X軸縮放得到

將式(7)代入式(6)，并將結(jié)果代入拉格朗日方程，把重力作為廣義力考慮，得到系統(tǒng)運動方程

值得一提的是，式(13)中C為梁彎曲變形與軸向運動產(chǎn)生的陀螺項，而非結(jié)構(gòu)阻尼.并且，Vx=與Ax=分別為梁軸向運動速度與加速度.Vx=V0+Axt，V0為初始速度.由于式(17)為一系列冪級數(shù)函數(shù)乘積的積分運算，很容易采用符號運算軟件進行推導(dǎo)得到其代數(shù)格式，并自動生成MATLAB程序.

3 模型驗證

本節(jié)采用4個算例對本文方法準確性進行驗證.前兩個算例源于文獻[36]，在該算例中，梁作軸向勻速運動.后兩個算例源于文獻[29]，算例中梁作軸向勻加速運動.4個算例中，梁截面積、楊氏模型、慣性矩、密度均一致，分別為1.4661×10?3m2,6.8335N/m2,1.1073×10?8m4,2738.6kg/m3.其余參數(shù)設(shè)置見表1.

計算得到的梁自由端動力學響應(yīng)與文獻[36]和文獻[29]對比見圖2和圖3.本文計算結(jié)果與文獻結(jié)果吻合較好，說明本文采取的方法以及建立的模型準確可靠.

圖2 勻速收縮、伸展運動Fig.2 The retracting and extruding motion with constant velocity

圖3 勻加速收縮、伸展運動Fig.3 The retracting and extruding motion with constant acceleration

表1 驗證算例物理參數(shù)Table 1 Physical parameters of the veri fi cation examples

4 冪級數(shù)基函數(shù)擬合階數(shù)的影響

第2節(jié)中，冪級數(shù)函數(shù)擬合階數(shù)n1的大小將直接決定冪級數(shù)基函數(shù)的擬合精度，并最終影響到動力學計算精度.本節(jié)主要討論冪級數(shù)擬合階數(shù)對動力學計算精度的影響.

根據(jù)式(11)，采用MATLAB軟件編程，在計算域內(nèi)對各振型函數(shù)均勻采樣，本文采用500個樣本點.然后根據(jù)得到的樣本點以及給定的擬合階數(shù)，采用poly fi t函數(shù)便可擬合得到冪級數(shù)函數(shù)系數(shù)αji.為了給出擬合階數(shù)的選取原則，采用冪級數(shù)函數(shù)的最大相對擬合誤差作為量化指標

需要說明的是，筆者在研究過程中發(fā)現(xiàn)在多數(shù)情況下，細長懸臂梁結(jié)構(gòu)動力響應(yīng)主要取決于其一階模態(tài).因此，本節(jié)只討論第一階模態(tài)擬合階數(shù)對動力學計算精度的影響.通過計算，在表2中給出了擬合階數(shù)為1至10范圍內(nèi)的最大相對擬合誤差.對比發(fā)現(xiàn)，隨著擬合階數(shù)的增大，冪級數(shù)函數(shù)的最大相對擬合誤差逐漸減小.

表2 最大相對擬合誤差Table 2 The maximum relative fi tting error

圖4 冪級數(shù)函數(shù)擬合階數(shù)對動力學計算精度的影響Fig.4 The e ff ect of fi tting order of the power series function on computational accuracy

隨后，以第3節(jié)算例2為例，分別采用不同擬合階數(shù)的冪級數(shù)函數(shù)進行動力學計算，得到梁自由端位移響應(yīng)曲線對比如圖4.以10階擬合精度的冪級數(shù)函數(shù)計算得到的動力學響應(yīng)曲線為參考進行對比研究.圖4(a)表明，當擬合階數(shù)在1至4階范圍內(nèi)時，動力學計算產(chǎn)生了較大的誤差.尤其是當擬合階數(shù)為1時，計算得到的響應(yīng)曲線根本無法反應(yīng)梁的真實振動情況.這是因為采用一階冪級數(shù)函數(shù)時，式(17)中N00對空間X求導(dǎo)兩次后為零，相當于結(jié)構(gòu)剛度消失.當擬合階數(shù)增加至5階及其以上后，響應(yīng)曲線與10階響應(yīng)曲線基本吻合，差異很小.因此，為了保證動力學計算精度，應(yīng)采用至少5階冪級數(shù)函數(shù)進行擬合.由表2可知，5階冪級數(shù)函數(shù)對應(yīng)的最大相對擬合誤差為0.0034%.因此，筆者建議在選取冪級數(shù)函數(shù)擬合階數(shù)時，應(yīng)將其最大相對擬合誤差控制在千分之一以內(nèi).

5 頻響分析

本節(jié)討論軸向振動頻率對梁橫向振動的影響.其工程應(yīng)用背景為武器浮動身管.很多身管武器在設(shè)計過程中，為了減小后坐力，經(jīng)常將其身管設(shè)計為浮動形式，使得身管在發(fā)射過程中能夠沿著槍(炮)身做前后相對運動，吸收一部分后坐動能.武器連續(xù)發(fā)射過程中，身管不停地做軸向往復(fù)運動，使得膛口橫向振動受到擾動.由于身管的射擊精度直接于彈頭出膛口瞬時膛口的振動狀態(tài).因此，研究此類振動系統(tǒng)的頻響特性對匹配身管武器系統(tǒng)動力學特性、避免結(jié)構(gòu)共振具有重要意義.

在研究頻響特性時，使軸向移動梁做軸向簡諧運動，其懸臂長度、軸向速度、加速度由下式給出

式中，L0表示梁平均懸臂長度；Amp和ω分別表示簡諧運動幅值與頻率.算例所用參數(shù)設(shè)置如表3.

表3 算例參數(shù)Table 3 Example parameters

5.1 考慮重力影響

圖5給出了不同軸向振動幅值下的梁自由端頻率響應(yīng)特性曲線.從圖5可以發(fā)現(xiàn)，軸向振動幅值越大，幅頻響應(yīng)曲線中特殊頻率點越多.當軸向運動幅值為10mm時(圖5(a))，幅頻響應(yīng)曲線中只有c和d兩個特殊頻率點.當軸向運動幅值增加至100mm時，幅頻響應(yīng)曲線中共出現(xiàn)了5個特殊頻率響應(yīng)點，如圖5(d)，各頻率值見表4.

圖5 軸向移動梁頻率響應(yīng)特性(考慮重力)Fig.5 Frequency response characteristic of the axially moving beam(considering gravity)

表4 各頻率響應(yīng)點頻率值(Hz)Table 4 Frequency of each response point(Hz)

當梁軸向振動頻率由頻率點c向頻率點d靠近過程中，梁橫向振動幅度逐漸減小，并在頻率點d達到最小值.因此，在往復(fù)軸向運動系統(tǒng)動力學設(shè)計過程中，匹配其軸向振動頻率使其盡可能接近d頻率點，能夠有效減小橫向振動，改善系統(tǒng)動力性能.

然而，需要注意的是d頻率點會因梁軸向振動幅度而變得不穩(wěn)定，見圖5(b)～圖5(d).隨著梁軸向振動幅度的增加，在頻率曲線中d頻率點附近會出現(xiàn)e頻率點，該頻率點為不穩(wěn)定頻率點.當梁在e頻率點做軸向振動時，梁橫向振動幅度迅速增加并失穩(wěn)，如圖5(d).并且，e頻率點頻帶寬度會隨著梁橫向振動幅度的增大而增大，對梁的橫向振動十分不利.避免e頻率點出現(xiàn)的一個有效方法就是調(diào)整梁軸向振動幅度.

筆者發(fā)現(xiàn)，在重力作用下，各頻率點頻率值分布有一定規(guī)律性：梁平均懸臂長度對應(yīng)的第一階模態(tài)頻率為14.07Hz，a,b,e頻率點頻率值分別分布在梁第一階模態(tài)頻率的1/3,1/2,2倍處附近.因此，在軸向運動系統(tǒng)動力學設(shè)計過程中，不僅要使軸向振動頻率避開一階模態(tài)頻率.在軸向振動幅值較大的情況下，避開一階模態(tài)頻率的1/3,1/2,2倍頻率值也是必要的.

5.2 忽略重力影響

將重力加速度g設(shè)為0，則可忽略重力對梁橫向振動的影響.同樣分析不同軸向振動幅值下的梁頻率響應(yīng)特性，結(jié)果見圖6.

對比圖5與圖6可以發(fā)現(xiàn)，忽略重力影響后，梁的頻率響應(yīng)特性變得相對簡單.當梁軸向振動幅度較小時，軸向運動頻率對梁的橫向振動影響很小，如圖6(a).隨著梁軸向振動幅度的增大，頻響曲線在28.36Hz附近出現(xiàn)共振峰(即5.1節(jié)中e頻率點)，使系統(tǒng)失穩(wěn).梁軸向振動幅度越大，e頻率點共振帶寬越大，這一現(xiàn)象與文獻[33]得出的結(jié)論一致.

關(guān)于軸向運動梁的穩(wěn)定性，主要是由于式(15)中的C為一陀螺項，其符號受到梁運動方向的影響，即：梁做伸展運動時，C符號為正，具有阻尼效應(yīng)，能夠耗散系統(tǒng)能量使系統(tǒng)獲得穩(wěn)定.當梁做收縮運動時，C符號為負，收縮運動過程中梁會從固定邊界處吸收能量，使系統(tǒng)處于不穩(wěn)定運動狀態(tài)，甚至失穩(wěn).關(guān)于軸向運動梁的穩(wěn)定性的詳細討論，可參考文獻[32].

圖6 軸向移動梁頻率響應(yīng)特性(忽略重力)Fig.6 Frequency response characteristic of the axially moving Beam(neglecting gravity)

6 結(jié)論

本文采用Rayleigh-Ritz法、拉格朗日方程推導(dǎo)了軸向移動懸臂梁時變動力學方程.與以往離散方法不同的是，本文選取冪級數(shù)函數(shù)構(gòu)造試函數(shù)對軸向移動系統(tǒng)動力問題進行求解.由于冪級數(shù)基函數(shù)具有較好的積分與微分性能，這一特點使積分運算變得非常迅速，并且很容易以矩陣的形式進行推導(dǎo)，并可以采用符號運算軟件直接生成MATLAB程序.由于MATLAB基本數(shù)據(jù)單位為矩陣，符號軟件生成的程序只需經(jīng)過簡單修改便可進行動力學計算.這樣大大縮短了軸向移動梁從建模到動力學分析的時間，過程十分高效.相關(guān)算例對比表明，本文采取的方法與建立的模型準確可靠.

通過建立的時變動力學方程對梁的軸向共振問題進行了研究，得到如下結(jié)論：

(1)在選取冪級數(shù)函數(shù)擬合階數(shù)時，應(yīng)將其最大相對擬合誤差控制在千分之一以內(nèi)，以保證較好的動力學計算精度.

(2)軸向振動幅值越大，幅頻響應(yīng)曲線中特殊頻率點越多，并且容易導(dǎo)致系統(tǒng)失穩(wěn).同時，失穩(wěn)頻率寬度越大.

(3)考慮重力影響時，會在梁一階模態(tài)頻率值的1/3、1/2、2倍處附近產(chǎn)生共振峰.共振峰峰值隨著軸向振動幅度的增加而增加.忽略重力影響，即梁在水平面內(nèi)的振動，梁的頻率響應(yīng)特性相對簡單，只在梁一階模態(tài)頻率值的2倍處附近出現(xiàn)系統(tǒng)失穩(wěn).

1 Gimm HI,Cha KU,Cho CK.Characterizations of gun barrel vibrations of during fi ring based on shock response analysis and shorttime Fourier transform.Journal of Mechanical Science and Technology,2012,26(5):1463-1470

2 Hua HL,Liao ZQ,Song JE.Vibration reduction and fi ring accuracy improvement by natural frequency optimization of a machine gun system.Journal of Mechanical Science and Technology,2015,29(9):3635-3643

3 劉寧,楊國來.移動質(zhì)量作用下軸向運動懸臂梁振動特性分析.振動與沖擊,2012,31(3):102-105(Liu Ning,Yang Guolai.Vibration property analysis of an axially moving cantilever beam with the effect of a moving mass.Journal of Vibration and Shock,2012,31(3):102-105(in Chinese))

4 SuwekenG,VanHorssenWT.Onthetransversalvibrationsofaconveyor belt with a low and time-varying velocity.Part II:The beamlike case.Journal of Sound and Vibration,2003,267(5):1007-1027

5 Pakdemirli M,Oz HR.In fi nite mode analysis and truncation to resonant modes of axially accelerated beam vibrations.Journal of Sound and Vibration,2008,311(3-5):1052-1074

6 Ponomareva SV,van Horssen WT.On the transversal vibrations of an axially moving continuum with a time-varying velocity:Transient from string to beam behavior.Journal of Sound and Vibration,2009,325(4):959-973

7 Al-Bedoor BO,Khulief YA.Vibrational motion of an elastic beam with prismatic and revolute joints.Journal of Sound and Vibration,1996,190(2):195-206

8 Al-Bedoor BO,Khulief YA.Finite element dynamic modeling of a translating and rotating fl exible link.Computer Methods in Applied Mechanics&Engineering,1996,131(1):173-189

9 Al-Bedoor BO,Khulief YA.General planar dynamics of a sliding flexible link.Journal of Sound&Vibration,1997,206(5):641-661

10 Lee U,Oh H.Dynamics of an axially moving viscoelastic beam subject to axial tension.International Journal of Solids and Structures,2005,42(8):2381-2398

11 丁虎,陳立群.軸向運動黏彈性梁平面耦合非線性受迫振動.上海大學學報(自然科學版),2009,15(6):649-652,658(Ding Hu,Chen Liqun.Coupled forced response of nonlinear axially moving viscoelastic beam.Journal of Shanghai University(Natural Science Edition),2009,15(6):649-652,658(in Chinese))

12 丁虎,胡超榮,陳立群等.軸向變速黏彈性Rayleigh梁非線性參數(shù)振動穩(wěn)態(tài)響應(yīng).振動與沖擊,2012,31(5):135-138(Ding Hu,Hu Chaorong,Chen Liqun,et al.Steady state response of nonlinear vibration of an axially accelerating viscoelastic Rayleigh beam.Journal of Vibration and Shock,2012,31(5):135-138(in Chinese))

13 Marynowski K.Dynamic analysis of an axially moving sandwich beam with viscoelastic core.Composite Structures,2012,94(9):2931-2936

14 唐有綺.軸向變速黏彈性Timoshenko梁的非線性振動.力學學報,2013,45(6):965-973(Tang Youqi.Nonlinear vibrations of axially accelerating viscoelastic Timoshenko beams.Chinese Journal of Theoretical and Applied Mechanics,2013,45(6):965-973(in Chinese))

15 Zhao FQ,Wang MM.Vibration analysis of an axially moving viscoelastic beam.Materials,Mechanical Engineering and Manufacture,Pts 1-3,2013,268(270):1177-1181

16 丁虎,嚴巧赟,陳立群等.軸向加速運動黏彈性梁受迫振動中的混沌動力學.物理學報,2013,(20):78-84(Ding Hu,Yan Qiaoyun,Chen Liqun,et al.Chaotic dynamics in the forced nonlinear vibration of an axially accelerating viscoelastic beam.Acta Physica Sinica,2013,(20):78-84(in Chinese))

18 王冬梅,張偉,李慕榮等.用微分求積法分析軸向移動黏彈性梁的非平面非線性振動.動力學與控制學報,2015,(1):23-27(Wang Dongmei,Zhang Wei,Li Murong,et al.Using DQM to analyze the nonplanar nonlinear vibrationa of an axially moving viscoelastic beam.Journal of Dynamics and Control,2015,(1):23-27(in Chinese))

19 Lin W,Qiao N.Vibration and stability of an axially moving beam immersed in fl uid.International Journal of Solids and Structures,2008,45(5):1445-1457

20 Ni Q,Li MW,Tang M,et al.Free vibration and stability of a cantilever beam attached to an axially moving base immersed in fl uid.Journal of Sound and Vibration,2014,333(9):2543-2555

21 Li MW,Ni Q,Wang L.Nonlinear dynamics of an underwater slender beam with two axially moving supports.Ocean Engineering,2015,108:402-415

22 Yan H,Ni Q,Dai HL,et al.Dynamics and stability of an extending beam attached to an axially moving base immersed in dense fl uid.Journal of Sound and Vibration,2016,383:364-383

23 Spelsberg-Korspeter G,Kirillov ON,Hagedorn P.Modeling and stability analysis of an axially moving beam with frictional contact.Journal of Applied Mechanics-Transactions of the ASME,2008,75(3):031001-10

24 Guo XX,Wang ZM.Thermoelastic Coupling vibration characteristics of the axially moving beam with frictional contact.Journal of Vibration and Acoustics-Transactions of the ASME,2010,132(5):051010-7

25 Zheng P,Yang TZ,Yang XD,et al.An approximate analytical solution of an axially moving beam subjected to harmonic and parametric excitations simultaneously.Proceedings of First International Conference of Modelling and Simulation,2008,5:340-345

26 Yan QY,Ding H,Chen LQ.Nonlinear dynamics of axially moving viscoelastic Timoshenkobeam under parametric and external excitations.Applied Mathematics and Mechanics-English Edition,2015,36(8):971-984

27 WangL,ChenHH,HeXD.ActiveH-in fi nitycontrol ofthe vibration of an axially moving cantilever beam by magnetic force.Mechanical Systems and Signal Processing,2011,25(8):2863-2878

28 Zhang YW,Zhang Z,Chen LQ,et al.Impulse-induced vibration suppression of an axially moving beam with parallel nonlinear energy sinks.Nonlinear Dynamics,2015,82(1):61-71

29 Ghayesh MH,Kazemirad S,Amabili M.Coupled longitudinaltransverse dynamics of an axially moving beam with an internal resonance.Mechanism and Machine Theory,2012,52:18-34

30 Park S,Yoo HH,Chung J.Vibrations of an axially moving beam with deployment or retraction.AIAA Journal,2013,51(3):686-696

31 Yang XD,Zhang W.Nonlinear dynamics of axially moving beam with coupled longitudinal-transversal vibrations.Nonlinear Dynamics,2014,78(4):2547-2556

32 Stylianou M,Tabarrok B.Finite element analysis of an axially moving beam,Part I:time integration.Journal of Sound and Vibration,1994,178(4):433-453

33 Chang JR,Lin WJ,Huang CJ,et al.Vibration and stability of an axially moving Rayleigh beam.Applied Mathematical Modelling,2010,34(6):1482-1497

34 Ghayesh MH,Amabili M.Nonlinear vibrations and stability of an axially moving Timoshenko beam with an intermediate spring support.Mechanism and Machine Theory,2013,67:1-16

35 Park S,Chung JT.Dynamic analysis of an axially moving fi nitelength beam with intermediate spring supports.Journal of Sound and Vibration,2014,333(24):6742-6759

36 Oz HR,Pakdemirli M,Boyaci H.Non-linear vibrations and stability of an axially moving beam with time-dependent velocity.International Journal of Non-Linear Mechanics,2001,36(1):107-115

37 Al-Bedoor BO,Khulief YA.An approximate analytical solution of beam vibrations during axial motion.Journal of Sound and Vibration,1996,192(1):159-171

38 Matsuzaki Y.Vibration of a cantilevered beam during deployment and retrieval:analysis and experiment.Smart Materials&Structures,1998,4(4):334-339

39 Zhang H,Chen L.Vibration of an axially moving string supported by a viscoelastic foundation.Acta Mechanica Solida Sinica,2016,29(3):221-231

AN EFFICIENT DYNAMIC MODELING METHOD OF AN AXIALLY MOVING CANTILEVER BEAM AND FREQUENCY RESPONSE ANALYSIS1)

Hua Hongliang Liao Zhenqiang2)Zhang Xiangyan
(School of Mechanical Engineering，Nanjing University of Science and Technology，Nanjing210094，China)

The dynamics of the axially moving beam has wide application in engineering,such as robot manipulators,machine tools and gun barrel,et al.Computing the dynamic response of axially moving beam is an important method to evaluate the dynamic performance and fi nally the structure design.The time-varying motion equations of the axially movingcantileverbeamarederivedusingtheRayleigh-RitzmethodandLagrange’sequation.Thepowerseriesfunctionis used to construct the trial function to solve the dynamic problem.Due to the good integral and di ff erential performance of power series function,the derivation is easy to be carried out in the form of matrix.In this way,the symbolic computation software can generate the MATLAB program directly.And the generated MATLAB program can be used to conduct the dynamic computation with few modi fi cations,because the basic data unit of MATLAB is matrix.The overall process is efficiency and the time from dynamic modeling to computation is greatly reduced.Through four sets of numerical examples,the computational accuracy of the presented method is validated by comparing the dynamic responses with those from previous literatures.Then,the e ff ects of fi tting order of the power series function on computational accuracy are discussed.And the principle to select the fi tting order of the power series function to achieve good convergence and computational accuracy is given.Based on the dynamic model,the e ff ects of axial motion frequency on transverse vibration are studied.The e ff ects of axial vibration amplitude on the frequency response characteristic are explored.And the di ff erence between considering gravity and neglecting gravity e ff ect are compared.

computational dynamics，power series basis function，frequency response analysis，axially moving beam

O313

A doi：10.6052/0459-1879-17-182

2017–05–16 收稿，2017–08–07 錄用,2017–08–10 網(wǎng)絡(luò)版發(fā)表.

1)國家自然科學基金資助項目(51375241).

2)廖振強，教授，主要研究方向：自動武器發(fā)射動力學.E-mail:huahl123@126.com

華洪良，廖振強，張相炎.軸向移動懸臂梁高效動力學建模及頻率響應(yīng)分析.力學學報,2017,49(6):1390-1398

Hua Hongliang，Liao Zhenqiang，Zhang Xiangyan.An efficient dynamic modeling method of an axially moving cantilever beam and frequency response analysis.Chinese Journal of Theoretical and Applied Mechanics,2017,49(6):1390-1398