基于损伤传播模型和数据分析的航空发动机剩余寿命预测的制作方法
【技术领域】
[0001] 本发明属于航空发动机预测与健康管理技术领域，特别涉及一种基于损伤传播模 型和数据分析的航空发动机剩余寿命预测方法。
【背景技术】
[0002] 航空发动机是个典型的复杂机械系统，由大量的热力旋转部件构成，长期处在高 温、高压、高转速等恶劣的工作环境下，且经常面临未知的外部环境;随着航空工业的迅猛 发展，航空设备不断的更新换代，飞机的综合性、复杂性和智能化程度不断提高，对飞机的 机动、安全、经济、环保等性能提出了更高的要求，这对航空发动机的要求变得日益苛刻；因 此，航空发动机在整个飞机系统中是个敏感、易损的故障多发部件。对航空发动机进行故障 诊断、预测和健康管理(PHM)具有十分重大的意义。
[0003] 航空发动机预测是指对部件、系统剩余可用的或安全的寿命的估计，即寿命预测， 是航空发动机健康管理中技术含量最高、难度最大的综合性技术。有研究指出，大部分装备 部件的故障率不服从浴盆曲线分布，由此提出了基于状态的维修，即以往对部件的监控是 及时发现并判断出故障，不具备预测的能力。通过预测，可以提高飞行的安全性，任务的完 成率;根据预测结果有效安排维护/检修工作，去除失效部件，根据其余部件的剩余使用寿 命做出合理决策，最大化部件使用价值的，可以大大降低维修成本。针对整台航空发动机的 寿命周期预测，可提前对航空发动机进行检修，防止航空发动机在空中出现故障，避免重大 安全事故。对民用航空发动机而言，通过预测航空发动机换发时间，可有效地进行整个航空 发动机机队的使用及大修计划管理。对军用航空发动机而言，预测航空发动机在未来某特 定时间和特定环境下完成任务的能力，也是非常重要的。
[0004] 飞机发动机从某个初始退化点开始，流量和效率会发生指数变化，对其他参数产 生影响，故障传播继续直到达到某个故障标准。损伤传播建模跟踪和预测飞机发动机涡轮 机械的损伤进程，损伤传播建模在不同应用领域使用的模型包括阿伦尼斯模型，Coffin-Manson机械裂纹增长模型，和艾林模型（超过三种应力或当前两项模型不符合要求时采 用）。这些模型具有一个共同点，即故障演变是种指数变化的过程。考虑宏观的退化特征，假 设磨损的传播遵循w=Ae B(t)。
[0005] 对于寿命预测的数据分析方法有递归神经网络，SVM方法，隐马尔科夫模型等。本 发明涉及的数据分析方法包括已公开的线性回归，非线性最小二乘法，相似性分析和K近邻 法等基础算法。
【发明内容】
[0006] 本发明的目的针对目前已有技术存在的问题，提供一种基于损伤传播模型和数据 分析的航空发动机剩余寿命预测方法，本发明方法能够给航空发动机的预测和健康管理提 供有效的技术支持，给航空发动机的后续维修决策提供有效的参考。
[0007] 为实现上述目的，本发明提供了一种基于损伤传播模型和数据分析的航空发动机 剩余寿命预测方法，其特征在于，该方法包括如下步骤：
[0008] 1)采集某型航空发动机的全寿命运行周期的传感器测量数据，从中选择构建损伤 传播模型的传感器参数；
[0009] 2)将选中的构建损伤模型健康指标h的传感器参数融合成描述航空发动机的健康 状况的健康指标，形成健康指标序列Ih(k)}，其中h(k)代表航空发动机在第k次运行时的健 康指标，ke [1，n]，n代表该航空发动机总共的运行次数；
[0010] 3)将每组健康指标序列{h(k)}拟合一个指数型损伤传播模型Mi，Mi代表第i组训练 数据对应的损伤传播模型并将各组健康指标序列{h(k)}拟合的指数型损伤传播模型组成 损伤传播模型库{Mi};
[0011] 4)对待预测的多台该型号航空发动机的历史运行的传感器测量数据，按照步骤1) 选取用于构建损伤传播模型的传感器参数，将这些传感器参数融合成健康指标，形成多台 待预测航空发动机的健康指标序列lh(k)};
[0012] 5)将每台待预测的航空发动机的健康指标序列{h(k)}与步骤3)中损伤传播模型 库{MJ中的模型进行相似性匹配；
[0013] 6)根据每台待预测的航空发动机的匹配结果预测该台航空发动机的剩余寿命，将 这些剩余寿命采用k近邻法进行融合获取最终的剩余寿命。
[0014] 本发明的特点及有益效果：
[0015] 本发明将损伤传播建模和数据分析方法应用到航空发动机领域，针对航空发动机 从健康运行到失效的全寿命周期的数据，采用单调变化的原则对传感器数据进行了选取， 将选中的传感器数据融合成航空发动机健康指标;根据损伤传播建模方法使用最小二乘法 进行了发动机失效模型库的搭建;将在役航空发动机的历史数据与模型库进行相似性匹 配，使用了标准化欧式距离进行相似度比较，求得每组测试数据在模型库下的剩余寿命库； 最后采用了 k近邻法对剩余寿命库进行了最终剩余寿命的求取。本发明能够给航空发动机 的预测和健康管理提供有效的技术支持，给航空发动机的后续维修决策提供有效的参考。
【附图说明】
[0016] 图1为本发明一种基于损伤传播模型和数据分析的航空发动机剩余寿命预测方法 的实施例流程框图；
【具体实施方式】
[0017] 下面结合附图及实施例对本发明进行详细说明。
[0018] 本发明的一种基于损伤传播模型和数据分析的航空发动机剩余寿命预测方法的 实施例流程，如图1所示，该方法包括如下步骤：
[0019] 1)采集某型航空发动机的全寿命运行周期的传感器测量数据，从中选择构建损伤 传播模
[0020] 型的传感器参数；
[0021] 具体包括以下步骤：
[0022] 11)本实施例中，利用美国国家航空航天局(NASA)卓越预测中心的涡扇发动机(航 空发动机的一种)退化仿真数据集作为分析的数据对象，该数据集来自NASA的民用模块化 航空推进系统仿真软件(CMAPSS)的仿真结果，CMAPSS通过改变涡扇发动机模型部件的流量 和效率（以指数率变化的形式)来仿真退化的故障，当故障达到磨损阈值(表示一个超过了 部件/子系统不能使用的操作限制)时认为该涡扇发动机模型失效。CMAPSS中涡扇发动机模 型的输出是传感器测量值的时间序列;本实施例选用CMAPSS数据集中的FD002数据(该数据 包括一共有6种运行状态，21个传感器测量参数）；该数据集包括训练数据，测试数据和测试 数据的剩余寿命；
[0023] 12)根据某些应用领域已有成熟的损伤传播模型（如应用于各种故障机械机理的 阿伦尼斯模型，应用于机械故障、材料疲劳的Coffin-Manson模型，应用于温度和应力的艾 琳模型），其故障演变具有指数变化的过程的特点；假设磨损的传播遵循w = AeB(t)(其中，w 是磨损程度，A是比例因子，B(t)是随时间变化的指数），磨损发展到某个阈值thw时，系统将 失效;设用以描述航空发动机的健康状况的健康指标h，在航空发动机初始运行时，该健康 指标h值为1，在航空发动机的磨损到达阈值1匕时，该健康指标h值退化为0;由于健康指标h 的值是随着时间变化的，为此设以时间为变量的健康方程为：
[0024] h(t) = l-AeB(t)/thw (1)；
[0025] 13)将Π )002训练数据的运行状态划分成多种运行状态，并将多个传感器测量值转 化为健康指标h，再从中选择出作为构建损伤模型健康指标h的传感器参数;本实施例划分 为6种运行状态；由于健康方程中健康指标h是个单变量，因此需要将21个传感器测量值转 化为健康指标h;本实施例选择在6种不同的运行状态下都符合该要求的传感器参数测量 值，最终选择了 21个传感器测量值中的5个测量值作为构建损伤模型健康指标h的传感器参 数(损伤传播模型中的健康指标h是指数型单调变化的，融合参数时只选择单调变化且受噪 声影响较小(参数变化量对参数本身而言较大）的参数），这5个传感器测量值对应的航空发 动机变量分别为:LPC出口总温，HPC出口总温，LPT出口总温，HPC出口静压，旁通比；
[0026] 2)将选中的构建损伤模型健康指标h的传感器参数融合成描述航空发动机的健康 状况的健康指标，形成健康指标序列lh(k)}，其中h(k)代表航空发动机在第k次运行时的健 康指标，ke [1，n]，n代表该航空发动机总共的运行次数；
[0027] 本实施例使用线性回归模型对选中的5个构建损伤模型的传感器参数进行融合， 转换成健康指标h:
[0028]
[0029] 其中，x = (XI，X2, . . .，xn)是N(N = 5)维特征向量，h为健康指标；（α，β) = (α ,βχ, 02，···，βΝ)是Ν+1维模型参数；
[0030] 从训练数据集中选取样本集Q={(x，h)}，用以求取线性回归模型：
[0031]
[0032] 式中，下标i代表训练数据集中第i组数据序列，Cf代表第i组数据的第r次运行； Ch〇代表健康指标h取0的次数;Chl代表健康指标h取1的次数；
[0033] 本实施例中ro〇02训练数据包含260组训练数据(每组训练数据代表一台航空发动 机的全寿命周期的运行过程），将260组训练数据按运行状态划分为6组状态数据(本实施例 划分了 6种运行状态，每种运行状态对应1组状态数据）；对每种运行状态对应的1组状态数 据使用公式(3)选取样本集Ω，将样本集Ω代入公式(2)中，使用非线性最小二乘法求取该 运行状态下公式(2)中的模型参数(α，β)(每种运行状态应对一组模型参数(α，β)，本实施例 划分了 6种运行状态，因