相交轴转子系统当量化建模与扭振响应分析

陈立满; 杨大炼; 蒋玲莉; 李学军; Chen Liman; Yang Dalian; Jiang Lingli; Li Xuejun

网刊加载中。。。

使用Chrome浏览器效果最佳，继续浏览，你可能不会看到最佳的展示效果，

确定继续浏览么?

复制成功，请在其他浏览器进行阅读

相交轴转子系统当量化建模与扭振响应分析

- ORCID：
陈立满 ¹
- ORCID：
杨大炼 ¹
- ORCID：
蒋玲莉 ²
- ORCID：
李学军 ²

1. 湖南科技大学机械设备健康维护湖南省重点实验室，湘潭 411201； 2. 佛山科学技术学院机电工程学院，佛山 528225

最近更新：2021-03-08

DOI：10.6052/1672-6553-2020-043

目录contents

摘要

齿轮联结的相交轴转子系统具有转子不同体、不同速、不同轴线等结构特点，使其动力学建模异于常规转子系统，论文采用当量化建模方法，将相交轴转子系统转化为同转速、同轴线的转子系统并建立动力学模型，通过与已有研究对比验证了文中基于DyRoBeS当量化建模方法的有效性.以相交轴齿轮联结转子系统典型结构——某型直升机尾传动系统为研究对象，建立了当量化转子系统动力学模型，分析了系统各阶扭振临界转速及其组成、联结齿轮啮合刚度对扭振的影响、联结齿轮转动惯量对系统临界转速的影响，结果表明：传动系统第1、2阶临界转速为派生临界转速，其余8阶均源于各组成轴；中、尾减啮合刚度值与扭振的各阶共振峰值存在映射关系，部分啮合刚度值会激发共振极大值，应尽量避免；中、尾减啮合刚度变化会引发系统临界转速的变化，尤其是第1、2阶临界转速；不同阶临界转速对各个联结齿轮转动惯量的敏感性亦不同.论文从转子动力学角度研究了齿轮联结的相交轴转子系统的动力学特性，可为此类转子系统的结构优化设计及运行维护提供理论依据.

关键词

转子系统; 临界转速; 齿轮联结; 当量化建模; 扭振响应

引言

转轴及轴上旋转部件和支撑部件共同构成转子系统.最典型的转子系统是Jeffcott转子，此类转子系统的转轴保持同体、同速、同轴线等特性.由于结构、传动等方面的需求，转子系统可以是不同体的两组轴或多组轴.一种情况是通过联轴器联接，该类转子系统转轴不同体、但同转速、同轴线；另一种是通过齿轮联结，又可分为圆柱齿轮联结的平行轴转子系统和锥齿轮联结的相交轴转子系统，其中相交轴转子系统由螺旋锥齿轮联结，广泛应用于航空、航海等领域.

由于齿轮副的存在，转子系统将受到明显的啮合激励，啮合刚度、阻尼不仅使转子系统建模困难，而且对动力学特性产生极大影响，且不同类型的齿轮副产生的啮合激励差异较大.目前圆柱齿轮联结的平行轴转子系统动力学特性的研究较详尽.洪清泉等^［

1］在对传统的齿轮副扭转振动模型进行动力学等效变换的基础上，研究了齿轮传动转子系统的动态响应.蒋庆磊等^{［参考文献 2

百度学术}2］建立了齿轮传动多转子耦合系统动力学分析的通用模型，给出了建立多转子耦合振动模型的方法以及振动特性分析方法.张义民等^{［参考文献 3

百度学术}3］推导了三平行轴齿轮转子的啮合刚度矩阵，通过建立有限元模型，发现了系统不平衡响应会在非临界转速处进行积累，出现最大振动峰值的现象.Tang等^{［参考文献 4

百度学术}4］建立了一种新的简化扭转振动动力学模型，研究了复合行星混合动力系统的扭转振动特性.Zhang等^{［参考文献 5

百度学术}5］建立了五平行轴齿轮转子有限元模型，研究了转子系统在外载荷作用下，轴承刚度变化引起的临界转速和共振特性的变化.

