陳宏威　向建華　劉　帥　吳斌輝（北京理工大學(xué)機(jī)械與車輛學(xué)院北京100081）

·振動·噪聲·

曲拐彎扭耦合振動影響規(guī)律研究

陳宏威向建華劉帥吳斌輝
（北京理工大學(xué)機(jī)械與車輛學(xué)院北京100081）

基于圓盤非對稱轉(zhuǎn)子簡化模型和達(dá)朗貝爾原理，建立了曲軸單拐彎-扭耦合非線性振動數(shù)學(xué)模型，并利用解析推導(dǎo)法和四階龍格庫塔的迭代數(shù)值求解方法研究了彎曲振動和扭轉(zhuǎn)振動間的相互影響規(guī)律。研究結(jié)果表明，彎曲振動（或扭轉(zhuǎn)振動）會耦合產(chǎn)生多種頻率成分的扭轉(zhuǎn)振動（或彎曲振動），耦合振動幅值隨著激勵頻率和偏心距的增加而加大，且彎曲振動對扭轉(zhuǎn)振動的耦合影響相對較大。

曲拐彎扭耦合振動非線性振動影響規(guī)律

引言

內(nèi)燃機(jī)曲軸應(yīng)用范圍廣泛，結(jié)構(gòu)復(fù)雜，功能繁多，是內(nèi)燃機(jī)系統(tǒng)關(guān)鍵部件［1］。曲軸承受周期性變化的切向力和法向力作用，常常產(chǎn)生彎曲振動、扭轉(zhuǎn)振動及二者的耦合振動形式，很容易引發(fā)多種振動問題，導(dǎo)致軸段的疲勞損壞、軸承磨損以及工作噪聲等［2］。以往人們關(guān)于內(nèi)燃機(jī)曲軸的研究主要集中于軸的扭轉(zhuǎn)振動，但近年來隨著認(rèn)識的深入，一些學(xué)者開始著重關(guān)注扭轉(zhuǎn)—彎曲耦合振動，形成了比較完備的理論分析、數(shù)值模擬與實驗研究相結(jié)合的體系［3-5］。由于曲軸高速回轉(zhuǎn)時的質(zhì)量不平衡或系統(tǒng)不對稱等非線性因素，使得數(shù)值仿真中曲軸彎曲振動的數(shù)學(xué)建模和求解比較復(fù)雜［6-8］，其相應(yīng)的曲軸彎扭耦合振動分析的相關(guān)文獻(xiàn)目前尚不多見。李慧珍［9］在有限元計算曲軸振動中將曲軸視為圓截面梁和變截面矩形梁的組合；梁興雨［10］等建立了由多個自由度組成的發(fā)動機(jī)剛?cè)狁詈隙囿w動力學(xué)系統(tǒng)模型；何成兵［11］推導(dǎo)了三類微元軸段的彎扭耦合振動微分方程組，提出數(shù)值計算解法。Nataraj C［12］利用攝動的方法得到了非線性彎扭耦合振動微分方程組的初步解；梁明軒［13］分別通過諧波平衡法與數(shù)值積分法得到系統(tǒng)幅頻特性，分析了質(zhì)量偏心和轉(zhuǎn)速對軸系非線性動力學(xué)響應(yīng)的影響；劉占生［14］得到了共振區(qū)彎曲振動和扭轉(zhuǎn)振動的頻率特征，分析了不平衡量對彎扭耦合共振區(qū)的影響。目前，對于內(nèi)燃機(jī)曲軸的振動耦合分析，大多是基于計算精度較差的離散質(zhì)量或分布參數(shù)模型，且彎扭耦合振動影響等理論研究不夠深入。

本文針對某內(nèi)燃機(jī)曲軸單拐，并基于圓盤非對稱轉(zhuǎn)子簡化模型和達(dá)朗貝爾原理，建立了曲拐彎-扭耦合非線性振動數(shù)學(xué)模型，并利用解析推導(dǎo)和龍格庫塔的迭代數(shù)值求解方法深入研究了彎曲振動和扭轉(zhuǎn)振動間的相互影響規(guī)律，從而為內(nèi)燃機(jī)曲軸彎扭耦合研究提供理論指導(dǎo)。

1　非線性彎扭耦合模型建立和求解

