一种用于低雷诺数翼型优化的方法技术

一种用于低雷诺数翼型优化的方法，根据任务要求，在Matlab软件上编译CST函数，用CST函数将翼型参数化，输入CST函数中的参数，得到CST函数对应的翼型的几何外形，将翼型的几何外形导入ICEM软件，利用ICEM宏录制功能，完成对翼型的有限元网格划分，得到网格文件；在Fluent软件中进行宏录制即对飞行工况、气动计算条件进行设置，将网格文件导入Fluent软件，完成气动计算，得到翼型升力系数、升阻比；对步骤(3)得到的翼型升力系数、升阻比进行判断，若翼型升力系数、升阻比均满足任务要求，则保留翼型的几何外形对应的CST函数；否则调整CST函数中的参数；得到了气动效率高的翼型，提高了飞行器的飞行性能。

A Method for Airfoil Optimization with Low Reynolds Number

A method for airfoil optimization with low Reynolds number is presented. According to the task requirement, CST function is compiled on Matlab software. The airfoil is parameterized by CST function, and the parameters in CST function are input. The airfoil geometry corresponding to CST function is obtained. The airfoil geometry is imported into ICEM software. The finite element meshing of airfoil is completed by using ICEM macro recording function, and the mesh file is obtained. Macro recording in Fluent software is to set the flight conditions and aerodynamic calculation conditions, import grid files into Fluent software, complete aerodynamic calculation, and get the lift coefficient and lift-drag ratio of airfoil; Judge the lift coefficient and lift-drag ratio of airfoil obtained in step (3). If the lift coefficient and lift-drag ratio of airfoil meet the mission requirements, the corresponding C of airfoil geometry will be retained. ST function; otherwise adjust the parameters of CST function; get the airfoil with high aerodynamic efficiency, improve the flight performance of aircraft.

