基于包络谱分析的航空发动机主轴轴承故障诊断
2021-08-09
来源:星星旅游
第30卷第4期 2 0 l 3年8月 沈阳航空航天大学学报 Journal of Shenyang Aerospace University VO1.30 No.4 Aug.2 0 1 3 文章编号:2095—1248(2013)04—0018—05 基于包络谱分析的航空发动机主轴轴承故障诊断 梁先芽,沙云东,栾孝驰,赵奉同,张强 (沈阳航空航天大学辽宁省航空推进系统先进测试技术重点实验室,沈阳110136) 摘要:针对某型航空发动机地面检测条件下采集到的机匣振动信号,利用小波变换和包络分析相 结合的方法提取主轴轴承故障信息并进行故障诊断。首先对所采集到的航空发动机振动信号进 行小波降噪,再利用小波变换提取降噪后的滚动轴承故障特征信号,然后对该故障特征频段进行 包络谱分析,以获得特征峰值频率。将该方法应用到试验数据中,分析发现该方法能够有效地诊 断出航空发动机主轴轴承故障及具体的故障位置。这为航空发动机主轴轴承故障诊断提供了重 要的判断依据,具有广阔的应用前景。 关键词:小波变换;包络谱;航空发动机;主轴轴承;故障诊断 中图分类号:V261.2 文献标志码:A doi:10.3969/j.issn.2095—1248.2013.04.004 Fault diagnosis of aero-engine main shaft bearings based on envelope spectrum analysis LIANG Xian—ya,SHA Yun—dong,LUAN Xiao—chi,ZHAO Feng-tong,ZHANG Qiang (Liaoning Key Laboratory of Advanced Test Technology for Aeronautical PropulsionSystem, Shenyang Aerospace University,Shenyang 1 10136) Abstract:A method combining wavelet transform and envelope spectrum is proposed for the collection of fault information and the fault diagnosis for the aero—engine main shaft bearings,taking advantage of the re— corded aero—engine vibration signals collected from the ground tests.First,the recorded aero—engine vibration signals are wavelet de—noised;then,the wavelet de—noised fault feature signals are extracted through wavelet ransfortm;lastly,the fault feature signals are analyzed through the Hilbert envelope spectrum to obtain the characteristic peak frequencies.After being applied in tests,this method is approved to be effective in the di— agnosis of the fault and the definite fault position of the aero—engine main shaft bearings,which may provide a virtual basis on the fault diagnosis of aero—-engine main shaft bearings and has a broad prospect of applica・・ tjon. Key words:wavelet transform;envelope spectrum;aero—engine;main shaft bearings;fault diagnosis 航空发动机作为飞机的心脏,其性能全面影 响到飞机的各项技术指标。支承传动系统是发动 理方法的发展,小波理论作为适用于研究非平稳 信号的时频分析工具广泛应用到轴承故障诊断之 中。文献3利用小波变换对信号的奇异性进行研 究,有效地诊断出滚动轴承的故障及其故障类型。 Fred K.Choy等 针对滚动轴承内圈和滚珠两种 机的重要组成部分,滚动轴承是航空发动机支承 传动系统中的关键环节。针对滚动轴承振动信号 的特性,人们相继引进了冲击脉冲 和涡流传感 器 的分析方法。近些年来,随着信号分析和处 收稿日期:2013—01—10 故障形式,对采集到的振动信号进行频谱分析、 作者简介:梁先芽(1987一),男,湖南耒阳人,在读硕士,主要研究方向:航空发动机强度、振动及噪声,E-mail:Liangxy24@yeah.net 沙云东(1966一),男,黑龙江阿城人,教授,主要研究方向:航空发动机强度,振动及噪声。 第4期 梁先芽,等:基于包络谱分析的航空发动机主轴轴承故障诊断 l9 Wigner.Ville和小波变换的时频分析以及庞加莱 映射的混沌振动分析。研究发现在轴旋转时,与 内圈缺陷相关的损伤区域周期性的出现。而对于 滚珠存在缺陷时,利用庞加莱映射可以有效地分 析研究由于轴承滚珠自身复杂的运动形式而造成 的振动信号的混沌现象。s.A.Mclnerny等 提 出了一套滚动轴承故障诊断专家系统。该系统能 够对轴承振动信号进行包络分析以及对故障轴承 特征和调制与解调的振幅进行相关分析。余光伟 等 利用离散小波变换和包络谱相结合的方法, 对采集到的轴承振动信号进行处理进而诊断出内 圈和滚球故障。赵美云等 提出了一种小波分 析和Hilbert包络分析相结合的方法对滚动轴承 声发射信号进行处理。通过5层小波分解和重构 以提取出故障特征频带,再对故障特征信息进行 包络谱分析,从功率谱中获得轴承外圈故障特征 频率的基频及倍频,以此诊断出滚动轴承外圈故 障。 本文主要研究的内容是某型涡扇发动机在地 面试验条件下,先对测得的发动机机匣振动信号 进行给定阀值技术的小波降噪处理,然后通过小 波分解与重构对降噪后的信号进行分析以提取到 主轴轴承的故障信息,再对该故障特征频段进行 包络谱分析。对仿真信号进行分析很好地提取出 淹没在噪声之中的周期信息。对试验数据分析后 发现该方法较好地判断出主轴轴承是否存在故障 以及故障的具体位置,这为主轴轴承故障诊断提 供了重要的依据,具有重要的研究价值。 1波变换和包络分析 1.1小波变换 小波变换是对一个函数在空问和时间上进行 局部化的一种数学变换,通过平移母小波获得信 号的时间信息,通过缩放母小波的尺度获得信号 的频率特性。对母小波的平移和缩放操作是为计 算小波的系数,这些系数代表局部信号和小波之 间的相互关系。常见的有连续小波变换和离散小 波变换 J。 (1)连续小波变换(CWT) 由基本小波或母小波 (t)通过伸缩a和b平 移产生一个函数族{ (t)}称为小波基,则有 )(f)= f4 、 1 “/ (1) 式中:a为尺度因子,a>0,b为时移因子。 设 ( )∈Z2( ),定义x(f)的连续小波变换 为 (a,b)=[ , (f)] ( (2)离散小波变换(DWT) 在计算机进行计算时需要将连续小波离散 化,最常见的离散化是二进制,将尺度因子a和时 移因子b离散化为 ra=2 ‘【b=2 k k,J∈Z (3) ,且口口]。 1.2小波降噪 一维信号小波降噪过程主要包括以下三个步 骤: (1)信号的小波分解。选择一个小波并确定 分解的尺度,再进行分解计算; (2)小波分解高频系数的阀值量化。对于各 个分解尺度下的高频系数选择好一个阀值进行软 阀值量化处理。该步骤为核心部分,常用的有默 认、给定以及强制消噪处理三类。考虑到给定阀 值消噪的可信度高以及一些经验公式,故采用给 定阀值消噪处理振动信号; (3)一维小波重构。根据小波分解的底层低 频系数和各层高频系数进行一维小波重构 j。 1.3小波分解与重构 离散二进制小波变换是一种常用的方法,因 为,因此a=2 , ∈Z,适合在计算机上进行高效地 运算。常采用Mallat塔式算法进行离散二进制 小波变换的计算,设信号 ( )的离散序列为 (n) (n=1,2,3,…,Ⅳ),其离散二进制小波变换为 C0k= ,、 cjt=,∑cj_1,.h I(七:0,l 2..,N一1)(4) =, ∑Cj-1,ng JI 式中: 为信号的离散采样数据,Ⅳ为离散采样点 数, 为分解的层数,h,g为滤波器的脉冲响应,用 于信号的分解,Cj, 为信号的逼近系数,d 为信号 的细节系数。 相应的信号重构算法为 cj =, ∑Cj一 h +∑ 一 g (k: 0,1,2,…,Ⅳ一1) (5) 式中: ,g 为滤波器的脉冲响应,用于信号的重 构 。 1.4包络分析 包络分析主要包括共振信号拾取和包络检波 20 沈阳航空航天大学学报 第30卷 两个过程,即拾取共振峰(传感器或机械系统固 有频率)附近处或高频谐振器谐振频率处的信 号,得到高信噪比的高频振动信号,再通过包络检 波的方法得到包络波形,包络波形中包含有与滚 动轴承故障诊断相关的特征信息。 通常使用Hilbert变换实现信号的解包络提 取故障信息。Hilbert变换定义是: +∞ 障诊断之中,首先利用该方法对仿真信号进行分 析。设,. 为仿真故障轴承冲击信号,1"2为干扰信 号,且仿真信号公式为: ,.1=sgn(sint) r2=rand(1,n) (7) (8) 仿真信号模拟滚动轴承在振源r1的振动下, 受到振源 的干扰,时域波形如图1(a),1(b)所 示。利用本文提出的小波包络方法对混合信号进 行分析。首先进行小波分解重构,结果如图1 / 、 (;): (t):一1:一1 f丛 t ! t—r (6) 利用Hilbert变换进行信号包络时的原理是 让测试信号产生一个90。的相移,从而与原信号 构成解析信号来进行包络,进行谱分析可以得到 清晰的故障信息"J。 (c)。再进行包络分析处理,结果如图1(d)。显 然频谱中出现方波的周期特征频率,由此可以说 明该方法的能够在噪声干扰的情况之下提取轴承 故障特征信息以用来进行故障诊断。 2仿真信号试验及分析 为验证提出的分析方法能够用于滚动轴承故 3 试验分析 滚动轴承在运转时会产生一些特有的振动 3 吕 、一 2 迥 1 1.0 0 罂 一1 l删E O l00 200 300 Time/ms (a)仿真信号时域图 (b)混合信号时域波形 Time/ms 频率f/Hz (d)小波分析包络谱图 (c)小波分析细节信号图 图1仿真信号的小波变换及其包络功率谱 频率,而这些频率在一定程度上表现了轴承故障 = 水平。某型航空发动机在地面检测条件下,主轴 ( +D bc。 )z (10) 轴承故障特征频率计算公式为 : 滚动体上某一固定点与外圈或内圈接触频率 为 z个滚动体与外圈上某一固定点接触频率为 =— 一 (1一/Y bc0 )z 、 一c , (9) z个滚动体与内圈上某一固定点接触频率为 第4期 梁先芽,等:基于包络谱分析的航空发动机主轴轴承故障诊断 21 ,c= ( 一 Dbc。 ) (12) 某型发动机高压涡轮后轴承的基本参数为: Z=34个直径为D =8 ITLII'I的滚棒组成的轴承, 其轴承节圆直径为Dc:148 nlin,接触角为p= 保持转速为.厂r=120 r/min的稳态过程,试验过程 中外圈转速视为零。计算得故障特征频率:外圈 32.2 Hz,内圈35.8 Hz,滚棒18.4 Hz。 对某次试验所测得数据进行小波降噪处理。 其中采样频率 =1 024 Hz,采样数据长度为N= 0。,外圈半径为140 iilin,内圈半径111 mid。地面 试验过程中,发动机在低压转子驱动器的驱动下 O.15 O.10 一0.05 0 一2 048。小波降噪过程中采用给定阀值消噪,消噪 处理前后的时域和频域如图2和图3。 ×10一 10 8 邑6 4 2 0.05 0 10 一-0.1 5 0.20 —土 O.5 1.0 1.5 2.O 一.. 0 5O 1OO 15O 200 250 300 350 400 450 500 时间(s) (a)原始信号时域图图2原始信号时域和频域波形 频率f/Hz (b)NNNg-N ̄ 1O 8 旦6 墨4 -3 42 0 O.5 1.O 1.5 2.O 一O 时间(S) (a)降噪后时域图 图3小波降噪后时域和频域波形 频率ffHz (b)降噪后频域图 对降噪后的信号进行4层小波分解,分解后 各尺度重构波形如图4(a)。其中d1一d4分别表 示1—4层细节信号时域波形。为了提取滚动轴 Hilbert包络变换并做谱分析,且采样频率 = 1 024 Hz,故频率分辨率为0.5 Hz。分析结果如 图4(b)。 承外圈故障特征频率,对第4层细节信号进行 名 5名 x 1=0==-3 一 — 州 二 脚 _. 忡静 _ ~~] { o.o;?…… QQ…… Q一一l …一 旦 一~ 0 ・穹o’050 一-05I 一一1 打f一一r …T …— i —— —南。 。 j6…~…… ~…… ¨r 曲“6………1 ~…… Jd …~ [00 o o.2 0 500 1 000 1 500 2 000 2 500 35 40 45 5O 55 60 65 … ] 采样点数/爪 频率f/Hz (a)小波分析细节系数波形 (b)细节系数的包络功率谱 图4 降噪后信号4层小波分析的细节系数波形及第4层信号的包络功率谱 对另一次试验所采集的数据进行相同处理。 其中采样频率 =1 024 Hz,采样数据长度为N= 小波分解的细节系数波形和第4层细节信号的包 络谱图如图5。 5 120。谱分析时的频率分辨率为0.2 Hz。4层 在实际试验过程中,发动机低压转子转速存 22 沈阳航空航天大学学报 第30卷 在一定的波动。两次试验数据的上述方法分析得 逐渐下降。由此可以初步判断该轴承外圈存在故 障,这和发动机台架试验过程中出现瓦基姆报警 结果相吻合。 到的峰值频率均为32 Hz,如图4和图5,其与轴 承故障特征频率32.2 Hz间的差异在误差范围之 内,而且谱线的幅值在远离各阶故障特征频率时 ×10一’ 名 一6 , 号0.0O5 0.05 0.05 勺 O 0.05 0.2 --毛—0 0.2 0 0.02 1 000 2 000 3 000 4 000 5 000 6 000 35 40 45 50 55 60 65 采样点数/个 频率f/Hz (a)小波分析细节系数波形 (b)第4层细节系数的包络功率谱 图5 降噪后信号4层小波分析的细节系数波形及第4层信号的包络功率谱 SAE l3th Proposion Conference,1977:l一7. 4结论 0 O O 0 O O 一 ∞\Ⅲ) 粤 ]jKim P Y.Proximity transducer technique for bearing health monitoringl J].J.Propulsion,1986,3(1):84— 89. Mallat S.A Theory for multi—resolution signal decom— ]j 1J];]J、j 鸺% 在航空发动机主轴轴承测试过程中,所测得 的信号为非平稳信号。通过分析发动机机匣振动 信号,采用小波变换和包络分析提取主轴轴承故 障特征信息,经过试验表明具有如下结论: (1)某型航空发动机在地面检测的条件下, 针对机匣振动信号进行小波降噪,能够有效地提 position:the wavelet representation signal decomposi— tion l J 1.IEEE Transon Pattem Analysis and Machine Intelligence。1989.11(2):257—263. Fred K Choy,Lei Wang.On—line vibration monitoring of ball bearing damage using all experimental test rig 1 J 1.Journal of Propulsion and Power,2007,23(3): 629—636. McInerny S A。Dai Y.Basic vibration signal processing 取表征主轴轴承故障的特征信息; (2)在一定条件下,针对小波降噪后的信号 进行小波分解与重构,对重构后的特征频段进行 包络谱分析可以较好地诊断主轴轴承故障; (3)针对发动机机匣振动信号,小波变换和 for bearing fault detectionI J I.IEEE Transactions on Education,2003。46(1):149~156. 包络分析相结合的信号分析方法在一定程度上适 用于航空发动机主轴轴承故障诊断。 参考文献(References): 余光伟,郑敏.小波变换在滚动轴承故障分析中的应 用[J].轴承,201l,7:37—40. 赵美云,李力,陈保家.滚动轴承故障声发射信号的 小波包络谱分析[J].轴承,2008,4:38—41. 范显峰,姜兴渭.基于小波包变换的信号去噪方法研 究[J].哈尔滨工业大学学报,2003,7:809—811. Jason R.Stack.Thomas G.Habetler.Fault classi丘cation nd faulat signature production for rolling element bear— ings in electric machines lJ 1.IEEE。Symposium on Di— agnostics for Electric Machines.Power Electronics and Drives,2003:172—176. [1]George J A,Andres R M.Helicopter bearing failure de— tection utilizing shock pulse techniques[C].AIAA/ (责任编辑:刘划英文审校:宋晓英)