任意邊界條件彈性桿結(jié)構(gòu)扭轉(zhuǎn)振動特性分析

許得水， 杜敬濤， 李文達， 楊鐵軍， 李玩幽

(哈爾濱工程大學 動力與能源工程學院, 哈爾濱 150001)

任意邊界條件彈性桿結(jié)構(gòu)扭轉(zhuǎn)振動特性分析

許得水， 杜敬濤， 李文達， 楊鐵軍， 李玩幽

(哈爾濱工程大學 動力與能源工程學院, 哈爾濱 150001)

采用改進傅里葉級數(shù)方法建立了任意邊界條件彈性桿扭轉(zhuǎn)振動特性預報模型。針對傳統(tǒng)傅里葉級數(shù)在扭振邊界處存在的位移導數(shù)不連續(xù)問題，通過改進傅里葉級數(shù)的方法改善解的收斂性和準確性。彈性桿結(jié)構(gòu)扭振微分方程與任意邊界條件方程進行聯(lián)合求解，得到彈性桿扭振問題的特征矩陣方程。數(shù)值算例分析結(jié)果充分驗證了本文模型的可行性與正確性。

彈性桿；扭轉(zhuǎn)振動；邊界條件；傅里葉級數(shù)

彈性桿(或軸)結(jié)構(gòu)廣泛應用于船舶推進、石油鉆探、傳動軸系等眾多工程場合，當主動力矩與負荷反力矩存在不平衡時，容易引起彈性桿結(jié)構(gòu)繞軸線進行的扭轉(zhuǎn)振動。全面理解彈性桿結(jié)構(gòu)扭振特性是對此類復雜結(jié)構(gòu)系統(tǒng)優(yōu)化設計及其振動控制的重要基礎與前提。為此，國內(nèi)外眾多學者圍繞彈性桿扭轉(zhuǎn)振動分析開展了大量的研究。

TEFFT等[1]通過實驗測試研究了方形截面均勻桿結(jié)構(gòu)的扭轉(zhuǎn)共振行為，扭振固有頻率采用比例因子進行表征，他們發(fā)現(xiàn)該參數(shù)隨彈性桿結(jié)構(gòu)截面積與長度比增大而增大。PAUL等[2]分析了彈性圓桿結(jié)構(gòu)扭振強迫響應問題，圓桿一端與彈性半空間自由表面相連，另一端受到簡諧扭矩作用，頻響曲線表明阻尼效應隨彈性桿結(jié)構(gòu)的長徑比減小而增強。RAO[3]考慮帶有彈性支撐的軸系與管路系統(tǒng)扭轉(zhuǎn)振動特性，得到了各種情況下的扭振固有頻率參數(shù)。劉清友等[4]利用差分法求解了鉆柱扭轉(zhuǎn)振動模型。

劉志先[5]采用Bessel函數(shù)法對楔形直桿的縱振和扭振固有頻率進行了推算。RAFIEE等[6]利用Bessel及Neumann特殊函數(shù)建立了變截面彈性桿扭轉(zhuǎn)自由振動解析解。楊立軍等[7]采用攝動法研究邊界面桿軸向和扭轉(zhuǎn)固有振動，結(jié)合邊界條件得到固有頻率的特征方程及頻率參數(shù)。CHONDROS[8]根據(jù)變分原理建立了含有周向裂紋的彈性圓桿扭轉(zhuǎn)振動分析模型。WEI等[9]研究了使用壓電式執(zhí)行機構(gòu)和LQG狀態(tài)反饋對彈性桿結(jié)構(gòu)扭轉(zhuǎn)振動進行主動抑制的可行性。BARR[10]研究了非圓形均勻彈性桿的扭轉(zhuǎn)波動行為特性。舒歌群[11]基于扭轉(zhuǎn)彈性波理論分析了船舶柴油機推進軸系的扭轉(zhuǎn)振動。

