一种基于非概率区间分析模型的翼型鲁棒优化设计方法技术

本发明专利技术公开了一种基于非概率区间分析模型的翼型鲁棒优化设计方法，属于优化设计技术领域。本发明专利技术充分考虑翼型设计中的不确定因素，在不确定参数概率密度未知的情况下，利用区间向量实现不确定参数的定量化表征。通过样本点试验，建立翼型设计变量、不确定参数及气动力系数之间的映射关系。在此基础上，利用非概率区间分析模型得到翼型气动力系数的区间上下界，建立鲁棒优化模型，进而应用遗传算法对翼型进行优化设计。数值预测表明，本发明专利技术方法在保持翼型升力系数以及几何形状约束的前提条件下，降低了所设计翼型的阻力系数，同时减小了阻力系数的变动范围，为翼型优化设计提供了新思路。

【技术实现步骤摘要】
一种基于非概率区间分析模型的翼型鲁棒优化设计方法
本专利技术涉及翼型优化设计领域，特别涉及一种基于非概率区间分析模型的翼型鲁棒优化设计方法。
技术介绍
对翼型进行气动优化时，存在确定性优化设计与鲁棒性优化设计两种区别较大的方法。确定性优化设计是在特定的飞行状态下，构造相应的翼型剖面，以最大限度提高翼型在该状态下的气动性能，是目前应用较广泛的翼型优化方法。然而，在实际巡航中，大气来流速度、气压、温度、突风等不确定因素都会使飞行器服役环境出现偏离设计预期的情况，从而使确定性优化方法得到的最优解对于这些不确定因素较为敏感。与确定性优化不同，鲁棒优化设计是一种寻求对各种不确定因素波动变化不敏感的设计方法，通过合理选择翼型参数，使翼型的气动性能稳定在目标值附近，提高翼型的可靠性。二维翼型鲁棒设计的实质是在满足约束条件的前提下，得到一个对不确定参数扰动并不敏感的翼型。对于翼型鲁棒设计中不确定参数的描述一般采用概率模型，利用随机变量对不确定参数定量化，然后通过MonteCarlo仿真等方法获得气动力系数的均值和方差。但是，在翼型设计中，有时很难获得较多的统计数据来描述不确定参数的概率分布，如大气来流速度等不确定参数的统计资料很少，不足以推断其分布于数字特征值。与概率模型相比，非概率区间分析模型仅需明确不确定参数的分布界限，能够在不确定参数的概率密度未知的情况下，对不确定参数进行定量化，并以此为基础计算响应所在的范围。非概率区间分析模型在结构的静、动力特性分析领域已经取得了一定成果，但在翼型外形优化中的应用还处于起步阶段，相关研究成果十分有限，在一定程度上限制了翼型优化设计技术的发展。
技术实现思路
本专利技术要解决技术问题为：针对传统翼型优化设计方法效率低、优化结果鲁棒性差等问题，提出一种基于非概率区间分析模型的翼型鲁棒优化设计方法。该方法在翼型设计初期就考虑了来流速度的不确定性，通过Kriging代理模型建立了包含翼型设计变量、不确定参数以及气动力系数之间的映射关系，利用非概率区间分析模型计算翼型气动力系数的上下界，建立一种针对二维翼型的非概率鲁棒优化设计方法。本专利技术解决上述技术问题采用的技术方案为：一种基于非概率区间分析模型的翼型鲁棒优化设计方法，包括以下步骤：(1)首先，通过类型函数/形状函数转换CST方法//KulfanBM.Universalparametricgeometryrepresentationmethod[J].JournalofAircraft,2008,45(1):142-158//设定参数化翼型曲线的表达式。翼型几何曲线可用下列函数表示：式中，x/c为翼型弦向的无量纲坐标值，y/c为翼型法向的无量纲坐标值，C(x/c)和S(x/c)为类型函数和形状函数，zte/c为翼型后缘点坐标，C(x/c)可表示为：对于翼型，取指数N1＝0.5,N2＝1，则S(x/c)可表示为：S(x/c)与翼型前缘半径Rle及后缘倾角β满足以下关系：在形状函数S(x/c)中取n＝3，提取其中8个参数作为翼型设计变量x，可表示为：x＝(x1,x2,…,x8)＝(Rle/c,β1,β2,zte/c,b1,b2,b1',b2')(5)其中，β1和β2表示上下翼面的后缘倾角，b1、b2、b1'、b2'为上下翼面形状函数多项式的加权系数。(2)设定翼型设计变量的区间上下界，其中下界记为：x＝(x1,x2,…,x8)，上界记为：(3)设定不确定设计参数的区间上下界，记为(4)采用均匀设计方法在由翼型设计变量和不确定参数组成的混合变量空间中生成样本点Ui；(5)对于每一个样本点Ui，根据样本点中的翼型设计变量，建立翼型的CAD模型；(6)由于翼型的CAD模型不能直接用于气动力系数计算，需要将其导入前处理软件GAMBIT中进行流场网格划分；(7)将步骤(6)中创建的网格文件导入FLUENT软件中，来流马赫数根据样本点Ui中的参数确定。采用S-A湍流模型、二阶迎风格式进行气动力系数的计算。(8)重复步骤(5)～(7)，完成所有样本点对应气动力系数的计算流程；(9)根据步骤(8)中的样本点输入参数及对应的气动力系数，通过Kriging代理模型建立翼型设计变量、不确定参数及气动力系数之间的映射关系；(10)选取翼型设计变量的初始值；(11)在步骤(9)中建立的Kriging模型基础上，引入非概率区间分析模型计算翼型气动力系数的区间上下界，进而得到翼型气动力系数的区间中心值和半径；(12)在保持翼型升力系数以及几何形状的约束条件下，以阻力系数的区间中心值和半径最小化为优化目标，采用遗传算法对翼型进行多目标鲁棒优化设计；(13)判断优化目标是否满足收敛条件，若不满足，转到步骤(10)，更新设计变量，重复步骤(11)～(12)；(14)直至设计目标相邻两次迭代值的变化小于设定容许偏差时，完成二维翼型的多目标鲁棒优化设计；(15)将鲁棒优化得到的翼型与初始翼型及进行传统确定性优化得到的翼型进行对比，分析比较三种翼型的几何形状及气动特性。其中，所述步骤(2)中，翼型设计变量的上下界由表2确定。表2翼型优化设计变量的取值范围其中，所述步骤(3)中，考虑不确定设计参数为来流马赫数Ma，其区间上下界为Ma∈[0.74,0.76]，中心值为Mac＝0.75，区间半径为ΔMa＝0.01。其中，所述步骤(9)中，通过Kriging模型建立了同时包含翼型设计变量、不确定参数及气动力系数之间的函数映射关系，即：其中，所述步骤(11)中，在不确定参数的概率密度未知的条件下，引入非概率区间分析模型，计算气动力系数的区间上下界，得到气动力系数的区间中心值和半径。根据区间数学理论，有界不确定参数向量属于某一区间向量，即：式中，和称为区间数的中心值和半径；δ＝[-1,1]。将由式(6)确定的气动响应函数在有界不确定参数的中心值处进行Taylor级数展开，可得：略去式(8)中的二阶及其以上的高阶小量，应用区间数学中的自然扩张原理，可以近似得到气动响应在不确定参数作用下的区间上下界：进一步计算可以得到气动力系数的区间中心值和半径为：其中，所述步骤(12)中，建立如下的多目标鲁棒优化模型：其中，表示对目标函数的区间中心值优化，而表示对不确定参数引起的偏差优化，从而降低目标函数对不确定参数的敏感程度，提高翼型的鲁棒性；为约束条件，n为约束条件的个数。本专利技术的有益效果是：本专利技术利用CST参数化建模方法实现翼型几何外形的快速表征，同时考虑翼型设计中存在的不确定参数，利用非概率区间分析模型计算翼型气动力系数的区间上下界，建立含不确定参数的多目标鲁棒优化模型，进而通过遗传算法对翼型进行鲁棒优化设计。与初始翼型和通过确定性优化得到翼型相比，本专利技术所得翼型的阻力系数显著减小，并且对于来流马赫数的敏感性显著降低，保证了翼型的安全可靠性。附图说明图1为翼型几何控制参数示意图；图2为翼型的CAD模型示意图；图3为翼型外流场结构网格示意图；图4为翼型几何形状对比图；图5为翼型阻力系数对比图；图6为本专利技术的方法实现流程图。具体实施方式以下将参照附图，对本专利技术的设计实例进行详细描述。应当理解，所选实例仅为了说明本专利技术，而不是限制本专利技术的保护范围。(1)以NACA0012翼型作为初始翼型，利用类型函数/本文档来自技高网...

【技术保护点】
一种基于非概率区间分析模型的翼型鲁棒优化设计方法，其特征在于实现步骤如下：(1)首先，利用类型函数/形状函数转换方法设定参数化翼型曲线的表达式，翼型的几何曲线可用下列函数表示：yc(xc)=C(xc)S(xc)+xcztec---(1)]]>式中，x/c为翼型弦向的无量纲坐标值，y/c为翼型法向的无量纲坐标值，C(x/c)和S(x/c)为类型函数和形状函数，zte/c为翼型后缘点坐标，C(x/c)可表示为：C(xc)=(xc)N1(1-xc)N2,0≤xc≤1---(2)]]>对于翼型，取指数N1＝0.5,N2＝1，则S(x/c)可表示为：S(xc)=Σi=0n[bi·n!i!(n-i)!·(xc)i·(1-xc)n-i],0≤xc≤1---(3)]]>S(x/c)与翼型前缘半径Rle及后缘倾角β满足以下关系：S(0)=b0=2Rlec,S(1)=bn=tanβ+ztec---(4)]]>在形状函数S(x/c)中取n＝3，提取其中8个参数作为翼型设计变量x，可表示为：x＝(x1,x2,…,x8)＝(Rle/c,β1,β2,zte/c,b1,b2,b1',b2')      (5)其中，β1和β2表示上下翼面的后缘倾角，b1、b2、b1'、b2'为上下翼面形状函数多项式的加权系数；(2)设定翼型设计变量的区间上下界，其中下界记为：x＝(x1,x2,…,x8)，上界记为：x‾=(x‾1,x‾2,...,x‾8);]]>(3)设定不确定设计参数的区间上下界，记为(4)采用均匀设计方法在由翼型设计变量和不确定参数组成的混合变量空间中生成样本点Ui；(5)对于每一个样本点Ui，根据样本点中的翼型设计变量，建立翼型的CAD模型；(6)由于翼型的CAD模型不能直接用于气动力系数计算，需要将其导入前处理软件GAMBIT中进行流场网格划分；(7)将步骤(6)中创建的网格文件导入FLUENT软件中，来流马赫数根据样本点Ui中的参数确定，采用S‑A湍流模型、二阶迎风格式进行气动力系数的计算；(8)重复步骤(5)～(7)，完成所有样本点对应气动力系数的计算流程；(9)根据步骤(8)中的样本点输入参数及对应的气动力系数，通过Kriging代理模型建立翼型设计变量、不确定参数及气动力系数之间的映射关系；(10)选取翼型设计变量的初始值；(11)在步骤(9)中建立的Kriging模型基础上，引入非概率区间分析模型计算翼型气动力系数的区间上下界，进而得到翼型气动力系数的区间中心值和半径；(12)在保持翼型升力系数以及几何形状的约束条件下，以阻力系数的区间中心值和半径最小化为优化目标，采用遗传算法对翼型进行多目标鲁棒优化设计；(13)判断优化目标是否满足收敛条件，若不满足，转到步骤(10)，更新设计变量，重复步骤(11)～(12)；(14)直至设计目标相邻两次迭代值的变化小于设定容许偏差时，完成二维翼型的多目标鲁棒优化设计；(15)将鲁棒优化得到的翼型与初始翼型及进行传统确定性优化得到的翼型进行对比，分析比较三种翼型的几何形状及气动特性。...

