基于改进高斯粒子滤波的涡扇发动机突变故障诊断方法技术

本发明专利技术公开了基于改进高斯粒子滤波的涡扇发动机突变故障诊断方法，涉及航空发动机故障诊断领域，能够实现突变故障的快速诊断，并且诊断结果噪声水平低，估测结果精度高。本发明专利技术包括：对发动机非线性数学模型注入气路部件突变故障；设计基于伪协方差的改进高斯粒子滤波算法；发动机非线性数学模型载入改进高斯粒子滤波算法，改进高斯粒子滤波算法对气路部件突变故障进行诊断,得到诊断结果。本发明专利技术采用伪协方差代替协方差、高斯采样代替重采样，减少了诊断时间、提高了诊断精度，能够实现发动机寿命周期内气路突变故障的快速诊断。

【技术实现步骤摘要】
基于改进高斯粒子滤波的涡扇发动机突变故障诊断方法
本专利技术涉及航空发动机故障诊断领域，尤其涉及了基于改进高斯粒子滤波的涡扇发动机突变故障诊断方法。
技术介绍
涡扇发动机结构复杂，工作环境恶劣，属于一种故障多发系统。据统计，涡扇发动机气路部件故障占涡扇发动机总体故障的90％以上，因此，实时检测发动机健康状况并进行气路性能分析是提高发动机安全可靠性的重要途径。发动机气路部件的效率变化系数、流量变化系数等健康参数是发动机气路故障的状态特征，而这些健康参数会直接导致转速、温度、压力等发动机测量参数的变化，因此发动机气路部件故障诊断主要是采用特定的算法根据测量参数的变化来估计出健康参数，从而分析发动机气路部件的健康状况。发动机故障诊断的方法主要有基于模型(如最小二乘方法，卡尔曼滤波方法)，基于数据(如神经网络，支持向量机)和基于知识(如专家系统)的方法。随着发动机部件级模型的精度以及计算机速度的提高，基于模型的故障诊断方法在工程实际中得到了广泛的研究和应用。卡尔曼滤波算法受限于线性高斯系统，一些改进的卡尔曼算法虽然适用于非线性问题，但依然依赖于高斯噪声的假设。而粒子滤波算法适用于含非高斯噪声的非线性系统，于是有学者提出使用粒子滤波算法进行气路部件健康参数的估计，并取得了较好的估计效果。标准粒子滤波用于发动机故障诊断主要存在以下两个方面的问题：(1)由于标准粒子滤波直接采用先验分布作为重要性密度函数，没有结合最新观测值，导致在突变故障下，不能实现突变故障的快速诊断；(2)由于发动机健康参数的维度高，滤波器的抽样率低，同时出于实时性考虑，粒子数少，所以有效粒子数低，导致诊断结果噪声水平高，估计精度低。因此，现有技术中，缺乏一种粒子滤波算法，能够实现突变故障的快速诊断，并且诊断结果噪声水平低，估测结果精度高。
技术实现思路
本专利技术提供了基于改进高斯粒子滤波的涡扇发动机突变故障诊断方法，能够实现突变故障的快速诊断，并且诊断结果噪声水平低，估测结果精度高。为达到上述目的，本专利技术采用如下技术方案：基于改进高斯粒子滤波的涡扇发动机故障诊断方法，包括：S1、对发动机非线性数学模型注入气路部件突变故障，发动机非线性数学模型，其中，根据发动机各部件间的流量连续、功率平衡、转子动力学原理建立共同方程，采用牛顿拉夫森迭代算法求解非线性方程组，获得发动机各工作截面的参数，依据发动机的气动热力学特性和部件特性数据，建立发动机非线性数学模型。S2、设计基于伪协方差的改进高斯粒子滤波算法；S3、发动机非线性数学模型载入改进高斯粒子滤波算法，改进高斯粒子滤波算法对气路部件突变故障进行诊断,得到诊断结果。进一步的，发动机非线性数学模型表征了以下故障量，效率变化程度系数SEi，流量变化程度系数SWi；效率变化程度系数和流量变化程度系数为健康参数，以反应发动机部件的健康状况。其中i＝1,2,3,4分别为风扇，高压压气机，高压涡轮，低压涡轮这四个旋转部件的序号。进一步的，气路部件突变故障包括：风扇部件突变故障，压气机部件突变故障，高压涡轮部件突变故障和低压涡轮部件突变故障。进一步的，S2具体包括：S21、k表示步数，在k＝0时，赋予初始值，再进入k＝1时刻，完成初始化，其中，为状态量估计的初始值，为粒子集伪协方差估计的初始值；S22、在k时刻，根据高斯采样从获得高斯采样粒子集其中为k-1时刻的等权重的高斯采样粒子，为k-1时刻的粒子加权均值，为k-1时刻的粒子集伪协方差；S23、在k时刻，从重要密度函数中依据重要性采样法抽样得到时刻更新粒子集其中，zk为k时刻的测量值，；S24、选取先验密度函数作为重要性密度函数，即先验密度函数为计算得到各粒子的重要性权重即为各粒子的后验似然概率分布，S25、将各粒子权值归一化，为归一化值，S26、估计状态值及粒子集伪协方差S27、判断是否达到最大步数kmax，否则令k＝k+1，返回S22，循环执行S22-S27。进一步的，发动机故障诊断中，状态估计值列向量由2个状态变量和8个健康参数组成，2个状态量分别为高压转速和低压转速，8个健康参数分别为风扇效率健康参数，风扇流量健康参数，压气机效率健康参数，压气机流量健康参数，高压涡轮效率健康参数，高压涡轮流量健康参数，低压涡轮效率健康参数，低压涡轮流量健康参数。在无突变故障发生的情况下，对于粒子加权均值向量和有其中为列向量第j行的值，j为列向量的行数。又因为分别表示发动机转速状态量，由于惯性作用，假设在相邻时刻的值不会突变，所以有于是在无突变故障发生的情况下，有粒子集的协方差∑k为粒子集伪协方差为当发动机无突变故障时伪协方差能够代替协方差；若第j,j＝3,4,…,10维的发动机健康参数发生突变故障且其他健康参数无突变故障，则相对第j行元素突变，一般情况下，归一化权值大的粒子必然聚集在附近，则的第j行元素突变，于是的对角线上第j个元素显著大于其他元素，下一时刻从中抽样获得的粒子比从中获得的粒子的第j个元素波动更加剧烈，而其他元素的波动程度基本保持为无突变故障情形。进一步的，在S3中，诊断包括：S31、将发动机的高度马赫数环境参数，燃油量、尾喷管面积控制参数，以及健康参数输入发动机非线性数学模型，发动机非线性数学模型计算得出预测观测参数；S32、发动机输出的观测参数为低压转子转速、高压转子转速、风扇出口温度、风扇出口压力、压气机出口温度、压气机出口压力、高压涡轮出口温度、高压涡轮出口压力。将观测参数和预测观测参数的残差，依据改进高斯粒子滤波算法，计算得出部件健康参数，实现发动机气路部件的突变故障诊断。本专利技术的有益效果是：本专利技术针对粒子滤波应用于突变故障诊断时诊断步数多和诊断结果噪声水平高的问题，在高斯粒子滤波框架下，提出了粒子集伪协方差的概念，通过考虑粒子与上一时刻估计值的偏差量，加剧突变故障时刻粒子在故障维度上的波动，保证采样得到的粒子更能反映真实突变情况，减少突变故障诊断步数；本专利技术采用高斯粒子滤波算法，通过高斯采样代替了粒子重采样，避免了粒子滤波重采样引起的样本贫化问题，从而提高滤波精度，保证了滤波稳定性；并且，伪协方差相比粒子集协方差，显著降低了算法的均方根误差，并且计算量基本不变，在保证计算速度的同时降低诊断结果噪声水平；综上，本专利技术适用于发动机突变故障诊断系统，能够有效提高突变故障的诊断精度和诊断速度。附图说明为了更清楚地说明本专利技术实施例中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图。图1是基于改进高斯粒子滤波的突变故障诊断原理图；图2是涡扇发动机模型气路工作截面示意图；图3(a)是高斯噪声下风扇效率突变8％标准粒子滤波器诊断效果；图3(b)是高斯噪声下风扇效率突变8％改进高斯粒子滤波器诊断效果；图4(a)是高斯噪声下风扇效率突变8％、压气机效率突变4％标准粒子滤波器诊断效果；图4(b)是高斯噪声下风扇效率突变8％、压气机效率突变4％改进粒子滤波器诊断效果；图5(a)是伽马噪声下风扇效率突变8％标准粒子滤波器诊断效果；图5(b)是伽马噪声下风扇效率突变8％改进高斯粒子滤波器诊断效果；图6(a)是伽马噪声下风本文档来自技高网...

