优化设计的数学模型及基本要素

1、第2章 优化设计的数学模型及基本要素Chapter 2 Mathematical Modeli ng for Optimizati on2-1 数学模型的建立（mathematical modeli ng）建立数学模型，就是把实际问题按照一定的格式转换成数学表达式的过程。数学模型 建立的合适、正确与否，直接影响到优化设计的最终结果。建立数学模型，通常是根据设计要求，应用相关基础和专业知识，建立若干个相应的 数学表达式。如机械结构的优化设计，主要是根据力学、 机械设计基础等专业基础知识及机械设备等专业知识来建立数学模型的。当然，要建立能够反映客观实际的、比较准确的数学模型并非容易之事。数学模型建

2、 的过于复杂，涉及的因素太多，数学求解时可能会遇到困难；而建的太简单，又不接近实际情况，解出来也无多大意义。因此，建立数学模型的原则：抓主要矛盾，尽量使问题合理简化。Principle : The problem is simplified as much as possible.由于设计对象千变万化，即使对同一个问题，由于看问题的角度不同，数学模型建的 可能也不一样。建立数学模型不可能遵循一个不变的规则，本课也不准备把大量的时间花在数学模型的建立上。仅想以几个例子来演示一下数学模型的建立过程，使学生从中得到一些启发。Exp. 2-1例2-1用宽度为24cm，长度100cm的薄 铁皮做成100

3、cm长的梯形槽，确定折边的尺寸 X和折角二（如图2-1所示），使槽的容积最大。解：由于槽的长度就是板的长度， 槽的梯形 截面积最大就意味着其容积最大。 因此，该问题 就由，求体积最大变成求截面积最大。 槽的梯形 截面积为：1S高（上底边+下底边）2其中，上底边=24 -2x ；下底边=24 - 2x - 2xcos ；高=xsin 二定义：该优化设计问题的目标函数是槽的梯形截面积S，设计变量为x门。问题可以简单地归结为：选择适当的设计变量 x门，在一定的限制条件下，使目标函数S达到最大，限制条件为：0, 0 ： x ： 122Exp. 2-2例2-2如图2-2所示是一根简化了的机床主轴。在设计

