机械优化设计实例

1、压杆的最优化设计I材料的惯性矩，1、设计变量 现以管状压杆的内径机械优化设计实例压杆是一根足够细长的直杆，以学号为p值，自定义有设计变量的尺寸限制值，求在p一定时di、d2和l分别取何值时管状压杆的体积或重量最小（内外直径分别为di、d2）两端承向轴向压力，并会因轴向压力达到临界值时而突然弯曲，失去稳定性，所以，设计时，应使压应力不超过材料的弹性极限，还必须使轴向压力小于压杆的临界载荷。解：根据欧拉压杆公式，两端钱支的压杆，其临界载荷为：EI为抗弯刚度di、外径d2和长度l作为设计变量2、目标函数以其体积或重量作为目标函数3、约束条件以压杆不产生屈服和不破坏轴向稳定性，以及尺寸限制为约束条件，

2、在外力为p的情况下建立优化模型:2）min=Enin尸（由,4）=/加二目式心.由）=产一下一=产一一市一之。*（X）=且4）=dg-dx<0口=g<*-d1UHK-U目耳（尤）=y2）=dn皿一M。目式X）nW6（心）二心一2mue-0g式灯二g式9二A一收。罚函数:口3卜心上，闻）=-&：）£+"训【Q可产71+讪。*一,篇也）_/4打0二一山）64/+minO/iF+mm0,4族+nunOJ2+解：1、设计变量:传递扭矩的等截面轴的优化设计2、目标函数0）-F（J0过一以轴的重量最轻作为目标函数：3、约束条件：r=<r1）要求扭矩应力小于许用

3、扭转应力，即：3r监式中：一一轴所传递的最大扭矩¥d7抗扭截面系数。对实心轴16小=幻二-订过07Z32）要求扭转变形小于许用变形。即:扭转角式中：G材料的剪切弹性模数Jp极惯性矩，对实心轴：rsaCr32一之03)结构尺寸要求的约束条件:若轴中间还要承受一个集中载荷,束条件、弯曲刚度的约束条件、则约束条件中要考虑：根据弯矩联合作用得出的强度与扭转约对于较重要的和转速较高可能引起疲劳损坏的轴,应采用疲劳强度校核的安全系数法，增加一项疲劳强度不低于许用值的约束条件。二级齿轮减速器的传动比分配二级齿轮减速器，总传动比i=4,求在中心距A最小下如何分配传动比设齿轮分度圆直径依次为di、d2

4、、d3、d4。第一、二级减速比分别为ii、i2。假设di=d3,则：11414公#4+或）=-1+“+”=5g+1+2七辐矫直实验罚函数法是一种对实际计算和理论研究都非常有价值的优化方法，广泛用来求解约束问题。其原理是将优化问题中的不等式约束和等式约束加权转换后，和原目标函数结合成新的目标函数，求解该新目标函数的无约束极小值，以期得到原问题的约束最优解。考虑到本优化程序要处理的是一个兼而有之的问题,故采用混合罚函数法。P7一设Ps7,|Pl一)、优化过程(1)、设计变量以试件通过各矫直辐量：X=MMP4MlP2作用在瓣直辐上的压力时所受到的弯矩为设计变mJ=上工(2)、目标函数以矫直机的驱动功

5、率为目标函数式中:矫直速度， mm/s一,矫直辐直径，mm-传动效率%作用在辐子上的总传动力矩轧件弯曲变形所需的转动力矩,%=（%+%+肛+&）克服轧件与辐子间滚动摩擦所需的转动力矩,明克服辐子轴承的摩擦及支承辐与工作辐间的滚动摩擦所需力矩,简化问题，可以将程序中的上式表明，编制程序时也可以把目标函数简化为求弯矩和的最小值。目标函数改为ZX=+*&+(3)、不等式约束首先，试件应满足咬入条件，即T九式中： ：、二辐之间的相对压下量其次，要保证试件每经过一个矫直单元，实现一次反向弯曲，且弯矩值在极限范围内，即式中：'八一一使试件产生反向弯曲的最小弯矩值，M1nBi=min

