约束优化设计.docx

第四章 约束优化设计l 概述l 约束坐标轮换法l 随机方向法l 罚函数法概述结构优化设计的问题,大多属于约束优化设计问题,其数学模型为:.min f (x)gu (x) £ 0hv (x) = 0x Î Rnu = 1, 2,..., mv = 1, 2,..., p < n行域 φ 内,选择一个初始点 X   然后确定一个可行得一个目标函数有所改善的可行的新点 X  即完成了根据求解方式的不同,可分为直接解法和间接解法两类。直接解法是在仅满足不等式约束的可行设计区域内直接求出问题的约束最优解。属于这类方法的有:随机实验法、随机方向搜索法、复合形法、可行方向法等。其基本思路:在由 m 个不等式约束条件 gu(x)≤0 所确定的可0搜索方向 S,且以适当的步长沿 S 方向进行搜索,取1一次迭代。以新点为起始点重复上述搜索过程,每次均按如下的基本迭代格式进行计算:X k+1=X k +a kSk(k=0,1,2,..) 逐步趋向最优解,直到满足终止准则才停止迭代。直接解法的原理简单,方法实用,其特点是:1) 由于整个过程在可行域内进行,因此,迭代计算不论何时终止,都可以获得比初始点好的设计点。2) 若目标函数为凸函数,可行域为凸集,则可获得全域最优解,否则,可能存在多个局部最优解,当选择的初始点不同,而搜索到不同的局部最优解。3) 要求可行域有界的非空集a) 可行域是凸集;b)可行域是非凸集间接解法间接解法是将约束优化问题转化为一系列无约束优化问题来解的一种方法。由于间接解法可以选用已研究比较成熟的无约束优化方法,并且容易处理同时具有不等式约束和等式约束的问题。因而在机械优化设计得到广泛的应用。间接解法中具有代表性的是惩罚函数法。将约束函数进行特殊的加权处理后,和目标函数结合起来,构成一个新的目标函数,即将原约束优化问题转化为一个或一系列的无约束优f  (X ,m 1,m 2 ) = F (X )+ åm 1G éë g j (X )ùû + åm 2H éëhk (X )ùû化问题。m                      lj=1 k=1新目标函数然后对新目标函数进行无约束极小化计算。加权因子间接法是结构优化设计中广泛使用的有效方法,其特点:1) 由于无约束优化方法的研究日趋成熟,为间接法提供可靠基础。这类算法的计算效率和数值计算的稳定性大有提高;2) 可以有效处理具有等式约束的约束优化问题;3) 目前存在的主要问题,选取加权因子较为困难,选取不当,不仅影响收敛速度和计算精度,甚至导致计算失败。在可行域任取一点 ,取一个初始步长 ,X 0a按 ,取得沿 坐标轴第一个1 1X X ea= + 1x如图所示,直到逼近最优点 。X从初始点 出发,沿 方向以一定X Sl 约束坐标轮换法约束坐标轮换法是在无约束坐标轮换法的基础上,再加上由约束条件构成的可行性逻辑判断而构成的方法,这样可以使搜索点保持在可行域内,求得最优解。迭代步长不是采用最优步长,而是加速步长。其基本思路:01 0迭代点,检查该点是否满足可行性和适用性:X1 Î D (可行性条件)F ( X1 ) < F ( X 0 ) (适用性条件)若两者均满足,步长加倍,迭代计算1X 2 = X 0 + 2a e1 ,只要迭代点满足条件,加倍增大步长,继续迭代获得新点;当迭代点不满足条件,取前一个迭代点,转而沿 x2 坐标轴方法搜索,不满足条件时,取负步长进行,*约束坐标轮换法虽然方法简单、算法明确,便于设计,但当维数较高时收敛速度慢,还会出现“死点”,导致出现伪最优点。l 约束随机方法随机方向法的基本思路:在可行域内选择一个初始点,利用随机数的概率特性,产生若干个随机方向,并从中选择一个能使目标函数值下降最快的随机方向作为搜索方向S。0步长进行搜索,得到新点 X ,新点X 应满足约束条件且f ( X ) < f ( X 0 ) ,至此完成一次迭代。随机方向法程序设计简单,搜索速度快,是解决小型机械优化问题的十分有效的算法。如图所示。1. 随机数的产生首先令 r1 = 2  , r2 = 2  , r3 = 2   取 r=2657863,按一下步骤计算:下面介绍一种常用的产生随机数的数学模型35 36 37令r ¬ 5rr ¬ r - r2r ³ r r ¬ r - r则   q = r / r1随机方向法的初始点 X   必须是一个可行点,既满足全部不等式约束条件。若 r ³ r3 则 r ¬ r - r3若 r ³ r2 则若 则(0,1)之间的随机数在任意(a,b)区间内的随机数x = a + q(b - a)2. 初始点的选择0初始点可以通过随机选择的方法产生。1)输入设计变量的下限值和上限值,即ai £ xi £ bi2)在区间(0,1)内产生 n