CN111046493A - 一种基于集中载荷的旋翼气动噪声计算方法 - Google Patents

一种基于集中载荷的旋翼气动噪声计算方法 Download PDF

Info

Publication number
CN111046493A
CN111046493A CN201911232693.6A CN201911232693A CN111046493A CN 111046493 A CN111046493 A CN 111046493A CN 201911232693 A CN201911232693 A CN 201911232693A CN 111046493 A CN111046493 A CN 111046493A
Authority
CN
China
Prior art keywords
blade
micro
segment
lift
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.)
Pending
Application number
CN201911232693.6A
Other languages
English (en)
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 CN201911232693.6A priority Critical patent/CN111046493A/zh
Publication of CN111046493A publication Critical patent/CN111046493A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Structures Of Non-Positive Displacement Pumps (AREA)

Abstract

本发明属于直升机空气动力学技术领域,公开了一种基于集中载荷的旋翼气动噪声计算方法。包括:S1,进行旋翼气动仿真得到多个桨叶剖面的集中载荷,每个桨叶剖面的集中载荷包括对应桨叶剖面的升力系数、阻力系数;S2,根据每个桨叶剖面的升力系数、阻力系数,计算每个桨叶微段的升力、阻力;其中,桨叶微段与桨叶剖面一一对应;S3,根据每个桨叶微段的升力、阻力分别计算对应桨叶微段的升力、阻力的离散载荷;S4,根据所有桨叶微段的升力、阻力的离散载荷计算旋翼气动噪声,将旋翼气动计算得到的集中载荷转换成沿弦向分布的离散载荷,高效率地计算出各个观测点处的旋翼噪声值。

Description