6、/为了使试件变形充分、均匀，在经过第一、第二个矫直单元时反弯曲率值与原始曲率值应尽量接近。也就是说，前几个矫直单元采用大变形矫直方案。试件从最后一个矫直辐中出来后，要满足质量要求，符合国家有关标准，即有式中：有关标准规定的残余曲率值，对于本试件cj=2.OxIO-3-1(4)、等式约束=尸声应满足:工,，.二式中：了门、734、了45、1比分别为相邻两辐之间的相对压下量二)、优化数据以原始曲率半径为2000mm的08F双层焊管为例，得到优化矫直力矩为X=马马x3匕三产=21274622321922045221983一221874,7)7优化矫直力为P=品招名片写"耳=(274.4S2

7、0.S1092.61079.510626783.42f分别对原始双层曲率半径为1500mm、2500mm和3000mm的断面系数相同的08F双层卷焊管进行矫直力矩优化，得到的优化值与原始曲率半径为2000mm的值相差不大。5级齿轮传动的传动比分配在指挥仪及精密仪器中，常用如图所示的多级减速器。为了提高运动精度，不仅要求减重，还要求转动惯量小。已知总减速比为i,假定各级小齿轮参数相同，各级减速比分别为“，也“工4，15,且有工=".匕工"5。下面推导转动惯量和中心距的表达式。解：1.总转动惯量/士T常产181a1巧1工彳q口二30北1+T+TT+-T-T+.J.2.2+-33

8、1瓦L瓦/%十%式中为小齿轮转动惯量,V国内。尸。若令 陌 等于下式右端括号内各项，即当一定时，对也求极小，则心必为极小，局称传动的转动惯量系数。2.中心距八乱云十式中式为小齿轮分度圆直径。令笔等于下式右端括号内各项，则一定下，对左2求极小，必使上亦极小，这意味着重量亦小。%称中心距系数，即%=i本问题要求减速器重量和转动惯量小，这是双目标优化问题，设计变量为，可用线形加权法构造双目标转化为单目标无约束优化问题，即匚式中叼=1、%!=1为权因子，表明设计者对跖（标志转动惯量）和/（标志重量）的重要程度而选定的系数，一般有我!=1,例如选的=1，%=1等。选取不同%=1、%=1求解，可了解它们对

9、传动比分配的影响，从中选设计者认为满意的解。曲柄摇杆机构运动规律的最优化设计（复合形法）当曲柄1】由其初始位置职口转到朝2时，要求摇杆由其极限角啊。开始按下列规律运动：并且其传动角/（机构的连杆,口一与摇杆片之间的夹角）的最大值及最小值应分别不大于、不小于其许用值，即了皿卬皿=13$；&之71dlJ=45口曲柄摇杆机构僮图曲柄播杆机构简图1.确定设计变量考虑到机构杆长按比例变化时不会改变其运动规律，因此在计算时常取曲柄为单位长度，即=1,而其他杆长则按比例取为1的倍数；机架长式：心十一/J十研二维八足一"工>2（11十，口），4g二arccos：上卬*因此，仅k,右.为

10、独立变量，是二维最优化设计问题，叱H：2.建立目标函数可根据期望与偏差为最小的要求，来建立目标函数：近幻二£3七尸i-a式中Wj期望输出角%*实际输出角/笈一以一房0金仍<:需舞“一吗+/郭<甲&<.2况暂由右图求得A¥“凝&摇杆右与BD连线夹角2,2_2%=arccosf-L）3七机架”与BD连线夹角，雕一.灯常由结构布置事先给定，分析上图得关系上UA.1/4g组lB吊塔肯轩也何EW叶量用田/?j=arccosfpt=+/J-2/Jjcos物BD连线长度匕"1414式中h（单位长度）、L为已定常量。3 .给定约束条件曲柄与机架处

