一种基于emd的低湍流度信号处理分析方法-朱博.pdf

《一种基于emd的低湍流度信号处理分析方法-朱博.pdf》由会员分享，可在线阅读，更多相关《一种基于emd的低湍流度信号处理分析方法-朱博.pdf（6页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、 第30卷 第5期实验流体力学Vol. 30,No. 5 2016年10月Journal of Experiments in Fluid M echanics Oct. , 2016 文章编号:1672- 9897(2016)05- 0074- 06 doi:10. 11729/syltlx20150148 一种基于EM D的低湍流度信号处理分析方法朱 博* ,彭 强,汤更生(中国空气动力研究与发展中心设备设计及测试技术研究所,四川绵阳 621000)摘要:采用热线风速仪在3座典型低速风洞中进行了流场湍流度测量,这3座低速风洞包括1个低湍流风洞、1个常规闭口风洞和1个开口射流风洞。针对湍流信号

2、通常受噪声干扰的问题,在湍流度值处理中引入了经验模式分解( EM D)自适应滤波和HHT时频谱分析方法。将EM D方法与其他几种湍流度值处理方法进行了比较,包括带通滤波方法( BP) 、电磁噪声解耦方法( ENC)和高通惯性衰滤波方法( HPIA) 。采用EM D方法测得低湍流风洞的湍流度值,在流场速度30 100 m/ s的范围内小于0. 05% 。采用HHT方法完成了脉动速度信号的时频分析,分析发现开口风洞试验段的脉动速度HHT时频谱有突出的低频信号。所构建的EM D自适应滤波器可以有效控制噪声对热线输出信号的影响,是一种有效的低湍流度信号处理方法。关键词:低湍流度;经验模式分解;低速风洞

3、;开口风洞;热线风速仪中图分类号:O357. 5+ 4 文献标识码:A收稿日期:2015- 12- 08;修订日期:2016- 02- 27*通信作者E- mail: bobjou 139. com引用格式:Zhu B , Peng Q , Tang G S . Digital signal process of low turbulence intensity based on EM D . Journal of Experiments in Fluid M echanics, 2016 ,30 ( 5 ) : 74- 79 .朱 博,彭 强,汤更生.一种基于EM D的低湍流度信号处理分析方

4、法.实验流体力学, 2016 , 30 ( 5 ) : 74- 79 .Digital signal process of low turbulence intensity based on EM DZhu Bo* , Peng Qiang, Tang Gengsheng( Facility Design and Instrumentation Institute, China Aerodynamics Research and DevelopmentCenter, M ianyang Sichuan 621000, China)Abstract: Hot wire anemometer is

5、 used in three typical low speed wind tunnels to measurethe flow field turbulence intensity, including a low turbulence wind tunnel, a general closed windtunnel and an open jet wind tunnel. To deal with the disturbance to the turbulence signal causedby noise, an EM D( Empirical M ode Decomposition)

6、self- adaptive filter and HHT( Hilbert- HuangTransform) time- frequency analysis methodology are introduced in turbulence intensity data pro-cessing. The EM D method is compared with several other turbulence intensity processing meth-ods, including BPF( Band Pass Filter) , ENC( Electronic Noise Corr

7、ection) and HPIA( High PassInertia Attenuation Filter) . Using the EM D filter, the turbulence intensity in the low turbulencewind tunnel test section is measured and found to be less than 0. 05% in the flow speed range of30 100m/ s. Using HHT method, we complete the time- frequency analysis of the

8、flow fluctua-tion signal, and the low frequency flow flucturation signal is found predominatnt in the open jetwind tunnel. The designed EM D self- adaptive filter is effective at controlling the influence of thenoise on the hot wire output signal, and therefore it is a practical method for low turbu

9、lence in-tensity signal processing.Keywords: low turbulence intensity; empirical mode decomposition; low speed wind tunnel;open jet wind tunnel; hot wire anemometer0 引 言由于湍流度对飞行器绕流的流动结构有较大影响,可直接影响翼型等气动部件的边界层转捩位置及分离位置,从而改变飞行器的升阻比和俯仰力矩特性 1- 2 ,因此,湍流度是风洞流场校测的重要指标,需要在飞行器模型的风洞试验中对风洞湍流度进行准确测量,并对其影响进行修正以提高

