机械优化设计实例压杆是一根足够细长的直杆，以学号为p值，自定义有设计变量的尺寸限制值，求在p 一定时d i、d2和丨分别取何值时管状压杆的体积或重量最小？（内外直径分别为d i、d2）两端承向轴向压力，并会因轴向压力达到临界值时而突然弯曲，失去稳定性，所以，设计时，应使压应力不超过材料的弹性极限，还必须使轴向压力小于压杆的临界载荷。

解：根据欧拉压杆公式，两端铰支的压杆，其临界载荷为:I ――材料的惯性矩，EI为抗弯刚度1、设计变量现以管状压杆的内径d i、外径d2和长度l作为设计变量2、目标函数以其体积或重量作为目标函数3、约束条件以压杆不产生屈服和不破坏轴向稳定性，以及尺寸限制为约束条件，在外力为p的情况下建立优化模型：min F〔孟）=oiia 右詡异）=—一和）J2） '压杆的最优化设计也㈤皿也2亟刍-皿。

4勿（忙）=韵佃1 dJ = P -卩匕型—瘩辽（茁—町）I2~M?嵐二（工）==止皿_£］玉o血（兀）=呂.SJ =右 ~ ^lmax —°3) .3 ■■' -J」j -工—二.g$ (光)~ & (£) —^2 2JHK—"」^W = ^W = U-/^ogO劭刘罚函数:反耐皿上严）二7寓-町）f +円｛［诡［o,［cr］-+ mm[ Q@『+nun[鸥『+min[ 0,?]3 + ■■■)传递扭矩的等截面轴的优化设计2、目标函数JT 讪(為4-d「)―^—胡解：1、设计变量:冈区I3、约束条件:式中：――轴所传递的最大扭矩勿（匿）二內⑴二兰拿-罔空J7EZ2）要求扭转变形小于许用变形。

即:式中：G――材料的剪切弹性模数32M T13）结构尺寸要求的约束条件:若轴中间还要承受一个集中载荷，则约束条件中要考虑：根据弯矩联合作用得出的强度与扭转约束条件、弯曲刚度的约束条件、对于较重要的和转速较高可能引起疲劳损坏的轴，疲劳强度校核的安全系数法，增加一项疲劳强度不低于许用值的约束条件。

Jp ――极惯性矩，对实心轴:以轴的重量最轻作为目标函数:a = = --- Ip41）要求扭矩应力小于许用扭转应力，即:抗扭截面系数。

对实心轴16应采用二级齿轮减速器的传动比分配二级齿轮减速器，总传动比i=4，求在中心距A最小下如何分配传动比？设齿轮分度圆直径依次为d i、d2、d3、d4。

第级减速比分别为i 1、i 2。

假设d i=d3,则:1 1 4 1 4」飞3化)+囲+血)］飞如+ 訂二护［卄廿2］七辊矫直实验罚函数法是一种对实际计算和理论研究都非常有价值的优化方法，广泛用来求解约束问题。

其原理是将优化问题中的不等式约束和等式约束加权转换后，和原目标函数结合成新的目标函数，求解该新目标函数的无约束极小值，以期得到原问题的约束最优解。

考虑到本优化程序要处理的是一个兼而有之的问题，故采用混合罚函数法。

-)、优化过程(1)、设计变量以试件通过各矫直辊时所受到的弯矩为设计变量:X -M3- [x1工4 x5]r(2)、目标函数住用在轿直辗上的压力以矫直机的驱动功率为目标函数式中：——矫直速度，mm/s匸一一矫直辊直径，mm耳一一传动效率亠」二一一作用在辊子上的总传动力矩叽=M t+皿+心: --- 轧件弯曲变形所需的转动力矩，N.mmM t=十十亞占+ M补丄一4%M" ' ――克服轧件与辊子间滚动摩擦所需的转动力矩，N.mm山――克服辊子轴承的摩擦及支承辊与工作辊间的滚动摩擦所需力矩，.mm上式表明，编制程序时也可以把目标函数简化为求弯矩和的最小值。

