CN110844115A - 一种桨涡干扰噪声和桨叶挥舞载荷数据有效性判别方法 - Google Patents
一种桨涡干扰噪声和桨叶挥舞载荷数据有效性判别方法 Download PDFInfo
- Publication number
- CN110844115A CN110844115A CN201910993739.XA CN201910993739A CN110844115A CN 110844115 A CN110844115 A CN 110844115A CN 201910993739 A CN201910993739 A CN 201910993739A CN 110844115 A CN110844115 A CN 110844115A
- Authority
- CN
- China
- Prior art keywords
- blade
- interference noise
- vortex interference
- load
- propeller
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
- B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
- B64F5/60—Testing or inspecting aircraft components or systems
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Transportation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明属直升机旋翼试验领域,涉及一种桨涡干扰噪声和桨叶挥舞载荷数据有效性判别方法。通过对同步采集得到的桨涡干扰噪声数据和桨叶挥舞载荷数据进行整周期平均,提取出桨涡干扰噪声和桨叶挥舞载荷最大峰值处的相位,换算成距离差,然后和理论计算的距离差进行对比,从而判别出桨涡干扰噪声数据和桨叶挥舞载荷数据是否有效,能够避免某些试验状态下因实测数据和理论计算有所偏差而认为数据失效的误判。
Description
技术领域
本发明属直升机模型旋翼试验技术领域,涉及一种桨涡干扰噪声和桨叶挥舞载荷数据有效性判别方法。
背景技术
旋翼桨涡干扰噪声是直升机旋翼产生的一种典型噪声,尤其在小速度前飞和斜下降状态,桨涡干扰噪声的能量占到直升机噪声总能量的80%。因此它是目前直升机降低噪声的主要研究领域。桨叶载荷反应的是桨叶承受多大的交变力,是检测并判断桨叶是否工作正常的一种常规手段。桨涡干扰是桨叶切割前一片桨叶脱落桨尖涡产生的气动扰动,该扰动同时也会对桨叶载荷尤其是挥舞载荷产生激励作用,两者之间存在内在的联系。
目前国内外对旋翼桨涡干扰噪声数据的有效判断都是仅从其自身数据的规律是否和计算相符来进行判断,桨叶挥舞载荷数据的判断也是如此,并没有将桨涡干扰噪声数据和桨叶载荷数据关联起来进行分析判断,会造成某些试验状态下因实测数据和理论计算有所偏差而认为数据失效的误判。
发明内容
本发明的目的:提供一种桨涡干扰噪声和桨叶挥舞载荷数据有效性判别方法,通过对同步采集得到的桨涡干扰噪声数据和桨叶挥舞载荷数据进行整周期平均,提取出桨涡干扰噪声和桨叶挥舞载荷最大峰值处的相位,换算成空间距离差,然后和理论计算的空间距离差进行对比,从而判别出桨涡干扰噪声数据和桨叶挥舞载荷数据是否有效。
本发明的技术方案:
第一方面,提供一种桨涡干扰噪声和桨叶挥舞载荷数据有效性判别方法,包括:
同步采集转速、桨涡干扰噪声和桨叶挥舞载荷;
根据转速对桨涡干扰噪声和桨叶挥舞载荷分别进行整周期平均得到整周期平均的桨涡干扰噪声和整周期平均的桨叶挥舞载荷;
分别提取整周期平均的桨涡干扰噪声的峰值相位和整周期平均的桨叶挥舞载荷的峰值相位得到桨涡干扰噪声的峰值相位和桨叶挥舞载荷的峰值相位;
利用桨涡干扰噪声的峰值相位和桨叶挥舞载荷挥舞峰值相位计算桨涡干扰噪声的测点试验距离差;
利用桨涡干扰噪声的测点试验距离差计算桨涡干扰噪声的测点与桨叶挥舞载荷的测点试验距离差得到第一试验距离差;
计算桨涡干扰噪声的测点与桨叶挥舞载荷测点的理论时间差;
根据桨涡干扰噪声的测点与桨叶挥舞载荷测点的理论时间差计算桨涡干扰噪声的测点与桨叶挥舞载荷测点的理论距离差得到第一理论距离差;
根据第一试验距离差计算第一试验距离差的标准误差;
根据第一理论距离差计算第一理论距离差的标准误差;
在第一试验距离差的标准误差与第一理论距离差的标准误差的差值的绝对值小于第一阈值的情况下,将桨涡干扰噪声和桨叶挥舞载荷数据判定为有效,否则将桨涡干扰噪声和桨叶挥舞载荷数据判定为无效。
可选地,同步采集转速、桨涡干扰噪声和桨叶挥舞载荷,具体包括:采集转速;根据转速设置采样率;按照所述采样率同步采集不少于100转的转速、桨涡干扰噪声和桨叶挥舞载荷。
可选地,根据转速对桨涡干扰噪声和桨叶挥舞载荷分别进行整周期平均得到整周期平均的桨涡干扰噪声和整周期平均的桨叶挥舞载荷,具体包括:
以转速的任意一个脉冲上升沿为起点,截取数据点数为360000的桨涡干扰噪声数据和桨叶挥舞载荷数据;
将桨涡干扰噪声数据和桨叶挥舞载荷数据依次按照数据点数为3600等间隔分割成100组数据;
将分割后的100组桨涡干扰噪声中的每组桨涡干扰噪声的第1个点、第2个点到第3600个点累加并除以100,得到整周期平均的桨涡干扰噪声;
将分割后的100组桨叶挥舞载荷中的每组桨叶挥舞载荷的第1个点、第2个点到第3600个点累加并除以100,得到整周期平均的桨叶挥舞载荷。
可选地,分别提取整周期平均的桨涡干扰噪声的峰值相位和整周期平均的桨叶挥舞载荷的峰值相位得到桨涡干扰噪声的峰值相位和桨叶挥舞载荷的峰值相位,具体包括:
将整周期平均的桨涡干扰噪声的第三个波峰的峰值相位确定为桨涡干扰噪声的峰值相位;
计算整周期平均的桨涡干扰噪声的最小峰值相位,根据整周期平均的桨涡干扰噪声的最小峰值相位确定整周期平均的桨叶挥舞载荷的峰值参考数据位置点;
根据整周期平均的桨叶挥舞载荷的峰值参考数据位置点确定桨叶挥舞载荷的峰值相位。
可选地,利用桨涡干扰噪声的峰值相位和桨叶挥舞载荷挥舞峰值相位计算桨涡干扰噪声的测点试验距离差,具体包括:
以桨叶挥舞载荷挥舞峰值相位为基准,计算桨涡干扰噪声和桨叶挥舞载荷的相位差;
根据桨涡干扰噪声和桨叶挥舞载荷的相位差,计算桨涡干扰噪声的测点试验距离差。
可选地,利用桨涡干扰噪声的测点试验距离差计算桨涡干扰噪声的测点与桨叶挥舞载荷的测点试验距离差得到第一试验距离差,具体包括:将桨涡干扰噪声的测点试验距离差乘以声音速度并除以采样率得到第一试验距离差。
可选地,计算桨涡干扰噪声的测点与桨叶挥舞载荷测点的理论时间差,具体包括:
将理论的旋翼桨尖到桨叶挥舞载荷测点声传播的时间差计算为第一理论时间差;
计算旋翼桨尖到桨涡干扰噪声测点的理论距离差,根据旋翼桨尖到桨涡干扰噪声测点的理论距离差得到第二理论时间差;
将第二理论时间差减去第一理论时间差得到桨涡干扰噪声的测点与桨叶挥舞载荷测点的理论时间差。
可选地,根据桨涡干扰噪声的测点与桨叶挥舞载荷测点的理论时间差计算桨涡干扰噪声的测点与桨叶挥舞载荷测点的理论距离差得到第一理论距离差,具体包括:将桨涡干扰噪声的测点与桨叶挥舞载荷测点的理论时间差乘以声音速度得到第一理论距离差。
可选地,第一阈值为0.1。
本发明的有益效果:
本发明提供的桨涡干扰噪声和桨叶挥舞载荷数据有效性判别方法,通过对同步采集得到的桨涡干扰噪声数据和桨叶挥舞载荷数据进行整周期平均,提取出桨涡干扰噪声和桨叶挥舞载荷最大峰值处的相位,换算成距离差,然后和理论计算的距离差进行对比,从而判别出桨涡干扰噪声数据和桨叶挥舞载荷数据是否有效,可避免某些试验状态下因实测数据和理论计算有所偏差而认为数据失效的误判。
附图说明
图1为根据本发明实施例的桨涡干扰噪声和桨叶挥舞载荷数据有效性判别方法的流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本方法通过对同步采集得到的桨涡干扰噪声数据和桨叶挥舞载荷数据进行整周期平均,提取出桨涡干扰噪声和桨叶挥舞载荷最大峰值处的相位,换算成距离差,然后和理论计算的距离差进行对比,从而判别出桨涡干扰噪声数据和桨叶挥舞载荷数据是否有效。
(1)转速和桨涡干扰噪声、桨叶挥舞载荷同步采集
转速和桨涡干扰噪声、桨叶挥舞载荷同步采集步骤如下:
第1步,将涡发生器安装在和桨盘平面同高的位置,涡发生器出涡口正对桨尖位置,距离为100mm;
第2步,已涡发生器位置为0度方位,以距离桨榖中心为圆心,以1.5倍桨叶长度Lr(Lr的单位为米)为半径,按桨叶旋转方向顺时针旋转45度,在该位置竖直布置一根传声器枝干,然后从桨盘平面高度6米开始,在传声器枝干上以0.5米等间隔布置6个点传声器,与桨盘平面相同高度的为1号传声器,低于桨盘平面0.5米的为2号传声器,依次类推最低的点为6号传声器;
第3步,将转速传感器安装在涡发生器的同一个方位,在0.8Lr剖面位置布有挥舞载荷应变片的桨叶(以下假定为第3号桨叶)对面的桨叶(以下假定为第1号桨叶)上贴上反光片;
第4步,根据转速Np(Np的单位为转/秒),设置采样率Sm=Np*3600,然后待试验状态稳定后同步采集不少于Tn=100/Np秒的转速和桨涡干扰噪声、桨叶挥舞载荷数据,既至少采集100圈的数据,同时要求试验时操纵总距在±0.1度范围以内,涡发生器的出涡口速度不小于11米/秒。
(2)转速桨涡干扰噪声、桨叶挥舞载荷试验数据整周期平均
整周期平均步骤如下:
第1步,以转速任意一个脉冲上升沿为起点,截取数据点长度为360000的桨涡干扰噪声数据和桨叶挥舞载荷数据;
第2步,然后将上述数据依次按照数据点长度为3600等间隔分割成100圈数据,每圈3600个点;
第3步,将分割后的100圈数据每圈的第1个点、第2个点到第3600个点分别提取出来累加后除100,从而得到1圈(即3600个点)平均数据,及整周期平均数据。
(3)桨涡干扰噪声、桨叶挥舞载荷试验数据峰值相位提取
峰值相位提取步骤如下:
第1步,从1号传声器的整周期平均数据中找到第3个桨涡干扰噪声波峰(对应第3号桨叶)的峰值点的位置L1v0,则对应的相位为:P1v0=L1v0/10;
第2步,重复第1步,依次从第2个到第6个传声器整周期平均数据中找到第3个桨涡干扰噪声波峰的峰值点的位置L2v0、L3v0、L4v0、L5v0、L6v0,则对应的相位分别为:P2v0=L2v0/10、P3v0=L3v0/10、P4v0=L4v0/10、P5v0=L5v0/10、P6v0=L6v0/10;
第3步,从P1vn、P2vn、P3vn、P4vn、P5vn、P6vnn中找到值最小的相位Pmin,从而计算得到在1圈平均数据中的桨叶挥舞载荷波峰参考数据位置点Lh0=Pmin×10;
第4步,以第3步计算得到的参考数据位置点Lh0为起点,从桨叶挥舞载荷整周期平均数据向前查找到出现第1个波峰的峰值点的位置Lh1,则对应的相位为:Ph1=Lh1/10。
(4)桨涡干扰噪声、桨叶挥舞载荷测点试验距离差计算
桨涡干扰噪声、桨叶挥舞载荷测点试验距离差计算步骤如下:
第1步,以桨叶挥舞载荷挥舞峰值相位Ph1为基准,分别计算得到1号到6号传声器的相位差P1v1=P1v0-Ph1、P2v1=P2v0-Ph1、P3v1=P3v0-Ph1、P4v1=P4v0-Ph1、P5v1=P5v0-Ph1、P6v1=P6v0-Ph1;
第2步,计算得到1号到6号传声器的数据点差分别是:L1v1=P1v1×10、L2v1=P2v1×10、L3v1=P3v1×10、L4v1=P4v1×10、L5v1=P5v1×10、L6v1=P6v1×10;
第3步,计算得到1号到6号传声器到桨叶挥舞载荷测点之间的第一试验距离差分别是:D1t=L1v1×340/Sm、D2t=L2v1×340/Sm、D3t=L3v1×340/Sm、D4t=L4v1×340/Sm、D5t=L5v1×340/Sm、D6t=L6v1×340/Sm。
(5)桨涡干扰噪声、桨叶挥舞载荷测点理论距离差计算
桨涡干扰噪声、桨叶挥舞载荷测点理论距离差计算
第1步,计算理论的旋翼桨尖到桨叶挥舞载荷测点声传播的第一理论时间差Tch=0.2×R/340;
第2步,计算旋翼桨尖到1号传声器测点处的第一理论距离差为L1c0=0.707×1.5×R,然后计算得到旋翼桨尖到1号传声器测点处的第二理论时间差T1c0=L1c0/340;
第3步,计算旋翼桨尖到2号传声器测点处的第一理论距离差为L2c0=SQRT[(0.707×1.5×R)^2+0.5×0.5],SQRT为开平方根函数,然后计算得到旋翼桨尖到2号传声器测点处的第二理论时间差T2c0=L2c0/340;
第4步,计算旋翼桨尖到3号传声器测点处的第一理论距离差为L3c0=SQRT[(0.707×1.5×R)^2+1.0×1.0],SQRT为开平方根函数,然后计算得到旋翼桨尖到3号传声器测点处的第二理论时间差T3c0=L3c0/340;
第5步,计算旋翼桨尖到4号传声器测点处的第一理论距离差为L4c0=SQRT[(0.707×1.5×R)^2+1.5×1.5],SQRT为开平方根函数,然后计算得到旋翼桨尖到4号传声器测点处的第二理论时间差T4c0=L4c0/340;
第6步,计算旋翼桨尖到5号传声器测点处的第一理论距离差为L5c0=SQRT[(0.707×1.5×R)^2+2.0×2.0],SQRT为开平方根函数,然后计算得到旋翼桨尖到5号传声器测点处的第二理论时间差T5c0=L5c0/340;
第7步,计算旋翼桨尖到6号传声器测点处的第一理论距离差为L6c0=SQRT[(0.707×1.5×R)^2+2.5×2.5],SQRT为开平方根函数,然后计算得到旋翼桨尖到6号传声器测点处的第二理论时间差T6c0=L6c0/340;
第8步,计算得到1号到6号传声器到桨叶挥舞载荷测点之间理论时间差分别是:T1c1=T1c0-Tch、T2c1=T2c0-Tch、T3c1=T3c0-Tch、T4c1=T4c0-Tch、T5c1=T5c0-Tch、T6c1=T6c0-Tch。
第9步,计算得到1号到6号传声器到桨叶挥舞载荷测点之间理论距离差分别是:D1c=T1c1×340、D2c=T2c1×340、D3c=T3c1×340、D4c=T4c1×340、D5c=T5c1×340、D6c=T6c1×340。
(6)试验距离差与理论距离差对比判别数据有效性
试验距离差与理论距离差对比判别数据有效性步骤如下:
第1步,计算得到1号到6号传声器的第一试验距离差和第二理论距离差的误差分别是:D1=ABS(D1c-D1t)、D2=ABS(D2c-D2t)、D3=ABS(D3c-D3t)、D4=ABS(D4c-D4t)、D5=ABS(D5c-D5t)、D6=ABS(D6c-D6t),其中ABS为去绝对值函数。
第2步,将计算得到的误差D1、D2、D3、D4、D5、D6和设定的第一阈值Ds=0.1进行对比,如果D1、D2、D3、D4、D5、D6全部小于该标准误差,则认为试验采集到的桨涡干扰噪声和桨叶挥舞载荷数据是有效的。
示例如下:
(1)转速和桨涡干扰噪声、桨叶挥舞载荷同步采集
第1步,将涡发生器安装在和桨盘平面同高的位置,涡发生器出涡口正对桨尖位置,距离为100mm;
第2步,已涡发生器位置为0度方位,以距离桨榖中心为圆心,以1.5倍桨叶长度Lr(Lr的单位为米)为半径,按桨叶旋转方向顺时针旋转45度,在该位置竖直布置一根传声器枝干,然后从桨盘平面高度6米开始,在传声器枝干上以0.5米等间隔布置6个点传声器,与桨盘平面相同高度的为1号传声器,低于桨盘平面0.5米的为2号传声器,依次类推最低的点为6号传声器;
第3步,将转速传感器安装在涡发生器的同一个方位,在0.8Lr剖面位置布有挥舞载荷应变片的桨叶(以下假定为第3号桨叶)对面的桨叶(以下假定为第1号桨叶)上贴上反光片;
第4步,根据转速Np=10(Np的单位为转/秒),设置采样率Sm=Np*3600=36000,然后待试验状态稳定后同步采集不少于Tn=100/Np=10秒的转速和桨涡干扰噪声、桨叶挥舞载荷数据,既至少采集100圈的数据,同时要求试验时操纵总距为0度,涡发生器的出涡口速度为15米/秒。
(2)转速桨涡干扰噪声、桨叶挥舞载荷试验数据整周期平均
整周期平均步骤如下:
第1步,以转速任意一个脉冲上升沿为起点,截取数据点长度为360000的桨涡干扰噪声数据和桨叶挥舞载荷数据;
第2步,然后将上述数据依次按照数据点长度为3600等间隔分割成100圈数据,每圈3600个点;
第3步,将分割后的100圈数据每圈的第1个点、第2个点到第3600个点分别提取出来累加后除100,从而得到1圈(即3600个点)平均数据,既整周期平均数据。
(3)桨涡干扰噪声、桨叶挥舞载荷试验数据峰值相位提取
峰值相位提取步骤如下:
第1步,从1号传声器的整周期平均数据中找到第3个桨涡干扰噪声(对应第3号桨叶)波峰的峰值点的位置L1v0=2364,则对应的相位为:P1v0=L1v0/10=236.4;
第2步,重复第1步,依次从第2个到第6个传声器整周期平均数据中找到第3个桨涡干扰噪声波峰的峰值点的位置L2v0=2369、L3v0=2390、L4v0=2413、L5v0=2450、L6v0=2485,则对应的相位分别为:P2v0=L2v0/10=236.9、P3v0=L3v0/10=239.0、P4v0=L4v0/10=241.3、P5v0=L5v0/10=245.0、P6v0=L6v0/10=248.5;
第3步,从P1vn、P2vn、P3vn、P4vn、P5vn、P6vnn中找到值最小的相位Pmin=236.4,从而计算得到在1圈平均数据中的桨叶挥舞载荷波峰参考数据位置点Lh0=Pmin×10=2364;
第4步,以第3步计算得到的参考数据位置点Lh0为起点,从桨叶挥舞载荷整周期平均数据向前查找到出现第1个波峰的峰值点的位置Lh1=2182,则对应的相位为:Ph1=Lh1/10=218.2。
(4)桨涡干扰噪声、桨叶挥舞载荷测点试验距离差计算
桨涡干扰噪声、桨叶挥舞载荷测点试验距离差计算步骤如下:
第1步,以桨叶挥舞载荷挥舞峰值相位Ph1为基准,分别计算得到1号到6号传声器的相位差P1v1=P1v0-Ph1=18.2、P2v1=P2v0-Ph1=18.7、P3v1=P3v0-Ph1=20.8、P4v1=P4v0-Ph1=23.1、P5v1=P5v0-Ph1=26.8、P6v1=P6v0-Ph1=30.3;
第2步,计算得到1号到6号传声器的数据点差分别是:L1v1=P1v1×10=182、L2v1=P2v1×10=187、L3v1=P3v1×10=208、L4v1=P4v1×10=231、L5v1=P5v1×10=268、L6v1=P6v1×10=303;
第3步,计算得到1号到6号传声器到桨叶挥舞载荷测点之间的试验距离差分别是:D1t=L1v1×340/Sm=1.7189、D2t=L2v1×340/Sm=1.7661、D3t=L3v1×340/Sm=1.9644、D4t=L4v1×340/Sm=2.1817、D5t=L5v1×340/Sm=2.5311、D6t=L6v1×340/Sm=2.8617。
(5)桨涡干扰噪声、桨叶挥舞载荷测点理论距离差计算
桨涡干扰噪声、桨叶挥舞载荷测点理论距离差计算
第1步,计算理论的旋翼桨尖到桨叶挥舞载荷测点声传播的第一理论时间差Tch=0.2×R/340=0.2×2/340=0.001176,其中R=2;
第2步,计算旋翼桨尖到1号传声器测点处的第一理论距离差为L1c0=0.707×1.5×R=2.1210,然后计算得到第二理论时间差T1c0=L1c0/340=0.006238;
第3步,计算旋翼桨尖到2号传声器测点处的第一理论距离差为L2c0=SQRT[(0.707×1.5×R)^2+0.5×0.5]=2.1791,SQRT为开平方根函数,然后计算得到第二理论时间差T2c0=L2c0/340=0.006409;
第4步,计算旋翼桨尖到3号传声器测点处的第一理论距离差为L3c0=SQRT[(0.707×1.5×R)^2+1.0×1.0]=2.3449,SQRT为开平方根函数,然后计算得到第二理论时间差T3c0=L3c0/340=0.006897;
第5步,计算旋翼桨尖到4号传声器测点处的第一理论距离差为L4c0=SQRT[(0.707×1.5×R)^2+1.5×1.5]=2.5978,SQRT为开平方根函数,然后计算得到第二理论时间差T4c0=L4c0/340=0.007641;
第6步,计算旋翼桨尖到5号传声器测点处的第一理论距离差为L5c0=SQRT[(0.707×1.5×R)^2+2.0×2.0]=2.9152,SQRT为开平方根函数,然后计算得到第二理论时间差T5c0=L5c0/340=0.008574;
第7步,计算旋翼桨尖到6号传声器测点处的第一理论距离差为L6c0=SQRT[(0.707×1.5×R)^2+2.5×2.5]=3.2785,SQRT为开平方根函数,然后计算得到第二理论时间差T6c0=L6c0/340=0.009643;
第8步,计算得到1号到6号传声器到桨叶挥舞载荷测点之间理论时间差分别是:T1c1=T1c0-Tch=0.005062、T2c1=T2c0-Tch=0.005233、T3c1=T3c0-Tch=0.005720、T4c1=T4c0-Tch=0.006464、T5c1=T5c0-Tch=0。007398、T6c1=T6c0-Tch=0.008466。
第9步,计算得到1号到6号传声器到桨叶挥舞载荷测点之间理论距离差分别是:D1c=T1c1×340=1.7210、D2c=T2c1×340=1.7791、D3c=T3c1×340=1.9449、D4c=T4c1×340=2.1978、D5c=T5c1×340=2.5152、D6c=T6c1×340=2.8785。
(6)试验距离差与理论距离差对比判别数据有效性
试验距离差与理论距离差对比判别数据有效性步骤如下:
第1步,计算得到1号到6号传声器的试验距离差和理论距离差的误差分别是:D1=ABS(D1c-D1t)=0.0021、D2=ABS(D2c-D2t)=0.0130、D3=ABS(D3c-D3t)=0.0195、D4=ABS(D4c-D4t)=0.0161、D5=ABS(D5c-D5t)=0.0159、D6=ABS(D6c-D6t)=0.0168,其中ABS为去绝对值函数。
第2步,将计算得到的误差D1、D2、D3、D4、D5、D6和设定的第一阈值Ds=0.025进行对比,因为D1、D2、D3、D4、D5、D6全部小于该标准误差,所以认为试验采集到的桨涡干扰噪声和桨叶挥舞载荷数据是有效的。
以上所述,仅为本发明的具体实施例,对本发明进行详细描述,未详尽部分为常规技术。但本发明的保护范围不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。本发明的保护范围应以所述权利要求的保护范围为准。
Claims (9)
1.一种桨涡干扰噪声和桨叶挥舞载荷数据有效性判别方法,其特征在于,包括:
同步采集转速、桨涡干扰噪声和桨叶挥舞载荷;
根据转速对桨涡干扰噪声和桨叶挥舞载荷分别进行整周期平均得到整周期平均的桨涡干扰噪声和整周期平均的桨叶挥舞载荷;
分别提取整周期平均的桨涡干扰噪声的峰值相位和整周期平均的桨叶挥舞载荷的峰值相位得到桨涡干扰噪声的峰值相位和桨叶挥舞载荷的峰值相位;
利用桨涡干扰噪声的峰值相位和桨叶挥舞载荷挥舞峰值相位计算桨涡干扰噪声的测点试验距离差;
利用桨涡干扰噪声的测点试验距离差计算桨涡干扰噪声的测点与桨叶挥舞载荷的测点试验距离差得到第一试验距离差;
计算桨涡干扰噪声的测点与桨叶挥舞载荷测点的理论时间差;
根据桨涡干扰噪声的测点与桨叶挥舞载荷测点的理论时间差计算桨涡干扰噪声的测点与桨叶挥舞载荷测点的理论距离差得到第一理论距离差;
根据第一试验距离差计算第一试验距离差的标准误差;
根据第一理论距离差计算第一理论距离差的标准误差;
在第一试验距离差的标准误差与第一理论距离差的标准误差的差值的绝对值小于第一阈值的情况下,将桨涡干扰噪声和桨叶挥舞载荷数据判定为有效,否则将桨涡干扰噪声和桨叶挥舞载荷数据判定为无效。
2.根据权利要求1所述的方法,其特征在于,同步采集转速、桨涡干扰噪声和桨叶挥舞载荷,具体包括:采集转速;根据转速设置采样率;按照所述采样率同步采集不少于100转的转速、桨涡干扰噪声和桨叶挥舞载荷。
3.根据权利要求1所述的方法,其特征在于,根据转速对桨涡干扰噪声和桨叶挥舞载荷分别进行整周期平均得到整周期平均的桨涡干扰噪声和整周期平均的桨叶挥舞载荷,具体包括:
以转速的任意一个脉冲上升沿为起点,截取数据点数为360000的桨涡干扰噪声数据和桨叶挥舞载荷数据;
将桨涡干扰噪声数据和桨叶挥舞载荷数据依次按照数据点数为3600等间隔分割成100组数据;
将分割后的100组桨涡干扰噪声中的每组桨涡干扰噪声的第1个点、第2个点到第3600个点累加并除以100,得到整周期平均的桨涡干扰噪声;
将分割后的100组桨叶挥舞载荷中的每组桨叶挥舞载荷的第1个点、第2个点到第3600个点累加并除以100,得到整周期平均的桨叶挥舞载荷。
4.根据权利要求1所述的方法,其特征在于,分别提取整周期平均的桨涡干扰噪声的峰值相位和整周期平均的桨叶挥舞载荷的峰值相位得到桨涡干扰噪声的峰值相位和桨叶挥舞载荷的峰值相位,具体包括:
将整周期平均的桨涡干扰噪声的第三个波峰的峰值相位确定为桨涡干扰噪声的峰值相位;
计算整周期平均的桨涡干扰噪声的最小峰值相位,根据整周期平均的桨涡干扰噪声的最小峰值相位确定整周期平均的桨叶挥舞载荷的峰值参考数据位置点;
根据整周期平均的桨叶挥舞载荷的峰值参考数据位置点确定桨叶挥舞载荷的峰值相位。
5.根据权利要求1所述的方法,其特征在于,利用桨涡干扰噪声的峰值相位和桨叶挥舞载荷挥舞峰值相位计算桨涡干扰噪声的测点试验距离差,具体包括:
以桨叶挥舞载荷挥舞峰值相位为基准,计算桨涡干扰噪声和桨叶挥舞载荷的相位差;
根据桨涡干扰噪声和桨叶挥舞载荷的相位差,计算桨涡干扰噪声的测点试验距离差。
6.根据权利要求1所述的方法,其特征在于,利用桨涡干扰噪声的测点试验距离差计算桨涡干扰噪声的测点与桨叶挥舞载荷的测点试验距离差得到第一试验距离差,具体包括:将桨涡干扰噪声的测点试验距离差乘以声音速度并除以采样率得到第一试验距离差。
7.根据权利要求1所述的方法,其特征在于,计算桨涡干扰噪声的测点与桨叶挥舞载荷测点的理论时间差,具体包括:
将理论的旋翼桨尖到桨叶挥舞载荷测点声传播的时间差计算为第一理论时间差;
计算旋翼桨尖到桨涡干扰噪声测点的理论距离差,根据旋翼桨尖到桨涡干扰噪声测点的理论距离差得到第二理论时间差;
将第二理论时间差减去第一理论时间差得到桨涡干扰噪声的测点与桨叶挥舞载荷测点的理论时间差。
8.根据权利要求1所述的方法,其特征在于,根据桨涡干扰噪声的测点与桨叶挥舞载荷测点的理论时间差计算桨涡干扰噪声的测点与桨叶挥舞载荷测点的理论距离差得到第一理论距离差,具体包括:将桨涡干扰噪声的测点与桨叶挥舞载荷测点的理论时间差乘以声音速度得到第一理论距离差。
9.根据权利要求1所述的方法,其特征在于,第一阈值为0.1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910993739.XA CN110844115B (zh) | 2019-10-18 | 2019-10-18 | 一种桨涡干扰噪声和桨叶挥舞载荷数据有效性判别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910993739.XA CN110844115B (zh) | 2019-10-18 | 2019-10-18 | 一种桨涡干扰噪声和桨叶挥舞载荷数据有效性判别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110844115A true CN110844115A (zh) | 2020-02-28 |
CN110844115B CN110844115B (zh) | 2022-04-12 |
Family
ID=69597557
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910993739.XA Active CN110844115B (zh) | 2019-10-18 | 2019-10-18 | 一种桨涡干扰噪声和桨叶挥舞载荷数据有效性判别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110844115B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112487360A (zh) * | 2020-11-20 | 2021-03-12 | 中国直升机设计研究所 | 一种基于统计分析的旋翼桨涡干扰噪声整周期平均方法 |
CN112504426A (zh) * | 2020-11-20 | 2021-03-16 | 中国直升机设计研究所 | 一种基于峰值搜索的旋翼桨涡干扰噪声整周期平均方法 |
CN114235360A (zh) * | 2021-11-19 | 2022-03-25 | 中国直升机设计研究所 | 一种直升机桨叶实测载荷相位关系的分析方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1103901A (en) * | 1965-07-30 | 1968-02-21 | Copperfield Corp | Helicopter rotor |
US4608650A (en) * | 1983-08-09 | 1986-08-26 | Becton Dickinson And Company | Imbalance measuring system and method |
JPH03258198A (ja) * | 1990-03-08 | 1991-11-18 | Takaoka Electric Mfg Co Ltd | 周期性設備音の雑音除去方式 |
FR2782307A1 (fr) * | 1998-08-17 | 2000-02-18 | Onera (Off Nat Aerospatiale) | Procede pour la reduction du bruit d'interaction pales-tourbillons engendre par une voilure tournante |
US20110280410A1 (en) * | 2009-01-21 | 2011-11-17 | Panasonic Corporation | Active noise control apparatus |
CN105092204A (zh) * | 2015-07-24 | 2015-11-25 | 南京航空航天大学 | 一种用于直升机旋翼桨-涡干扰噪声实验的装置 |
WO2017031581A1 (en) * | 2015-08-24 | 2017-03-02 | UNIVERSITé LAVAL | System and method for detecting spikes in noisy signals |
CN108100300A (zh) * | 2017-12-03 | 2018-06-01 | 中国直升机设计研究所 | 一种用于全消声室桨-涡干扰试验的涡发生装置 |
CN108216617A (zh) * | 2017-12-29 | 2018-06-29 | 厦门大学 | 一种抑制直升机桨-涡干扰噪声的方法 |
CN108910079A (zh) * | 2018-06-26 | 2018-11-30 | 中国直升机设计研究所 | 一种旋翼试验台动平衡调整相位识别方法 |
-
2019
- 2019-10-18 CN CN201910993739.XA patent/CN110844115B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1103901A (en) * | 1965-07-30 | 1968-02-21 | Copperfield Corp | Helicopter rotor |
US4608650A (en) * | 1983-08-09 | 1986-08-26 | Becton Dickinson And Company | Imbalance measuring system and method |
JPH03258198A (ja) * | 1990-03-08 | 1991-11-18 | Takaoka Electric Mfg Co Ltd | 周期性設備音の雑音除去方式 |
FR2782307A1 (fr) * | 1998-08-17 | 2000-02-18 | Onera (Off Nat Aerospatiale) | Procede pour la reduction du bruit d'interaction pales-tourbillons engendre par une voilure tournante |
US20110280410A1 (en) * | 2009-01-21 | 2011-11-17 | Panasonic Corporation | Active noise control apparatus |
CN105092204A (zh) * | 2015-07-24 | 2015-11-25 | 南京航空航天大学 | 一种用于直升机旋翼桨-涡干扰噪声实验的装置 |
WO2017031581A1 (en) * | 2015-08-24 | 2017-03-02 | UNIVERSITé LAVAL | System and method for detecting spikes in noisy signals |
CN108100300A (zh) * | 2017-12-03 | 2018-06-01 | 中国直升机设计研究所 | 一种用于全消声室桨-涡干扰试验的涡发生装置 |
CN108216617A (zh) * | 2017-12-29 | 2018-06-29 | 厦门大学 | 一种抑制直升机桨-涡干扰噪声的方法 |
CN108910079A (zh) * | 2018-06-26 | 2018-11-30 | 中国直升机设计研究所 | 一种旋翼试验台动平衡调整相位识别方法 |
Non-Patent Citations (3)
Title |
---|
刘正江等: "旋翼桨涡干扰噪声特性试验技术研究", 《直升机技术》 * |
史勇杰等: "旋翼平行桨-涡干扰噪声的参数影响", 《南京航空航天大学学报》 * |
王亮权等: "基于气弹耦合的旋翼桨涡干扰噪声预测方法", 《航空动力学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112487360A (zh) * | 2020-11-20 | 2021-03-12 | 中国直升机设计研究所 | 一种基于统计分析的旋翼桨涡干扰噪声整周期平均方法 |
CN112504426A (zh) * | 2020-11-20 | 2021-03-16 | 中国直升机设计研究所 | 一种基于峰值搜索的旋翼桨涡干扰噪声整周期平均方法 |
CN112504426B (zh) * | 2020-11-20 | 2022-10-18 | 中国直升机设计研究所 | 一种基于峰值搜索的旋翼桨涡干扰噪声整周期平均方法 |
CN112487360B (zh) * | 2020-11-20 | 2022-12-30 | 中国直升机设计研究所 | 一种基于统计分析的旋翼桨涡干扰噪声整周期平均方法 |
CN114235360A (zh) * | 2021-11-19 | 2022-03-25 | 中国直升机设计研究所 | 一种直升机桨叶实测载荷相位关系的分析方法 |
CN114235360B (zh) * | 2021-11-19 | 2023-10-27 | 中国直升机设计研究所 | 一种直升机桨叶实测载荷相位关系的分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110844115B (zh) | 2022-04-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110844115B (zh) | 一种桨涡干扰噪声和桨叶挥舞载荷数据有效性判别方法 | |
US10591519B2 (en) | Detecting faults in wind turbines | |
CN103890383B (zh) | 控制风力涡轮机的方法及相关系统 | |
EP3557050B1 (en) | Doubly-fed wind turbine blade imbalance detection method based on stator current data | |
JP5033033B2 (ja) | 水平軸風車の乱流強度計測方法 | |
ES2764974T3 (es) | Mejoras relacionadas con la determinación de desequilibrios del rotor en una turbina eólica | |
KR101974293B1 (ko) | 회전 기계의 상태 감시 장치, 회전 기계, 및 회전 기계의 상태 감시 방법 | |
ES2911511T3 (es) | Sistema y método de diagnóstico de un desequilibrio del rotor de un aerogenerador | |
TW201741561A (zh) | 風力渦輪機偏航失準估測 | |
CN109083815B (zh) | 用于风力涡轮操作的偏移切换方法 | |
CN105654239B (zh) | 用于风力发电机组极端风况的识别方法、装置及系统 | |
CN109724802B (zh) | 一种基于频谱图评价与优选的电机轴承微弱故障诊断方法 | |
CN105756864A (zh) | 基于双馈风力发电机组定子电流的桨叶不平衡故障诊断方法 | |
CN102749478A (zh) | 测量风力涡轮机转子尾流中的风向的方法 | |
CN108700492B (zh) | 异常诊断装置和异常诊断方法 | |
JP6184771B2 (ja) | タービンブレードの状態監視方法及び装置 | |
SA522432997B1 (ar) | طريقة ونظام لمراقبة صلاحية مُثبت جذر الريشة | |
CN110794367B (zh) | 一种桨涡干扰声源定位方法 | |
CN113883014B (zh) | 风电机组叶轮不平衡检测方法、装置、设备及存储介质 | |
EP3619427A1 (en) | Determining a wind speed value | |
CN212454697U (zh) | 一种风力发电机组叶片运行状态检测装置 | |
CN111173687B (zh) | 一种风电风机叶片裂纹损伤在线监测装置及方法 | |
CN110920931A (zh) | 一种旋转部件飞行测试静态载荷、动态载荷分离方法 | |
CN113504310B (zh) | 基于单个叶端定时传感器的叶片固有频率识别方法 | |
CN111103451B (zh) | 核电站汽轮机转子划痕测量信号处理方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |