CN107220412A - 一种基于可变面积阻力片的旋翼厚度噪声控制方法 - Google Patents
一种基于可变面积阻力片的旋翼厚度噪声控制方法 Download PDFInfo
- Publication number
- CN107220412A CN107220412A CN201710307028.3A CN201710307028A CN107220412A CN 107220412 A CN107220412 A CN 107220412A CN 201710307028 A CN201710307028 A CN 201710307028A CN 107220412 A CN107220412 A CN 107220412A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- blade
- rotor
- resisting force
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/10—Noise analysis or noise optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Toys (AREA)
- Wind Motors (AREA)
Abstract
本发明实施例公开了一种基于可变面积阻力片的旋翼厚度噪声控制方法,涉及气动声学技术领域,能够针对性地对旋翼厚度噪声进行大幅度降噪。本发明包括:读取旋翼尺寸和飞行参数,计算给定飞行状态下旋翼厚度噪声P′T;在桨叶端部四分之一弦线处施加控制力F,计算给定飞行状态下控制力F产生的控制源噪声p′F;根据旋翼厚度噪声降噪目标,设定抵消因子f,计算桨叶端部四分之一弦线处施加的控制力F;在桨叶端部四分之一弦线处添加可变面积阻力片,基于CFD方法计算不同阻力片面积与控制力的关系,得到桨叶端部四分之一弦线处阻力片的面积;改变控制力展向分布位置,计算不同位置的阻力片面积;对桨叶不同位置添加相应面积的阻力片。
Description
技术领域
本发明涉及气动声学技术领域,尤其涉及一种基于可变面积阻力片的旋翼厚度噪声控制方法。
背景技术
随着直升机的广泛应用,直升机噪声大的缺点越来越引起人们的重视,使得直升机噪声已经成为直升机设计过程中需要着重考虑的问题。旋翼噪声是直升机噪声中影响最重要的部分,旋翼噪声可分为厚度噪声、载荷噪声、桨-涡干扰噪声和高速脉冲噪声等。旋翼厚度噪声是由于桨叶周期性地排开空气而产生的,具有低频及桨盘面内传播的特性,传播距离较远,是声探测的主要对象,探索能有效抑制、甚至消除旋翼厚度噪声不利影响的控制措施不仅具有重要科学研究意义,更有十分广阔的应用前景,可为未来绿色、安静直升机的发展提供重要技术支撑。
直升机旋翼厚度噪声与桨叶外形密切相关。在过去的数十年间,国内外研究人员主要通过对旋翼外形进行优化设计,从而降低旋翼噪声。美国RAH-66直升机旋翼桨叶采用后掠桨尖,使噪声水平降低2dB左右;黑鹰直升机经过历次改进设计,旋翼桨尖由最初的矩形变为下反尖削形式,以用于减小噪声水平;欧洲研制的Blue-edge(蓝色刀锋)旋翼桨叶采用了双掠型桨叶(前掠-后掠组合),使得平飞状态下的旋翼气动噪声可以降低2-3dB。这类基于桨叶气动外形设计的被动噪声控制方法虽然能在一定程度上降低旋翼噪声,但是存在以下几个方面的缺点:1)旋翼气动性能直接决定于桨叶外形,因此为了保证旋翼具有优良的气动性能,在进行低噪声旋翼设计时需要在噪声和性能设计之间进行折衷,从而使得噪声控制效果往往不明显;2)没有针对性地进行旋翼厚度噪声降噪,导致厚度噪声的降噪效果较差;3)被动噪声控制方法存在着对于不同飞行状态的适应性问题,不能实现全包线飞行状态下的直升机降噪。
发明内容
本发明的实施例提供一种基于可变面积阻力片的旋翼厚度噪声控制方法,在旋翼桨叶端部设置可变面积阻力片,通过控制***使阻力片面积发生非定常的脉动,从而诱发出一个与被控噪声相位相同、幅值相反的可控声波,通过原旋翼噪声声波与所述可控声波的叠加,直接抵消某一方向上的声压脉冲,进而从根本上降低直升机飞行中的噪声辐射。
为达到上述目的,本发明的实施例采用如下技术方案:
第一方面,本发明实施例提供一种基于可变面积阻力片的旋翼厚度噪声控制方法,包括以下步骤:
S1,读取旋翼尺寸和飞行参数,建立基于声学类比法的旋翼厚度噪声计算模型,计算给定飞行状态下旋翼厚度噪声PT';
S2,在桨叶端部四分之一弦线处施加控制力F,建立基于声学类比法的控制源噪声计算模型,计算给定飞行状态下控制力F产生的控制源噪声p′F;
S3,根据旋翼厚度噪声降噪目标,设定抵消因子f,根据所述旋翼厚度噪声PT'和控制源噪声PT',计算桨叶端部四分之一弦线处施加的控制力F;
S4,在桨叶端部四分之一弦线处添加可变面积阻力片,基于CFD方法计算不同阻力片面积与控制力的关系,得到桨叶端部四分之一弦线处阻力片的面积;
S5,改变控制力展向分布位置,根据所述不同阻力片面积与控制力的关系,计算不同位置的阻力片面积;
S6,根据步骤S5的结果,对桨叶不同位置添加相应面积的阻力片。
作为一个实施例,还包括:S7,对带有相应面积阻力片的旋翼,根据步骤S1计算添加阻力片之后的旋翼厚度噪声,进行检验和完善。
作为一个实施例,所述步骤S1中,选取桨盘平面内5R处的点为远场目标观测点,计算给定飞行状态下控制力旋翼的厚度噪声PT';其中,R是桨叶半径。
作为一个实施例,所述控制源噪声p′F是:
其中,MH为旋翼悬停时的桨尖马赫数,R是桨叶半径,Ψ是桨叶方位角,F为控制力,和分别表示F和ψ的偏导数,a0是声速,d是声源到远场目标观测点的距离,xobs是远场目标观测点的位置,τ是源时间,Ω是旋翼转速。
作为一个实施例,所述步骤S3远场目标观测点处的厚度噪声为:
f·p′T+p′F≈0
其中,抵消因子f=1表示声场完全对消,0<f<1代表部分抵消。
作为一个实施例,步骤S3中,所述计算桨叶端部四分之一弦线处施加的控制力F是:
其中,MH为旋翼悬停时的桨尖马赫数,As为桨盘面积,ρ0是空气密度,a0是声速,Ψ是桨叶方位角。
作为一个实施例,步骤S4中所述基于CFD方法计算不同阻力片面积与控制力的关系,包括:
S101,读取桨叶几何参数、运动特点和旋翼所处流场特征划分旋翼桨叶的贴体网格和背景网格;
S102,对桨叶贴体网格和背景网格进行前处理,并对带有阻力片部分的桨叶进行网格加密;
S103,对桨叶网格在背景网格中进行挖洞和贡献单元搜索,并将得到的桨叶网格和背景网格的嵌套关系进行储存;
S104,初始化桨叶流场;
S105,更新所述桨叶网格坐标;
S106,更新所述桨叶网格和背景网格的嵌套关系;
S107,计算桨叶网格上的流场,更新背景网格上对应洞边界单元的流场信息。
作为一个实施例,步骤S5中所述改变控制力展向分布位置的改变范围是:0.6R至R,R是桨叶半径。
本发明实施例提供的一种基于可变面积阻力片的旋翼厚度噪声控制方法,通过在桨叶不同位置添加可变面积阻力片,可以针对性地对旋翼厚度噪声进行大幅度降噪,旋翼厚度噪声声压可降低6dB以上,噪声声压级降低50%以上。本发明实施例提供的方法属于主动控制方法,解决了被动控制方法中存在的不同飞行状态的适应性问题,实现了全包线飞行状态下的旋翼厚度噪声降噪。同时,相比现有的旋翼主动控制方法,实现方式简单有效,可显著降低降噪成本。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1为本发明实施例的流程示意图;
图2为本发明实施例提供的旋翼厚度噪声的控制原理图;
图3为本发明实施例提供的带有可变面积阻力片的桨叶示意图;
图4为本发明实施例提供的控制过程中同一位置处阻力片的面积变化示意图;
图5为采用本发明实施例后的旋翼厚度噪声控制效果;
图6为本发明实施例提供的点对点噪声控制图;
图7为本发明实施例提供的控制力随桨叶方位角变化示意图。
具体实施方式
为使本领域技术人员更好地理解本发明的技术方案,下面结合附图和具体实施方式对本发明作进一步详细描述。下文中将详细描述本发明的实施方式,所述实施方式的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
本技术领域技术人员可以理解,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。
本发明的实施例提供一种基于可变面积阻力发生装置的旋翼厚度噪声控制方法,在旋翼桨叶端部设置可变面积阻力片,通过控制***使阻力片面积发生非定常的脉动,如图2所示,从而诱发出一个与被控噪声相位相同、幅值相反的可控声波,通过原旋翼噪声声波与所述可控声波的叠加,直接抵消某一方向上的声压脉冲,进而从根本上降低直升机飞行中的噪声辐射。
为达到上述目的,本发明实施例以一副两米直径旋翼为例,对本发明实施例作进一步的详细说明。旋翼的具体参数如下表所示:
旋翼翼型 | NACA0012 | 旋翼半径(m) | 1 |
旋翼根切(m) | 0.12 | 旋翼总距(°) | 8 |
桨尖马赫数 | 0.6 | 桨叶扭度(°) | 0 |
桨叶弦长(m) | 0.1905 | 桨叶片数 | 2 |
第一方面,本发明实施例提供一种基于可变面积阻力片的旋翼厚度噪声控制方法,如图1所示,包括以下步骤:
S1,读取旋翼尺寸和飞行参数,建立基于声学类比法的旋翼厚度噪声计算模型,计算给定飞行状态下旋翼厚度噪声PT';
S2,在桨叶端部四分之一弦线处施加控制力F,建立基于声学类比法的控制源噪声计算模型,计算给定飞行状态下控制力F产生的控制源噪声p′F;
S3,根据旋翼厚度噪声降噪目标,设定抵消因子f,根据所述旋翼厚度噪声PT'和控制源噪声PT',计算桨叶端部四分之一弦线处施加的控制力F;
S4,在桨叶端部四分之一弦线处添加可变面积阻力片,如图3所示,基于CFD方法计算不同阻力片面积与控制力的关系,得到桨叶端部四分之一弦线处阻力片的面积;
S5,改变控制力展向分布位置,如图4所示,是本发明实施例提供的同一位置阻力片面积变化示意图,根据所述不同阻力片面积与控制力的关系,计算不同位置的阻力片面积;
S6,根据步骤S5的结果,对桨叶不同位置添加相应面积的阻力片。如图5所示,是采用本发明实施例方法后位于桨叶端部四分之一弦线处的阻力片对厚度噪声降噪效果图。
作为一个实施例,还包括:S7,对带有相应面积阻力片后的旋翼,根据步骤S1计算添加阻力片之后的旋翼厚度噪声,进行检验和完善。
作为一个实施例,如图6所示,所述步骤S1中,考虑90°方位角处桨尖外端的声源对远场目标观测点的噪声控制,选取桨盘平面内5R(R是桨叶半径)处的点为远场目标观测点,计算给定飞行状态下控制力旋翼的厚度噪声PT'。其中,所述桨盘是桨叶旋转构成的平面。
作为一个实施例,所述控制源噪声p′F是:
其中,MH为旋翼悬停时的桨尖马赫数,R是桨叶半径,Ψ是桨叶方位角,F为控制力,和分别表示F和ψ的偏导数,a0是声速,d是声源到远场目标观测点的距离,xobs是远场目标观测点的位置,τ是源时间,Ω是旋翼转速。
作为一个实施例,所述步骤S3远场目标观测点处的厚度噪声为:
f·p′T+p′F≈0
其中,抵消因子f=1表示声场完全对消,0<f<1代表部分抵消。
作为一个实施例,步骤S3中,结合设定的抵消因子f,将所述旋翼的厚度噪声PT'和所述控制源噪声p′F代入,计算桨叶端部四分之一弦线处控制力F是:
其中,MH为旋翼悬停时的桨尖马赫数,As为桨盘面积,ρ0是空气密度,a0是声速,Ψ是桨叶方位角。如图7所示,是本发明实施例中控制力随桨叶方位角变化示意图。
作为一个实施例,步骤S4中所述基于CFD方法计算不同阻力片面积与控制力的关系,包括:
S101,读取桨叶几何参数、运动特点和旋翼所处流场特征划分旋翼桨叶的贴体网格和背景网格;
S102,对桨叶贴体网格和背景网格进行前处理,并对带有阻力片部分的桨叶进行网格加密;
S103,对桨叶网格在背景网格中进行挖洞和贡献单元搜索,并将得到的桨叶网格和背景网格的嵌套关系进行储存;
S104,初始化桨叶流场;
S105,更新所述桨叶网格坐标;
S106,更新所述桨叶网格和背景网格的嵌套关系;
S107,计算桨叶网格上的流场,更新背景网格上对应洞边界单元的流场信息。
作为一个实施例,步骤S5中所述改变控制力展向分布位置的改变范围是:0.6R至R,R是桨叶半径。
本发明实施例提供的一种基于可变面积阻力片的旋翼厚度噪声控制方法,通过在桨叶不同位置添加可变面积阻力片,可以针对性地对旋翼厚度噪声进行大幅度降噪,旋翼厚度噪声声压可降低6dB以上,噪声声压级降低50%以上。本发明实施例提供的方法属于主动控制方法,解决了被动控制方法中存在的不同飞行状态的适应性问题,实现了全包线飞行状态下的旋翼厚度噪声降噪。同时,相比现有的旋翼主动控制方法,实现方式简单有效、可显著降低降噪成本。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于设备实施例而言,由于其基本相似于方法实施例,所以描述得比较简单,相关之处参见方法实施例的部分说明即可。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (8)
1.一种基于可变面积阻力片的旋翼厚度噪声控制方法,其特征在于,包括以下步骤:
S1,读取旋翼尺寸和飞行参数,建立基于声学类比法的旋翼厚度噪声计算模型,计算给定飞行状态下旋翼厚度噪声P'T;
S2,在桨叶端部四分之一弦线处施加控制力F,建立基于声学类比法的控制源噪声计算模型,计算给定飞行状态下控制力F产生的控制源噪声p′F;
S3,根据旋翼厚度噪声降噪目标,设定抵消因子f,根据所述旋翼厚度噪声P'T和控制源噪声P'T,计算桨叶端部四分之一弦线处施加的控制力F;
S4,在桨叶端部四分之一弦线处添加可变面积阻力片,基于CFD方法计算不同阻力片面积与控制力的关系,得到桨叶端部四分之一弦线处阻力片的面积;
S5,改变控制力展向分布位置,根据所述不同阻力片面积与控制力的关系,计算不同位置的阻力片面积;
S6,根据步骤S5的结果,在桨叶不同位置添加相应面积的阻力片。
2.根据权利要求1所述的一种基于可变面积阻力片的旋翼厚度噪声控制方法,其特征在于:还包括,
S7,对带有相应面积阻力片的旋翼,根据步骤S1计算添加阻力片之后的旋翼厚度噪声,进行检验和完善。
3.根据权利要求1所述的一种基于可变面积阻力片的旋翼厚度噪声控制方法,其特征在于:
所述步骤S1中,选取桨盘平面内5R处的点为远场目标观测点,计算给定飞行状态下控制力旋翼的厚度噪声P'T;其中,R是桨叶半径。
4.根据权利要求1所述的一种基于可变面积阻力片的旋翼厚度噪声控制方法,其特征在于,所述控制源噪声p′F是:
<mrow>
<msubsup>
<mi>p</mi>
<mi>F</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>o</mi>
<mi>b</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>,</mo>
<mi>&tau;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mi>&Omega;</mi>
<mrow>
<mn>4</mn>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mi>d</mi>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>&psi;</mi>
</mrow>
</mfrac>
<mfrac>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>M</mi>
<mi>H</mi>
</msub>
<mi>sin</mi>
<mi>&psi;</mi>
</mrow>
<mo>)</mo>
</mrow>
</mfrac>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>M</mi>
<mi>H</mi>
</msub>
<mi>sin</mi>
<mi>&psi;</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>RM</mi>
<mi>H</mi>
</msub>
</mrow>
<mrow>
<mn>4</mn>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mi>d</mi>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>F</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>&psi;</mi>
</mrow>
</mfrac>
<mi>sin</mi>
<mi>&psi;</mi>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>M</mi>
<mi>H</mi>
</msub>
<mi>sin</mi>
<mi>&psi;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>F</mi>
<mi> </mi>
<mi>cos</mi>
<mi>&psi;</mi>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>M</mi>
<mi>H</mi>
</msub>
<mi>sin</mi>
<mi>&psi;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mfrac>
</mrow>
其中,MH为旋翼悬停时的桨尖马赫数,R是桨叶半径,Ψ是桨叶方位角,F为控制力,和分别表示F和ψ的偏导数,a0是声速,d是声源到远场目标观测点的距离,xobs是远场目标观测点的位置,τ是源时间,Ω是旋翼转速。
5.根据权利要求1所述的一种基于可变面积阻力片的旋翼厚度噪声控制方法,其特征在于,所述步骤S3远场目标观测点处的厚度噪声为:
f·p′T+p′F≈0
其中,抵消因子f=1表示声场完全对消,0<f<1代表部分抵消。
6.根据权利要求1所述的一种基于可变面积阻力片的旋翼厚度噪声控制方法,其特征在于,步骤S3中,所述计算桨叶端部四分之一弦线处施加的控制力F是:
<mrow>
<mi>F</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A</mi>
<mi>s</mi>
</msub>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
<msup>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mn>2</mn>
</msup>
<msup>
<msub>
<mi>M</mi>
<mi>H</mi>
</msub>
<mn>3</mn>
</msup>
<mi>cos</mi>
<mi>&psi;</mi>
</mrow>
<mrow>
<mn>60</mn>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>M</mi>
<mi>H</mi>
</msub>
<mi>sin</mi>
<mi>&psi;</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mrow>
<mn>50</mn>
<mo>+</mo>
<mn>39</mn>
<msup>
<msub>
<mi>M</mi>
<mi>H</mi>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>45</mn>
<msub>
<mi>M</mi>
<mi>H</mi>
</msub>
<mi>sin</mi>
<mi>&psi;</mi>
<mo>-</mo>
<mn>11</mn>
<msup>
<msub>
<mi>M</mi>
<mi>H</mi>
</msub>
<mn>2</mn>
</msup>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&psi;</mi>
<mo>+</mo>
<mn>12</mn>
<msup>
<msub>
<mi>M</mi>
<mi>H</mi>
</msub>
<mn>3</mn>
</msup>
<mi>sin</mi>
<mi>&psi;</mi>
<mo>-</mo>
<mn>18</mn>
<msup>
<msub>
<mi>M</mi>
<mi>H</mi>
</msub>
<mn>3</mn>
</msup>
<msup>
<mi>sin</mi>
<mn>3</mn>
</msup>
<mi>&psi;</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
其中,MH为旋翼悬停时的桨尖马赫数,As为桨盘面积,ρ0是空气密度,a0是声速,Ψ是桨叶方位角。
7.根据权利要求1所述的一种基于可变面积阻力片的旋翼厚度噪声控制方法,其特征在于,步骤S4中所述基于CFD方法计算不同阻力片面积与控制力的关系,包括:
S101,读取桨叶几何参数、运动特点和旋翼所处流场特征划分旋翼桨叶的贴体网格和背景网格;
S102,对桨叶贴体网格和背景网格进行前处理,并对带有阻力片部分的桨叶进行网格加密;
S103,对桨叶网格在背景网格中进行挖洞和贡献单元搜索,并将得到的桨叶网格和背景网格的嵌套关系进行储存;
S104,初始化桨叶流场;
S105,更新所述桨叶网格坐标;
S106,更新所述桨叶网格和背景网格的嵌套关系;
S107,计算桨叶网格上的流场,更新背景网格上对应洞边界单元的流场信息。
8.根据权利要求1所述的一种基于可变面积阻力片的旋翼厚度噪声控制方法,其特征在于,步骤S5中所述改变控制力展向分布位置的改变范围是:
0.6R至R,R是桨叶半径。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710307028.3A CN107220412B (zh) | 2017-05-04 | 2017-05-04 | 一种基于可变面积阻力片的旋翼厚度噪声控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710307028.3A CN107220412B (zh) | 2017-05-04 | 2017-05-04 | 一种基于可变面积阻力片的旋翼厚度噪声控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107220412A true CN107220412A (zh) | 2017-09-29 |
CN107220412B CN107220412B (zh) | 2020-06-23 |
Family
ID=59944822
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710307028.3A Active CN107220412B (zh) | 2017-05-04 | 2017-05-04 | 一种基于可变面积阻力片的旋翼厚度噪声控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107220412B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108021772A (zh) * | 2017-12-26 | 2018-05-11 | 南京航空航天大学 | 一种基于射流的降低旋翼噪声的方法及*** |
CN108170939A (zh) * | 2017-12-26 | 2018-06-15 | 南京航空航天大学 | 一种基于后缘襟翼的降低旋翼噪声的方法及*** |
CN108382606A (zh) * | 2018-02-23 | 2018-08-10 | 南京航空航天大学 | 一种用于降低直升机旋翼厚度噪声的实验装置 |
CN108446419A (zh) * | 2018-01-17 | 2018-08-24 | 天津大学 | 一种超声速边界层特征厚度估算方法 |
CN112214835A (zh) * | 2020-09-25 | 2021-01-12 | 中国直升机设计研究所 | 一种旋翼悬停状态气动噪声工程估算方法 |
CN113051666A (zh) * | 2021-03-25 | 2021-06-29 | 南京航空航天大学 | 一种旋翼飞行器噪声数字化分析方法及*** |
CN113449460A (zh) * | 2021-05-20 | 2021-09-28 | 北京理工大学 | 破片阻力系数的测量方法、装置及计算机可读存储介质 |
CN113670559A (zh) * | 2021-10-19 | 2021-11-19 | 中国空气动力研究与发展中心低速空气动力研究所 | 基于后缘小翼的直升机旋翼噪声主动控制风洞试验方法 |
CN113887090A (zh) * | 2021-11-18 | 2022-01-04 | 南京航空航天大学 | 一种旋翼气动噪声的计算方法及*** |
CN114275158A (zh) * | 2022-01-06 | 2022-04-05 | 南京航空航天大学 | 一种基于可收放阻力片机构的旋翼噪声控制装置 |
WO2023103168A1 (zh) * | 2021-12-06 | 2023-06-15 | 南京航空航天大学 | 一种融合声阵列和在桨控制的直升机主动噪声抑制装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015089404A2 (en) * | 2013-12-13 | 2015-06-18 | Lord Corporation | Redundant active vibration and noise control systems and methods |
CN105468818A (zh) * | 2015-11-12 | 2016-04-06 | 南京航空航天大学 | 一种用于直升机旋翼地面声场预测的方法 |
CN106114821A (zh) * | 2015-05-27 | 2016-11-16 | 深圳市高巨创新科技开发有限公司 | 一种低噪声飞行器螺旋桨的设计方法及螺旋桨构型 |
-
2017
- 2017-05-04 CN CN201710307028.3A patent/CN107220412B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015089404A2 (en) * | 2013-12-13 | 2015-06-18 | Lord Corporation | Redundant active vibration and noise control systems and methods |
CN106114821A (zh) * | 2015-05-27 | 2016-11-16 | 深圳市高巨创新科技开发有限公司 | 一种低噪声飞行器螺旋桨的设计方法及螺旋桨构型 |
CN105468818A (zh) * | 2015-11-12 | 2016-04-06 | 南京航空航天大学 | 一种用于直升机旋翼地面声场预测的方法 |
Non-Patent Citations (2)
Title |
---|
GAURAV GOPALAN 等: "Helicopter Thickness Noise Reduction Possibilities Through Active On-Blade Acoustic Control", 《JOURNAL OF AIRCRAFT》 * |
史勇杰 等: "桨叶气动外形对直升机桨_涡干扰噪声影响研究", 《南京航空航天大学学报》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108170939A (zh) * | 2017-12-26 | 2018-06-15 | 南京航空航天大学 | 一种基于后缘襟翼的降低旋翼噪声的方法及*** |
CN108021772A (zh) * | 2017-12-26 | 2018-05-11 | 南京航空航天大学 | 一种基于射流的降低旋翼噪声的方法及*** |
CN108446419B (zh) * | 2018-01-17 | 2022-02-22 | 天津大学 | 一种超声速边界层特征厚度估算方法 |
CN108446419A (zh) * | 2018-01-17 | 2018-08-24 | 天津大学 | 一种超声速边界层特征厚度估算方法 |
CN108382606A (zh) * | 2018-02-23 | 2018-08-10 | 南京航空航天大学 | 一种用于降低直升机旋翼厚度噪声的实验装置 |
CN112214835B (zh) * | 2020-09-25 | 2022-05-06 | 中国直升机设计研究所 | 一种旋翼悬停状态气动噪声工程估算方法 |
CN112214835A (zh) * | 2020-09-25 | 2021-01-12 | 中国直升机设计研究所 | 一种旋翼悬停状态气动噪声工程估算方法 |
CN113051666A (zh) * | 2021-03-25 | 2021-06-29 | 南京航空航天大学 | 一种旋翼飞行器噪声数字化分析方法及*** |
CN113051666B (zh) * | 2021-03-25 | 2022-05-31 | 南京航空航天大学 | 一种旋翼飞行器噪声数字化分析方法及*** |
CN113449460A (zh) * | 2021-05-20 | 2021-09-28 | 北京理工大学 | 破片阻力系数的测量方法、装置及计算机可读存储介质 |
CN113670559A (zh) * | 2021-10-19 | 2021-11-19 | 中国空气动力研究与发展中心低速空气动力研究所 | 基于后缘小翼的直升机旋翼噪声主动控制风洞试验方法 |
CN113670559B (zh) * | 2021-10-19 | 2022-01-25 | 中国空气动力研究与发展中心低速空气动力研究所 | 基于后缘小翼的直升机旋翼噪声主动控制风洞试验方法 |
CN113887090A (zh) * | 2021-11-18 | 2022-01-04 | 南京航空航天大学 | 一种旋翼气动噪声的计算方法及*** |
CN113887090B (zh) * | 2021-11-18 | 2024-01-30 | 南京航空航天大学 | 一种旋翼气动噪声的计算方法及*** |
WO2023103168A1 (zh) * | 2021-12-06 | 2023-06-15 | 南京航空航天大学 | 一种融合声阵列和在桨控制的直升机主动噪声抑制装置 |
JP7465039B2 (ja) | 2021-12-06 | 2024-04-10 | 南京航空航天大学 | 音響アレイ及びプロペラでの制御を統合したヘリコプター用能動騒音抑制装置 |
CN114275158A (zh) * | 2022-01-06 | 2022-04-05 | 南京航空航天大学 | 一种基于可收放阻力片机构的旋翼噪声控制装置 |
CN114275158B (zh) * | 2022-01-06 | 2023-01-03 | 南京航空航天大学 | 一种基于可收放阻力片机构的旋翼噪声控制装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107220412B (zh) | 2020-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107220412A (zh) | 一种基于可变面积阻力片的旋翼厚度噪声控制方法 | |
Glegg | The response of a swept blade row to a three-dimensional gust | |
Mangler et al. | The induced velocity field of a rotor | |
Gur et al. | Propeller performance at low advance ratio | |
CN108170939B (zh) | 一种基于后缘襟翼的降低旋翼噪声的方法及*** | |
CN103530482A (zh) | 一种非均匀入流中螺旋桨噪声数值预报方法 | |
Lee et al. | Numerical modeling of wind turbine aerodynamic noise in the time domain | |
Huang et al. | Helicopter rotor noise prediction using a convected FW-H equation in the frequency domain | |
CN104192302A (zh) | 基于绕尖头冯卡门曲线回转体基准流场的乘波体设计方法 | |
US4168939A (en) | Acoustically swept rotor | |
Chaviaropoulos | Flap/lead–lag aeroelastic stability of wind turbine blade sections | |
CN108021772A (zh) | 一种基于射流的降低旋翼噪声的方法及*** | |
Kopiev et al. | Investigation of airframe noise for a large-scale wing model with high-lift devices | |
Stuermer et al. | Progress in aerodynamic and aeroacoustic integration of CROR propulsion systems | |
CN109747818A (zh) | 基于桨尖质量射流的直升机旋翼气动干扰控制方法 | |
Li et al. | Extensions and Applications of Lyu and Ayton's Serrated Trailing-Edge Noise Model to Rotorcraft | |
Spurlock et al. | Cartesian mesh simulations for the third AIAA sonic boom prediction workshop | |
Stoll | Design of quiet UAV propellers | |
Roger et al. | Vortex-shedding noise and potential-interaction noise modeling by a reversed sears' problem | |
Nakamura | Prediction of blade-vortex interaction noise from measured blade pressure | |
Mao et al. | Computation of instantaneous and time-averaged active acoustic intensity field around rotating source | |
Xu et al. | Investigation on improving efficiency of high-altitude propeller with tandem configuration | |
Cramer et al. | Investigation of Anti-Phase Asymmetric Quiet Rotor Technology | |
Arun Prakash et al. | Experimental investigation of stepped aerofoil using propeller test rig | |
CN206367596U (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 |