直升机尾传动系统是典型的锥齿轮联结的相交轴转子系统，朱自冰等^［

6］、许兆棠等^{［参考文献 7

百度学术}7，8］使用集中质量和模态叠加法建立了直升机尾传动系统动力学模型，并获得传动系统的临界转速和模态振型.柴山，范磊等^{［参考文献 9

百度学术}9，10］推导了考虑弯扭耦合的齿轮啮合单元的传递矩阵，将传递矩阵法应用于研究存在啮合耦合的多级齿轮传动系统弯扭耦合振动.为了建立更精确的相交轴齿轮传动系统模型，许多学者将多种建模方法整合在一起处理由于螺旋锥齿轮的存在而产生的问题.Li等^{［参考文献 11

百度学术}11］首先建立离散的集中质量模型，然后根据锥齿轮啮合边界条件利用拉格朗日方程建立弯扭轴耦合动力学模型.虽然集中质量法计算速度快，计算效率高，但它通常使用有限的坐标，不能完全准确地反映传动系统的结构特性，Hua^{［参考文献 12

百度学术}12］将集中质量法和有限元法结合起来，提出一种改进的集中质量建模方法.锥齿轮联结的相交轴转子系统，是一类不同体、不同转速、不同轴线的转子系统，在空间中改变动力传递的方向，不同传动轴具有不同的局部坐标系，角位移、角速度和转动惯量等无法在全局坐标系下度量；尾传动系统的三组传动轴由两级减速器联结，两级减速器是螺旋锥齿轮副，产生的啮合激励更加复杂，且其装配精度更难控制，导致建模、临界转速和振动响应分析更加困难.国内外学者在研究齿轮转子系统的扭振动力学特性时，常对转子系统做当量化处理.王琪^{［参考文献 13

百度学术}13］、Chen等^{［参考文献 14

百度学术}14］提出了齿轮转子的当量化建模思想和理论.Wu等^{［参考文献 15

百度学术}15］推导了多分支齿轮转子系统当量化建模公式.Tang等^{［参考文献 16

百度学术}16］应用当量化建模方法建立了混合动力汽车传动系统扭振等效模型，Senjanović等^{［参考文献 17

百度学术}17］、Murawski等^{［参考文献 18

百度学术}18］将当量化建模方法应用到舰船推进系统建模中.

齿轮转子系统在工作过程中不可避免的产生扭振，因此进行传动系统扭振响应研究很有必要.Han等^［

19］通过参数化研究，对实际船舶发生高扭转振动的原因进行了分析.Zhao等^{［参考文献 20

百度学术}20］研究了具有两级行星齿轮级和一级平行齿轮级的风力发电机组齿轮箱的扭振.王飞鹏等^{［参考文献 21

百度学术}21］改进了基于扭矩激励的多支承转子系统力学模型，分析了无扭矩激励和暂态扭矩激励作用下，多支承转子系统的振动响应.廖鹏飞等^{［参考文献 22

百度学术}22］采用经典的频域集总参数多质点模型研究了冰载荷冲击下船舶推进轴系扭振响应.

本文以某型直升机尾传动系统作为研究对象，采用当量转化建模方法对相交轴转子系统做等效处理，进而使用DyRoBes软件建立相交轴齿轮转子系统有限元模型，从转子动力学角度研究相交轴齿轮转子系统扭振响应规律，可为齿轮联结相交轴转子系统的机构优化设计和运行维护提供参考.

1 齿轮联结相交轴转子系统有限元模型

论文以典型相交轴齿轮联结转子系统——直升机尾传动系统为研究对象，直升机尾传动系统由水平轴（HS）、中间减速器、尾斜轴（OS）、尾减速器和尾桨轴（TRS）组成，见图1.其中，中间减速器与尾减速器均由螺旋锥齿轮副构成^［

23］.表1给出了中、尾减螺旋锥齿轮的主要参数，其中三根组成轴的长度分别为

L_{H S} = 2400

mm，

L_{O S} = L_{T R S} = 800

mm，外径

D = 90

mm，内径

d = 80

mm，材料为结构钢，弾性模量

