发动机试验动态数据处理方法综述.pdf
《发动机试验动态数据处理方法综述.pdf》由会员分享,可在线阅读,更多相关《发动机试验动态数据处理方法综述.pdf(9页珍藏版)》请在咨信网上搜索。
1、宇航计测技术2023 年2023 年 8 月第 43 卷 第 4 期宇航计测技术Journal of Astronautic Metrology and MeasurementAug.2023Vol.43No.4文章编号:1000-7202(2023)04-0066-09DOI:10.12060/j.issn.1000-7202.2023.04.11发动机试验动态数据处理方法综述刘子星1,高 俊2,宋炯亮2,林 萌2,杨 懿2(1.北京动力机械研究所,北京 10074;2.北京航天试验技术研究所,北京 100074)摘要:发动机试验动态数据通常包含关键部位振动、脉动压力等重要信息,是评判发动机
2、性能、零部件特性和优化设计的重要参考指标。总结了快速傅里叶变换、小波分析、小波包分析、Hilbert-Huang 变换等几种适用于动态数据处理方法的优缺点,分别介绍了上述几种方法以及多种方法相结合后在动态数据分析和处理方面的应用。关键词:发动机试验;动态数据;FFT;小波分析;小波包分析;希尔伯特-黄变换中图分类号:V438+.4文献标识码:AReview of Dynamic Data Processing Methods for Engine TestLIU Zi-xing1,Gao Jun2,SONG Jiong-liang2,LIN Meng2,YANG Yi2(1.Beijing P
3、ower Machinery Institute,Beijing 100074,China;2.Beijing Institute of Aerospace Testing Technology,Beijing 100074,China)Abstract:Engine test dynamic data usually contains important information such as vibration and pulsating pressure ofkey parts,which is an important reference index for evaluating en
4、gine performance,component characteristics and optimaldesign.The advantages and disadvantages of several methods suitable for dynamic data processing are summarized,such asfast Fourier transform,wavelet analysis,wavelet packet analysis,Hilbert-Huang transform,and the application of thesemethods and
5、their combination in dynamic data analysis and processing are introduced.Key words:Engine test;Dynamic data FFT;Wavelet analysis;Wavelet packet analysis;Hilbert-Huang transform收稿日期:2023-01-11,修回日期:2023-03-22作者简介:刘子星(1985-),男,高级工程师,本科,主要研究方向:航空发动机控制技术。1 引 言发动机试验是解决发动机研制过程所遇到的问题、推动理论研究和完善发展发动机设计理论的基本途
6、径,也是分析发动机技术参数和评估发动机性能和评判发动机能否定型的重要方法1。在发动机试验中,各种测量数据是研究发动机性能、定位各类型故障、计算可靠性等的重要依据。数据处理是将反映物体状态的物理量随时间变化的原始信号转换成某些适合于工程设计、试验和分析的简化量。数据处理贯穿于发动机设计、试验和生产的各个阶段,其目的是揭示动力学环境、制定试验条件、设计准则以及诊断故障和提供信息2。对发动机试车的动态参数进行分析是研究发动机各项性能指标的重要依据。例如,分析推进剂的流量数据能够精确计算发动机比冲和混合比。分析推力和稳态压力能够研究发动机实际工作性能以及发动机与试车台之间的协调性。推力室脉第 4 期发
7、动机试验动态数据处理方法综述动压力是研究发动机启动特性与不稳定燃烧、管路特性和定位发动机故障原因的重要依据。振动参数是研究推力室、燃气发生器、涡轮泵等关键部件结构性能、不稳定燃烧状态等非常重要的参考依据。因此,研究动态数据的分析和处理方法,不仅关系到发动机零部件的性能评价,也是准确评估发动机总体性能、发动机和试车台匹配程度等的重要依据。目前,随着数字信号分析技术的发展,业内研究人员对动态数据分析已经由最初单一的傅里叶变换分析发展到小波分析、小波包分析、Hilbert-Huang 变换、人工神经网以及多种分析方法相结合的灵活运用和推广3,4。在本文中,总结了上述多种方法的理论特征,梳理了该方法在
8、动态数据分析领域的应用情况,为深度挖掘试验数据信息、提高发动机性能和提高发动机设计能力提供了参考依据。2 发动机试验动态数据处理方法2.1 FFT 方法快速傅里叶变换(FFT,1965)的基本思路是将时域采样点数分解为组合因子,按照组合因子将整个数据序列分割为若干子序列。先计算各子序列的离散傅里叶变换,然后将子序列变换结果加以合成,从而得到整个序列的离散傅里叶变换2。该方法的局限性在于缺乏分析时域内局部频率数据的能力。在工程应用中,采用计算机对时域连续傅里叶变换进行离散化时会不可避免地产生一些问题:(1)时域离散化引起频域的周期延拓,很可能产生频率混淆;(2)时域的截断引起泄漏;(3)频域的离
9、散化引起时域的周期延拓。为解决上述问题,诸多研究人员制定了相应的解决办法,如加入抗混滤波器、窗函数、提高采样率等2。目前,FFT 方法已经广泛应用于航空、航天发动机试验振动数据和脉动压力数据等的分析处理中5-8。胡海峰等根据 R 语言在统计分析方面的优势,在 VC 软件平台结合 R 语言的编程实现发动机数据分析处理。采用 FFT 方法分析点火启动及拖尾段的含噪波动数据,根据时域内振动信号的频谱图,识别信号高低频的有用信息9。孙百红等系统归纳了 FFT 方法对液体火箭发动机试验振动数据进行分析的步骤和路线。采用FFT 分析方法对某型号液体火箭发动机试验出现故障的稳定段振动数据进行频谱分析,通过对
10、比正常数据的幅频特性,成功定位某试验件结构破损的故障。通过对振动数据 RMS(Root-Mean-Square,均方根)特性分析,绘制振动数据特征频段全程瀑布图,能够有效识别试车中的故障10。牛伟等将航空发动机转速高齿信号作为分析发动机振动状态的基础,提出了基于 FPGA 的振动信号采样设计流程:采样-滤波抽样-三次样条插值-信号抽取-傅里叶变换,制定了振动数据频域特征提取处理流程。通过对试验数据进行分析,证明该方法不仅运算简单、分析精度高、实时性强,而且有效解决了 FFT 方法分析过程中的频谱泄露和栅栏效应等问题11。杨懿等针对液体火箭发动机脉动压力数据采样率高和系统响应频率快的特点,采用
11、FFT 方法对脉动压力数据进行分析。根据脉动压力数据的频谱特性,结合推力室轴、径向和切向三个振动测点的数据,有效定位发动机不稳定燃烧的声振类型12。陈子豪等为了验证固体火箭发动机试车过程中因不稳定燃烧所产生的压强异常升高和发动机推力值异常振荡的问题,采用 FFT 方法分析推力室压强的幅频特性和燃烧器的试验数据,得到压强振荡的倍频关系,结合声学仿真的结果,证明该方法能够有效分析发动机不稳定燃烧的机理13。洪亮等采用 FFT 方法,通过统计分析转子在不同转速下的倍频特征和振动量级之间的关系以及转子临界试验中的幅频特性,得出了转子发生故障的特征和原因,为优化和提高发动机设计提供了重要的参考依据14。
12、钟永恒等针对航空发动机外部管路流固耦合振动是导致发动机外部管路故障的重要因素之一,根据所推导的直管管路流固耦合偏微分方程建立管路系统物理模型,采用 C+语言开发了仿真分析软件,研究 MOC-FFT 方法分析计算外部管路时频响应特性和固有频率。通过动态压力特性试验,验证了该模型的有效性和准确性15。发动机试车过程中,不同系统之间往往存在关联,因此,通常将不同测点的脉动压力数据、振动数据等结合起来进行分析,便于分析发动机的性能和故障定位。以某型号发动机试验为例,将振动测点数据和推76宇航计测技术2023 年力室脉动压力数据12的分析相结合。脉动数据采样率为10 kHz,稳定段0.02 s 的数据如
13、图1 所示。图 1 推力室脉动压力数据图Fig.1 Fluctuating pressure data diagram of thrust chamber采用 FFT 方法对脉动压力数据进行频谱分析,0 11 000 Hz 频率段的频率-幅值特性曲线如图 2所示。图 2 推力室脉动压力数据 FFT 频谱图Fig.2 FFT spectrum of fluctuating pressure data分析图 2 的幅频曲线,推力室脉动压力的幅值峰值主要集中在 05 500 Hz 的频率段内,5 50010 000 Hz的高频段幅值曲线较为平坦。在中、低频段,5003 500 Hz 幅值存在多个小幅
14、度波动,可能与发动机工作倍频有关系。5 000 Hz 左右存在一个高峰值波峰,可能与发动机中频共振和不稳定燃烧有关。对于脉动压力数据频谱分析需要结合发动机振动数据进一步分析。采用 FFT 方法对推力室轴向、径向和切向三个振动测点数据进行分析,其中径向测点 ab1-2 的频谱分析结果如图 3 所示。图 3 推力室径向振动测点 FFT 频谱图Fig.3 FFT spectrum of radial vibration measuring point从图 3 可知,径向振动测点在 1 0005 000 Hz低频率段中存在多个幅值峰值,5 000 Hz 左右也存在一个高峰值波峰,说明脉动压力数据与径向
15、振动测点存在关联。对脉动压力、推力室轴向、径向和切向振动测点数据 FFT 频谱分析结果如表 1 所示。表 1 推力室振动测点与脉动压力幅频数据表Tab.1 Data of vibration measuring points andfluctuating pressure amplitude frequency of thrust chamberPfs 幅值峰值对应的频率/Hzab1-1(轴向)ab1-2(径向)ab1-3(切向)904-1 611-2 243-2 884-3 685-4 822-9 643-注:表示在振动数据 FFT 分析的幅值峰值的频率值与 Pfs吻合。从表 1 可知,在1
16、 611 Hz、2 884 Hz、4 822 Hz 和9 643 Hz 时,振动径向测点的幅频峰值特性与脉动压力测点的幅频峰值特性存在耦合关系。轴向和切向振动测点与脉动压力耦合点较少,由此可以推断推力室发生高频横向声振。从上述分析过程可知,采用 FFT 分析方法,结合脉动压力、振动等多类型的动态数据能够为判断发动机不稳定燃烧状态和评价发动机性能提供有86第 4 期发动机试验动态数据处理方法综述力的参考依据。2.2 小波分析方法小波分析方法以傅里叶变换为基础,具有良好的自适应优势。比利时数学家 Ingrid Daubechies 系统总结了上世纪 80 年代以来诸多研究人员以及自己对小波的先进研
17、究成果,编写了文献16,对小波分析的基本理论、性质、适用领域和使用规则进行了详细介绍,使得小波分析成为一门真正的应用学科,并成为国际研究热点。小波分析在发动机动态数据处理中的应用主要集中在数据去噪、奇异点检测分析等方面。符娆等针对航空发动机转、静子发生碰摩时,振动信号频率产生急剧变化的特点,利用小波分析在局部信息中具有高分辨率的优势,采用阈值法先对原始含噪振动信号进行去噪处理,再提取有价值信号中的幅值特征和能量信息。通过对试验数据分析表明,降噪效果较好,能够准确跟踪信号中的突变信息17。文壁等重点分析了涡扇发动机起动过程中转子做功脉动压力、噪声脉动压力、失速团引起的脉动压力和燃烧波动引起的脉动
18、压力四种成分,利用小波分析在时频分析领域的优势,对含噪信号进行分解,对多次起动过程的频谱图,结合起动过程的压力脉动现象进行分析,结果显示小波分析方法能够清晰识别压力测点不同转速下的脉动频率,定位脉动压力异常的原因18。耿纪洲等运用小波分析和分形算法相结合的方法对发动机轴承磨损的含噪数据进行去噪处理,分析结果显示该方法能够清晰提取轴承的故障特征以及判断最佳的故障诊断部位19。肖云魁等针对小波分析理论中 Lipschitz 指数不能清晰定位汽车发动机曲轴轴承工作异常过程中的奇异点和对奇异点能量大小的描述不够精确的问题,先对试验数据进行去噪,再采用模极大值多尺度分析方法,对轴承故障数据进行分析。仿真
19、和试验验证的结果显示,该方法能够精确定位试验过程中的奇异点位置20。戴屹梅总结了液体火箭发动机涡轮泵低温轴承研制以及试验中的性能参数和故障模式特点,介绍了傅里叶变换、小波分析等多种数据分析方法,建立了仿真模型,开发集成了短时傅里叶变换、小波分析等多种方法的数据分析系统,对行业内发动机试验涡轮泵轴承数据分析提供了重要的参考依据21。杨懿等介绍了小波去噪的基本原理、四种阈值去噪的基本方法以及具体实施的基本步骤。将信噪比作为选择小波基函数的参考标准,将相关系数作为去噪结果的评价标准。通过对发动机试验含噪流量数据进行去噪分析,证明小波分析方法能够有效消除流量数据中的含噪分量22。白泉等介绍了小波分析在
20、实际应用中的基本步骤和方法。小波分析中所用到的小波基函数具有非唯一性,采用不同小波基函数分析同样的数据时将会产生不同的结果。因此,研究人员列举了 L-P 小波、Meyer 小波、Db 小波等多种小波基函数在支撑长度、对称性和正则性等方面的特征性质,根据小波基函数性质和被分析数据的特征选择恰当的小波基函数23。采用小波分析方法对某型号发动机试验关机段推力室脉动压力数据的奇异点特征进行分析。数据样本数量为4 000 个,采样率为10 kHz,选择DB9 小波小波基函数,分解层数为7 层,分解结果如图4 所示。图 4 DB9 小波 7 层分解结果Fig.4 DB9 wavelet 7-level d
21、ecomposition results对图 4 中 d1d7 的分解结果进行分析可知,小波分析能够准确定位数据奇异点在时域内的信息。d1d5 层的分解重构效果较 d6d7 好,能够清晰观测到时域内不连续点和幅值变化信息。d1d4 低层96宇航计测技术2023 年高频信号中突变能量较 d5d7 层高,包含了原始信号中大部分能量突变成分。随着分解层数的增加,小波分解过滤的信息较多,d6d7 层存在一定程度的失真,曲线较为平缓,对分析结果的影响不大。在运用小波分析方法对数据进行分析的过程中,可以将流量、振动等其他数据有效结合,为判断不稳定燃烧状态及发动机整体和零部件性能提供重要的参考依据。2.3
22、小波包分析方法小波包分析方法是小波概念的优化和提升。该方法克服了小波变换在中、高频段频率局部性分析能力差的缺陷,对信号中的高、低频部分进行二次分解,从而进一步提高了信号的时频分辨率24,25。小波包分析在信号去噪、奇异点定位、数据压缩、图形处理和故障检测等领域有着广泛的应用。杨懿等介绍了小波包分析的基本原理和小波包阈值去噪的基本步骤。在对氢氧火箭发动机试验流量数据进行去噪过程中,将信噪比和相关系数作为去噪效果的评价标准,取得了良好的去噪效果26。杨懿等采用小波包阈值去噪与分解重构算法分析液体火箭发动机试验脉动压力数据。分析了脉动压力数据在高、低频段内的系数特征和能量分布情况。结合相应振动测点
23、FFT 频谱特性对发动机在多个工况内的工作状态进行分析27。崔建国等针对主燃油泵轴承损坏以及调节器故障引发的喘振问题,将小波包多层分解后的能量信号作为原始信号的特征向量,结合神经网络中的极限学习算法,建立主燃油泵的故障诊断模型,并将该模型与 BP 神经网络模型进行对比验证。结果显示,该模型能够更加准确、迅速地实现燃油泵故障诊断28。任学平等采用小波包分析方法对汽车发动机运行的含噪振动信号进行去噪处理。根据计算机仿真结果选定合适的小波基函数和分解层数,对降噪后的信号进行残差分析和统计量分析。降噪结果显示该方法能有效区分信号中的突变和噪声成分,具有良好的信噪分离效果29。马建仓等为有效判断航空发动
24、机故障类型,统计了 5 种常见的航空发动机工作故障特征,将峭度指标、裕度指标、脉冲指标和有效值等时域指标作为评判发动机是否正常工作的参考标准。根据分解结果中高低频带的能量信息,绘制多层频带的相对能量分布图,并结合故障特征和时域指标,实现航空发动机故障特征的准确定位30。小波包分析方法在动态数据分析中的使用方法和流程与小波分析基本一致。在采用 MATLAB 实现编程分析中,调用的分析程序和基函数稍有不同。以某型号发动机试验流量数据为例26,采用小波包分析方法对含噪流量数据进行去噪分析。样本点数为 8 000 个,数据采集系统采样率为 10 kHz。原始数据如图 5 所示。图 5 含噪流量数据图F
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 发动机 试验 动态 数据处理 方法 综述
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。