工程优化报告

华中科技大学1/15工程优化方法综述一、优化建模1概述建立数学模型是优化设计中的一个重要的组成部分。优化结构是否可用，主要取决于所建立的数学模型是否能够确切而简介地反映工程问题的客观实际。在建立数学模型时，片面地强调确切，往往会使数学模型十分冗长、复杂，增加求解问题的困难程度，有时甚至会使问题无法求解；片面强调简洁，则可能使数学模型过分失真，以致失去了求解的意义。合理的做法是在能够确切反映工程实际问题的基础上力求简洁。设计变量、目标函数和约束条件是组成优化设计模型的三要素，下面分别予以讨论。2设计变量的选择机械设计中的所以参数都是可变的，但是将所以的设计参数都列为设计变量不仅会使问题复杂化，而且是没有必要的。例如材料的机械性能由材料的种类决定，在机械设计中常用材料的种类有限，通常可根据需要和经验事先选定，因此诸如弹性模量、泊松比、许用应力等参数按选定材料赋以常量更为合理；另一类状态参数，如功率、温度、应力、应变、绕度、压力、速度、加速度等则通常可由设计对象的尺寸、载荷以及各构件间的运动关系等计算得出，多数情况下也没有必要作为设计变量。因此，在充分了解设计要求的基础上，应根据各设计参数对目标函数的影响程度认真分析其主次，尽量减少设计变量的数目，以简化优化设计问题。另外还应注意设计变量应当相当独立，否则会使目标函数出现“山脊”或“沟谷”，给优化带来困难。3.目标函数的确定目标函数是一项设计所追求的指标的数学反映，因此对它最基本的要求是能够用来评价设计的优劣，同时必须是设计变量的可计算函数。选择目标函数是整个优化设计过程中最重要的决策之一。有些问题存在明显的目标函数，例如一个没有特殊要求的承受静载的梁，自然希望它越轻越好，因此选择其自重作为目标函数是没有异议的。但涉及一台复杂的机器，追求的目标往往较多，华中科技大学2/15就目前使用较成熟的优化方法来说，还不能把所有要追求的指标都列为目标函数，因为这样做并不一定能有效地求解。因此应当对所追求的各项指标进行细致的分析，从中选择最重要最具有代表性的指标作为设计追求的目标。例如一架好的飞机，应该具有自重轻、静载重量大、航程长、使用经济、价格便宜、跑道长度合理等性能，显然这些都是设计时追求的指标。但并不需要把它们都列为目标函数，在这些指标中最重要的指标是飞机的自重。因为采用轻的零部件建造的自身重量最轻的飞机只会促进其它几项指标，而不会损害其中任何一项。因此选择飞机自重作为优化的目标函数应该是最合适了。若一项工程设计中追求的目标是互相矛盾的，这时常常取其中最主要的指标作为目标函数，而其余的指标列为约束条件。也就是说，不指望这些次要参数都达到最优，只要它们不至于过劣就可以了。4约束条件的确定约束条件是就工程设计本身而提出的对设计变量取值范围的限制条件。和目标函数一样，它们也是设计变量的可计算函数。如前所述，约束条件可分为性能约束和边界约束两大类。性能约束通常与设计原理有关，有时非常简单，如设计曲柄连杆机构时，按曲柄存在条件而写出的约束函数均为设计变量的线性显函数；有时却相当复杂，如对一个复杂的机构系统，要计算其中各构件的应力和位移，常采用有限元法，这时相应的约束函数为设计变量的隐函数，计算这样的约束函数往往要花费很大的计算量。在选择约束条件时应当特别注意避免出现相互矛盾的约束。因为相互矛盾的约束必然导致可行域为一空集，使问题的解不存在。另外应当尽量减少不必要的约束，不必要的约束不仅增加优化设计的计算量，而且可能使可行域缩小，影响优化结果。二、优化理论1综述优化理论归根到底是优化数学理论，机械优化设计是建立在多元函数的极值理论上的。主要包括以下几个方面的内容：多元函数的方向导数与梯度多元函数的泰勒展开华中科技大学3/15无约束优化问题的极值条件凸集、凸函数与凸规划等式约束优化问题的极值条件不等式约束优化问题的极值条件2K-T条件在以上众多理论中，K-T条件应用最为广泛，下面重点讨论一下K-T条件。不等式约束多元函数极值条件：用起作用约束的下标集合表示用梯度形式表示，可得K-T条件的几何意义：在约束极小点处，函数的负梯度一定能表示成所有起作用约束在该点梯度的非负线性组合。三、优化算法1.无约束优化算法1.1概述无约束优化方法是优化技术中基本的也是非常重要的内容。无约束优化问题的数学模型如下：**101,2,...,01,2,...,01,2,...,mjjjiijjjdfxdgxindxdxgxjmjm***01,2,...,00jjjJiijjdfxdgxindxdxgxjJjJ**jjjJfxgx12min()[]nnFxxxxxR　　华中科技大学4/15使用无约束优化方法，不仅可以直接求无约束优化设计问题的最优解，而且通过对无约束优化方法的研究给约束优化方法建立明确的概念、提供良好的基础，某些优化设计方法就是先把约束优化设计问题转化为无约束问题后，再直接用无约束优化方法求解。无约束优化问题一般模型为给定x(0)，置k=1a)确定局部下降方向d(k),使▽f(x(k))Td(k)0；b)确定步长tk0,使f(x(k)+tkd(k))f(x(k))；c)令x(k+1)=x(k)+tkd(k)；d)判断是否终止，终止则输出，否则k:=k+1，返回(a)从上面的分析我们看到一个算法的搜索方向成为该优化方法的基本标志，分析、确定搜索方向d(k)是研究优化方法的愚根本的任务之一。根据构成搜索方向使用的信息性质不同，无约束优化可以分为两类。一类是利用目标函数的一阶导数或二阶导数的无约束优化方法，如最速下降法、共轭梯度法、牛顿法等。另一类是只利用目标函数值得无约束优化算法，如坐标轮换法、单纯形法和鲍威尔法。下面将分别讨论。1.2最速下降法优化设计是追求目标函数值最小，因此，自然可以设想从某点出发，其搜索方向取该点的负梯度方向，使函数值在该点附近下降最快。这种方法也称为最速下降法。1.2.1基本原理梯度法的迭代公式为：(k1)(k)(k)(k(k)(k))(111)()()()()()()(      `0 ())()))0 0(kkkkkkkkkkkkkTkkTkkTgFxxfxfxminaminfxagfxagafxaggfxgxgxgg其中是函数在迭代点处的梯度一般采用一维搜索的最优步长，即据一元函数极值条件和多元复合函数求导公式，得即（或华中科技大学5/15此式表明，相邻的两个迭代点的梯度是彼此正交的。也即在梯度的迭代过程中，相邻的搜索方向相互垂直。梯度法向极小点的逼近路径是锯齿形路线，越接近极小点，锯齿越细，前进速度越慢。这是因为，梯度是函数的局部性质，从局部上看，在该点附近函数的下降最快，但从总体上看则走了许多弯路，因此函数值的下降并不快。1.2.2算法流程1.3共轭梯度法共轭梯度法是共轭方向法的一种，因为该方法中每一个共轭向量都是依赖于迭代点处的负梯度而构造出来的，所以称作共轭梯度法。1.3.1基本原理共轭梯度法的搜索方向采用梯度法基础上的共轭方向，如图所示，目标函数F(x)在迭代点x(k+1)处的负梯度为(k1)(x)f，该方向与前一搜索方向(k)互为正交，在此基础上构造一种具有较高收敛速度的算法，该算法的搜索方向要满足以下两个条件：华中科技大学6/15以x(k+1)点出发的搜索方向(k1)是(k1)(x)f与(k)的线性组合。即(k1)(k1)(k)f(x)k以与为基底的子空间中，矢量与相共轭，即满足(k1)T(k)[]0G1.3.2算法流程1.4DFP变尺度法变尺度法也称拟牛顿法，它是基于牛顿法的思想而又作了重大改进的一类方法。我们所介绍的变尺度法是由Davidon于1959年提出又经Fletcher和Powell加以发展和完善的一种变尺度法，故称为DFP变尺度法。1.4.1基本思想变尺度法的基本思想与牛顿法和梯度法有密切联系。观察梯度法华中科技大学7/15和牛顿法的迭代公式和(k1)(k)(k)(k)(k1)(k)(k)1(k)kxgxHgxx分析比较这两种方法可知：梯度法的搜索方向为-g(k)，只需计算函数的一阶偏导数，计算量小，当迭代点远离最优点时，函数值下降很快，但当迭代点接近最优点时收敛速度极慢。牛顿法的搜索方向为-Hk-1g(k)，不仅需要计算一阶偏导数，而且要计算二阶偏导数及其逆阵，计算量很大，但牛顿法具有二次收敛性，当迭代点接近最优点时，收敛速度很快。若迭代过程先用梯度法，后用牛顿法并避开牛顿法的海赛矩阵的逆矩阵的烦琐计算，则可以得到一种较好的优化方法，这就是“变尺度法”产生的基本构想。为此，综合梯度法和牛顿法的优点，提出变尺度法的基本思想。变尺度法的基本迭代公式写为下面的形式：1(k1)(k)(k)kkx     Agx kkkkkkkkAnnAIAHgxAg式中的为构造的阶对称矩阵，它是随迭代点的位置的变化而变化的。若，上式为梯度法的迭代公式，若，上式为阻尼牛顿法的迭代公式。变尺度法的搜索方向称为拟牛顿方向。1.4.2算法流程华中科技大学8/151.5小结下面讨论本章所介绍的几种优化方法的适用范围：可靠性：牛顿法较差，因为它对目标函数要求太高，解题成功率较低。有效性：坐标变换法和梯度法的计算效率较低，因为它们从理论上不具有二次收敛性。简便性：牛顿法和变尺度法的程序编制较复杂，牛顿法还占用较多的存储单元。由以上讨论可知，在选用无约束优化方法时，一方面要考虑优化方法的特点，另一方面要考虑目标函数的情况。1、一般而言，对于维数较低或者很难求得导数的目标函数，使用坐标轮换法或鲍威尔法较合适。2、对于二次性较强的目标函数，使用牛顿法效果好。3、对于一阶偏导数易求的目标函数，使用梯度法可使程序编制简单，但精度不宜过高。4、综合而言，鲍威尔法和DFP法具有较好的性能。2.约束优化算法华中科技大学9/152.1惩罚函数法2.1.1基本思想将有约束条件的优化问题转化为无约束优化问题来求解。前提：一是不能破坏约束问题的约束条件，二是使其最优解归结到原约束问题的同一最优解上去。nkmin(x),xR.(x)0j1,2,...mh(x)0k1,2,...ljfstg构成一个新的目标函数，称为惩罚函数根据惩罚项得不同构成形式，惩罚函数法又可分为外点惩罚函数法、内点惩罚函数法和混合惩罚函数法。2.1.2算法流程2.2拉格朗日乘子法2.2.1PH算法效用函数为华中科技大学10/15()()min(,,)()()()()kkTTkkMxufxuhxhxhx判断函数为()kkhx当()()kkx时迭代停止。步骤1：选定初始点(0)x，初始拉格朗日乘子向量(1)u，初始罚因子1及其放大系数1c，控制误差0与常数(0,1)，令1k。步骤2：以(1)kx为初始点，求解无约束问题：()()min(,,)()()()()kkTTkkMxufxuhxhxhx得到无约束问题最优解()kx步骤3：当()khx时，()kx为所求的最优解，停；否则转步骤4.步骤4：当()()/kkhxhx时，转步骤5；否则令k1kc，转步骤5.步骤5：令(1)()()(),1kkkkuuhxkk，转步骤1。2.2.2PHR算法松弛变量法：12222111(,,)()max0,()2()()2iiilljjjjjMuvfxugxuvhxhx乘子的修正公式为：(1)()()(1)()()(),1,...,max0,(),1,...,kkkjjjkkkiiivvhxjluugxim华中科技大学11/15判断函数为：1/