E = 206

GPa，

ρ = 7800

kg/m³.

图1 某型直升机尾传动系统三维实体模型

Fig.1 3D model of a helicopter tail drive system

表1 尾传动系统各齿轮主要参数

Table1 Main parameters of each gear in the tail drive system

	Input of system	Middle reducer		Tail reducer
Modulus	9	8		6
Number of teeth	34	32	40	23	64
Angle of shafts(°)	90	90		90
Tooth width(mm)	110	40		50

1.1　齿轮联结相交轴转子系统当量化有限元建模方法

有限元法的基本思想是将连续体看成有限个基本单元在节点处彼此相连接的组合体，使问题变成有限自由度的力学问题，从而借助线性方程组求解.对于扭振而言，有限元法将一个连续扭振问题离散化成有限自由度的振动问题，扭振系统的离散化力学模型包括质量元件、弹性元件和阻尼元件.对于螺旋锥齿轮联结，假设齿轮采用刚性支撑，啮合的螺旋锥齿轮副等效为由弹性单元和阻尼单元联结的等效圆盘，建立两自由度齿轮副动力学模型，见图2.

图2 螺旋锥齿轮副动力学模型

Fig.2 Dynamic model of spiral bevel gear pair

$I_{p}, I_{g}, T_{p}, T_{g}, θ_{p}, θ_{g}, R_{p}, R_{g}$ 分别是主从动齿轮的转动惯量、扭矩、扭转角位移和基圆半径， $C_{m}, K_{m}$ 是啮合阻尼和啮合刚度， $e (t)$ 为静态传动误差， $f (x)$ 为齿侧间隙函数，根据动力学模型，参照文献［

24］给出螺旋锥齿轮副扭振微分方程：

\begin{matrix} \begin{array}{l} I_{P} \overset{. .}{θ_{p}} + C_{m} R_{p} (\overset{\cdot}{R_{p} θ_{p}} - \dot{R_{g} θ_{g}} - \dot{e (t)}) + \\ K_{m} R_{p} f (R_{p} θ_{p} - R_{g} θ_{g} - e (t)) = T_{p} \end{array} \\ \begin{array}{l} I_{g} {\overset{. .}{θ}}_{g} - C_{m} R_{g} (\overset{\cdot}{R_{p} θ_{p}} - \dot{R_{g} θ_{g}} - \dot{e (t)}) + \\ K_{m} R_{g} f (R_{p} θ_{p} - R_{g} θ_{g} - e (t)) = - T_{g} \end{array} \end{matrix}

（1）

采用DyRoBes有限元软件，进行相交轴转子系统有限元建模.由于相交轴齿轮联结转子系统各级输入轴、联结齿轮、输出轴各部分的转速不同，因此采用当量化建模方法，以输入轴为等效构件，输出轴的惯性、弹性和角位移等使用水平轴上的等效量来表征.当量转化的原则是使当量化后的动能及弹性变形能与当量化前相同^［

12］.按照式（2）进行计算.

