CN112199772A - 一种直升机旋翼气动噪声计算方法 - Google Patents
一种直升机旋翼气动噪声计算方法 Download PDFInfo
- Publication number
- CN112199772A CN112199772A CN202011020831.7A CN202011020831A CN112199772A CN 112199772 A CN112199772 A CN 112199772A CN 202011020831 A CN202011020831 A CN 202011020831A CN 112199772 A CN112199772 A CN 112199772A
- Authority
- CN
- China
- Prior art keywords
- blade
- rotor
- dimensional
- spanwise
- chord length
- 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
Images
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
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Turbine Rotor Nozzle Sealing (AREA)
Abstract
本发明属于直升机旋翼气动噪声分析与控制技术领域,具体涉及一种直升机旋翼气动噪声快速计算方法。基于刚性桨叶假设,将旋翼桨叶分解为二维翼型和一维梁模型;首先生成二维翼型网格,然后沿展向进行插值,最终得到三维桨叶表面网格;流场输入参数代入旋翼飞行器综合性分析软件计算旋翼展向分布的截面升力;噪声计算输入参数、桨叶表面网格信息和展向截面升力用来求解旋翼噪声声场。
Description
技术领域
本发明属于升机旋翼气动噪声分析与控制技术领域,具体涉及一种直升机旋翼气动噪声计算方法。
背景技术
直升机旋翼产生的气动噪声主要由厚度噪声和载荷噪声两项叠加,厚度噪声取决于桨叶的几何结构和运动,载荷噪声与桨叶表面的压力分布有关。目前直升机旋翼气动噪声的快速计算方法根据几何函数来拟合桨叶外形,生成的桨叶表面网格与真实桨叶几何结构存在较大差异,导致噪声计算存在误差。
发明内容
针对背景技术里存在的问题,本发明提供了如下方案:
一种直升机旋翼气动噪声计算方法,所述方法包括:
第一步,采用翼型生成软件生成二维翼型表面网格;
第二步,将生成的二维翼型表面网格根据弦长分布规律进行缩放得到沿桨叶展向各剖面的网格;
第三步,对桨根和桨尖部分进行处理以生成O型桨尖;
第四步,根据桨叶展向各个位置的剖面二维翼型表面网格生成旋翼桨叶三维桨叶表面网格;
第五步,流场输入参数代入旋翼飞行器综合性分析软件计算旋翼展向分布的截面升力;
第六步,噪声计算输入参数、桨叶表面网格信息和展向截面升力求解旋翼噪声声场。
进一步,所述第一步中,还包括:生成沿桨叶展向分布的一维站位点。
进一步,所述第二步,将生成的桨叶二维翼型表面网格按照桨叶展向一维站位点分布规律进行插值并根据弦长分布规律进行缩放以得到沿桨叶展向各剖面的网格。
进一步,所述第三步中,对桨根和桨尖部分进行翻折处理以生成O型翻折网格。
进一步,所述第四步中,根据旋翼桨叶的扭转角和1/4弦长位置沿桨叶展向的分布规律,分别对桨叶展向各个位置的剖面二维翼型表面网格及翻折网格进行坐标变换,最终生成旋翼桨叶三维桨叶表面网格。
进一步,根据旋翼桨叶的扭转角对桨叶展向各个位置的剖面二维翼型表面网格及翻折网格进行旋转;根据旋翼桨叶1/4弦长位置对桨叶展向各个位置的剖面二维翼型表面网格及翻折网格进行平移。
进一步,以旋翼桨叶1/4弦长位置为旋转中心,将桨叶展向各个位置的剖面二维翼型表面网格及翻折网格进行旋转,旋转度数为旋翼桨叶的扭转角。
本发明的有益效果:本发明对背景技术进行改进,能够精确刻画桨叶表面网格信息,提供了一种计算更加精确的直升机旋翼气动噪声计算方法。
附图说明
图1是二维翼型网格示意图;
图2是三维桨叶表面网格示意图;
图3是本发明计算得到的厚度噪声云图;
图4是本发明计算得到的总噪声云图。
具体实施方式
下面结合附图对本发明所涉及的直升机旋翼气动噪声计算方法做进一步详细说明。
(1)基于刚性桨叶假设,将旋翼桨叶分解为二维翼型和一维梁模型,首先生成二维翼型网格,如图1所示,然后沿展向进行插值,最终得到三维桨叶表面网格,如图2所示。
(2)流场输入参数代入旋翼飞行器综合性分析软件计算旋翼展向分布的截面升力;
(3)噪声计算输入参数、桨叶表面网格信息和展向截面升力代入噪声计算公式求解旋翼噪声声场。
根据桨叶结构特点,具体步骤如下:
第一步,采用翼型生成软件生成翼型表面二维坐标点,将坐标点连接生成二维翼型网格;并生成沿桨叶展向分布的一维站位点;
第二步,将生成的桨叶二维翼型表面网格按照桨叶展向一维站位点分布规律进行插值并根据弦长分布规律进行缩放以得到沿桨叶展向各剖面的网格;
第三步,对桨根和桨尖部分进行翻折处理以生成O型桨尖;
第四步,根据旋翼桨叶的扭转角(旋转)和1/4弦长位置(平移)沿桨叶展向的分布规律,分别对桨叶展向各个位置的剖面二维翼型表面网格及翻折网格进行坐标变换,最终生成旋翼桨叶三维桨叶表面网格。1/4弦长位置由水平平移量和竖直平移量组成,桨叶展向各个位置的剖面二维翼型表面网格及翻折网格在桨叶剖面所在平面内根据平移量进行平移;
第五步,流场输入参数代入旋翼飞行器综合性分析软件计算旋翼展向分布的截面升力;
第六步,噪声计算输入参数、桨叶表面网格信息和展向截面升力代入以下公式求解旋翼噪声声场。
其中,p′T代表厚度噪声,p′L代表载荷噪声,a0代表无扰动介质中的音速,ρ0为密度,f=0为桨叶表面方程,v为速度,[*]ret表示括号内的变量都在延迟时间下取值;r为源项到观察点之间的距离,L为桨叶截面升力,角标r和M分别是沿传播方向和源面运动马赫数方向的投影,M为源面运动马赫数,dS表示桨叶表面的面元,dl(y)是相邻紧致源点的距离,y是紧致源点径向位置。
特别地,针对具有复杂桨尖外形的旋翼,以弦长为权重系数,采用如下插值公式获取桨尖的旋翼外形;
θ1=|r-r1|/|r2-r1|
θ2=|r-r2|/|r2-r1|
CH(r)=C1·θ2+C2·θ1
Coe1=|CH(r)-CH(r1)|/|CH(r1)-CH(r2)|
Coe2=|CH(r)-CH(r2)|/|CH(r1)-CH(r2)|
其中,θ1、θ2为展向权重系数,r1、r2、r为翼型展向位置,C1、C2为展向剖面翼型弦长,CH为插值后的弦长,Coe1、Coe2为弦长权重系数。
采用该网格生成方法能够包含前后掠桨尖、多翼型配置、非线性负扭转及尖削等现代直升机桨叶的多种特征。
厚度噪声计算结果如图3所示。总噪声云图计算结果如图4所示。
Claims (8)
1.一种直升机旋翼气动噪声计算方法,其特征在于:所述方法包括:
第一步,采用翼型生成软件生成二维翼型表面网格;
第二步,将生成的二维翼型表面网格根据弦长分布规律进行缩放得到沿桨叶展向各剖面的网格;
第三步,对桨根和桨尖部分进行处理以生成O型桨尖;
第四步,根据沿桨叶展向各剖面的网格生成旋翼桨叶三维桨叶表面网格;
第五步,计算旋翼展向分布的截面升力;
第六步,根据旋翼桨叶三维桨叶表面网格和旋翼展向分布的截面升力求解旋翼噪声声场。
2.根据权利要求1所述的一种直升机旋翼气动噪声计算方法,其特征在于:所述第一步中,还包括:生成沿桨叶展向分布的一维站位点。
3.根据权利要求2所述的一种直升机旋翼气动噪声计算方法,其特征在于:所述第二步,将生成的桨叶二维翼型表面网格按照桨叶展向一维站位点分布规律进行插值并根据弦长分布规律进行缩放以得到沿桨叶展向各剖面的网格。
4.根据权利要求3所述的一种直升机旋翼气动噪声计算方法,其特征在于:所述第三步中,对桨根和桨尖部分进行翻折处理以生成O型翻折网格。
5.根据权利要求4所述的一种直升机旋翼气动噪声计算方法,其特征在于:所述第四步中,
根据旋翼桨叶的扭转角和1/4弦长位置沿桨叶展向的分布规律,分别对桨叶展向各个位置的剖面二维翼型表面网格及翻折网格进行坐标变换,最终生成旋翼桨叶三维桨叶表面网格。
6.根据权利要求5所述的一种直升机旋翼气动噪声计算方法,其特征在于:坐标变换具体为:根据旋翼桨叶的扭转角对桨叶展向各个位置的剖面二维翼型表面网格及翻折网格进行旋转;根据旋翼桨叶1/4弦长位置对桨叶展向各个位置的剖面二维翼型表面网格及翻折网格进行平移。
7.根据权利要求6所述的一种直升机旋翼气动噪声计算方法,其特征在于:以旋翼桨叶1/4弦长位置为旋转中心,将桨叶展向各个位置的剖面二维翼型表面网格及翻折网格进行旋转,旋转度数为旋翼桨叶的扭转角。
8.根据权利要求3所述的一种直升机旋翼气动噪声计算方法,其特征在于:以弦长为权重系数,采用如下插值公式获取桨叶某个展向位置的剖面网格;
θ1=|r-r1|/|r2-r1|
θ2=|r-r2|/|r2-r1|
CH(r)=C1·θ2+C2·θ1
Coe1=|CH(r)-CH(r1)|/|CH(r1)-CH(r2)|
Coe2=|CH(r)-CH(r2)|/|CH(r1)-CH(r2)|
其中θ1、θ2为展向权重系数,r1、r2、r为翼型展向位置,C1、C2为展向剖面翼型弦长,CH为插值后的弦长,Coe1、Coe2为弦长权重系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011020831.7A CN112199772B (zh) | 2020-09-25 | 2020-09-25 | 一种直升机旋翼气动噪声计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011020831.7A CN112199772B (zh) | 2020-09-25 | 2020-09-25 | 一种直升机旋翼气动噪声计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112199772A true CN112199772A (zh) | 2021-01-08 |
CN112199772B CN112199772B (zh) | 2022-11-04 |
Family
ID=74008266
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011020831.7A Active CN112199772B (zh) | 2020-09-25 | 2020-09-25 | 一种直升机旋翼气动噪声计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112199772B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115774903A (zh) * | 2023-02-13 | 2023-03-10 | 中国空气动力研究与发展中心计算空气动力研究所 | 直升机旋翼桨盘载荷实时生成方法、装置、设备及介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104881510A (zh) * | 2015-02-13 | 2015-09-02 | 南京航空航天大学 | 一种直升机旋翼/尾桨气动干扰数值仿真方法 |
CN108021772A (zh) * | 2017-12-26 | 2018-05-11 | 南京航空航天大学 | 一种基于射流的降低旋翼噪声的方法及系统 |
CN111046493A (zh) * | 2019-12-04 | 2020-04-21 | 中国直升机设计研究所 | 一种基于集中载荷的旋翼气动噪声计算方法 |
CN111563299A (zh) * | 2020-04-30 | 2020-08-21 | 南京航空航天大学 | 一种旋翼噪声确定方法及系统 |
-
2020
- 2020-09-25 CN CN202011020831.7A patent/CN112199772B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104881510A (zh) * | 2015-02-13 | 2015-09-02 | 南京航空航天大学 | 一种直升机旋翼/尾桨气动干扰数值仿真方法 |
CN108021772A (zh) * | 2017-12-26 | 2018-05-11 | 南京航空航天大学 | 一种基于射流的降低旋翼噪声的方法及系统 |
CN111046493A (zh) * | 2019-12-04 | 2020-04-21 | 中国直升机设计研究所 | 一种基于集中载荷的旋翼气动噪声计算方法 |
CN111563299A (zh) * | 2020-04-30 | 2020-08-21 | 南京航空航天大学 | 一种旋翼噪声确定方法及系统 |
Non-Patent Citations (7)
Title |
---|
KENNETH S. BRENTNER 等: "AN ANALYTICAL COMPARISON OF THE ACOUSTIC ANALOGY AND KIRCHHOFF FORMULATION FOR MOVING SURFACES", 《AIAA JOURNAL》 * |
QIJUN ZHAO 等: "Robust Navier-Stokes method for predicting unsteady flowfield and aerodynamic characteristics of helicopter rotor", 《CHINESE JOURNAL OF AERONAUTICS》 * |
招启军 等: "基于CFD方法的倾转旋翼/螺旋桨气动优化分析", 《空气动力学学报》 * |
曹亚雄 等: "带下反桨尖旋翼气动噪声数值分析", 《直升机技术》 * |
朱正 等: "低HSI噪声旋翼桨尖外形优化设计方法", 《航空学报》 * |
王阳 等: "基于非结构网格CFD技术的旋翼气动噪声计算方法研究", 《空气动力学学报》 * |
陈丝雨 等: "基于HHC方法的旋翼噪声抑制机理及参数影响", 《航空学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115774903A (zh) * | 2023-02-13 | 2023-03-10 | 中国空气动力研究与发展中心计算空气动力研究所 | 直升机旋翼桨盘载荷实时生成方法、装置、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN112199772B (zh) | 2022-11-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109711048B (zh) | 一种考虑气动和结构多类型设计变量的机翼设计方法 | |
Brentner | An efficient and robust method for predicting helicopter high-speed impulsive noise | |
CN109878721B (zh) | 一种微小型旋翼无人飞行器旋翼翼型的设计方法及产品 | |
Shukla et al. | An inviscid model for vortex shedding from a deforming body | |
CN114036638A (zh) | 一种直升机旋翼桨-涡干扰流场和噪声高精度仿真方法 | |
CN108416075A (zh) | 基于cfd技术的风力机气动计算模型建模方法 | |
CN113051666B (zh) | 一种旋翼飞行器噪声数字化分析方法及系统 | |
CN111444661B (zh) | 交互式棱柱网格生成中翘曲现象消除方法 | |
CN113569360B (zh) | 一种风力机叶片抗颤振翼型簇设计方法 | |
CN112199772B (zh) | 一种直升机旋翼气动噪声计算方法 | |
CN109977526B (zh) | 一种基于三维cst技术的调整机翼有限元模型的方法 | |
CN109747818B (zh) | 基于桨尖质量射流的直升机旋翼气动干扰控制方法 | |
CN111563299A (zh) | 一种旋翼噪声确定方法及系统 | |
Bird et al. | A geometrically non-linear time-domain unsteady lifting-line theory | |
Cayiroglu et al. | Wing aerodynamic optimization by using genetic algoritm and Ansys | |
CN117313526A (zh) | 一种漂浮式海上风力机的耦合协同仿真计算方法 | |
CN115525988A (zh) | 一种风电机组自主载荷仿真计算与修正系统 | |
CN117272870B (zh) | 基于自适应重叠网格的动态编队飞行数值仿真方法 | |
CN112199777A (zh) | 一种适用于模化仿生前缘流场特征的方法 | |
CN115828421A (zh) | 一种直升机噪声预警方法 | |
CN114756952A (zh) | 一种悬停状态直升机旋翼气动噪声计算方法 | |
CN113935112A (zh) | 一种直升机旋翼厚度噪声的预测方法及系统 | |
CN111597741B (zh) | 一种基于改进Hicks-henne算法的仿生蟹滑翔姿态下翼型优化设计方法 | |
Kumar et al. | Computational Analysis of Morphing Geometry Inspired From Butterfly Wings in the Application Micro Aerial Vehicles | |
CN117688863A (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 |