一种基于集中载荷的旋翼气动噪声计算方法
技术领域
本发明属于直升机空气动力学技术领域,尤其涉及一种基于集中载荷的旋翼气动噪声计算方法。
背景技术
外部噪声大是直升机的固有缺点之一,使得直升机会对周围环境以及地面人员产生较严重的噪声污染和影响,不论在民用领域还是军事领域,这都限制了直升机的进一步使用。而旋翼和尾桨的气动噪声是直升机远场噪声最主要的来源,因此,开展直升机旋翼/尾桨气动噪声研究,建立可靠的噪声预估算法,是实现直升机降噪技术的基础。
直升机旋翼气动噪声计算领域常用的方法是固体边界发声的气动声学基本方程,即FW-H方程,它的实用形式是Farassat 1A公式:
Figure BDA0002303120440000011
Figure BDA0002303120440000012
Figure BDA0002303120440000013
Figure BDA0002303120440000014
Figure BDA0002303120440000015
式中,p’是扰动声压,p’T代表旋翼旋转噪声中的厚度噪声成分,p’L代表旋翼旋转噪声中的载荷噪声成分;a0和ρ0分别表示无扰动介质中的音速和密度,vn为桨叶表面的法向运动速度;r为积分面到观察点之间的距离,
Figure BDA0002303120440000016
为积分面指向观察点的单位矢量,[·]ret表示括号内的变量都在延迟时间下取值;桨叶载荷定义为
Figure BDA0002303120440000021
Pij是压缩应力张量,其包含表面压强和粘性应力;
Figure BDA0002303120440000022
和t分别表示观测点所在位置和声波到达观测点的时刻;字母头上圆点“·”表示时间导数,n表示积分面的单位法向矢量,M为积分面运动马赫数。
上述噪声计算公式中,需要确定lr等参数。这些是与桨叶气动载荷相关的量,在计算旋翼噪声之前,应确定桨叶表面的气动载荷分布。通过旋翼流场计算可以得到桨叶的载荷分布,常用的方法之一是基于RANS方程的CFD(计算流体力学)方法,该方法将桨叶离散成一个个网格,通过求解N-S方程得到每个网格点处的载荷。相比于CFD方法,自由尾迹方法计算效率高,并且能模拟出桨尖涡轨迹,特别适合于旋翼桨-涡干扰状态的计算。它输出的是桨叶沿径向分布剖面的集中载荷,还不能直接用于噪声计算。
发明内容
针对上述背景技术中的问题,本发明的目的在于提供一种基于集中载荷的旋翼气动噪声计算方法,根据非紧致源假设,将旋翼气动计算得到的集中载荷转换成沿弦向分布的离散载荷,高效率地计算出各个观测点处的旋翼噪声值。
为达到上述目的,本发明采用如下技术方案予以实现。
一种基于集中载荷的旋翼气动噪声计算方法,所述方法包括:
S1,进行旋翼气动仿真得到多个桨叶剖面的集中载荷,每个桨叶剖面的集中载荷包括对应桨叶剖面的升力系数、阻力系数;
S2,根据每个桨叶剖面的升力系数、阻力系数,计算每个桨叶微段的升力、阻力;其中,桨叶微段与桨叶剖面一一对应;
S3,根据每个桨叶微段的升力、阻力分别计算对应桨叶微段的升力、阻力的离散载荷;
S4,根据所有桨叶微段的升力、阻力的离散载荷计算旋翼气动噪声。
本发明技术方案的特点和进一步的改进为:
(1)S1中,采用自由尾迹方法获取多个桨叶剖面的集中载荷。
(2)S2中,采用叶素理论计算每个桨叶微段的升力、阻力。
(3)S2具体为:
对于展向位置为r,宽度为dr的微段,得到该桨叶微段的升力dL(ψ,r)为:
Figure BDA0002303120440000031
式中,Cl(ψ,r)表示桨叶剖面的升力系数,Cd(ψ,r)表示桨叶剖面的阻力系数,c表示桨叶微段的弦长,U表示桨叶微段的来流速度,ρ表示空气密度,ψ表示桨叶的方位角。
(4)S2具体为:
对于展向位置为r,宽度为dr的微段,得到该桨叶微段的阻力dD(ψ,r)为:
Figure BDA0002303120440000032
式中,Cl(ψ,r)表示桨叶剖面的升力系数,Cd(ψ,r)表示桨叶剖面的阻力系数,c表示桨叶微段的弦长,U表示桨叶微段的来流速度,ρ表示空气密度,ψ表示桨叶的方位角。
(5)S3,具体为:
在非紧致源假设下,根据每个桨叶微段的升力计算对应桨叶微段的升力的离散载荷
Figure BDA0002303120440000033
为:
Figure BDA0002303120440000034
式中,x为沿弦向的位置,c表示桨叶微段的弦长。
(6)S3,具体为:
在非紧致源假设下,根据每个桨叶微段的升力计算对应桨叶微段的阻力的离散载荷
Figure BDA0002303120440000035
为:
Figure BDA0002303120440000036
式中,x为沿弦向的位置,c表示桨叶微段的弦长。
(7)S4中,根据FW-H方程推导的Farassat 1A公式,采用所有桨叶微段的升、阻力的离散载荷计算旋翼气动噪声。
本发明的有益效果:本发明技术方案基于集中载荷计算旋翼气动噪声的方法,可以利用自由尾迹模拟的旋翼气动载荷结果作为输入,计算效率高,而且非常适用于计算旋翼的桨-涡干扰噪声。
附图说明
图1为本发明实施例提供的旋翼气动计算结果示意图;
图2为本发明实施例提供的旋翼噪声评估结果示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明中,采用旋翼气动仿真得到的集中载荷作为噪声计算输入。利用自由尾迹方法计算得到桨叶剖面升力系数、阻力系数和来流马赫数,对于展向位置为r,宽度为dr的叶素,由叶素理论可知,桨叶微段升、阻力分别为:
Figure BDA0002303120440000041
Figure BDA0002303120440000042
式中,Cl(ψ,r),Cd(ψ,r)分别表示桨叶升力系数和阻力系数,c表示桨叶该段弦长,U表示桨叶微段来流速度。
利用Farassat 1A公式进行气动噪声计算需要对桨叶的弦向和展向展开积分计算,而桨叶的集中载荷位于在桨叶弦向1/4点处,因此,需将桨叶剖面集中载荷沿弦向进行离散分布。在非紧致源假设下,桨叶载荷离散公式为:
Figure BDA0002303120440000043
Figure BDA0002303120440000044
已知桨叶沿径向的集中载荷,通过公式(6),就可以得到每个剖面载荷沿弦向的分布,从而进行旋翼气动噪声计算。
本发明技术方案的基于集中载荷计算旋翼气动噪声算法,可以利用自由尾迹模拟的旋翼气动载荷结果作为输入,计算效率高,而且非常适用于计算旋翼的桨-涡干扰噪声。图1和图2展示了某旋翼斜下降状态的计算结果,斜下降状态是典型的桨-涡干扰状态,该飞行状态下会发生平行桨-涡干扰,从而产生严重的桨-涡干扰噪声。从图1可以看出,剖面的法向力系数在90°方位角(前行侧)和270°方位角(后行侧)附近剧烈振荡,说明发生了桨-涡干扰。图2是利用该算法得到的旋翼噪声分布云图,可以看出在前行侧和后行侧斜下方各有一处噪声集中区域,为桨-涡干扰噪声,与气动计算结果一致。
以上所述,仅为本发明的具体实施例,对本发明进行详细描述,未详尽部分为常规技术。但本发明的保护范围不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。本发明的保护范围应以所述权利要求的保护范围为准。