【技术保护点】
1.基于改进高斯粒子滤波的涡扇发动机故障诊断方法，其特征在于，包括：S1、对发动机非线性数学模型注入气路部件突变故障；S2、设计基于伪协方差的改进高斯粒子滤波算法；S3、所述发动机非线性数学模型载入所述改进高斯粒子滤波算法，所述改进高斯粒子滤波算法对所述气路部件突变故障进行诊断,得到诊断结果。

【技术特征摘要】
1.基于改进高斯粒子滤波的涡扇发动机故障诊断方法，其特征在于，包括：S1、对发动机非线性数学模型注入气路部件突变故障；S2、设计基于伪协方差的改进高斯粒子滤波算法；S3、所述发动机非线性数学模型载入所述改进高斯粒子滤波算法，所述改进高斯粒子滤波算法对所述气路部件突变故障进行诊断,得到诊断结果。2.根据权利要求1所述的方法，其特征在于，所述发动机非线性数学模型表征了以下故障量:SEi为效率变化程度系数，SWi流量变化程度系数；所述效率变化程度系数和所述流量变化程度系数为健康参数。3.根据权利要求1所述的方法，其特征在于，所述气路部件突变故障包括：风扇部件突变故障，压气机部件突变故障，高压涡轮部件突变故障和低压涡轮部件突变故障。4.根据权利要求1所述的方法，其特征在于，S2具体包括：S21、k表示步数，在k＝0时，赋予初始值，再进入k＝1时刻，完成初始化，其中，为状态量估计的初始值，为粒子集伪协方差估计的初始值；S22、在k时刻，根据高斯采样从获得高斯采样粒子集其中为k-1时刻的等权重的高斯采样粒子，为k-1时刻的粒子加权均值，为k-1时刻的粒子集伪协方差；S23、在k时刻，从重要密度函数中依据重要性采样法抽样得到时刻更新粒子集其中，zk为k时刻的测量值；S24、选取先验密度函数作为重要性密度函数，即所述先验密度函数为计算得到各粒子的重要性权重即为所述各粒子的后验似然概率分布，S25、将所述各粒子权值归一化，为归一化值，S26、估计状态值及粒子集伪协方差S27、判断是否达到最大步数kmax，否则令k＝k+1，返回S22，循环执行S22-S27。5.根据权利要求1所述的方法，其特征在于，发动机故障诊...

【专利技术属性】
技术研发人员：黄金泉，卢俊杰，鲁峰，刘宸闻，王启航，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：江苏,32