CN106599353A - 一种飞机多段翼型外流场动态数值模拟方法 - Google Patents
一种飞机多段翼型外流场动态数值模拟方法 Download PDFInfo
- Publication number
- CN106599353A CN106599353A CN201610972223.3A CN201610972223A CN106599353A CN 106599353 A CN106599353 A CN 106599353A CN 201610972223 A CN201610972223 A CN 201610972223A CN 106599353 A CN106599353 A CN 106599353A
- Authority
- CN
- China
- Prior art keywords
- model
- aircraft
- flow field
- parameter
- airfoils
- 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
Classifications
-
- 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
-
- 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)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明涉及一种飞机多段翼型外流场动态数值模拟方法,该方法包括如下步骤:(1)建立带有前缘缝翼后缘襟翼多段翼型飞机模型并导入fluent软件;(2)构建飞机模型的垂直下落运动函数,并导入fluent软件;(3)选取湍流模型;(4)设置仿真参数,包括流体参数、翼型材料参数、工作环境参数、边界条件参数、动网格更新模型、求解方法和求解精度;(5)运行fluent软件,获取飞机垂直下落过程中翼型周围流体流动参数变化图,包括速度、压力和温度变化图。与现有技术相比,本发明能有效获取飞机降落过程中的外流场动态数值,便于对飞机起降性能评估以及增升装置的设计提供参考。
Description
技术领域
本发明涉及一种飞机外流场数值模拟方法,尤其是涉及一种飞机多段翼型外流场动态数值模拟方法。
背景技术
飞机起飞着陆过程中,因为距离地面较近,有很大可能受到地面效应的影响。所谓地面效应,就是一种使飞行器诱导阻力减小,同时能获得比空中飞行更高升阻比的流体力学效应:当运动的飞行器掉到距地面(或水面)很近时,整个飞行器体的上下压力差增大,升力会陡然增加。飞机通常利用增升装置来增加低速时飞机的升力、增加载重、增加航程、减轻飞机重量等,增升装置多由前缘缝翼、主翼及后缘襟翼组成,飞机的安全性与经济性受到增升装置的很大影响。
由于模型尺寸的限制,风洞试验很难真实观察到机翼中如缝道、空腔和拐角等处的流动细节,使人们很难对整个多段翼型周围的流动状态有完整的理解。对多段翼型进行仿真数值模拟,可高效率、低成本、快速获取最佳的翼型周围流动状态,为后期的参数优化提供完整的技术支持。
传统的机翼近地面仿真数值模拟领域,大多数采用定常状态下进行模拟计算,通过改变迎角等方式进行模拟,这种方式对于高度固定的状态下研究飞行迎角、展弦比和后掠角对地面效应的流动机理以及气动性能有较大帮助。
然而,对于飞机下落过程中翼型位移运动问题,不能使用定常的方法进行模拟。模拟运动过程中的状态更靠近降落中的真实状态,才能直接体现出飞机降落过程中地面效应对于多段机翼周围流动的影响,因此研究一种飞机降落过程中飞机多段翼型外流场动态数值模拟方法十分必要。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种飞机多段翼型外流场动态数值模拟方法。
本发明的目的可以通过以下技术方案来实现:
一种飞机多段翼型外流场动态数值模拟方法,该方法包括如下步骤:
(1)建立带有前缘缝翼后缘襟翼多段翼型飞机模型并导入fluent软件;
(2)构建飞机模型的垂直下落运动函数,并导入fluent软件;
(3)选取湍流模型;
(4)设置仿真参数,包括流体参数、翼型材料参数、工作环境参数、边界条件参数、动网格更新模型、求解方法和求解精度;
(5)运行fluent软件,获取飞机垂直下落过程中翼型周围流体流动参数变化图,包括速度、压力和温度变化图。
步骤(1)中采用ANSYS ICEM软件进行带有前缘缝翼后缘襟翼多段翼型飞机模型的建立。
所述的湍流模型采用SST k-ω模型。
所述的流体参数包括密度和粘度。
所述的工作环境参数包括大气压大小以及域环境,其中大气压设置为零,域环境设置为流体域。
所述的边界条件参数包括入口边界的边界类型和壁面边界的边界类型。
设置的动网格更新模型具体为:前缘襟翼、主翼及后缘襟翼部位的动网格模型采用局部重划模型,其余区域采用弹簧近似光滑模型。
求解方法仿真参数设定包括压力速度耦合方程模式设定、压力方程的松弛因子设定、压力修正方程模式设定。
求解精度包括迭代参数步长设置为0.0005s,步数设置为10000。
与现有技术相比,本发明具有如下优点:
(1)采用fluent软件进行仿真,建立带有前缘缝翼后缘襟翼多段翼型飞机模型能够有效模拟飞机实际构造,同时通过垂直下落运动函数来构建飞机下落过程,进而通过fluent软件来仿真飞机下落过程,对翼型周围流体流动情况,包括速度、压力、温度分布数据进行测量,很好地反应飞机降落过程中飞机多段翼型外流场动态数值,便于对飞机起降性能评估以及增升装置的设计提供参考;
(2)本发明提供了一种优于现有飞机机翼气动性能数值模拟,传统机翼气动性能模拟方法一般采用定常的固定飞行高度的方法进行流体力学仿真,对比不同高度下周围流体的流动情况,这种方法对于翼型动态下落过程会产生较大误差;
(3)围绕着前缘襟翼、主翼及后缘襟翼部位需要做下降的刚体运动,由于运动边界的位移远远大于网格尺寸,采用弹簧近似光滑模型会导致网格质量下降,造成网格畸变,导致计算不收敛,因此采用局部重划模型,使重新划分的网格满足畸变率和尺寸的要求,其余部位采用弹簧近似光滑模型,从而使得最后仿真结果更加精可靠。
附图说明
图1为本发明的流程框图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例
如图1所示,一种飞机多段翼型外流场动态数值模拟方法,该方法包括如下步骤:
(1)建立带有前缘缝翼后缘襟翼多段翼型飞机模型并导入fluent软件;
(2)构建飞机模型的垂直下落运动函数,并导入fluent软件;
(3)选取湍流模型;
(4)设置仿真参数,包括流体参数、翼型材料参数、工作环境参数、边界条件参数、动网格更新模型、求解方法和求解精度;
(5)边界初始化,设置监视窗口,设置迭代参数,设置自动截图,运行fluent软件,进行自动计算,从而获取飞机垂直下落过程中翼型周围流体流动参数变化图,包括速度、压力和温度变化图。
步骤(1)中采用ANSYS ICEM软件进行带有前缘缝翼后缘襟翼多段翼型飞机模型的建立。具体地:采用ANSYS ICEM软件进行带有前缘缝翼后缘襟翼多段翼型飞机模型,本实施例中翼型为麦道航空30p-30n,根据最终多段翼型模型的结构特点,选取网格类型划分网格,网格单位数为16万,设置边界条件,输出网格文件;将输出的网格文件导入fluent软件中;检测网络check,用于检测步骤一中模型建立后网格划分是否正确;根据多段翼型的真实尺寸设定单位。
所述的湍流模型采用SST k-ω模型。
所述的流体参数包括密度和粘度。
所述的工作环境参数包括大气压大小以及域环境,其中大气压设置为零,域环境设置为流体域。
所述的边界条件参数包括入口边界的边界类型和壁面边界的边界类型,本实施例中设置入口边界的边界类型为“Pressure Far-Field”,设置壁面边界的边界类型为“Wall”。边界条件为设置表压为11111Pa,设置马赫数为0.2,设置速度向量为直角坐标方式。设置方向向量为(0.96126,0.27563)。该向量为通过攻角16°计算获得,Cos16°=0.96126,sin16°=0.27563。
设置的动网格更新模型具体为:前缘襟翼、主翼及后缘襟翼部位的动网格模型采用局部重划模型,其余区域采用弹簧近似光滑模型。弹簧近似光滑模型中的参数设定:弹簧倔强系数设定为0.05,边界节点松弛设定为0.5;收敛判据保持默认的0.001;迭代次数保持默认的20。局部重划模型中的参数设定:将最小长度比例及最大长度比例均设置为0,最大单元畸变设定为0.4,依照尺寸标准重新划分的间隔设定为1,尺寸函数变量设置为0.3,尺寸函数变化率保持默认的0.3。
求解方法仿真参数设定包括压力速度耦合方程模式设定、压力方程的松弛因子设定、压力修正方程模式设定。压力速度耦合方程模式选择“SIMPLE模式”,压力方程的松弛因子设为1,压力修正方程模式选择“Presto模式”,其他方程的差分格式选择二阶迎风模式(Second Order Upwind)。求解精度包括迭代参数步长设置为0.0005s,步数设置为10000。
Claims (9)
1.一种飞机多段翼型外流场动态数值模拟方法,其特征在于,该方法包括如下步骤:
(1)建立带有前缘缝翼后缘襟翼多段翼型飞机模型并导入fluent软件;
(2)构建飞机模型的垂直下落运动函数,并导入fluent软件;
(3)选取湍流模型;
(4)设置仿真参数,包括流体参数、翼型材料参数、工作环境参数、边界条件参数、动网格更新模型、求解方法和求解精度;
(5)运行fluent软件,获取飞机垂直下落过程中翼型周围流体流动参数变化图,包括速度、压力和温度变化图。
2.根据权利要求1所述的一种飞机多段翼型外流场动态数值模拟方法,其特征在于,步骤(1)中采用ANSYS ICEM软件进行带有前缘缝翼后缘襟翼多段翼型飞机模型的建立。
3.根据权利要求1所述的一种飞机多段翼型外流场动态数值模拟方法,其特征在于,所述的湍流模型采用SST k-ω模型。
4.根据权利要求1所述的一种飞机多段翼型外流场动态数值模拟方法,其特征在于,所述的流体参数包括密度和粘度。
5.根据权利要求1所述的一种飞机多段翼型外流场动态数值模拟方法,其特征在于,所述的工作环境参数包括大气压大小以及域环境,其中大气压设置为零,域环境设置为流体域。
6.根据权利要求1所述的一种飞机多段翼型外流场动态数值模拟方法,其特征在于,所述的边界条件参数包括入口边界的边界类型和壁面边界的边界类型。
7.根据权利要求1所述的一种飞机多段翼型外流场动态数值模拟方法,其特征在于,设置的动网格更新模型具体为:前缘襟翼、主翼及后缘襟翼部位的动网格模型采用局部重划模型,其余区域采用弹簧近似光滑模型。
8.根据权利要求1所述的一种飞机多段翼型外流场动态数值模拟方法,其特征在于,求解方法仿真参数设定包括压力速度耦合方程模式设定、压力方程的松弛因子设定、压力修正方程模式设定。
9.根据权利要求1所述的一种飞机多段翼型外流场动态数值模拟方法,其特征在于,求解精度包括迭代参数步长设置为0.0005s,步数设置为10000。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610972223.3A CN106599353A (zh) | 2016-11-04 | 2016-11-04 | 一种飞机多段翼型外流场动态数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610972223.3A CN106599353A (zh) | 2016-11-04 | 2016-11-04 | 一种飞机多段翼型外流场动态数值模拟方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106599353A true CN106599353A (zh) | 2017-04-26 |
Family
ID=58590918
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610972223.3A Pending CN106599353A (zh) | 2016-11-04 | 2016-11-04 | 一种飞机多段翼型外流场动态数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106599353A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107729638A (zh) * | 2017-10-09 | 2018-02-23 | 中国民航大学 | 各向异性大气湍流场数值模拟方法 |
CN107748052A (zh) * | 2017-11-25 | 2018-03-02 | 南京航空航天大学 | 一种基于piv方法的测量襟翼缝道流动的装置 |
CN107818197A (zh) * | 2017-09-15 | 2018-03-20 | 南京航空航天大学 | 一种基于piv技术的超声速翼型的测力方法和装置 |
CN109614644A (zh) * | 2018-11-02 | 2019-04-12 | 中国航空工业集团公司西安飞机设计研究所 | 一种外吹式襟翼布局飞机动力增升效果评估方法 |
CN111339672A (zh) * | 2020-03-02 | 2020-06-26 | 上海索辰信息科技有限公司 | 进气道前缘激波气动热仿真分析方法 |
CN113569440A (zh) * | 2021-06-10 | 2021-10-29 | 上海工程技术大学 | 一种模拟机场区域飞机起降污染物排放扩散的方法 |
CN113704886A (zh) * | 2021-08-16 | 2021-11-26 | 成都飞机工业(集团)有限责任公司 | 一种快速、择优的缝道翼型设计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104765927A (zh) * | 2015-04-16 | 2015-07-08 | 西北工业大学 | 基于多学科的飞机增升装置高低速综合优化设计方法 |
CN104978449A (zh) * | 2015-08-17 | 2015-10-14 | 北京航空航天大学 | 二维三段翼型的前缘缝翼及后缘襟翼位置的气动优化方法 |
CN105775159A (zh) * | 2016-03-07 | 2016-07-20 | 南京航空航天大学 | 具有抑制机翼分离流功能的吹气口设计方法 |
-
2016
- 2016-11-04 CN CN201610972223.3A patent/CN106599353A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104765927A (zh) * | 2015-04-16 | 2015-07-08 | 西北工业大学 | 基于多学科的飞机增升装置高低速综合优化设计方法 |
CN104978449A (zh) * | 2015-08-17 | 2015-10-14 | 北京航空航天大学 | 二维三段翼型的前缘缝翼及后缘襟翼位置的气动优化方法 |
CN105775159A (zh) * | 2016-03-07 | 2016-07-20 | 南京航空航天大学 | 具有抑制机翼分离流功能的吹气口设计方法 |
Non-Patent Citations (2)
Title |
---|
佟增军 等: "二维多段翼型襟翼吹吸气增升控制的数值模拟", 《中国力学学会学术大会 2009论文集》 * |
朱一西 等: "多段翼型非定常地面效应数值模拟", 《空气动力学学报》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107818197A (zh) * | 2017-09-15 | 2018-03-20 | 南京航空航天大学 | 一种基于piv技术的超声速翼型的测力方法和装置 |
CN107818197B (zh) * | 2017-09-15 | 2021-05-25 | 南京航空航天大学 | 一种基于piv技术的超声速翼型的测力方法和装置 |
CN107729638A (zh) * | 2017-10-09 | 2018-02-23 | 中国民航大学 | 各向异性大气湍流场数值模拟方法 |
CN107748052A (zh) * | 2017-11-25 | 2018-03-02 | 南京航空航天大学 | 一种基于piv方法的测量襟翼缝道流动的装置 |
CN107748052B (zh) * | 2017-11-25 | 2018-09-21 | 南京航空航天大学 | 一种基于piv测速方式的测量襟翼缝道流动的装置 |
CN109614644A (zh) * | 2018-11-02 | 2019-04-12 | 中国航空工业集团公司西安飞机设计研究所 | 一种外吹式襟翼布局飞机动力增升效果评估方法 |
CN109614644B (zh) * | 2018-11-02 | 2023-03-14 | 中国航空工业集团公司西安飞机设计研究所 | 一种外吹式襟翼布局飞机动力增升效果评估方法 |
CN111339672A (zh) * | 2020-03-02 | 2020-06-26 | 上海索辰信息科技有限公司 | 进气道前缘激波气动热仿真分析方法 |
CN113569440A (zh) * | 2021-06-10 | 2021-10-29 | 上海工程技术大学 | 一种模拟机场区域飞机起降污染物排放扩散的方法 |
CN113704886A (zh) * | 2021-08-16 | 2021-11-26 | 成都飞机工业(集团)有限责任公司 | 一种快速、择优的缝道翼型设计方法 |
CN113704886B (zh) * | 2021-08-16 | 2023-10-03 | 成都飞机工业(集团)有限责任公司 | 一种快速、择优的缝道翼型设计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106599353A (zh) | 一种飞机多段翼型外流场动态数值模拟方法 | |
Patel et al. | CFD Analysis of an Aerofoil | |
CN106840572B (zh) | 一种临近空间大展弦比柔性飞行器风洞试验数据修正方法 | |
CN102066196B (zh) | 主动地使空气动力学轮廓变形的方法 | |
CN105183996A (zh) | 面元修正与网格预先自适应计算方法 | |
EP3509944A1 (en) | Blade or wing | |
CN104317980A (zh) | 一种共轴刚性旋翼气动弹性响应方法 | |
Mohd et al. | COMPUTATIONALAERODYNAMICS OF HOVERING HELICOPTER ROTORS | |
Sitorus et al. | The role of elytra in beetle flight: I. Generation of quasi-static aerodynamic forces | |
Fisher | The effect of freestream turbulence on fixed and flapping micro air vehicle wings | |
Qin et al. | Numerical study of aerodynamic forces and flow physics of a delta wing in dynamic ground effect | |
Ryu et al. | Aerodynamic Analysis of the Ducted Fan for a VTOL UAV in Crosswinds | |
Folkersma et al. | Steady-state aeroelasticity of a ram-air wing for airborne wind energy applications | |
CN105129071B (zh) | 太阳能飞机翼型设计方法及太阳能飞机翼型 | |
Datta et al. | Validation of Structural and Aerodynamic Modeling Using UH‐60A Airloads Program Data | |
Reveles et al. | Rotor/propeller model fidelity effects on propeller-wing interactional aerodynamics | |
Jain et al. | Investigation on 3-d wing of commercial aeroplane with aerofoil naca 2415 using cfd fluent | |
Willems et al. | Verification, validation, and application of shear stress transport transitional model to a R/C aircraft | |
Andan et al. | Aerodynamics of a hybrid airship | |
Boffadossi et al. | Analysis on aerodynamic characteristics of a paraglider airfoil | |
Sitorus et al. | Hydrodynamic characteristics of cambered NACA0012 for flexible-wing application of a flapping-type tidal stream energy harvesting system | |
CN109871562B (zh) | 一种基于Catia二次开发计算飞机水上迫降静水面漂浮特性的方法 | |
Khalil et al. | Computational analyses of aerodynamic characteristics of naca653218airfoil | |
Roy et al. | Delay or control of flow separation for enhanced aerodynamic performance using an effective morphed surface | |
Wild et al. | Comparison of Instantaneous Aerodynamic Loads on Sharp and Blunt Trailing-Edged Blades of High Advance Ratio Rotors |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170426 |