近年來，隨著材料科學的發(fā)展，具有各種材料特性的彈性桿扭轉(zhuǎn)振動也受到了一定關注。YANG[12]針對大變形情況下的圓柱形壓電桿結(jié)構(gòu)，采用單模態(tài)Galerkin近似，對其非線性扭轉(zhuǎn)振動行為進行了分析。李立群等[13]基于非局部理論研究了具有宏觀尺度的非局部彈性直桿的扭轉(zhuǎn)振動問題，采用傳遞函數(shù)法對常見邊界條件下非局部彈性直桿扭轉(zhuǎn)振動的固有頻率及振型進行了求解。NARENDAR等[14]分析了微納米尺度桿應變梯度扭轉(zhuǎn)振動特性。ANSARI等[15]從分子動力學角度研究了氮化硼納米管扭轉(zhuǎn)行為。郁殿龍等[16]分析了聲子晶體桿結(jié)構(gòu)扭轉(zhuǎn)振動帶隙特性。

在現(xiàn)有文獻中，針對彈性桿結(jié)構(gòu)振動特性分析多集中為理想邊界條件，如自由或固定，所得到的結(jié)論具有較大的局限性。事實上，邊界條件作為重要的結(jié)構(gòu)參數(shù)，對彈性結(jié)構(gòu)振動特性具有顯著影響。楊天智等[17]采用微分求積法求解了兩端彈性支承管道橫向振動固有頻率。劉向堯等[18]運用傳遞矩陣法建立了復雜邊界條件下Timoshenko梁結(jié)構(gòu)自由振動模型。本文采用一種改進傅里葉級數(shù)方法對彈性桿扭轉(zhuǎn)振動控制微分方程及任意邊界條件進行聯(lián)合求解，建立了任意邊界條件下彈性桿扭振特性分析模型。通過與現(xiàn)有解析解和有限元法得到的結(jié)果對比，充分驗證了本模型的正確性。在此基礎上，探討了邊界約束剛度系統(tǒng)對彈性桿扭振特性的影響。

1 理論模型

1.1 彈性扭轉(zhuǎn)約束邊界

如圖1所示，考慮彈性桿結(jié)構(gòu)自由扭轉(zhuǎn)振動問題，其中，桿長為L，橫截面半徑為R，扭轉(zhuǎn)角位移用θ(x,t)表示。通過在彈性桿結(jié)構(gòu)左右兩端引入扭轉(zhuǎn)約束彈簧來模擬任意邊界條件，對應的剛度系數(shù)分別為K0和KL。相應地，經(jīng)典邊界條件如自由和固支便可以通過將彈性剛度系統(tǒng)設置為零或無窮大而統(tǒng)一得到，當剛度系數(shù)取某有限值時，則可以表示任意彈性約束邊界條件。

彈性桿扭轉(zhuǎn)振動θ(x,t)可以表示為：

(1)

式中:Θ(x)為彈性桿扭振位移分布函數(shù)，ejωt為簡諧時間因子，ω是圓頻率，在后續(xù)推導過程中為了簡化，將忽略該時間因子。

圖1 任意邊界條件彈性桿扭轉(zhuǎn)振動分析模型Fig.1 Torsional vibration analysis model of elastic rod structure with arbitrary boundary conditions

根據(jù)彈性桿結(jié)構(gòu)兩端的力矩平衡關系，邊界條件方程可以寫為如下形式：

(2a)

(2b)

經(jīng)典邊界條件下彈性桿扭轉(zhuǎn)位移函數(shù)通常可以采用傅里葉級數(shù)進行展開：

(3)

這里，λn=nπ/L。

根據(jù)傅里葉級數(shù)的求導性質(zhì)，可以發(fā)現(xiàn)，傳統(tǒng)的級數(shù)展開形式將使得位移函數(shù)的一階導數(shù)在端點位置始終為零，即：

(4)

(5)

