CN110794367A - 一种桨涡干扰声源定位方法 - Google Patents

一种桨涡干扰声源定位方法 Download PDF

Info

Publication number
CN110794367A
CN110794367A CN201910971201.9A CN201910971201A CN110794367A CN 110794367 A CN110794367 A CN 110794367A CN 201910971201 A CN201910971201 A CN 201910971201A CN 110794367 A CN110794367 A CN 110794367A
Authority
CN
China
Prior art keywords
blade
peak
interference noise
vortex interference
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
Application number
CN201910971201.9A
Other languages
English (en)
Other versions
CN110794367B (zh
Inventor
刘正江
陈焕
陈卫星
廖会生
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Helicopter Research and Development Institute
Original Assignee
China Helicopter Research and Development Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Helicopter Research and Development Institute filed Critical China Helicopter Research and Development Institute
Priority to CN201910971201.9A priority Critical patent/CN110794367B/zh
Publication of CN110794367A publication Critical patent/CN110794367A/zh
Application granted granted Critical
Publication of CN110794367B publication Critical patent/CN110794367B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
    • 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)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明属直升机模型旋翼试验技术领域，涉及一种桨涡干扰声源定位方法。该方法包括：在一片桨叶的N个剖面处贴挥舞载荷片,其中N为不小于5的正整数；同步采集转速、桨涡干扰噪声和所述N个剖面位置处的桨叶挥舞载荷；根据转速、桨涡干扰噪声和桨叶挥舞载荷计算桨叶上的桨涡干扰噪声具体声源位置。本方法能够计算出桨叶上出现桨涡干扰噪声的具体声源位置而不需要专门的传声器阵列，不需要消声室环境。

Description

