基于遗传算法的平面叶栅多目标优化设计 免费.docVIP

基于遗传算法的平面叶栅多目标优化设计 1??数学模型????奇点法[2]的基本出发点是用一系列分布在翼型骨线上的奇点来代替叶栅中的翼型对水流的作用，将叶栅绕流的计算转化为基本势流的叠加计算，利用绕流无分离的条件来绘制翼型的形状。其前提是假定来流为无旋有势流动、叶片无限薄。在设计过程中，所求的骨 线可先假设一个翼型的骨线形状，计算出骨线上各点的合成速度W，由于骨线 是假定的，W并不能和骨线相切。根据骨线和速度W相切的条件修改第一次假设的骨 线形状，得到第二次近似骨线。重复上述计算，直至逼近为止。1.1??目标函数??优化模型为： (1) 式中：ξ为损失系数，σ为气蚀系数；βi(i=1,2,3,…6)为奇点法设计骨线时得到的6个等分计算点的夹角(如图1)，奇点法就是通过这六个夹角和骨线总长来得到骨线形状。????得到骨线后，为得到翼型，必须以骨线为中线进行加厚，继而得到平面叶栅。本文先根据强度要求来确定翼型的最大厚度，而后以NACA-0012型空气动力翼型的厚度分布规律为基础，对所得骨线进行双边加厚。1.2??约束条件??叶栅设计时的约束条件有几何约束、流动约束、能量约束等。在确定约束条件时应考虑在允许的情况下，尽量减少约束条件的数目。本课题提出约束条件有环量约束，以此来满足给定的环量，除此之外，还有速度分布约束、边界层无分离约束等。1.3??多目标优化设计方法??设优化的多目标优化模型[3]为：  (2) 式中：X为优化个体；x1,x2,…,xk为优化变量；f1(X),f2(X)，…，fn(X)为各个优化目标。????在多目标优化过程中，通常引入一个“有效解”，即Pareto最优解的概念。本文采用以下形式构造统一的目标函数： (3) 这里，P为正整数，通常取P=2;X*i为第i个目标对应的最优个体，fi(X*i)为其最优解。如此构造的统一目标函数可同时可虑多个优化目标。1.4??个体适应值的计算??要得到个体适应值，需要对已知翼型的叶栅进行流场分析。取与水轮机同时转动的相对坐标系，则所求解问题的相对流函数满足拉普拉斯方程。鉴于我们关心的是翼型表面的速度分布这一特点，同时由于拉普拉斯方程的基本解已知，参照边界元[4]的特点，采用边界元方法来计算由翼型组成叶栅的平面无旋流场。流场求出之后，即可得到翼型表面的流速分布，从而可求得翼型的气蚀系数；而后，假设叶栅中的总损失正比于叶栅表面边界层中的流动损失，采用积分法[5]计算平面叶栅边界层，从而得到损失系数。在气蚀系数和损失系数已知的条件下，即可求得个体的适应值。 2? 数值算例2.1??数值算例的给定设计参数??本文以ZZ440轴流式水轮机为研究对象，对轮毂处的平面叶栅进行优化设计。基本参数[6]为：Q10=0.8m3/s，n10=115r/min,ns=440mkw，Z=6,Dh=0.5给定设计参数为：nd=1.3n10,Qd=1.45Q10,H=1m,D=1m。流场计算区域如图2所示。2.2??优化变量取值范围??在设计出初始无厚翼型骨线时，即可得到的骨线上六个骨线等分计算点的夹角值。各优化变量βi在初始骨线角的基础上左右波动20%，即可得到优化变量的取值范围。如在第i个优化点的初始骨线角度为β0i，则在该点的优化变量βi的变化范围[0.8β0i，1.2β0i]。如此确定的设计变量取值的上下限可在编码方案中自动被考虑，无须再作为优化约束另行处理。 2.3??遗传算法参数的选择??在本例中，采用可以较好处理多峰问题的基于共享机制的小生境技术遗传算法[7,8]。对各个优化参量采用实数编码。权衡考察结果，当种群规模N=40，杂交概率Pc=0.7，变异概率Pm=0.3，海明距离阀值Ls=0.5时，就本例而言，可望得到较好的算法形态。2.4??数值计算结果及其分析??上述思想用Fortran PowerStation语言编程实现，先分别以叶栅损失和气蚀系数最小为单目标进行优化，得出单目标优化结果：σmin=0.36，ξmin=0.050。而后，运用距离法进行多目标优化。遗传算法运行至第32代时达到最优，其优化设计的计算结果如图所示。 2