顯然，通過對比式(1)可知，傳統(tǒng)傅里葉級數(shù)展開形式，由于在端點位置上一階導數(shù)不連續(xù)，不具備表示任意彈性約束情況下內(nèi)力矩的能力，將很大程度上影響級數(shù)解的預報精度和收斂速度。

為了使傅里葉級數(shù)在包含端點在內(nèi)的整個求解域[0,L]上所有相關導數(shù)足夠連續(xù)(或函數(shù)足夠光滑)，我們在傳統(tǒng)傅里葉級數(shù)基礎上引入兩項附加多項式函數(shù)，以去除邊界位置的導數(shù)跳躍，相應的改進傅里葉級數(shù)形式如下：

(6)

式中:ζ1(x)和ζ2(x)是彈性桿扭振角位移輔助函數(shù)。這種情況下，位移場函數(shù)在端點位置上一階導數(shù)不再恒等于零，從而具備表示任意彈性約束情況下內(nèi)力矩的能力，有助于改善級數(shù)解的精確性和收斂性[19]。此處，相應的輔助函數(shù)構(gòu)造如下：

(7)

可以證明：

(8)

(9)

將位移表達式(6)代入上述方程并整理可得：

(10)

在實際計算中，可以采用級數(shù)截斷為有限維，從而上述線性方程組可以進一步寫成矩陣形式：

HA=QB

(11)

式中:

(12)

A={a0a1…aN}T

(13)

(14)

B={b1b2}T

(15)

從式(11)可以看出，改進傅里葉級數(shù)旋轉(zhuǎn)角位移形式(6)中傳統(tǒng)傅里葉級數(shù)系數(shù)向量A與輔助多項式系數(shù)B之間并非任意關系，而是取決于扭轉(zhuǎn)桿結(jié)構(gòu)兩端彈性約束剛度情況。

1.2 彈性桿扭轉(zhuǎn)振動方程求解

彈性桿扭轉(zhuǎn)振動偏微分方程為：

(16)

式中:β2=ω2J/GIp，ω為彈性桿扭振圓頻率，J為單位長度轉(zhuǎn)動慣量，GIp為扭轉(zhuǎn)剛度。

將彈性桿扭轉(zhuǎn)角位移改進傅里葉級數(shù)形式(6)代入上述扭轉(zhuǎn)振動控制微分方程，可以得到：

(17)

將式(17)中輔助函數(shù)及其二階空間導數(shù)進一步采用傅里葉余弦級數(shù)展開，利用三角函數(shù)系數(shù)相等可以得到：

(18)

這里，α1n,α2n,ε1n,ε2n分別為所構(gòu)建輔助函數(shù)及其二階空間導數(shù)形式傅里葉展開系數(shù)，相應定義如下：

(19)

(21)

(22)

(23)

(24)

(25)

(26)

對方程(17)中級數(shù)進行截斷，可以將其寫為更為簡潔的矩陣形式：

[CA+BD]+β2[EA+FB]=0

(27)

式中，

(28)

(29)

(30)

(31)

結(jié)合彈性約束邊界方程(11)，任意邊界條件彈性桿結(jié)構(gòu)扭轉(zhuǎn)振動特征方程為：

(K+β2M)A=0

(32)

式中:K=C+DTQ-1H，M=E+FTQ-1H。

通過求解式(32)這一標準的矩陣特征值問題可以得到彈性約束邊界條件桿結(jié)構(gòu)扭轉(zhuǎn)振動模態(tài)參數(shù)，將所對應的特征向量代入至彈性桿扭振角位移改進傅里葉級數(shù)形式(6)中，即可得到各種邊界條件下彈性桿扭轉(zhuǎn)模態(tài)振型。

2 數(shù)值算例分析

2.1 兩端固定邊界