對于截取的曲拐，在考慮扭轉(zhuǎn)振動的前提下，采用離散化模型分析彎曲振動運動狀況。將曲軸單拐簡化為一個剛性圓盤并安裝在具有彎曲和扭轉(zhuǎn)彈性而無質(zhì)量的軸上。假定軸是各向同性，兩端軸承為剛性支撐，忽略了圓盤的陀螺效應(yīng)，進(jìn)而將整個曲拐簡化為一個單質(zhì)量圓盤轉(zhuǎn)子模型，如圖1所示。

圖1　曲拐和簡化后的單質(zhì)量圓盤轉(zhuǎn)子模型

建立考慮扭轉(zhuǎn)振動耦合影響和偏心距的單質(zhì)量圓盤轉(zhuǎn)子示意圖如圖2所示。圓盤旋轉(zhuǎn)中心為O，質(zhì)心為O1，偏心距為e，φ為質(zhì)心旋轉(zhuǎn)的角度，ω為轉(zhuǎn)動角頻率，即φ=ωt+θ，θ為扭轉(zhuǎn)角。

圖2　圓盤轉(zhuǎn)子示意圖

由于偏心距e的存在，圓盤具有回轉(zhuǎn)速度φ˙和回轉(zhuǎn)加速度，使得質(zhì)心上有切向加速度eφ¨和法向加速度eφ˙2φ¨，如圖3所示。由此產(chǎn)生的慣性力為：

圖3　慣性力分析示意圖

同理，質(zhì)心具有彎曲加速度y¨和z¨，產(chǎn)生的慣性力為-my¨和-mz¨，由此施加在回轉(zhuǎn)中心的額外扭矩為：

在假設(shè)圓盤為剛體且左右支撐軸彎扭各向同性的前提下，根據(jù)達(dá)朗伯原理得到彎扭耦合振動平衡方程組：式中：T為作用在曲柄銷處的切向力產(chǎn)生的扭矩；Fy、Fz為作用在曲柄銷處的徑向力在y、z方向的投影；c和k為彎曲阻尼和剛度，此處不再考慮y和z方向的各向異性；ct和kt為扭轉(zhuǎn)阻尼和剛度。

從式（3）定性分析可以看出，耦合振動的產(chǎn)生是由于偏心距e的存在，e越大耦合作用就越大。而當(dāng)e=0時，圓盤的彎扭耦合振動平衡方程組變成：

即當(dāng)偏心距不存在時，彎曲扭轉(zhuǎn)振動間不存在明顯的耦合關(guān)系。

式（3）所示的彎扭耦合振動平衡方程組中變量z、y、θ均為時間的函數(shù)，呈現(xiàn)高度非線性。本文利用Matlab/Simulink搭建式（3）所對應(yīng)的彎扭耦合振動平衡方程組仿真模型（如圖4所示），并基于四階-龍格庫塔法（Runge-Kutta）迭代求解，獲得相應(yīng)的數(shù)值解。

圖4中，左右兩邊是模型輸入和輸出模塊，淺灰色的是兩個彎曲振動和1個扭轉(zhuǎn)振動計算模塊，深灰色的是彎扭和扭彎耦合計算模塊。仿真求解時，將初始值（作用力或力矩）賦給對應(yīng)的振動計算模塊，得到的值經(jīng)由各耦合模塊計算獲得新的耦合力或是力矩，再一并輸入作為總激勵來計算下一輪的值。以此反復(fù)迭代求解，直至滿足收斂條件，從而實現(xiàn)彎扭耦合振動平衡方程組的數(shù)值求解。

由于彎扭耦合影響是相互的，本文在分析彎扭耦合振動作用時，采用單獨控制一種振動而激發(fā)另一種振動的辦法，即通過單獨施加力或是力矩，觀察一種振動產(chǎn)生時另一種振動頻域下的幅值變化，從而獲得彎扭振動之間的耦合作用規(guī)律。

圖4　彎扭耦合振動求解模型

2　彎曲振動對扭轉(zhuǎn)振動的影響研究

由式（3）可以看出，由于偏心距e的存在，彎曲振動在圓盤上產(chǎn)生了額外扭矩M，即M=me（y¨sinφz¨cosφ），而重力作用產(chǎn)生的額外扭矩大小為megsin（tω+θ）。對M進(jìn)行傅立葉變換，將y、z方向的彎曲振動分解為若干簡諧振動的疊加，得到：

式中：ωi為第i個彎曲振動分量的振動角頻率；Yi、Zi分別為y、z兩個方向上振動的振幅；αi、βi為各自的相位。

