专利名称：航空发动机转子系统工况振动信号特征信息提取方法
技术领域：
本发明涉及一种信号特征信息提取方法，尤其涉及一种航空发动机转子系统工况 振动信号特征信息提取方法。
背景技术：
航空发动机是飞机的心脏，是航空机械设备中故障率最高、调整最复杂、维护工作 量最大的系统，其工作状态的好坏直接影响到飞机的安全可靠与运行。一旦发动机在运行 中出现了故障，就可能导致飞机无法正常飞行，危及人身安全，造成重大的社会、经济损失。据统计，航空发动机的故障发生率约占整个飞机故障的30%。飞机因机械原因发 生的重大飞行事故中40%左右是由于发动机故障导致的。发动机由于自身特点不同于一般 机械，即使非常轻微的一些机械缺陷或损伤都会引起转子系统的振动，包括从几乎小得微 不足道的振动到大得足以导致发动机破坏的振动，由于航空发动机转子系统的复杂性，多 种故障表现为相同的振动，同一故障表现为不同的振动，往往导致漏报和误报，严重影响了 航空发动机转子系统故障诊断的可靠性。为了提高航空发动机转子系统故障诊断水平，需 要发展实用、高效的特征信息提取方法。目前，航空发动机振动故障诊断和信号特征提取的方法很多。例如，侯胜利、李应 红等在《基于主元核相似度免疫机制的故障诊断方法及应用》一文中基于主元核理论和免 疫系统机制，提出了基于主元核相似度免疫机制的故障诊断方法。文振华、左洪福在《基于 粗糙集_集成神经网络的航空发动机磨损故障诊断方法》一文中将粗糙集理论和神经网 络相结合并应用到航空发动机磨损故障诊断中，依据属性的重要性和决策表的相容性，采 用粗糙集理论对征兆信息进行属性约简，获取征兆的主要特征。徐启华、师军在《基于支持 向量机的航空发动机故障诊断》一文中提出一种基于支持向量机的航空发动机故障诊断方 法，应用该方法成功地对发动机气路部件的几种典型故障进行了正确诊断。王威、侯胜利 在《一种基于人工免疫理论的性能监控与故障诊断方法》一文中针对航空发动机故障样本 获取比较困难等问题，提出了一种基于人工免疫理论的航空发动机性能监控与故障诊断方 法。蔡开龙，谢寿生在《航空发动机的模糊故障诊断方法研究》一文中提出了一种基于T-S 模糊模型的故障诊断方法，将基于T-S模糊模型的故障诊断方法应用于航空发动机。这些方法都对航空发动机振动信号特征信息的提取做出了贡献，但诊断手段单 一，不能高效、准确、快捷的提取航空发动机转子系统的工况特征信息，全面反映航空发动 机转子系统的运行状态。