10、风洞试验数据精准度。目前,测量流场湍流度的最理想仪器是热线风速仪(简称热线) ,但是由于测量目标为动态信号,导致测量结果容易受到干扰,特别对于低湍流度测量而言,更容易受到电磁干扰 3- 4和支架干扰 5 。因此,在低速风洞低湍流度测量中,需要对数据进行有效的信号处理 6- 11 。万方数据 自20世纪30年代末至今,为推动湍流机理、边界层转捩和翼型层流化研究的不断发展,国内外已建造约30座低湍流度风洞 4并开展了大量研究工作,研究提出多种信号处理方法,降低干扰噪声的影响从而提高了低湍流测量的准确性。文献 6采用0. 1Hz 1kHz带通滤波BP( Band Pass)方法,测得试验段湍流度0.

11、 069% ( 5m/ s) 、 0. 086% ( 15m/ s) ;文献 7采用10Hz 1kHz的BP方法,测得开口试验段喷口位置X/ D= 0. 042的核心流湍流度0. 02% ( 33ft/ s) ;文献 8采用1Hz 10kHz的BP方法,测得试验段湍流度0. 05% (动压200m/ s) ;文献 9采用1Hz4kHz的BP方法测得试验段湍流度0. 05% ( 150 ft/ s) ;文献 10采用电磁噪声解耦ENC( Electronic NoiseCorrection)方法,测得试验段湍流度0. 025% ( M a0. 05) 、 0. 3% ( M a0. 2) ;文献

12、 11采用声波干扰去相关方法,测得试验段湍流度0. 02% ( 18 53m/ s) ;文献 12提出了用尺度间变化的门限值来抑制带噪湍流信号在不同尺度上的噪声子波系数,从而实现了在重构湍流信号中消除噪声;文献 3采用快速收敛迭代法对测量数据进行信号处理并考察处理前后统计湍流量的变化,发现噪声的存在不仅对湍流小尺度统计量的估算影响显著,也对与大尺度关联的量产生一定的影响;文献 13采用高通惯性衰减滤波方法HPIA( High Pass Inertia Attenuation Filter)测得试验段湍流度0. 05% ( 30 100m/ s) 。经验模式分解EM D( Empirical M

13、 ode Decompo-sition)方法是美国学者Huang等在1998年提出的一种信号处理技术 14 ,成功应用于海洋、地震、机械振动分析、声音处理、纹理分析、影像滤波以及气象等领域的数据分析中 15 。通过EM D分解得到的基本分量称为本征模函数( Intrinsic M ode Function,IM F)分量,对所有的IM F分量进行HHT( Hilbert-Huang Transform)变换,可以得到整个信号的HHT时频谱。通过HHT时频谱不但能够知道信号具有哪些频率成分以及对应成分的大小,更关键的是还能够知道各个频率成分出现的时刻,从而将时域和频域特征很好地联系起来。本文根据

14、低速风洞低湍流度信号的特征,通过实验构建了一种基于EM D的自适应滤波器进行信号处理,采用同一阈值处理3种低速风洞的湍流度数据,避免了由于不同低速风洞的干扰频率差异,导致的滤波参数不统一的问题,便于比较不同低速风洞的湍流度值;对热线信号进行了HHT时频谱分析,分析发现开口风洞试验段存在明显的低频速度脉动信号。1 风洞设备及热线仪1. 1 风洞设备本文主要试验在低湍流风洞、常规闭口风洞和开口射流风洞等3种典型低速风洞上进行。低湍流风洞稳定段设置6层阻尼网,流场湍流度设计指标为0. 05% 。开口风洞流场湍流度设计指标为0. 2% 。常规低速风洞稳定段设置3层阻尼网,流场湍流度设计指标为0. 5%