简化问题，可以将程序中的目标函数改为工X =兀+兀$十工4 +屯+兀&（3）、不等式约束首先，试件应满足咬入条件，即式中：丁二一一一、二辊之间的相对压下量口一一试件与矫直辊之间的滑动摩擦系数其次，要保证试件每经过一个矫直单元，实现一次反向弯曲，且弯矩值在极限范围内，即式中：门』一一使试件产生反向弯曲的最小弯矩值，N.mm为了使试件变形充分、均匀，在经过第一、第二个矫直单元时反弯曲率值与原始曲率值应尽量接近。

也就是说，前几个矫直单元采用大变形矫直方案。

试件从最后一个矫直辊中出来后, 要满足质量要求，符合国家有关标准，即有C式中：L门——有关标准规定的残余曲率值，对于本试件[cj= Z.OxlO'5mm~x(4)、等式约束式中：*--、•-＞、・’匸、•’工——分别为相邻两辊之间的相对压下量二)、优化数据以原始曲率半径为2000mm的08F双层焊管为例，得到优化矫直力矩为X 二｛百花花石x.y =(21274.6 22321 9 22045.Z 21983.2 21874.7)7优化矫直力为P={P l马E 耳耳死耳f 二{274.4 820.8 1088.6 1079.5 1062.8 783.4 分别对原始双层曲率半径为1500mm 2500mm和3000mmB勺断面系数相同的5级齿轮传动的传动比分配在指挥仪及精密仪器中，常用如图所示的多级减速器。

为了提高运动精度，不仅要求减重，还要求转动惯量小。

已知总减速比为i,假定各级小齿轮参数相同，各级减速比分别为,且有一：。

下面推导转动惯量和中心距的表达式。

解：i.总转动惯量’」式中一1为小齿轮转动惯量，■- ■-。

若令"-等于下式右端括号内各项，即当1 一定时，对'1求极小，则必为极小，「称传动的转动惯量系数。

2•中心距-进行矫直力矩优化，得到的优化值与原始曲率半径为2000mm的值相差不大。

£二㊁哄力十弓+%+九+258.2)r08F 双层卷焊管式中T为小齿轮分度圆直径。

令-7等于下式右端括号内各项，则T一定下，对-八求极小，必使二亦极小，这意味着重量亦小。

-称中心距系数，即二」-本问题要求减速器重量和转动惯量小，这是双目标优化问题，设计变量为•- ■- ■-，可用线形加权法构造双目标转化为单目标无约束优化问题，即匚:」-式中的=1、笑=1 为权因子，表明设计者对禺（标志转动惯量）和疋2 （标志重量）的重要程度而选定的系数，一般有空一-，例如选二」-，二」- 等。

选取不同□一 - 、二」- 求解，可了解它们对传动比分配的影响，从中选设计者认为满意的解。

曲柄摇杆机构运动规律的最优化设计（复合形法）并且其传动角「（机构的连杆I•与摇杆之间的夹角）的最大值及最小值应分别不大于、不小于其许用值，即『叱站；込］=4罗。

