CN111046493A - 一种基于集中载荷的旋翼气动噪声计算方法 - Google Patents
一种基于集中载荷的旋翼气动噪声计算方法 Download PDFInfo
- 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
Links
Images
Classifications
-
- 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
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本发明属于直升机空气动力学技术领域,公开了一种基于集中载荷的旋翼气动噪声计算方法。包括:S1,进行旋翼气动仿真得到多个桨叶剖面的集中载荷,每个桨叶剖面的集中载荷包括对应桨叶剖面的升力系数、阻力系数;S2,根据每个桨叶剖面的升力系数、阻力系数,计算每个桨叶微段的升力、阻力;其中,桨叶微段与桨叶剖面一一对应;S3,根据每个桨叶微段的升力、阻力分别计算对应桨叶微段的升力、阻力的离散载荷;S4,根据所有桨叶微段的升力、阻力的离散载荷计算旋翼气动噪声,将旋翼气动计算得到的集中载荷转换成沿弦向分布的离散载荷,高效率地计算出各个观测点处的旋翼噪声值。
Description
技术领域
本发明属于直升机空气动力学技术领域,尤其涉及一种基于集中载荷的旋翼气动噪声计算方法。
背景技术
外部噪声大是直升机的固有缺点之一,使得直升机会对周围环境以及地面人员产生较严重的噪声污染和影响,不论在民用领域还是军事领域,这都限制了直升机的进一步使用。而旋翼和尾桨的气动噪声是直升机远场噪声最主要的来源,因此,开展直升机旋翼/尾桨气动噪声研究,建立可靠的噪声预估算法,是实现直升机降噪技术的基础。
直升机旋翼气动噪声计算领域常用的方法是固体边界发声的气动声学基本方程,即FW-H方程,它的实用形式是Farassat 1A公式:
式中,p’是扰动声压,p’T代表旋翼旋转噪声中的厚度噪声成分,p’L代表旋翼旋转噪声中的载荷噪声成分;a0和ρ0分别表示无扰动介质中的音速和密度,vn为桨叶表面的法向运动速度;r为积分面到观察点之间的距离,为积分面指向观察点的单位矢量,[·]ret表示括号内的变量都在延迟时间下取值;桨叶载荷定义为Pij是压缩应力张量,其包含表面压强和粘性应力;和t分别表示观测点所在位置和声波到达观测点的时刻;字母头上圆点“·”表示时间导数,n表示积分面的单位法向矢量,M为积分面运动马赫数。
上述噪声计算公式中,需要确定lr等参数。这些是与桨叶气动载荷相关的量,在计算旋翼噪声之前,应确定桨叶表面的气动载荷分布。通过旋翼流场计算可以得到桨叶的载荷分布,常用的方法之一是基于RANS方程的CFD(计算流体力学)方法,该方法将桨叶离散成一个个网格,通过求解N-S方程得到每个网格点处的载荷。相比于CFD方法,自由尾迹方法计算效率高,并且能模拟出桨尖涡轨迹,特别适合于旋翼桨-涡干扰状态的计算。它输出的是桨叶沿径向分布剖面的集中载荷,还不能直接用于噪声计算。
发明内容
针对上述背景技术中的问题,本发明的目的在于提供一种基于集中载荷的旋翼气动噪声计算方法,根据非紧致源假设,将旋翼气动计算得到的集中载荷转换成沿弦向分布的离散载荷,高效率地计算出各个观测点处的旋翼噪声值。
为达到上述目的,本发明采用如下技术方案予以实现。
一种基于集中载荷的旋翼气动噪声计算方法,所述方法包括:
S1,进行旋翼气动仿真得到多个桨叶剖面的集中载荷,每个桨叶剖面的集中载荷包括对应桨叶剖面的升力系数、阻力系数;
S2,根据每个桨叶剖面的升力系数、阻力系数,计算每个桨叶微段的升力、阻力;其中,桨叶微段与桨叶剖面一一对应;
S3,根据每个桨叶微段的升力、阻力分别计算对应桨叶微段的升力、阻力的离散载荷;
S4,根据所有桨叶微段的升力、阻力的离散载荷计算旋翼气动噪声。
本发明技术方案的特点和进一步的改进为:
(1)S1中,采用自由尾迹方法获取多个桨叶剖面的集中载荷。
(2)S2中,采用叶素理论计算每个桨叶微段的升力、阻力。
(3)S2具体为:
对于展向位置为r,宽度为dr的微段,得到该桨叶微段的升力dL(ψ,r)为:
式中,Cl(ψ,r)表示桨叶剖面的升力系数,Cd(ψ,r)表示桨叶剖面的阻力系数,c表示桨叶微段的弦长,U表示桨叶微段的来流速度,ρ表示空气密度,ψ表示桨叶的方位角。
(4)S2具体为:
对于展向位置为r,宽度为dr的微段,得到该桨叶微段的阻力dD(ψ,r)为:
式中,Cl(ψ,r)表示桨叶剖面的升力系数,Cd(ψ,r)表示桨叶剖面的阻力系数,c表示桨叶微段的弦长,U表示桨叶微段的来流速度,ρ表示空气密度,ψ表示桨叶的方位角。
(5)S3,具体为:
式中,x为沿弦向的位置,c表示桨叶微段的弦长。
(6)S3,具体为:
式中,x为沿弦向的位置,c表示桨叶微段的弦长。
(7)S4中,根据FW-H方程推导的Farassat 1A公式,采用所有桨叶微段的升、阻力的离散载荷计算旋翼气动噪声。
本发明的有益效果:本发明技术方案基于集中载荷计算旋翼气动噪声的方法,可以利用自由尾迹模拟的旋翼气动载荷结果作为输入,计算效率高,而且非常适用于计算旋翼的桨-涡干扰噪声。
附图说明
图1为本发明实施例提供的旋翼气动计算结果示意图;
图2为本发明实施例提供的旋翼噪声评估结果示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明中,采用旋翼气动仿真得到的集中载荷作为噪声计算输入。利用自由尾迹方法计算得到桨叶剖面升力系数、阻力系数和来流马赫数,对于展向位置为r,宽度为dr的叶素,由叶素理论可知,桨叶微段升、阻力分别为:
式中,Cl(ψ,r),Cd(ψ,r)分别表示桨叶升力系数和阻力系数,c表示桨叶该段弦长,U表示桨叶微段来流速度。
利用Farassat 1A公式进行气动噪声计算需要对桨叶的弦向和展向展开积分计算,而桨叶的集中载荷位于在桨叶弦向1/4点处,因此,需将桨叶剖面集中载荷沿弦向进行离散分布。在非紧致源假设下,桨叶载荷离散公式为:
已知桨叶沿径向的集中载荷,通过公式(6),就可以得到每个剖面载荷沿弦向的分布,从而进行旋翼气动噪声计算。
本发明技术方案的基于集中载荷计算旋翼气动噪声算法,可以利用自由尾迹模拟的旋翼气动载荷结果作为输入,计算效率高,而且非常适用于计算旋翼的桨-涡干扰噪声。图1和图2展示了某旋翼斜下降状态的计算结果,斜下降状态是典型的桨-涡干扰状态,该飞行状态下会发生平行桨-涡干扰,从而产生严重的桨-涡干扰噪声。从图1可以看出,剖面的法向力系数在90°方位角(前行侧)和270°方位角(后行侧)附近剧烈振荡,说明发生了桨-涡干扰。图2是利用该算法得到的旋翼噪声分布云图,可以看出在前行侧和后行侧斜下方各有一处噪声集中区域,为桨-涡干扰噪声,与气动计算结果一致。
以上所述,仅为本发明的具体实施例,对本发明进行详细描述,未详尽部分为常规技术。但本发明的保护范围不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。本发明的保护范围应以所述权利要求的保护范围为准。
Claims (8)
1.一种基于集中载荷的旋翼气动噪声计算方法,其特征在于,所述方法包括:
S1,进行旋翼气动仿真得到多个桨叶剖面的集中载荷,每个桨叶剖面的集中载荷包括对应桨叶剖面的升力系数、阻力系数;
S2,根据每个桨叶剖面的升力系数、阻力系数,计算每个桨叶微段的升力、阻力;其中,桨叶微段与桨叶剖面一一对应;
S3,根据每个桨叶微段的升力、阻力分别计算对应桨叶微段的升力、阻力的离散载荷;
S4,根据所有桨叶微段的升力、阻力的离散载荷计算旋翼气动噪声。
2.根据权利要求1所述的一种基于集中载荷的旋翼气动噪声计算方法,其特征在于,S1中,采用自由尾迹方法获取多个桨叶剖面的集中载荷。
3.根据权利要求1所述的一种基于集中载荷的旋翼气动噪声计算方法,其特征在于,S2中,采用叶素理论计算每个桨叶微段的升力、阻力。
8.根据权利要求1所述的一种基于集中载荷的旋翼气动噪声计算方法,其特征在于,S4中,根据FW-H方程推导的Farassat 1A公式,采用所有桨叶微段的升、阻力的离散载荷计算旋翼气动噪声。
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 (10)
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 | 中国直升机设计研究所 | 一种基于峰值搜索的旋翼桨涡干扰噪声整周期平均方法 |
CN113051666A (zh) * | 2021-03-25 | 2021-06-29 | 南京航空航天大学 | 一种旋翼飞行器噪声数字化分析方法及系统 |
CN113051668A (zh) * | 2021-04-14 | 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 | 中国飞行试验研究院 | 航空器低空辐射声场获取方法、装置、设备及存储介质 |
Citations (3)
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 | 中国科学院合肥物质科学研究院 | 一种低噪声的高效便携式无人机旋翼结构及其设计方法 |
-
2019
- 2019-12-04 CN CN201911232693.6A patent/CN111046493A/zh active Pending
Patent Citations (3)
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)
Title |
---|
宋辰瑶: "基于尾迹分析的旋翼旋转噪声计算及桨-涡干扰噪声研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (16)
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) | 一种基于集中载荷的旋翼气动噪声计算方法 | |
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 | |
Li et al. | Prediction of rotorcraft broadband trailing-edge noise and parameter sensitivity study | |
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 | |
CN113051666A (zh) | 一种旋翼飞行器噪声数字化分析方法及系统 | |
CN110077575A (zh) | 飞机系统和方法 | |
Carroll et al. | Comparison of a blade element momentum model to 3D CFD simulations for small scale propellers | |
Lind et al. | Ducted rotor performance calculations and comparison with experimental data | |
Singh et al. | Dynamic stall modeling using viscous vortex particle method for coaxial rotors | |
Rajakumar et al. | Computational fluid dynamics of wind turbine blade at various angles of attack and low Reynolds number | |
Wu et al. | Using the combined flow control accessory to the aerodynamic performance enhancement of bio-inspired seagull airfoils | |
Agostinelli et al. | Propeller-Wing Interaction using Rapid Computational Methods | |
Peng et al. | A novel composite calculation model for power coefficient and flapping moment coefficient of wind turbine | |
Stuermer et al. | Aerodynamic and aeroacoustic analysis of contra-rotating open rotor propulsion systems at low-speed flight conditions | |
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 | |
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 | |
Tejero E et al. | Numerical simulation of the tip aerodynamics and acoustics test | |
Butt et al. | Numerical Anslysis of Cropped Delta Wing with Reflex Camber Cross-Section at Transonic Speed with Leading-Edge Tubercles | |
Rasam et al. | Aerodynamic noise prediction for a wind turbine using numerical flow simulations and semi-empirical modelling approaches |
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 |