Claims (8)

1.一种基于集中载荷的旋翼气动噪声计算方法,其特征在于,所述方法包括:
S1,进行旋翼气动仿真得到多个桨叶剖面的集中载荷,每个桨叶剖面的集中载荷包括对应桨叶剖面的升力系数、阻力系数;
S2,根据每个桨叶剖面的升力系数、阻力系数,计算每个桨叶微段的升力、阻力;其中,桨叶微段与桨叶剖面一一对应;
S3,根据每个桨叶微段的升力、阻力分别计算对应桨叶微段的升力、阻力的离散载荷;
S4,根据所有桨叶微段的升力、阻力的离散载荷计算旋翼气动噪声。
2.根据权利要求1所述的一种基于集中载荷的旋翼气动噪声计算方法,其特征在于,S1中,采用自由尾迹方法获取多个桨叶剖面的集中载荷。
3.根据权利要求1所述的一种基于集中载荷的旋翼气动噪声计算方法,其特征在于,S2中,采用叶素理论计算每个桨叶微段的升力、阻力。
4.根据权利要求1所述的一种基于集中载荷的旋翼气动噪声计算方法,其特征在于,S2具体为:
对于展向位置为r,宽度为dr的微段,得到该桨叶微段的升力dL(ψ,r)为:
Figure FDA0002303120430000011
式中,Cl(ψ,r)表示桨叶剖面的升力系数,Cd(ψ,r)表示桨叶剖面的阻力系数,c表示桨叶微段的弦长,U表示桨叶微段的来流速度,ρ表示空气密度,ψ表示桨叶的方位角。
5.根据权利要求1所述的一种基于集中载荷的旋翼气动噪声计算方法,其特征在于,S2具体为:
对于展向位置为r,宽度为dr的微段,得到该桨叶微段的阻力dD(ψ,r)为:
Figure FDA0002303120430000021
式中,Cl(ψ,r)表示桨叶剖面的升力系数,Cd(ψ,r)表示桨叶剖面的阻力系数,c表示桨叶微段的弦长,U表示桨叶微段的来流速度,ρ表示空气密度,ψ表示桨叶的方位角。
6.根据权利要求5所述的一种基于集中载荷的旋翼气动噪声计算方法,其特征在于,S3,具体为:
在非紧致源假设下,根据每个桨叶微段的升力计算对应桨叶微段的升力的离散载荷
Figure FDA0002303120430000022
为:
Figure FDA0002303120430000023
式中,x为沿弦向的位置,c表示桨叶微段的弦长。
7.根据权利要求5所述的一种基于集中载荷的旋翼气动噪声计算方法,其特征在于,S3,具体为:
在非紧致源假设下,根据每个桨叶微段的升力计算对应桨叶微段的阻力的离散载荷
Figure FDA0002303120430000024
为:
Figure FDA0002303120430000025
式中,x为沿弦向的位置,c表示桨叶微段的弦长。
8.根据权利要求1所述的一种基于集中载荷的旋翼气动噪声计算方法,其特征在于,S4中,根据FW-H方程推导的Farassat 1A公式,采用所有桨叶微段的升、阻力的离散载荷计算旋翼气动噪声。
CN201911232693.6A 2019-12-04 2019-12-04 一种基于集中载荷的旋翼气动噪声计算方法 Pending CN111046493A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911232693.6A CN111046493A (zh) 2019-12-04 2019-12-04 一种基于集中载荷的旋翼气动噪声计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911232693.6A CN111046493A (zh) 2019-12-04 2019-12-04 一种基于集中载荷的旋翼气动噪声计算方法