15、 。1. 2 热线仪热线仪采用丹麦丹迪公司的Streamline系统,测量探头主要有55P11一维探头和55P61二维探头,测点布置于风洞试验段的中心,系统配置的采集卡分辨率为16位。数据采样频率为20kHz,采样时间为30s。2 低湍流度信号的EM D处理和自适应滤波根据EM D方法 14- 17对采样信号x( t)进行分解,可表示为:x( t) = ni= 1ci ( t) + r( t) ( 1)式中: ci ( t)表示分解的第i个本征模函数分量,r( t)表示残余信号分量。对ci ( t)进行HHT变换可以表示为:H( ci ( t) ) = 1+ - ci ( )t - d ( 2

16、)构造解析函数:zi ( t) = ci ( t) + jH ( ci ( t) ) = Ai ( t) ej i ( t) = Ai ( t) ej fi ( t) d( t)( 3)式中:瞬时幅值为Ai ( t) = ( ci ( t) ) 2 + ( H( ci ( t) ) ) 2 ( 4)瞬时相位为 i ( t) = arctan H( ci ( t) )ci ( t) ( 5)瞬时频率为fi ( t) = 12 d( i ( t) )dt ( 6)瞬时幅值Ai ( t)和瞬时频率fi ( t)同时为时间t的函数, Ai ( t) , fi ( t) ,( )t构成了HHT时频谱。由

17、于EM D方法是按照高频到低频分解信号,因此,只要对低频分解信号进行重构,就可以去掉低湍流度信号中的高频干扰:x ( t) = x( t) - m- 1i= 1ci ( t) ( 7)57 第5期 朱 博等:一种基于EM D的低湍流度信号处理分析方法 万方数据 又由于高斯白噪声信号的IM F分量的能量密度与其平均周期的乘积为一常量 16 ,即Em Tm = const ( 8)因此,可以通过查验相邻IM F的const值变化情况来判定分解信号是否受到了白噪声干扰,从而确定分解信号重构的IM F分量的级数m, m值可以是EM D从高频到低频分解过程中满足下式的第一个解(试验中发现,热线信号容易受

18、低频白噪声影响,可导致低频部分存在第2个解,若将第2个解也作为判定级数,从而滤除低频干扰信号,可以进一步构建一个自适应的带通滤波器) :Rm = ( Em Tm - Em- 1 Tm- 1 ) / ( 1m - 1m- 1i= 1Ei Ti )Rm C ( 9)能量密度为Em = 1N Ni= 1cm( i ) 2 ( 10)平均周期为Tm = N 2Pm( 11)上述通过判定级数m值从而进行IM F分量的重构过程,可以理解为一种自适应滤波方法。其中,N为分解信号的数据总量, Pm为第m个分解信号的极值点总数, C值为常数, C直接影响到滤波的质量。表1是本文3种低速风洞中的能量密度与其平均周

19、期的乘积变化情况,可见, Rm值在2 3的范围内是一个突变区域, Rm值在这个范围内迅速增大。采用文献 13对热线高频干扰信号的自谱分析方法,认为Rm值的突变是干扰引起的,即Rm值大于2 3以上的信号受到了电磁干扰和支架干扰的影响,因此,对C取值为2作为EM D滤波器的阈值,采用这一阈值滤波后的信号频谱如图1所示,在图中,有效信号得到最大限度保留,高频干扰信号被抑制。表1 3种低速风洞中的能量密度与其平均周期的乘积Table 1 A list of Em Tm in the 3 kinds of low speed wind tunnel风洞类型R1 R2 R3 R4 R5 R6 R7 R8

20、R9 R10 R11 R12 R13 R14 R15 R16低湍流0. 44 0. 05 0. 10 0. 04 0. 22 0. 21 0. 03 0. 01 1. 05 33. 79 37. 20 47. 05 27. 88 27. 62 7. 76 3. 96常规低速0. 34 0. 21 0. 03 0. 17 0. 07 0. 01 0. 53 0. 20 0. 50 1. 99 16. 60 25. 63 36. 38 11. 32 20. 11 4. 31开口风洞0. 65 0. 25 0. 03 0. 01 0. 02 0. 11 0. 03 0. 20 0. 09 0. 79