一种桨涡干扰声源定位方法
技术领域
本发明属直升机模型旋翼试验技术领域，涉及一种桨涡干扰声源定位方法。
背景技术
旋翼桨涡干扰噪声是直升机旋翼产生的一种典型噪声，尤其在小速度前飞和斜下降状态，桨涡干扰噪声的能量占到直升机噪声总能量的80％。因此它是目前直升机降低噪声的主要研究领域。桨叶载荷反应的是桨叶承受多大的交变力，是检测并判断桨叶是否工作正常的一种常规手段。桨涡干扰是桨叶切割前一片桨叶脱落桨尖涡产生的气动扰动，该扰动同时也会对桨叶载荷尤其是挥舞载荷产生激励作用，两者之间存在内在的联系。
目前国内外对旋翼桨涡干扰噪声声源定位主要采用传声器阵列的方法来采集声源数据，然后进行空间计算生成定位云图，该方法需要较多的传声器来组成声阵列，同时声阵列支架需要安装在试验台附近尽量没有紊流场的位置，而且每次试验前需要对阵列进行校准检查，因此声阵列只能用于特殊处理过的消声环境下的旋翼试验台进行声源定位试验才能获得理想的结果。
本方法通过同步采集转速方位信号、桨涡干扰信号和桨叶多个剖面的挥舞载荷数据进行整周期平均，提取出桨涡干扰噪声和多个剖面桨叶挥舞载荷最大峰值处的相位，换算成桨涡干扰声源位置到多个桨叶挥舞载荷剖面的距离，从而反向计算出桨叶剖面出现桨涡干扰噪声的具体声源位置。
发明内容
本发明的目的：提供一种利用桨叶挥舞载荷数据进行桨涡干扰声源方法，通过同步采集转速方位信号、桨涡干扰信号和桨叶多个剖面的挥舞载荷数据进行整周期平均，提取出桨涡干扰噪声和多个剖面桨叶挥舞载荷最大峰值处的相位，换算成桨涡干扰声源位置到多个桨叶挥舞载荷剖面的距离，从而计算出桨叶剖面出现桨涡干扰噪声的具体声源位置。
本发明的技术方案：
第一方面，提供了一种桨涡干扰声源定位方法，包括：
在一片桨叶的N个剖面处贴挥舞载荷片,其中N为不小于5的正整数；
同步采集转速、桨涡干扰噪声和所述N个剖面位置处的桨叶挥舞载荷；
根据转速、桨涡干扰噪声和桨叶挥舞载荷计算桨叶上的桨涡干扰噪声具体声源位置。
可选地，在一片桨叶的N个剖面位置处贴挥舞载荷片，具体包括：在一片桨叶的0.95Lr剖面、0.90Lr剖面、0.85Lr剖面、0.80Lr剖面、0.75Lr剖面、0.70Lr剖面、0.65Lr剖面、0.60Lr剖面、0.55Lr剖面位置和0.5Lr剖面位置处粘贴挥舞载荷应变片，其中Lr是桨叶长度。
可选地，同步采集转速、桨涡干扰噪声和桨叶挥舞载荷，具体包括：采集转速；根据转速设置采样率；按照所述采样率同步采集不少于100转的转速、桨涡干扰噪声和桨叶挥舞载荷。
可选地，根据转速、桨涡干扰噪声和桨叶挥舞载荷计算桨叶上的桨涡干扰噪声具体声源位置，具体包括：
根据转速、桨涡干扰噪声和桨叶挥舞载荷计算得到桨涡干扰噪声的峰值位置和桨叶挥舞载荷的峰值位置；
根据桨涡干扰噪声的峰值位置和桨叶挥舞载荷的峰值位置计算桨叶上的桨涡干扰噪声具体声源位置。
可选地，根据转速、桨涡干扰噪声和桨叶挥舞载荷计算得到桨涡干扰噪声的峰值位置和桨叶挥舞载荷的峰值位置，具体包括：
根据转速对桨涡干扰噪声和桨叶挥舞载荷分别进行整周期平均得到整周期平均的桨涡干扰噪声和整周期平均的桨叶挥舞载荷；
分别提取整周期平均的桨涡干扰噪声的峰值相位和整周期平均的桨叶挥舞载荷的峰值相位计算得到桨涡干扰噪声的峰值位置和桨叶挥舞载荷的峰值位置。
可选地，根据转速对桨涡干扰噪声和桨叶挥舞载荷分别进行整周期平均得到整周期平均的桨涡干扰噪声和整周期平均的桨叶挥舞载荷，具体包括：
以转速的任意一个脉冲上升沿为起点，截取数据点数为360000的桨涡干扰噪声数据和桨叶挥舞载荷数据；
将桨涡干扰噪声数据和桨叶挥舞载荷数据依次按照数据点数为3600等间隔分割成100组数据；
将分割后的100组桨涡干扰噪声中的每组桨涡干扰噪声的第1个点、第2个点到第3600个点累加并除以100，得到整周期平均的桨涡干扰噪声；
将分割后的100组桨叶挥舞载荷中的每组桨叶挥舞载荷的第1个点、第2个点到第3600个点累加并除以100，得到整周期平均的桨叶挥舞载荷。
可选地，分别提取整周期平均的桨涡干扰噪声的峰值相位和整周期平均的桨叶挥舞载荷的峰值相位得到桨涡干扰噪声的峰值相位和桨叶挥舞载荷的峰值相位，具体包括：
将整周期平均的桨涡干扰噪声的第三个波峰的峰值相位确定为桨涡干扰噪声的峰值相位；
计算整周期平均的桨涡干扰噪声的最小峰值相位，根据整周期平均的桨涡干扰噪声的最小峰值相位确定整周期平均的桨叶挥舞载荷的峰值参考数据位置点；
根据整周期平均的桨叶挥舞载荷的峰值参考数据位置点确定桨叶挥舞载荷的峰值相位。
可选地，根据桨涡干扰噪声的峰值位置和桨叶挥舞载荷的峰值位置计算桨叶上的桨涡干扰噪声具体声源位置，具体包括：
将多个剖面处的桨叶挥舞载荷的峰值位置组成峰值位置数组；
将峰值位置数组划分成第一峰值位置数组Y1和第二峰值位置数组Y2；
根据第一峰值位置数组Y1和第二峰值位置数组Y2得到对应的第一剖面位置数组X1和第二剖面位置数组X2；
利用第一峰值位置数组Y1和第一剖面位置数组X1以及第二峰值位置数组Y2和第二剖面位置数组X2，计算得到桨涡干扰噪声源位置。
可选地，将峰值位置数组划分成第一峰值位置数组和第二峰值位置数组，具体包括：
确定峰值位置数组中的最小值；
将峰值位置数组中的最小值及其左侧的数据组成第一峰值位置数组Y1；将所述数组中的最小值右侧的数据组成第二峰值位置数组Y2。
可选地，利用第一峰值位置数组Y1和第一剖面位置数组X1以及第二峰值位置数组Y2和第二剖面位置数组X2，计算得到桨涡干扰噪声源位置，具体包括：
利用第一峰值位置数组Y1和第一剖面位置数组X1，利用最小二乘法拟合得到相应的第一方程Y1＝K1*X1+B1；
利用第二峰值位置数组Y2和第二剖面位置数组X2，利用最小二乘法拟合得到相应的第二方程Y2＝K2*X2+B2；
将第一方程和第二方程联立求解得到桨涡干扰噪声源位置。
本发明的有益效果：
本发明的桨涡干扰声源定位方法，通过同步采集转速方位信号、桨涡干扰信号和桨叶多个剖面的挥舞载荷数据进行整周期平均，提取出桨涡干扰噪声和多个剖面桨叶挥舞载荷最大峰值处的相位，换算成桨涡干扰声源位置到多个桨叶挥舞载荷剖面的距离，从而计算出桨叶上出现桨涡干扰噪声的具体声源位置，不需要专门的传声器阵列，不需要消声室环境，具有低成本且便捷、高效的优点。
附图说明
图1为根据本发明实施例的桨涡干扰声源定位方法的流程图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得
的所有其他实施例，都属于本发明保护的范围。
下面对本发明做进一步详细说明。
(1)选择一片桨叶N个剖面贴挥舞载荷片
第1步，选择一片桨叶(假设为第3号桨叶)，在其0.95Lr剖面(Lr是桨叶长度，单位为米)、0.90Lr剖面、0.85Lr剖面、0.80Lr剖面、0.75Lr剖面、0.70Lr剖面、0.65Lr剖面、0.60Lr剖面、0.55Lr剖面位置和0.5Lr剖面位置粘贴挥舞载荷应变片
(2)转速和桨涡干扰噪声、桨叶挥舞载荷同步采集
转速和桨涡干扰噪声、桨叶挥舞载荷同步采集步骤如下：
第1步，将涡发生器安装在和桨盘平面同高的位置，涡发生器出涡口正对桨尖位置，距离为100mm；
第2步，已涡发生器位置为0度方位，以距离桨榖中心为圆心，以1.5倍桨叶长度Lr(Lr的单位为米)为半径，按桨叶旋转方向顺时针旋转45度，在该位置竖直布置一根传声器枝干，然后从桨盘平面高度6米开始，在传声器枝干上以0.5米等间隔布置6个点传声器，与桨盘平面相同高度的为1号传声器，低于桨盘平面0.5米的为2号传声器，依次类推最低的点为6号传声器；
第3步，将转速传感器安装在涡发生器的同一个方位，在贴有挥舞载荷应变片的桨叶(第3号桨叶)对面的桨叶(以下假定为第1号桨叶)上贴上反光片；
第4步，根据转速Np(Np的单位为转/秒)，设置采样率Sm＝Np*3600，然后待试验状态稳定后同步采集不少于Tn＝100/Np秒的转速和桨涡干扰噪声、桨叶挥舞载荷数据，既至少采集100圈的数据，同时要求试验时操纵总距在±0.1度范围以内，涡发生器的出涡口速度不小于11米/秒。
(3)转速桨涡干扰噪声、桨叶挥舞载荷试验数据整周期平均
整周期平均步骤如下：
第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为起点(即0点)，分别从在0.95Lr剖面、0.90Lr剖面、0.85Lr剖面、0.80Lr剖面、0.75Lr剖面、0.70Lr剖面、0.65Lr剖面、0.60Lr剖面、0.55Lr剖面位置和0.5Lr剖面位置的挥舞载荷整周期平均数据向前查找到出现第1个波峰的峰值点的位置分别对应为：Lh1-0.95、Lh1-0.90、Lh1-0.85、Lh1-0.80、Lh1-0.75、Lh1-0.70、Lh1-0.65、Lh1-0.60、Lh1-0.55、Lh1-0.50，这10个数组成一个数组。
(4)桨叶上桨涡干扰噪声具体声源位置计算
桨叶上的桨涡干扰噪声源位置计算如下：
第1步，从数组[Lh1-0.95，Lh1-0.90，Lh1-0.85，Lh1-0.80，Lh1-0.75，Lh1-0.70，Lh1-0.65，Lh1-0.60，Lh1-0.55，Lh1-0.50]中找到数值最小数据，假设为Lh1-0.80，该数组分割成2个数组，分别是第一峰值位置数组Y1＝[Lh1-0.95，Lh1-0.90，Lh1-0.85，Lh1-0.80]和第二峰值位置数组Y2＝[Lh1-0.80，Lh1-0.75，Lh1-0.70，Lh1-0.65，Lh1-0.60，Lh1-0.55，Lh1-0.50]，同时得到其剖面位置数组，分别是第一剖面位置数组X1＝[0.95，0.90，0.85，0.80]和第二剖面位置数组X2＝[0.80，0.75，0.70，0.65，0.60,0.55,0.50]；
第2步,由第一剖面位置数组X1＝[0.95，0.90，0.85，0.80]，第一峰值位置数组Y1＝[Lh1-0.95，Lh1-0.90，Lh1-0.85，Lh1-0.80]，按照最小二乘法拟合计算得到对应的一次方程Y1＝K1*X1+B1；
第3步,由第一剖面位置数组X2＝[0.80，0.75，0.70，0.65，0.60,0.55,0.50]，第一峰值位置数组Y1＝[Lh1-0.80，Lh1-0.75，Lh1-0.70，Lh1-0.65，Lh1-0.60，Lh1-0.55，Lh1-0.50]，按照最小二乘法拟合计算得到对应的一次方程Y2＝K2*X2+B2；
第4步，对两个方程取相同的自变量和应变量，即Y＝K1*X+B1和Y＝K2*X+B2，将两个一元一次方程联立求解得到X＝(B1-B2)/(K1-K2)、Y＝K1*[(B1-B2)/(K1-K2)]+B1。其中的X就是桨叶上的桨涡干扰噪声源位置。
示例如下：
(1)选择一片桨叶N个剖面贴挥舞载荷片
第1步，选择一片桨叶(假设为第3号桨叶)，在其0.95Lr剖面(Lr是桨叶长度，单位为米)、0.90Lr剖面、0.85Lr剖面、0.80Lr剖面、0.75Lr剖面、0.70Lr剖面、0.65Lr剖面、0.60Lr剖面、0.55Lr剖面位置和0.5Lr剖面位置粘贴挥舞载荷应变片
(2)转速和桨涡干扰噪声、桨叶挥舞载荷同步采集
转速和桨涡干扰噪声、桨叶挥舞载荷同步采集步骤如下：
第1步，将涡发生器安装在和桨盘平面同高的位置，涡发生器出涡口正对桨尖位置，距离为100mm；
第2步，已涡发生器位置为0度方位，以距离桨榖中心为圆心，以1.5倍桨叶长度Lr(Lr的单位为米)为半径，按桨叶旋转方向顺时针旋转45度，在该位置竖直布置一根传声器枝干，然后从桨盘平面高度6米开始，在传声器枝干上以0.5米等间隔布置6个点传声器，与桨盘平面相同高度的为1号传声器，低于桨盘平面0.5米的为2号传声器，依次类推最低的点为6号传声器；
第3步，将转速传感器安装在涡发生器的同一个方位，在贴有挥舞载荷应变片的桨叶(第3号桨叶)对面的桨叶(以下假定为第1号桨叶)上贴上反光片；
第4步，根据转速Np＝10(Np的单位为转/秒)，设置采样率Sm＝Np*3600＝36000，然后待试验状态稳定后同步采集不少于Tn＝100/Np＝10秒的转速和桨涡干扰噪声、桨叶挥舞载荷数据，既至少采集100圈的数据，同时要求试验时操纵总距为0度，涡发生器的出涡口速度为15米/秒。
(3)转速桨涡干扰噪声、桨叶挥舞载荷试验数据整周期平均
整周期平均步骤如下：
第1步，以转速任意一个脉冲上升沿为起点，截取数据点长度为360000的桨涡干扰噪声数据和桨叶挥舞载荷数据；
第2步，然后将上述数据依次按照数据点长度为3600等间隔分割成100圈数据，每圈3600个点；
第3步，将分割后的100圈数据每圈的第1个点、第2个点到第3600个点分别提取出来累加后除100，从而得到1圈(即3600个点)平均数据，既整周期平均数据。
(4)桨涡干扰噪声、桨叶挥舞载荷试验数据峰值相位提取
峰值相位提取步骤如下：
第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＝247.9，从而计算得到在1圈平均数据中的桨叶挥舞载荷波峰参考数据位置点Lh0＝Pmin×10＝2479；
第4步，以第3步计算得到的参考数据位置点Lh0为起点(即0点)，别从在0.95Lr剖、0.90Lr剖面、0.85Lr剖面、0.80Lr剖面、0.75Lr剖面、0.70Lr剖面、0.65Lr剖面、0.60Lr剖面、0.55Lr剖面位置和0.5Lr剖面位置的挥舞载荷整周期平均数据向前查找到出现第1个波峰的峰值点的位置分别对应为：549，327，108，3，29，300，354，423，517，759。
(4)桨叶上桨涡干扰噪声的具体声源位置计算
桨叶上的桨涡干扰噪声源位置计算如下：
第1步，从数组[549，327，108，3，29，300，354，423，517，759]中找到数值最小数据为3，该数组分割成2个数组，分别是第一峰值位置数组Y1＝[549，327，108，3]和第二峰值位置数组Y2＝[3，29，300，354，423，517，759]，同时得到其剖面位置数组，分别是第一剖面位置数组X1＝[0.95，0.90，0.85，0.80]和第二剖面位置数组X2＝[0.80，0.75，0.70，0.65，0.60,0.55,0.50]；
第2步,由第一峰值位置数组X1＝[0.95，0.90，0.85，0.80]，第一峰值位置数组Y1＝[549，327，108，3]，按照最小二乘法拟合计算得到对应的一次方程Y1＝3714*X1-3003；
第3步,由第一峰值位置数组X2＝[0.80，0.75，0.70，0.65，0.60,0.55,0.50]，第二峰值位置数组Y2＝[3，29，300，354，423，517，759]，按照最小二乘法拟合计算得到对应的一次方程Y2＝-2405*X2+1904；
第4步，对两个方程取相同的自变量和应变量，即Y＝3714*X-3003和Y＝-2405*X+1904，将两个一元一次方程联立求解得到X＝0.8019,X就是桨叶上的桨涡干扰噪声源位置0.8019Lr。
以上所述，仅为本发明的具体实施例，对本发明进行详细描述，未详尽部分为常规技术。但本发明的保护范围不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。本发明的保护范围应以所述权利要求的保护范围为准。