11、于共线位置时:了咏:arced+勺7H卬蚪卜充匚丫miinarccosfXA>r_ = 450cosy之上佻产叩及1'-*-得约束条件：gl（©=7；一冠-应工】工”&+"<0g式X）=+岩-五西电一圈一，1404 .设计变量的可行域曲柄存在的条件：'口之'口4之1，4之k+。汕+AA+&T）之八幻WO=1一'M0公阳”松工04（X）=区+4）-a-/三。4（x）=©-14）+/一巧M。得约束条件：；1:一'-j取单位长度即4=?4=$，只有&（为和自式工）为起作用约束，它们是两个椭圆方程

12、。minf(X)=min£(巴-%丫2-0&t.g(X)=一工；一工;一41x2+36<0最后的数学模型为g2(Y)=xf+Xj-&x1A16M0这是一个带有不等式约束具有两个设计变量的小型最优化设计问题，可采用直接法来求解曲柄摇杆机构如图所示，当曲柄AB整周转动时，连杆BC上一点M实现给定轨迹，轨迹坐标方程为xg=0.4+sin2*（川-中）/2"一Cl.16=2+09sin（窿一甲j要求确定各构件的长度，使点M的实际轨迹与给定轨迹间的偏差最小解：1.确定设计变量：取各构件的杆长为设计变量，即N=勺/3,/，之4.35，/=/">&a

13、mp;J4匕。F2.建立目标函数：按连杆上点M实际实现的轨迹与要求实现的轨迹的均方偏差最小的原则建立目标函数，其函数形式可表达为，''''1i = £ （心-工学+尸j式中， （s表示将曲柄转角2厘分成s等份）；,豆?廿为连杆上点m第j个位置所要求实现的坐标值；勺Vj为点M实际实现的坐标值。看与可按下列公式计算得到：%=可cos（Pj-Fz5丫金-z6siri匕=为sin甲j+句sin门+/cosy,cosy=（z*+3aosr-zlcos（p）f22sinV=（町sin-一勺脸f名a一8±（£"产gs审=2金达二不+财B

14、=2KLC=Q-守33-22*K.=-2a十%十"$-2zt24COSL-27374-2zxzcostp小岫=2马句sinP其中第三式中的根号前的正负号应参照运动连续性来选取。3 .确定约束条件：首先要保证曲柄AB能作整周回转。若结构上要求飞，则可取下列约束条件（曲柄存在条件）：务（Z）=4-十0M式Z）=向+马一勺MQ其次考虑机构具有较好的传力性能，按最小传动角大于45°写成的约束条件为宫<2)=nJ-(z4百产-2%cos450<04(2)=-z2a-eJ+(z4+zj”-2与与cof45qW0综上所述，这一优化设计问题的数学模型可表示为求设计变量N=L，与

15、，%，二4必，去F,使miny(z)X3一%)-y)一满足于约束条件'广-：4 .计算结果：以邙°)二口4.5453引为例，计算结果为：简化的机床主轴当跨距1,外伸端a求轴的内、外径d、D取何值的时候，主轴重量最轻解：1、确定设计变量由于材料一定时，主轴内孔只与机床型号有关，所以设计变量为：x=x1X2X3T=1DaT2、目标函数：考虑主轴最轻，所以manPg=三(0。/)(£+0),p二阻J：/)5十七)44-一式中：p材料密度3、约束条件：由于主轴刚度是一个重要的性能指标，其外伸端的挠度y不得超过规定值y0所以可依此建立性能的约束：.二1=二（口4_1）在外力F给定的情况下，挠度3EI,64/、孑（无】+不壬）.g=4-/-A之0所以：丁''-'强度有富裕，所以这里可以不考虑应力约由于机床主轴对刚度要求比较高，当满足刚度要求时,束条件。边界的约束条件由具体机床结构与主轴材料而定：1.<1<