CN115092414A - 一种环量控制翼型气动与电磁隐身联合优化方法 - Google Patents
一种环量控制翼型气动与电磁隐身联合优化方法 Download PDFInfo
- Publication number
- CN115092414A CN115092414A CN202210785250.5A CN202210785250A CN115092414A CN 115092414 A CN115092414 A CN 115092414A CN 202210785250 A CN202210785250 A CN 202210785250A CN 115092414 A CN115092414 A CN 115092414A
- Authority
- CN
- China
- Prior art keywords
- airfoil
- wing
- pneumatic
- stealth
- dimensional
- 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
- 238000005457 optimization Methods 0.000 title claims abstract description 66
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000008569 process Effects 0.000 claims abstract description 22
- 238000004088 simulation Methods 0.000 claims abstract description 5
- 238000013461 design Methods 0.000 claims description 15
- 238000009826 distribution Methods 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000004422 calculation algorithm Methods 0.000 claims description 5
- 230000002068 genetic effect Effects 0.000 claims description 5
- 125000004122 cyclic group Chemical group 0.000 claims 4
- 238000012795 verification Methods 0.000 abstract description 8
- 230000007246 mechanism Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000003050 experimental design method Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 238000009966 trimming Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
- B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64C—AEROPLANES; HELICOPTERS
- B64C3/00—Wings
- B64C3/10—Shape of wings
- B64C3/14—Aerofoil profile
-
- 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)
- Aviation & Aerospace Engineering (AREA)
- Manufacturing & Machinery (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本发明公开了一种环量控制翼型气动与电磁隐身联合优化方法，包括：获取翼型气动与电磁隐身联合优化的Pareto前沿解；在Pareto前沿解中选择气动性能占优、隐身性能占优以及气动隐身性能均衡的三个二维翼型；将选择的二维翼型作为飞翼布局飞行器外翼段的翼型1和翼型2，形成三种三维飞翼布局构型；对各三维飞翼布局构型进行三维的气动与隐身仿真分析，并根据三维飞翼布局构型的气动和隐身特性数据确定最终翼型选取。本发明应用于飞翼布局优化技术领域，提出了一种二维优化+三维验证的2.5D优化概念，在常规二维翼型优化基础上增加三维飞翼布局的验证流程，对优化结果继续优选，获得适用于三维飞翼布局的环量控制翼型。
Description
技术领域
本发明涉及飞翼布局优化技术领域，具体是一种环量控制翼型气动与电磁隐身联合优化方法。
背景技术
飞翼布局通常指没有尾翼且翼身融合，以整个翼身作为升力面的布局。采用飞翼布局的飞行器通常具有较高的升阻比和良好的电磁隐身性能，能够穿破对方防空系统，执行持久的监视和情报收集工作。
飞翼布局相比常规飞行器有着如下优势：其一，飞翼布局没有尾翼，翼身融合可以消除翼-尾涡流和激波的影响，减小飞行阻力；其二，整个机身作为升力面，翼面积较大因而翼载荷较低，即使是在空气稀薄的高空仍能提供足够升力，相比常规大展弦比高空长航飞行器减少了大翼展所需的结构重量；其三，取消尾翼的同时也取消掉了用于支撑尾翼的支撑结构，使得全机结构重量大大减轻，增大飞行器装载能力；最后，平滑的机身和取消尾翼消除了一些强散射源，使得全机的雷达反射面积显著减少，提高了飞行器的隐身性能。
飞翼布局飞行器的上述优势，使它受到世界范围内的广泛重视。但随着雷达预警探测系统的快速发展，飞翼布局使用传统舵面配平或机动时，仍将暴露于对方雷达之下，因而不能满足日益激烈的“探测与反探测”的生存需求。传统飞翼布局无人机的舵面偏转会破坏飞翼布局的平滑外形，同时，舵面偏转机构会在飞行器外表面上留下开孔和狭缝，形成电磁波反射源。为了消除舵面机构对飞行器隐身性能的影响，利用后缘环量控制设备替换舵面来进行无人机的飞行控制成为研究前沿。安装后缘环量控制设备的翼型被称为环量控制翼型(CCA)，其部件名称和参数定义如图1所示。
图1中，坐标系原点O固定在翼型前缘点，翼型的弦向和法向分别为x轴和y轴方向。环量控制翼型的弦长为在x轴方向从翼型前缘点到射流出口位置的距离，定义为c。后缘射流机构包含Coanda曲面、射流出口和高压气室。Coanda曲面的半径定义为R。气体经过高压气室出口(射流出口)形成射流，射流夹带边界层流动让翼型产生气动型弯曲，进而改变翼型环量和升力，环量控制的原理如图2所示。在两侧机翼后缘分别使用环量控制设备，可以进一步产生姿态控制力矩。
使用环量控制设备代替传统舵面虽让飞机外表面减少了开孔和缝隙，同时避免了偏转舵面对雷达波的散射。但环量控制设备的安装将增大机翼后缘厚度，必然增强后向电磁散射强度。机翼上下表面和Coanda曲面的几何形状决定了环量控制能力和翼型的电磁散射特性，因此有必要开展气动和隐身联合优化，同时提升翼型的气动和隐身特性。另一方面，环量控制翼型主要用于飞翼布局，其后掠角较大，机翼表面的展向流动显著。常规的二维翼型优化方法得到的优化翼型在三维飞翼布局上往往表现不佳。若直接开展三维飞翼布局整机优化，气动和隐身计算量又难以承受。
发明内容
针对上述现有技术中的不足，本发明提供一种环量控制翼型气动与电磁隐身联合优化方法，提出了二维优化+三维验证的2.5D优化概念，在常规二维翼型优化基础上增加三维飞翼布局的验证流程，对优化结果继续优选，获得适用于三维飞翼布局的环量控制翼型。
为实现上述目的，本发明提供一种环量控制翼型气动与电磁隐身联合优化方法，包括如下步骤：
步骤1，获取翼型气动与电磁隐身联合优化的Pareto前沿解；
步骤2，在Pareto前沿解中选择气动性能占优、隐身性能占优以及气动隐身性能均衡的三个二维翼型；
步骤3，将选择的二维翼型作为飞翼布局飞行器外翼段的翼型1和翼型2，并放样生成翼型1和翼型2之间的外翼段，形成三种三维飞翼布局构型；
步骤4，对各三维飞翼布局构型进行三维的气动与隐身仿真分析，并根据三维飞翼布局构型的气动和隐身特性数据确定最终翼型选取。
在其中一个实施例，步骤1的过程具体为：
步骤1.1，绘制初始CCA翼型网格，并在翼型上下表面以及Coanda曲面上设置若干控制点，进行翼型参数化建模；
步骤1.2，基于FFD的翼型几何和网格变形，在设计空间内选取若干样本点，并收集每个样本点的气动数据与隐身数据，构建翼型气动与电磁隐身联合优化的代理模型；
步骤1.3，基于代理模型，以最小化前、后向±20°RCS均值、最大化升阻比为优化目标，以设计空间与俯仰力矩系数绝对值小于预设值为约束条件，经过遗传算法优化，得到翼型气动与电磁隐身联合优化的Pareto前沿解。
在其中一个实施例，步骤1.1中，在翼型上下表面设置控制点的过程为：
在翼型上表面，沿弦向间隔设置若干第一控制点，其中，第一控制点在翼型前缘的分布密度大于在翼型后缘的分布密度；
在翼型下表面，沿弦向间隔设置若干第二控制点，其中，第二控制点在翼型前缘的分布密度大于在翼型后缘的分布密度。
在其中一个实施例，所述第一控制点的数量与所述第二控制点的数量相同且一一对应，相互对应的第一控制点与第二控制点位于翼型的同一法向上。
在其中一个实施例，所述第一控制点、所述第二控制点沿弦向分布位置为：
式中，x表示第n个第一控制点或第n个第二控制点到翼型前缘的距离，c为弦长，a、b为计算系数，N为第一控制点或第二控制点的数量，n＝1,2,…,N。
在其中一个实施例，步骤1.1中，在Coanda曲面上设置控制点的过程为：
在半圆形Coanda曲面上，设置若干沿周向均匀分布的第三控制点。
在其中一个实施例，步骤1.2与步骤1.3中，所述设计空间为：
对于翼型上下表面的控制点，其沿翼型法向的纵坐标变化范围为原始翼型厚度的20％；
对于Coanda曲面上的控制点，其沿Coanda曲面径向移动范围为20％R，其中，R为Coanda曲面的半径。
在其中一个实施例，步骤1.3中，所述俯仰力矩系数绝对值小于预设值，具体为：
俯仰力矩系数绝对值小于或等于0.01。
本发明提供的一种环量控制翼型气动与电磁隐身联合优化方法，具有如下有益技术效果：
1、基于FFD方法的翼型表面和Coanda曲面联合参数化建模：
对于翼型上下表面，FFD方法的控制点在翼型前缘分布较密，可更准确地控制小曲率翼型前缘的几何变形；对于半圆形Coanda曲面，控制点沿周向均匀分布，可更全面地控制整个Coanda曲面的几何变形；这样的联合参数化建模方式产生的优化设计空间范围更大，分布更加合理。
2、建立二维优化+三维验证的翼型2.5D优化流程：
常规翼型优化终止于二维优化结果，但包含后缘射流机构的翼型主要应用平台为飞翼布局飞行器，飞翼布局常常具有较大的后掠角，机翼展向流动带来显著的三维效应。常规二维翼型优化结果往往无法在飞翼布局飞行器上表现出最佳气动和电磁隐身性能。同时在飞行器上，射流与绕流的相互作用方式与二维翼型的情况差异较大，因此射流在增升或飞行控制方面的效果将与二维翼型优化结果不同。
本发明提出二维优化+三维验证的翼型2.5D优化流程，在常规二维翼型气动与电磁隐身多目标优化的基础上，增加在三维飞行器上的验证环节。在Pareto前沿上选择气动性能占优、隐身性能占优和气动隐身性能均衡的三个二维翼型，替换飞翼布局飞行器外翼段原始翼型，形成三种构型，进行三维的气动与隐身模拟。根据三种构型的气动与隐身性能数据开展对比优选工作，从而获得面向特定飞翼布局飞行器的翼型优化结果。整个优化流程结合了传统二维翼型优化流程与新的三维飞行器平台验证流程，因此可视为一种2.5D优化流程，能够保证翼型在飞行器上能够展示出气动和电磁隐身特性优势。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图示出的结构获得其他的附图。
图1为背景技术中环量控制翼型部件和参数定义的示意图；
图2为背景技术中环量控制的原理图；
图3为本发明实施例中翼型气动与电磁隐身联合优化的流程图；
图4为本发明实施例中初始CCA翼型网格示意图；
图5为本发明实施例中翼型变形控制点分布示意图；
图6为本发明实施例中翼型几何变形范围示意图；
图7为本发明实施例中CCA翼型几何和网格变形示意图，其中，(a)为变形前的示意图，(b)为变形后的示意图；
图8为本发明实施例中遗传算法优化解分布与Pareto前沿解的示意图；
图9为本发明实施例中优化翼型替换飞翼布局外翼段的示意图，其中，(a)为替换前的示意图，(b)为替换后的示意图。
本发明目的的实现、功能特点及优点将结合实施例，参照附图做进一步说明。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
需要说明，本发明实施例中所有方向性指示(诸如上、下、左、右、前、后……)仅用于解释在某一特定姿态(如附图所示)下各部件之间的相对位置关系、运动情况等，如果该特定姿态发生改变时，则该方向性指示也相应地随之改变。
另外，在本发明中如涉及“第一”、“第二”等的描述仅用于描述目的，而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中，“多个”的含义是至少两个，例如两个，三个等，除非另有明确具体的限定。
另外，本发明各个实施例之间的技术方案可以相互结合，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在，也不在本发明要求的保护范围之内。
如图3所示为本实施例公开的一种环量控制翼型气动与电磁隐身联合优化方法，其主要包括如下步骤1-4。
步骤1，获取翼型气动与电磁隐身联合优化的Pareto前沿解，其具体实施过程为：
步骤1.1，利用网格绘制软件Pointwise生成图4所示的初始CCA翼型网格，然后在几何和网格变形软件Sculptor中，在翼型上下表面以及Coanda曲面上设置若干控制点，进行翼型参数化建模。
本实施例中，在翼型上下表面以及Coanda曲面上设置若干控制点的过程为：
在翼型上表面，沿弦向间隔设置若干第一控制点，其中，第一控制点在翼型前缘的分布密度大于在翼型后缘的分布密度，即图5中的u1～u9；在翼型下表面，沿弦向间隔设置若干第二控制点，其中，第二控制点在翼型前缘的分布密度大于在翼型后缘的分布密度，即图5中的d1～d9；在半圆形Coanda曲面上，设置若干沿周向均匀分布的第三控制点，即图5中的ru1～ru4、r0、rd1～rd4。其中，控制点(u1～u9,r1～d9,ru1～ru4,r0,rd1～rd4)即作为优化设计变量，共27个。第一控制点与第二控制点在翼型前缘分布较密，可更准确地控制小曲率前缘的几何变形。第三控制点沿Coanda曲面周向均匀分布可更全面地控制整个Coanda曲面的几何变形。该联合参数化建模方式产生的优化设计空间范围更大，分布更加合理。
需要注意的是，虽然本实施例中图示的第一控制点、第二控制点、第三控制点的数量均为9个，但在实际应用过程中根据需要自适应的调整各控制点的数量。
本实施例中，第一控制点与第二控制点一一对应，且相互对应的第一控制点与第二控制点位于翼型的同一法向上。例如，第一控制点、第二控制点沿弦向分布位置为：
式中，x表示第n个第一控制点或第n个第二控制点到翼型前缘的距离，c为弦长，a＝1.0392、b＝0.0372为计算系数，N为第一控制点或第二控制点的数量，n＝1,2,…,N。通过该公式可以将有限的控制点(控制点多了，求解空间的维度就大，不容易找到最优解，陷入“维度灾难”。)多分给前缘，后缘少分配控制点，进而有效地减少优化问题的求解维度。
步骤1.2，基于FFD的翼型几何和网格变形，采用实验设计法(DOE)在设计空间内选取若干样本点，并收集每个样本点的气动数据与隐身数据，构建翼型气动与电磁隐身联合优化的代理模型。
参考图6，本实施例中的设计空间具体为：对于翼型上下表面的控制点，其沿翼型法向的纵坐标变化范围为原始翼型厚度的20％，翼型表面几何变化范围(即翼型表面的设计空间)被最大厚度和最小厚度翼型表面包络；对于Coanda曲面上的控制点，Coanda表面控制点r0沿Coanda曲面径向移动范围为20％R,其它控制点按照比例变化，令Coanda曲面几何变化范围(即Coanda曲面的设计空间)被最大椭圆Coanda曲面和最小椭圆Coanda曲面包络。
参考图7，对于所选取的每个样本点，可以通过在初始CCA翼型网格的基础上移动控制点得到，依靠自由曲面变形技术(FFD)实现翼型几何和网格整体变形，从而实现对CCA翼型的参数化表示。本实施例中代理模型的构建过程为：
变形后的网格文件可直接用于气动计算，调用Fluent开展气动计算，完成后可计算得到样本点对于翼型的升阻比K和俯仰力矩系数Cm，即得到样本点的气动数据；
变形后的CCA翼型几何文件，使用Hypermesh绘制网格用于隐身计算。使用FEKO可计算得到样本点对于翼型的前向、后向±20°RCS，即得到样本点的隐身数据；
在Isight中构建优化流程，收集每个样本点的气动和隐身数据，即可构建翼型气动与电磁隐身联合优化的代理模型。
步骤1.3，基于代理模型，以最小化前、后向±20°RCS均值
最大化升阻比K为优化目标，以设计空间与俯仰力矩系数Cm的绝对值小于或等于0.01为约束条件，经过遗传算法优化，得到翼型气动与电磁隐身联合优化的Pareto前沿解，即图8所示。
需要注意的是，关于气动数据、隐身数据的计算过程，以及Isight构建代理模型的过程与遗传算法的优化过程均为所属领域的常规技术手段，因此本实施例中不再对其进行赘述。
步骤2，在Pareto前沿解中选择气动性能占优、隐身性能占优以及气动隐身性能均衡的三个二维翼型，其中，气动性能占优指的是Pareto前沿解中升阻比K最大的优化解，即图8找那个点C代表的翼型C；隐身性能占优指的是Pareto前沿解中前、后向±20°RCS均值
最小的优化解，即图8找那个点A代表的翼型A；气动隐身性能均衡指的是Pareto前沿解中点C与电A连线中间区域的优化解，即图8找那个点B代表的翼型B。与对初始翼型相比，优化翼型A、优化翼型B、优化翼型C的优化结果如表1-2所示：
表1优化翼型与初始翼型气动特性参数对比
表2优化翼型与初始翼型RCS均值对比
从表1-2可知，相对原始翼型，优化翼型A、优化翼型B、优化翼型C的升阻比有显著提升，RCS明显减小。
步骤3，将选择的二维翼型作为飞翼布局飞行器外翼段的翼型1和翼型2(即将同一二维翼型分别作为翼型1和翼型2)，并放样生成翼型1和翼型2之间的外翼段，形成三种三维飞翼布局构型，即图9所示；
步骤4，对各三维飞翼布局构型进行三维的气动与隐身仿真分析，并收集气动与隐身数据如表3-4所示：
表3优化翼型与初始翼型应用于飞翼布局的气动数据
表4优化翼型与初始翼型应用于飞翼布局的隐身数据
各三维飞翼布局构型的气动和隐身特性数据确定最终翼型选取。例如，根据表3-4中三种优化翼型和初始翼型的飞翼布局气动和隐身特性数据可以发现，优化翼型A和优化翼型B布局的升阻比小于原始翼型布局。因此二维翼型优化结果应用到三维飞行器上并不一定能够展示出优越性。但优化翼型C布局的升阻比显著大于初始翼型布局，同时其各个方位的RCS均小于初始翼型。因此翼型C可被优选为最终翼型选取。需要注意的是，若是出现三种优化翼型都没有同时满足升阻比最大、RCS最小的情况，则可根据当前翼型设计的应用需求选择气动性能占优或隐身性能占优或气动隐身性能均衡的优化翼型作为最终翼型选取，也可以在Pareto前沿中继续选取对应的二维翼型，拓展至三维翼型中后，分析气动与隐身性能，并作出取舍。
以上所述仅为本发明的优选实施例，并非因此限制本发明的专利范围，凡是在本发明的发明构思下，利用本发明说明书及附图内容所作的等效结构变换，或直接/间接运用在其他相关的技术领域均包括在本发明的专利保护范围内。
Claims (8)
1.一种环量控制翼型气动与电磁隐身联合优化方法，其特征在于，包括如下步骤：
步骤1，获取翼型气动与电磁隐身联合优化的Pareto前沿解；
步骤2，在Pareto前沿解中选择气动性能占优、隐身性能占优以及气动隐身性能均衡的三个二维翼型；
步骤3，将选择的二维翼型作为飞翼布局飞行器外翼段的翼型1和翼型2，并放样生成翼型1和翼型2之间的外翼段，形成三种三维飞翼布局构型；
步骤4，对各三维飞翼布局构型进行三维的气动与隐身仿真分析，并根据三维飞翼布局构型的气动和隐身特性数据确定最终翼型选取。
2.根据权利要求1所述的环量控制翼型气动与电磁隐身联合优化方法，其特征在于，步骤1的过程具体为：
步骤1.1，绘制初始CCA翼型网格，并在翼型上下表面以及Coanda曲面上设置若干控制点，进行翼型参数化建模；
步骤1.2，基于FFD的翼型几何和网格变形，在设计空间内选取若干样本点，并收集每个样本点的气动数据与隐身数据，构建翼型气动与电磁隐身联合优化的代理模型；
步骤1.3，基于代理模型，以最小化前、后向±20°RCS均值、最大化升阻比为优化目标，以设计空间与俯仰力矩系数绝对值小于预设值为约束条件，经过遗传算法优化，得到翼型气动与电磁隐身联合优化的Pareto前沿解。
3.根据权利要求2所述的环量控制翼型气动与电磁隐身联合优化方法，其特征在于，步骤1.1中，在翼型上下表面设置控制点的过程为：
在翼型上表面，沿弦向间隔设置若干第一控制点，其中，第一控制点在翼型前缘的分布密度大于在翼型后缘的分布密度；
在翼型下表面，沿弦向间隔设置若干第二控制点，其中，第二控制点在翼型前缘的分布密度大于在翼型后缘的分布密度。
4.根据权利要求3所述的环量控制翼型气动与电磁隐身联合优化方法，其特征在于，所述第一控制点的数量与所述第二控制点的数量相同且一一对应，相互对应的第一控制点与第二控制点位于翼型的同一法向上。
5.根据权利要求3所述的环量控制翼型气动与电磁隐身联合优化方法，其特征在于，所述第一控制点、所述第二控制点沿弦向分布位置为：
式中，x表示第n个第一控制点或第n个第二控制点到翼型前缘的距离，c为弦长，a、b为计算系数，N为第一控制点或第二控制点的数量，n＝1,2,…,N。
6.根据权利要求2所述的环量控制翼型气动与电磁隐身联合优化方法，其特征在于，步骤1.1中，在Coanda曲面上设置控制点的过程为：
在半圆形Coanda曲面上，设置若干沿周向均匀分布的第三控制点。
7.根据权利要求2至6任一项所述的环量控制翼型气动与电磁隐身联合优化方法，其特征在于，步骤1.2与步骤1.3中，所述设计空间为：
对于翼型上下表面的控制点，其沿翼型法向的纵坐标变化范围为原始翼型厚度的20％；
对于Coanda曲面上的控制点，其沿Coanda曲面径向移动范围为20％R，其中，R为Coanda曲面的半径。
8.根据权利要求2至6任一项所述的环量控制翼型气动与电磁隐身联合优化方法，其特征在于，步骤1.3中，所述俯仰力矩系数绝对值小于预设值，具体为：
俯仰力矩系数绝对值小于或等于0.01。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210785250.5A CN115092414B (zh) | 2022-07-05 | 2022-07-05 | 一种环量控制翼型气动与电磁隐身联合优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210785250.5A CN115092414B (zh) | 2022-07-05 | 2022-07-05 | 一种环量控制翼型气动与电磁隐身联合优化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115092414A true CN115092414A (zh) | 2022-09-23 |
| CN115092414B CN115092414B (zh) | 2024-02-02 |
Family
ID=83297069
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210785250.5A Active CN115092414B (zh) | 2022-07-05 | 2022-07-05 | 一种环量控制翼型气动与电磁隐身联合优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115092414B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115649417A (zh) * | 2022-12-14 | 2023-01-31 | 中国空气动力研究与发展中心空天技术研究所 | 一种高亚声速自配平高隐身翼型 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060226295A1 (en) * | 2005-03-25 | 2006-10-12 | Gulick Douglas S | Lifting body tuned for passive re-entry |
| US20190031322A1 (en) * | 2013-10-25 | 2019-01-31 | Ioannis Micros | Stealth design with multi-faceted dihedral planform and insufflation mechanism |
| CN111959816A (zh) * | 2020-07-15 | 2020-11-20 | 成都飞机工业（集团）有限责任公司 | 一种提高飞翼布局飞机高低速性能的气动设计方法 |
| CN113415412A (zh) * | 2021-06-25 | 2021-09-21 | 中国人民解放军国防科技大学 | 宽速域射流控制飞行器 |
| CN114528634A (zh) * | 2021-12-31 | 2022-05-24 | 中国航空工业集团公司沈阳飞机设计研究所 | 高隐身高机动布局飞机的弹性机翼气动隐身优化设计方法 |
-
2022
- 2022-07-05 CN CN202210785250.5A patent/CN115092414B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060226295A1 (en) * | 2005-03-25 | 2006-10-12 | Gulick Douglas S | Lifting body tuned for passive re-entry |
| US20190031322A1 (en) * | 2013-10-25 | 2019-01-31 | Ioannis Micros | Stealth design with multi-faceted dihedral planform and insufflation mechanism |
| CN111959816A (zh) * | 2020-07-15 | 2020-11-20 | 成都飞机工业（集团）有限责任公司 | 一种提高飞翼布局飞机高低速性能的气动设计方法 |
| CN113415412A (zh) * | 2021-06-25 | 2021-09-21 | 中国人民解放军国防科技大学 | 宽速域射流控制飞行器 |
| CN114528634A (zh) * | 2021-12-31 | 2022-05-24 | 中国航空工业集团公司沈阳飞机设计研究所 | 高隐身高机动布局飞机的弹性机翼气动隐身优化设计方法 |
Non-Patent Citations (7)
| Title |
|---|
| 夏露;张欣;杨梅花;米百刚;: "飞翼布局翼型气动隐身综合设计", 西北工业大学学报, no. 05, pages 821 - 826 * |
| 孙奕捷;申功璋;: "飞翼布局飞机控制/气动/隐身多学科优化设计", 北京航空航天大学学报, no. 11, pages 1357 - 1360 * |
| 张彬乾: "飞翼布局隐身翼型优化设计", 航空学报, vol. 35, no. 4, pages 957 - 967 * |
| 樊华羽;詹浩;程诗信;米百刚;姚会勤;: "高效粒子群算法研究及飞翼无人机气动隐身优化设计", 航空工程进展, no. 06, pages 735 - 743 * |
| 王荣;闫溟;白鹏;杨云军;徐国武;: "飞翼无人机平面外形气动隐身优化设计", 航空学报, no. 1, pages 1 - 8 * |
| 陈保;白俊强;黎明;: "基于分解策略的飞行器气动隐身优化设计研究", 气体物理, no. 06, pages 40 - 49 * |
| 陈曦;白俊强;李权;: "某飞翼布局隐身飞行器的翼型优化", 航空计算技术, no. 06, pages 46 - 53 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115649417A (zh) * | 2022-12-14 | 2023-01-31 | 中国空气动力研究与发展中心空天技术研究所 | 一种高亚声速自配平高隐身翼型 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115092414B (zh) | 2024-02-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109895992B (zh) | 设计为与偏转状态下的襟翼密封的预变形飞机扰流板和下垂板 | |
| JP2018127202A (ja) | 翼フラップ偏位制御除去 | |
| CN115092414A (zh) | 一种环量控制翼型气动与电磁隐身联合优化方法 | |
| CN111278728A (zh) | 用于飞行器的圆顶整流罩及其制造方法 | |
| Elelwi et al. | Wing component allocation for a morphing variable span of tapered wing using finite element method and topology optimisation–application to the UAS-S4 | |
| Marks et al. | Investigation into the effect of shape deviation on variable camber compliant wing performance | |
| CN105117541B (zh) | 一种正向型架外形优化设计方法 | |
| KR101472388B1 (ko) | 외부 무장에 따른 경무장 헬기의 공력특성 및 성능예측장치 | |
| CN113361017B (zh) | 一种考虑机构约束的管射无人机二次折叠翼气动设计方法 | |
| Chekkal et al. | Design of a morphing wing tip | |
| Deslich et al. | Effects of spanloading and slew angle on an Oblique Flying Wing | |
| Schmitz et al. | Acoustic Phasing, Directionality and Amplification Effects of Helicopter Blade‐Vortex Interactions | |
| Obayashi et al. | Transonic wing design by inverse optimization using MOGA | |
| McAndrew et al. | Twin-wing design options used for unmanned aerial vehicles to achieve high altitudes at low speeds | |
| CN117634334A (zh) | 兼顾气动/隐身的战斗机翼型优化设计方法及宽速域翼型族 | |
| Lange et al. | Wing Shape Improvement in the Presence of an Over-the-Wing Mounted UHBR Engine for a Short Range Transport Aircraft | |
| Ding et al. | Structural design and optimization of a morphing wing trailing edge flap | |
| CN113536458B (zh) | 一种基于类别形状函数变换的可变形机翼的快速几何建模方法 | |
| CN113626934B (zh) | 飞翼布局无人机机翼多构型协调设计方法 | |
| EP1694560B1 (en) | Area ruling for supersonic vehicles | |
| ŞAHİN et al. | Experimental and numerical investigation of different wing profiles | |
| CN113665837B (zh) | 基于前缘线等激波强度尖头冯卡门曲线舵机鼓包设计方法 | |
| Jestus | Aerodynamic Characterization of Multiple Wing-Wing Interactions for Distributed Lift Applications | |
| Krenz | Transonic configuration design | |
| Gabor et al. | Aerodynamic analysis of upper surface wing morphing efficiency for the S4 Éhecatl unmanned aerial system |
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 |