Publications (1)

Publication Number Publication Date
CN111046493A true CN111046493A (zh) 2020-04-21

Family

ID=70234929

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911232693.6A Pending CN111046493A (zh) 2019-12-04 2019-12-04 一种基于集中载荷的旋翼气动噪声计算方法

Country Status (1)

Country Link
CN (1) CN111046493A (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112145371A (zh) * 2020-09-08 2020-12-29 扬州大学 一种风力机气动噪声快速预测方法
CN112173075A (zh) * 2020-09-25 2021-01-05 中国直升机设计研究所 一种直升机低噪声旋翼桨叶气动外形
CN112199772A (zh) * 2020-09-25 2021-01-08 中国直升机设计研究所 一种直升机旋翼气动噪声计算方法
CN112487730A (zh) * 2020-10-30 2021-03-12 南京航空航天大学 一种基于相角控制的多旋翼飞行器噪声抑制方法
CN112504426A (zh) * 2020-11-20 2021-03-16 中国直升机设计研究所 一种基于峰值搜索的旋翼桨涡干扰噪声整周期平均方法
CN113051668A (zh) * 2021-04-14 2021-06-29 长沙神弓信息科技有限公司 一种旋翼类无人机等效阻力系数的测算方法
CN113051666A (zh) * 2021-03-25 2021-06-29 南京航空航天大学 一种旋翼飞行器噪声数字化分析方法及系统
CN113887090A (zh) * 2021-11-18 2022-01-04 南京航空航天大学 一种旋翼气动噪声的计算方法及系统
CN115339619A (zh) * 2021-05-14 2022-11-15 国家直升机中心米尔&卡莫夫股份公司 飞行器支承构件空气动力剖面
CN116086599A (zh) * 2023-04-10 2023-05-09 中国飞行试验研究院 航空器低空辐射声场获取方法、装置、设备及存储介质
CN112487730B (zh) * 2020-10-30 2024-05-28 南京航空航天大学 一种基于相角控制的多旋翼飞行器噪声抑制方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106275423A (zh) * 2016-08-23 2017-01-04 成都翼高九天科技有限公司 一种新型无人机旋翼空气动力学结构
CN108021772A (zh) * 2017-12-26 2018-05-11 南京航空航天大学 一种基于射流的降低旋翼噪声的方法及系统
CN110466751A (zh) * 2019-09-03 2019-11-19 中国科学院合肥物质科学研究院 一种低噪声的高效便携式无人机旋翼结构及其设计方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106275423A (zh) * 2016-08-23 2017-01-04 成都翼高九天科技有限公司 一种新型无人机旋翼空气动力学结构
CN108021772A (zh) * 2017-12-26 2018-05-11 南京航空航天大学 一种基于射流的降低旋翼噪声的方法及系统
CN110466751A (zh) * 2019-09-03 2019-11-19 中国科学院合肥物质科学研究院 一种低噪声的高效便携式无人机旋翼结构及其设计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
宋辰瑶: "基于尾迹分析的旋翼旋转噪声计算及桨-涡干扰噪声研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112145371A (zh) * 2020-09-08 2020-12-29 扬州大学 一种风力机气动噪声快速预测方法
CN112173075A (zh) * 2020-09-25 2021-01-05 中国直升机设计研究所 一种直升机低噪声旋翼桨叶气动外形
CN112199772A (zh) * 2020-09-25 2021-01-08 中国直升机设计研究所 一种直升机旋翼气动噪声计算方法
CN112173075B (zh) * 2020-09-25 2022-12-30 中国直升机设计研究所 一种直升机低噪声旋翼桨叶气动外形
CN112199772B (zh) * 2020-09-25 2022-11-04 中国直升机设计研究所 一种直升机旋翼气动噪声计算方法
CN112487730A (zh) * 2020-10-30 2021-03-12 南京航空航天大学 一种基于相角控制的多旋翼飞行器噪声抑制方法
CN112487730B (zh) * 2020-10-30 2024-05-28 南京航空航天大学 一种基于相角控制的多旋翼飞行器噪声抑制方法
CN112504426A (zh) * 2020-11-20 2021-03-16 中国直升机设计研究所 一种基于峰值搜索的旋翼桨涡干扰噪声整周期平均方法
CN113051666B (zh) * 2021-03-25 2022-05-31 南京航空航天大学 一种旋翼飞行器噪声数字化分析方法及系统
CN113051666A (zh) * 2021-03-25 2021-06-29 南京航空航天大学 一种旋翼飞行器噪声数字化分析方法及系统
CN113051668A (zh) * 2021-04-14 2021-06-29 长沙神弓信息科技有限公司 一种旋翼类无人机等效阻力系数的测算方法
CN115339619A (zh) * 2021-05-14 2022-11-15 国家直升机中心米尔&卡莫夫股份公司 飞行器支承构件空气动力剖面
CN113887090A (zh) * 2021-11-18 2022-01-04 南京航空航天大学 一种旋翼气动噪声的计算方法及系统
CN113887090B (zh) * 2021-11-18 2024-01-30 南京航空航天大学 一种旋翼气动噪声的计算方法及系统
CN116086599A (zh) * 2023-04-10 2023-05-09 中国飞行试验研究院 航空器低空辐射声场获取方法、装置、设备及存储介质
CN116086599B (zh) * 2023-04-10 2023-10-20 中国飞行试验研究院 航空器低空辐射声场获取方法、装置、设备及存储介质

Similar Documents

Publication Publication Date Title
CN111046493A (zh) 一种基于集中载荷的旋翼气动噪声计算方法
Baltazar et al. Open-water thrust and torque predictions of a ducted propeller system with a panel method
Lee et al. Numerical investigation of the aerodynamics and wake structures of horizontal axis wind turbines by using nonlinear vortex lattice method
Hall et al. Power requirements for large-amplitude flapping flight
Jain et al. An assessment of RCAS performance prediction for conventional and advanced rotor configurations
Gupta et al. CFD analysis of wind turbine blade with winglets
Carroll et al. Comparison of a blade element momentum model to 3D CFD simulations for small scale propellers
Park et al. Loosely coupled multibody dynamics–CFD analysis for a rotor in descending flight
CN110077575A (zh) 飞机系统和方法
Singh et al. Dynamic stall modeling using viscous vortex particle method for coaxial rotors
CN113051666A (zh) 一种旋翼飞行器噪声数字化分析方法及系统
Rajakumar et al. Computational fluid dynamics of wind turbine blade at various angles of attack and low Reynolds number
Lind et al. Ducted rotor performance calculations and comparison with experimental data
Stuermer et al. Aerodynamic and aeroacoustic analysis of contra-rotating open rotor propulsion systems at low-speed flight conditions
Agostinelli et al. Propeller-Wing Interaction using Rapid Computational Methods
CN112283051B (zh) 一种基于升力线模型的振动信号特征优化方法及系统
Wie et al. Helicopter rotor noise in the merged tip-vortex and blade interaction condition
Balakrishnan et al. Effects of different thickness of UAV airfoil on aerodynamics performance using CFD techniques
Peng et al. A novel composite calculation model for power coefficient and flapping moment coefficient of wind turbine
Tuncer et al. Computational study of subsonic flow over a delta canard-wing-body configuration
Filotas Vortex induced helicopter blade loads and noise
Rozehnal et al. Hysteresis assumptions for propeller in yaw in close stall conditions
Bain et al. Computational modeling of variable-droop leading edge in forward flight
Tejero E et al. Numerical simulation of the tip aerodynamics and acoustics test
Jayasundara et al. CFD and Aeroacoustic Analysis of Wingtip-Mounted Propellers

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20200421

RJ01 Rejection of invention patent application after publication