CN104598734B - 滚动轴承集成期望最大化和粒子滤波的寿命预测方法 - Google Patents
滚动轴承集成期望最大化和粒子滤波的寿命预测方法 Download PDFInfo
- Publication number
- CN104598734B CN104598734B CN201510033397.9A CN201510033397A CN104598734B CN 104598734 B CN104598734 B CN 104598734B CN 201510033397 A CN201510033397 A CN 201510033397A CN 104598734 B CN104598734 B CN 104598734B
- Authority
- CN
- China
- Prior art keywords
- sigma
- beta
- theta
- rolling bearing
- delta
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
滚动轴承集成期望最大化和粒子滤波的寿命预测模型,首先采用峭度指标对轴承健康状态进行实时监测,确定寿命预测起始时刻;当满足预测起始条件后,采用有效值对轴承剩余寿命进行预测;在预测阶段,采用期望最大化方法对模型参数进行评估,同时采用粒子滤波方法对轴承状态进行评估,通过对模型参数和轴承状态的准确评估,提高剩余寿命预测精度,本发明能够实现对模型参数和滚动轴承状态的准确评估,并且在滚动轴承寿命预测中表现出了比传统指数模型更好的预测效果。
Description
技术领域
本发明涉及滚动轴承剩余寿命预测技术领域,具体涉及滚动轴承集成期望最大化和粒子滤波的寿命预测模型。
背景技术
滚动轴承广泛应用于机械设备中,由于复杂多变的工况,滚动轴承故障时有发生。为了保证设备的安全运行,传统定期维修策略需要投入大量人力物力对滚动轴承进行定期检修,对存在故障或安全隐患的滚动轴承进行更换。滚动轴承从发生故障到完全失效往往经历很长的衰退期,对存在早期故障的滚动轴承进行更换势必会大大缩短滚动轴承的有效服役周期,造成资源的浪费。采用预防性维修策略,对滚动轴承的剩余寿命进行有效预测,可以延长滚动轴承的有效服役周期,提高经济效益。因此如何对滚动轴承进行有效剩余寿命预测正受到国内外学者的广泛关注。
基于模型的滚动轴承剩余寿命预测方法,试图采用物理或数学模型对滚动轴承的衰退趋势进行描述,并根据观测数据对模型参数进行评估,以预测滚动轴承未来时刻健康衰退趋势和剩余寿命。美国普渡大学Gebraeel等人提出的指数预测模型与滚动轴承加速衰退趋势相适应,因此在滚动轴承剩余寿命预测中得到广泛应用。在指数预测模型中,模型参数和状态评估的准确性是影响模型预测精度的两个关键因素。Gebraeel最先提出的指数预测模型采用贝叶斯框架对模型参数进行评估。为了提高参数评估的准确性,清华大学司小胜等人对指数模型进行改进,采用期望最大化和贝叶斯更新相结合的方法对模型参数进行评估,得到了较好的参数评估效果。但是以上研究工作都存在一个弊端,即直接将滚动轴承观测值作为状态值进行寿命预测,没有对滚动轴承的真实健康状态进行评估。实际上,观测值与滚动轴承真实状态之间存在一定差异,而且观测值中有大量随机噪声的干扰。因此,传统指数预测模型只对模型参数进行评估而忽略了对滚动轴承状态的评估,导致模型预测精度降低。
发明内容
为了克服上述现有技术存在的缺点,本发明的目的在于提供滚动轴承集成期望最大化和粒子滤波的寿命预测模型,采用期望最大化方法对模型参数进行评估,采用粒子滤波对滚动轴承状态进行评估,以提高指数模型的预测精度。
为了达到上述目的,本发明采取的技术方案为:
滚动轴承集成期望最大化和粒子滤波的寿命预测模型,包括以下步骤:
1)实时监测并采集滚动轴承振动信号,计算滚动轴承峭度指标和有效值;
2)计算滚动轴承健康时刻峭度指标的均值μ和标准差σ,以确定健康状态下峭度指标的3σ区间[μ-3σ,μ+3σ];
3)判断滚动轴承连续l+1个时刻的峭度指标{mp+i}i=0:l是否满足{|mp+i-μ|>3σ}i=0:l,如果满足该条件,则确定mp所对应的时刻tp为寿命预测的起始时刻;
4)从寿命预测起始时刻开始,将滚动轴承振动信号有效值带入衰退模型:
其中,xi为ti时刻状态值,是已知常数,θ′,β′和σ为三个未知参数,σB(ti)~N(0,σ2ti)服从布朗运动,对上式求对数得到以下变形形式:
其中,θ=ln(θ′)服从分布β=β′-σ2/2服从分布且规定以简化计算;
对模型参数进行初始化,得初始参数μ0=μθ,0,μ1=μβ,0和
5)从概率密度函数中进行随机采样,得到Ns个初始粒子粒子权值为其中Δt=tj-tj-1为时间间隔;
6)tk时刻得到观测序列S1:k={s1,s2,...,sk},在参数已知条件下,观测序列的条件概率密度为:
根据贝叶斯理论得到参数θ和β的联合概率密度函数为:
由此得到参数更新公式如下:
采用以上公式对参数μθ,k,μβ,k和进行更新;
7)为待评估的模型参数,计算似然函数:
似然函数的期望值为:
其中为第i次评估的结果,计算满足条件的参数为第i+1次评估的结果,
8)将衰退模型改写为以下形式:
sk=sk-1+βΔt+σW(Δt), (9)
其中W(Δt)=B(tk)-B(tk-1),由此得到重要密度函数为:
从以上重要密度函数中进行重要性采样,得到粒子集
9)采用tk时刻的观测值Sk对粒子权值进行更新,
采用下式计算有效粒子数,
如果有效粒子数小于阈值NT,需要根据粒子权值大小进行重采样,得到新的粒子集使其满足粒子权值重置为
10)采用粒子集对滚动轴承状态进行评估,
然后将状态评估结果带入下式对滚动轴承剩余寿命的概率密度函数进行预测,
其中lk为tk时刻的剩余寿命,γ为滚动轴承失效阈值。
本发明利用期望最大化方法在参数评估方面的优势和粒子滤波方法在状态评估方面的优势,将两种方法集成到指数预测模型中,对滚动轴承的剩余寿命进行预测,克服了原有指数模型只注重参数评估而忽视了状态评估的弊端,采用粒子滤波方法对滚动轴承的健康状态进评估,状态评估结果能更好地反映滚动轴承的真实衰退趋势,通过同时对模型参数和滚动轴承状态的准确评估,提高了剩余寿命的预测精度。
附图说明
图1为滚动轴承集成期望最大化和粒子滤波的寿命预测模型流程图。
图2为PRONOSTIA实验台整体结构图。
图3为两个测试滚动轴承全寿命期内的振动信号。
图4为两个测试滚动轴承峭度指标和预测起始时刻选择结果图。
图5为两个测试滚动轴承有效值和对应的预测起始时刻。
图6为采用改进指数模型对两个测试滚动轴承有效值从预测起始时刻进行评估的结果。
图7为采用传统指数模型和改进指数模型对两个测试滚动轴承剩余寿命从预测起始时刻开始进行预测的结果对比图。
具体实施方式
下面结合附图和实施例对本发明做进一步详细描述。
如图1所示,滚动轴承集成期望最大化和粒子滤波的寿命预测模型,包括以下步骤:
1)实时监测并采集滚动轴承振动信号,计算滚动轴承峭度指标和有效值;
2)计算滚动轴承健康时刻峭度指标的均值μ和标准差σ,以确定健康状态下峭度指标的3σ区间[μ-3σ,μ+3σ];
3)判断滚动轴承连续l+1个时刻的峭度指标{mp+i}i=0:l是否满足{|mp+i-μ|>3σ}i=0:l,如果满足该条件,则确定mp所对应的时刻tp为寿命预测的起始时刻;
4)从寿命预测起始时刻开始,将滚动轴承振动信号有效值带入衰退模型:
其中,xi为ti时刻状态值,是已知常数,θ′,β′和σ为三个未知参数,σB(ti)~N(0,σ2ti)服从布朗运动,对上式求对数得到以下变形形式:
其中,θ=ln(θ′)服从分布β=β′-σ2/2服从分布且规定以简化计算;
对模型参数进行初始化,得初始参数μ0=μθ,0,μ1=μβ,0和
5)从概率密度函数中进行随机采样,得到Ns个初始粒子粒子权值为其中Δt=tj-tj-1为时间间隔;
6)tk时刻得到观测序列S1:k={s1,s2,...,sk},在参数已知条件下,观测序列的条件概率密度为:
根据贝叶斯理论得到参数θ和β的联合概率密度函数为:
由此得到参数更新公式如下:
采用以上公式对参数μθ,k,μβ,k和进行更新;
7)为待评估的模型参数,计算似然函数:
似然函数的期望值为:
其中为第i次评估的结果,计算满足条件的参数为第i+1次评估的结果,
8)将衰退模型改写为以下形式:
sk=sk-1+βΔt+σW(Δt), (9)
其中W(Δt)=B(tk)-B(tk-1),由此得到重要密度函数为:
从以上重要密度函数中进行重要性采样,得到粒子集
9)采用tk时刻的观测值Sk对粒子权值进行更新,
采用下式计算有效粒子数,
如果有效粒子数小于阈值NT,需要根据粒子权值大小进行重采样,得到新的粒子集使其满足粒子权值重置为
10)采用粒子集对滚动轴承状态进行评估,
然后将状态评估结果带入下式对滚动轴承剩余寿命的概率密度函数进行预测,
其中lk为tk时刻的剩余寿命,γ为滚动轴承失效阈值。
下面结合实施例对本发明做详细描述。
实施例:采用PRONOSTIA实验台上采集的滚动轴承加速寿命实验数据对本发明进行验证。
PRONOSTIA实验台如图2所示,该实验台通过对滚动轴承气压加载,使滚动轴承在高负荷条件下工作,可以在数小时内实现滚动轴承从正常状态退化到完全失效。实验过程中,滚动轴承转速为1800rpm,负载为4kN。采用加速度传感器对滚动轴承振动信号进行采样,采样频率为25.6kHz,数据长度为2560,每次采样持续时间为0.1s,采样间隔为10s。当振动幅值超过20m/s2时,滚动轴承完全失效。两组实验滚动轴承的全寿命振动信号如图3所示。
从振动信号中提取峭度指标,对寿命预测起始时刻进行选择,选择过程如图4所示,由健康历史数据确定滚动轴承峭度指标的3σ区间,然后对滚动轴承峭度指标进行实时监测,当峭度指标mp超出3σ区间即|mp-μ|>3σ,说明此时滚动轴承状态异常,滚动轴承异常状态可能由滚动轴承故障引起,也可能由随机噪声引起,为了排除随机噪声对峭度指标带来的干扰,连续观测tp时刻之后的l个时刻的峭度指标{mp+i}i=1:l是否也超出了3σ区间,即判断{|mp+i-μ|>3σ}i=1:l,如果满足以上条件,说明滚动轴承的异常状态确实由故障引起。因此选择tp时刻为寿命预测的起始时刻。
确定预测起始时刻之后,从该时刻开始,将滚动轴承有效值带入衰退模型,对滚动轴承剩余寿命进行预测,两滚动轴承的有效值如图5所示,可以看出,所选的时刻与滚动轴承加速衰退的起始时刻相近,作为滚动轴承的预测起始时刻较为合适,本发明采用粒子滤波方法对滚动轴承有效值进行评估,评估结果如图6中所示,从图中可以看出,由于随机噪声的影响,有效值局部波动性较强,而粒子滤波方法评估结果可以有效抑制随机噪声的干扰,反映出滚动轴承真实的加速衰退趋势,因此本发明通过采用粒子滤波方法,实现了对滚动轴承健康状态的有效评估。
图7显示了滚动轴承剩余寿命预测结果,从图中可以看出,在寿命预测前期,传统指数模型和改进指数模型的预测效果都不理想,随着时间的推移,两种方法都逐渐收敛到真实值,但是改进指数模型的收敛速度更快,预测精度更高,其原因为,在刚开始预测时,由于观测数据点数较少,模型参数无法得到有效评估,因此引起较大的预测误差,当观测点数逐渐增多,模型参数获得较准确的评估结果,此时状态评估结果便成为影响寿命预测结果的主要因素,由于改进指数模型采用粒子滤波方法实现了对滚动轴承健康状态的准确评估,因此其预测效果明显好于传统指数模型。
通过采用滚动轴承加速寿命实验数据验证了本发明在滚动轴承剩余寿命预测中的优势,由于传统指数模型只是对模型参数进行评估,而滚动轴承观测值中包含大量的随机噪声,影响了剩余预测结果的精度,本发明对传统指数模型进行改进,采用粒子滤波方法对滚动轴承健康状态进行有效评估,从而提高了指数模型对滚动轴承剩余寿命的预测精度。
本发明所提出的滚动轴承集成期望最大化和粒子滤波的寿命预测模型,并不只局限于滚动轴承的剩余寿命预测,还可以应用于其他机械电子产品的剩余寿命预测问题。大量研究工作证明,本方法适用于各类具有指数衰退形式的机电产品的剩余寿命预测。实施者只需对本方法相应步骤进行适当调整,以适应不同产品的应用需求。另外,本发明提供了一种同时对模型参数和产品健康状态进行准确评估的思路,应当指出,在不脱离本发明构思的前提下,所做的调整和变形,也应视为本发明的保护范围。
Claims (1)
1.滚动轴承集成期望最大化和粒子滤波的寿命预测方法,其特征在于,包括以下步骤:
1)实时监测并采集滚动轴承振动信号,计算滚动轴承峭度指标和有效值;
2)计算滚动轴承健康时刻峭度指标的均值μ和标准差v,以确定健康状态下峭度指标的3v区间[μ-3v,μ+3v];
3)判断滚动轴承连续l+1个时刻的峭度指标{mp+i}i=0:l是否满足{|mp+i-μ|>3v}i=0:l,如果满足该条件,则确定mp所对应的时刻tp为寿命预测的起始时刻;
4)从寿命预测起始时刻开始,将滚动轴承振动信号有效值代入衰退模型:
其中,xi为ti时刻状态值,是已知常数,θ′,β′和σ为三个未知参数,σB(ti)~N(0,σ2ti)服从布朗运动,对上式求对数得到以下变形形式:
其中,θ=ln(θ′)服从分布β=β′-σ2/2服从分布且规定以简化计算;
对模型参数进行初始化,得初始参数μ0=μθ,0,μ1=μβ,0和
5)从概率密度函数中进行随机采样,得到Ns个初始粒子粒子权值为其中Δt=tj-tj-1为时间间隔;
6)tk时刻得到观测序列S1:k={s1,s2,...,sk},在参数已知条件下,观测序列的条件概率密度为:
根据贝叶斯理论得到参数θ和β的联合概率密度函数为:
由此得到参数更新公式如下:
采用以上公式对参数μθ,k,μβ,k和进行更新;
7)为待评估的模型参数,计算似然函数:
似然函数的期望值为:
其中为第i次评估的结果,计算满足条件的参数为第i+1次评估的结果,
8)将衰退模型改写为以下形式:
sk=sk-1+βΔt+σW(Δt), (9)
其中W(Δt)=B(tk)-B(tk-1),由此得到重要密度函数为:
从以上重要密度函数中进行重要性采样,得到粒子集
9)采用tk时刻的观测值Sk对粒子权值进行更新,
采用下式计算有效粒子数,
如果有效粒子数小于阈值NT,需要根据粒子权值大小进行重采样,得到新的粒子集使其满足粒子权值重置为
10)采用粒子集对滚动轴承状态进行评估,
然后将状态评估结果代入下式对滚动轴承剩余寿命的概率密度函数进行预测,
其中lk为tk时刻的剩余寿命,γ为滚动轴承失效阈值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510033397.9A CN104598734B (zh) | 2015-01-22 | 2015-01-22 | 滚动轴承集成期望最大化和粒子滤波的寿命预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510033397.9A CN104598734B (zh) | 2015-01-22 | 2015-01-22 | 滚动轴承集成期望最大化和粒子滤波的寿命预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104598734A CN104598734A (zh) | 2015-05-06 |
CN104598734B true CN104598734B (zh) | 2017-05-17 |
Family
ID=53124513
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510033397.9A Active CN104598734B (zh) | 2015-01-22 | 2015-01-22 | 滚动轴承集成期望最大化和粒子滤波的寿命预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104598734B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104881580B (zh) * | 2015-05-27 | 2017-08-29 | 中国人民解放军国防科学技术大学 | 一种卫星驱动机构在轨健康状态指标提取及寿命预测方法 |
CN105224792B (zh) * | 2015-09-21 | 2018-04-06 | 河南科技大学 | 一种滚动轴承性能保持可靠性的预测方法 |
CN105653851B (zh) * | 2015-12-27 | 2018-09-21 | 北京化工大学 | 基于分阶段物理模型及粒子滤波的滚动轴承剩余寿命预测方法 |
CN105740625B (zh) * | 2016-01-31 | 2018-02-23 | 太原科技大学 | 一种齿轮的实时剩余寿命预测方法 |
CN105956236B (zh) * | 2016-04-22 | 2019-03-12 | 西安交通大学 | 一种双重更新的四因素随机退化模型齿轮寿命预测方法 |
CN106021719A (zh) * | 2016-05-19 | 2016-10-12 | 浙江工业大学 | 一种基于无迹卡尔曼滤波算法的轴承剩余寿命预测方法 |
CN106934126B (zh) * | 2017-02-28 | 2020-03-31 | 西安交通大学 | 基于循环神经网络融合的机械零部件健康指标构造方法 |
CN106934125B (zh) * | 2017-02-28 | 2020-02-18 | 西安交通大学 | 一种梯形噪声分布的指数模型机械设备剩余寿命预测方法 |
CN108304348B (zh) * | 2018-01-31 | 2021-06-18 | 浙江工业大学 | 一种基于二元维纳过程的轴承剩余寿命预测方法 |
CN109212966B (zh) * | 2018-08-14 | 2020-04-10 | 西安交通大学 | 一种多工况动态基准化的机械设备剩余寿命预测方法 |
CN110231165B (zh) * | 2019-05-31 | 2020-03-31 | 西安交通大学 | 一种基于期望差异约束置信网络的机械设备故障诊断方法 |
CN112683535B (zh) * | 2021-01-14 | 2022-04-12 | 大连理工大学 | 基于多阶段维纳过程的轴承寿命预测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103955750A (zh) * | 2014-04-04 | 2014-07-30 | 西安交通大学 | 基于特征融合和粒子滤波的滚动轴承剩余寿命预测方法 |
CN103954450A (zh) * | 2014-05-19 | 2014-07-30 | 重庆交通大学 | 基于主成分分析的轴承寿命退化性能评估指标构建方法 |
-
2015
- 2015-01-22 CN CN201510033397.9A patent/CN104598734B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103955750A (zh) * | 2014-04-04 | 2014-07-30 | 西安交通大学 | 基于特征融合和粒子滤波的滚动轴承剩余寿命预测方法 |
CN103954450A (zh) * | 2014-05-19 | 2014-07-30 | 重庆交通大学 | 基于主成分分析的轴承寿命退化性能评估指标构建方法 |
Non-Patent Citations (3)
Title |
---|
A particle filtering-based approach for remaining useful life predication of rolling element bearings;Naipeng Li et al.;《IEEE Conference on Prognostics and Health Management》;20141231;第1-8页 * |
Monte Carlo-based filtering for fatigue crack growth estimation;F. Cadini et al.;《Probabilistic Engineering Mechanics》;20081120;第367-373页 * |
基于粒子滤波的滚动轴承寿命预测方法;雷亚国等;《2014年全国机械行业可靠性技术学术交流会论文集》;20140802;第194-199页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104598734A (zh) | 2015-05-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104598734B (zh) | 滚动轴承集成期望最大化和粒子滤波的寿命预测方法 | |
CN110414155B (zh) | 一种带有单测点的风机部件温度异常检测和报警方法 | |
CN112304613B (zh) | 一种基于特征融合的风电机组发电机轴承预警方法 | |
CN105241660B (zh) | 基于健康监测数据的高铁大型桥梁性能测试方法 | |
CN110414154B (zh) | 一种带有双测点的风机部件温度异常检测和报警方法 | |
US20190285517A1 (en) | Method for evaluating health status of mechanical equipment | |
CN105300692B (zh) | 一种基于扩展卡尔曼滤波算法的轴承故障诊断及预测方法 | |
CN104792529A (zh) | 基于状态空间模型的滚动轴承寿命预测方法 | |
Hu et al. | Online performance assessment method for a model-based prognostic approach | |
CN103955750B (zh) | 基于特征融合和粒子滤波的滚动轴承剩余寿命预测方法 | |
EP2384971B1 (en) | Method of determining a maneuver performed by an aircraft | |
CN105834835A (zh) | 一种基于多尺度主元分析的刀具磨损在线监测方法 | |
CN107016235A (zh) | 基于多特征自适应融合的设备运行状态健康度评估方法 | |
EP2989705B1 (en) | Method and apparatus for defect pre-warning of power device | |
CN104598736B (zh) | 一种自适应多核组合相关向量机的滚动轴承寿命预测模型 | |
CN101799367A (zh) | 一种机电设备神经网络故障趋势预测方法 | |
CN105653851B (zh) | 基于分阶段物理模型及粒子滤波的滚动轴承剩余寿命预测方法 | |
CN110688617B (zh) | 风机振动异常检测方法及装置 | |
CN103217280A (zh) | 航空发动机转子剩余寿命的多变量支持向量机预测方法 | |
CN108896299A (zh) | 一种齿轮箱故障检测方法 | |
CN107132454A (zh) | 基于随机矩阵谱半径法的电网异常快速检测方法 | |
EP3480455B1 (en) | Wind turbine monitoring device, wind turbine monitoring method, wind turbine monitoring program, and storage medium | |
CN103438794B (zh) | 一种基于压电传感网络的工程结构冲击位置确定方法 | |
CN112966630B (zh) | 一种超/特高压gis开关设备运行状态评估方法、***及设备 | |
DE102016013406A1 (de) | Verfahren zum Betrieb eines Zustandsüberwachungssystems einer Schwingmaschine und Zustandsüberwachungssystem |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |