《翼型多目标气动优化设计方法.pptx》由会员分享,可在线阅读,更多相关《翼型多目标气动优化设计方法.pptx(23页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、总述总述 导言导言什么是优化什么是优化 优化设计策略优化设计策略优化设计的工作流程优化设计的工作流程翼型表示方法翼型表示方法优化设计算法优化设计算法数值模拟方法数值模拟方法 优化算例优化算例(NACA0012)多目标遗传算法算例多目标遗传算法算例多目标模拟退火算法算例多目标模拟退火算法算例 结论结论第1页/共23页导言导言 什么是优化什么是优化第2页/共23页优化包含的要素 设计变量设计变量 计算过程计算过程 计算结果计算结果 优化目标优化目标 优化算法优化算法第3页/共23页优化设计策略优化设计策略优化设计的工作流程优化设计的工作流程翼型表示方法翼型表示方法优化设计算法优化设计算法数值模拟方
2、法数值模拟方法第4页/共23页优优化化算算法法2、计算翼型曲线3、gambit生成网格4、Fluent计算6、是否满足终止条件1、生成初始参数程序终止7、由 modelFRONTIER的优化算法生成一组新的参数5、得到Cl,Cd,Cl/Cd否是设计变量设计变量计算过程计算过程结果结果目标目标第5页/共23页翼型表示方法翼型表示方法 解析函数形状扰动方法解析函数形状扰动方法 为为Hicks-HenneHicks-Henne函数簇函数簇:m、n为参数为参数i i即为设计变量,改变即为设计变量,改变i i的值即可得到的值即可得到一系列光滑翼型一系列光滑翼型第6页/共23页Hicks-HenneHic
3、ks-Henne函数簇函数簇:函数满足的条件函数满足的条件 参数参数m m对应函数极值点的位置,函数在对应函数极值点的位置,函数在x=mx=m点处取到最大值点处取到最大值1 1,向两侧迅,向两侧迅速减小;速减小;函数在函数在0 0、1 1点对应取点对应取0 0值,保证了上下一面结合处的连续性;值,保证了上下一面结合处的连续性;参数参数n n对应函数的形状,对应函数的形状,n n的值越大,函数峰值两侧下降的速度越快;的值越大,函数峰值两侧下降的速度越快;函数在函数在x=0 x=0点导数为零,这就保证了上下翼面在点导数为零,这就保证了上下翼面在0 0点结合处的光滑性。点结合处的光滑性。第7页/共2
4、3页Hicks-HenneHicks-Henne函数簇函数簇:函数图形函数图形 m mi i=0.1=0.1、0.30.3、0.50.5、0.70.7、0.90.9 n ni i=3 =3 、3 3、3 3、3 3、1 (i=11 (i=1、2 2、3 3、4 4、5)5)第8页/共23页优化设计算法优化设计算法:遗传算法遗传算法 简简单单遗遗传传算算法法的的遗遗传传操操作作主主要要有有三三种种:选选择择(selection)selection)、杂杂交交(cross-over)cross-over)、变变异异(mutation)mutation)。选选择择操操作作根根据据个个体体的的适适应应
5、度度函函数数值值所所度度量量的的优优劣劣程度决定它在下一代是被淘汰还是被遗传。程度决定它在下一代是被淘汰还是被遗传。杂交算法杂交算法交换随机挑选的两个个体的某些位,交换随机挑选的两个个体的某些位,变变异异算算子子则则直直接接对对一一个个个个体体中中的的随随机机挑挑选选的的某某一一位进行突变。位进行突变。第9页/共23页优化设计算法:优化设计算法:模拟退火算法模拟退火算法 模拟退火算法是一种基于模拟退火算法是一种基于Monte Carlo Monte Carlo 迭代求解的迭代求解的启发式随机搜索算法启发式随机搜索算法,它源于,它源于固体退火原理,首先将固体加温至充分高,固体退火原理,首先将固体
6、加温至充分高,使固体内部的粒子变为无序状,内能增大,使固体内部的粒子变为无序状,内能增大,再让其徐徐冷却,粒子渐趋有序,最后在再让其徐徐冷却,粒子渐趋有序,最后在常温时达到基态,内能减为最小。模拟退常温时达到基态,内能减为最小。模拟退火算法用内能火算法用内能E E模拟模拟目标函数目标函数,用温度,用温度T T作作为为控制参数控制参数。第10页/共23页优化设计算法优化设计算法:多目标算法及约束条件处理多目标算法及约束条件处理 ParetoPareto边界边界的处理方法:仅对的处理方法:仅对ParetoPareto边界进行操作,这就同时考虑到了边界进行操作,这就同时考虑到了每个目标的进展,并且能
7、大大简化运算。每个目标的进展,并且能大大简化运算。权系数权系数w wi i处理为函数:处理为函数:用一定的算法计算当前用一定的算法计算当前w wi i的值以消除操作中的的值以消除操作中的主观因素,使得对每个目标的处理更有效。主观因素,使得对每个目标的处理更有效。第11页/共23页 利用利用GambitGambit生成生成结构化网格能大大提结构化网格能大大提高高FLUENTFLUENT的计算速度的计算速度和精度。和精度。数值模拟方法数值模拟方法:网格生成网格生成第12页/共23页数值模拟方法数值模拟方法:流场计算流场计算 本文选择本文选择N-S方程方程作为流场作为流场计算的计算的主控方程主控方程
8、,湍流模式,湍流模式使用使用Spalart-Allmaras(S-A)模型)模型,边界采用,边界采用远场压远场压力条件力条件采用采用一阶迎风差分格一阶迎风差分格式式离散控制方程。离散控制方程。我们首先验证了翼型在我们首先验证了翼型在不同攻角不同攻角下的下的气动性能气动性能,其,其中,中,2o 攻角,攻角,2.63马赫下,马赫下,计算值与实验值的压力系数计算值与实验值的压力系数曲线对比如图。证明了流场曲线对比如图。证明了流场计算模型和参数的选择是相计算模型和参数的选择是相当有效的。当有效的。第13页/共23页优化算例及结果优化算例及结果 多目标遗传算法(多目标遗传算法(MOGA)算例算例 多目标
9、模拟退火算法多目标模拟退火算法(MOSA)算例算例 结果讨论结果讨论第14页/共23页 优化算例的设计要求优化算例的设计要求 初始翼型选用NACA0012,Ma=0.75,迎角2.57o;优化后翼型的最大厚度不小于原始翼型的厚度;要求提高升力、降低阻力、提高升阻比。第15页/共23页优化算例优化算例:翼型及压力系数曲线对比翼型及压力系数曲线对比优化前后机翼形状及压力系数对比曲线优化前后机翼形状及压力系数对比曲线第16页/共23页MOGA算例:算例:结果表格结果表格 ClCdCl/Cd初始翼型(2.57o)0.43630.024218.03中间翼型(2.57o)0.57800.017832.47
10、优化翼型(2.57o)0.68120.013450.84优化翼型(1.829 o)0.52370.011545.54表1、多目标遗传算法优化结果表格第17页/共23页MOSA算例:算例:结果表格结果表格 ClCdCl/Cd初始翼型(2.57o)0.43630.024218.03中间翼型(2.57o)0.49680.016130.86优化翼型(2.57o)0.52060.010947.76优化翼型(2.5969 o)0.52730.0109648.11表2、多目标模拟退火算法优化结果第18页/共23页MOGA算例:算例:图形对比图形对比MOGA优化前后流场等马赫数曲线及翼面附近速度矢量图 第19页/共23页MOSA算例:算例:图形对比图形对比MOSA优化前后流场等马赫数曲线图 第20页/共23页 MOGA优化后翼型升阻比提高182.0%,而其中形状变化带来的提高是152.6%,其余部分则源自有效攻角的增大。MOSA优化后翼型升阻比提高164.9%,但单纯形状变化带来的提高是166.8%,翼型形状优化效果更加明显。结果比较结果比较第21页/共23页Questions?Thank you第22页/共23页感谢您的观看!第23页/共23页