CN115649417B - 一种高亚声速自配平高隐身翼型 - Google Patents
一种高亚声速自配平高隐身翼型 Download PDFInfo
- Publication number
- CN115649417B CN115649417B CN202211601613.1A CN202211601613A CN115649417B CN 115649417 B CN115649417 B CN 115649417B CN 202211601613 A CN202211601613 A CN 202211601613A CN 115649417 B CN115649417 B CN 115649417B
- Authority
- CN
- China
- Prior art keywords
- airfoil
- stealth
- wing
- profile
- self
- 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.)
- Active
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
- Turbine Rotor Nozzle Sealing (AREA)
Abstract
本发明属于翼型设计领域,公开了一种高亚声速自配平高隐身翼型,所述翼型的最大厚度为8.38%,所述翼型的最大厚度位置为32.6%,所述翼型的最大弯度为1.57%,所述翼型最大弯度位置为30.5%,所述翼型后缘厚度为1.1%。本发明解决了现有的飞翼布局飞行器无法兼顾自配平和高隐身性能以满足高亚音速飞行的问题,适用于飞翼布局飞行器高亚音速巡航。
Description
技术领域
本发明涉及翼型设计领域,尤其涉及一种高亚声速自配平高隐身翼型。
背景技术
翼型是飞行器气动布局设计的基础,直接影响飞行器的飞行性能,特别是高亚音速飞翼布局飞行器对翼型的设计要求更加严格。为了提高飞翼飞行器的高亚音速巡航性能和隐身性能,飞翼布局飞行器通常取消平垂尾,整个飞行器气动布局非常干净,具有极高的气动性能和隐身性能。翼型弯度、前缘半径、尾缘角和翼型厚度等几何参数都会对其隐身特性有重要影响,通过优化设计翼型结构可以改善翼型的隐身特性,翼型隐身重点需要抑制其前向和后向雷达散射截面。需要注意的是:某些翼型参数对气动和隐身性能均有较大影响,但作用相反,应作为综合优化设计的主要设计参数,并采用不同的优化设计策略寻求总体最优结果。
对于飞翼布局飞行器,由于机身长度通常较短,舵面尺寸有限,导致飞翼布局的机动性能比较差,因此为了提高飞行器的操纵能力,翼型设计中通常要求具有自配平能力,减小俯仰舵面配平偏转导致的气动性能和隐身性能下降,同时增加巡航状态下的舵面有效偏转范围。针对高亚音速飞翼布局对翼型的特殊需求,从气动性能、隐身性能和自配平等方面设计出发,急需设计了一种能够实现自配平的高亚音速高隐身翼型。
发明内容
本发明意在提供一种高亚声速自配平高隐身翼型,以解决现有的飞翼布局飞行器无法兼顾自配平和高隐身性能以满足高亚音速飞行的问题。
为了实现上述目的,本发明提供如下技术方案:
一种高亚声速自配平高隐身翼型,所述翼型的最大厚度为8.38%,所述翼型的最大厚度位置为32.6%,所述翼型的最大弯度为1.57%,所述翼型最大弯度位置为30.5%,所述翼型后缘厚度为1.1%。
优选地,所述翼型在以0.8马赫数下巡航时,最大升阻比为14.5,对应的攻角为2°。
优选地,所述翼型在以0.8马赫数下巡航时,失速迎角为16°以上,翼型最大升阻比超过14。
优选地,所述翼型在以0.7马赫数下巡航时,最大升阻比超过54。
优选地,所述翼型的设计方法包括如下步骤:
S1、首先将翼型进行参数化,选择FFD自由变形技术进行翼型表面的参数化;
S2、针对FFD自由变形技术的特点,在均匀控制格点的基础上实现非均匀控制点的网格变形控制方法;
S3、结合CFD数值计算和电磁隐身计算工具得到样本点的气动参数和隐身性能;
S4、基于样本点构建代理模型,根据目标函数选择多目标优化算法;
S5、通过多轮优化得到所需的高亚音速翼型。
优选地,所述网格变形控制方法为:根据翼型型面分布设计翼型非均匀控制点,建立控制点与网格表面节点的对应关系:
其中是翼型表面网格在局部坐标系中的局部坐标,采用最小二乘法得到翼型表面坐标的局部坐标/>,/>表示FFD框的节点个数,/>是FFD控制点的坐标,/>是翼型表面的原始坐标,/>、/>是贝塞尔基函数,/>分别是x方向和y方向的节点数,移动翼型附近的控制点/>到新的位置/>,由于局部坐标不变,从而可以得到变形后的网格坐标为:
本发明的有益效果:
本方案能够实现自配平的高亚声速高隐身翼型,M=0.8对应的设计升力系数为0.26,M=0.7对应的设计升力系数为0.5,该翼型具有更小的阻力系数和雷达散射截面。优化得到的翼型M=0.8时失速迎角超过16°,具有较缓和的失速特性,最大升阻比为14.5,M=0.7对应的最大升阻比为54,可以满足飞翼布局飞行器高亚音速下的高性能设计需求。
附图说明
图1为本发明一种高亚声速自配平高隐身翼型的几何图;
图2为本发明一种高亚声速自配平高隐身翼型的升力系数曲线;
图3为本发明一种高亚声速自配平高隐身翼型的俯仰力矩特性曲线;
图4为本发明一种高亚声速自配平高隐身翼型的升阻比特性曲线;
图5为本发明一种高亚声速自配平高隐身翼型的设计状态的压力系数图;
图6为本发明一种高亚声速自配平高隐身翼型的设计状态的压力云图;
图7为本发明一种高亚声速自配平高隐身翼型在M=0.7时的升阻比曲线;
图8为本发明一种高亚声速自配平高隐身翼型在M=0.7时的俯仰力矩特性曲线。
注:图中X、Y分别指的是翼型沿着x、y方向的坐标,单位为m。AOA指的是翼型的攻角,CL指的是翼型升力系数,CM指的是翼型俯仰力矩系数,K指的是翼型升阻比,cp指的是翼型表面压力系数。
具体实施方式
下面结合附图和实施方式对本发明作进一步的详细说明:
如图1所示,一种高亚声速自配平高隐身翼型,翼型的最大厚度为8.38%,翼型的最大厚度位置为32.6%,翼型的最大弯度为1.57%,翼型最大弯度位置为30.5%,翼型后缘厚度为1.1%。翼型上表面相对平缓,在翼型后缘采用后卸载的设计,本设计相当于在翼型后缘产生负升力增加抬头力矩,使得翼型在巡航状态下实现力矩自配平,通过调整翼型表面型线的分布改善翼型的隐身特性。
翼型能够实现在M=0.8时,雷诺数百万量级具有较好的气动性能,翼型下表面前半部分相对平缓,翼型前缘曲率半径较小,一方面是减小波阻,一方面是减小雷达散射面积。翼型前半部分设计有一定的弯度,使得翼型在小迎角下具有较大的升力系数,整个翼型的升力贡献主要是由前半部分产生,由于零升迎角的翼型产生的是抬头力矩。高亚音速飞机机翼设计一般有一定的后掠角,同时保证一定的静稳定裕度,因此翼型自身的抬头力矩可以通过静稳定裕度的设计使得飞机达到力矩自配平。
飞翼布局在高亚音速下巡航飞行,如果迎角过大,翼型容易上表面出现较大范围内的超音速区域和表面激波,同时激波也会诱导翼型上表面分离,导致飞翼布局的气动性能明显下降,另外翼型的曲率分布对隐身性能影响较大,特别是翼型前缘半径、翼型弯度和翼型厚度等参数。
本翼型的设计方法包括如下步骤:
S1、选择某基准翼型进行设计,首先将翼型进行参数化,选择FFD自由变形技术进行翼型表面的参数化;对基准翼型进行网格划分,然后根据参数化方法确定FFD框,并选择一定(即从设计的FFD控制节点中选择一部分作为变量)的控制点作为控制变量。
S2、针对FFD自由变形技术的特点,在均匀控制格点的基础上实现非均匀控制点的网格变形控制方法。首先根据翼型型面分布设计翼型非均匀控制点,建立控制点与网格表面节点的对应关系:
其中是翼型表面网格在局部坐标系中的局部坐标,采用最小二乘法得到翼型表面坐标的局部坐标/>,/>表示FFD框的节点个数,/>是FFD控制点的坐标,/>是翼型表面的原始坐标,一个/>对应唯一的 />组合,/>、/>是贝塞尔基函数,/>分别是x方向和y方向的节点数,移动翼型附近的控制点/>到新的位置/>,由于局部坐标不变,从而可以得到变形后的网格坐标为:
S3、结合CFD数值计算和电磁隐身计算工具得到样本点(即所有控制点变形状态的组合)的气动参数和隐身性能。
S4、基于样本点构建代理模型,根据目标函数选择多目标优化算法,以巡航状态点升力系数为气动约束,阻力最小和雷达散射截面最小为优化目标,建立翼型优化设计流程。
S5、通过多轮优化得到所需的高亚音速翼型。
图1中X和Y为无量纲的数值,x,y是翼型沿横轴纵轴的坐标。
图2给出的是本翼型在M=0.8时的升力系数,最大升力系数达到0.8以上,当迎角较小时升力线斜率较大,当迎角继续增大时升力线斜率变小,但是翼型仍没有失速,最大失速迎角在16°以上。
图3可以看出在小迎角范围内翼型的力矩基本保持不变,翼型产生的是抬头力矩。
图4可以看出最大升阻比出现在2°左右,最大升阻比为14.5左右。
图5是M=0.8设计点下的翼型上下表面压力系数对比,可以看出在设计状态下,翼型上表面50%区域出现了明显的压力增加过程,说明翼型表面出现了激波,由于升力产生的主要区域是翼型前半部分,翼型尾部区域出现了负升力,负升力产生的是抬头力矩,前半部分压力作用和尾部共同产生了翼型的抬头力矩。
图6是翼型设计点下的压力云图,翼型下表面的压力分布相对平缓,翼型上表面出现了明显的激波,这也是高亚音速翼型阻力增加的主要原因。
图7和图8给出了翼型在马赫数为0.7时的气动性能,可以看出本翼型在M=0.7也具有较好的气动性能,且在一定范围内翼型能够实现力矩自配平,升阻比也出现了明显的增加,原因主要是在M=0.7时翼型表面出现的超音速区域很小且激波较弱,证明了本专利设计的翼型具有较好的宽速度适应范围。
以上所述的仅是本发明的实施例,实施方案中公知的具体技术方案或特性等常识在此未作过多描述。应当指出,对于本领域的技术人员来说,在不脱离本发明技术方案的前提下,还可以作出若干变形和改进,这些也应该视为本发明的保护范围,这些都不会影响本发明实施的效果和专利的实用性。本申请要求的保护范围应当以其权利要求的内容为准,说明书中的具体实施方式等记载可以用于解释权利要求的内容。
Claims (4)
1.一种高亚声速自配平高隐身翼型,其特征在于:所述翼型的最大厚度为8.38%,所述翼型的最大厚度位置为32.6%,所述翼型的最大弯度为1.57%,所述翼型最大弯度位置为30.5%,所述翼型后缘厚度为1.1%;
所述翼型的设计方法包括如下步骤:
S1、首先将翼型进行参数化,选择FFD自由变形技术进行翼型表面的参数化;
S2、针对FFD自由变形技术的特点,在均匀控制格点的基础上实现非均匀控制点的网格变形控制方法;
S3、结合CFD数值计算和电磁隐身计算工具得到样本点的气动参数和隐身性能;
S4、基于样本点构建代理模型,根据目标函数选择多目标优化算法;
S5、通过多轮优化得到所需的高亚音速翼型;
所述网格变形控制方法为:根据翼型型面分布设计翼型非均控制点,建立控制点与网格表面节点的对应关系:
其中s,t是翼型表面网格在局部坐标系中的局部坐标,采用最小二乘法得到翼型表面坐标的局部坐标s,t,l,m表示FFD框的节点个数,是FFD控制点的坐标,X1是翼型表面的原始坐标,/>是贝塞尔基函数,i,j分别是x方向和y方向的节点数,移动翼型附近的控制点/>到新的位置/>由于局部坐标不变,从而可以得到变形后的网格坐标为:
2.根据权利要求1所述的一种高亚声速自配平高隐身翼型,其特征在于:所述翼型在以0.8马赫数下巡航时,最大升阻比为14.5,对应的攻角为2°。
3.根据权利要求1所述的一种高亚声速自配平高隐身翼型,其特征在于:所述翼型在以0.8马赫数下巡航时,失速迎角为16°以上,翼型最大升阻比超过14。
4.根据权利要求1所述的一种高亚声速自配平高隐身翼型,其特征在于:所述翼型在以0.7马赫数下巡航时,最大升阻比超过54。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211601613.1A CN115649417B (zh) | 2022-12-14 | 2022-12-14 | 一种高亚声速自配平高隐身翼型 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211601613.1A CN115649417B (zh) | 2022-12-14 | 2022-12-14 | 一种高亚声速自配平高隐身翼型 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115649417A CN115649417A (zh) | 2023-01-31 |
CN115649417B true CN115649417B (zh) | 2023-06-30 |
Family
ID=85022814
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211601613.1A Active CN115649417B (zh) | 2022-12-14 | 2022-12-14 | 一种高亚声速自配平高隐身翼型 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115649417B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105752314A (zh) * | 2016-03-22 | 2016-07-13 | 西北工业大学 | 一种高空低速自然层流高升力翼型 |
CN114684348A (zh) * | 2022-04-01 | 2022-07-01 | 西北工业大学 | 一种适用于轻中型高速公务机的高阻力发散马赫数层流翼型 |
Family Cites Families (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7070146B2 (en) * | 2003-08-29 | 2006-07-04 | Supersonic Aerospace International, Llc | Aircraft thickness/camber control device for low sonic boom |
US8473254B2 (en) * | 2003-11-14 | 2013-06-25 | Supersonic Aerospace International, Llc | Methods for configuring aircraft to meet performance goals and shock wave disturbance constraints |
US20050197811A1 (en) * | 2004-02-02 | 2005-09-08 | Ariyur Kartik B. | Prediction of dynamic ground effect forces for fixed wing aircraft |
US8186616B2 (en) * | 2004-12-21 | 2012-05-29 | Israel Aerospace Industries Ltd. | Hybrid transonic-subsonic aerofoils |
US7793884B2 (en) * | 2008-12-31 | 2010-09-14 | Faruk Dizdarevic | Deltoid main wing aerodynamic configurations |
WO2010129210A2 (en) * | 2009-04-27 | 2010-11-11 | University Of Miami | Supersonic flying wing |
US8448893B2 (en) * | 2009-10-26 | 2013-05-28 | Aerion Corporation | Laminar flow wing optimized for transonic cruise aircraft |
CN203306221U (zh) * | 2012-09-13 | 2013-11-27 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种兼顾增升结构的超音速飞机机翼翼型 |
CN104118556B (zh) * | 2014-08-07 | 2015-03-18 | 西北工业大学 | 一种极低雷诺数高升阻比低速特殊勺型翼型 |
CN204310031U (zh) * | 2014-12-05 | 2015-05-06 | 成都飞机设计研究所 | 一种模块化飞翼布局 |
CN106043688A (zh) * | 2016-06-08 | 2016-10-26 | 南京航空航天大学 | 一种直升机旋翼翼型 |
CN106335624B (zh) * | 2016-08-30 | 2018-07-03 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种采用双尖锥鱼尾型前缘气动布局的飞行器 |
CN106864726B (zh) * | 2016-12-27 | 2019-02-26 | 西北工业大学 | 一种适用于通用飞行器偏重巡航特性的层流翼型 |
US11126758B2 (en) * | 2018-03-22 | 2021-09-21 | Facebook, Inc. | Automatic airfoil and wing design based on dynamic modeling of structural and aerodynamic performance |
EP3543108A1 (en) * | 2018-03-22 | 2019-09-25 | Facebook, Inc. | Automatic airfoil and wing design based on dynamic modeling of structural and aerodynamic performance |
CN109204777B (zh) * | 2018-10-31 | 2023-12-15 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种直升机翼型 |
CN109808913B (zh) * | 2019-01-29 | 2022-05-17 | 西北工业大学 | 一种带有可偏转翼梢小翼的无人机设计方法 |
CN109808869A (zh) * | 2019-03-10 | 2019-05-28 | 西北工业大学 | 一种超声速环形翼设计方法 |
CN109969381A (zh) * | 2019-04-12 | 2019-07-05 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种具有尖前缘大弯度的低雷诺数高升阻比翼型 |
CN111498084A (zh) * | 2020-04-15 | 2020-08-07 | 成都飞机工业(集团)有限责任公司 | 一种应用于高空高速长航时无人机的低阻层流翼型 |
CN111611652B (zh) * | 2020-06-04 | 2023-06-09 | 电子科技大学 | 一种基于主动柔性变形的扑翼气动特性分析方法 |
CN112329140B (zh) * | 2020-10-30 | 2022-07-15 | 北京理工大学 | 基于改进位置向量期望改善度的变体飞行器气动优化方法 |
CN113962025B (zh) * | 2021-10-26 | 2024-05-14 | 成都飞机工业(集团)有限责任公司 | 超扁平无尾布局超音速飞机的翼型优化方法及装置 |
CN114021265A (zh) * | 2021-12-08 | 2022-02-08 | 辽宁科技大学 | 一种基于对数螺旋线的飞机用翼型 |
CN114658595A (zh) * | 2022-03-22 | 2022-06-24 | 华能湖北新能源有限责任公司 | 一种具有前缘延伸增功装置的风力机叶片及其设计方法 |
CN115092414B (zh) * | 2022-07-05 | 2024-02-02 | 中国人民解放军国防科技大学 | 一种环量控制翼型气动与电磁隐身联合优化方法 |
CN114936430B (zh) * | 2022-07-25 | 2022-10-25 | 西北工业大学 | 一种宽速域高超声速飞行器气动布局设计方法及系统 |
CN115320827B (zh) * | 2022-10-14 | 2023-01-31 | 中国航空工业集团公司沈阳空气动力研究所 | 一种高亚音速低雷诺数流动的高升阻比翼型 |
-
2022
- 2022-12-14 CN CN202211601613.1A patent/CN115649417B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105752314A (zh) * | 2016-03-22 | 2016-07-13 | 西北工业大学 | 一种高空低速自然层流高升力翼型 |
CN114684348A (zh) * | 2022-04-01 | 2022-07-01 | 西北工业大学 | 一种适用于轻中型高速公务机的高阻力发散马赫数层流翼型 |
Also Published As
Publication number | Publication date |
---|---|
CN115649417A (zh) | 2023-01-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7092421B2 (ja) | ヘリコプター回転翼の翼型の決定方法及びシステム | |
Sasaki et al. | Navier-Stokes optimization of supersonic wings with four objectives using evolutionary algorithm | |
CN112572761B (zh) | 一种具有低阻高发散马赫数高升力的跨空域稳健层流翼型 | |
US20090084904A1 (en) | Wingtip Feathers, Including Paired, Fixed Feathers, and Associated Systems and Methods | |
US20120166148A1 (en) | Method of designing natural laminar flow wing for reynolds numbers equivalent to actual supersonic aircraft | |
CN114528634A (zh) | 高隐身高机动布局飞机的弹性机翼气动隐身优化设计方法 | |
CN110104164B (zh) | 一种用于跨声速机翼的前加载-吸气组合流动控制方法 | |
CN115649417B (zh) | 一种高亚声速自配平高隐身翼型 | |
CN114154434A (zh) | 一种无尾飞翼布局的多约束精细化气动优化设计方法 | |
CN116443245B (zh) | 一种基于前行桨叶概念的高速特性显著提高的超临界旋翼翼型 | |
CN109677630B (zh) | 基准流场激波形状可控的强几何约束下的乘波体设计方法 | |
CN116306044A (zh) | 全湍流构型的不确定性分析方法及其梯度优化设计方法 | |
CN115618499A (zh) | 一种考虑三维影响的宽速域薄翼型多目标优化设计方法 | |
CN111498084A (zh) | 一种应用于高空高速长航时无人机的低阻层流翼型 | |
CN115092414B (zh) | 一种环量控制翼型气动与电磁隐身联合优化方法 | |
CN114906343A (zh) | 一种适用于扁平融合体布局飞行器的腹部襟翼设计方法 | |
CN105205208A (zh) | 一种飞机机翼型架外形设计方法 | |
CN113626935B (zh) | 一种高巡航效率跨音速月形机翼设计方法 | |
CN111597741B (zh) | 一种基于改进Hicks-henne算法的仿生蟹滑翔姿态下翼型优化设计方法 | |
CN211076318U (zh) | 一种用于机身减阻及延缓气流分离的涡流发生器 | |
CN115593612B (zh) | 一种自配平抗失速高性能翼型 | |
CN209852562U (zh) | 一种翼型及无人机 | |
CN117236231B (zh) | 融合工况与几何不确定性的翼型梯度优化设计方法 | |
CN113536458B (zh) | 一种基于类别形状函数变换的可变形机翼的快速几何建模方法 | |
CN112948976B (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 |