【技术实现步骤摘要】
一种用于低雷诺数翼型优化的方法
本专利技术涉及一种用于低雷诺数翼型优化的方法，属于临近空间飞行器领域。
技术介绍
美国在上世纪90年代先后推出“捕食者”、“全球鹰”无人机，我国也紧跟步伐，连续研制出“彩虹”、“翼龙”等型号无人机，卓有成效地完成了包括打击、侦察和预警等多种军事任务。就目前的发展趋势来看，高空、长航时、多用途是未来无人机的主要发展方向。高空无人机的巡航高度位于平流层，该区域距离地面一般在20-30km，空气稀薄，密度低，运动粘性系数较大，因此飞行雷诺数较常规而言更低(Re≈2×105)，所以提升飞机巡航时的升力系数和升阻比是至关重要的。在目前的技术条件下，高空长航时无人机的飞行速度一般较大，如果强行降低巡航速度，飞机会有失速的危险。
技术实现思路
本专利技术解决的技术问题为：克服现有技术不足，提供一种用于低雷诺数翼型优化的方法，通过翼型优化，实现了高空无人机的巡航高度位于平流层，该区域距离地面一般在20-30km，空气稀薄，密度低，运动粘性系数较大的情况下，提升飞机巡航时的升力系数和升阻比，实现高空无人机的巡航作业。本专利技术解决的技术方案为：一种用于低雷诺数翼型优化的方法，步骤如下：(1)根据任务要求，在Matlab软件上编译CST函数，用CST函数将翼型参数化，输入CST函数中的参数，得到CST函数对应的翼型的几何外形；(2)将翼型的几何外形导入ICEM软件，利用ICEM宏录制功能，完成对翼型的有限元网格划分，得到网格文件；(3)在Fluent软件中进行宏录制即对飞行工况、气动计算求解方程进行设置，将网格文件导入Fluent软件，完成气动计算，得到翼型巡航时升力系数、升阻比；(4)对步骤(3)得到的翼型巡航时升力系数、升阻比进行判断，若翼型巡航时升力系数、升阻比均满足任务要求，则保留步骤(1)中翼型的几何外形对应的CST函数；否则，返回步骤(1)调整CST函数中的参数。几何外形，包括翼型的厚度、弯度和前缘半径。编译的CST函数，包括类别函数和形状函数两部分相乘得到CST函数，CST函数表示的曲线能够准确描述低雷诺数翼型的外形。输入CST函数中的参数，具体为：CST函数中的部分中的参数输入后不变，S(ψ)中的参数能够调整。利用ICEM宏录制功能，完成对翼型的有限元网格划分，具体为：网格划分时采用全结构网格，设定翼型的远场前方、上方、下方距离、后方距离、网格高度。网格文件为.mesh格式。对飞行工况、气动计算求解方程进行设置，具体为：飞行工况设置包括：任务要求的飞行高度、飞行速度、飞行姿态；气动计算求解方程包括模拟翼型层流流动、转捩预测、湍流流动的方程。完成气动计算，得到翼型升力系数、升阻比，具体为：利用模拟翼型层流流动、转捩预测、湍流流动的方程，得到翼型巡航时升力系数、升阻比。步骤(4)对步骤(3)得到的翼型巡航时升力系数、升阻比进行判断时，任务要求包括巡航升力系数的最小值和升阻比的最小值，步骤(3)得到的翼型巡航时升力系数分别大于任务要求的巡航升力系数的最小值和升阻比的最小值判定为满足任务要求。步骤(4)返回步骤(1)调整CST函数中的参数，具体为：将CST函数中的部分中的参数保持不变，调整中的参数。本专利技术与现有技术相比的优点在于：(1)本专利技术一种用于低雷诺数翼型优化的方法，翼型几何建模简单，有利于优化输入参数的提取；(2)本专利技术一种用于低雷诺数翼型优化的方法，在优化时充分利用计算流体力学分析软件，计算结果精度高，可直接作为工程设计的参考数据；(3)本专利技术一种用于低雷诺数翼型优化的方法，优化目标明确，迭代周期短，可缩短工程研制周期。(4)本专利技术创新性的将翼型建模、网格划分、气动计算搭建在ISIGHT平台上，根据优化目标，实现了翼型的自动选取。(5)本专利技术通过优化设计出的翼型升力系数和升阻比高，有利于高空无人机的低速巡航，增长飞机的留空时间。附图说明图1是E387翼型优化流程示意图；图2是CST参数化E387翼型示意图；图3是CST参数化E387翼型误差示意图；图4是ICEM划分网格示意图；图5是翼型优化前后外形对比图；图6是翼型优化前后升力系数对比图；图7是翼型优化前后升阻比对比图。具体实施方式下面结合附图和具体实施例对本专利技术做进一步详细描述。本专利技术优化原始翼型为E387翼型，该翼型为典型低雷诺数翼型，翼型主体呈前圆后尖的形状，最大厚度为弦长的9％，最大弯度为弦长的3.7％，前缘半径为弦长的1.6％。(1)根据任务要求，在Matlab软件上编译CST函数，用CST函数将翼型参数化，输入CST函数中的参数，得到CST函数对应的翼型的几何外形。(2)将翼型的几何外形导入ICEM软件，利用ICEM宏录制功能，完成对翼型的有限元网格划分，得到网格文件；(3)在Fluent软件中进行宏录制即对飞行工况、气动计算求解方程进行设置，将网格文件导入Fluent软件，完成气动计算，得到翼型巡航时升力系数、升阻比；(4)对步骤(3)得到的翼型巡航时升力系数、升阻比进行判断，若翼型巡航时升力系数、升阻比均满足任务要求，则保留步骤(1)中翼型的几何外形对应的CST函数；否则，返回步骤(1)调整CST函数中的参数；在参数化翼型时，需要通过CST函数准确描述翼型的几何外形特征，包括翼型的厚度、弯度和前缘半径。该函数包括类别函数和形状函数S(ψ)，两部分相乘得到CST函数，CST函数表示的曲线能够准确描述低雷诺数翼型的外形。类别函数其中，N1，N2为翼型几何外形类别函数的特征参数，优选N1为0.5，优选N2为1，定义原点在翼型的前缘的最前端，x轴从翼型的前缘的最前端指向翼型尾缘的中心；，Ψ为翼型横坐标(横坐标记为x轴的坐标)，n优选取4,B(i)为控制翼型外形曲率的参数,控制翼型外形曲率的参数,以优化的基准翼型(即E387翼型)为例，B(0)取0.306、B(1)取0.123、B(2)取0.457B(2)取0.277(B3麻烦帮我改一下，我这里编辑不了)和B(4)取0.454,为函数的n阶多项式展开式，通过以上公式的表述，可以清晰的理解CST参数化二维翼型的基本原理。i＝0,1,2,..,n；在整体优化流程设计时，既要考虑翼型计算过程中每一步的计算精度，同时要兼顾计算周期，因此，本专利通过搭建自动优化平台实现了翼型参数建模、有限元分析以及比对筛选的功能，该专利可以根据任务要求，快速完成迭代，实现优化目标。如图1所述，本专利技术优选以E387低雷诺数翼型作为基础翼型，参数化翼型与原始翼型的对比图和误差图如图2、图3所示，图2横坐标为翼型的横坐标数值与弦长的比值，纵坐标为翼型的纵坐标数值与弦长的比值，图3的横坐标为翼型的横坐标数值与弦长的比值，总坐标为CST函数拟合误差，由图中数据可以看到，CST方法拟合的翼型与原翼型曲线基本吻合，误差精度控制在10-3量级，满足工程设计需要；第二部分为ICEM网格自动生成模块，如图4所示，网格划分时采用全结构网格，远场前方、上方、下方距离为30倍弦长，后方距离为40倍弦长，第一层网格高度为10-5倍弦长，y+小于1，该部分主要完成翼型点、线、面映射的任务，实现自动划分不同翼型网格的功能；第三部分是Fluent计算宏录制，该部分主要完成翼型气动力计算的工作，在进行宏录制时要设定好飞行工况以本文档来自技高网...

