CN107945266A - 一种机翼三维结冰模型的生成方法 - Google Patents
一种机翼三维结冰模型的生成方法 Download PDFInfo
- Publication number
- CN107945266A CN107945266A CN201711229972.8A CN201711229972A CN107945266A CN 107945266 A CN107945266 A CN 107945266A CN 201711229972 A CN201711229972 A CN 201711229972A CN 107945266 A CN107945266 A CN 107945266A
- Authority
- CN
- China
- Prior art keywords
- ice type
- wing
- dimensional
- ice
- icing
- 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
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Software Systems (AREA)
- Automation & Control Theory (AREA)
- Computer Graphics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供一种机翼三维结冰模型的生成方法,步骤如下:利用机翼结冰风洞试验得到机翼某一状态时x个截面冰型;提取截面冰型的特征长度,生成二维冰型轮廓;在二维冰型轮廓基础上根据所需厚度b,拉伸出三维冰型单元,并对其进行编号;根据厚度b,结合需要模拟的机翼三维冰型的展向长度L,获得需要组合三维冰型单元的个数n;生成n个1至x之间的随机整数数列A(n);对应数列A(n),按照展向某个方向,将每一项对应编号的冰型单元组合到机翼前缘上,生成机翼三维冰型外形模型。本发明所提供的方法,三维冰型保留了结冰试验的更多冰型特性,而且在排列冰型单元时,采用了程序生成随机数的方法,能够生成不同的三维冰型。
Description
技术领域
本发明属于飞机结冰试验领域,具体涉及一种机翼三维结冰模型的生成方法。
背景技术
在飞行中,飞机表面可能出现结冰,飞机表面的冰会对飞机气动特性产生严重的影响,如阻力增大,俯仰力矩特性剧烈变化等,严重影响飞行安全。
飞机的适航性要求飞机具有防除冰的能力以及在结冰状态下还具有安全飞行的能力。因此,研究不同气候条件对飞机结冰的影响,以及飞机结冰后对全机气动力特性的影响就非常重要。首先在结冰风洞中进行结冰试验,即在不同环境,不同大气条件下飞机的结冰情况,得到不同的截面冰型。然后进行结冰影响试验,研究不同条件结冰后对全机气动特性的影响。
飞机结冰是一个非常复杂的过程,受温度、湿度、大气中水滴情况和飞机表面条件,以及结冰时间等因素的影响,而且结冰时,冰型的生成具有很大的随机性。这给结冰过程的数值模拟,以及结冰试验结果的分析带来了很大的难度。
目前根据结冰试验截面冰型生成机翼三维冰型的方法,通用方法为:1)对比结冰试验得到的不同冰型,分析选取结冰情况最恶劣(截面冰型最复杂,结冰尺寸最大)的冰型截面,认为利用该截面冰型进行结冰影响试验,对全机气动特性影响最大;2)将选取的冰型直接沿机翼展向拉伸,获得具有二维特征的冰型模型,见图1和图2。这种冰型生成方法忽略了结冰外形在展向上的三维形状特征,并且无法体现出结冰过程中随机性,不均匀性。根据相关研究资料,研究结冰对全机气动力特性的影响时,直接拉伸获得的二维冰型,与接近真实的三维冰型的试验结果有一定差异,而且三维冰型的试验结果更加接近试飞结果。
发明内容
本发明的目的在于解决以上问题,发明一种能够模拟出机翼结冰过程随即性和快捷生成机翼三维结冰模型的生成方法。
本发明的目的通过如下技术方案实现:一种机翼三维结冰模型的生成方法,包括如下步骤,
步骤一:利用机翼结冰风洞试验得到机翼某一状态时x个截面冰型;
步骤二:提取截面冰型的特征长度,生成二维冰型轮廓;
步骤三:在二维冰型轮廓基础上根据所需厚度b,拉伸出三维冰型单元,并对三维冰型单元进行编号;
步骤四:根据步骤三中三维冰型单元的厚度b,结合需要模拟的机翼三维冰型的展向长度L,获得需要组合三维冰型单元的个数n;
步骤五:生成n个1至x之间的随机整数数列A(n);
步骤六:对应数列A(n),按照展向某个方向,将每一项对应编号的冰型单元组合到机翼前缘上,生成机翼三维冰型外形模型。
优选地是,所述步骤一中,截面冰型数量x为3~5个。
本发明所提供的一种机翼三维结冰模型的生成方法的有益效果在于,由于在生成三维冰型时应用了结冰试验得到的所有截面冰型,因此三维冰型保留了结冰试验的更多冰型特性,而且在排列冰型单元时,采用了程序生成随机数的方法,即利用相同的冰型截面,可以生成不同的三维冰型,这与机翼结冰过程的随机性相似。使用该方法,可以快捷地生成具有三维特征的机翼冰型模型;另外,由于利用冰型截面生成冰型单元,在设计和加工三维冰型时可以使用模块化的方法,使模型加工、装配过程更加简便。
附图说明
图1为现有冰型拉伸得到的机翼冰型示意图;
图2为图1的局部放大图;
图3至图6为本发明一实施例某一状态时的截面冰型示意图;
图7至图10为本发明一实施例根据二维冰型轮廓拉伸得到的三维冰型单元示意图;
图11为本发明一实施例按照随机整数数列排列得到的机翼三维冰型外形示意图;
图12为图11的局部放大图。
具体实施方式
为使本发明实施的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行更加详细的描述。在附图中,自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本发明一部分实施例,而不是全部的实施例。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
下面结合附图对本发明的机翼三维结冰模型的生成方法做进一步详细说明。
如图3至12所示,一种机翼三维结冰模型的生成方法,通过如下步骤实施。
首先利用机翼结冰风洞试验得到某一状态x个截面冰型,一般每个结冰风洞试验状态取截面冰型数量为3至5个,本实施例选取4个(见图3至图6)所示,在应用时根据具体情况进行选取,不限于4个。再提取截面冰型的特征长度,以图6中的羊角冰说明,羊角冰即冰型截面有两个较大的突起类似羊角,羊角冰的截面冰型的特征长度为高度、宽度等,适当简化,忽略较小的曲线特征,生成二维冰型轮廓。
然后利用二维冰型轮廓,根据所需厚度b,拉伸出三维冰型单元,每个截面冰型生成的三维冰型单元编号为1至4,例如b=10mm,见图7至图10所示。若需要模拟的机翼三维冰型的展向长度L,则需要组合n个三维冰型单元,对n,有n×b=L。例如,机翼三维冰型的展向长度L=500mm,则n=500/10=50,可以利用程序生成50个1至4之间的随机整数,即数列A(50):
对应数列A(50),按照展向某个方向(本发明采用从机翼根部向机翼翼尖方向,相反方向也可),将每一项对应编号的冰型单元组合到机翼前缘上,则生成机翼三维冰型外形模型,见图11和图12,使用本发明方法生成的三维冰型模型可用于加工冰型模型,进行冰型对机翼气动力影响风洞试验。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (2)
1.一种机翼三维结冰模型的生成方法,其特征在于,包括如下步骤,
步骤一:利用机翼结冰风洞试验得到机翼某一状态时x个截面冰型;
步骤二:提取截面冰型的特征长度,生成二维冰型轮廓;
步骤三:在二维冰型轮廓基础上根据所需厚度b,拉伸出三维冰型单元,并对三维冰型单元进行编号;
步骤四:根据步骤三中三维冰型单元的厚度b,结合需要模拟的机翼三维冰型的展向长度L,获得需要组合三维冰型单元的个数n;
步骤五:生成n个1至x之间的随机整数数列A(n);
步骤六:对应数列A(n),按照展向某个方向,将每一项对应编号的冰型单元组合到机翼前缘上,生成机翼三维冰型外形模型。
2.根据权利要求1所述的机翼三维结冰模型的生成方法,其特征在于,所述步骤一中,截面冰型数量x为3~5个。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711229972.8A CN107945266B (zh) | 2017-11-29 | 2017-11-29 | 一种机翼三维结冰模型的生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711229972.8A CN107945266B (zh) | 2017-11-29 | 2017-11-29 | 一种机翼三维结冰模型的生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107945266A true CN107945266A (zh) | 2018-04-20 |
CN107945266B CN107945266B (zh) | 2021-08-17 |
Family
ID=61946782
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711229972.8A Active CN107945266B (zh) | 2017-11-29 | 2017-11-29 | 一种机翼三维结冰模型的生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107945266B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110889233A (zh) * | 2019-12-04 | 2020-03-17 | 中国空气动力研究与发展中心 | 结冰风洞试验冰形三维实体生成方法及三维实物冰 |
CN111310278A (zh) * | 2020-01-17 | 2020-06-19 | 智慧航海(青岛)科技有限公司 | 一种基于仿真的船舶自动化建模方法 |
CN111323200A (zh) * | 2020-05-11 | 2020-06-23 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种结冰风洞试验冰形面积计算方法 |
CN111351630A (zh) * | 2020-03-17 | 2020-06-30 | 南京航空航天大学 | 一种3d打印夹芯复合材料模拟冰型 |
CN111368367A (zh) * | 2020-02-25 | 2020-07-03 | 空气动力学国家重点实验室 | 一种冰风洞结冰试验冰形的参数化方法 |
CN114701661A (zh) * | 2022-05-07 | 2022-07-05 | 中山大学 | 一种结冰表面原位制备共形可控冰型方法 |
CN114701661B (zh) * | 2022-05-07 | 2024-06-11 | 中山大学 | 一种结冰表面原位制备共形可控冰型方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101872373A (zh) * | 2010-05-14 | 2010-10-27 | 南京航空航天大学 | 三维机翼结冰后翼稍网格重构方法 |
CN105508152A (zh) * | 2015-12-31 | 2016-04-20 | 北京金风科创风电设备有限公司 | 叶片结冰模型的构建方法、结冰状态的监测方法和装置 |
-
2017
- 2017-11-29 CN CN201711229972.8A patent/CN107945266B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101872373A (zh) * | 2010-05-14 | 2010-10-27 | 南京航空航天大学 | 三维机翼结冰后翼稍网格重构方法 |
CN105508152A (zh) * | 2015-12-31 | 2016-04-20 | 北京金风科创风电设备有限公司 | 叶片结冰模型的构建方法、结冰状态的监测方法和装置 |
Non-Patent Citations (3)
Title |
---|
孔满昭: "翼型结冰过程数值模拟验算与分析", 《航空工程进展》 * |
常士楠: "三维机翼结冰模拟", 《航空学报》 * |
李寒清: "三维机翼结冰外形生成及流场分析", 《复旦学报(自然科学版)》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110889233A (zh) * | 2019-12-04 | 2020-03-17 | 中国空气动力研究与发展中心 | 结冰风洞试验冰形三维实体生成方法及三维实物冰 |
CN110889233B (zh) * | 2019-12-04 | 2023-03-28 | 中国空气动力研究与发展中心 | 结冰风洞试验冰形三维实体生成方法及三维实物冰 |
CN111310278A (zh) * | 2020-01-17 | 2020-06-19 | 智慧航海(青岛)科技有限公司 | 一种基于仿真的船舶自动化建模方法 |
CN111310278B (zh) * | 2020-01-17 | 2023-05-02 | 智慧航海(青岛)科技有限公司 | 一种基于仿真的船舶自动化建模方法 |
CN111368367A (zh) * | 2020-02-25 | 2020-07-03 | 空气动力学国家重点实验室 | 一种冰风洞结冰试验冰形的参数化方法 |
CN111368367B (zh) * | 2020-02-25 | 2023-08-25 | 空气动力学国家重点实验室 | 一种冰风洞结冰试验冰形的参数化方法 |
CN111351630A (zh) * | 2020-03-17 | 2020-06-30 | 南京航空航天大学 | 一种3d打印夹芯复合材料模拟冰型 |
CN111323200A (zh) * | 2020-05-11 | 2020-06-23 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种结冰风洞试验冰形面积计算方法 |
CN114701661A (zh) * | 2022-05-07 | 2022-07-05 | 中山大学 | 一种结冰表面原位制备共形可控冰型方法 |
CN114701661B (zh) * | 2022-05-07 | 2024-06-11 | 中山大学 | 一种结冰表面原位制备共形可控冰型方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107945266B (zh) | 2021-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107945266A (zh) | 一种机翼三维结冰模型的生成方法 | |
CN105203291B (zh) | 一种用于矢量推进飞机模型的风洞试验系统 | |
CN105115692B (zh) | 一种cfd数值模拟与风洞试验相结合的气动预测方法 | |
Reid et al. | FENSAP-ICE simulation of icing on wind turbine blades, part 1: performance degradation | |
CN108304597B (zh) | 一种高速飞行器头部前缘转捩预测装置及方法 | |
CN107103123A (zh) | 基于Bagging‑FSE算法的风洞马赫数混合模型建模方法 | |
Christodoulou et al. | Aerodynamic Analysis of a Quadcopter Drone Propeller with the Use of Computational Fluid Dynamics. | |
Fenercioglu et al. | Effect of unequal flapping frequencies on flow structures | |
CN115356071A (zh) | 一种对转桨扇缩尺试验件风洞试验数据修正方法和系统 | |
CN108341073A (zh) | 飞机颤振高温试验图像测试方法 | |
Ara et al. | Performance study of winglets on tapered wing with curved trailing edge | |
Potapczuk et al. | LEWICE modeling of swept wing ice accretions | |
Ni et al. | A numerical study of the anti-icing heat load for a three-dimensional aircraft wing | |
Venkataramani et al. | Wind tunnel tests and modeling studies of ice accretion relevant to aircraft engines | |
Laban et al. | Reshaping engine nacelles for testing in wind tunnels with turbofan propulsion simulators | |
Maneia et al. | Aerodynamics of a rigid curved kite wing | |
Pena et al. | Parallel computing of icing on threedimensional airfoils | |
Gong et al. | Numerical simulation of powered high-lift flow | |
Zohary et al. | A successful microblowing strategy for airfoils | |
Sun et al. | Ice accretion prediction for single airfoil and multi-element airfoil | |
Viguier | Wind tunnel test of a non-linear flutter model | |
Shan et al. | Numerical Simulation Research on Aerodynamic Characteristics of the Fixed Transition Airfoil | |
Parlindungan et al. | An Experimental Analysis on the Effects of Stagger on the Aerodynamic Forces of Tandem Wings | |
Griffin et al. | Numerical and experimental study on the ability of dynamic roughness to alter the development of a leading edge vortex | |
Gates et al. | Servo-flap rotor power improvement through wind tunnel testing |
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 |