含磨损故障的齿轮传动系统非线性动力学特性

本资料包含pdf文件1个，下载需要1积分

Nonlinear dynamics analysis of gear transmission system with wear faultWANG Xiao sun ， Shi-jing ，ZHOU Xu-hui ，HU Ji-cai(1.School of Power and Mechanical Engineering，Wuhan University，Wuhan 430072，China；2.Wuhan 2 Shipbuilding Institute，Wuhan 430064，China)Abstract： The Achard and W eber-Banaschek formula were employed to calculate the dynamic surface wear and themesh stifnesses of worn surface gear pair respectively. The three-degree-of-freedom translational-rotational couplednonlinear dynamic equations of gear transmission system，including the wom teeth pair's time-varying mesh stifness，piece-wise backlash and internal error excitation，were presented to do in-depth investigation of the surface wearS effect on geartransmissionS nonlinear vibration characteristics. The varying step GILL integration method was employed to performnumerical simulation of the dynamic mode1.The internal excitation frequency was selected as the parameter to calculatethe bifurcation diagram.The orthonormalization treatment of the systemS Jacobie matrix was caried out by utilizingGRAM-SCHMIDT method in order to obtain the Lyapunov exponents.The Poincar6 section and power spectrum were usedto verify the results of bifurcation and Lyapunov exponents of the system under some special parameter settings.Under theinfluence of systemS strong nonlinearities，a rich variety of bifurcation phenomena were illustrated.The classic periodic-doubling routes to chaos，intermitent routes to chaos and abundant different quasi-periodic routes were revealed bybifurcation diagram and Lyapunov exponents. One ordinary periodic-doubling route and two particular quasi-periodicroutes were demonstrated in detail with the aid of Poincar6 maps plotted in the phase plane.Ahernant quasi-period andphase-locking were observed in the systemS quasi-route to chaos.In addition，it was observed that the frequencies ofquasi-periodic motion satisfy the familiar Farey sequence.All the results indicate that the dynamic characteristics of geartransmission with wear fault are very complex and the systemS routes to chaos are abundant and diverse。

Key words：bifurcation；chaos；wear fault；mesh stiffness；gear transmission system齿轮传动系统是目前机械系统中使用最为广泛的基金项 目：国家自然科学基金资助项目(51005170)；武汉市学科带头人计划(20105I73O545)收稿 日期：2012-08-14 修改稿收到日期：2012-09-26第-作者 王晓笋 男 ，博士，副教授，1979年生传动方式，其通过轮齿啮合实现载荷和运动的传递 11。

而由于其内部存在的齿侧间隙、时变啮合刚度和内部误差激励等非线性因素，使得齿轮传动系统的负载运行的过程中呈现复杂的非线性动力学特征。Kahraman等-zI4 考虑了齿侧间隙和轴承间隙等非线性因素，建立单级齿轮传动系统单自由度、二自由度和三自由度振 动 与 冲 击 2013年第 32卷非线性动力学模型，应用解析法和数值法对模型进行求解，发现系统存在的跳跃不连续和混沌等典型非线性现象。Theodossiades等 和 Raghothama等。 分别采用多尺度法和增量谐波平衡法对齿轮传动系统非线性动力学模型进行了求解，发现了系统存在的亚谐波共振和超谐波共振现象。刘梦军等 。采用数值积分的方法研究了单级齿轮传动系统所具有的分岔与混沌现象。Kahraman等 和 Parker等 分别建立了齿轮传动系统非线性动力学试验系统，通过试验测试获得的幅频特性曲线验证了系统存在的跳跃不连续现象。Cai等 考虑了齿轮传动系统的非线性支撑力，通过数值积分的方法研究了啮合阻尼和激励频率对于齿轮传动系统振动特性的影响。