【技术保护点】
1.一种用于低雷诺数翼型优化的方法，其特征在于步骤如下：(1)根据任务要求，在Matlab软件上编译CST函数，用CST函数将翼型参数化，输入CST函数中的参数，得到CST函数对应的翼型的几何外形；(2)将翼型的几何外形导入ICEM软件，利用ICEM宏录制功能，完成对翼型的有限元网格划分，得到网格文件；(3)在Fluent软件中进行宏录制即对飞行工况、气动计算求解方程进行设置，将网格文件导入Fluent软件，完成气动计算，得到翼型巡航时升力系数、升阻比；(4)对步骤(3)得到的翼型巡航时升力系数、升阻比进行判断，若翼型巡航时升力系数、升阻比均满足任务要求，则保留步骤(1)中翼型的几何外形对应的CST函数；否则，返回步骤(1)调整CST函数中的参数。

【技术特征摘要】
1.一种用于低雷诺数翼型优化的方法，其特征在于步骤如下：(1)根据任务要求，在Matlab软件上编译CST函数，用CST函数将翼型参数化，输入CST函数中的参数，得到CST函数对应的翼型的几何外形；(2)将翼型的几何外形导入ICEM软件，利用ICEM宏录制功能，完成对翼型的有限元网格划分，得到网格文件；(3)在Fluent软件中进行宏录制即对飞行工况、气动计算求解方程进行设置，将网格文件导入Fluent软件，完成气动计算，得到翼型巡航时升力系数、升阻比；(4)对步骤(3)得到的翼型巡航时升力系数、升阻比进行判断，若翼型巡航时升力系数、升阻比均满足任务要求，则保留步骤(1)中翼型的几何外形对应的CST函数；否则，返回步骤(1)调整CST函数中的参数。2.根据权利要求1所述的一种用于低雷诺数翼型优化的方法，其特征在于：几何外形，包括翼型的厚度、弯度和前缘半径。3.根据权利要求1所述的一种用于低雷诺数翼型优化的方法，其特征在于：编译的CST函数，包括类别函数和形状函数两部分相乘得到CST函数，CST函数表示的曲线能够准确描述低雷诺数翼型的外形。4.根据权利要求3所述的一种用于低雷诺数翼型优化的方法，其特征在于：输入CST函数中的参数，具体为：CST函数中的部分中的参数输入后不变，S(ψ)中的参数能够调整。5.根据权利要求1所述的一种用于低雷诺数翼型...

【专利技术属性】
技术研发人员：王超，姚爱霞，吴尚卿，丁志超，
申请(专利权)人：航天神舟飞行器有限公司，
类型：发明
国别省市：天津,12