一种航空发动机传感器故障诊断方法

一种航空发动机传感器故障诊断方法黄强;王健;张桂刚【摘要】Inordertodiagnoseaero—enginesensors'faultsaccuratelyandquickly,theenginesensors'hardfaultdi—agnosisofconstantdeviation,constantgain,stuckareachievedwiththeresidualtestfaultdiagnosismodulebasedonKalmanfilter.ThemoduleisdesignedinSimulink.Thereisabigdelayindetectionofsoftfaultwithtraditionalre一sidualtestmethodbasedonKalmanfilter,sotheimprovedmethodisadoptedtodetectfaultandcomparedwiththetraditionalmethod.SimulationresultshowthatthetraditionalresidualtestmethodbasedonKalmanfiltercaneffec—tivelydiagnoseseveraltypicalaero—enginesensors'hardfaults,butitisinsensitivetosoftfaultdetection.Theim一provedmethodcanacceleratethespeedofdetectionandachieveagoodperformance.%为了对航空发动机传感器故障作出准确快速的诊断，在Simulink中设计了基于卡尔曼滤波器的残差检验的故障诊断模块，实现发动机传感器恒偏差、恒增益、卡死故障等硬故障的诊断。针对传统卡尔曼滤波器残差检验法检测缓变故障时存在较大延时的缺点，采用改进卡尔曼滤波器的残差检验法来检测故障，并与传统卡尔曼滤波器法进行比较。仿真结果表明，传统卡尔曼滤波器的残差检验法能有效地诊断航空发动机传感器的几种典型硬故障，但对于缓变故障却不敏感。改进卡尔曼滤波器的残差检验法能在一定程度上加快检测速度，取得较好的效果。【期刊名称】《传感技术学报》【年(卷)，期】2014(000)010【总页数】6页(P1315-1320)【关键词】测试计量技术及仪器;故障诊断;卡尔曼滤波器;残差检验;航空发动机传感【作者】黄强;王健;张桂刚【作者单位】中国科学院自动化研究所，北京100190;北京科技大学自动化系，北京100083;中国科学院自动化研究所，北京100190;中国科学院自动化研究所，北京100190【正文语种】中文【中图分类】TP212.6随着全权限数字电子控制FADEC(FullAuthorityDigitalEngineControl)系统的快速发展，人们对航空发动机的可靠性提出了更高的要求，因此航空发动机上部署了各种各样的传感器以监测发动机的工作状况。大量发动机传感器都在强震动,高温，高压等恶劣的环境中工作，因此，业界越来越重视传感器故障的及时诊断与排除。故障诊断的方法有很多种，但大致可分为两类:基于数据的故障诊断方法和基于模型的故障诊断方法。在基于数据的故障诊断方法中，故障诊断系统可以直接从各种各样的传感器得到表征系统运行状态的测量数据,并且对数据进行处理，从提取的特征数据中获得系统当前的健康状态，从而进行故障诊断，得出结果。常见的基于数据的故障诊断方法有模糊神经网络法［1］、径向基神经网络法［2］、相关向量机法［3］、支持向量机法［4］等等。在基于模型的故障诊断方法中，故障诊断系统并不像基于数据的故障诊断方法那样从传感器获取表征运行状态的数据，因为这时系统的模型是已知的，并且系统模型已经确定了系统的状态、参数和输入/输出等之间的关系，从而可以通过模型分析、变化监测、滤波器、观测器以及分类与推理的方法来实现系统的故障诊断功能。常见的基于模型的故障诊断方法有:观测器法［5］、等价空间法［6］、故障树分析法［7］、因果图分析法［8］、卡尔曼滤波器法［9］等等。采用卡尔曼滤波方法来实现航空发动机传感器的故障诊断是目前比较成熟的做法，因此，本文采用卡尔曼滤波器法来诊断传感器故障。航空发动机传感器有恒偏差故障，恒增益故障，卡死故障，漂移故障，短路故障，冲激故障，周期性干扰故障等硬故障，本文选用恒偏差，恒增益,卡死这3种常见硬故障［10］。传感器在运行过程中，随着运行时间的增长，会出现性能退化，这就是传感器的软故障。恒增益故障当系统运行到某一时刻t1发生恒增益故障，在之后的时间里yout=a・yin,得恒偏差故障当系统运行到某一时刻t1发生恒偏差故障，在之后的时间里yout=yin+&得卡死故障当系统的输出为P时发生卡死故障，在之后的时间里yi=B，得⑷软故障(性能下降)当系统运行到某一时刻t1发生软故障，得卡尔曼滤波器进行故障诊断的基本思想是:先用卡尔曼滤波器对系统的状态进行估计,接着从估计的状态得到系统的估计输出值，再把系统输出真实值与系统输出估计值进行比较，得到残差。当系统正常工作时，残差的均值近似为0,这是因为卡尔曼滤波值是对x(k)的线性最小方差估计,但当系统发现故障时，残差会迅速变大，因此，可以通过分析残差从而知道系统是否发生了故障。2.1基于卡尔曼滤波器的残差生产方法假设系统如下：式中3(k)是估计噪声,v(k)是测量噪声。根据式(5),可以得到卡尔曼滤波器法的残差生成原理图，如图1所示：所以计算残差的过程为：计算一步最佳估计值初始值计算预报误差即为残差r(k);式中是测量估计值,y(k)是测量值。计算预报误差方差阵P(k|k-1);式中P(0)=Var[x(0)],为对角矩阵。计算滤波增益矩阵K(k);⑸计算滤波方差矩阵P(k);(6)计算k时刻x(k)的最优估计值重复上面的过程就可以进行实时状态估计，并且可以得出下面用于故障诊断的残差r(k)。2.2基于卡尔曼滤波器的故障诊断策略当系统出于正常状态时，卡尔曼滤波器产生的残差r(k)近似为高斯白噪声，平均值近似为0,但在故障状态时,r(k)会明显增大，这样就可以判断故障。本文采取WSSR(WeightedSum-SquaredResidual)检测方法[10],该方法如下：式中,N是人为选择的数据窗口长度。N取得较大时能降低故障的误报率，但同时也减慢了故障的检测速度,因此需要均衡考虑故障的准确率和快速性来选择N。因此可以把故障检测器表示为式中M是人为选择的阈值。因此可以取Lk+1大于某一门限M的概率为a,即:p{LK+1>M}=a，其中a为允许的虚警概率。此时检测过程为：Lk+1>M系统异常Lk+1<M系统正常门限M的选择要兼顾虚警概率a和遗漏概率p两个因素。其中，虚警概率a为实际系统故障而判为故障的概率;遗漏概率P为实际系统故障而判为无故障的概率[11]。对航空发动机传感器进行故障诊断之前，必须要建立发动机传感器模型(状态方程)，得到仿真数据，再在某一时间点上加入故障，以便对其进行故障的检测。因此，首先要做的是用Simulink建立发动机传感器模型，这个模型是为了产生仿真数据，再在某一时间点人为注入故障。然后对待检测的数据进行卡尔曼滤波，与原输出相比较得到残差，最后进行残差分析，确定故障。3.1传统卡尔曼滤波器残差检验法故障诊断3.1.1仿真数据产生模块的实现选取一个典型传感器的状态方程用来产生数据，分别在输入端和输出端加上两个独立的高斯白噪声。传感器的状态方程如式(5)所示,其中：;H(k+1)=[10]。3.1.2卡尔曼滤波器及残差检验模块的实现卡尔曼滤波器的实现卡尔曼滤波器算法根据式(6)~式(11),在Simulink中用S函数编写。残差检测模块的实现根据前文叙述的残差检测(WSSR)方法。本文取N=3,因此得到残差检测模块如图2所示：图2中UnitDelay是延时模块，产生0.001s的延时,MathFunction是矩阵转置模块,Product是乘法模块,Add是加法模块，将相乘后的结果相加。图3中RelationalOperator1是判断模块，判断残差加权平方和是否超过了阈值,Constant1中的0.5是阈值,TriggeredSubsystem1模块的作用是当残差超过阈值时就输出该时间点的残差加权平方和及故障时间点,Display2显示残差加权平方和,Display3显示故障时间点。3.1.3故障诊断连接各个模块即可得到整个故障诊断系统模型图，如图4所示。故障诊断计算流程图如图5所示。在Simulink中分别设置恒增益故障，恒偏差故障，卡死故障,软故障，得到结果如下(残差加权平方和超过0.5时即可认为发生了故障)。恒增益故障在5s时加入故障,Simulink显示5.001s时检测到故障，因此满足故障检测准确快速的要求。恒偏差故障在5s加入故障,Simulink显示5.001s时检测到故障，因此满足故障检测准确快速的要求。卡死故障Simulink显示在5.025s时检测到故障满足要求。⑷软故障(性能下降)在5s加入故障,Simulink显示在5.224s时检测到故障，但检测到故障却用了224个采样周期(本文设置采样周期为0.001s),不满足故障检测的快速性要求。从上面的实验结果可以看出，卡尔曼滤波器法对于诊断传感器恒增益，恒偏差，卡死等硬故障具有准确快速的特点，这满足了故障诊断的性能要求。可是面对性能退化这一软故障时，卡尔曼滤波器虽然能正确地判断出故障，但是延时了224个采样周期，这延时相对来说比较大，不满足系统故障诊断的快速性，因此，本文接着采用了一种改进的算法来诊断软故障。3.2改进卡尔曼滤波器残差检验法故障诊断3.2.1改进算法介绍若在k时刻之前的某一个时刻(k-m)开始发生传感器的缓变故障，但该故障在k时刻之前一直未能被发现，那么系统状态估计值中含有故障信息，则量测预报值为中该预报值也包含了故障信息,因此与yi(k)的差值不够大，以至于故障不能被识别［12］。若系统预测状态的先验值取(k-n)时刻的状态估计值当n的值比较恰当时，则可以认为不包含系统的故障信息或者被污染程度很小，该状态可以被称为伪正常状态。这么做后就可以用状态递推器来检测故障。将k时刻的量测预报值改为：式中中(k/k-1)。则k时刻的残差为：当无故障发生时，可以通过证明得rs(k)是零均值白噪声［11］。只要n取的恰当，故障基本不会影响到量测预报值，因此判断传感器是否发生了软故障就可以通过检查残差来实现。对于n的选取需要注意两方面会随着n的增大，受故障的影响程度也减小，也许会有可能根本不受故障影响；(2)随着n的增大的预报方差也会增大，会降低检查效率，增大计算量。由此可以看出,n的选择需要权衡上面两方面的影响。3.2.2故障诊断在Simulink设置软故障,得到结果如下(本文取n=200,5s时加入故障)。在5s时加入故障,Simulink显示在5.15s时检测的故障，虽然延时了150个采样周期，但是比未改进前检测速度加快了。基于卡尔曼滤波器的故障检测方法对软故障的检测不敏感，导致故障积累到一定程度也不能及时检测出。基于改进的卡尔曼滤波器故障检测法检测速度明显比未改进前快，在一定程度上保证了故障检测的快速性。因此，对于检测传感器软故障,基于改进的卡尔曼滤波器故障检测法比未改进前更为有效。本文首先使用了传统卡尔曼滤波器的残差检验法来诊断航空发动机传感器的几种典型故障，对硬故障的诊断收到了很好的效果。针对传统卡尔曼滤波器的残差检验法诊断传感器软故障存在较大延时的缺点，本文使用了改进卡尔曼滤波器(基于移动伪正常状态)的卡尔曼滤波器法来诊断故障检测速度明显加快。因此，本文所采取的算法有效、正确、可靠，这是对基于模型的PHM算法的一次有益的探索。黄强(1991-),男，硕士研究生，主要研究方向为智能控制与计算智能；王健(1969-),男，中国科学院自动化研究所副研究员，硕士研究生导师,主要研究方向为智能计算与智能控制、飞机综合健康管理,.cn;张桂刚(1978-),男，博士，中国科学院自动化研究所副教授，硕士研究生导师，主要研究方向为大数据、大飞机综合健康管理,.cn 。【相关文献】刘宜平，沈毅，刘志言.一种基于模糊神经网络的故障分类器及其在多传感器故障诊断中的应用[J].传感技术学报,2000,13(1):38-43.房方，魏乐.传感器故障的神经网络信息融合诊断方法[J].传感技术学报,2000,13(4):272-276.TippingME.SparseBayesianLearningandRelevanceVectorMachine[J].JournalofMachineLearingResearch,2001,24(1):211-244.曾小军，黄宜坚.基于AR模型和支持向量机的故障诊断法[J].机械科学与技术,2010,29⑺:972-976.PattonR,ChenJ.Observer-BasedFaultDetectionandIsolation:RobustnessandApplications[J].ControlEngineeringPractice,1997,5(5):671-682.FrankP.On-LineFaultDetectioninUncertainNonlinearSystemsUsingDiagnosticObserver: