[理学]第5章 约束优化方法2010.ppt

例1 用内点法求 的约束最优解。 解: 用内点法求解该问题时，首先构造内点惩罚函数: 用解析法求函数的极小值，运用极值条件： 联立求解得： 时不满足约束条件 应舍去 。 无约束极值点为 当 ——惩罚函数法 1）? 初始点x0的选取 使用内点法时，初始点应选择一个离约束边界较远的可行点。如太靠近某一约束边界，构造的惩罚函数可能由于障碍项的值很大而变得畸形，使求解无约束优化问题发生困难。 2）? 惩罚因子初值r0的选取 惩罚因子的初值应适当，否则会影响迭代计算的正常进行。一般而言，太大，将增加迭代次数；太小，会使惩罚函数的性态变坏，甚至难以收敛到极值点。无一般性的有效方法。对于不同的问题，都要经过多次试算，才能决定一个适当 r0。 2 、参数选择及收敛条件 一般的看法是，c值的大小在迭代过程中不起决定性作用，通常的取值范围在0.1~0.7之间。 4) 收敛条件 ? 3) 惩罚因子的缩减系数c的选取 惩罚因子r是一个逐次递减到0的数列，相邻两次迭代的惩罚因子的关系为 : 3、算法步骤： 1）选择可行域内初始点X(0); 2）选取初始罚因子r(0)与罚因子缩减系数c，并置K←0； 3）求minφ(x(K),r(K))解出最优点xK*； 4）按终止准则判别，若满足转步骤5），否则令K←K+1，r(K+1)←r(K)， xK+10←xK*转步骤3）； 5）输出最优解（X*，F*），停止计算。 内点法程序流程图 = = = + = = 优化设计实例: 欲设计一空心传动轴, D和d分别为空心轴的外径和内径, 轴长L=2m。轴的材料密度?=7.8×103kg/m3，杨氏模量E=2 × 105MPa，剪切模量G=80GPa，许用剪应力[?]=50MPa, 单位长度许用扭转角[?]=0.5?/m , 轴所传递的功率P=10KW，转速n= 400r/min。要求在满足使用性能条件和结构尺寸限制的前提下使其重量最小。 1.设计变量和目标函数 该传动轴的力学模型是一个受扭转的圆柱轴。其外径和内径是决定圆轴的主要独立参数, 故将其作为设计变量。写成向量形式为X=[x1 x2]T=[D d]T。依题意, 取重量最小为优化目标。空心圆轴的重量可按下式计算： w = ??gL(D2-d2)/4 2.约束条件 (1)扭转强度 根据扭转强度, 要求扭转剪应力需满足： (2)扭转刚度 单位长度的最大扭转角不超过规定的许用值: (3)扭转稳定性 扭转剪应力不得超过临界剪应力： (4)结构尺寸 由结构尺寸要求决定的约束条件为： 3.数学模型 将所有函数表达式规范化并代人已知数据, 可得传动轴优化设计的数学模型为： 优化结果: D = 43.2mm d = 8.7mm wmin = 215.25N 二、外点惩罚函数法 外点法是从可行域的外部构造一个点序列去逼近原约束问题的最优解。外点法可以用来求解含不等式和等式约束的优化问题。 外点惩罚函数的形式为： r是惩罚因子 , 外点法的迭代过程在可行域之外进行，惩罚项的作用是迫使迭代点逼近约束边界或等式约束曲面。由惩罚项的形式可知，当迭代点x 不可行时，惩罚项的值大于0。 ——惩罚函数法 用外点法求下列问题的约束最优解 ——惩罚函数法 例2 用外点法求解下列有约束优化问题 解：惩罚函数为： 对上式求偏导，得 无约束目标函数极小化问题的最优解系列为： 当惩罚因子渐增时，由下表可看出收敛情况。 r 0.01 -0.80975 -50.00000 -24.9650 -49.9977 0.1 -0.45969 -5.00000 -2.2344 -4.9474 1 0.23607 -0.50000 0.9631 0.1295 10 0.83216 -0.05000 2.3068 2.0001 1000 0.99800 -0.00050 2.6624 2.6582 ∞ 1 0 8/3 8/3 内点法和外点法的简单比较 内点法的特点： ??? （1）始点必须为严格内点 ??? （2）不适于具有等式约束的数学模型 ??? （3）迭代过程中各个点均为可行设计方案 ??? （4）一般收敛较慢 ??? （5）初始罚因子要选择得当 ??? （6）罚因子为递减，递减率c有0c1 ? ?????外点法的特点： ??? （1）初始点可以任选 ??? （2）对等式约束和不等式约束均可适用 ??? （3）仅最优解为可行设计方案 ??? （4）一般收敛较快 ??? （5）初始罚因子要选择得当 ??? （6）罚因子为递增，递增率c’有c’1 三、混合惩罚函数法 1、混合惩罚函数法