Illi柄掘朴机构简图当曲柄人由其初始位置册转到% 2列规律运动：时，要求摇杆由其极限角” 1开始按下曲柄揺杆机构简图1. 确定设计变量考虑到机构杆长按比例变化时不会改变其运动规律，因此在计算时常取曲柄为单位长度， 即打一 1,而其他杆长则按比例取为 U 的倍数；机架长 二常由结构布置事先给定，分析上图得关 系式：陆]-arccos[啊-arccos[因此，仅为独立变量，是二维最优化设计问题,可根据期望与偏差为最小的要求，来建立目标函数:弘口二旳-瓯尸i-0式中*——期望输岀角1'――实际输出角他 H [ V! i :岬’t h 工网:勺H■禹）盘€俐< 2^由右图求得摇杆匚与BD连线夹角3 . 2 _ 2CL= arccosf—--- --- )S心机架7与BD连线夹角2 ； _ 1 2毘=心七产BD连线长度诩曲将杆『杓EL什楚押聞・h U15A■:J炯空口沁'沆;式中I'' '■ I （单位长度）、二为已定常量。

3.给定约束条件曲柄与机架处于共线位置时:+叮我鮎2=留=arccos[ ]>[r^] = 45°込煜之［8£怙］及^ ' 得约束条件:Si（蜀二一彳-琉-忑心心+&+几r < o自（K）= £ +请-忑心心-（存一许卩< 04.设计变量的可行域曲柄存在的条件:十厶A 斗+/丄 厶+4-4冋右+他一切幼殆©0 = I —丐 <0 刘3)=1-尢沙0务(X)= 71 +.) ■耳■屯兰°S6(^) = 7i -几)+可-码盖o 场 UD - (fl -ZJ- +A 3 <0取单位长度即'1 '■ ^'，只有二、'二:和-J-':为起作用约束，它们是两个椭圆方程。

min f(X) = min £ 仏-陽J_ 2-0 .s,i.g L (X) = _兀：一就- Ji 兀*] +36 <0 乞(刃=兀十理-72%^2-16 £0这是一个带有不等式约束具有两个设计变量的小型最优化设计问题，可采用直接法来求解曲柄摇杆机构如图所示，当曲柄 AB 整周转动时，连杆 BC 上一点M 实现给定轨迹，轨迹坐标方程为可=0.4 + an{ 2jr[(nr - ^I2ir- 0 16]}=■ 2 + 0 5sin(疋-肋要求确定各构件的长度，使点 M 的实际轨迹与给定轨迹间的偏差最小。

2. 建立目标函数:取各构件的杆长1 :为设计变量，即解：1.确定设计变量：T" = *1丿2丿*4人山『得约束条件:最后的数学模型为按连杆上点M实际实现的轨迹与要求实现的轨迹的均方偏差最小的原则建立目标函数，其函数形式可表达为，'"'(町一心尸+ 3-丿』一式中，(s表示将曲柄转角■分成s等份)；一'丁，，山为连杆上点M第j个位置所要求实现的坐标值；*： ^为点M实际实现的坐标值与I；可按下列公式计算得到：卞：=CGS密』+ '二g ¥$一%sin 丫$yj= z L sin 巒+■?」tin y,4-^cosy = +^3 cos^r - s L cossin y = (z l sin 爭一巧sin 叭"】-丘±(肝-4加严CQSkT = --------------2Z4 = 2?+以B=2KLy C^K^-M2K工2^ -23J j - 22-154 COS (pL = 2広样斗一2S］巧cos tp其中第三式中的根号前的正负号应参照运动连续性来选取。

3.确定约束条件：首先要保证曲柄AB能作整周回转。

若结构上要求'' ，则可取下列约束条件(曲柄存在条件)：gj® 二 6 一勾一屯+^4< 0也(刃=习+石一习-2A<0其次考虑机构具有较好的传力性能，按最小传动角大于」：“写成的约束条件为gj(Z)二习'+空：—(习一习尸—Ng屯cos45° < 0g^(Z) =-z23- r53+ (z4+ Z-!)3- 2z a z3co.'45° 兰Q综上所述，这一优化设计问题的数学模型可表示为求设计变量 - 一-一一亠一一一_，使min y（z）- £〔勺・兀» +Cx广打）满足于约束条件■- -4.计算结果:以u - -为例，计算结果为:简化的机床主轴当跨距1,夕卜伸端a求轴的内、外径d、D取何值的时候，主轴重量最轻?解：1、确定设计变量由于材料一定时，主轴内孔只与机床型号有关，所以设计变量为:x = [ x i X2 X3]T =[l D a]T 2、目标函数：考虑主轴最轻，所以min二—（D 2 一沪）（f 十盘）’p -（xf —d ‘）（厲十与）44式中：p ——材料密度3、约束条件：由于主轴刚度是一个重要的性能指标，其外伸端的挠度y 不得超过规定值 y 。