4、这根轴时，有二个重要因素需要 考虑，主轴的重量和外伸端的扰度。对于加工精度要求不高的普通机车而言，以选取主轴重量最轻为优化设计的目标，外伸端的扰度可以作为限制条件来考虑。图2-2解：当主轴的材料选定后，其重量仅与四个量有关。轴的内经d ,外经D，支撑间的跨距丨及外伸端a。由于机床主轴的内孔是用来通过待加工的棒料，其大小由机床型号决定，不能选作设计变量。因此，该问题的设计变量取l, D, a ；目标函数，即主轴的重量为f = P(l a)(D2 -d2)；主轴的限制条件，取它的刚度条件，即外伸端的扰度小于某Fa(| + a3EJ规定值 yc乞力及尺寸。在外力f作用下，外伸端的扰度为纸其中,J (

5、D4-d4)。因此，主轴的刚度约束为Fa越乞y。它的尺寸约束为643EJl0 < I乞h , D0乞D乞Dj , a0乞a乞印。Exp. 2-3 (p8)例2-3如图2-3所示，钢梁C的一端与刚性支撑 B焊接在一起，另一端承受作用力 6000N。最优的设计钢梁尺寸，使梁的重量最轻。2#图2-3解：钢梁包括梁本身及焊缝，选择独立的设计变量为尺寸h,l,t和b，并给定长度L=1.4m。用X - X X2 X3 x - h l t b A表示设计变量。钢梁的总重量，即 目标函数为v =乂 Vw其中，Vc-梁C的体积，立方英寸；Vw-焊缝的体积，立方英寸。从图上看，它们的体积分别是Vc 二tb（

6、L l）所以，总重量为V =tb(L l) h2l对于焊接钢梁的限制条件有f （X ） =X3X4（L X2） X12X2(1)焊接应力.(X)焊接应力由二部分组成,t(X )=才+才'，其中,F”_MR.2 片x2'JM - F产生的扭矩,M =FL +区/2） ； J -极惯性矩,-I-X；X3 ：:;凶R-X!.24212J(2) 弯曲应力二(X)最大的弯曲应力为；-(X)空X4X3(3) 失稳临界载荷PC(X)当t/b=X3X4值变大，即梁变薄时，会出现失稳的趋势。对于矩形梁，失稳临界载荷近似地表示成Pc（X） =±011其中，E -杨氏模量；I m/x3 ；

7、 ? =Gx3x4， G -剪切模量123(4 )梁的变形、")4FL3假定钢梁是长L的简支梁，其变形是、：(X)二3EX3X4上面四种约束，加上尺寸约束表示如下g1（X）-d -（X）_0 g2（x）h, Y（X） 一0 g3（X）二 x -人 _0 g4（X） =X2 _0 g5（X） = Xa _0 g6（X） =R（X） _F _0 g7（X）订-0.125 _0 g8（X） =0.25 -、（X） _0Exp. 2-4例2-4 某工厂生产代B二种产品。产品 A每件需用材料9kg，3个工时和4kw.h电，产值为60元；产品B每件需用材料4kg，10个工时和5kw.h电，产值为

8、120元；若每天可提 供材料300kg，300个工时和200kw.h电，问每天生产 A, B产品各多少件，获得的总产值 才能最大？解：这是一个生产计划的优化问题。假设每天生产A产品x1件，B产品x2件，在材料、工时和电力供应量的限制下，求x1 ,x2的值，使总产值最大。该优化问题的设计变量为X1和X2 ；目标函数为 f =60xr 120x2r max满足限制条件材料9为 4x2 < 360工时3x110x2 乞 300电力4为 5x2 乞 2002-2数学模型的三要素及一般形式无论什么样的优化设计问题，尽管其物理概念不同，但数学模型一般均由设计变量、 目标函数和约束条件组成，称其为三要

9、素。2-2-1 设计变量（Design Variable ）1）设计变量在机械设计中，一个零件、部件或是一台设备的设计方案，通常是由一组基本参数来 确定和表示的。在设计中，选择哪些参数表示一个设计方案，需要根据各种设计问题的性质来定。有的可以用几何参数，如零件的外形尺寸、截面尺寸、机构的运动学尺寸等；有的可 用某些物理量，如构件的重量、惯性距、频率、力和力矩等；还有的可用一些代表工作性能 的导出量，如应力、扰度、冲击系数等。总之，这些基本参数是对设计指标性能好坏有直接 影响的量。在设计中，有些基本参数可以根据设计要求事先给定，称为设计常数，如弹模、许用 应力等材料特性等。而有些则需要通过在设计

10、过程中进行调整、优选来定，如尺寸等。对于需要优选的参数，在设计过程中均把它们看作是变化的量，称为设计变量。应注意，设计变量一定是独立参数 （Variables must be in depe nde nt），任何导出量不能作为设计变量（如式iL中只能取三个量中的二个作为设计变量）。设计变量有连续变量和离散变量二种形Z2式（Continuous & Dispersive Variable ）。大多数机械优化设计中的设计变量都是连续变 量，可以用常规的优化算法来求解。而对于像齿轮的齿数、模数、钢板的厚度等只能在一定的数集里取值的离散变量的优化设计问题，则需用特殊的优化算法。2）设计变量的表

11、示对于一个优化问题，设计变量的个数则称为该问题的维数（Dimension ），用一数组X或向量表示： （n-dime nsional vectors）L峙X1X = x： = * X2xn 厂xn以n个设计变量为坐标轴张成的实空间称为设计空间（Design Space ），用Rn表示。设计空间中的每一个点都对应着一个设计方案。二个设计变量（n = 2）对应的设计空间是一个平面（plan ），三个设计变量（n = 3）对应的设计空间是一个三维立体空间（space ）（如图2-4所示），Xi r图2-4当维数大于三时(n .3)，设计空间就无法用图来表示，称为超越空间(tran see nd sp

12、ace)。3) 设计变量的选取设计空间的维数表示设计的自由度数。设计变量越多，设计自由度就越大，可供选择 的方案就越多，容易得到比较理性的结果。但随着设计变量数目的增多，必然会使问题复杂化，给寻优带来更大的困难。因此，在满足基本设计要求的前提下，应尽量减少设计变量的个数，把对目标函数影响较大的那些参数选作设计变量。但也应注意实用性，如为了选择一种最合适的材料，将材料的某些性能取为设计变量，但这样求得的最优值，从材料供应方面往往难以实现 (The variables are chose n as a few as possible)。2-2-2 目标函数(Objective Function)1

13、 )目标函数的表示在优化设计中用于评价设计方案好坏的衡量标准(Criteri on )，称为目标函数或评价函数。它是设计变量的函数，记作f(X) = f (x-i x xn)。在工程实际中，优化设计问题的目标函数有二种形式，目标函数的极小化或极大化, 即 (Maximization & Minimization)f(X)； min 或f(X)； m ax其实，目标函数f(X)的极大化就等价于 -f(X)的极小化，为了统一优化算法和程序，以 后最优化均指目标函数的极小化。建立目标函数是整个优化设计中的重要环节。在机械设计中，目标函数主要根据设计 准则来建立的。对于机构的优化设计，这个准则

