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

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

Chapter 2 Mathematical Modeling for Optimization

2-1 数学模型的建立 (mathematical modeling)

建立数学模型，就是把实际问题按照一定的格式转换成数学表达式的过程。数学模型建立的合适、正确与否，直接影响到优化设计的最终结果。

建立数学模型，通常是根据设计要求，应用相关基础和专业知识，建立若干个相应的数学表达式。如机械结构的优化设计，主要是根据力学、机械设计基础等专业基础知识及机械设备等专业知识来建立数学模型的。

当然，要建立能够反映客观实际的、比较准确的数学模型并非容易之事。数学模型建的过于复杂，涉及的因素太多，数学求解时可能会遇到困难；而建的太简单，又不接近实际情况，解出来也无多大意义。因此，建立数学模型的原则：抓主要矛盾，尽量使问题合理简化。Principle ：The problem is simplified as much as possible.

由于设计对象千变万化，即使对同一个问题，由于看问题的角度不同，数学模型建的可能也不一样。建立数学模型不可能遵循一个不变的规则，本课也不准备把大量的时间花在数学模型的建立上。仅想以几个例子来演示一下数学模型的建立过程，使学生从中得到一些启发。

Exp. 2-1

例2-1 用宽度为cm 24，长度cm 100的薄

铁皮做成cm 100长的梯形槽，确定折边的尺寸

x 和折角θ（如图 2-1所示），使槽的容积最大。

解: 由于槽的长度就是板的长度，槽的梯形

截面积最大就意味着其容积最大。因此，该问题

就由，求体积最大变成求截面积最大。槽的梯形

截面积为： 图 2-1

?=2

1S 高 ?（上底边+下底边） 其中，上底边=x 224-；下底边=θcos 2224x x +-；高=θsin x 定义：该优化设计问题的目标函数是槽的梯形截面积S ，设计变量为θ,x 。问题可以简单地归结为：选择适当的设计变量θ,x ，在一定的限制条件下，使目标函数S 达到最大，限制条件为: 120,20<<<

θ

Exp. 2-2

例2-2 如图 2-2所示是一根简化了的机床主轴。在设计这根轴时，有二个重要因素需要考虑，主轴的重量和外伸端的扰度。对于加工精度要求不高的普通机车而言，以选取主轴重量最轻为优化设计的目标，外伸端的扰度可以作为限制条件来考虑。

图 2-2

解： 当主轴的材料选定后，其重量仅与四个量有关。轴的内经d ，外经D ，支撑间的跨距l 及外伸端a 。由于机床主轴的内孔是用来通过待加工的棒料，其大小由机床型号决定，不能选作设计变量。因此，该问题的设计变量取 a D l ,,；目标函数，即主轴的重量为

))((4122d D a l f -+=πρ ；主轴的限制条件，取它的刚度条件，即外伸端的扰度小于某一规定值 }[y y c ≤及尺寸。 在外力F 作用下，外伸端的扰度为EJ

a l Fa y c 3}(2+= 其中，)(644

4d D J -=π。因此，主轴的刚度约束为][3}(2y EJ a l Fa ≤+。它的尺寸约束为101010,,a a a D D D l l l ≤≤≤≤≤≤。

Exp. 2-3 (p8)

例2-3 如图 2-3所示，钢梁C 的一端与刚性支撑B 焊接在一起,另一端承受作用力6000N 。最优的设计钢梁尺寸，使梁的重量最轻。

图 2-3

解： 钢梁包括梁本身及焊缝，选择独立的设计变量为尺寸,,h l t 和b ，并给定长度 1.4L m =。

用{}{}1234T T

X x x x x h l t b ==表示设计变量。

钢梁的总重量，即目标函数为 c w V V V =+

其中，C V -- 梁C 的体积，立方英寸；W V -- 焊缝的体积，立方英寸。

从图上看，它们的体积分别是

()C V tb L l =+

2212()2W V h l h l ==

所以，总重量为 2()V tb L l h l =++ → 234212()()f X x x L x x x =++

对于焊接钢梁的限制条件有

（1）焊接应力 ()X τ

焊接应力由二部分组成，()'''X τττ=+

，其中，'''MR J ττ= M --F 产生的扭矩，2[(/2)]M F L x =+；J -- 极惯性矩，222311220.707[()122x x x J x x ??+=+????

； 122

2231[]42x x x R ??+=+????

（2）弯曲应力()X σ

最大的弯曲应力为 243

6()FL X x x σ= （3）失稳临界载荷 ()C P X

当 34/t b x x = 值变大，即梁变薄时，会出现失稳的趋势。对于矩形梁，失稳临界载荷

近似地表示成

()C P X 其中，E -- 杨氏模量；334112I x x =

；33413

Gx x α=，G -- 剪切模量 （4）梁的变形()X δ

假定钢梁是长L 的简支梁，其变形是 33344()FL X Ex x δ=

上面四种约束，加上尺寸约束表示如下 1234142536718()()0

()()0

()0

()0

()0

()()0

()0.1250

()0.25()0

d d c g X X g X X g X x x g X x g X x g X P X F g X x g X X ττσσδ=-≥=-≥=-≥=≥=≥=-≥=-≥=-≥

Exp. 2-4

例2-4 某工厂生产B A ,二种产品。产品A 每件需用材料kg 9，3个工时和h kw .4电，产值为60元；产品B 每件需用材料kg 4，10个工时和h kw .5电，产值为120元；若每天可提供材料kg 300，300个工时和h kw .200电，问每天生产B A ,产品各多少件，获得的总产值才能最大？

解： 这是一个生产计划的优化问题。假设每天生产A 产品1x 件，B 产品2x 件，在材料、工时和电力供应量的限制下，求21,x x 的值，使总产值最大。

该优化问题的设计变量为 1x 和2x ；

目标函数为 m ax 1206021→+=x x f

满足限制条件 材料 3604921≤+x x

工时 30010321≤+x x

电力 2005421≤+x x

2-2 数学模型的三要素及一般形式

无论什么样的优化设计问题，尽管其物理概念不同，但数学模型一般均由设计变量、目标函数和约束条件组成，称其为三要素。

2-2-1 设计变量 （Design Variable ）

1) 设计变量

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

在设计中，有些基本参数可以根据设计要求事先给定，称为设计常数，如弹模、许用应力等材料特性等。而有些则需要通过在设计过程中进行调整、优选来定，如尺寸等。对于需要优选的参数，在设计过程中均把它们看作是变化的量，称为设计变量。应注意，设计变量一定是独立参数（Variables must be independent ），任何导出量不能作为设计变量（如式2

1z z i =中只能取三个量中的二个作为设计变量）。设计变量有连续变量和离散变量二种形式（Continuous & Dispersive Variable ）。大多数机械优化设计中的设计变量都是连续变量，可以用常规的优化算法来求解。而对于像齿轮的齿数、模数、钢板的厚度等只能在一定的数集里取值的离散变量的优化设计问题，则需用特殊的优化算法。

2) 设计变量的表示

对于一个优化问题，设计变量的个数则称为该问题的维数（Dimension ），用一数组X 或向量表示：（n-dimensional vectors ）

{}T n n x x x x x x X K M 2121=???

???????????=

以n 个设计变量为坐标轴张成的实空间称为设计空间 （Design Space ），用n R 表示。设计空间中的每一个点都对应着一个设计方案。二个设计变量（2=n ）对应的设计空间是一个平面（plan ），三个设计变量（3=n ）对应的设计空间是一个三维立体空间（space ）（如图 2-4所示），