核估计和随机滤波集成的齿轮剩余寿命预测方法技术

一种基于核估计和随机滤波理论的齿轮剩余寿命预测方法，属于机械可靠性技术领域，具体实施步骤如下：1、利用加速度传感器对主试齿轮箱内齿轮退化状态实时监测；2、对齿轮退化状态进行特征提取；3、利用核函数对数据的分布不做任何假定,从数据样本本身出发的特点对齿轮连续退化状态的概率密度函数进行非参数估计，得到齿轮基于实时状态监测数据的退化状态概率密度函数；4、利用实时状态监测数据来更新随机滤波递推模型参数，建立基于核估计和随机滤波相结合的预测模型；5、通过核估计和随机滤波相结合的预测模型预测齿轮剩余寿命。优点是可以有效地预测齿轮退化状态及实时剩余寿命，为齿轮预防性维修提供依据。

【技术实现步骤摘要】
核估计和随机滤波集成的齿轮剩余寿命预测方法
本专利技术属于机械可靠性设计领域，具体涉及一种齿轮剩余寿命预测方法，
技术介绍
齿轮是机械工业中应用广泛的机械设备传动系统中的关键部件。当齿轮发生断齿、齿面疲劳、胶合等故障时，常常会引起整个机械设备灾难性的破坏。以风力发电机组为例，齿轮故障率是整个风电机组中最高的，约占60％，而且其维护费用也较高，约占40％，因此，对齿轮提出合理有效的维修方案已成为风电行业急需解决的问题，而在整个维修方案制定过程中，齿轮的剩余寿命预测是重点和难点，随着信息传感设备的发展，对齿轮的运行状态进行实时监测，利用接收到的大量实时监测信息准确预测系统的退化状态及其剩余寿命，提供齿轮健康状态的关键信息，进而识别和管理故障的发生、规划维修活动，为更合理地制定基于状态的维护维修策略提供依据，目前，齿轮剩余寿命的预测方法分为四类：基于物理模型的预测方法、基于统计经验的预测方法、基于知识的预测方法和基于数据驱动的预测方法。现有预测方法存在下列问题：首先，现有预测方法需要进行状态退化模型结构假设，需要假设其作为判断依据的样本符合某种特定的模型结构，这些模型结构的假设与实际的物理模型之间常常存在较大的差距；其次，预测模型中所涉及到的参数估计问题，大多不能保证全局收敛；最后，由于齿轮处于变化的环境中，它的状态退化模型会发生改变，单一的预测模型不能适应环境的变化，需要多种预测模型相结合，
技术实现思路
本专利技术目的是提供一种核估计和随机滤波集成的齿轮剩余寿命预测方法，利用实时状态监测数据来更新随机滤波递推模型参数，可准确地预测齿轮的实时剩余寿命，本专利技术是这样实现的，其特征是包括以下实施步骤：步骤1、通过试验台架获取表征主试齿轮箱内齿轮状态的实时监测数据：采用如图2所示的试验台架，试验台架的中心距为a＝150mm。试验采用机械杠杆4加载，扭矩采用转矩传感器13#进行测量。试验过程中对齿轮的振动、加速度、温度和噪声通过接收传感器信号进行监测。如图3所示，主试齿轮箱1中为正反面交错搭接啮合的一对齿轮，齿轮的断齿状态等效为齿轮的失效。试验对主试齿轮箱1进行数据监测，试验过程中对齿轮的振动、加速度、温度和噪声通过接收传感器信号进行监测；试验共布置八个加速度传感器1#～8#,两个声音传感器9#和10#，一个温度传感器11#,一个转速传感器12#，一个转矩传感器13#。有四个加速度传感器1#～4#分别布置在主试齿轮箱轴承座的径向位置，有两个加速度传感器7#和8#分别布置在主试齿轮箱的轴向位置，有两个加速度传感器5#和6#分别布置在陪试齿轮箱轴承座的径向位置；有两个声音传感器9#和10#分别悬挂在主试齿轮箱和陪试齿轮箱正上方40cm处；一个温度传感器11#布置在主试齿轮箱内部；一个转速传感器12#布置在驱动电机3与陪试齿轮箱2联接轴的中部；一个转矩传感器13#布置在主试齿轮箱1与陪试齿轮箱2联接轴的中部。试验中依次加载八级载荷，八级载荷的大小在330～850N.m之间，每级载荷的加载时间为9～12个小时，在第八级载荷发生断齿。通过用加速度传感器4#记录齿轮的加速度数据，剩余寿命预测选取从第八级加载开始到断齿的加速度测试点整个时域信号进行分析。采样信息如下：采样频率为20～30kHz，每次采样持续50～70秒，每隔8～10分钟记录一次采样文本；步骤2、对主试齿轮箱内齿轮的退化状态进行特征提取，利用均方幅值对齿轮磨损退化性能进行衰退评估，对于在每次采样时间长度内，采样信号的均方幅值特征值表示为：式中：∑为求和号，m为每个采样周期的采样点数,yi为ti时刻齿轮的状态信息均方幅值，yj为齿轮每个采样周期的数据；步骤3、根据主试齿轮箱内齿轮的初始故障数据，利用核密度估计方法进行非参数估计，得出：式中：∑为求和号，f0(x0)为t0时刻的剩余寿命概率密度函数，h为平滑参数,K(u)为核函数，n为随机变量X的样本数，xi为i时刻齿轮的故障数据，选用高斯核函数：其中exp是以自然常数e为底的指数函数，判断窗宽的好坏通过积分均方误差来检验，遵循的原则为极小化渐进积分均方误差法AMISE，其中极小化渐进积分均方误差AMISE值为：式中：f(x)是真实值服从正态分布，f″(x)为对f(x)求二次导，为f(x)的估计值；要使渐进积分均方误差达到最小，必须选择恰当的h，即用求导的方法来求得最优窗宽h的值。令则可以得到最优窗宽如下：在取K(u)为高斯核函数时，利用Silverman大拇指法则，得到最优窗宽如下：为样本标准差(6)其中为样本标准差，选择核密度估计方法对fi(yi|xi)进行估计，根据齿轮的状态信息即：采样数据的特征值和剩余寿命的预测值去估计得：式中：fi(yi|xi)为ti时刻，yi在xi条件下的概率密度函数，α、β为平滑参数,K1(a,b)、K2(u)为核函数，选用高斯核函数：其中exp是以自然常数e为底的指数函数，ρ是随机变量a和b的相关系数、α、β为式fi(yi|xi)中平滑参数，由式(6)得出；步骤4、利用实时状态监测数据来更新随机滤波递推模型参数，建立基于核估计和随机滤波相结合的预测模型，由贝叶斯定理得：又因为：fi(yi|xi,Yi-1)＝fi(yi|xi)(9)所以：在ti时刻的剩余寿命等于ti-1时刻的剩余寿命减去时刻ti和时刻ti-1之间的间隔，即：根据式(11)和式(12)，得到随机滤波方程为：通过非参数估计得到初始寿命概率密度f0(x0)及状态退化概率密度fi(yi|xi)，将其代入递推公式(13)，得到剩余寿命概率密度函数fi(xi|Yi)；步骤5、通过核估计和随机滤波相结合的预测模型预测齿轮剩余寿命：最后得出齿轮平均剩余寿命为：本专利技术优点及积极效果是：本专利技术应用到齿轮的剩余寿命预测中，建立预测模型，可以对齿轮剩余寿命做出有效预测，附图说明图1为本专利技术实施例中齿轮剩余寿命预测方法流程图；图2为本专利技术实施例中试验台架示意图；图3为本专利技术实施例中齿轮正反面交错搭接啮合的示意图；图4为本专利技术实施例中特征值随监测时间变化曲线图；图5为本专利技术实施例中f0(x0)的概率密度曲线图；图6为本专利技术实施例中f1(y1|x1)的概率密度曲线图；图7为本专利技术实施例中f2(y2|x2)的概率密度曲线图；图8为本专利技术实施例中剩余寿命的概率密度函数fi(xi|Yi)图；图中：1-主试齿轮箱，2-陪试齿轮箱，3-驱动电机，4-机械杠杆，5-联轴器，a-主试齿轮箱与陪试齿轮箱的中心距；1#～8#-加速度传感器，9#、10#-声音传感器，11#-温度传感器，12#-转速传感器，13#-转矩传感器；具体实施方式下面结合附图对本专利技术实施例做进一步说明：如图2所示为本测试应用的试验台架示意图，测试方法步骤如下：步骤1、通过试验台架获取表征主试齿轮箱内齿轮状态的实时监测数据：采用如图2所示的试验台架，试验台架的中心距为a＝150mm。试验采用机械杠杆4加载，扭矩采用转矩传感器13#进行测量。试验过程中对齿轮的振动、加速度、温度和噪声通过接收传感器信号进行监测。如图3所示，主试齿轮箱1中为正反面交错搭接啮合的一对齿轮，齿轮的断齿状态等效为齿轮的失效。试验对主试齿轮箱1进行数据监测，试验过程中对齿轮的振动、加速度、温度和噪声通过接收传感器信号进行监本文档来自技高网...