而轮齿在运行的过程中由于受到载荷和齿面润滑状态的影响，齿面会出现不间断的磨损现象。齿面磨损的直接结果是导致齿面金属材料的非均匀去除 ，但同时也会进-步引起齿轮传动系统振动与噪声的加剧 J，由此可能使得轮齿的动载荷急剧增加，并由此可能引起轮齿受到更大的冲击载荷作用，最终将导致轮齿折断。分析含有磨损故障特征的齿轮传动系统非线性动力学特性是揭示齿轮传动系统在出现磨损之后会出现早期失效的有效途径之-。为此，本文考虑了齿面磨损可能引起的直齿轮传动系统的特性参数的变化，建立了含有齿面磨损故障的齿轮传动系统非线性动力学模型并对其进行数值特性分析，以变步长积分方法计算获取了系统的分岔图和对应的李雅普诺夫指数谱，由此分析系统在某个关键参数变化时的全局非线性动力学特性。同时，通过 Poincar6映射、功率谱对系统出现的周期、拟周期和混沌运动进行深入细致地研究，揭示系统由周期经分岔进入混沌的途径以及非线性系统内部所蕴含的-些普适性规律。

1 三自由度动力学的模型- 对直齿轮的三自由度平移 -扭转耦合的动力学模型如图 1所示，模型中0表示角位移，Y表示竖直方向位移，k表示刚度，C表示阻尼，表示转动惯量，m表示质量，R为基圆半径， 表示扭矩，F表示力，e表示内部误差激励。下标 gl， ，61，62分别代表主从动齿轮和主从动齿轮的支撑轴承，h表示啮合齿轮对，eq代表等效。符号-”表示带有量纲的变量。为消除系统内存在的刚体位移， jI入轮齿之间动静态传递误差之差 Rl 0 1-Rg20 Y 61-Y62-e( )，由此建立齿轮传动系统三自由度动力学方程如式(1)所示。

m 1-)，t6lc6l歹 lc6 6 ( 61) (了 (6)-F 1- - Im ：2c62-y 612-Ch占 kbzfb2( 62)-k ( ) (6)F62 (1)- ～ - lm (6 --yt61l )Ch占 (了 ( )F - ”(了)ch22图 1 三 自由度 直齿 轮传 动系统动力学模型Fig.1 Three-degree-of-freedom dynamic modelof a spur gear pair transmission system其中：m叼，g /( R 尺 )， ( )Esin( h)，F 1/R /Rg2。

引入标称尺度 b。和时间尺度 tm t，无量纲广义坐标表示为Y 6 /b (i1，2)，6：6/b ，可以得到方程式(1)的无量纲形式如式(2)所示。

l2( 1ly l(136 )kl 1(y62)(t)fh(6)-Fy 2( 2 y 2- 3占 )k2 (Y6 J)-：。(t)fh(6)Fy 1-y：2 ”2 36 3(t)fh(6)F F H(2)其中： √ /m ， / /m (i1，2)， Cbi/2m 09 (i1，2)， bCh/2m ∞ (i1，2)， 3Ch/2m 。 ， (t)k(t)/m (cJ (i1，2)， ∞2 /(i1，2)，F F /m bc(J：(i1，2)，F ：F /m6 ：，F H Esin(Dt )/b。， ， 3(t)k (t)/k k 为平均啮合刚度。

Kahraman等 通过研究发现，轴承的间隙与齿轮的齿侧间隙均可以用分段非线性函数描述，其中主动轴轴承和从动轴轴承的间隙可以描述为：，)，6 -bbi/b Y > bbi/b。(Y6 ) 0 I Y l≤bbi/b i1，2 (3)y 66 /b Y6 <-6/b齿侧间隙可以描述为：r - 6 /b > bh/b(6)： 0 ≤bh/b (4)L 6/b 6 <-bh/b第 16期 王晓笋等：含磨损故障的齿轮传动系统非线性动力学特性 392 齿面磨损故障特征的引入2.1 齿面累积磨损量的计算为研究齿面磨损对于系统动力学特性的影响，首先采用 Archard公式1 计算齿面动态磨损特性，主要考虑轮齿啮合的动载荷、啮合点滑移速度等计算轮齿负载运行-定周期之后的累积磨损量，其中主从动齿轮齿面的累积磨损量可以描述为：，： ， -l2kp0 l(M- 2)/u ( 1，2)(5)式中：u，( 1，2)为主从动齿轮的滚动速度， 和P分别为当前啮合点的动态磨损系数和最大接触力， 为当前啮合点处 Hertz接触宽度。

2.2 轮齿啮合刚度的计算根 据 Weber.Banaschek公式m ，单个轮齿 中线和啮合线的交点(图2中P点)在啮合线方向的变形量包括弯曲剪切变形量 ，基础部分变形量 和皇 载 作图2轮齿的几何参数 用点到P点的接近量。 ， Fi二和 可以分别采用式(6a)，。 m e。ters ofGear To。th(6b)和(6c)计算获龋-对轮齿啮合时，轮齿在啮合线方向总体变形量与啮合刚度分别用式(6e)和(6f)表示。

FNs2卢 923.1(10.294tan Jf dy) (6a)FNs2晔2筹 1.4(10.294tan2 ) (6b)-o.58 ( n 2 hl -0.429)(6c)- 2√6 6 6Z2 尺1 艘 6。 (6e)kFⅣ/3 (6f)式(6(a)～(d))各参数如图 2所示，E为等效弹性模量，b为齿宽，P 和p：为接触点的当量圆柱体半径。

齿面累积磨损造成渐开线齿廓的变化，为此在式 6((a)-(C))中引入从动齿轮的齿面累积磨损量 H 和H2，由此计算出现磨损之后轮齿的啮合刚度。

2.3 算例取主动齿轮齿数z。30，从动齿轮齿数 z：40，模数 m4，齿宽 b0.03，标准压力角 20。，输入力矩 60 N·m，输出负载为 ：80 N·m，主从动齿轮弹性模量 EE2207 GPa，a1.33×10 (m /N)，磨损系数 k5×10 。(m /N) lZ]，通过式(5)计算获得轮齿运行 l0 次之后的累积磨损量沿啮合线的分布如图3所示，在节点附近的单齿啮合区间内磨损量最小，因为单齿啮合区间内虽然载荷由单个轮齿承担，但滑移速度较小，因而最终的累积磨损量相对较小；在靠近齿顶处的双齿啮合区间 1和靠近齿根处的双齿啮合区间2的磨损量则相对较大，特别在齿顶处累积磨损量是整个轮齿上最大的。

呈40翟zo盥 O图 3 啮合齿轮对齿面累积磨损量Fig.3 Accumulated surface wear of gear pair采用式(6(a)-(f))计算获得齿轮对啮合刚度如图4所示，齿轮对啮合刚度与磨损量-样分为双齿啮合区间1、单齿啮合区间和双齿啮合区间2〖虑啮合齿轮对长期负载运行之后出现较为严重的磨损，对比未磨损的齿轮对啮合刚度和出现磨损之后的啮合刚度可以发现，沿渐开线齿廓不均匀分布(如图4所示)磨损量使得单个轮齿的啮合刚度在齿顶和齿根位置减小，但在节点处变化相对较校髓压力角0Degree)图4 未磨损与出现磨损的轮齿时变啮合刚度Fig.4 Time-varying mesh stifnessof worn and unworn gear pair在建立齿轮传动系统动力学方程时，时变啮合刚度通厨行-定的简化处理，将图4所示的啮合刚度以-系列三角函数和的形式进行描述，从而便于求解处理。本文采用 Raghothama等 提出的式(7)表达时变啮合刚度，由于齿面磨损引起的啮合刚度的变化通过分析发现主要体现在 在上。通过对图4的磨损齿轮对的啮合刚度进行傅里叶级数展开，可以得到-系列的s 在值 ，本文取 r3。

k3(t)k(t)/k 1-∑ cos(rOt) (7)r l3 非线性系统的数值特性分析系统的诸多非线性因素会引起系统可能出现周期、拟周期和混沌响应，为全面获得系统的动力学响应第 16期 王晓笋等：含磨损故障的齿轮传动系统非线性动力学特性 4l(a) 力1.4，周期 1运动 (b) 门1.358，周期 2运动 (C) n1.347 5，周期 4运动(d) 1.346 2，周期 8运动 (e) 力1.344 4，周期 16运动 (f) 力1.344，混沌运动图7 系统的功率谱Fig.7 The power spectrums系统从周期 1运动开始(图6(a))，经过不多的周期加倍，经历周期2(图6(b))，周期 4(图6(C))，周期8(图6(d))和周期 16(图6(e))，最终系统进入混沌运动(图6(f))。在 Poincar6映射上，周期 n运动对应这rt个不动的离散点，而功率谱则对应着在 1/n基频及m/n(m为正整数)处存在尖峰(图7(a)～(e))。当系统进入混沌运动之后，Poincar6映射上出现了奇怪吸引子(图6(f))，功率谱曲线则表现为具有噪声背景的连续谱，同时在 1/8及 m/8基频处存在尖峰，表明了系统访问奇怪吸引子 8个不 同区域 的严格周期性 (图7(f))。通过相轨线、Poincar6映射和功率谱的计算结果验证了分岔图与李雅普诺夫指数谱的正确性。

3.2 拟周期分岔途径通过拟周期分岔进入混沌是系统存在的主要途径之-，在分岔图上可以观察到多个经 由拟周期分岔途径进入混沌的区域，而各个不同区域其具体的分岔途径又存在差异，下面选取两个具有代表意义的拟周期途径进行分析。

3.2.1 拟周期分岔途径 1在区间/2∈[1.2，1.22]存在着-个典型的拟周期进入分岔的途径，当系统的激励频率 1.22时，系统处于周期2运动状态，如图8(a)所示，此时 Poinear6映射上是两个稳定的焦点，功率谱上对应m/2倍基频处存在着尖峰图9(a)。随着激励频率的变化，当 1.217 2时，稳定的焦点失去稳定性，并经由 Hopf分岔变成两个吸引环图8(b)，此时，除了在m/2倍基频处存在着尖峰，分析其频域特性，如图9(b)所示，系统的频域上存在1/7，2/9，3/11，9/23等-系列不可约频率分量，而通过研究发现这些频率分量之间存在着以下关系：2 1l 3 12 X l 9 18 X 19 - 7 211 - 7 2×223 - 7 8×2这表明在系统的拟周期频域区间的频率分量之问满足周期相加法则 16j(period-adding law)。继续改变激励频率，系统Poincar6映射上的吸引环呈现面积不断增大的趋势 (图 8(b))，再次观察功率谱 曲线 (图9(b))可以发现系统的响应上出现了更多的频率分量。

随后 Poincar6映射上的吸引环开始发生变形(图8(C))和环面上的自缠绕(图8(d))，系统的功率谱上出现了更多不可公约的特征频率(图9(C)～(d))。环面上的拟周期运动逐渐失去稳定性 ，此时稍有扰动就会过渡到另外-个周期解，拟周期运动中的两频率之比将从无理数变成与之接近的有理数，当力1.207现，系统出现了周期 8响应(图8(e))，频率谱上(图9(e))也出现了 m/8基频处出现了尖峰。最终通过多周期锁相进入混沌状态(图7(f))，系统的功率谱变成了具有噪声背景的连续谱(图 9(f))，对应在 Poincar6映射上也出现了具有分维数的奇怪吸引子。特别的，在功率谱上面对应 m/2基频处也存在着尖峰，代表系统周期性的访问 Poincar6映射上的两个混沌区间。这种经过周期运动-拟周期-锁相-混沌是系统进入混沌的常规途径之- 。

3.2.2 拟周期分岔途径 2在区间 ∈[1.8，1.82]也存在拟周期分岔进人混沌的途径，随着激励频率 的改变，Poincar6映射上首先出现的是多周期离散点(图 10(a))；此时系统处于周期 13运动状态，接着周期运动通过 Hopf分岔被拟周期运动替代，Poincar6映射出现了具有扭曲缠绕特性的环面(图 10(b))；不稳定 的拟周期运动再次进入P107多周期运动状态，即锁相(图10(c))；随后是进-步扭曲变形的 环面(图 10(d))；不稳定拟周期运动再次被 P94多周期运动所替代(图 10(e))；从图 10可以发现，系统经历了多次类似拟周期和多周期锁相的交替出现，即具有扭曲 环面的拟周期运动(图10(f))- - - P81多周期锁相(图 10(g))-扭曲 环面的拟周期42 振 动 与 冲 击 2013年第 32卷运动 (图 10(h))-P136多周期锁相(图 10(i))-扭曲 环面的拟周期运动(图10(j))-P55多周期锁相(图10(k))。最终系统在经历上述多次拟周期-多周期锁相之后，在 Poincar6映射出现了具有连续功率谱的占(a) n1.22，周期2运动混沌运动(图 10(1))。这种经历多次拟周期-多周期锁相的通向混沌的演化路径是非常规通路，说明了含有磨损故障齿轮传动系统经由拟周期运动通向混沌的复杂性和多样性。

(b) 1.4，拟周期运动 (c) 1.2172，拟周期运动(d) ：1.214，拟周期运动 (e) 1.207，周期 8运动 (f) 门1.204，混沌运动图8 系统的Poincar6映射(d) 1.21，拟周期运动(b)力1.2172，拟周期运动 (c) 1.214，拟周期运动(e)力1.207，周期 8运动 (f)力：1.204，混沌运动图9 系统的功率谱Fig.9 Power speetrus同时，还可以发现嵌套在拟周期运动区间内相邻的多周期运动之间周期数存在-个有趣的现象，相邻三个周期运动中间- -个的周期数等于第-个和第三个周期运动周期数之和，当力1.807 5，系统处于周期 l3运动(图 10(a))；而当 1.814 8，系统处于周期 94运动(图 10(e))，在这两个周期窗 口之间，当 力1.813 6，存在着-个周期 107运动区间(图10(C))，而1071394。同时，系统另外三个嵌套在拟周期运动区问的多周期窗口也满足上述规律，当 1.816 8，系统处于周期81运动(图 10(g))；而当 1.821，系统处于周期 55运动(图 10(k))，在这两个周期窗 口之间，当 1.818 8，存在着-个周期 136运动区间(图lo(i))，而 1368155。从分岔图和李雅普诺夫指数谱上还可以发现，在系统进入混沌运动之后，混沌振动区间内依然间或出现了大量的多周期窗口，这种交替出现的混沌、多周期运动最终导致了超级混沌运动的出现，对应的李雅普诺夫指数谱具有[，，-，-，-， - ]的形式。

第 16期 王晓笋等：含磨损故障的齿轮传动系统非线性动力学特性 43动(e) 1.814 8，周期94运动 (f) n1.816，拟周期运动 (g) n1.816 8，周期 81运动 (h) 1.818，拟周期运动(i) 力1.818 8，周期 136运动 (j) 力1.820 4，拟周期运动 (k) 1.821，周期 55运动 (1) 力1.821 4，混沌图 10 系统的Poincar6映射Fig.10 Poincar6 maps4 结 论本文结合 Archard和 Weber-Banaschek公式分别研究了齿面动态磨损特性和由此引起的轮齿时变啮合刚度的动态变化，从而在齿轮传动系统的动力学方程中可以引入齿面磨损故障特性。

建立了考虑磨损的三自由度齿轮传动系统的多间隙非线性动力学方程，通过分岔图、李雅普诺夫指数谱、Poincar6映射和功率谱等定性和定量分析方法对系统的非线性特性进行了研究。通过研究发现了多间隙的齿轮动力学模型内部所蕴含的丰富非线性现象，发现以倍周期分岔为代表的常规进入混沌的途径，尧现了以交替出现的拟周期-锁相为途径的非常规途径，说明了含有磨损故障的多间隙齿轮传动系统非线性特性的复杂性和多样性。

发现了系统在拟周期运动时功率谱上的不可约频率分量之间满足频率相加法则，尧现了嵌套在拟周期运动窗口内部的三个相邻多周期运动之间满足两边的周期运动周期数之和与中间的周期运动周期数相等的规律。这都表明在齿轮传动系统非线性动力学模型也满足非线性系统的诸多普适性规律。

通过分岔图和对应的李雅普诺夫指数谱还可以发现具有[，-，-，-，-，-]的混沌振动区间和具有[，，-，-，-，-]的超级混沌振动区间总共大致存在7个。而-般齿轮传动系统的工作区间选取在亚临界和超I临界区间，本文研究系统非线性动态特性随系统激励频率的变化规律，能够为齿轮传动系统动态设计与动态性能优化提供依据。

本文在考虑齿面磨损对于动力学方程的影响只是考虑了对于时变啮合刚度的变化幅值的影响，而实际上齿面磨损可能还会引起齿侧间隙、齿面误差激励幅值等其它动态特性参数的变化，为此需要定量的研究齿侧间隙和齿面误差激励等随磨损的动态变化规律，今后将进行更加深入和系统的研究。