21、 9. 01 8. 68 451. 32 18. 75 13. 57 19. 573 湍流度测量结果和分析3. 1 几种滤波方法获得的脉动速度频谱比较从脉动速度的频谱图可以直观地比较不同方法对湍流度值的处理效果 13 。图1是低湍流风洞试验段风速为90m/ s时的脉动速度幅值谱,采用EM D自适应滤波、高通惯性滤波和频域BP方法 18等3种方法进行了信号处理。图2是图1采用3种方法获得的脉动速度的时域数据。由于低频干扰在风洞测试数据中普遍存在,低频脉动速度对转捩影响未知 8 ,一般数据处理都对低频干扰进行了高通滤波,高通滤波参数在0. 1 10Hz之间 6- 9 ,因此,本文考虑热线在较长时间

22、工作中的加热氧化和温度漂移等低频干扰影响,在上述3种方法中也同时采用了0. 5Hz高通滤波。从图1可见, 3种信号处理方法在0. 5 40Hz之间的处理结果接近,由于40Hz以下的幅值比40Hz以上的幅值大1至2个数量级,因此,这3种信号处理方法获得的结果也基本接近(见图2) ,但是,高频部分的处理对测量结果略有影响, EM D自适应滤波结果比0. 5 40Hz的BP结果更接近于高通惯性滤波结果。图1 低湍流风洞试验段脉动速度幅值谱比较Fig. 1 Comparition of flow velocity fluctuation amplitudespectrum in the low tur

23、bulence W T图2 脉动速度时域数据比较Fig. 2 Comparision of flow velocity fluctuation time data67实 验 流 体 力 学 ( 2016)第30卷 万方数据3. 2 几种方法测得的湍流度结果比较图3为采用4种处理方法对开口风洞试验段数据进行处理获得的湍流度结果比较。由图可见,因为开口试验段低频压力脉动较大, 0. 5Hz高通滤波对试验数据的影响较大,高频干扰信号对测量结果的影响相对较小,因此ENC方法 10测得的湍流度结果与原始值接近, 0. 5Hz 5kHz的BP、 0. 5Hz的HPIA和EM D自适应滤波结果接近。 BP、

24、 HPIA和EM D方法测得的开口风洞试验段湍流度值在除60m/ s外的大部分风速段都满足设计指标0. 2% 。从本文的使用情况看, BP方法最简洁,但不具备自适应性,而且截止频率以外的有效信号有可能被抑制。 ENC方法经过初始电磁环境标定后,原则上可用于类似的电磁环境,但仅对电磁干扰的抑制有作用,无法抑制支架干扰信号。 HPIA可较大限度地保留原始信号的高频特征,但不属于自适应滤波方法。 EM D方法采用统一阈值滤波,对不同环境中的干扰信号抑制效果较好,具有自适应性。图3 开口风洞试验段湍流度不同信号处理结果的比较Fig. 3 Comparison of results obtained b

25、y different kinds of dataprocessing methods in the open jet W T图4和5分别是低湍流风洞试验段和常规低速风洞的湍流度信号处理结果。从图中可见,采用EM D和HPIA方法测得的结果比较接近, EM D比HPIA的处理结果略小。低湍流风洞试验段湍流度值在流场速度100m/ s以下优于设计指标0. 05% ,常规低速风洞试验段湍流度值在所测试风速范围内优于设计指标0. 5% 。图4 低湍流风洞试验段湍流度信号处理结果Fig. 4 Turbulence intensity results in the low turbulence W T图

26、5 常规低速风洞闭口试验段湍流度信号处理结果Fig. 5 Turbulence intensity results in the general low speed W T3. 3 脉动速度信号的HHT时频分析图6是对图1中原始数据的所有IM F分量进行变换获得的0 200Hz信号的低湍流风洞试验段90m/ s风速HHT时频谱。从图中可见,在100Hz以下,气流脉动速度的频率和幅值随时间呈现一定的连续变化。在100 200Hz之间,这种连续变化逐渐变成离散状态。 EM D自适应滤波器在100 200Hz之间根据高斯白噪声的特点,确定了分解重构的IM F级数m并进行滤波。从图1的FFT幅值谱也可

27、见,EM D滤波获得的信号幅值在100Hz以后逐渐衰减。图6 低湍流风洞试验段脉动速度0 200Hz信号HHT时频谱Fig. 6 HHT spectrum of 0 200Hz flow velocity fluctu-ation in the low turbulence W T图7是常规低速风洞90m/ s风速时的0 200Hz脉动速度HHT时频谱,图8是开口风洞试验段90m/ s风速的0 200Hz脉动速度HHT时频谱。由图可见,常规低速风洞脉动速度的幅值和频率变化在50 200Hz范围内,比低湍流风洞试验段略丰富些。开口风洞试验段脉动速度幅值较大,而且主要是50Hz以下的低频信号,这种