正如前文所述，在計算仿真中需要對改進傅里葉級數(shù)項數(shù)進行截斷，表2中所列數(shù)據(jù)為固定邊界條件彈性桿前12階固有頻率在不同截斷項數(shù)截斷數(shù)計算所得結(jié)果。從此表能夠看出，隨著截斷數(shù)增加彈性桿扭振模態(tài)頻率計算精度逐漸提高，當截斷數(shù)分別取35和45時，頻率計算結(jié)果已基本收斂，為此在隨后數(shù)值仿真計算中，所有傅里葉級數(shù)均取前35項。

2.2 固定-自由邊界

對于固定-自由邊界條件，我們只需要將彈性桿結(jié)構(gòu)左端邊界約束剛度系數(shù)設置為1×106，而右端剛度系數(shù)取為零即可。

表1 模擬固定邊界條件無限大彈簧約束剛度系數(shù)收斂情況

表2 改進傅里葉級數(shù)截斷項數(shù)收斂情況

通過將特征方程求解所得到特征向量系統(tǒng)代入式(6)和(11)便可得到固定-自由邊界條件彈性桿扭轉(zhuǎn)振動模態(tài)振型。圖2繪制出本文模態(tài)振型分布與解析解比較，可以看出兩者能夠很好地吻合。對于一端固定一端自由邊界，解析解模態(tài)函數(shù)為[20]：

(33)

式中:i為模態(tài)序數(shù)。

圖2 固支-自由邊界條件彈性桿扭振模態(tài)振型圖Fig.2 Torsional mode shapes of elastic rod with clamped and free boundary conditions

2.3 彈性-自由邊界

為了更好地揭示此種彈性-自由邊界條件彈性桿扭轉(zhuǎn)振動固有頻率隨無量剛彈性剛度系數(shù)的變化，表3給出了不同邊界約束剛度下彈性-自由邊界條件本文模型扭轉(zhuǎn)振動固有頻率變化情況?？梢园l(fā)現(xiàn)，隨著彈性約束剛度逐漸增加，各階固有頻率逐漸增大，圖3更為直觀地反映了這一變化趨勢，當邊界約束剛度增大到一定程度時，固有頻率取值趨于穩(wěn)定。

表3 不同無量綱剛度系數(shù)下彈性桿扭轉(zhuǎn)振動固有頻率

圖3 前五階彈性桿扭振固有頻率隨彈性約束剛度變化情況Fig.3 Variation of first five torsional natural frequencies of elastic rod with different non-dimensional boundary restraining stiffnesses

2.4 任意彈性邊界

模態(tài)階數(shù)ANSYS/Hz本方法/Hz偏差/%187.16887.090.0942290.69290.640.0183539.81539.770.0084797.02796.960.00751056.71056.560.01361317.31317.160.01171578.51578.280.01481840.11839.690.02292101.92101.290.029102363.92363.010.038112626.02624.820.045122888.32886.700.055

圖4 邊界約束剛度對彈性桿扭轉(zhuǎn)振動特性的影響Fig4 Influence of boundary restraining stiffness on the torsional vibration characteristics of elastic rod

3 結(jié) 論

本文采用改進傅里葉級數(shù)方法建立了任意邊界條件彈性桿扭轉(zhuǎn)振動特性預報模型。通過在彈性桿兩端引入扭轉(zhuǎn)邊界約束彈簧，所有邊界條件都可以通過設置相應的剛度系數(shù)獲得。針對傳統(tǒng)傅里葉級數(shù)在扭振邊界處存在的位移導數(shù)不連續(xù)問題，通過改進傅里葉級數(shù)的方法改善解的收斂性和精確性。結(jié)合彈性邊界條件，實現(xiàn)對彈性桿控制微分方程的精確求解，通過求解一個標準的矩陣特征值問題得到系統(tǒng)的模態(tài)參數(shù)。