Claims (10)

1.一种桨涡干扰声源定位方法，其特征在于，包括：
在一片桨叶的N个剖面处贴挥舞载荷片,其中N为不小于5的正整数；
同步采集转速、桨涡干扰噪声和所述N个剖面位置处的桨叶挥舞载荷；
根据转速、桨涡干扰噪声和桨叶挥舞载荷计算桨叶上的桨涡干扰噪声具体声源位置。
2.根据权利要求1所述的方法，其特征在于，在一片桨叶的N个剖面位置处贴挥舞载荷片，具体包括：在一片桨叶的0.95Lr剖面、0.90Lr剖面、0.85Lr剖面、0.80Lr剖面、0.75Lr剖面、0.70Lr剖面、0.65Lr剖面、0.60Lr剖面、0.55Lr剖面位置和0.5Lr剖面位置处粘贴挥舞载荷应变片，其中Lr是桨叶长度。
3.根据权利要求1所述的方法，其特征在于，同步采集转速、桨涡干扰噪声和桨叶挥舞载荷，具体包括：采集转速；根据转速设置采样率；按照所述采样率同步采集不少于100转的转速、桨涡干扰噪声和桨叶挥舞载荷。
4.根据权利要求1所述的方法，其特征在于，根据转速、桨涡干扰噪声和桨叶挥舞载荷计算桨叶上的桨涡干扰噪声具体声源位置，具体包括：
根据转速、桨涡干扰噪声和桨叶挥舞载荷计算得到桨涡干扰噪声的峰值位置和桨叶挥舞载荷的峰值位置；
根据桨涡干扰噪声的峰值位置和桨叶挥舞载荷的峰值位置计算桨叶上的桨涡干扰噪声具体声源位置。
5.根据权利要求4所述的方法，其特征在于，根据转速、桨涡干扰噪声和桨叶挥舞载荷计算得到桨涡干扰噪声的峰值位置和桨叶挥舞载荷的峰值位置，具体包括：
根据转速对桨涡干扰噪声和桨叶挥舞载荷分别进行整周期平均得到整周期平均的桨涡干扰噪声和整周期平均的桨叶挥舞载荷；
分别提取整周期平均的桨涡干扰噪声的峰值相位和整周期平均的桨叶挥舞载荷的峰值相位计算得到桨涡干扰噪声的峰值位置和桨叶挥舞载荷的峰值位置。
6.根据权利要求5所述的方法，其特征在于，根据转速对桨涡干扰噪声和桨叶挥舞载荷分别进行整周期平均得到整周期平均的桨涡干扰噪声和整周期平均的桨叶挥舞载荷，具体包括：
以转速的任意一个脉冲上升沿为起点，截取数据点数为360000的桨涡干扰噪声数据和桨叶挥舞载荷数据；
将桨涡干扰噪声数据和桨叶挥舞载荷数据依次按照数据点数为3600等间隔分割成100组数据；
将分割后的100组桨涡干扰噪声中的每组桨涡干扰噪声的第1个点、第2个点到第3600个点累加并除以100，得到整周期平均的桨涡干扰噪声；
将分割后的100组桨叶挥舞载荷中的每组桨叶挥舞载荷的第1个点、第2个点到第3600个点累加并除以100，得到整周期平均的桨叶挥舞载荷。
7.根据权利要求5所述的方法，其特征在于，分别提取整周期平均的桨涡干扰噪声的峰值相位和整周期平均的桨叶挥舞载荷的峰值相位得到桨涡干扰噪声的峰值相位和桨叶挥舞载荷的峰值相位，具体包括：
将整周期平均的桨涡干扰噪声的第三个波峰的峰值相位确定为桨涡干扰噪声的峰值相位；
计算整周期平均的桨涡干扰噪声的最小峰值相位，根据整周期平均的桨涡干扰噪声的最小峰值相位确定整周期平均的桨叶挥舞载荷的峰值参考数据位置点；
根据整周期平均的桨叶挥舞载荷的峰值参考数据位置点确定桨叶挥舞载荷的峰值相位。
8.根据权利要求4所述的方法，其特征在于，根据桨涡干扰噪声的峰值位置和桨叶挥舞载荷的峰值位置计算桨叶上的桨涡干扰噪声具体声源位置，具体包括：
将多个剖面处的桨叶挥舞载荷的峰值位置组成峰值位置数组；
将峰值位置数组划分成第一峰值位置数组Y1和第二峰值位置数组Y2；
根据第一峰值位置数组Y1和第二峰值位置数组Y2得到对应的第一剖面位置数组X1和第二剖面位置数组X2；
利用第一峰值位置数组Y1和第一剖面位置数组X1以及第二峰值位置数组Y2和第二剖面位置数组X2，计算得到桨涡干扰噪声源位置。
9.根据权利要求1所述的方法，其特征在于，将峰值位置数组划分成第一峰值位置数组和第二峰值位置数组，具体包括：
确定峰值位置数组中的最小值；
将峰值位置数组中的最小值及其左侧的数据组成第一峰值位置数组Y1；将所述数组中的最小值右侧的数据组成第二峰值位置数组Y2。
10.根据权利要求1所述的方法，其特征在于，利用第一峰值位置数组Y1和第一剖面位置数组X1以及第二峰值位置数组Y2和第二剖面位置数组X2，计算得到桨涡干扰噪声源位置，具体包括：
利用第一峰值位置数组Y1和第一剖面位置数组X1，利用最小二乘法拟合得到相应的第一方程Y1＝K1*X1+B1；
利用第二峰值位置数组Y2和第二剖面位置数组X2，利用最小二乘法拟合得到相应的第二方程Y2＝K2*X2+B2；
将第一方程和第二方程联立求解得到桨涡干扰噪声源位置。
CN201910971201.9A 2019-10-12 2019-10-12 一种桨涡干扰声源定位方法 Active CN110794367B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910971201.9A CN110794367B (zh) 2019-10-12 2019-10-12 一种桨涡干扰声源定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910971201.9A CN110794367B (zh) 2019-10-12 2019-10-12 一种桨涡干扰声源定位方法