28、信号一般是开口射流中的大涡结构引起的低频压力脉动和动压波动,在时频图中显示的幅值较为突出,是开口风洞调试和优化的重要指标。图9是图7的200 1000Hz高频部分,图10是图8的200 1000Hz高频部分。在高频部分可见, 2座风洞的速度脉动幅值和频率特征都比较接近,频率越高,信号被白噪声干扰越严重,信号幅值变化不连77 第5期 朱 博等:一种基于EM D的低湍流度信号处理分析方法 万方数据续,呈现为离散特征。图7 常规低速风洞试验段脉动速度0 200Hz信号HHT时频谱Fig. 7 HHT spectrum of 0 200Hz flow velocityfluctuation in th

29、e general low speed W T图8 开口风洞试验段脉动速度0 200Hz信号HHT时频谱Fig. 8 HHT spectrum of 0 200Hz flow velocityfluctuation in the open jet W T图9 常规低速风洞试验段脉动速度200 1000Hz信号HHT时频谱Fig. 9 HHT spectrum of 200 1000Hz flow velocityfluctuation in the general low speed W T图10 开口风洞试验段脉动速度200 1000Hz信号HHT时频谱Fig. 10 HHT spectru

30、m of 200 1000Hz flow velocityfluctuation in the open jet W T4 结 论本文给出了一种基于EM D的低湍流度信号处理分析方法,采用这种方法对3种低速风洞流场进行了测试,并与ENC、 HPIA和BP等方法进行了比较。( 1) EM D自适应滤波方法可根据信号中的白噪声参数确定分解重构信号的级数,其低湍流度信号处理结果略低于ENC、 HPIA和BP等方法测得的结果。( 2)采用同一EM D自适应滤波器分别测得3种低速风洞的流场湍流度,测量结果表明该方法适应性较好。低湍流风洞试验段流场的湍流度值在流场速度为30 100m/ s时优于设计指标0

31、. 05% ,开口风洞试验段的湍流度值在除60m/ s外的大部分风速段都满足设计指标0. 2% ,常规低速风洞各风速段流场湍流度低于0. 2% 。( 3)采用EM D方法获得的HHT时频谱可以分析脉动速度信号的时频特征。分析发现开口风洞试验段的HHT时频谱有明显较大的低频脉动信号,可以为风洞调试提供分析依据。参考文献: 1何克敏,白存儒,郭渠渝,等.较低湍流度范围湍流度对风洞实验结果的影响 J .流体力学实验与测量, 1997, 11( 3) : 11-17.He K M , Bai C R, Guo Q Y, et al. The effect of turbulence onwind tu

32、nnel results in the range of low turbulence J . Experi-ments and M easurements in Fluid M echanics, 1997, 11 ( 3) :11- 17. 2李峰,白存儒,周伟,等.湍流度对飞行器模型大迎角气动特性影响的初步研究 J .实验流体力学, 2006, 20( 3) : 45- 52.Li F, Bai C R, Zhou W , et al. Primal research of the effect offlow turbulence on aerodynamic characteristi

33、cs of a aircraftmodel at high angles of attack J . Journal of Experiments inFluid M echanics, 2006, 20( 3) : 45- 52. 3米建春,冯宝平.平面射流沿轴线的特征尺度及其对测量信号过滤程度的依赖 J .物理学报, 2010, 59( 7) : 4748- 4755.M i J C, Feng B P. Centerline characteristic scales of a turbulentplane jet and their dependence on filtration o

34、f measured signals J . Acta Physica Sinica, 2010, 59( 7) : 4748- 4755. 4侯志勇,王连泽,周建和,等.低(变)湍流度风洞设计再探讨 J .实验流体力学, 2011, 25( 1) : 92- 96.Hou Z Y, W ang L Z, Zhou J H, et al. The further research onthe design of low ( varying) turbulence wind tunnel J . Journalof Experiments in Fluid M echanics, 2011, 2

