魚雷渦輪轉子-滾動軸承系統(tǒng)非線性動力學特性分析

宋 文, 楊赪石, 陳志偉, 孫 巖, 任 鑫, 李 剛, 柴曉帥

魚雷渦輪轉子-滾動軸承系統(tǒng)非線性動力學特性分析

宋 文1, 楊赪石1, 陳志偉2, 孫 巖2, 任 鑫2, 李 剛2, 柴曉帥2

(1. 中國船舶集團有限公司 第705研究所, 陜西 西安, 710077; 2. 山西平陽重工機械有限責任公司, 山西 臨汾, 043003)

針對魚雷發(fā)動機渦輪轉子系統(tǒng), 考慮軸承徑向游隙, 非線性赫茲接觸剛度以及轉子不平衡量的基礎上, 建立了懸臂轉子-滾動軸承系統(tǒng)的動力學模型, 并基于變步長的龍格庫塔方法對所建立的動力學模型進行了數(shù)值仿真。計算了雙參數(shù)平面最大Lyapunov指數(shù)譜圖, 得到了轉子系統(tǒng)隨雙參數(shù)變化時的周期和非周期運動的分布規(guī)律, 并進行了系統(tǒng)在特定參數(shù)條件下的時域特性、頻域特性及周期特性分析。分析表明, 隨著軸承徑向游隙增加, 進入混沌區(qū)域的最低轉速在不斷減低。該研究可為魚雷渦輪機系統(tǒng)設計提供參考。

魚雷; 發(fā)動機; 渦輪轉子; 非線性動力學; 接觸剛度; 雙參數(shù)

0 引言

渦輪機是熱動力魚雷的重要推進形式之一, 相對于魚雷活塞發(fā)動機, 魚雷渦輪發(fā)動機工質的焓降大、功率密度大, 能達到活塞發(fā)動機無法實現(xiàn)的高航速和遠航程[1-2]。基于魚雷渦輪機的結構所限以及工作時長等特點, 渦輪盤轉子多設計為懸臂結構和滾動軸承支撐的結構形式。試驗數(shù)據(jù)表明, 在魚雷渦輪機的結構振動頻譜中, 常常出現(xiàn)軸承滾珠尺寸差及其橢圓度引起的振動分量[3]。同時, 軸承徑向游隙的存在會使軸承內(nèi)、外圈與滾動體在工作中碰撞產(chǎn)生沖擊, 導致出現(xiàn)非線性振動。研究軸承徑向游隙引發(fā)的魚雷渦輪轉子系統(tǒng)運動穩(wěn)定性問題具有重要意義。

針對軸承振動非線性機理, 采用關聯(lián)維數(shù)、最大Lyapunov指數(shù)、Kolmogorov熵及相關復雜性測度等方法對振動的混沌特性進行分類, 能夠較好地解釋軸承振動的復雜運動現(xiàn)象[4-5]。

但目前相關的多參數(shù)研究大多是其他參數(shù)為定值, 只針對單一參數(shù)為變量的運動狀態(tài)分類研究。董文凱等[6]采用Floquet理論對給定軸承徑向間隙條件下的解的穩(wěn)定性進行了研究。甄滿等[7]采用了變轉速分岔圖研究了不對中-碰磨條件下滾動軸承-轉子系統(tǒng)非線性動力學行為。

文中在考慮軸承內(nèi)部游隙、非線性赫茲接觸剛度及轉子不平衡量的基礎上, 建立了懸臂轉子-滾動軸承系統(tǒng)的動力學模型。重點研究了該轉子在特有工況和典型結構下, 滾動軸承徑向內(nèi)部游隙雙參數(shù)變化的系統(tǒng)振動失穩(wěn)的參數(shù)區(qū)域。

1 數(shù)學模型

由于魚雷渦輪機的空間結構特點, 采用懸臂梁式單盤轉子結構, 其系統(tǒng)動力學模型如圖1所示。

1.1 有限元模型

發(fā)動機轉子的動力學方程為

式中:為總質量陣;為總陀螺陣;為總阻尼陣;為總剛度陣;(,)為各節(jié)點軸承力、不平衡力和其他外載荷陣;為轉子各節(jié)點的位移轉角向量, 即

圖1 滾動軸承支承的魚雷渦輪轉子模型

首先, 對轉子系統(tǒng)按照軸段關鍵位置進行劃分, 然后進行轉子系統(tǒng)的質量陣、阻尼陣、剛度陣、陀螺陣和載荷陣計算, 并按照節(jié)點的依次順序對各矩陣進行組裝, 獲得轉子系統(tǒng)的總質量陣、總陀螺陣、總阻尼陣、總載荷陣以及總剛度陣。