Publications (2)

Publication Number Publication Date
CN110794367A true CN110794367A (zh) 2020-02-14
CN110794367B CN110794367B (zh) 2022-06-21

Family

ID=69440257

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910971201.9A Active CN110794367B (zh) 2019-10-12 2019-10-12 一种桨涡干扰声源定位方法

Country Status (1)

Country Link
CN (1) CN110794367B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112504426A (zh) * 2020-11-20 2021-03-16 中国直升机设计研究所 一种基于峰值搜索的旋翼桨涡干扰噪声整周期平均方法
CN114598983A (zh) * 2022-01-24 2022-06-07 北京航空航天大学 一种民机增升装置噪声麦克风阵列测试方法

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
BR9712021A (pt) * 1996-09-10 2000-01-18 Allison Engine Company Aparelho e método para redução de ruìdo de turbo propulsor.
JP2002333367A (ja) * 2001-05-10 2002-11-22 Mitsuba Corp 音源探査方法および装置
CN102588211A (zh) * 2012-02-29 2012-07-18 沈阳华人风电科技有限公司 一种风力发电机组全工况模型预测控制方法及系统
CN102722606A (zh) * 2012-05-24 2012-10-10 北京航空航天大学 一种降低直升机旋翼桨毂振动载荷的方法
US20140112103A1 (en) * 2012-10-22 2014-04-24 U.S. Army Research Laboratory Attn: Rdrl-Loc-I Method and system for motion compensated target detection using acoustical focusing
WO2014124643A1 (en) * 2013-02-14 2014-08-21 Vestas Wind Systems A/S Detecting blade structure abnormalities
CN104655081A (zh) * 2013-11-25 2015-05-27 中国直升机设计研究所 一种桨叶结构预扭角测量方法
KR20180018164A (ko) * 2016-08-12 2018-02-21 한국항공우주연구원 무인 항공기 및 이의 제어 방법
CN108120592A (zh) * 2017-11-29 2018-06-05 中国直升机设计研究所 一种直升机桨叶静强度的试验方法
CN108910079A (zh) * 2018-06-26 2018-11-30 中国直升机设计研究所 一种旋翼试验台动平衡调整相位识别方法
US20190035288A1 (en) * 2017-07-27 2019-01-31 Intel Corporation Methods, systems, and apparatus for drone collision avoidance and acoustic detection

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
BR9712021A (pt) * 1996-09-10 2000-01-18 Allison Engine Company Aparelho e método para redução de ruìdo de turbo propulsor.
JP2002333367A (ja) * 2001-05-10 2002-11-22 Mitsuba Corp 音源探査方法および装置
CN102588211A (zh) * 2012-02-29 2012-07-18 沈阳华人风电科技有限公司 一种风力发电机组全工况模型预测控制方法及系统
CN102722606A (zh) * 2012-05-24 2012-10-10 北京航空航天大学 一种降低直升机旋翼桨毂振动载荷的方法
US20140112103A1 (en) * 2012-10-22 2014-04-24 U.S. Army Research Laboratory Attn: Rdrl-Loc-I Method and system for motion compensated target detection using acoustical focusing
WO2014124643A1 (en) * 2013-02-14 2014-08-21 Vestas Wind Systems A/S Detecting blade structure abnormalities
CN104655081A (zh) * 2013-11-25 2015-05-27 中国直升机设计研究所 一种桨叶结构预扭角测量方法
KR20180018164A (ko) * 2016-08-12 2018-02-21 한국항공우주연구원 무인 항공기 및 이의 제어 방법
US20190035288A1 (en) * 2017-07-27 2019-01-31 Intel Corporation Methods, systems, and apparatus for drone collision avoidance and acoustic detection
CN108120592A (zh) * 2017-11-29 2018-06-05 中国直升机设计研究所 一种直升机桨叶静强度的试验方法
CN108910079A (zh) * 2018-06-26 2018-11-30 中国直升机设计研究所 一种旋翼试验台动平衡调整相位识别方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
刘正江 等: "旋翼桨涡干扰噪声特性试验技术研究", 《直升机技术》 *
刘正江 等: "旋翼桨涡干扰噪声特性试验技术研究", 《直升机技术》, no. 1, 31 March 2019 (2019-03-31), pages 43 - 47 *
徐国华 等: "直升机旋翼气动噪声的研究新进展", 《航空学报》 *
徐国华 等: "直升机旋翼气动噪声的研究新进展", 《航空学报》, vol. 38, no. 7, 31 December 2017 (2017-12-31), pages 23 - 38 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112504426A (zh) * 2020-11-20 2021-03-16 中国直升机设计研究所 一种基于峰值搜索的旋翼桨涡干扰噪声整周期平均方法
CN112504426B (zh) * 2020-11-20 2022-10-18 中国直升机设计研究所 一种基于峰值搜索的旋翼桨涡干扰噪声整周期平均方法
CN114598983A (zh) * 2022-01-24 2022-06-07 北京航空航天大学 一种民机增升装置噪声麦克风阵列测试方法