將第i分量代入M=me（y¨sinφ-z¨cosφ），得到頻率為ωi的額外扭矩：

從式（6）可以看出，角頻率為ωi的彎曲振動可以引起角頻率為ωi±ω的扭轉(zhuǎn)振動；彎曲振動產(chǎn)生的額外扭矩的幅值與不平衡量me成正比，與彎曲振幅（Yi，Zi）成正比，與該諧分量下角頻率的平方成正比，彎曲振動相位角（αi，βi）對扭矩的幅值影響也較大。此外，由于重力作用也會產(chǎn)生額外扭矩，其大小為meg sin（ωt+θ），故重力會引起角頻率為ω的扭轉(zhuǎn)振動。

取單位曲拐模型仿真計算的基本參數(shù)為：轉(zhuǎn)速為3 000 r/min，即轉(zhuǎn)動頻率為50 Hz，k=1.736×1010N/m，J=0.045 kg·m2，kt=1.285×106（N·m）/rad，e= 0.02 m，m=9.07 kg，彎曲激勵頻率為35 Hz，幅值大小為Fy=1.0×107N，F(xiàn)z=5×105N，c=80（N·s）/m，Ct=100（N·m）/rad，通過圖4的數(shù)值迭代計算，結(jié)果如圖5所示。

從圖5可以看出，在激勵力35 Hz作用下，彎曲振動在35 Hz產(chǎn)生了明顯的波峰。同時，在扭轉(zhuǎn)振動的15 Hz、50 Hz、85 Hz附近產(chǎn)生較為明顯的峰值；耦合作用產(chǎn)生的扭轉(zhuǎn)振動幅值和彎曲振動幅值比（無量綱）為1∶40。扭轉(zhuǎn)振動在15 Hz和85 Hz下產(chǎn)生的峰值符合│ω±ωi│的規(guī)律，即為激勵力頻率與轉(zhuǎn)動頻率之差或之和的絕對值。50 Hz處的峰值產(chǎn)生是因為重力所導(dǎo)致的不平衡引起的，重力可以產(chǎn)生meg sin（ωt+θ）的額外扭矩進(jìn)而引起頻率為ω的扭轉(zhuǎn)振動。值得注意的是，雖然由彎曲振動產(chǎn)生的扭矩M通過傅里葉變換可以得出N個不同角頻率的簡諧振動，進(jìn)而引起相應(yīng)的一系列的扭轉(zhuǎn)振動；但只有峰值較為明顯（通常是低階）的彎曲振動，耦合產(chǎn)生的扭轉(zhuǎn)振動才較為明顯。其余情況下產(chǎn)生的扭轉(zhuǎn)振動可以忽略。

此外，在轉(zhuǎn)速一定的情況下，隨著彎曲振動激勵力頻率從55 Hz到90 Hz的變化，扭轉(zhuǎn)振動幅值也呈一定的變化。其中圖6顯示了│ω-ωi│扭轉(zhuǎn)振動隨著激勵力ω變化關(guān)系。

從圖6中可以看出，隨著彎曲振動激勵力頻率從55 Hz到90 Hz的變化，頻率│ω-ωi│扭轉(zhuǎn)振動幅值逐漸增加。因為根據(jù)公式（6）可知，彎曲振動耦合產(chǎn)生的額外扭矩大小的決定因素中有激勵力頻率，且與其平方成正比，同時受重力產(chǎn)生的旋轉(zhuǎn)扭矩的影響，但總體上扭轉(zhuǎn)振動會隨著激勵力頻率的增加而相應(yīng)的增加。特別指出，當(dāng)扭轉(zhuǎn)固有頻率和彎曲固有頻率都不大時，會出現(xiàn)由于彎曲振動耦合作用所導(dǎo)致的扭轉(zhuǎn)共振的情況發(fā)生。

圖5　彎曲振動和耦合產(chǎn)生的扭轉(zhuǎn)振動頻域圖

圖6　│ω-ωi│扭轉(zhuǎn)振動隨著激勵力ω變化關(guān)系

3　扭轉(zhuǎn)振動對彎曲振動的影響研究

由式（3）可以看出，由于存在偏心距和扭轉(zhuǎn)角的原因，在圓盤部產(chǎn)生了切向力和法向力，這些激勵力作用在彎曲振動方向，使其產(chǎn)生相應(yīng)的彎曲振動。這里，以Z方向彎曲振動為研究對象；為進(jìn)一步分析微分方程，可以把扭轉(zhuǎn)角θ分解為若干諧振動的和，即

扭對彎耦合項Z方向力me（φ˙2sinφ-φ¨cosφ）進(jìn)行分析，即：

式中：第1項是由于旋轉(zhuǎn)不平衡引起的，為傳統(tǒng)的不平衡彎曲振動響應(yīng)激勵力，它將引起同頻率的彎曲振動。后邊各項均是由于扭轉(zhuǎn)振動引起的離心激勵力，頻率為：

1）ω±ωi，幅值與、不平衡me、轉(zhuǎn)速ω成正比，與ωi或是成正比。

3）ω±ωi±ωj，幅值與乘積、不平衡me、ωiωj乘積成正比。

4）2ωi±ω，幅值與、不平衡me成正比。

通過對比可知，這些激勵力由于幅值引入了扭振角，扭轉(zhuǎn)振動引起的慣性力相比不平衡直接引起的同頻率彎曲振動幅值要小很多；而其中ω±ωi頻率下對應(yīng)的振幅相對大一些。

同樣，取單位曲拐模型仿真計算的基本參數(shù)為：轉(zhuǎn)速為3 000 r/min，即轉(zhuǎn)動頻率為50 Hz，k=1.736× 1010N/m，J=0.045 kg·m2，kt=1.285×106（N·m）/rad，e=0.02 m，m=9.07 kg，扭矩頻率為35 Hz，幅值大小為T=2e4N·m，c=80（N·s）·m-1，ct=100（N·m·s）·rad-1。通過圖4的數(shù)值迭代計算，結(jié)果如圖7～圖9所示。

圖7　扭轉(zhuǎn)振動頻域圖

圖8　耦合引起的同頻率彎曲振動幅值

圖9　耦合引起的ω±ωi頻率彎曲振動幅值

從圖7可以看出，在35 Hz扭矩的作用下，產(chǎn)生了比較明顯的波峰；圖9中彎曲振動在15 Hz、50 Hz和85Hz附近產(chǎn)生了明顯的波峰；耦合作用產(chǎn)生的彎曲振動幅值和扭轉(zhuǎn)振動幅值比（無量綱）為1∶12 000。圖8中在彎曲振動50 Hz處產(chǎn)生明顯峰值的原因是傳統(tǒng)的同頻彎曲振動作用的結(jié)果，其幅值與不平衡成正比，與轉(zhuǎn)速ω平方成正比。圖9中彎曲振動在15 Hz和85 Hz處產(chǎn)生峰值的原因與│ω±ωi│相符合，即是轉(zhuǎn)動角頻率與扭轉(zhuǎn)振動頻率之差或之和，其幅值與不平衡me成正比，與轉(zhuǎn)速ω成正比，與ωi或是ω2i成正比，與θi成正比。由于這里的扭轉(zhuǎn)剛度很大，故θi一般比較小，使得幅值要遠(yuǎn)小于同頻率引起的彎曲振動；而其它頻率成分對應(yīng)的幅值更小，耦合作用不明顯，可以忽略。

圖10和圖11描述了在軸系轉(zhuǎn)速頻率為50 Hz，外扭矩頻率從55 Hz逐漸增加到90 Hz時，扭轉(zhuǎn)振動引起的│ω±ωi│頻率的彎曲振動隨著外扭矩頻率的變化趨勢。隨著扭轉(zhuǎn)激勵力頻率的增加，彎曲振動逐漸增加，這與彎曲振動對扭轉(zhuǎn)振動的結(jié)論相似。

圖10　Y向彎曲振動隨扭振頻率的變化

圖11　Z向彎曲振動隨扭振頻率的變化

4　結(jié)論

本文建立了曲軸單拐彎-扭耦合振動的數(shù)學(xué)模型，并針對彎-扭耦合振動和扭-彎耦合振動的影響規(guī)律進(jìn)行了解析推導(dǎo)和數(shù)值迭代求解，計算結(jié)果具有很好的一致性。具體結(jié)論如下：

1）彎曲振動（或扭轉(zhuǎn)振動）會耦合產(chǎn)生多種頻率成分的扭轉(zhuǎn)振動（或彎曲振動），都包括頻率為ω±ωi的振動成分。

