主权项 |
基于退化模型匹配的涡轮发动机剩余寿命的预测方法,它包括具体步骤如下:步骤一、数据预处理:从采集到的原始数据中提取出决定涡轮发动机运行状态变量: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>ζ</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>ζ</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>}。 |