一种双重更新的四因素随机退化模型齿轮寿命预测方法
【专利摘要】一种双重更新的四因素随机退化模型齿轮寿命预测方法，充分考虑实际齿轮退化过程中广泛存在的时变性、个体差异性、非线性和测量不确定性四方面随机因素，建立贴合其实际退化过程的四因素随机退化模型，并推导出该模型下齿轮剩余寿命基于观测值的概率密度函数；采用MLE方法对该模型参数进行初始化，避免参数初始化的盲目性；引入双重更新方法对模型状态及参数进行更新，获得更准确的更新结果，进而很大程度上提高了齿轮剩余寿命的预测精度。
【专利说明】
-种双重更新的四因素随机退化模型齿轮寿命预测方法
技术领域
[0001] 本发明属于机械设备剩余寿命预测技术领域，具体设及一种双重更新的四因素随 机退化模型齿轮寿命预测方法。
【背景技术】
[0002] 剩余寿命预测在齿轮状态监测数据的基础上为齿轮维护决策提供直接的指导和 依据，成为齿轮预知性维护策略不可或缺的重要环节。基于随机退化模型的齿轮寿命预测 方法将齿轮的剩余寿命视作随机变量，从监测信号中提取齿轮的健康状态信息，继而从概 率密度角度建立其剩余寿命的预测分布。该方法能够充分挖掘齿轮状态监测数据的动态特 性和随机特性，契合齿轮的失效机理，对齿轮的退化过程进行准确描述，因而受到众多学者 的广泛认同。在实际应用中，齿轮退化的随机性往往表现在四个方面:时变性、个体差异性、 非线性和测量不确定性，应用随机退化模型对齿轮的退化过程进行建模描述时需要同时考 虑运四种因素，然而现有模型大都仅考虑部分随机性因素，加之模型参数未进行初始化，模 型参数和系统状态更新方法表现欠佳，致使获得的齿轮剩余寿命概率分布偏离实际，剩余 寿命预测精度不尽人意。
【发明内容】
[0003] 为了克服上述现有技术的上述缺点，本发明提供一种双重更新的四因素随机退化 模型齿轮寿命预测方法，使得齿轮剩余寿命预测结果精度更高。
[0004] 为了达到上述目的，本发明采取的技术方案为：
[0005] -种双重更新的四因素随机退化模型齿轮寿命预测方法，包括W下步骤：
[0006] 1)使用监测系统采集齿轮振动信号，提取健康指标形成观测值向量y=[yi，...， ys ]T，S为独立观测值个数；
[0007] 2)构建四因素随机退化模型：
[000引
（1)
[0009] 其中，
[0010] Xk-tk时刻齿轮的真实退化状态；
[001U 0k-tk时刻衰退系数，表征个体差异性，为随机变量，服从分布WC%，请)；
[0012] b-模型阶数，表征非线性；
[00U] Ok-I-过程噪声，服从分布N(0,〇2(tk-tk-i))，《k-i = 〇(B(tk)-B(tk-i))，B(tk)为标 准布朗运动，地(tk)表征时变性；
[0014] yk-tk时刻齿轮健康指标观测值；
[001引 v_观测噪声，服从分布N(0, 丫 2);
[0016] 3)采用极大似然估计(MLE)方法对式（1)四因素随机退化模型的参数进行初始化， 具体步骤为：
[0017] 3.1)tk时刻齿轮健康指标观测值yk~N(xk，丫 2)，根据步骤1)所得观测值向量y计 算获得对数似然函数：
[001引
尚
[0019]其中，贷=[任2,你，请点rY，为四因素随机退化模型的未知参数向量，
Q = 02m+ 丫 2ls，M=[min{ti,tj}]i《i,j《s，Is为S阶单位矩阵；
[0020] 3.2)计算式(2)关于叫和成的偏导，极大化似然函数：
[0021]
货
[0022] 其中
I求解方程组(3)得到：
[0023]
W
[0024] 3.3)将式(4)计算结果回代至式(2)，简化对数似然函数为：
[00巧]
稱
[00%] 3.4)通过Matlab提供的寻优命令对式(5)进行寻优获得t和a，回带至式(4)求解 出山和，进而根据关系游=打2Ari，r =fV诚，-玲=鸣咕，Q =O2M+ 丫 2Is,求解出未知参 数向量@ = 与，的极大似然估计，作为式（1)四因素随机退化模型的初始参数 0;>=片2，,"。,(,，(请:,，6，^勺'；
[0027] 4)将式(1)四因素随机退化模型表示为式(6)形式，并带入初始参数0 0:
[002引 佩
[0029] 5：
[0030] 5)引入双重更新方法，对四因素随机退化模型的状态和参数进行实时更新，该方 法由两次相互补充、嵌套进行的更新过程组成：
[003。 5.1)外层更新开始，令k=0,对正态分布成沾A:)，即W(4，P。)进行随机采样，获得由N个初始 粒子构成粒子集柄},_1 V，对应的粒子权值集合为扣a=iv，其牛
[0032] 5.2)里层更新开始，令k = k+l，将Pk-IIk-I和N个未更新粒子扣1，/ = 1,2,…,妒，带入式 (7)迭代更新，计算获得Pk I k和4，Z' = I，2,...，W ;
[0033]
[0034]
[0039]
斷
[003引5.3)里层更新结束，由步骤5.2)所得？^1<和苗|*，<' = 1，2，...，^^形成正态分布列 W(4,巧如=1，2，...办，对每个分布进行一次采样获得新的粒子集挺^1、，各粒子权值更新 如下：
[0036] ㈱
[0037]
[0038] 5 - 4)由击驢5 - 3)所得粒子及权值计算tk时刻齿轮更新状态和参数：
[0040] 外层更新结束；
[0041] 6)将步骤5.4)的更新结果代入式（10)，获得tk时刻基于观测值的齿轮剩余寿命概 率密度扎k|yl:kak|yl:k)，式（10)为tk时刻基于观测值的齿轮剩余寿命概率密度理论推导结 果，
[00创
（I巧
[0043] 其中
，lk为tk时刻齿轮剩余寿命， 入为失效阔值
.
[0044] 7)将步骤6)所得打kivi:k(lk yi:k)代入式（11)，计算tk时刻齿轮预测剩余寿命为：
[0045]
…）
[0046] 8)判断齿轮更新状态知是否达到失效阔值A，如未达到则获取其下一时刻的健康 指标观测值，令k = k+l，重复步骤5.2)至步骤7)，直到％达到失效阔值停止。
[0047] 本发明的有益效果:充分考虑实际齿轮退化过程中广泛存在的时变性、个体差异 性、非线性和测量不确定性四方面随机因素，建立贴合其实际退化过程的四因素随机退化 模型，并推导出该模型下齿轮剩余寿命基于观测值的概率密度函数;采用MLE方法对该模型 参数进行初始化，避免参数初始化的盲目性；引入双重更新方法对模型状态及参数进行更 新，获得更准确的更新结果，进而很大程度上提高了齿轮剩余寿命的预测精度。
【附图说明】
[004引图1为本发明流程图。
[0049] 图2为齿轮加速退化寿命实验台结构图。
[0050] 图3为测试齿轮全寿振动信号：图3(a)为第一组测试齿轮全寿振动信号；图3(b)为 第二组测试齿轮全寿振动信号。
[0051] 图4为测试齿轮RMS指标和起始预测时刻：图4(a)为第一组测试齿轮的RMS指标和 起始预测时刻；图4(b)为第二组测试齿轮的RMS指标和起始预测时刻。
[0052] 图5为五种方法对测试齿轮健康状态的估计结果：图5(a)为五种方法对第一组测 试齿轮健康状态估计结果;图5 (b)为五种方法对第二组测试齿轮健康状态估计结果。
[0053] 图6为五种方法对测试齿轮剩余寿命的预测结果：图6(a)为五种方法对第一组测 试齿轮剩余寿命预测结果;图6(b)为五种方法对第二组测试齿轮剩余寿命预测结果。
【具体实施方式】
[0054] 下面结合附图和实施例对本发明做进一步详细描述。
[0055] 如图1所示，一种双重更新的四因素随机退化模型齿轮寿命预测方法，包括W下步 骤：
[0056] 1)使用监测系统采集齿轮振动信号，提取健康指标形成观测值向量y=[yi，...， ys ]T，S为独立观测值个数；
[0057] 2)构建四因素随机退化模型：
[005引
（1)
[0059] 其中，
[0060] Xk-tk时刻齿轮的真实退化状态；
[0061] 0k-tk时刻衰退系数，表征个体差异性，为随机变量，服从分布\ (AV讀）；
[0062] b-模型阶数，表征非线性；
[00创 Wk-I-过程噪声，服从分布N(0,02(tk-tk-i))，W k-i = 0(B(tk)-B(tk-i))，B(tk)为柄； 准布朗运动，地(tk)表征时变性；
[0064] Yk-tk时刻齿轮健康指标观测值；
[0065] V-观测噪声，服从分布N(0, 丫 2);
[0066] 3)采用极大似然估计(MLE)方法对式（1)四因素随机退化模型的参数进行初始化， 具体步骤包括：
[0067] 3.1)*1^寸刻齿轮健康指标观测值71<~^^，丫2)，根据步骤1)所得观测值向量7计 算获得对数似然函数：
[0068]
(2)
[0069] 其中，@ =片.A,.如/，yf，为四因素随机退化模型的未知参数向量
Q = 〇2m+ 丫 2ls，M=[min{ti,tj}]i《i,j《s，Is为S阶单位矩阵；
[0070] 3.2)计算式(2)关于叫和〇|的偏导，极大化似然函数：
[0071] 。）
[0072] 程组(3)得到：
[007；3] (令
[0074] 3.3)将式(4)计算结果回代至式(2)，简化对数似然函数为：
[0075]
贷
[0076] 3.4)通过Matlab提供的寻优命令对式(5)进行寻优获得t和O，回带至式(4)求解 出叫和诚，进而根据关系沪=打7请，F2=对封，揖: = 砖.,Q =O2M+丫 2Is,求解出未知参数 向量0 =护，心诚6,巧f的极大似然估计，作为式（1)四因素随机退化模型的初始参数
[0077] 4)将式(1)四因素随机退化模型表示为式(6)形式，并带入初始参数0 0:
[007引 （乂）
[0079]
[0080] 5)引入双重更新方法，对四因素随机退化模型的状态和参数进行实时更新，该方 法由两次相互补充、嵌套进行的更新过程组成：
[00川 5.1)外层更新开始，令k=0,对正态分布WUiA)，即WUiPu)进行随机采样，获得由N个初 始粒子构成粒子集柄U并娜应的粒子权值集合为;<Uv，其中如，iV「户1三；]，、培二1;
[0082] 5.2)里层更新开始，令k二k+1，将Pk-I I k-i和N个未更新粒子= 1,2,,: ? 入式 (7)迭代更新，计算获得Pk I k和論，i = 1，2,….W ;
[0083]
[0084] '，'
[0085] 5.3)里层更新结束，由步骤5.2)所得？^1<和4*，^ = 1，2，...，^^形成正态分布列 W技|&，巧|山=：1，2，...，^，对每个分布进行一次采样获得新的粒子集切心,、，，各粒子权值更新 如下：
[0089]
役)
[0086] 御
[0087]
[0088] 5.4)由步骤5.3)所得粒子及权值计算tk时刻齿轮更新状态和参数：
[0090] 外层更新结束；
、方矢义,' -____ ............ .............-[0094] 7)将步骤6)所得fLk|yi:k(lk|yi:k)代入式（11)，计算tk时刻齿轮预测剩余寿命为：[0095]
帅
[0091] 6)将步骤5.4)的更新结果代入式（10)，获得tk时刻基于观测值的齿轮剩余寿命概 率密度扎k|yl:kak|yl:k)，式（10)为tk时刻基于观测值的齿轮剩余寿命概率密度理论推导结 果，
[oow] ㈱
[0093] :剩余寿命，
[0096] 8)判断齿轮更新状态％是否达到失效阔值A，如未达到则获取其下一时刻的健康 指标观测值，令k = k+l，重复步骤5.2)至步骤7)，直到4达到失效阔值停止。
[0097] 实施例:从齿轮加速退化寿命实验台上采集齿轮加速退化实验数据对本发明进行 验证，实验台结构图如图2所示，该实验台采用与FZG标准减速机相似的背靠背式结构，主体 包括一个驱动电机和两个齿轮箱，齿轮箱分别为陪侍齿轮箱和测试齿轮箱，前者安装斜齿 轮，后者安装直齿圆柱齿轮。测试齿轮安装位置为测试齿轮箱齿轮晒合对中小齿轮，小齿轮 轴上安装有联轴器实现扭矩加载。实验过程中，驱动电机转速为1460rpm，负载扭矩约为 1300Nm。加速度传感器安装在测试齿轮箱箱体上，采用间歇采样的形式，采样频率为20曲Z， 采样持续时间为6s，即数据长度为120000。实验共进行两组，两组测试齿轮上均未提前加工 故障，实验结束后在第一组测试齿轮上观察到磨损、点蚀和裂纹故障，在第二组测试齿轮上 观察到严重裂纹故障。加速度传感器采集得到的两组测试齿轮的全寿命振动信号如图3所 /J、- O
[0098] 分别提取两组振动测试信号的RMS值作为健康指标，并选择起始预测时刻如图4所 示，分别为2088s和2098s。为获得当前时刻的剩余寿命预测结果，将自起始预测时刻至当前 时刻之间的RMS值作为观测值，输入本发明四因素随机退化模型，使用MLE方法对四因素随 机退化模型参数进行初始化，循环进行状态和参数双重更新、计算剩余寿命近似概率分布、 计算当前时刻齿轮预测寿命、提取下一个健康指标，直到齿轮的状态估计值达到预先设定 的失效阔值。
[0099] 本发明方法概括为"非线性退化模型+KPF更新"方法(M4)，为验证其有效性，引入 其他四种寿命预测方法与之进行比较，分别为："线性退化模型+KF更新"方法(Ml)、"非线性 退化模型+KF更新"方法(M2 )、"非线性退化模型+PF更新"方法(M3)和叩ari S-Erdogan模型+ KPF更新"方法(M5)。分别应用五种方法对两组齿轮进行剩余寿命预测，状态估计结果和寿 命预测结果如图5和图6所示。从图5中可W看出，M3及本发明方法的状态估计结果曲线更加 平滑，说明两方法较其他=种方法具有更加抗噪声干扰能力；从图6中可W看出，随着预测 时刻接近齿轮寿命终点，本发明方法的剩余寿命预测结果能够更快速的收敛于真实寿命。
[0100] 为定量验证本发明方法较于W上四种方法有更优良的剩余寿命预测表现，分别计 算五种方法的累计相对精度(CRA)，结果如表1所示，可W看出，本发明方法在两组齿轮实例 中均取得了最大值，验证了其更高的寿命预测精确度。
[0101] 亲1
[0102]
[0103]本发明提出一种双重更新的四因素随机退化模型齿轮寿命预测方法，适用对象为 齿轮，但齿轮与其他常见旋转零部件退化机理相似，因此本发明方法作出调整后同样适用 于与本发明模型描述随机退化特征相似的其他零部件或机械设备。此外，由于加速退化实 验通过增大负载加速齿轮退化并未改变其退化机理，本发明方法同样适用于实际应用中缓 慢随机退化齿轮的剩余寿命预测，通过更换模型输入指标、修改模型参数或拓展模型适用 条件等都未脱离本发明的构思，也应视为本发明的保护范围。
【主权项】
1. 一种双重更新的四因素随机退化模型齿轮寿命预测方法，其特征在于，包括以下步 骤： 1) 使用监测系统采集齿轮振动信号，提取健康指标形成观测值向量y=[yi，...，ys]T，s 为独立观测值个数； 2) 构建四因素随机退化模型：(1) 其中， Xk-tk时刻齿轮的真实退化状态； 0k一tk时刻衰退系数，表征个体差异性，为随机变量，服从分布〃(/%_，1 b-模型阶数，表征非线性； ω k-1-过程噪声，服从分布Ν(0，σ2(tk-tk-1))，ω k-：l = σ (B(tk)-B(tk-1))，B(tk)为标准布 朗运动，〇B(tk)表征时变性； yk-tk时刻齿轮健康指标观测值； v-观测噪声，服从分布Ν(0, γ2); 3) 采用极大似然估计(MLE)方法对式（1)四因素随机退化模型的参数进行初始化，具体 步骤为： 3.1) tk时刻齿轮健康指标观测值yk~N(Xk，γ2)，根据步骤1)所得观测值向量y计算获得 对数似然函数：(°) Z Z Z 其中，(^[^外^太/^^为四因素随机退化模型的未知参数向量乂吨丨tbJ 5 Ω = σ2Μ+ γ 2Is，M=[min{ti，tj}]issi,jsSs，Is为S阶单位矩阵； 3.2) 计算式(2)关于的偏导，极大化似然函数：(3) 其中；W/<，f=r2/W，β = 求解方程组(3)得到：(4) 3.3) 将式(4)计算结果回代至式(2)，简化对数似然函数为：(5) 3.4) 通过1^〖1&13提供的寻优命令对式(5)进行寻优获得〖和石,回带至式(4)求解出以(3和 ，进而根据关系及2 =σ2/σ〗，产=r74，Ω = Ω/σ〗，Ω =σ2Μ+ γ 2IS，求解出未知参数向量 的极大似然估计，作为式（1)四因素随机退化模型的初始参数4) 将式(1)四因素随机退化模型表示为式(6)形式，并带入初始参数Θ0:(6) 5) 引入双重更新方法，对四因素随机退化模型的状态和参数进行实时更新，该方法由 两次相互补充、嵌套进行的更新过程组成： 5 · 1)外层更新开始，令k=0,对正态分布，即#(4Pe)进行随机采样，获得由Ν个初始粒 子构成粒子集i<} 、，对应的粒子权值集合为，其中：5.2)里层更新开始，令k = k+Ι，将Ρη |k-dPN个未更新粒子U = 1，2,…，，带入式(7) 迭代更新，计算获得Pk 11<和％丨*；.，2'.= 1，2,…，況?5. 3)里层更新结束，由步骤5.2)所得Pk|dP = 2,...,7V，形成正态分布列 iV('，/^，^^..，〃，对每个分布进行一次采样获得新的粒子集以^^各粒子权值更新 如下：5.4)由步骤5.3)所得粒子及权值计算tk时刻齿轮更新状态和参数：(9) 外层更新结束； 6) 将步骤5.4)的更新结果代入式（10)，获得tk时刻基于观测值的齿轮剩余寿命概率密 度I凡^，式(10)为t k时刻基于观测值的齿轮剩余寿命概率密度理论推导结果， (10) .^yukK，k\u；kij-ek ~T ^k) L , 」.，. 其1i1 ^ik为tk时刻齿轮剩余寿命，λ为失 效阈值7) 将步骤6)所得(4 |崧)代入式(11)，计算tk时刻齿轮预测剩余寿命为：(Π) 8) 判断齿轮更新状态毛是否达到失效阈值λ，如未达到则获取其下一时刻的健康指标观 测值，令k = k+l，重复步骤5.2)至步骤7)，直到4达到失效阈值停止。
【文档编号】G06F17/50GK105956236SQ201610257102
【公开日】2016年9月21日
【申请日】2016年4月22日
【发明人】雷亚国, 林京, 周昕, 李乃鹏
【申请人】西安交通大学