1.2 滾動軸承模型

建立滾動軸承模型示意如圖2所示。滾動軸承由于滾珠和內(nèi)外滾道的相互作用, 產(chǎn)生總剛度的周期性變化, 形成變?nèi)岫?varying compliance, VC)振動。VC振動頻率即為滾珠的通過頻率, 即

式中: 為轉子角速度; 為外滾道半徑; 為內(nèi)滾道半徑; 為滾珠個數(shù)。

得出滾動軸承產(chǎn)生的時變軸承力, 并將其作用到方程(1)的右項, 即

1.3 渦輪機轉子雙參數(shù)Lyapunov指數(shù)計算

Lyapunov指數(shù)可以定量的表征相平面中2條相鄰軌線間的距離隨時間的平均指數(shù)發(fā)散率。對于魚雷渦輪機轉子, 定義節(jié)點位移為連續(xù)系統(tǒng), 即

2 計算與分析

采用4階龍格庫塔方法對動力學方程進行數(shù)值求解, 研究滾動軸承徑向內(nèi)部游隙隨轉速變化對轉子系統(tǒng)動力學特性的影響。

圖3 轉子圓盤節(jié)點x向位移分岔圖

由圖3可知, 伴隨著軸承徑向游隙從0增大到15μm, 振動位移的不穩(wěn)定域向低轉速區(qū)移動。從初始游隙時只存在擬周期狀態(tài), 到游隙增大后出現(xiàn)倍周期分岔及混沌, 軸承徑向游隙對魚雷渦輪轉子系統(tǒng)的運動穩(wěn)定性有著重要影響, 因此需在設計過程中合理控制系統(tǒng)中軸承的徑向游隙。

圖4 參數(shù)平面上的最大Lyapunov值分布圖

圖5 轉速4000 r/min, =8×10–7 m時向振動特性

當系統(tǒng)處于圖7的參數(shù)時, 由于不可約頻率及連續(xù)譜線的出現(xiàn), Poincare截面圖表現(xiàn)出云狀的離散點, 系統(tǒng)步入混沌。

圖6 轉速11 500 r/min, =1.1×10–6時x向振動特性圖

圖7 轉速10500 r/min, =1.65×10–6時向振動特性

3 結論

針對魚雷渦輪發(fā)動機, 在考慮軸承非線性赫茲接觸剛度、轉子不平衡量的基礎上, 建立了滾動軸承支撐的剛性水平轉子系統(tǒng)的動力學模型。通過雙參數(shù)平面的動力學分布圖, 研究分析轉子轉速、軸承內(nèi)部游隙對系統(tǒng)非線性動力學特性的影響規(guī)律。得出以下結論:

1) 轉速較低的情況下,VC所占的比例比較大, VC振動占主導地位, 轉速增大時, 系統(tǒng)不平衡激勵所導致的旋轉頻率的強迫振動分量增強;

2) 通過雙參數(shù)平面的動力學分布圖可以看出, 隨著徑向游隙增加, 進入混沌區(qū)域的最低轉速在不斷減低;

3) 由于概周期和混沌運動狀態(tài)下系統(tǒng)將出現(xiàn)連續(xù)的頻譜特征和漸變失穩(wěn)的時域特性, 說明在不同轉速下, 轉子的徑向游隙穩(wěn)定區(qū)域是不斷變化的, 因此在工程上應該盡量避免這樣的不穩(wěn)定參數(shù)條件。

[1] 查志武, 史小峰, 錢志博. 魚雷熱動力技術[M]. 北京: 國防工業(yè)出版社, 2006．

[2] 趙寅生, 錢志博. 魚雷渦輪機原理[M]. 西安: 西北工業(yè)大學出版社, 2002.

[3] 王儼凱, 王紅紅. 基于振動分析的魚雷渦輪機故障診斷技術初探[J]. 魚雷技術, 2009, 17(3): 34-38.

Wang Yan-kai, Wang Hong-hong. Fault Diagnosis for Torpedo Turbine Based on Vibration Analysis[J]. Torpedo Technology, 2009, 17(3): 34-38.

[4] 崔立, 鄭建榮. 柔性轉子滾動軸承系統(tǒng)混沌行為研究[J]. 中國機械工程, 2014, 25(3): 393-398.

Cui Li, Zheng Jian-rong. Study on Chaos Behavior of Flexible Rotor Rolling Bearing System[J]. Chinese Jou- rnal of Mechanical Engineering, 2014, 25(3): 393-398.

[5] 徐永智, 夏新濤, 南翔. 基于混沌理論滾動軸承振動穩(wěn)健化試驗數(shù)據(jù)的動態(tài)分析[J]. 航空動力學報, 2015, 30(8): 1958-1966.

Xu Yong-zhi, Xia Xin-tao, Nan Xiang. Dynamic Analysis of the Robust Test Data on Rolling Bearing Vibration Based on Chaos Theory[J]. Journal of Aerospace Power, 2015, 30(8): 1958-1966.

[6] 董文凱, 張智勇, 芮筱亭. 球軸承非線性動態(tài)剛度特性及其復雜共振行為研究[J]. 動力學與控制學報, 2018, 16(6): 483-487.

Dong Wen-kai, Zhang Zhi-yong, Rui Xiao-ting. Nonlinear Dynamic Stiffness Characteristics and Complicated Resonance Behaviors of Ball Bearings[J]. Journal of Dyna- mics and Control, 2018, 16(6): 483-487

[7] 甄滿, 孫濤, 田擁勝, 等. 滾動軸承轉子系統(tǒng)不對中-碰摩耦合故障非線性動力學分析[J]. 振動與沖擊, 2020, 7: 140-147.

Zhen Man, Sun Tao, Tian Yong-sheng, et al. Nonlinear Dynamic Analysis of Misalignment-Rub-impact Coupling Fault ofRolling Bearing Rotor System[J]. Journal of Vibration and Shock, 2020, 7: 140-147.

[8] Ivana Atanasovska. The Mathematical Phenomenological Mapping in Nonlinear Dynamics of Spur Gear Pair and Radial Ball Bearing due to the Variable Stiffness[J]. International Journal of Non-Linear Mechanics, 2015, 73(7): 114-120.

[9] Wahyu, Caesarendra, Kosasih B. Application of the Largest Lyapunov Exponent Algorithm for Feature Extraction in Low Speed Slew Bearing Condition Monitoring[J]. Mechanical Systems and Signal Processing, 2016, 50(1): 116-138.

[10] Hong W, Han Q, Luo R, et al. Dynamic Modeling of Moment Wheel Assemblies with Nonlinear Rolling Bearing Supports[J]. Journal of Sound and Vibration, 2017, 406: 124-145.

[11] Sharma A, Muniyappa A, Kankar P K. Nonlinear Dynamic Analysis of Defective Rolling Element Bearing Using Higuchi’s Fractal Dimension[J]. Sadhana, 2019, 44: 76.

[12] Aditya, Sharma, Nitin Upadhyay, et al. Nonlinear Dynamic Investigations on Rolling Element Bearings: A Review[J]. Advances in Mechanical Engineering, 2018, 10(3): 1-15.

Nonlinear Dynamic Characteristics Analysis of Torpedo Turbine Rotor-Rolling Bearing System

SONG Wen1, YANG Cheng-shi1, CHEN Zhi-wei2, SUN Yan2, REN Xin2, LI Gang2, CHAI Xiao-shuai2

(1. The 705 Research Institute, China State Shipbuilding Corporation Limited, Xi’an 710077, China; 2 ShanXi PingYang Industry Machinery Co.LTD, Linfen 043003, China)

To investigate the dynamic response of a torpedo turbine rotor, a dynamic model of a cantilever rotor-rolling bearing system is established, and the effect of the radial clearance of the bearing, nonlinear Hertz contact stiffness, and unbalance of the rotor are considered. The numerical simulation of the dynamic model of the system is carried out based on the variable step-size Runge-Kutta method. The maximum Lyapunov exponential spectra of the two-parameter plane are calculated, and the distribution of periodic and aperiodic motions of the rotor system with the change in the two parameters is obtained. The time-domain, frequency-domain, and periodic characteristics of the system under specific parameters are also studied. The analysis results show that the minimum speed in the chaotic region decreases with an increase in the bearing radial clearance. This study provides a useful reference for the design of torpedo-turbine systems.

torpedo; engine; turbine rotor; nonlinear dynamics; contact stiffness; two-parameter

TJ630.32; O194

A

2096-3920(2021)06-0690-05

10.11993/j.issn.2096-3920.2021.06.007

宋文, 楊赪石, 陳志偉, 等. 魚雷渦輪轉子-滾動軸承系統(tǒng)非線性動力學特性分析[J]. 水下無人系統(tǒng)學報, 2021, 29(6): 690-694.

2021-12-06;

2021-01-04.

國家自然科學基金項目(61403306).

宋 文(1986-), 男, 在讀博士, 高級工程師, 主要研究方向為結構動力學設計及振動控制.

(責任編輯: 許 妍)