【技术保护点】
1.一种核估计和随机滤波集成的齿轮剩余寿命预测方法，其特征在于实施步骤如下：步骤1、通过试验台架获取表征主试齿轮箱内齿轮状态的实时监测数据：采用试验台架的中心距为a＝150mm，试验采用机械杠杆(4)加载，扭矩采用转矩传感器(13#)进行测量，试验过程中对齿轮的振动、加速度、温度和噪声通过接收传感器信号进行监测，主试齿轮箱(1)中为正、反面交错搭接啮合的一对齿轮，齿轮的断齿状态等效为齿轮的失效，试验对主试齿轮箱(1)进行数据监测，试验过程中对齿轮的振动、加速度、温度和噪声通过接收传感器信号进行监测；试验共布置八个加速度传感器(1#～8#),两个声音传感器(9#、10#)，一个温度传感器(11#),一个转速传感器(12#)，一个转矩传感器(13#)，有四个加速度传感器(1#～4#)分别布置在主试齿轮箱轴承座的径向位置，有两个加速度传感器(7#、8#)分别布置在主试齿轮箱的轴向位置，有两个加速度传感器(5#、6#)分别布置在陪试齿轮箱轴承座的径向位置；有两个声音传感器(9#、10#)分别悬挂在主试齿轮箱和陪试齿轮箱正上方40cm处；一个温度传感器(11#)布置在主试齿轮箱内部；一个转速传感器(12#)布置在驱动电机(3)与陪试齿轮箱(2)联接轴的中部；一个转矩传感器(13#)布置在主试齿轮箱(1)与陪试齿轮箱(2)联接轴的中部，试验中依次加载八级载荷，八级载荷的大小在330～850N.m之间，每级载荷的加载时间为9～12个小时，在第八级载荷发生断齿，通过用加速度传感器(4#)记录齿轮的加速度数据，剩余寿命预测选取从第八级加载开始到断齿的加速度测试点整个时域信号进行分析；采样信息如下：采样频率为20～30kHz，每次采样持续50～70秒，每隔8～10分钟记录一次采样文本；步骤2、对主试齿轮箱内齿轮的退化状态进行特征提取，利用均方幅值对齿轮磨损退化性能进行衰退评估，对于在每次采样时间长度内，采样信号的均方幅值特征值表示为：...

【技术特征摘要】
1.一种核估计和随机滤波集成的齿轮剩余寿命预测方法，其特征在于实施步骤如下：步骤1、通过试验台架获取表征主试齿轮箱内齿轮状态的实时监测数据：采用试验台架的中心距为a＝150mm，试验采用机械杠杆(4)加载，扭矩采用转矩传感器(13#)进行测量，试验过程中对齿轮的振动、加速度、温度和噪声通过接收传感器信号进行监测，主试齿轮箱(1)中为正、反面交错搭接啮合的一对齿轮，齿轮的断齿状态等效为齿轮的失效，试验对主试齿轮箱(1)进行数据监测，试验过程中对齿轮的振动、加速度、温度和噪声通过接收传感器信号进行监测；试验共布置八个加速度传感器(1#～8#),两个声音传感器(9#、10#)，一个温度传感器(11#),一个转速传感器(12#)，一个转矩传感器(13#)，有四个加速度传感器(1#～4#)分别布置在主试齿轮箱轴承座的径向位置，有两个加速度传感器(7#、8#)分别布置在主试齿轮箱的轴向位置，有两个加速度传感器(5#、6#)分别布置在陪试齿轮箱轴承座的径向位置；有两个声音传感器(9#、10#)分别悬挂在主试齿轮箱和陪试齿轮箱正上方40cm处；一个温度传感器(11#)布置在主试齿轮箱内部；一个转速传感器(12#)布置在驱动电机(3)与陪试齿轮箱(2)联接轴的中部；一个转矩传感器(13#)布置在主试齿轮箱(1)与陪试齿轮箱(2)联接轴的中部，试验中依次加载八级载荷，八级载荷的大小在330～850N.m之间，每级载荷的加载时间为9～12个小时，在第八级载荷发生断齿，通过用加速度传感器(4#)记录齿轮的加速度数据，剩余寿命预测选取从第八级加载开始到断齿的加速度测试点整个时域信号进行分析；采样信息如下：采样频率为20～30kHz，每次采样持续50～70秒，每隔8～10分钟记录一次采样文本；步骤2、对主试齿轮箱内齿轮的退化状态进行特征提取，利用均方幅值对齿轮磨损退化性能进行衰退评估，对于在每次采样时间长度内，采样信号的均方幅值特征值表示为：式中：∑为求和号，m为每个采样周期的采样点数,yi为ti时刻齿...

【专利技术属性】
技术研发人员：石慧，白尧，张岩，
申请(专利权)人：太原科技大学，
类型：发明
国别省市：山西,14