MATLAB第12章工程优化设计实例
MATLAB第12章工程优化设计实例
第12章工程优化设讣实例优化设讣的数学模型
优化设讣的数学模型是描述实际优化问题的设计内容、变量关系、有关设计条件和意图的数学表达式，它反映了物理现象各主要因素的内在联系，是进行优化设计的基础。优化设计数学模型的三大要素：
?设计变量
?约束条件
? LI标函数
1.设计变量
一个设计方案可以用一组基本参数的数值来表示，这些基本参数可以是构件尺寸等儿何量，也可以是质量等物理量，还可以是应力、变形等表示丄作性能的导出量。设计变量:在设讣过程中进行选择并最终必须确定的各项独立的基本参数，乂叫做优化参数。
设计变量的全体实际上是一组变量，可用一个列向量表示
2.约束条件
设计空间是所有设计方案的集合,但这些设计方案有些是工程上所不能接受
的。如一个设计满足所有对它提出的要求，就称为可行设计。
一个可行设讣必须满足某些设讣限制条件，这些限制条件称作约束条件,简称约
束。3.目标函数
为了对设计进行定量评价，必须构造包含设讣变量的评价函数，它是优化的口标, 称为LI标函数，以F(X)表示。
FxFxxx ()(),, , , 12n
在优化过程中，通过设计变量的不断向F(X)值改善的方向自动调整，最后求得 F(X)值最好或最满意的X值。在构造LI标函数时,应注意：訂标函数必须包含全部设计变量，所有的设计变量必须包含在约束函数中。
模型输入时需要注意的问题
使用优化工具箱时，山于优化函数要求口标函数和约束条件满足一定的格式，所以需要用户在进行模型输入时注意以下儿个问题：
1.目标函数最小化
优化函数 fminbnd、fminsearchx fminunc、fmincon、fgoalattdin、fminmax 和1 sqnonlin都要求LI标函数最小化，如果优化问题要求LI标函数最大化，可以通过使该LI标函数的负值最小化即-f(x)最小化来实现。近似地，对于quadprog函数提供-H 和-f,对于1 inprog函数提供-f。
2.约束非正
优化工具箱要求非线性不等式约束的形式为Ci(x)?0,通过对不等式取负可以达到使大于零的约束形式变为小于零的不等式约束形式的U的，如Ci (x)?0形式的约束等价于- Ci(x)?0;Ci(x)?b形式的约束等价于- Ci(x)+b?0。
3.避免使用全局变量
Fmincon是mat lab最主要内置的求解约束最优化的函数，该函数的优化问题的标准形式为：1.数学模型标准形式：
min f , X,
s. t. AX?b ,线性不等式约束，
AeqX二beq ,线性等式约束，
C(X)?0 ,非线性不等式约束条件，
Ceq(X)二0,非线性等式约束,
Lb ?X ?Ub ,边界约束条件，
Fmincon函数使用的约束优化算法都是訂前比较普适的有效算法：
,1,对于中等的约束优化问题fmincon使用序列?.次规划,SQP sequential quadratic programming,算法；
,2,对于大规模约束优化问题fmincon使用基于内点反射牛顿法的信赖域算法(subspace trust region method and is based on the interior-reflective Newton method) ; , 3,对于大规模的线性系统使用共觇梯度算法，PCG preconditioned conjugate gradients,。山于这些算法都具有一定的复杂性，具体算法这里不再详述。
2?使用格式:
|vhal.cAhfl4EJBUl卩uh fknlnronf a
ftmaO.
12. 1螺栓组连接的优化设讣
螺栓组联接的设计计算，主要根据被联接机械设备的载荷大小、功能要求和结构特点，确定螺栓的个数，规格和布置方式。螺栓组联接的优化设计,可以在保证机械设备的可黑性和提高寿命的前提下，达到降低成本的LI的。
/
当掾迄的材钳.丘度和割
ifi工艺等丙滦用同时. 出
栓组的成:< 亏耀性牛我町
和直住〃成正比
例如茁所示.有一卜压力客罷內部气休任
强P = g川比?客器内ga r240;w?.螺栓臣中上?(！宣
Q=MO畑!?壕哇村蚪癌5Cr.抗扭强曳
为届眼强良为打=320?”“?挤乐硕劳强魔力=#如伽.许用錢劳安
全翊H?皈住余藏紧力KH采用钢皮石帽密封英二?埶对刚
殘:为K =0?3?占设il戎玄最任的喋治组廡接方廉
一越逢组连撞的设计些甩?H标酉
数約亲条件
T J tfkitx*
?由于集怜纽砂虧*取很于1?怡直径〃和个数打.Kit. ???■???<
uiin f (A ) = (7.=胡=\x
3)蜂遂组联篌的肉康氓件
?"学用洵琪Efh満足痕劳强度年件
_ 26严(K.-叭2丽."1
产冇圧纭亍⑻
■跖甘用条件为*(、?>?[?]-亠」
?②劇劃豹W乐斗为了保证鑼栓之何的SfJfffi力均匀. 砺止局部備气?根齬经复.稣栓的阿距应小干M
:因此，钳到约点条FI ￡小戶泌一附="竺Js/o
- “ *
?巾扶手工岸空耐的義举件亍
为丁悍还《性型拱的滾配工乞陸?细崔之斶的伺》?不虢小F”
?因此.無到旳東条件取(妙=2—竺二2—^^:；0
&心)=_询"
也(”)=一比 V。
nx-f
丿
y
"V
/
Q2”
a却
gm
IMS
U0tl
?Y =
2> U标宙数
_?饶桧述上的优化散学樱掣综上ffWh谊问隊的0学模住为I
心[和订
“II 标乐J n]in = j??]
WMR *&(4闯W4D0zr
—帆"
■ J，
&J T22X：------- S U
自找)=7； gD=“7 <° 三.优化方法和计算结果 调用、IATLAB优化工具箱函数如加如谜行优化计算 1)编制目标函数的文件(Screw_f) %成本目标函数 function f^§cr z ew f(x) f-x(l)*x(2); %工(1)个数，*(2)直径 2)编制约束因数的、1文件(Screwy) 忆约束函数 function [gxeq]-Screw_g(i); p=S.5*Ie6: % 压魂 Ksigm@=4?4;psai=0?23;signi3_l=330;D2=240; %参数 P^l/4*pl*D2A2*1.0e-6*p % 工作压力 F0二(1.6*P+P) ；x(l)：o心单个嫖栓最大工作戦荷 Fl=F0?0?旷PM⑴；％单个螺栓最小工作载荷 d 1 =0.85*x(2); % 螺栓小径 A=l/4*pi*dl A2; %螺栓危险截面积 sisiua nia5：=FO/A; sigma min=Fl/A; §赅贝2_9=($短1|12_|11臥工二gigm合npin〉/2： Sa=(2*sigma 14fKsisma- Dsai)*s^2ma m in)' (fKsis m a + p sa i)* (2 * sig m a a+sisma min)); g(l)=l .7-Sa: %疲劳强度约束条件 g(2>=400*pi/l(l )-H*s(2);味密封约束条件g(3y=-400*pi/x(l>+2*xp)； %装配约束条件 g(4>-16-x(l); %螺栓个数约束条件 g(5)--x(2>; %非负约束条件 .?ca=o； 3)在命令窗口编写调用主(Screw^commond) x0=[7,20]； optioDS-optimsttCdispliiy\ ?iter f); fniincQn(@ Scnm _Am(]打(]?w@ $cr輕岛QJ W QM) disp **********螺栓组优化设计最优解*********** fRJdnir(l ?'螺栓组的个数n-H3.4r<，Ml)) fpiintf(l. ?螺栓的直径d=°o3.4fn'.x(2)) 利用MATLAB解决工程中的实际问题，其具体步骤如下：1)根据设讣要求和U的 定义优化设讣问题，判断优化问题的类型。分析时要区分:单訂标与多口标问题;线 性与非线性问题;是否为线性规划问题等儿种情况； 2)根据优化问题的类型建立数学模型，选定优化函数。3)确定必要的数据和设 计初始点； 4)根据U标函数的性态，预设优化选项； 5)在所有的输入参数定义后，调用优化函数进行优化程序调试；6)根据优化过 程的具体提示信息，修改优化选项的设置,直到达到满足优化函数所需的优化条件为止； 7)对所得数据进行合理性和适应性分析。