发明名称 基于退化模型匹配的涡轮发动机剩余寿命的预测方法
摘要 基于退化模型匹配的涡轮发动机剩余寿命的预测方法,涉及涡轮发动机剩余寿命预测方法,它为了解决现有涡轮发动机剩余寿命的预测采用通用RUL预测模型的预测效果无法达到预测要求的问题,它包括具体步骤如下:步骤一、数据预处理:从采集到的数据提取运行状态变量;从传感器采集到特征向量;由运行状态变量与特征向量融合得到健康因子;步骤二、建立退化模型库:利用健康因子建立退化模型;多个退化模型组成退化模型库;步骤三、相似性评估:将退化轨迹与模型库中的模型匹配,每个模型给出一个RUL估计;步骤四、RUL融合:根据测试涡轮发动机与模型匹配程度,采用相似度加权,融合得到最终的剩余寿命预测值。适用于涡轮发动机剩余寿命预测。
申请公布号 CN102789545B 申请公布日期 2015.08.19
申请号 CN201210240981.8 申请日期 2012.07.12
申请人 哈尔滨工业大学 发明人 刘大同;李君宝;徐勇;罗悦;庞景月;王红;彭宇
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 张宏威
主权项 基于退化模型匹配的涡轮发动机剩余寿命的预测方法,它包括具体步骤如下:步骤一、数据预处理:从采集到的原始数据中提取出决定涡轮发动机运行状态变量:c<sub>1</sub>为高度、c<sub>2</sub>为马赫数和c<sub>3</sub>为油门角度;所述运行状态变量的数值组成集合c=(c<sub>1</sub>,c<sub>2</sub>,c<sub>3</sub>)<sup>T</sup>;c<sub>i</sub>表示c在时刻t<sub>i</sub>对应的高度、马赫数和油门角度的值,c<sub>i</sub>=(c<sub>1i</sub>,c<sub>2i</sub>,c<sub>3i</sub>)<sup>T</sup>,其中i为正整数;从21个位于涡轮发动机不同位置的传感器采集到能表征涡轮发动机退化过程的数值,所述数值组成21维特征向量:x=(x<sub>1</sub>,x<sub>2</sub>,...,x<sub>21</sub>)<sup>T</sup>,x<sub>i</sub>表示x在时刻t<sub>i</sub>对应21个传感器采集到的发动机不同点的温度、压力和速度值;在时刻t<sub>i</sub>由所述运行状态变量c<sub>i</sub>与特征向量x<sub>i</sub>,得到测量周期的三元组(c<sub>i</sub>,x<sub>i</sub>,t<sub>i</sub>);由c<sub>i</sub>和x<sub>i</sub>融合得到健康因子z<sub>i</sub>;步骤二、建立退化模型库:利用融合而成能表征系统退化趋势的健康因子z<sub>i</sub>为相应的训练实体建立退化模型;所述建立的多个退化模型组成退化模型库;步骤三、相似性评估:将某测试涡轮发动机的健康因子序列与模型库中所有的模型进行相似匹配,每个模型都给出一个对该测试涡轮发动机的RUL估计;步骤四、RUL融合:根据测试涡轮发动机与模型库中模型匹配程度的高低,采用相似度加权的方法融合得到最终的剩余使用寿命的RUL预测值,其特征在于,步骤一数据预处理具体过程为:步骤一一、对决定涡轮发动机运行状态变量的集合c=(c<sub>1</sub>,c<sub>2</sub>,c<sub>3</sub>)<sup>T</sup>采用K‑means算法聚类得到p类运行状态Ω={O<sub>1</sub>,O<sub>2</sub>,...,O<sub>p</sub>},p为正整数;步骤一二、对提取的每个训练实体中传感器的读数按p类运行状态Ω进行划分,共划分为p组;步骤一三、检测所有训练实体中每个传感器的读数在每个运行状态下随时间变化的情况,选择出在所有运行状态下读数均具有一致性单调趋势的传感器集合Α={X<sub>1</sub>,X<sub>2</sub>,...,X<sub>m</sub>},其中m为正整数;步骤一四、某组训练实体对选出的传感器集合Α的读数按有限的运行状态Ω进行划分,并按照<img file="FDA0000667183430000011.GIF" wi="714" he="101" />进行线性回归得到p组健康因子时间序列;步骤一五、p组健康因子时间序列按状态划分前的时间排序后,还原成完整的一组健康因子时间序列X<sub>0</sub>={X<sub>0</sub>(k)|k=1,2,……,n},其中n为正整数;步骤一六、确定反映系统退化行为特征的健康因子时间序列X<sub>0</sub>为参考序列,影响系统行为的传感器时间序列X<sub>i</sub>={X<sub>i</sub>(k)|}为比较序列,i=1,2,……,m;步骤一七、采用初值化对参考序列和比较序列进行无量纲化处理,用原始数列的所有数据除以序列的第一个值,即<img file="FDA0000667183430000021.GIF" wi="307" he="143" />步骤一八、计算经过无量纲化处理后的参考序列x<sub>0</sub>和比较序列x<sub>i</sub>之间的点关联系数为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>r</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munder><mi>min</mi><mi>i</mi></munder><munder><mi>min</mi><mi>k</mi></munder><mo>|</mo><msub><mi>x</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>x</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo><mo>+</mo><mi>&zeta;</mi><munder><mi>max</mi><mi>i</mi></munder><munder><mi>max</mi><mi>k</mi></munder><mo>|</mo><msub><mi>x</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>x</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo></mrow><mrow><mo>|</mo><msub><mi>x</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>x</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo><mo>+</mo><mi>&zeta;</mi><munder><mi>max</mi><mi>i</mi></munder><munder><mi>max</mi><mi>k</mi></munder><mo>|</mo><msub><mi>x</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>x</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000667183430000022.GIF" wi="1357" he="182" /></maths>其中ζ∈(0,∞)称为分辨系数,取ζ=0.5;步骤一九、计算参考序列和比较序列间的关联度,用参考序列和比较序列在各个时刻的关联系数之平均值表示,即:<img file="FDA0000667183430000023.GIF" wi="618" he="141" />步骤一十、计算点关联序数序列的稳定度<img file="FDA0000667183430000024.GIF" wi="778" he="157" />步骤一十一、计算改进后的灰色关联度<img file="FDA0000667183430000025.GIF" wi="451" he="143" />步骤一十二、对计算得到的关联度从大到小进行排序;步骤一十三、对所有训练实体重复步骤一四至步骤一十二,得到所有训练涡轮发动机中各个传感器读数与相应健康因子关联度排序的统计信息,选出关联度高的传感器子集B={X<sub>1</sub>,X<sub>2</sub>,...,X<sub>t</sub>}。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号