滚动轴承复合故障振动建模及Lempel-Ziv复杂度评价

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

Vibration modeling of rolling bearings considering compoundmulti-defect and appraisal with Lempel-Ziv complexityZHU -sheng，YUAN ，ZHANG You-yun，ZHU Chuan-feng，ZHOU Zhi(Theory of Lubrication＆Bearing Institute，Xian Jiaotong University，Xian 710049，China)Abstract： The vibration modeling is a key link to explore the cognitive rules in fault diagnosis.In order to foreseethe vibration characteristics accurately and eficiently，the transfer path and dynamic contact among joint interfaces werewell studied.A novel methodology for dynamic modeling of 6DOF systems consisting of beating inner race，outer race andhousing was presented based on the theories of Hertzian contact.The progressive processes of single and compound multi-defect were provided according to the practical kinematics of rolling bearings.The model describes the vibrations in timeand frequency domains.Computed results from the model were according to the validated with experimental results，whichwere generated on defective deep groove ball bearings.The complexities of vibration signal change with the emergence andgrowth of faults.Based on nonlinear dynamics theory，the Lempel-Ziv complexity measure was used to characterize thecomplexity of single and compound multi-defect signals. The results show that the complexity measure used as aquantitative criterion can effectively appraise the running condition of rolling bearings。

Key words：rolling bearings；compound multi-defect；vibration modeling；Lempel·Ziv complexity滚动轴承被广泛应用于机械工业，是基逮也是关键核心件。由于其通常工作在高速、高温、重载的恶劣环境下，滚动轴承疲劳导致的失效往往在整机故障中占大部分比率。因此故障诊断对于滚动轴承的安全、连续运行和预知更新至关重要1j。故障建模研究故障状态下动力学参数与响应征兆的内在联系，籍此获取服役正常与异常的证据信息 J，其 目的在于根据机械设备性能退化响应特点来识别潜在故障和构建故基金项 目：国家 自然科学基金资助项 目(51035007)；国家重点基础研究发展计划(973计划)资助项 目(2011CB706606)；国家科技重大专项资助项 目(2010ZX04001-021)；国家科技 重大专 项资助 项 目(2010ZX04014-015)收稿 日期：2012-06-26 修改稿收到日期：2012-09-19第-作者 朱永生 男，博士，副教授，1973年 1O月生障库。

滚动轴承动力学建模-直为学者所重视，Fukata3 J较早分析了滚动轴承的变柔性振动。研究结论显示，滚动轴承支承下的不平衡转子具有两个激励源，即来自不平衡的转频激励和来 自轴承总体刚度周期变化的内部参数激励。Tiwari等 采用两自由度模型研究轴承游隙造成的转子系统的非线性行为和稳定性。然后Lioulios等 探讨了速度波动对轴承非线性动力学的影响。Patil等 模拟内、外圈局部故障，研究单点缺陷所激励的时、频域特性。Liew等 进-步研究非线性Hertzian接触、游隙、转速等因素对动态行为的影响机制，并有效模 拟 了轴 承元 件 的振 动形态。Sopanen等 提出-种轴承模型用于分析局部及分布式故障，模型考虑了非线性 Hertzian接触，但是缺乏外圈 -24 振 动 与 冲 击 2013年第 32卷轴承座 的耦合效应。澳大利亚 UNSW 大学 Sawalhi等 采用附加自由度的方法建立 5DOF模型，从而得到 r相对精确的振动响应。-些新的建模手段也被引入，Mohsen等 用键合图功率流构建滚动轴承振动模型，分析了轴承元件的动态接触行为。机械设备出现故障时，故障的演变存在-个由轻微到严重、由单-到复合的过程 。对 F长期在变载工况下服役的机床轴承、机车轴承尤为显著，监测到的故障往往以多故障耦合的形态出现，研究滚动轴承复合故障建模有着重要工程意义 。然而，现有的动力学模型均以单点缺陷为研究对象，复合故障建模的研究鲜见报道。

为能够有效预示多耦合故障诱发的振动特征，基于轴承运动学与 Hertzian接触理论，建立-种计及内圈- 外圈 -轴承座的6DOF振动模型。通过渐进过程描述缺陷区域，以深沟球轴承 SKF6205-2RS为实例计算振动响应。与实验数据对比表明，利用该方法可深入揭示轴承局部缺陷动力特性与宏观动态性能之间的映射关系▲-步着重研究复合故障建模方法及其响应特点。在此基础上引入非线性动力系统中复杂性刻画疗法--I empel-Ziv复杂度指标用于评价不同故障状态动力系统的变化。研究表明，这种特征量是表征不同故障信息的有效参数。