14、可以是运动学或动力学的特性，如运动误差、振动特性等；对于另部件的设计，这个准则可以是重量、体积、效率等；对于产品设计，也 可以将成本、价格、寿命等作为设计追求的目标。2) 单目标和多目标优化问题(Single- or Multi- Objective Function)在优化设计中，数学模型中仅包含一项设计准则，即目标函数的称为单目标优化问题。同时包含若干个设计准则的就是多目标优化问题。一般来说，目标函数越多，对设计的评价就越周全，设计的综合效果就应该越好， 但对问题的求解就会越复杂。 本课主要解决单目标 优化问题，在最后介绍一些多目标问题的求解方法。3) 目标函数等值线(Level Curv

15、es,Isoline )目标函数f(X)是设计变量x的函数。一组设计变量.x?xnF就代表一个设计方案，在设计空间就确定了一个设计点xk，就有确定的目标函数 f (xk)与之对应。但反过来,定值的目标函数 f(X)=C，却有无穷多个设计点与之对应。这无穷多个目标函数相同的 设计点的集合，就称为目标函数的等值线，二维函数是等值线，三维函数是等值面，三维以 上是等值超曲面。2-2-3 约束条件(Constraint )如前所述，设计空间是所有设计方案的集合，但从工程实用角度上来说，不是所有的设计方案都能接受，如负面积等。为了得到可以接受的(可行的)设计方案(feasible project )，

16、必须根据实际情况和要求，对设计变量的取值加以限制，这种限制就称为约束条件。1)约束种类:等式约束和不等式(Equality & In equality )gu（X）二|凶）乞 0 （u =1,2,川m）hv(X) =0(% X2xj =0 (v =1,2,p n)其中，m, p分别表示不等式和等式约束的数目，注意 p必须小于n,即p ： n。因为，从 理论上讲，一个等式可以消去一个变量，若p = n，则可由n个等式约束中求出唯一的一组设计变量x2xn 丁，就没有优化的余地了。丄边界约束和性能约束 (Boundary& Performanee)边界约束是用于考虑设计变量允许的变化

17、范围的，如一个尺寸应满足Imin乞I乞Imax时可建立不等式约束方程为9 (X ) = 1 min - 1 - 0g2(X) = 1 - Imax - 0性能约束是由某种设计性能或指标推导出来的一种约束条件，如曲柄摇杆机构中曲柄存在的条件等Exp.2-4例2-4如图2-5所示，曲柄摇杆机构，构件的长度分别为，曲柄l1，连杆l2，摇杆l3，机架I 4。设计曲柄摇杆机构必须满足曲柄存在条件，即曲柄为最短杆，最短杆与最长杆之和必须小于其余二杆之和。于是，其性能约束条件为gi（X）=li1 （21,2,3） g4（X） “3 I4 -h I2 g5（x）“2 I4 -h I3 g6（x）J I3 -h

18、 J图2-5应注意在本书中，不等式约束都写成gu (X) _ 0的形式，对于gu (X) < 0的约束，可以写成-gMX) _0的形式。2)可行域与非可行域(Feasible Set )由于引入约束以后，设计点在n维设计空间内被分成二部分。满足约束条件的设计点称为可行设计点，可行设计点的集合称为可行域，位于可行域边界上的设计点亦是可行点(FeasiblePoint )，过该点的约束为起作用约束(Active Constraint )，否则，为不起作用约束(In active Con strai nt);不满足约束条件的设计点的集合为非可行域。下面以二维问题为例说明之。四个不等式约束，一个

19、等式约束，可行、非可行点，起作用、不起作用约束(如图2-6所示)。7#图2-62-2-4数学模型的一般形式(Ge neral Form )由设计变量、目标函数和约束条件组成的数学模型实际上就是优化问题的数学抽象。用文字可以表述为：在满足一定的约束条件下，寻找一组设计变量X =& XJH，使目标函数f(X)达到最优值。其数学表达式为min f (X)二 f 化 x2 xn) X Rns.t.gu(X) =gu(x xj|xn0 (um)hv(X) =1%(为 X2 x=0 (v=1,2, p n)在数学模型中，如果目标函数和约束函数f (X), gu(X), hv(X)都是设计变量 X的

20、线性函数，则称线性规划问题(Li near Programmi ng);否则为非线性规划问题(No nli nearProgramming)。当 m = p = 0时，则称为无约束优化问题 (UneonstrainedOptimization)，当m, p中有一个不为零， 即为约束优化问题(Constrained Optimization) 在工程实际中，不加任何限制的设计问题是很少遇到的，研究它仅有理论意义。而绝大多数的工程优化问题都属于非线性约束优化问题。8#The n -dime nsionalvectors X =xiX2III Xn，whose values are restricted to#satisfy a nu mber equatio nshv(X)=0a nd a set of in equalitiesgu(X)亠0, aresearched to get the minimum of a functionf(X).Review 2Mathematical Modeli ng