\{\begin{array}{l} \frac{1}{2} I_{e} ω_{e}^{2} = \frac{1}{2} I_{p} ω_{p}^{2} \\ \frac{1}{2} K_{e} (φ_{m} - φ_{m + 1})^{2} = \frac{1}{2} K_{p} (i φ_{m} - i φ_{m + 1})^{2} \end{array}

（2）

式中， $I_{p}, I_{e}$ 分别为原始转动惯量和当量转动惯量， $ω_{p}, ω_{e}$ 为原始角速度和当量角速度， $K_{p}, K_{e}$ 分别为原始刚度和当量刚度， $φ_{m}, φ_{m + 1}$ 分别为第 $m$ 和 $m + 1$ 质量的扭转角， $i$ 是速比， $i = \frac{ω_{e}}{ω_{p}}$

当量转化后相交轴转子系统解决了输入、输出轴不同轴线、转速等问题，为建立有限元模型创造条件.根据已建立的螺旋锥齿轮副动力学模型，结合当量转化建模思想，建立的一级相交轴齿轮转子系统当量化模型如图3所示，图中 $K_{m}$ 是以扭转刚度表示的综合啮合刚度， $C_{m}$ 是啮合阻尼， $I_{e}$ 表示当量化后从动轮的转动惯量.相交轴转子系统转化为了同速、同轴线转子系统，为从转子动力学角度研究相交轴齿轮传动系统动力学特性提供了便利.

图3 相交轴齿轮转子系统当量化模型

Fig.3 Idealized model of intersecting shaft rotor system

为验证本文基于DyRoBes的相交轴齿轮联结的转子系统当量化建模方法的有效性，与文献［

8］基于解析法建立的相交轴多级齿轮传动系统分析案列进行对比.根据文献［8］提供的数据，建立基于DyRoBes的当量化有限元模型，计算各阶扭振临界转速，并与文献［8］的分析结果进行对比，如表2所示.表中可见，两种方法误差不大，尤其是低阶误差.由此说明本文提出的基于DyRoBes的齿轮联结的相交轴转子系统当量化建模方法的有效性，可用于直升机尾传动系统动力学建模分析.

表2 文中方法与现有文献临界转速分析结果对比

Table 2 Comparison between the method and the results of natural frequency analysis in existing literatures

	Result in this paper (rad/s)	Result in[8](rad/s)	Error(%)
1^st order	0	0	0
2^ndorder	108	109	0.92
3^rd order	319	316	0.94
4^th order	2723	2285	17.49
5^th order	4776	4514	5.64

1.2　典型相交轴转子系统直升机尾传动系统有限元模型

应用基于DyRoBes的相交轴齿轮联结的转子系统当量转化建模方法，对直升机尾传动系统进行简化处理，把系统各部分向尾传动力输入轴处等效，以水平传动轴为等效构件，将尾斜轴的惯性、弹性、力矩和角位移都用等效构件上的等效量来代替，根据表1中的数据和图1所示的轴系基本结构，采用当量化建模方法，建立图4所示的尾传动系统等效扭振有限元模型，为便于分析，设置29个节点，关键节点位置见表3.将轴段离散化成n个微元轴段，有限元模型具有n个扭转自由度，其广义坐标为：

Θ = {θ_{1}, θ_{2}, θ_{3} \dots θ_{27}, θ_{28}, θ_{29}}

（3）

其中 $θ_{i}$ 为第i微元处的扭转角位移.根据转子动力学理论，尾传动系统的扭振微分方程为：

I \overset{. .}{Θ} + C \dot{Θ} + K Θ = Q

（4）

其中， $I$ 为转动惯量矩阵， $C$ 为扭转阻尼矩阵， $K$ 为扭转刚度矩阵， $Q$ 为扭矩激励矩阵.

图4 尾传动系统等效扭转振动模型

Fig.4 Equivalent torsional vibration model of tail drive system

表3 尾传动系统当量化模型关键节点位置

Table 3 Important nodes of idealized model of tail drive system

Node	Node 1	Node2/7/12	Node 17	Node 18	Node 23	Node 24	Node 29
Location	Input gear	Bearings	Input gear of middle reducer	Ouput gear of middle reducer	Input gear of tail reducer	Ouput gear of tail reducer	Tail rotor

2 齿轮联结的相交轴转子系统扭振特性分析

2.1　扭振临界转速和扭振模态

直升机尾传动系统由三根传动轴和联结转轴的两级减速器组成，齿轮联结副使系统存在明显的啮合耦合，将系统级和各组成轴前五阶临界转速进行对比分析，如表4所示，转子系统的临界转速远远低于各组成轴.进而求解了尾传动系统的前10阶临界转速，如表5所示，分析发现，系统第1、2阶临界转速是由耦合作用产生，其余8阶系统级临界转速均源于各组成轴，其1、2阶扭振模态振型见图5，由图5可知，尾斜轴和尾桨轴的扭振幅值较大，需重点关注.

表4 前5阶扭振临界转速对比

Table 4 Comparison of critical speed of first 5 orders torsional vibration

	1^st order (r/min)	2^nd order (r/min)	3^rd order (r/min)	4^th order (r/min)	5^th order (r/min)
TRS	7627	120881	253766	104352	776594
OS	41369	139513	263865	405581	736241
HS	6716	31247	59621	91479	125834
The transmission system	774	1721	6898	7947	31251

表5 前10阶系统级临界转速

Table 5 The first 10 orders critical speed of tail drive system

Order	Speed(R/Min)	Source	Order	Speed(R/Min)	Source
1	774	Couple	6	41065	OS(1)
2	1721	Couple	7	59622	HS(3)
3	6898	HS(1)	8	91497	HS(4)
4	7947	TRS(1)	9	120906	TRS(2)
5	31251	HS(2)	10	125834	HS(5)

图5 前两阶扭振模态振型

Fig.5 The first two torsional modes

2.2　联结齿轮啮合刚度对扭振的影响研究

相交轴转子系统的激励源主要来自于联结齿轮副，齿轮副产生啮合激励将严重影响传动系统的工作状态，造成断轴、断齿等严重机械故障.

齿轮啮合过程中，由于参与轮齿啮合的齿对数是周期性变化的，因此，啮合过程中的刚度也是周期性变化的，轮齿的时变刚度可以看作啮合线方向的时变弹性元件产生周期性变化.这种因啮合刚度的时变性产生的动态啮合力并对齿轮传动系统进行激励的现象就是刚度激励.时变啮合刚度采用式（5）计算：

k_{n} (t) = \frac{F_{n} (t) / B}{δ (t)}

（5）

其中， $k_{n} (t)$ 为时变啮合刚度， $F_{n} (t)$ 为法向动态啮合力，B为螺旋锥齿轮齿宽， $δ (t)$ 为齿面动态变形量.由于螺旋锥齿轮齿形复杂，传统解析法不适用于计算时变啮合刚度，因此使用Kisssoft软件计算时变啮合刚度，为简化计算，使用平均啮合刚度代替时变啮合刚度，同时忽略啮合阻尼对临界转速的影响.

基于已建立的有限元模型，研究相交轴转子系统的扭振响应，中、尾减平均啮合刚度为 $k_{1} = 1.6 \times 10^{6}$ N/m， $k_{2} = 3 \times 10^{6}$ N/m，

图6给出各节点处的扭振频谱，图中存在四个明显的共振转速，其数值分别对应于系统前四阶临界转速，各阶临界转速处不同节点的扭振角位移峰值存在显著差异.

图6 各节点扭振频谱图

Fig.6 Spectrum of angular displacement at all nodes

进一步研究前四阶共振角位移随中、尾减齿轮副啮合刚度的变化规律.中减齿轮副啮合刚度 $k_{1}$ 取0.6×10⁶N/m，依次递增0.2×10⁶N/m，至2.6×10⁶N/m，尾减齿轮副啮合刚度k₂=3.0×10⁶N/m保持不变，分析关键节点17、18和29，即中减主动齿轮、中减从动齿轮和尾旋翼，扭转角位移随啮合刚度改变的变化规律，如图7所示.

图7 不同中减啮合刚度下的代表性节点各阶共振峰值17

Fig. 7 The resonant peaks of the representative nodes with different middle meshing stiffness

(a)节点 (b)节点18 (c)节点29

(a)node17 (b)node18 (c)node29

图7中，第1、2阶共振角位移对中减齿轮副啮合刚度变化响应丰富，中减啮合刚度取 $k_{1} = 1.0 \times 10^{6}$ N/m， $2.0 \times 10^{6}$ N/m时，各节点一阶共振角位移出现极大值；中减啮合刚度取 $k_{1} = 1.4 \times 10^{6}$ N/m时，节点18、29的二阶共振角位移出现极大值.由此说明，在系统结构设计时，上述中减啮合刚度值要尽量避免；在系统运行监测中，扭振的变大有可能是啮合刚度变化引发，为系统结构优化设计及运行维护提供依据.第四阶共振峰值对啮合刚度变化不敏感，但节点29的四阶共振峰值偏大.

尾减啮合刚度 $K_{2}$ 取 $2.0 \times 10^{6}$ N/m，依次递增 $0.2 \times 10^{6}$ N/m，至 $4.0 \times 10^{6}$ N/m，中减齿轮副啮合刚度 $K_{1} = 1.6 \times 10^{6}$ N/m保持不变.取关键节点23、24和29，即尾减主动齿轮、尾减从动齿轮和尾旋翼，分析扭转角位移随啮合刚度改变的变化规律，如图8所示.

图8 不同尾减啮合刚度下的代表性节点各阶共振峰值

Fig. 8 The resonance peaks of the representative nodes with different tail meshing stiffness

(a)节点23 (b)节点24 (c)节点29

(a)node23 (b)node24 (c)node29

图8中，第一、二、四阶共振扭转角位移对尾减齿轮副啮合刚度变化响应丰富，尾减啮合刚度 $K_{2} = 3.8 \times 10^{6}$ N/m时，各节点一阶共振角位移出现极大值；尾减啮合刚度 $K_{2} = 3.0 \times 10^{6}$ N/m时，各节点二阶共振角位移出现极大值；尾减啮合刚度 $K_{2} = 2.6 \times 10^{6}$ N/m时，各节点四阶共振角位移出现极大值；在系统结构设计时，要尽量避免尾减的上述啮合刚度值.相较于图7，尾减齿轮副啮合刚度变化对各节点各阶共振的影响远大于中减齿轮副，在系统结构设计时，尾减的啮合刚度选用要更加谨慎.

中、尾减齿轮副啮合刚度变化对各阶临界转速同样会产生影响，图9描述了系统前四阶临界转速随中、尾减齿轮副啮合刚度变化的规律.由图可知，第一、二阶临界转速对啮合刚度变化较敏感；尾减速器齿轮副啮合刚度的变化对临界转速的影响大于中减速器.

图9 系统前四阶临界转速随联结齿轮副变化趋势图

Fig.9 The trend chart of the first four critical speed of the system with the coupling gear pair

(a)中减齿轮副啮合刚度变化 (b)尾减齿轮副啮合刚度变化

(a)meshing stiffness of middle gear pair (b)meshing stiffness of tail gear pair

2.3　联结齿轮转动惯量对系统扭振临界转速的影响

在尾传动系统中，两级减速器中的螺旋锥齿轮具有较大的转动惯量，对传动系统的扭振临界转速产生重要影响，将各部件的转动惯量按比例增减，令 $I_{i} = (δ I_{p})$ kg·m⁴，（ $i$ =1，2，3，4，）分别代表中减主动轮（node17）、从动轮（node18）、尾减主动轮（node23）、从动轮（node24）； $I_{p}$ 为各部件的原始转动惯量， ${I_{p}}_{c}$ 为各部件改变后的转动惯量； $δ = ({I_{p}}_{c} / I_{p}) \times 100 %$ =25%，50%，75%，100%，125%，150%，175%，200%.为便于表现各阶临界转速随各部分（节点）转动惯量变化，以各 $δ$ 值对应的临界转速与 $δ = 100 %$ 时临界转速的差值为纵坐标，以表征转动惯量变化的 $δ$ 值为横坐标，绘制尾传动系统各阶扭振临界转速对转动惯量变化的响应曲线.

由图10可知，尾传动系统扭振临界转速对不同部件转动惯量变化的敏感性不同，尾减速器主动轮（node23）对系统级临界转速的影响很小；尾传动系统不同阶临界转速对同一部件转动惯量变化的敏感性亦不相同.中减主动轮（node17）转动惯量的变化对尾传动系统第三、四、五阶临界转速影响较大.

图10 系统前四阶临界转速随转动惯量变化趋势图

Fig.10 The trend chart of the first four critical rotational speeds of the system with the moment of inertia

3 结论

（1）本文提出一种基于Dyrobes的当量化建模方法建立齿轮联结相交轴转子系统的有限元模型，与已有研究成果对比验证了本方法的正确性.应用该方法建立了典型齿轮联结相交轴传动系统——直升机尾传动系统的有限元模型.

（2）基于已建立的直升机尾传动系统有限元模型，获得了传动系统扭转振动临界转速和模态振型；对比分析了各组成轴和传动系统的扭振临界转速，发现第1、2阶为系统派生临界转速，其余各阶源于各组成轴，并且轴系的扭振临界转速较各组成轴显著降低.

（3）中、尾减啮合刚度值与扭振的各阶角位移峰值存在映射关系，部分啮合刚度值会激发角位移极大值，保持尾减啮合刚度 $k_{2}$ 不变，取 $k_{1} = 1.0 \times 10^{6}$ ， $1.4 \times 10^{6}$ ， $2.0 \times 10^{6}$ N/m时产生共振角位移极大值；保持中减啮合刚度 $k_{1}$ 不变，取 $k_{2} = 2.6 \times 10^{6}, 3.0 \times 10^{6}, 3.8 \times 10^{6}$ N/m时产生共振角位移极大值；在系统结构设计及运行中应尽量避免上述啮合刚度值.中、尾减啮合刚度变化会引发系统临界转速的变化，尤其是第1、2阶临界转速.

（4）系统级临界转速对不同部件转动惯量改变的敏感性不同，中间减速器主动轮对系统级临界转速影响最大，尾减速器主动轮对系统级临界转速的影响最小；系统级各阶临界转速对同一部件转动惯量改变的敏感性也不同.

参考文献

洪清泉，程颖. 基于ADAMS的多级齿轮传动系统动力学仿真. 北京理工大学学报， 2003（6）：690~693（ Hong Q Q，Cheng Y. Dynamic simulation of multistage gear train system in ADAMS. Journal of Beijing Institute of Technology， 2003（6）：690~693（in Chinese）） [百度学术]

蒋庆磊，吴大传，谭善光，等.齿轮传动多转子耦合系统振动特性研究.振动工程学报，2010，23（3）：254~259 [百度学术]

Jiang Q L，Wu D Z，Tan S G，et al. Development and application of a model for coupling geared rotors system. Journal of Vibration Engineering， 2010，23（3）：254~259（in Chinese） [百度学术]

张义民，何永慧，朱丽莎，等.多平行轴齿轮耦合转子系统的振动响应. 振动、测试与诊断， 2012，32（4）：527~531 [百度学术]

Zhang Y M，He Y H，Zhu L S，et al. Vibration response of multi-shaft rotor system with coupled gear mesh. Journal of Vibration， Measurement & Diagnosis，2012，32（4）：527~531（in Chinese） [百度学术]

Tang X， Yang W， Hu X， et al. A novel simplified model for torsional vibration analysis of a series-parallel hybrid electric vehicle. Mechanical Systems & Signal Processing， 2017， 85：329~338 [百度学术]

Zhang H ， Zhai J ， Han Q ， et al. Dynamics of a geared parallel-rotor system subjected to changing oil-bearing stiffness due to external loads. Finite Elements in Analysis and Design， 2015， 106：32~40 [百度学术]

朱自冰，朱如鹏，鲍和云.直升机尾传动系统扭转振动建模与特性.航空动力学报，2013，28（2）：432-438 [百度学术]

Zhu Z B， Zhu R P，Bao H Y. Torsional vibration characteristics analyzing and modeling for helicopter tail drive system. Journal of Aerospace Power，2013，28（2）：432~438（in Chinese） [百度学术]

许兆棠.直升机传动系统扭转振动的分析.工程力学，2012，29（9）：330~336 [百度学术]

Xu Z T. Analysis of torsional vibration of helicopter transmission system. Engineering Mechanics，2012，29（9）：330~336（in Chinese） [百度学术]

许兆棠，朱如鹏.直升机尾传动系扭转振动的分析.航空学报，2007，2：425~431（Xu Zhaotang， Zhu Rupeng. Analysis of torsional vibration of helicopter tail drive system. Acta Aeronauta Sinica，2007，2：425~431（in Chinese）） [百度学术]

柴山，高连勇. 多转子系统弯扭耦合振动分析的整体传递矩阵法. 机械设计， 2005， 22（10）：8~10 [百度学术]

Chai S， Gao L Y. Integral transfer matrix method for flexural-torsional coupling vibration analysis of multi-rotor systems . Journal of Mechanical Design， 2005，22（10）：8~10 （in Chinese） [百度学术]

范磊，王少萍. 高速内啮合齿轮轴减速器扭转振动传递矩阵法分析. 液压与气动， 2015，290（10）：18~22 [百度学术]

Fan L，Wang S P. Torsional vibration analysis of high speed internal gear shaft speed reducer by the transfer matrix method. Chinese Hydraulics &Pneumatics，2015，290（10）：18~22（in Chinese） [百度学术]

Li M， Hu H Y， Jiang P L， et al. Coupled axial-lateral-torsional dynamics of a rotor-bearing system geared by spur bevel gears. Journal of Sound and Vibration， 2002， 254（3）： 427~446 [百度学术]

Hua X， Lim T C ，Peng T ， e al. Dynamic analysis of spiral bevel geared rotor systems applying finite elements and enhanced lumped parameters. International Journal of Automotive Technology， 2012， 13（1）：97~107 [百度学术]

王祺. 内燃机轴系扭转振动. 北京：国防工业出版社，1984 （Wang Q. Torsional vibration of internal combustion engine shafting. Beijn g：National Defense Industry Press，1984（in Chinese）） [百度学术]

Chen W J， Gunter E J. Introduction to dynamics of rotor-bearing systems. Victoria， BC： Trafford， 2007 [百度学术]

Wu J S， Chen C H. Torsional vibration analysis of gear-branched systems by finite element method. Journal of Sound and Vibration， 2001， 240（1）： 159~182 [百度学术]

Tang X， Yang W， Hu X， et al. A novel simplified model for torsional vibration analysis of a series-parallel hybrid electric vehicle. Mechanical Systems and Signal Processing， 2017， 85： 329~338 [百度学术]

Senjanović I， Hadžić N， Murawski L， et al. Analytical procedures for torsional vibration analysis of ship power transmission system. Engineering Structures， 2019， 178： 227~244 [百度学术]

Murawski L， Charchalis A. Simplified method of torsional vibration calculation of marine power transmission system. Marine Structures， 2014， 39： 335~349 [百度学术]

Han H S， Lee K H， Park S H. Parametric study to identify the cause of high torsional vibration of the propulsion shaft in the ship. Engineering Failure Analysis， 2015：S1350630715301187 [百度学术]

Zhao M， Ji J C. Nonlinear torsional vibrations of a wind turbine gearbox. Applied Mathematical Modelling， 2015， 39（16）：4928~4950 [百度学术]

王飞鹏，杨兆建，宋高峰，等. 一种改进的基于扭矩激励的多支承转子系统模型. 机械科学与技术， 2014， 33（5）： 697~701 [百度学术]

Wang F P，Yang Z J，Song G F，et al. An Improved multi-support rotor system with torque excitation. Mechanical Science and Technology for Aerospace Engineering， 2014， 33（5）： 697~701（in Chinese） [百度学术]

廖鹏飞，周瑞平，李江. 冰区船舶电力推进轴系瞬态扭振计算研究. 推进技术， 2018， 39（8）： 1889~1896 [百度学术]

Liao P F，Zhou R P，Li J， et al. Research on calculation of transient torsional vibration of electric propulsion shafting in ice area. Journal of Propulsion Technology， 2018， 39（8）： 1889~1896（in Chinese） [百度学术]

蒋玲莉，华登荣，李学军，等.某型直升机传动系统动力学虚拟仿真研究.机械传动，2016，40（12）：12~16 [百度学术]

Jiang L L，Hua D R，Li X J，et al. Study on the dynamics virtual simulation of the transmission system of a type helicopter. Journal of Mechanical Transmission， 2016，40（12）：12~16（in Chinese） [百度学术]

李润方，王建军. 齿轮系统动力学—振动、冲击、噪声.北京：科学出版社，1997（F， Wang J J. Dynamics of gear systems-vibration， impact， noise. Beijing：Science Press， 1997（in Chinese）） [百度学术]

相交轴转子系统当量化建模与扭振响应分析

摘要

关键词

引言