1 滚动轴承接触建模1.1 非线性弹性恢复力如图 1，第 i个滚动体 -套圈弹性趋近最6 是套圈的径向位移( ，Y )，滚动体位置角 0 和游隙 的函数，由运动关系可知：6 cos0 y sin0 -i：1，2，。n (1)图1 滚动轴承运动示意图，- Fig.1 Schematic diagram滚动体位置角 0 。f mm g b。 ring fi。

60～t其中： 为公转速度，n为滚动体个数。设 为相配轴的转速，D 为滚动体直径，D 为轴承节圆直径， (1- )。设轴承座运动坐标为(Xp )，则内圈- 轴承座的综合弹性趋近量就为： ( - )cos0 (Y -Yp)sin0 -y (2)需要说明的是，由于载荷分布区具有角度范围， 只能取正值。

1.2 载荷 -位移关系H根据 Hertzian理论，点接触载荷-位移关系为：F 硒 (3)式巾：艿为弹性趋近量，K为Hertzian接触刚度。为了描述两个啮合表面之间的接触状况，定义 K、 ，、∑JD分别为滚动体与内、外国的接触刚度与曲率和，则：， ( )(∑ ㈤/xi1 (4)Kou，er ( )(∑pJ)m ( 。。 ㈩其中：E为弹性模量， 为泊松比， 为变形系数。这样，-旦轴承几何尺寸确定， ，、 就容易得到▲而求出轴承整体等效接触刚度为：K 1 (6)滚动轴承运动副的等效约束力可由每个滚动体所受弹性接触力求得，故总的非线性弹性恢复力可表示为： ∑K[( - )cos0。

(Y -Y )sin0 - ] (7)1.3 接触阻尼滚动轴承接触阻尼随平均负荷变化，其值很小，Kramer通过动态测试方法给出了经验计算公式 C (0.25-2.5)×10 K Ns/mm (8)2 振动系统模型滚动轴承振动性能是其内部运动关系在外部载荷作用下与传递路径的综合作用。内圈 -外圈 -轴承座为非线性强耦合关系，图 2表达了轴承接触界面的多自由度振动形态。其中，c。、 。、C。、 分别为外圈、轴承座的阻尼和刚度。基于多 自由度振动理论，构建多体系统环境下的轴承精细化模型如图2所示。由于滚动轴承接触问题的复杂性及各种非线性因素的影响，无疑会增加系统建模及求解的难度。鉴于此，本研究做以下基本假设。① 忽略滚动体滑移现象。② 忽略制造误差，如表面波纹度的影响。③ 假设滚动体与内、外圈之问只存在局部弹性接触变形。由第-节的分析结合图2～3，建立振动方程如下：图2 滚动轴承振动示意图Fig.2 Schematic diagram of roling beating vibration第 16期 朱永生等：滚动轴承复合故障振动建模及 Lempel-Ziv复杂度评价 25，vc/c∑ ( 1Ⅳc夕∑ (1- p)cos0 ( -Yp)sin0 - ] ～c。s0 mseto2eOStot- 3fp)cos0 (y -Yp)sin0 - ] ～sin0 m g m e sintot- ∑K[( -Xp)COs0 ( -Ye)sin - c。s0 0i1m。 。 (c。c )夕。( 。ke)y。-- -kpYe-∑K[( - )cos0 (y -Ye)sin0 - sin0。01p p P p( P- 。)0r凡P p Cpyp ke(Yp-Y。)：mPg式中：m 为内圈(相配轴)质量，m。、m 分别为外圈和轴承座质量，e为不平衡距离 为外载荷。

3 局部故障嵌入滚动轴承的失效大部分表现为局部缺陷，例如疲劳剥落、凹痕、点蚀等。缺陷的存在会形成-个有限的区域，具有长度 z，如图3所示。当滚动体运行至缺陷位置时，弹性趋近量瞬间增加，设 d为相应的实际增加值，即式(2)中的游隙 增加 d，这将导致冲击产生和麓F，。

图 3 表面缺陷区域Fig.3 The surface defects region能量突变。在外圈固定内圈旋转的情况下，外圈缺陷所对应的中心角 Ot2arcsin(1/2R)，R为外圈半径。

设缺陷区域相对初始滚动体的位置角为 ，则滚动体通过外圈缺陷区域的渐进过程可描述为：c。s( )≥c。s(c。 )≥c。s( - )1 。

sin(tO n)≤0 J内圈出现故障时，缺陷区域通过-个滚动体后每次经过 的角度与另-个滚动体相遇。内圈缺陷所对应的中心角 O/2arcsin(1/2r)，r为内圈半径。由于缺陷要 随 内 圈 同 时 旋 转，缺 陷 角 位 置 为 to t±arcsin(1/2r)，因此滚动体通过内圈缺陷区域的渐进过程为：to tarcsin(I/2r)≥to t2-7r .≥ ∞ t- aresin(I/2r) (11)滚动体出现故障时，缺陷区域以转速to 自转-周与(9)内、外圈分别接触-次，每-次接触弹性趋近量均增加d，此时判别函数由两部分组成。设滚动体缺陷区域与外圈初始相对角位置为 ，则缺陷区域通过外圈的渐进过程为：0.52aresin(1/D )≥sin(∞6t)≥0.5-、2arcsin(1/O ) (12)c。s( )≥0 J同理，滚动体缺陷区域通过内圈的渐进过程为：- 0·52arcsin(1/De)≥ 、sin(∞6t)≥-0.5-2aresin(1/D ) (13)。 (tot ≤0 J复合故障指同时发生两个或多个互相关联、交叉影响的故障。以内、外圈复合故障为例，动态建模表现为式(10)～(11)的级联，即同时运行内、外圈的渐进过程。

4 计算分析及验证式(9)所代表的状态方程组非线性强，难以解析，本文采用数值方法获取系统响应。计算方法为定步长Runge-Kuta法，步长设为 1/12 000 S，初始值为 10in。算例采用美国凯斯西储大学(CWRU)轴承研究中心的深沟球轴承(SKF6205-2RS)，转速 1 772 r/min，采样频率 12 000 Hz。故障轴承采用电火花加工制成，内、外圈及滚动体的故障频率依次为 5.415 2、3.584 8、4.713 5乘以工频，CWRU提供了详细设计参数可供参考 。外圈发生故障时频谱结构相对简单，主频为故障特征频率，如图4所示∩以看出，仿真信号包络谱体现了不平衡激励的转频fr29.66 Hz，内圈故障频率fo106 Hz及其谐波 、 、 成分，这与实验信号非常吻合。与外圈故障的区别在于内圈缺陷连同转轴-起运行，故障频率必然受到工频的调制，因此包络谱的主频以f/159.9 Hz 189.5 Hz - 130.4Hz及其组合谐波成分 十 、 - 、 、 ～ 等为主。图5的仿真信号包络谱成功突出了主频成分，并且频率预测达到了很高的精度。仿真和实验数据的对比表明，振动模型对显著频率成分实现了良好的预26 振 动 与 冲 击 2013年第 32卷测。为验证模型的鲁棒性，考察转速为 1 797 r/min时滚动体故障的振动响应。从图6中清晰看出滚动体缺陷激励引起的振动在包络谱中表现为保持架转频 11.9 Hz对故障频率. 141.2 Hz的调制。这是由于滚动体上缺陷位置随保持架旋转不断改变，缺陷引起的冲击响应幅值必然受到保持架旋转频率调制。以上分析充分体现了本文建模方法的可靠性。

g6- ×j磐I- x 29 66 L 1.上。3.1 11 6 24；， .色 4主l，，H(a)仿真信号包络谱l从 L 王 。

58 i 1flHz(b)实验信号包络谱图4 外圈故障仿真与实验信号包络谱Fig.4 Simulation and experimentalenvelope spectrum of outer race defect50 ，∞ 50 200 260 300 350 400 450 5∞ 550l，/Hz(a)仿真信号包络谱I x 1甜 9× 蚰 l Y 0054071 X：1 89 5 T 。

：l， ： fl ， 36。l ： 上工.L： iji ：ij .IJ L 上。 猢撩z(b)实验信号包络谱图5 内圈故障仿真与实验信号包络谱Fig.5 Simulation and experimentalenvelope spectrum of inner race defect设置相同的缺陷长度 f0.5 mm和深度 d0.1mm，转速 1 772 r/rain时内、外圈复合故障振动响应如图7所示。滚动轴承复合故障的信号特征比单-故障相对复杂。从刚度损失的角度可以解释为复合故障导致轴承接触刚度发生更大变化，在不同缺陷位置引入局部柔度，使得动力学性能改变。当内、外圈复合故障出现时，局部柔度同时增加，激励系统出现两种故障频率和. 。包络谱低频部分 及其谐波比较明显， 的幅值相对较猩以解释为复合故障由单点故障呈矢量叠加方式构成，以内、外圈复合故障为例，频域总幅值为：4 20 0×趔粤 .4f 曼b- ×骠I量6- X迥86查 4X 20//H(b)实验信号包络图6 滚动体故障仿真与实验信号包络谱Fig.6 Simulation and experimental envelopespectrum of roling element defect4邑 2旦 0Xj型 .2孽.4吕6- X：鲤山 . lI.L I JTr rr， r 。

t/24 i2 3 4 5/x 10 /HZ(d)包络谱图图7 内、外圈复合故障信号及包络谱Fig.7 Vibration response of innerand outer race compound defectsI 4( I：I ai. r( 2( I丁(14) 。 121A (1厂)1A ( )1cosqb 1式中：lA ( I、l ，(-厂)I、 分别为内、外圈故障频率幅值和相位角。总幅值是单点故障引起幅值的矢量和，多个频率分量之问呈现-定角度并具有方向性。

相位-致时，振幅增大；相位反向时，振幅减小，表现为强度不平衡。这是由故障接触区域的不同，通过缺陷区域的渐进过程不同造成的。复合故障同样具备单点故障的所有频率特征， 、 - 、 、 等频率成分在包络谱中均得以体现。图8～9分别为滚动体与内、外圈的复合故障响应。从中可以看出滚动体故障频率. 不突出，而其倍频成分比较显著，论证了复合故障具有强度的不平衡性及特征耦合性。

图10为内圈、外圈、滚动体三种故障复合时的振动响应，包络谱清晰显示了三种特征频率的2倍频，5 4 3 2 O 5 4 3 2 O - .Ⅲ-I0-× - .山- oI×第 16期 朱永生等：滚动轴承复合故障振动建模及 Lempel-Ziv复杂度评价15二 1052 0吕6- ×趔馨兰0-2罂 .d罾6- ×：理罂譬要篓-量篓t/s(b)细节信号l 2 3 4 5x 10VHz(d)包络谱图图8 外圈、滚动体复合故障信号及包络谱Fig.8 Vibration response of outer raceand roling element compound defects-4毒 2言 0 A L 厶 I f r r图9 内圈、滚动体复合故障信号及包络谱Fig.9 Vibration response of inner raceand roling element compound defects139.2 Hz相对微弱，并且各频率的幅值均不同。计算结果充分表明滚动轴承复合故障具有相互耦合、相互制约和强度不平衡性，也说明这些故障的产生具有内在确定性机制。同样显示，包络谱所提取的故障频率倍频分量可作为识别复合故障的有效依据。

5 复合故障Lempel-Ziv复杂度评价滚动轴承复合故障往往使其频率组成及各成分幅、频特性复杂化，从而使时域波形发生不同程度的畸变，这种畸变也是轴承故障在时域信号中的直接反映。

事实上，信号的波形、频谱结构发生变化时，信号的复杂程度尧生了变化。人们开始探索采用-些指标定量地描述轴承振动信号的复杂性 J-建东等 用Lempe1.Ziv复杂性测度定量评估大机组运行状态。

图10 内圈、外圈、滚动体复合故障信号及包络谱Fig.10 Vibration response of inner race，out race and roling element compound defectsHong等 研究了 Lempe1.Ziv复杂度指标用于评估故障轴承损伤状态。实践表明，这～指标是衡量有限时间序列复杂程度的高效工具。

序列S ，s2，，s 的Lempel-Ziv复杂度可按c (r)(r Ⅳ)经过 Ⅳ次循环得到，计算过程如下 J引：(1)初始化S删 ，Q。 ，C (0)0，r1。令 Q Q ，由于 Q 不属于 s - ，则 C (r)C (r-1)1，Q ，rr1。

(2)令 Q Q ，判断 Q，是否属于 S - S rs - ；若属于，则 C (r)C (r-1)，rr1，重复过程(2)。

(3)若不属于，则 c (r)C (r-1)1，Q， ，r，1，重复过程(2)。

例如序列0000的分段形式为0·000，C 2。符号序列 010101的分段形式为 0·1·0101 ，C 就为 3。出于对数据严格的可比性考虑，Lem。

pel和Ziv进-步提出了-种归-化基准公式，将 C 值界定在 [O 1]之间，详细算法为：0≤C ： ≤1 (15) Ⅳ u UL. NC，Ⅳ。l imCⅣ(Ⅳ)南 (16)其中：k为s 中元素的个数，对二进制序列k：2。根据已有的经验，-般 Ⅳ大于 3 600可保证上式成立 。

根据上述 Lempel-Ziv复杂度算法和文献[2，16]可知，- 个序列的复杂度越大，说明它的添加操作越多，序列的内涵分量也越多，描述给定符号序列所需最少的、互不相同的分段形式越多，给定的序列周期性越弱，出现新内涵模式的速度也就越快。Lempe1.Ziv复杂度的物28 振 动 与 冲 击 2013年第 32卷理意义在于它反映了-个时间序列随着长度的增加出现新模式的速率。复杂度越大，说明数据在窗口长度时期内随时间出现的新变化越多，发生新变化的速度越快，表明这-时期的数据变化是无序而复杂的。反之复杂度越小，则说明发生新变化的速度越慢，数据变化是规则的，周期性越强。因此，振动信号的Lempel-Ziv复杂度指标能够客观反映出系统状态发生变化的情况。

为进-步完善轴承故障评估方法，本文借助上述Lempel-Ziv复杂度指标的建立思想，并将之推广到复合故障诊断领域。同样设置缺陷长度 f0.5 mm和深度d0.1 mm的微弱故障进行数值计算。图 11从左往右依次对应正常、外圈、内圈、滚动体、内 -外圈复合、外圈 -滚动体复合、内圈 -滚动体复合、内外圈 -滚动体复合故障的 Lempel-Ziv复杂度变化趋势，表 1为其具体数值。

- X.三N ＆昌状态类别图 11 不同状态的 Lempel-Ziv数值变化趋势Fig.1 1 Lempel- Ziv complexity VS.diferent fault表 1 不同状态的 Lempel-Ziv复杂度Tab.1 Num erical value of Lem pel-Ziv complexity从结果中可以看出，正常状态的振动最为简单。

滚动体故障虽然激励冲击能量不大，但是滚动轴承中滚动体的运动比较复杂，既有公转又有 自转。在其运动过程中，故障区域既要与内圈接触又要与外圈接触，渐进过程由两部分组成(式(12)～(13))。不断变化的接触行为会引起强烈的调制作用，反映在频谱中会出现越来越多的谱峰及其谐波，波形中的成分越来越杂乱，振动的周期性减弱，非线性增强，相对于内、外圈更加倾向于较高的Lempd-z 值。复合故障破坏了单- 故障的周期冲击振动结构，故障冲击加剧，直观上可以看出，信号的复杂度增加了。不同形式的复合故障其响应的复杂度也不同，这是因为不同的复合故障形成机制不同，导致信号的冲击振动结构有所差异。内外圈-滚动体三种缺陷复合故障的振动响应与两种缺陷复合故障的 Lempel-Ziv指标差别不太大，体现了故障之间既相互干扰又相互制约。对复合故障的复杂度评价与我们的工程经验相吻合，再次证明了本文故障建模方法的可靠性。由上述分析可知，振动信号的复杂度较为敏感地反映了轴承状态演变过程中动力系统的变化历程，是-种有效的特征量。实际应用中整个信号的 Lempe1.Ziv指标单独作为诊断特征量稍微少了- 些。基于波形匹配思想，兼备优良性质的多小波拥有多个尺度函数和小波函数，可以很好地匹配复合故障的多种特征波形 。利用多小波变换消除噪声和分解信号，从而构造不同尺度下的 Lempe1.Ziv复杂度向量，可以同时完成单点及复合故障的监测诊断。此方法的实验研究加之工程应用，亦是下-步的研究课题，将另文给出。

6 结 论(1)建立-种计及内圈 -外圈 -轴承座的滚动轴承6DOF振动模型。从轴承几何学、运动学的角度推导了局部故障的渐进过程，通过与实验数据的比较表明，该方法能有效预测振动频谱的主要成分，为后续复合故障的研究奠定了基矗(2)进-步研究了复合故障的建模方法和振动形态。研究结果显示，复合故障具有相互耦合性、相互制约性和强度不平衡性。包络谱或共振解调谱可同时解耦复合故障特征频率及其倍频分量。-定条件下，包络谱的倍频分量可以作为诊断复合故障的依据。

(3)滚动轴承局部故障的出现相当于局部柔度的引入，进而在不同渐进过程瞬间改变系统动力学特性。

不同故障信号的Lempel-Ziv复杂度指标反映了不同故障的动力学形成机制。轴承振动信号的复杂度能有效地评价轴承运行状态，灵敏地反映其状态的变迁过程。