在數(shù)值算例中，應用MATLAB語言編程仿真，計算了各種邊界條件下彈性桿扭轉(zhuǎn)振動模態(tài)參數(shù)，通過與文獻中解析解和有限元ANSYS結(jié)果進行對比，充分驗證了本文方法對于分析不同邊界條件桿結(jié)構(gòu)扭轉(zhuǎn)振動特性的可靠性和有效性。當邊界條件發(fā)生改變時，本文模型不需要像有限元那樣對模型進行修正，或解析解那樣進行重新推導，僅須設置邊界約束剛度系數(shù)即可。在此基礎上，探討了邊界約束剛度對彈性桿扭振模態(tài)頻率的影響規(guī)律，結(jié)果表明彈性桿兩端邊界約束剛度對扭振固有頻率的影響存在一個敏感區(qū)域，當邊界約束剛度越過此區(qū)域時，將不再具有顯著影響；在此范圍內(nèi)，合適改變邊界約束剛度能夠調(diào)整系統(tǒng)某階固有頻率，進而可以避免工程實際中外部激勵載荷下發(fā)生共振。

[ 1 ] TEFFT W E, SPINNER S. Torsional resonance vibrations of uniform bars of square cross section [J]. Journal of Research of the National Bureau of Standards-A. Physics and Chemistry, 1961, 65(3): 167-171.

[ 2 ] PAUL H S, SRIDHARAN K. Forced torsional vibration of an elastic circular rod on an elastic half space and on anelastic stratum [J]. Journal of the Acoustical Society of America, 1979, 65(2): 391-398.

[ 3 ] RAO C K. Torsional frequencies and mode shapes of generally constrained shafts and piping [J]. Journal of Sound and Vibration, 1988, 125(1): 115-121.

[ 4 ] 劉清友, 馬德坤, 鐘青. 鉆柱扭轉(zhuǎn)振動模型的建立及求解[J]. 石油學報, 2000, 21(2): 78-82. LIU Qingyou, MA Dekun, ZHONG Qing. A drilling string torsional vibration model and its solution[J]. ACTA, 2000, 21(2): 78-82.

[ 5 ] 劉志先. 推算楔形直桿縱振與扭振自主頻率的一個方法[J]. 中國科學, 1977, 7(6): 536-546. LIU Zhixian. A method for calculating the frequency of the longitudinal and torsional vibration of a wedge shaped rod [J]. Science China, 1977, 7(6): 536-546.

[ 6 ] RAFIEE M, MEHRABADI S J, RASEKH-SALEH N. Analytical solutions for the torsional vibrations of variable cross-section rods [J]. International Journal of Engineering and Applied Science, 2010, 2(4): 64-71.

[ 7 ] 楊立軍, 鄧志恒, 葉柏龍, 等. 變截面桿軸向和扭轉(zhuǎn)振動分析的攝動法[J]. 中南大學學報, 2014, 45(3): 855-861. YANG Lijun, LIU Zhiheng, YE Bailong, et al. Perturbation method for longitudinal and torsional vibration analysis of bar with variable cross-section[J]. Journal of Central South University, 2014, 45(3): 855-861.

[ 8 ] CHONDROS T G. Variational formulation of a rod under torsional vibration for crack identification [J]. Theoretical and Applied Fracture Mechanics, 2005, 44: 95-104.

[ 9 ] WEI Y D, Lü Y G, CHEN Z C.Active control of torsional vibration of flexible bars [J]. International Journal of Information and Systems Sciences, 2005, 1(3): 347-354.

[10] BARR A D S. Torsional waves in uniform rods of non-circular section [J]. Journal of Mechanical Engineering Science, 1962, 4(2): 127-135.

[11] 舒歌群. 基于扭轉(zhuǎn)彈性波理論的船舶柴油機推進軸系扭振研究[J]. 船舶力學, 2005, 9(5): 137-142. SHU Gequn. Calculation of crankshaft torsional vibration response in diesel engine by torsional elastic wave theory[J]. Journal of Ship Mechanics, 2005, 9(5): 137-142.