35、5( 1) : 92- 96. 5 Harald Q, J rgen Q, Christian B. Hot- W ire measurements incryogenic environment R . AIAA- 2011- 880. 6 Saric W S, Takagi S, M ousseux M . The ASU unsteady windtunnel and fundamental requirements for freestream turbulencemeasurements R . AIAA- 1988- 0053, 1988.87实 验 流 体 力 学 ( 2016)

36、第30卷 万方数据 7 M ueller T J, Scharpf D F, Batill S M , et al. The design of asubsonic low- noise, low- turbulence wind tunnel for acousticmeasurements R . AIAA- 1992- 3883, 1992. 8 Neuhart Dan H, M cGinley Catherine B. Free- stream turbu-lence intensity in the Langley 14- by 22- foot subsonic tunnel R

37、. NASA/ TP- 2004- 213247, 2004. 9 Premi Amandeep, M aughmer M ark D, Brophy Christopher.Flow- quality measurements and qualification of thePennsylvaniaState University low- Speed, low- Turbulence wind tunnel R .AIAA- 2012- 1214, 2012. 10 Stainback P Calvin, Owen F Kevin. Dynamic flow qualitymeasurem

38、ents in the Langley low turbulence pressure tunnel R . AIAA- 1984- 0621, 1984. 11 Kohama Y, Kobayashi R, Ito H. Tohoku university low- tur-bulence wind tunnel R . AIAA- 1992- 3913, 1992. 12李士心,刘鲁源,舒玮.子波去噪技术在湍流信号处理中的应用 J .实验力学, 2001, 16( 4) : 433- 437.Li S X, Liu L Y, Shu W . Turbulent signal process

39、ing based onwavelet transform de- noising technique J . Journal of Experi-mental M echanics, 2001, 16( 4) : 433- 437. 13朱博,汤更生.声学风洞流场低湍流度及频谱测量研究 J .实验流体力学, 2015, 29( 4) : 58- 64.Zhu B, Tang G S. Low turbulence intensity and spectrummeasurement research in aeroacoustic wind tunnel J . Journalof Exper

40、iments in Fluid M echanics, 2015, 29( 4) : 58- 64. 14 Huang N E, Shen Z, Long S R. The empirical mode decompo-sition and hilbert spectrum for nonlinear and nonstationary timeseries analysis J . Proceedings of the Royal Society of LondonSERIESA, 1998, 454: 903- 995. 15戴吾蛟,丁晓利,朱建军,等.基于经验模式分解的滤波去噪法及其在G

41、PS多路径效应中的应用 J .测绘学报, 2006, 35( 11) : 321- 327.Dai W J, Ding X L, Zhu J J, et al. EM D filter method and itsapplication in GPS multipath J . Acta Geodaetica et Carto-graphica Sinica, 2006, 35( 11) : 321- 327. 16楼文娟,卢旦,杨毅,等.开孔建筑屋盖风振响应中的气动阻尼识别 J .空气动力学学报, 2007, 25( 4) : 419- 424.Lou W J, Lu D, Yang Y,

42、et al. Identification of aerodynamicdamping for roof wind induced response of opening building J .Acta Aerodynamica Sinica, 2007, 25( 4) : 419- 424. 17 Flandrin P, Rilling G. Empirical mode decomposition as filterbank J . IEEE Signal Process Letters, 2003, 11: 112- 114. 18朱博,王元兴,余永生.风洞模型-支撑系统涡激振动测量与

43、分析 J .实验流体力学, 2014, 28( 5) : 59- 64.Zhu B, W ang Y X, Yu Y S. Vortex- induced vibration measure-ment and analysis of model- sting system in wind tunnel J .Journal of Experiments in Fluid M echanics, 2014, 28( 5) : 59-64.作者简介:朱 博( 1973- ) ,男,广西百色人,高级工程师。研究方向:风洞测试、声学测试和热线技术研究。通信地址:四川省绵阳市二环路南段6号14信箱401分箱( 621000) 。 E-mail: bobjou 139. com(编辑:杨 娟)97 第5期 朱 博等:一种基于EM D的低湍流度信号处理分析方法 万方数据