发明内容
本发明的目的在于提供一种航空发动机转子系统工况振动信号特征信息提取方 法，利用传感器测取振动位移信号，采用小波分析与聚类分析相结合的方法，从大量振动信 号中挖掘隐含的特征信息，实现工况特征的准确、高效提取。本发明是这样来实现的，其特征是方法为
4
1)特征向量获取通过电涡流位移传感器，按确定的时间间隔或采样频率测取确 定数量个航空发动机转子系统在某个采样时段的振动位移信号，对振动位移信号进行小波 变换，并获取该信号沿尺度(频率)轴上的能量模分布情况，然后将其按尺度顺序排列成向 量，作为特征向量；将确定数量个特征向量组成初始样本集Q，求取样本集Q特征向量作为 待诊断样本；2)聚类分析将Q按距离准则逐步聚类，类别由多到少(聚集法)，直到满足合适 的分类要求为止，通过比较待诊断样本与各个初始样本之间的距离，最终获得转子系统工 况振动信号特征信息。所述的特征向量获取的具体步骤如下(1)通过电涡流位移传感器，按确定的时间间隔或采样频率测取确定数量个航空 发动机转子系统在某个采样时段的振动位移信号；(2)选择小波基函数，并用该基函数对信号进行连续小波变换；(3)确定尺度参数的取值范围；(4)求信号的连续小波变换的系数的模，由式iVTfia.7) - {At). ΨβΛ( }) - -L {fit )ψ i 印求得连续小波变换的系数，其中 y/j ^ j
ftα ■\ a }
为炉j Ll j的共轭函数，再通过式乞》⑷印2得到其尺度-能量谱，然后用式w
转化成小波变换系数的模(5)按照尺度因子a由小到大将信号连续小波变换的能量模顺序排列，即得到与 转子运行状态相对应的特征向量，本向量为基于尺度_能量模的特征向量；(6)重复步骤(1) (4)，得到η个表征转子运行状态的特征向量，将η个特征向 量组成初始样本集，求取该样本集特征向量为待诊断样本。所述的聚类分析的具体步骤为(1)初始化令η个特征向量作为初始样本自成一类，即建立η个子集.^f5，
,计算各类之间的距离，可得到一个ηΧη维的距离矩阵D(a)，其中右上角标号 (0)表示聚类开始运算前的状态，设定迭代计算器为b = 0 (2)求距离矩阵D(b)中的最小元素，对角线元素除外，如果该最小元素为^>和-斤5 两类之间的距离，则将对)和4”合并为一类并由此建立新的分类(3)采用最短距离法和最长距离法计算合并后的新类别之间的距离，得距离矩阵假设A和B是两个聚类，则两类间的最短距离定义为《=版对，最长 距离为=e ^be B]；式中，&表示A类中的、样本和B类中的样本vt间的距离，巧Jt表示A类中的所有 样本与B类中所有样本之间的最小距离，O。表示A类中的所有样本与B类中所有样本之间 的最长距离，如果B类是由E类和F类合并而成的，则有D4j s 咖队 > ^A (4)令b = b+Ι，跳转到步骤(2)，重复计算及合并，直到D(b)满足条件，算法停止，所得到的分类结果即为聚类结果，将聚类结果合并得到距离矩阵；(5)比较待诊断样本与各个初始样本之间的距离，最终获得航空发动机转子系统 工况振动信号的特征信息。本发明的优点是将小波分析和聚类分析有机结合，既运用了小波分析良好的时频 特性和对突变信号和非平稳信号突出的处理能力，又利用了聚类分析能提高数据集中数据对 象相似度的特性。既克服了传统特征提取方法，如傅里叶变换，只能获取信号整体频谱，不能 获取局部频谱，只能对平稳信号分析，不能对突变信号和非平稳吸信号进行较良好分析的缺 陷，又克服了单一信号分析特征信息不明显的不足。大大提高了航空发动机运行状态识别的 准确性和效率，能应用于航空发电机转子系统的工况监测、智能诊断和智能监测。
具体实施例方式第一步，获取特征向量。1)选取采样频率为2560Hz，采样点数为512个，通过电涡流位移传感器分别测取 转子在不平衡、不对中、碰摩和油膜涡动四种故障运行状态下的振动位移信号各3组；(2)选Daubechies小波基函数对信号进行连续小波变换；(3)取尺度因子的范围为[2，30]之间的间隔为4的整数；(4)求信号的连续小波变换的系数的模；(5)按照尺度因子由小到大将信号连续小波变换的能量模顺序排列，得到与转子 运行状态相对应的特征向量。
(6)重复步骤(1) (4)，得到转子在不平衡、不对中、碰摩和油膜涡动四种 故障运行状态下的特征向量各三组，作为用于聚类分析的初始样本集Λ； (>. .\f,^K Af1 如表 L 用上述同样的方法提

取初始样本集特征向量作为待诊断的样本，得x^i=U. 0031,3.5982,6.4125,9. 8615, 13. 8290 17. 2595，19. 6139,2L 2796).表1基于尺度_能量模的转子特征向量(初始样本集)
尺瘴2610”ts222630不 平 mXfi1.03346. 49739,0^9U. 914 414. 616 717. 3725ia 9344Xf1,0328ζ. miOTJl11.799 714.606 S17. 36S9s a.Xfi.oecf3-STOi .0S6r11. 9Il βκ· 117. §1242α 5679不 φXf0.99273.590S6.39765. $82313. 317. 342 2IiL 698421. 3468Xfi1.00273.61086.414 9.870913, 85 617, 310ISL 676121. 3377χΤ0.膽3.5921δ. 4C629.862513.84517.2SS 31<l 646821.蘭逋 Xfo,$m3.41 5.δ2517.949210. 42 312.922 715. 144216. 9151Xfo.dm3.4004s. eiiti7.9602 α <54 3U. SiT i15. 163216. 9601Xfi0.$7333.3917S.C3097.9780 α 470 912. 9ω 315. 178616s 9829
VtOl Λ JO1.36064.3228δ. 637012.轮 8417.935 924 112 330.乡37. 9684m Pi1.55 4.93308.6723IZ 66117. m S24. 15θ 931. 024$37. 9961jtI 21.55694.95018.793312. 9CS418. Old 221城 31. C63S38, 0207第二步，聚类分析(1)初始化令12个特征向量自成一类，即建立12个子集《6)，…,砣5，即表 1.计算各类之间的距离，可得到一个12、12维的距离矩阵D(a)，如表2，其中右上角标号(0) 表示聚类开始运算前的状态，设定迭代计算器为b = 0。由于在这里各个信号之间的差异很 微弱，故采用欧氏距离能够较好地区分开。表2初始样本集构成的距离矩阵
7jkfΛΓ；**Jf;.,χ·,'·X"jr"Χ'Jir;*Χ*K'Χ“·.ΛΓ;，Jf,,*.00.0260.8104.45 i. 4054.35 4. 5Τ94.JH34.50525.64 925. 72 25,80 Iχγ·0.02600.8204.t7S4. 4264, 378 .5624.526i. 48825.66 525.7* 325,80 1X*'0.8100.82003.8255.2885.2525,21524.89 $24.97 525.04 $X'4.4574. 478S. 92200. 0580.1008.7188.67Τβ. 63921.94 022. Cl 422.08 4X '4. 4 5i. 426z.m0,05800. 0588.67i8.S328. 5$521.97 222. 04 622. H 6X'·4.3S74.378S. 8250. IW0.058O8.6208.5788.54122.02 622.10 022.17 0ΛΓ-"'i.57S4.5625.2^8.718β. SHS. 62000,0450,08330.16 230.2* ι30. 31 6Xt4.5434. 52&5.2528.S77S. 6328.5780.04500. 04230.12 530.20 530. 28 0X；*·4.5054,4885.2158.6S98. S蚝8.5410.沏30.042030.08 930.16 830.24 3-C2S.64 θ25.66 S2*. 89 621.94 02J.9T 222.02 630.16 250.12 530.08 900. 096α 189■C25.72 725.74 32t.於 S22.01 422. CH $22.10 030. 24 iJO. 20 530.16 ε0.09600. 093 25.80 25,80 125.04 922 08 422.11 622. 17 030.31 β30.28 030.24 30.1890.0930(2)求距离矩阵D(b)中的最小元素(对角线元素除外)，如果该最小元素为於丨和 两类之间的距离，则将广和jrf合并为一类并由此建立新的分类^"，....(3)采用最短距离法和最长距离法来确定合并后的新类别之间的距离，得距离矩 阵产.(4)令b = b+Ι，跳转到步骤(2)，重复计算及合并，直到D(b)满足条件，算法停止， 所得到的分类结果即为聚类结果，经过8次聚类合并得到距离矩阵，如表3.表3经计算合并的距离矩阵
^tilΑ'；；ν *ΛΙ111 !0-0,8203.8S&-4. +784. 48 Κ5. 23824- S96-2S. mΛ ** S. 82&-<. 4 80-0. ΟΟ8. S41-8. 71821. 940-22.170Λ t 4. 488-5.2888.541-8. Πδ0-0.083m. 0W-30, 316,.rmm,24. m-m. m21.940-22, JTO30.089-30. 316该算法通过对数据集的一次扫描就可以生产质量比较好的聚类，并且可以通过
追加扫描进一步提高聚类质量，该算法的时间复杂度是0，N, (N是数据集中的数据对象数 目)。通过上述聚类算法合并后的距离矩阵，我们可以看出 该算法能够很好地把故障初始样本集分类，再将待诊断的样本 广=α 0031,3. 5982,6. 4125，9. 3615，13. 8290, 17. 259S, 19. 6139,21. 2796)^P入到样本集，然后求取待诊断样本和各个初始样本之间的欧式距离，得表4，
表4待诊断样本与各初始样本的距离
权利要求
一种航空发动机转子系统工况振动信号特征信息提取方法，其特征是方法为1)特征向量获取通过电涡流位移传感器，按确定的时间间隔或采样频率测取确定数量个航空发动机转子系统在某个采样时段的振动位移信号，对振动位移信号进行小波变换，并获取该信号沿尺度(频率)轴上的能量模分布情况，然后将其按尺度顺序排列成向量，作为特征向量；将确定数量个特征向量组成初始样本集Q，求取样本集Q特征向量作为待诊断样本；2)聚类分析将Q按距离准则逐步聚类，类别由多到少(聚集法)，直到满足合适的分类要求为止，通过比较待诊断样本与各个初始样本之间的距离，最终获得转子系统工况振动信号特征信息。
2.根据权利要求1所述的航空发动机转子系统工况振动信号特征信息提取方法，其特 征是所述的特征向量获取的具体步骤如下(1)通过电涡流位移传感器，按确定的时间间隔或采样频率测取确定数量个航空发动 机转子系统在某个采样时段的振动位移信号；(2)选择小波基函数，并用该基函数对信号进行连续小波变换；(3)确定尺度参数的取值范围；(4)求信号的连续小波变换的系数的模，由式巧阶·卜(/^. 《^ + .〖邝冲‘丨^^卜求得连续小波变换的系数，其中…^丨为Ve Λο1 α ‘ν!-11的共轭函数，再通过式民=g㈣力收作1得到其尺度-能量谱，然后用式& Μ。,转化成小波变换系数的模(5)按照尺度因子a由小到大将信号连续小波变换的能量模顺序排列，即得到与转子 运行状态相对应的特征向量，本向量为基于尺度_能量模的特征向量；(6)重复步骤(1) (4)，得到η个表征转子运行状态的特征向量，将η个特征向量组 成初始样本集，求取该样本集特征向量为待诊断样本。
3.根据权利要求1所述的航空发动机转子系统工况振动信号特征信息提取方法，其特 征是所述的聚类分析的具体步骤为(1)初始化令η个特征向量作为初始样本自成一类，即建立η个子集Aw, .…，Λ'Γ，计算各类之间的距离，可得到一个ηΧη维的距离矩阵D(a)，其中右上角标号(O)表示聚类开始运算前的状态，设定迭代计算器为b = 0 ；(2)求距离矩阵D(b)中的最小元素，对角线元素除外，如果该最小元素为A^和-两类 之间的距离，则将对合并为一类并由此建立新的分类-^0,....(3)采用最短距离法和最长距离法计算合并后的新类别之间的距离，得距离矩阵 假设A和B是两个聚类，则两类间的最短距离定义为成=^ B\ ,最长距离为= Ι1 Χ e式中，心表示A类中的\样本和B类中的样本&之间的距离，表示A类中的所有样本 与B类中所有样本之间的最小距离，&表示A类中的所有样本与B类中所有样本之间的最 长距离，如果B类是由E类和F类合并而成的，则有(4)令b= b+1，跳转到步骤(2)，重复计算及合并，直到满足条件，算法停止，所得到 的分类结果即为聚类结果，将聚类结果合并得到距离矩阵；(5)比较待诊断样本与各个初始样本之间的距离，最终获得航空发动机转子系统工况 振动信号的特征信息
全文摘要
一种航空发动机转子系统工况振动信号特征信息提取方法，其特征是方法为1)通过电涡流位移传感器测取振动位移信号；2)对信号进行小波变换，获取该信号沿尺度轴上的能量模分布情况，将其按尺度顺序排列成向量，作为特征向量，将多个特征向量组成初始样本集Q，求取样本集Q特征向量作为待诊断样本，3)将Q按距离准则逐步聚类得到聚类结果，比较待诊断样本与初始样本之间的距离，获得转子系统工况振动信号特征信息。本发明的优点是克服了单一方法信号特征信息提取方式复杂、特征信息不明显等不足，提高了航空发动机转子系统运行状态识别的准确性和效率,能应用于航空发电机转子系统的工况监测、智能诊断和监测。
文档编号G01H11/06GK101968379SQ20101029838
公开日2011年2月9日 申请日期2010年9月30日 优先权日2010年9月30日
发明者刘晓波, 孙康, 沈亮霓, 王志华 申请人:南昌航空大学