[12] YANG J S. Nonlinear torsional vibration of a circular cylindrical piezoelectric rod with relatively large shear deformation [J]. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, 2007, 54(7): 1482-1485.

[13] 李立群, 段靜波, 李家文, 等. 基于傳遞函數(shù)法的非局部彈性直桿扭轉(zhuǎn)振動研究[J]. 力學季刊, 2012, 33(2): 181-187. LI Liqun, DUAN Jingbo, LI Jiawen, et al. Research on torsional vibration of nonlocal elastic bar using transfer function method[J]. Chinese Quarterly of Mechanics, 2012, 33(2): 181-187.

[14] NARENDAR S, RAVINDER S,Gopalakrishnan S. Strain gradient torsional vibration analysis of micro/nano rods [J]. International Journal of Nano Dimension, 2012, 3(1): 1-17.

[15] ANSARI R, AJORI S.Molecular dynamics study of the torsional vibration characteristics of boron-nitride nanotubes [J]. Physics Letters A, 2014, 378: 2876-2880.

[16] 郁殿龍, 劉耀宗, 王剛, 等. 一維桿狀結(jié)構(gòu)聲子晶體扭轉(zhuǎn)振動帶隙研究[J]. 振動與沖擊, 2006, 25(1): 104-106. YU Dianlong, LIU Yaozong, WANG Gang, et al. Research on torsional vibration band gaps of one dimensional phononic crystals composed of rod structures[J]. Journal of Vibration and Shock, 2006, 25(1): 104-106.

[17] 楊天智，方勃，楊曉東，等. 復雜約束管道固有頻率的快速算法[J]. 哈爾濱工業(yè)大學學報, 2007, 39(5): 771-773. YANG Tianzhi, FANG Bo, YANG Xiaodong, et al. A fast algorithm for natural frequencies of pipe with complex restrictions[J]. Journal of Harbin Institute of Technology, 2007, 39(5): 771-773.

[18] 劉向堯，聶宏，魏小輝. 復雜邊界條件下的多跨梁的振動模型[J]. 北京航空航天大學學報, 2015, 41(5): 841-846. LIU Xiangyao, NIE Hong, WEI Xiaohui, et al. Vibration model for multi-span beam with arbitrary complex boundary conditions[J]. Journal of Beijing University of Aeronautics and Astronautics, 2015, 41(5): 841-846.

[19] DU J T, LI W L, JIN G Y,et al. An analytical method for the in-plane vibration analysis of rectangular plates with elastically restrained edges[J]. Journal of Sound and Vibration, 2007, 306: 908-927.

[20] HARRIS C M, PIERSOL A G. Harris’ Shock and Vibration Handbook 5thEd[A]. 劉樹林等譯. 沖擊與振動手冊[M].5版.北京：中國石化出版社, 2008：1-8.

Torsional vibration characteristics of an elastic rod structure under arbitrary boundary conditions

XU Deshui, DU Jingtao, LI Wenda, YANG Tiejun, LI Wanyou

(College of Power and Energy Engineering, Harbin Engineering University, Harbin 150001, China)

An improved Fourier series method was employed to model torsional vibration of an elastic rod under arbitrary boundary conditions. In order to overcome discontinuities of displacement derivatives of traditional Fourier series at boundary points, an improved method was constructed to improve the convergence and correctness of the series solution. The system characteristic equation was obtained through exactly solving the torsional vibration governing equations and boundary condition equations of the rod. Various numerical examples were presented to validate the feasibility and correctness of the proposed model.

elastic rod; torsional vibration; boundary condition; Fourier series

國家自然科學基金資助項目(11202056)

2015-05-29 修改稿收到日期：2015-11-17

許得水 男，博士生，1989年生

杜敬濤 男，教授，博士生導師，1981年生 E-mail: dujingtao@hrbeu.edu.cn

O327

A

10.13465/j.cnki.jvs.2017.01.024