Also Published As

Publication number Publication date
CN110794367B (zh) 2022-06-21

Similar Documents

Publication Publication Date Title
DK2169218T3 (en) System and method for estimating wind conditions for wind turbines
CN110794367B (zh) 一种桨涡干扰声源定位方法
CN102749478A (zh) 测量风力涡轮机转子尾流中的风向的方法
Barlas et al. Consistent modelling of wind turbine noise propagation from source to receiver
CN100390406C (zh) 涡轮机
CN206290378U (zh) 风力发电机组监测系统及风力发电机组
CN103245291B (zh) 叶片类零件装配精度检测方法
WO2018107726A1 (zh) 基于定子电流数据驱动的双馈风电机组桨叶不平衡检测方法
WO2012083958A3 (en) Control method for a wind turbine
Oerlemans et al. Experimental demonstration of wind turbine noise reduction through optimized airfoil shape and trailing-edge serrations
CN110844115B (zh) 一种桨涡干扰噪声和桨叶挥舞载荷数据有效性判别方法
CN105138845B (zh) 获得风力发电机风速值的方法
Castelli et al. Effect of dynamic stall, finite aspect ratio and streamtube expansion on VAWT performance prediction using the BE-M model
Baskaran et al. Aerodynamic and aeroacoustic characteristics of propellers with different blade numbers
Bedon et al. Experimental tests of a vertical-axis wind turbine with twisted blades
CN209354294U (zh) 一种用于实海况测量的潮流能水轮机叶尖比速测量装置
NO325790B1 (no) Vindturbin med vertikal akse
CN103452766B (zh) 一种叶片叶尖带v型开槽结构的风力机
CN117131608A (zh) 一种基于最佳环量分布的激励盘方法
CN110844110B (zh) 一种确定桨叶运动参数的相位的方法
CN102661787A (zh) 一种变转速风电机组转子振动信号的变换分析方法
CN107304749B (zh) 一种风力发电机组机舱位置测量的方法、装置和系统
EP3504426A1 (en) Wind data detecting system and method for wind turbines
CN105332856A (zh) 一种基于测量固定坐标系下主轴载荷的风电机组独立变桨控制方法
JP2014047742A (ja) 風力発電装置、及び、風力発電装置の制御方法

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