SERVICE LIFE EVALUATION OF TRACTION MOTOR BEARINGS IN COMPLICATED VIBRATION ENVIRONMENT OF A LOCOMOTIVE
-
摘要: 牵引电机是铁路机车的动力源, 其关键零部件(支承轴承等)的服役性能直接影响机车传动系统的稳定性和可靠性. 对于重载机车, 传统轴承服役寿命评估方法主要基于定载荷工况, 难以准确评估轨道不平顺等复杂外部激励作用下电机轴承的服役寿命. 因此, 本文根据车辆-轨道耦合动力学理论, 考虑轨道车辆运行过程中的轮轨相互作用和齿轮啮合作用, 建立了具有牵引动力传动系统的机车-轨道耦合动力学模型; 采用线性损伤累积准则和ISO 281标准计算方法, 评估了复杂机车振动环境下牵引电机轴承的服役寿命. 结果表明, 在轨道随机不平顺激励下, 机车轮轨垂向力、齿轮啮合力、牵引电机内部转子离心力、不平衡磁拉力等明显增大; 在复杂机车振动环境中, 电机轴承内部滚子-滚道相互作用加剧, 传动端与非传动端轴承的疲劳寿命缩短; 随着线路状态的不断恶化和机车运行速度的提高, 牵引电机轴承的预测寿命里程不断减小; 由于传动端轴承承受较大的外部动态载荷, 传动端轴承的服役寿命明显低于非传动端轴承. 本文提出的评估方法可为机车牵引电机轴承的设计、选型和寿命评估提供理论指导.Abstract: Traction motor is the power source of the railway locomotive, therefore, the stability and reliability of the locomotive transmission system are directly affected by the service lives of its key components (e.g., the support bearings, etc). For the heavy haul locomotive, the traditional evaluation methods are mainly based on the fixed load conditions, therefore, it is difficult to accurately evaluate the service lives of motor bearings under complex external excitations such as track random irregularities. Based on the vehicle-track coupling dynamics, a locomotive-track coupled dynamics model with traction power transmission system considering the specific structures of the traction motor and rolling bearing is established. In addition, the effect of the impact wheel-rail interaction and gear engagement is considered in this study. Besides, the service lives of traction motor bearings in the complex vibration environment of a locomotive are assessed based on the linear damage accumulation criterion and ISO 281 standard. Results indicate that the vertical wheel-rail interactive force, gear meshing force, the centrifugal force of rotor, and unbalanced magnetic pull increase obviously under the external excitation induced by the track random irregularity. In the complex vibration environment of the locomotive, the interactions between the roller and race in motor bearings are intensified, and the corresponding fatigue lives of the driving end and non-driving end motor bearings are shortened. Furthermore, with the continuous deterioration of the line state and the improvement of the locomotive running speed, the predicted life mileages of the motor bearings decrease in the increasingly harsh service environment. Because of the large external dynamic loads, the service life of the driving end bearing is shorter than that of the non-driving end one. This evaluation method proposed in this study can provide a theoretical guidance for the design, selection, and service life prediction of motor bearings used in the locomotive.
-
Keywords:
- railway vehicle /
- rolling bearing /
- life evaluation /
- traction motor /
- gear engagement
-
引言
随着我国轨道交通事业不断发展, 为满足“高速化”、“重载化”的发展需求, 铁路机车的牵引功率不断提高, 为机车牵引动力传动系统的稳定性和可靠性带来严峻的挑战. 作为铁路机车的动力源, 牵引电机的动态特性和服役寿命直接影响机车的运行品质和安全性. 据统计, 轴承故障约占牵引电机故障的44%[1]. 滚动轴承故障将影响电机传动效率, 产生异常振动[2], 甚至卡死导致切轴, 严重威胁机车的运行安全. 因此, 为保证牵引电机的服役性能和使用寿命, 必须探明电机轴承在复杂机车振动环境中的动力学特性, 准确评估其疲劳寿命.
作为极易失效的部件, 滚动轴承的动力学特性研究及其疲劳寿命预测获得了研究学者们的广泛关注, 其中, 滚动轴承内部动态载荷分布是轴承动力学研究的核心问题. 例如, 涂文兵等[3]分析了加速工况下滚动轴承的动态载荷特性, 研究了角加速度、径向游隙等外部激励及结构参数对轴承内部动态载荷的影响. Liu等[4]研究了保持架柔性对滚动轴承打滑特性的影响, 并根据滚子与滚道间的切向摩擦分析了径向载荷与内圈角加速度间的匹配关系. Liu等[5]建立了关于角接触球轴承的声振耦合动力学模型, 研究了存在滚道表面缺陷时轴承的声学特性. 马帅军等[6]基于ADAMS软件编写了Gfosub作用力子程序, 建立了角接触球轴承的动力学模型. 李志农等[7]提出了一种滚动体表面点蚀变形渐变释放模型, 分析了滚动体局部点蚀时的故障特征. 由于滚动轴承工作环境密闭, 现有技术手段无法直接观测其内部零部件的服役状态, 上述研究充分考虑了轴承各部件间的相互作用, 通过建立可靠的滚动轴承动力学模型, 获得了较为准确的轴承内部作用力. 但是, 恒定外部载荷的理论假设对于工业装备, 尤其是铁路机车、汽车等受到复杂、强烈外部激扰的移动装备, 不足以反映其支承轴承的真实动力学响应. 因此, 部分学者已经开始关注外部振动环境对滚动轴承服役状态的影响. Wang等[8]和Liu等[9]分别研究了轨道车辆轴箱轴承和电机轴承的动力响应, 分析了车轮多边形、轨道随机不平顺等外部激扰对滚动轴承服役性能的影响. 耿自林等[10]分析了轨道不平顺激扰下高速列车轴箱轴承的摩擦功耗, 研究了车轮多边形、轴承内、外圈故障对轴箱轴承温升的影响.
在工业领域, 滚动轴承的服役寿命是装备设计、寿命分析的重要依据. 1947年, Lundberg和Palmgren[11-12]根据Weibull金属疲劳概率分布提出了轴承寿命计算标准理论, 即L-P理论. 经过众多研究学者的不断完善和修正[13-14], 形成了ISO标准[15], 并得到广泛应用, 因此, 我国于2011年将其引入成为新的国家标准[16]. 然而, 该标准基于额定外部载荷, 通过修正系数对复杂外部激扰下的轴承服役寿命进行修正, 评估结果与工程实践仍存在较大差距[17]. 为克服该问题, 众多学者开展了大量的理论研究. 例如, 刘德昆等[18]基于实测轴箱动载荷, 参考ISO标准计算方法, 提出了一种高速列车轴箱轴承疲劳寿命的评估方法. 赵礼辉等[19]基于实测轮心载荷谱获得了汽车轮毂轴承动态外载荷, 分析了轮毂轴承的服役寿命. 然而, 基于实测外部载荷的寿命评估方法因试验周期长、成本高、影响因素复杂等难以广泛应用. Li等[20]采用有限元方法研究了时变外载荷对轴承疲劳寿命的影响. 但有限元方法计算效率低, 不适用于长时间、大里程的动力学仿真. 随着计算机技术的不断发展, 众多以数据驱动的新型预测方法[21-23]被不断提出, 为滚动轴承疲劳寿命预测和评估指出了新的路径.
综上所述, 为准确评估电机轴承的服役寿命, 应充分考虑其工作环境的影响. 对于架悬式牵引电机, 其一端悬吊于构架上, 一端通过抱轴承与轮对连接. 牵引力矩通过齿轮传动转置传递至轮对, 产生纵向蠕滑力驱动车辆克服运行阻力前进. 由轨道随机不平顺引起的轮轨作用力通过机车结构振动传递和齿轮啮合作用传递至牵引电机. 基于车辆与轨道耦合动力学理论[24], Chen等[25]建立了考虑齿轮动态啮合作用的机车与轨道耦合动力学模型, 齿轮传动转置与车辆-轨道系统通过齿轮啮合和振动传递耦合起来. Zhou等[26]进一步考虑机车牵引电机电系统与机械系统间的耦合效应, 通过电机电流信号分析车辆系统的振动特性. 在这些模型中, 牵引电机被视作一个质量块悬吊于构架上, 无法获得电机内部的动力学响应. 实际上, 牵引电机由转子、定子、支承轴承、小齿轮组成, 电机轴承被安装于牵引电机转子与定子之间. 因此, Liu等[9,27-28]根据牵引电机及其支承轴承的结构特点, 建立了更精细的牵引电机模型, 并将其与机车−轨道系统耦合为一个整体. 该模型可以准确获得在转子自身重力、由转子动态偏心导致的离心力和不平衡磁拉力、齿轮啮合力等复杂外部激励作用下的电机轴承内部载荷, 为电机轴承的疲劳寿命评估提供了理论基础.
为准确评估机车牵引电机轴承的服役寿命, 本文建立了具有牵引动力传动系统的机车与轨道耦合动力学模型, 考虑了牵引电机及其支承轴承的结构特点, 通过轴承、电机、车辆的耦合作用, 计算了轨道随机不平顺和齿轮时变啮合刚度激励下的电机轴承动态接触载荷. 基于ISO标准计算方法和线性损伤累积准则, 分析了机车振动环境中滚动轴承的瞬时疲劳损伤, 评估了其疲劳寿命, 揭示了线路状态、车辆运行速度对电机轴承寿命的影响.
1. 牵引电机轴承服役寿命评估方法
1.1 机车-轨道耦合动力学模型
本文采用Liu等[2,9]建立的具有牵引动力传动系统的机车与轨道耦合动力学模型. 如图1所示, 耦合动力学模型包括车辆系统, 轨下结构和牵引动力传动系统. 其中, 车辆系统由车体, 构架, 轮对和轴箱组成, 被视为一个多刚体系统. 一、二系悬挂系统、牵引拉杆等机车各部件间的连接单元等效为弹簧(K)-阻尼(C)单元.
牵引动力传动系统由牵引电机、齿轮箱和轮对组成, 该系统通过振动传递和齿轮啮合与车辆系统耦合为一体. 牵引力矩通过齿轮副传递至轮轨界面, 产生纵向蠕滑力驱动机车克服运行阻力前进. 牵引电机与齿轮箱通过螺栓刚性连接, 可视作一个整体. 齿轮箱小齿轮装配或加工在电机转轴上, 大齿轮与轮对固结. 对于架悬式牵引电机, 其一端通过电机吊杆悬吊于构架上, 另一端由轮对通过抱轴承支承. 电机轴承安装于转子与定子之间, 起到支承转子、减小转子与定子间的摩擦、保证转子回转精度的作用. 根据电机轴承与小齿轮横向距离的远近, 电机轴承被分别命名为非传动端轴承和传动端轴承. 考虑仿真计算效率, 轨下结构仅考虑钢轨, 由钢轨胶垫和扣件提供垂向和横向支承. 钢轨柔性采用Euler梁进行模拟. 该耦合动力学模型的动力学方程、数值积分方法等详细内容可参考文献[9].
1.2 牵引传动系统作用力
机车牵引传动系统包括: 牵引电机、齿轮箱和轮对, 承受多种复杂的内、外部激励. 轮轨相互作用是连接车辆系统与轨下结构的纽带. 机车车轮与钢轨接触区域的法向非线性接触力P(t)可表示为[29]
$$ P\left( t \right) = {\left[ {\frac{1}{G}\delta Z\left( t \right)} \right]^{3/2}} $$ (1) 式中, G为轮轨接触常数, δZ(t)为车轮与钢轨间的瞬时法向弹性压缩变形量.
齿轮时变啮合刚度是齿轮传动系统的固有内部激励[29-30]. 在机车的复杂振动环境中, 轮轨激扰将显著影响大、小齿轮的运动状态, 甚至发生反向啮合, 此时, 时变齿轮啮合力Fm可表示为[9]
$$ {F_{\text{m}}} = \left\{ {\begin{array}{*{20}{l}} {{K_{{\text{mf}}}}\delta + {C_{{\text{mf}}}}\dot \delta }, & {{\text{for forward - side contact}}} \\ 0,&{{\text{for contact loss}}} \\ {{K_{{\text{mb}}}}\delta + {C_{{\text{mb}}}}\dot \delta },& {{\text{for backward - side contact}}} \end{array}} \right. $$ (2) 式中, Km为齿轮副的时变啮合刚度, Cm为齿轮啮合阻尼. 本文采用Chen等[30-31]基于势能原理提出的一种改进的直齿轮啮合刚度计算方法, 计算机车齿轮传动系统中齿轮副时变啮合刚度. 相应的啮合刚度计算公式可参考文献[30-31].
齿轮副的动态传递误差δ可表示为[9]
$$\begin{split} & {\delta _i} = \left\{ {\begin{array}{*{20}{l}} \left.\begin{gathered} {\theta _{{\text{p}}i}}{R_{\text{p}}} + {\theta _{{\text{g}}i}}{R_{\text{g}}}+ \hfill \\ {\left( { - 1} \right)^i}\left( {{Z_{{\text{p}}i}} - {Z_{{\text{g}}i}}} \right)\cos {\alpha _{\text{m}}} -\hfill \\ {\left( { - 1} \right)^i}\left( {{X_{{\text{p}}i}} - {X_{{\text{g}}i}}} \right)\sin {\alpha _{\text{m}}}- \hfill \\ {b_0} - {e_i} \hfill \\ \end{gathered} \right\}&\begin{gathered} {\text{for forward - side}} \hfill \\ {\text{contact}} \hfill \\ \end{gathered} \\ \left.\begin{gathered} {\theta _{{\text{p}}i}}{R_{\text{p}}} + {\theta _{{\text{g}}i}}{R_{\text{g}}} + \hfill \\ {\left( { - 1} \right)^i}\left( {{Z_{{\text{p}}i}} - {Z_{{\text{g}}i}}} \right)\cos {\alpha _{\text{m}}}+ \hfill \\ {\left( { - 1} \right)^i}\left( {{X_{{\text{p}}i}} - {X_{{\text{g}}i}}} \right)\sin {\alpha _{\text{m}}} -\hfill \\ {b_0} - {e_{{\text{r}}i}} \hfill \\ \end{gathered}\right\} &\begin{gathered} {\text{for backward - side}} \hfill \\ {\text{contact}} \hfill \\ \end{gathered} \end{array}} \right.\\& \left( {i = 1,2,3,4} \right)\\[-15pt] \end{split}$$ (3) 式中, 脚标p和g分别表示小齿轮和大齿轮, θ为旋转角位移, R为节圆半径, X和Z分别表示纵向和横向, αm为啮合角, b0为齿间间隙, e和er为由加工和装配导致的误差.
牵引电机是铁路机车的动力源, 考虑加工精度、装配误差、转轴及支承轴承的弹性变形等因素的影响, 牵引电机因定转子间气隙分布不均和转子质量偏心而诱发产生的内部动态作用力包括: 转子离心力, 不平衡磁拉力和定子与转子碰摩力. 转子碰摩现象通常在电机发生严重故障时产生, 本文研究机车正常运行状态下的电机轴承的服役寿命, 因此定-转子偏磨力暂未考虑. 转子离心力Fω可表示为
$$ {F_{\text{ω }}} = {M_{{\text{rot}}}}\dot \theta _{{\text{rot}}}^2r(t) $$ (4) 式中, Mrot为转子质量, θrot为转子旋转角位移, r(t)为转子动态偏心距.
轨道机车通常采用4极鼠笼式三相异步电机, 此类电机负载运行时的不平衡磁拉力FUMP在纵向和垂向分力可表示为[32]
$$ \begin{split} {F_{{\text{UMPx}}}} =& {f_1}\cos \gamma + {f_{3{\text{c}}}}\cos (2{{\dot \theta }_{{\text{rot}}}}t - 3\gamma )+ \hfill \\ & {\text{ }} {f_{{\text{3s}}}}\sin (2{{\dot \theta }_{{\text{rot}}}}t - 3\gamma ) + {f_{{\text{4c}}}}\cos (2{{\dot \theta }_{{\text{rot}}}}t - 5\gamma ) +\hfill \\ & {\text{ }} {f_{{\text{4s}}}}\sin (2{{\dot \theta }_{{\text{rot}}}}t - 5\gamma ) \end{split} $$ (5) $$ \begin{split} {F_{{\text{UMPz}}}} = &{f_1}\sin \gamma + {f_{{\text{3c}}}}\sin (2{{\dot \theta }_{{\text{rot}}}}t - 3\gamma ) -\hfill \\ & {\text{ }} {f_{{\text{3s}}}}\cos (2{{\dot \theta }_{{\text{rot}}}}t - 3\gamma ) - {f_{{\text{4c}}}}\sin (2{{\dot \theta }_{{\text{rot}}}}t - 5\gamma )+ \hfill \\ & {\text{ }} {f_{{\text{4s}}}}\cos (2{{\dot \theta }_{{\text{rot}}}}t - 5\gamma ) \end{split} $$ (6) 式中, 脚标x和z分别表示纵向和垂向, γ为转子中心偏移的方位角, 各分量幅值f1, f3c, f3s, f4c, f4s的取值可参考文献[32].
滚动轴承被安装在牵引电机转子和定子之间, 定-转子通过电机轴承传递相互作用. 轴承内、外圈分别装配在转子、定子上, 因此, 如图2所示, 滚动轴承内、外圈在第j个滚子处的径向相对位移δ可表示为
$$ {\delta _j} = Z\sin {\sigma _j} + X\cos {\sigma _j} - \dfrac{e}{2} - {h_{{\rm{mi}}j}} - {h_{{\rm{mo}}j}}$$ (7) 式中, X和Z分别为轴承内、外圈的纵向和垂向相对位移, σ为滚子角位置, e为径向游隙, hm为滚子与滚道间的中心油膜厚度, 脚标i和o分别表示内圈和外圈.
根据Hertz接触理论, 考虑滚子绕轴承形心公转过程中的离心效应, 滚动轴承内滚子与内、外圈滚道间的径向非线性接触力可表示为[3]
$$ \left[ {\begin{array}{*{20}{c}} {{N_{\text{i}}}} \\ {{N_{\text{o}}}} \end{array}} \right]{\text{ = }}\left[ {\begin{array}{*{20}{c}} {{K_{\text{e}}}\chi \delta _{}^{10/9}} \\ {{K_{\text{e}}}\chi \delta _{}^{10/9} + {M_{\rm{r}}}\omega _{{\rm{cr}}}^{^2}{R_{\rm{m}}}} \end{array}} \right] $$ (8) 式中, Ke为内、外圈滚道在滚子处的等效接触刚度, 可由滚子与内、外圈滚道间的接触刚度Ki, Ko获得; χ为接触参数, 当径向位移δ≥0时, 取值为1, 当δ<0时, 取值为0; ωcr为滚子的公转转速, 可根据滚动轴承结构参数通过内圈转速获得; Mr为滚子质量.
因此, 滚动轴承的纵向、垂向的支承合力可表示为
$$ \left. {\begin{array}{*{20}{l}} {{F_{{\text{ix}}}} = \displaystyle\sum\limits_{j = 1}^{{N_{b} }} {\left( {{N_{{\rm{i}} j}}\sin {\sigma _j} + {F_{{i} j}}\cos {\sigma _j}} \right)} + {C_{\text{x}}}\dot X} \\ {{F_{{\text{iz}}}} = \displaystyle\sum\limits_{j = 1}^{{N_{b} }} {\left( {{N_{{\rm{i}} j}}\cos {\sigma _j} - {F_{{i} j}}\sin {\sigma _j}} \right)} + {C_{\text{z}}}\dot Z} \\ {{F_{{\text{ox}}}} = \displaystyle\sum\limits_{j = 1}^{{N_{b} }} {\left( {{N_{{\rm{o}} j}}\sin {\sigma _j} + {F_{{o} j}}\cos {\sigma _j}} \right)} + {C_{\text{x}}}\dot X} \\ {{F_{{\text{oz}}}} = \displaystyle\sum\limits_{j = 1}^{{N_{b} }} {\left( {{N_{{\rm{o}} j}}\cos {\sigma _j} - {F_{{o} j}}\sin {\sigma _j}} \right)} + {C_{\text{z}}}\dot Z} \end{array}} \right\} $$ (9) 式中, j为滚子编号, Nb为滚子个数, C为滚动轴承等效阻尼.
1.3 滚动轴承寿命评估方法
目前, 为评估滚动轴承的疲劳寿命, 工程领域普遍采用基于L-P理论改进的ISO标准计算方法. 该标准引入了寿命修正系数a1和aISO, 滚动轴承的修正额定寿命LISO可表示为[15-16]
$$ {L_{{\rm{ISO}}}} = {a_1}{a_{{\rm{ISO}}}}{\left( {\dfrac{{{C_{{\rm{ISO}}}}}}{{{P_{{\rm{ISO}}}}}}} \right)^{10/3}} $$ (10) 式中, a1为可靠性寿命修正系数; aISO为基于寿命计算系统方法的寿命修正系数, 该系数考虑了润滑、污染等轴承服役条件的影响, 其取值可参考文献[15-16]; PISO为当量动载荷, 可根据式(7)和式(8)获得; 基本额定动载荷CISO可表示为[15-16]
$$ {C_{{\text{ISO}}}} = {b_{\text{m}}}{f_{\text{c}}}{\left( {k{L_{{\text{e}}}}\cos \alpha } \right)^{7/9}}{N_{\text{b}}}^{3/4}{\left( {2{R_{\text{r}}}} \right)^{29/27}} $$ (11) 式中, bm, fc为滚动轴承材料、加工、结构相关参数, k为滚子列数, α为公称接触角, Rr为滚子半径.
根据Palmgren-Miner线性损伤累积准则, 滚动轴承的疲劳失效可视为服役过程中动载荷作用下轴承疲劳损伤的累积过程. 假设滚动轴承服役过程中, 每一微小时刻内承受恒定的外部载荷, 滚动轴承的累积疲劳损伤可表示为
$$ D = \sum\limits_{i = 1}^m {{\varepsilon _i}} = \sum\limits_{i = 1}^m {\frac{{{l_i}}}{{{L_i}}}} $$ (12) 式中, εi为第i个积分步内滚动轴承的瞬时疲劳损伤值, m为积分步个数, li为第i个积分步内轴承内圈的回转圈数, Li为第i个积分步内轴承的额定内圈回转圈数. 当累积疲劳损伤值D大于1时, 滚动轴承发生疲劳失效.
相应的, 滚动轴承的预测寿命里程Xpre可表示为
$$ {X_{{\text{pre}}}} = \frac{{{X_{\text{c}}}}}{D} = \frac{{{X_{\text{c}}}}}{{\displaystyle\sum\limits_{i = 1}^m {{\varepsilon _i}} }} $$ (13) 式中, Xc为车体纵向位移, 即车辆行驶距离.
综上所述, 考虑轨道随机不平顺、齿轮时变啮合刚度、电机内部动态激励、轴承内部相互作用等复杂内、外激扰的影响, 本文提出了机车牵引电机轴承服役寿命的评估方法, 其流程如图3所示.
2. 机车牵引电机轴承服役寿命评估
本文以某HX型电力机车为例, 采用美国5级谱作为车辆系统外部激励, 运行速度为80 km/h, 仿真时长为5 s, 详细的机车主要动力学参数、牵引特性曲线、齿轮时变啮合刚度等详见文献[9]. 牵引电机传动端与非传动端轴承结构参数如表1所示. 相应的轨道不平顺和齿轮时变啮合刚度曲线如图4所示. 为突出提出的寿命评估方法的准确性, 提取了轮轨垂向力、齿轮啮合力、转子离心力和不平衡磁拉力等系统内、外部激励, 同时, 为不失一般性, 提取了传动端和非传动端电机轴承的第一个滚动体与内圈滚道的接触力, 以均方根值和最大值为统计指标, 对比分析了轨道随机不平顺对牵引电机服役寿命的影响.
表 1 牵引电机轴承参数Table 1. Main parameters of the traction motor bearingsValue Driving end bearing Non-driving end bearing radius of inner race/mm 53.5 45 radius of outer race/mm 72.5 62 radius of roller/mm 9.5 8.5 length of roller/mm 40 24 number of rollers 17 14 number of columns 1 1 bm 1.1 1.1 fc 84.3 84.3 在机车运行过程中, 轮轨接触、齿轮啮合、电机定与转子相互作用是机车及其牵引动力传动系统的主要激励, 其时程曲线如图5所示. 当不考虑轨道随机不平顺激扰时, 轮轨垂向力、齿轮啮合力、转子离心力和不平衡磁拉力由于齿轮时变啮合刚度和滚动轴承时变支承刚度等内部激励存在微小的波动, 但基本保持恒定. 此时, 电机轴承外部载荷与理论推导出的静载荷一致, 可视作采用ISO标准评估额定静载荷下的电机轴承服役寿命. 在轨道随机不平顺的影响下, 轮轨垂向力、齿轮啮合力、转子偏心力和不平衡磁拉力等内、外部激励均产生剧烈波动, 其均方根值各增大0.94, 0.35, 0.01和0.01 kN, 最大值各增大2.13, 6.48, 0.08和0.08 kN, 牵引电机轴承服役环境明显恶化.
如图6所示, 轮轨激扰将加剧电机轴承滚子和滚道间的接触作用, 传动端电机轴承的滚子与滚道接触力的均方根值和最大值分别增大0.09和2.19 kN. 预测寿命里程从192.34万公里降低至142.71万公里. 在轨道随机不平顺的影响下, 非传动端轴承滚子与滚道间接触力仅为传动轴承的1/10, 其预测寿命里程远高于传动端轴承. 因此, 对于机车牵引电机检修和维护应重点关注传动端轴承.
3. 运行条件对电机轴承服役寿命的影响
3.1 线路状态对电机轴承服役寿命的影响
轨道随机不平顺是车辆−轨道系统的重要激励. 随着服役时间的增加, 线路状态的不断恶化, 轮轨相互作用将逐渐加剧, 引起机车系统的剧烈振动, 加剧齿轮传动装置的啮合冲击, 恶化牵引电机及其支承轴承的服役环境. 本节选用4种轨道随机不平顺载荷谱(工况1: 无轨道不平顺, 工况2: 美国6级谱, 工况3: 美国5级谱, 工况4: 美国4级谱)模拟不同恶化程度的线路状态, 机车运行速度为80 km/h.
机车系统内、外部动态激励及牵引电机轴承的内部接触载荷、预测服役寿命的统计指标(均方根值和最大值)随线路状态的变化趋势如图7所示. 随着线路状态的不断恶化, 不断增大的轮轨垂向作用力导致齿轮传动装置振动冲击加剧, 在转子定子机电动态相互作用和齿轮时变啮合力影响下, 电机传动端和非传动端支承轴承滚子和滚道接触力明显增大. 相较于工况1, 仿真工况2、工况3、工况4的传动端电机轴承滚子与滚道间的接触力最大值依次增大2.51 kN, 4.19 kN和8.13 kN; 非传动端轴承最大值依次增大0.97 kN, 1.14 kN和1.9 kN. 工况4传动端电机轴承的预测寿命里程仅为理论寿命里程的62.03%, 服役寿命严重缩短(请注意图7(d)中传动端与非传动端对应的左右纵坐标范围). 线路状态严重威胁机车牵引电机轴承的服役寿命.
3.2 相对恒定车速对电机轴承服役寿命的影响
行车速度的提升可以有效提高运营效率, 但也将提高轮轨间的作用频率. 此外, 更高的牵引功率将直接加剧齿轮间的啮合冲击, 对于牵引电机及其支承轴承的可靠运行十分不利. 本节采用美国5级谱, 在不同牵引力矩下使机车保持50, 60, 70, 80和90 km/h的恒定车速, 分析不同车速下机车的内、外部激励和牵引电机轴承服役寿命的变化规律.
机车系统动态激励、牵引电机轴承内部载荷的统计指标及电机轴承的预测寿命里程如图8所示. 随着机车车速的不断提升, 轮轨相互作用增大, 齿轮啮合、牵引电机定子与转子间机、电相互作用加剧, 电机轴承承受的动载荷不断增大. 与车速为50 km/h工况相比, 当机车运行速度达到90 km/h时, 传动端和非传动端电机轴承滚子与滚道间的接触力最大值分别增大49.93%和145.45%, 预测寿命里程分别减少92.48%和72.08% (请注意图8(d)中传动端与非传动端对应的左右纵坐标范围).
4. 结论
为准确评估复杂机车振动环境下牵引电机轴承的服役寿命, 本文建立了基于牵引动力传动系统的机车−轨道耦合动力学模型, 考虑了轨道随机不平顺、时变齿轮啮合刚度、牵引电机内部动态激励、滚动轴承时变支承刚度等复杂内、外激扰. 基于本文分析条件的研究结果表明:
(1)在轨道随机不平顺的影响下, 轮轨力、齿轮啮合力、转子偏心力和不平衡磁拉力等内、外部激励产生剧烈波动, 瞬时冲击载荷导致电机轴承的服役环境显著恶化, 疲劳寿命明显降低, 此外, 牵引电机传动端轴承疲劳寿命远低于非传动端轴承;
(2)随着线路状态的不断恶化, 轮轨相互作用不断加剧, 加剧了齿轮啮合装置的啮合冲击, 缩短了牵引电机轴承的服役寿命;
(3)相同线路条件下, 机车运行速度的提升将加剧轮轨相互作用和齿轮啮合冲击, 影响电机轴承的服役状态, 缩短服役寿命.
此外, 限于篇幅, 本文仅评估了轨道随机不平顺影响下的电机轴承疲劳寿命, 可采用本文提出的方法进一步开展车轮踏面损伤(车轮多边形、扁疤和剥离等)及轨道结构激励(波磨、钢轨接头和道岔等)对牵引电机轴承的服役寿命的影响研究.
-
表 1 牵引电机轴承参数
Table 1 Main parameters of the traction motor bearings
Value Driving end bearing Non-driving end bearing radius of inner race/mm 53.5 45 radius of outer race/mm 72.5 62 radius of roller/mm 9.5 8.5 length of roller/mm 40 24 number of rollers 17 14 number of columns 1 1 bm 1.1 1.1 fc 84.3 84.3 -
[1] Nandi S, Toliyat HA, Li X. Condition monitoring and fault diagnosis of electrical motors—A review. IEEE Trans. Energy Convers., 2005, 20(4): 719-729 doi: 10.1109/TEC.2005.847955
[2] Liu YQ, Chen ZG, Li W, et al. Dynamic analysis of traction motor in a locomotive considering surface waviness on races of a motor bearing. Railway Engineering Science, 2021. doi: 10.1007 s40534-021-00246-x
[3] 涂文兵, 梁杰, 罗丫等. 加速工况下滚动轴承动态载荷特性研究. 振动与冲击, 2021, 40(9): 152-159 (Tu Wenbing, Liang Jie, Luo Ya, et al. Dynamic load characteristics of rolling bearing under acceleration condition. Journal of Vibration and Shock, 2021, 40(9): 152-159 (in Chinese) Tu Wenbing, Liang Jie, Luo Ya, et al. Dynamic load characteristics of rolling bearing under acceleration condition. Journal of vibration and shock, 2021, 40(09): 152-159(in Chinese)
[4] Liu YQ, Chen ZG, Tang L, et al. Skidding dynamic performance of rolling bearing with cage flexibility under accelerating conditions. Mechanical Systems and Signal Processing, 2021, 150: 107257 doi: 10.1016/j.ymssp.2020.107257
[5] Liu J, Xu YJ, Pan G. A combined acoustic and dynamic model of a defective ball bearing. Journal of Sound and Vibration, 2021, 501: 116029 doi: 10.1016/j.jsv.2021.116029
[6] 马帅军, 闫柯, 张晓红等. 滚动轴承运动-接触精确建模与二次开发. 西安交通大学学报, 2021, 55(8): 33-41 (Ma Shuaijun, Yan Ke, Zhang Xiaohong, et al. An accurate modeling method and secondary development of Bearing contact-motion. Journal of Xi’an Jiaotong University, 2021, 55(8): 33-41 (in Chinese) Ma Shuaijun, Yan Ke, Zhang Xiaohong, et al. An accurate modeling method and secondary development of Bearing contact-motion. Journal of Xi’an Jiaotong University, 2021, 55(08): 33-41(in Chinese)
[7] 李志农, 李云龙, 任帅等. 局部点蚀故障滚动体的滚动轴承动力学模型研究. 振动工程学报, 2020, 33(3): 597-603 (Li Zhinong, Li Yunlong, Ren Shuai, et al. Research on dynamic model of rolling bearing with local pitting fault in rolling bearing element. Journal of Vibration and Shock, 2020, 33(3): 597-603 (in Chinese) Li Zhinong, Li Yunlong, Ren Shuai, et al. Research on dynamic model of rolling bearing with local pitting fault in rolling bearing element. Journal of vibration and shock, 2020, 33(03): 597-603. (in Chinese)
[8] Wang ZW, Mei GM, Xiong Q, et al. Motor car-track spatial coupled dynamics model of a high-speed train with traction transmission systems. Mechanism and Machine Theory, 2019, 137: 386-403 doi: 10.1016/j.mechmachtheory.2019.03.032
[9] Liu YQ, Chen ZG, Wang KY, et al. Dynamic modelling of traction motor bearings in locomotive-track spatially coupled dynamics system. Vehicle System Dynamics, 2021. doi: 10.1080/00423114.2021.1918728
[10] 耿自林, 宋冬利, 张卫华等. 高速列车轴箱轴承动力学行为及温度分析. 机械, 2020, 47(11): 54-62 (Geng Zilin, Song Dongli, Zhang Weihua, et al. Dynamic and temperature analysis of axle box bearing of high-speed train. Machinery, 2020, 47(11): 54-62 (in Chinese) doi: 10.3969/j.issn.1006-0316.2020.11.009 Geng Zilin, Song Dongli, Zhang Weihua, et al. Dynamic and temperature analysis of axle box bearing of high-speed train. Machinery, 2020, 47(11): 54-62(in Chinese) doi: 10.3969/j.issn.1006-0316.2020.11.009
[11] Lundberg G, Palmgren A. Dynamic capacity of rolling bearings. Acta Polytech. Scand., Mech. Eng. Ser 1., 1947, 1(3): 7
[12] Lundberg G, Palmgren A. Dynamic capacity of rolling bearings. J. Appl. Mech., 1949, 16(2): 165-172
[13] Ioannides E, Harris T. A new fatigue life model for rolling bearings. Journal of Tribology, 1985, 107: 367-378 doi: 10.1115/1.3261081
[14] Tallian T. Simplified contact fatigue life prediction model-Part I: Review of published models. Journal of Tribology, 1992, 114: 207-213 doi: 10.1115/1.2920875
[15] ISO. ISO 281 method of calculating dynamic load rating and rating life of rolling bearings. Geneva: International Organization for Standardization, 2007
[16] GB/T 6391—2010, 滚动轴承额定动载荷和额定寿命 [17] 赵美茹. HX_D3C型和HX_N3型机车牵引电机轴承故障原因分析与解决措施. 铁道机车与动车, 2018(4): 35-38, 6 (Zhao Meiru. Cause analysis and solution of typical bearing failure of the traction motors for HXD3C and HXN3 locomotives. Railway Locomotive and Motor Car, 2018(4): 35-38, 6 (in Chinese) Zhao Meiru. Cause analysis and solution of typical bearing failure of the traction motors for HXD3 C and HXN3 locomotives. Railway Locomotive and Motor Car, 2018(04): 35-38 + 6(in Chinese)
[18] 刘德昆, 李强, 王曦等. 动车组轴箱轴承基于实测载荷的寿命预测方法. 机械工程学报, 2016, 52(22): 45-54 (Liu Dekun, Li Qiang, Wang Xi, et al. Life prediction method for EMU axle box bearings based on actual measured loadings. Journal of Mechanical Engineering, 2016, 52(22): 45-54 (in Chinese) doi: 10.3901/JME.2016.22.045 Liu Dekun, Li Qiang, Wang Xi, et al. Life Prediction method for EMU axle box bearings based on actual measured loadings. Journal of Mechanical Engineering, 2016, 52(22): 45-54(in Chinese) doi: 10.3901/JME.2016.22.045
[19] 赵礼辉, 李其宸, 冯金芝等. 随机道路载荷下轮毂轴承服役寿命预测方法研究. 机械工程学报, 2021, 57(10): 77-86 (Zhao Lihui, Li Qichen, Feng Jinzhi, et al. Research on prediction method of hub-bearing service life under random road load. Journal of Mechanical Engineering, 2021, 57(10): 77-86 (in Chinese) Zhao Lihui, Li Qichen, Feng Jinzhi, et al. Research on prediction method of hub-bearing service life under random road load. Journal of Mechanical Engineering, 2021, 57(10): 77-86(in Chinese)
[20] Li FK, Hu WP, Meng QC, et al. A new damage-mechanics-based model for rolling contact fatigue analysis of cylindrical roller bearing. Tribology International, 2018, 120: 105-114 doi: 10.1016/j.triboint.2017.12.001
[21] 周哲韬, 刘路, 宋晓等. 基于Transformer模型的滚动轴承剩余使用寿命预测方法. 北京航空航天大学学报, 2021. doi: 10.13700/j.bh.1001-5965.2021.0247 Zhou Zhetao, Liu Lu, Song Xiao, et al. Remaining useful life prediction method of rolling bearing based on Transformer model. Journal of Beijing University of Aeronautics and Astronautics, 2021. doi: 10.13700/j.bh.1001-5965.2021.0247 (in Chinese)
[22] 王新刚, 韩凯忠, 王超等. 基于迁移学习的轴承剩余使用寿命预测方法. 东北大学学报(自然科学版), 2021, 42(5): 665-672 (Wang Xingang, Han Kaizhong, Wang Chao, et al. Bearing remaining useful life prediction method based on transfer learning. Journal of Northeastern University (Natural Science) , 2021, 42(5): 665-672 (in Chinese) Wang Xingang, Han Kaizhong, Wang Chao, et al. Bearing remaining useful life prediction method based on transfer learning. Journal of Northeastern University (Natural Science), 2021, 42(05): 665-672(in Chinese)
[23] 董绍江, 吴文亮, 贺坤等. 基于性能衰退评估的轴承寿命状态识别方法研究. 振动与冲击, 2021, 40(5): 186-192, 210 (Dong Shaojiang, Wu Wenliang, He Kun, et al. Bearing life state recognition method based on performance degradation evaluation. Journal of Vibration and Shock, 2021, 40(5): 186-192, 210 (in Chinese) Dong Shaojiang, Wu Wenliang, He Kun, et al. Bearing life state recognition method based on performance degradation evaluation. Journal of vibration and shock, 2021, 40(05): 186-192 + 210(in Chinese)
[24] Zhai WM, Wang KY, Cai CB. Fundamentals of vehicle-track coupled dynamics. Vehicle System Dynamics, 2009, 47(11): 1349-1376 doi: 10.1080/00423110802621561
[25] Chen ZG, Zhai WM, Wang KY. Vibration feature evolution of locomotive with tooth root crack propagation of gear transmission system. Mechanical Systems and Signal Processing, 2019, 115: 29-44 doi: 10.1016/j.ymssp.2018.05.038
[26] Zhou ZW, Chen ZG, Spiryagin M, et al. Dynamic response feature of electromechanical coupled drive subsystem in a locomotive excited by wheel flat. Engineering Failure Analysis, 2021, 122: 105248 doi: 10.1016/j.engfailanal.2021.105248
[27] Liu YQ, Chen ZG, Zhai WM, et al. Dynamic investigation of traction motor bearing in a locomotive under excitation from track random geometry irregularity. International Journal of Rail Transportation, 2021. doi: 10.1080/23248378.2020.1867658
[28] Liu YQ, Chen ZG, Wang KY, et al. Surface wear evolution of traction motor bearings in vibration environment of a locomotive during operation. Science China Technological Sciences, 2021. doi: 10.1007/s11431-021-1939-3
[29] Sui SQ, Wang KY, Ling L, et al. Effect of wheel diameter difference on tread wear of freight wagons. Engineering Failure Analysis, 2021. doi: 10.1016/j.engfailanal.2021.105501
[30] Chen ZG, Zhou ZW, Zhai WM, et al. Improved analytical calculation model of spur gear mesh excitations with tooth profile deviations. Mechanism and Machine Theory, 2020, 149: 103838 doi: 10.1016/j.mechmachtheory.2020.103838
[31] Chen ZG, Ning JY, Zhou ZW, et al. An improved dynamic model of spur gear transmission considering coupling effect between gear neighboring teeth. Nonlinear Dynamics, 2021. doi: 10.1007/s11071-021-06852-y
[32] 周生通, 朱经纬, 肖乾等. 初始静偏心与重力载荷对动车牵引电机转子轴心轨迹的影响. 机械工程学报, 2020, 56(17): 145-154 (Zhou Shengtong, Zhu Jingwei, Xiao Qian, et al. Initial static eccentricity and gravity load on rotor orbit of EMU traction motor. Journal of Mechanical Engineering, 2020, 56(17): 145-154 (in Chinese) doi: 10.3901/JME.2020.17.145 Zhou Shengtong, Zhu Jingwei, Xiao Qian, et al. Initial static eccentricity and gravity load on rotor orbit of EMU traction motor. Journal of Mechanical Engineering, 2020, 56(17): 145-154 (in Chinese) doi: 10.3901/JME.2020.17.145
-
期刊类型引用(4)
1. 史秋磊. 地铁牵引电动机轴承复合缺陷振动检测技术研究. 中国机械. 2024(30): 88-91 . 百度学术
2. 丁杰. HX_D1C型电力机车牵引变流器功率模块的振动测试与分析. 湖南文理学院学报(自然科学版). 2023(02): 42-47 . 百度学术
3. 丁杰,蒋莉,徐庆怡. HX_D1B型电力机车设备的线路振动测试及分析. 机电工程技术. 2023(06): 63-66+165 . 百度学术
4. 王家鑫,张汉辰,吴志强,刘禹清,陈再刚. 基于正交试验的轴箱内置式高速动车驱动系统悬挂刚度优化方法. 交通信息与安全. 2023(06): 12-19 . 百度学术
其他类型引用(7)