2）耦合產(chǎn)生的振動幅值隨著激勵頻率和偏心距的增加而加大。

3）無量綱下的彎扭耦合幅值比遠(yuǎn)大于扭彎耦合幅值比，故本曲軸單拐的彎曲振動對扭轉(zhuǎn)振動的耦合影響相對較大。

1Fuchun Ren，Kun Yang，Zhong Guang Fu，et al.Study on the coupling of the bending and vibrations of the unbalanced rotor［C］：Xi'an：Proceedings of the 2nd International Conference on Hydrodynamic Bearing-rotor System Dynamics，March 20-22，1997

2李震，桂長林，孫軍.內(nèi)燃機(jī)曲軸軸系振動分析研究的現(xiàn)狀、討論與展望［J］.內(nèi)燃機(jī)學(xué)報，2002，20（05）：469-474

3Nevsan S.Performance increase in turbo molecular pumps with curved type blades［J］.Vacuum，2012，86（11）：1764-1769

4韓建鑫，王煒，張琪昌.內(nèi)燃機(jī)曲軸系統(tǒng)變慣量公式的修正［J］.內(nèi)燃機(jī)學(xué)報，2014，32（03）：271-275

5Munawar I，Abdul W，Dimitri B，et al.Design modification in rotor blade of turbo molecular pump［J］.Nuclear Instruments and Methods in Physics Research A，2012，678：88-90

6張勇，蔣滋康.軸系彎扭耦合振動的數(shù)學(xué)模型［J］.清華大學(xué)學(xué)報（自然科學(xué)版），1998（08）：115-118

7周桐，徐健學(xué)，傅衛(wèi)平.彎扭耦合振動對次同步諧振響應(yīng)的影響［J］.應(yīng)用力學(xué)學(xué)報，2000，17（1）：12-17

8王雪嬌.內(nèi)燃機(jī)曲軸彎扭耦合非線性振動分析［D］.天津：天津大學(xué)，2009

9李慧珍，張德平.用有限元進(jìn)行曲軸扭轉(zhuǎn)振動計算［J］.內(nèi)燃機(jī)學(xué)報，1991，9（2）：157-162

10梁興雨，舒歌群，李東海，等.基于柔體曲軸多體動力學(xué)的軸系扭振響應(yīng)分析［J］.內(nèi)燃機(jī)工程，2007（04）：46-49

11何成兵，楊昆，顧煜炯.質(zhì)量偏心對碰磨轉(zhuǎn)子彎曲振動和扭轉(zhuǎn)振動特性的影響［J］.中國電機(jī)工程學(xué)報，2002，22（7）：105-110

12 Nataraj C.On the Interaction of Torsion and Bending in Rotating Shafts［J］.Journal of Applied Mechanics，1993，60（3）：239-241

13梁明軒，袁惠群，趙天宇，等.內(nèi)燃機(jī)曲軸軸系彎扭耦合非線性振動響應(yīng)［J］.東北大學(xué)學(xué)報（自然科學(xué)版），2014（03）：402-405，414

14劉占生，崔穎，黃文虎，等.轉(zhuǎn)子彎扭耦合非線性動力學(xué)特性研究［J］.中國機(jī)械工程，2003，14（7）：603-605

The Nonlinear Study on Coupled Bending and Torsion Vibration of Crank Throw

Chen Hongwei，Xiang Jianhua，Liu Shuai，Wu Binhui
School of Mechanical Engineering，Beijing Institute of Technology（Beijing，100081，China）

Based on the simplified model of asymmetric disc rotor and D'alembert's principle，the nonlinear mathematical model of coupled bending and torsion vibration is established，and then the mutual influence between the bending vibration and torsion vibration is studied through the theoretical derivation and the solving methods of iterative numerical of fourth-order Runge-Kutta.The results show that the bending vibration（or torsion vibration）will produce various frequency of torsion vibration（or bending vibration）if coupled.The coupling vibration amplitude will raise with the increasing of the excitation frequency and eccentricity，the bending vibration has a relatively large coupling influence on torsion vibration.

Crank throw，Coupled bending and torsion vibration，Nonlinear vibration，Influence rule

TK422

A

2095-8234（2016）04-0065-06

2016-02-03）

陳宏威（1991-），男，碩士研究生，主要研究方向為內(nèi)燃機(jī)結(jié)構(gòu)振動。

向建華（1976-），男，博士，副教授。