【技术特征摘要】
1.一种基于非概率区间分析模型的翼型鲁棒优化设计方法，其特征在于实现步骤如下：(1)首先，利用类型函数和形状函数转换方法设定参数化翼型曲线的表达式，翼型的几何曲线可用下列函数表示：式中，x/c为翼型弦向的无量纲坐标值，y/c为翼型法向的无量纲坐标值，C(x/c)和S(x/c)为类型函数和形状函数，zte/c为翼型后缘点坐标，C(x/c)可表示为：对于翼型，取指数N1＝0.5,N2＝1，则S(x/c)可表示为：式中，n为多项式阶数，取n＝3；bi为形状函数多项式的加权系数，其中加权系数b0和bn与翼型前缘半径Rle/c及后缘倾角β满足以下关系：c为单位翼型弦长，取c＝1m；在形状函数S(x/c)中取n＝3，提取其中8个参数作为翼型设计变量x，可表示为：x＝(x1,x2,…,x8)＝(Rle/c,β1,β2,zte/c,b1,b2,b1',b2')(5)其中，β1和β2表示上下翼面的后缘倾角，b1、b2、b1'、b2'为上下翼面形状函数多项式的加权系数；(2)设定翼型设计变量的区间上下界，其中下界记为：x＝(x1,x2,…,x8)，上界记为：(3)设定不确定设计参数的区间上下界，记为(4)采用均匀设计方法在由翼型设计变量和不确定参数组成的混合变量空间中生成样本点Ui；(5)对于每一个样本点Ui，根据样本点中的翼型设计变量，建立翼型的CAD模型；(6)由于翼型的CAD模型不能直接用于气动力系数计算，需要将其导入前处理软件GAMBIT中进行流场网格划分；(7)将步骤(6)中创建的网格文件导入FLUENT软件中，来流马赫数根据样本点Ui中的参数...

【专利技术属性】
技术研发人员：邱志平，郑宁宁，王晓军，王睿星，许孟辉，李云龙，王冲，王磊，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京;11