机械优化设计三个案例.

机械优化设计案例11.题目对一对单级圆柱齿轮减速器，以体积最小为目标进行优化设 计。2已知条件已知数输入功p=58kw，输入转速ni=1000r/min，齿数比u=5, 齿轮的许用应力、：H=550Mpa，许用弯曲应力:F=400Mpa。 3建立优化模型3.1问题分析及设计变量的确定由已知条件得求在满足零件刚度和强度条件下，使减速器体 积最小的各项设计参数。由于齿轮和轴的尺寸(即壳体内的零件) 是决定减速器体积的依据，故可按它们的体积之和最小的原则建 立目标函数。单机圆柱齿轮减速器的齿轮和轴的体积可近似的表示为：2 2 2 2 2 2v =0.25*0 -dz1) 0.25二b(d2 -dz2) -0.25(b -c)(Dg2 -dg2)-d0c 0.25二l(d； d；2) 7二d； 8：d；22 2 2 2 2 2 2 2= 0.25叫m z b _dz初+m z u b _dz2b-0.8b(mzu10m) +2 2 2 22.05bdZ2 -0.05b(mzu -10m -1.6dz2)+dzd + 28dz1 +32dz2】式中符号意义由结构图给出，其计算公式为d1 = mz, d2 = mz2Dg2 二 u m1z10mdg2 =1.6dz2,d0 =0.25(u m1z10m-1.6dz2)c =0.2b由上式知，齿数比给定之后，体积取决于b、Z1、m、l、dz1和dz2六个参数，则设计变量可取为x 二X1 X2 X3 X4 X 冷丁 =b 乙 m I dz1 dz2T3.2目标函数为2 2 2 2 2 2 f (x 0.785398(4.75x1x2x385x1x2x3 -85x1x30.92x1xx1x52222220.8x1x2x3x -1.6x1x3x6x4x5x4x628x532x6)；min3.3约束条件的建立1)为避免发生根切，应有 乙-為山=17，得gi（x） =17 _X2 乞 0:? .： b .:min 二 T max mCOCD2 ）齿宽应满足d , Fin和Fax为齿宽系数d的最大值和最小值，一般取；：min =0.9, max=1.4,得g2（x） =0.9 -为（X2X3） _0g3（x） =x1;（X2X3） -1.4 乞 03）动力传递的齿轮模数应大于 2mm,得g4（x） =2-X3 乞 04）为了限制大齿轮的直径不至过大，小齿轮的直径不能大于d 1 max彳得g5（x） 7x3 -3 0 005） 齿轮轴直径的范围：dzmin - dz /x得g6（x） =1 0 0X5 _0g7（x） =X5 -1 5 00g8（x） =1 3 0X6 岂0g9（x）=冷 一2 0 0 06） 轴的支撑距离1按结构关系，应满足条件：1 - b 2：伽 0.5dz2（可取比min =20），得g10 （x）二 X1 0.5x6 - X4 - 40 _ 0 7）齿轮的接触应力和弯曲应力应不大于许用值，得g11 （x） =1468250. （x2x3 _ % ） -550 - 0q12（x） =70982q-4 2 一400 012X1X2X3 （0.169 0.6666 10 x 0.854 10 X2）g13（x）=7098224 2 -400 空01X1X2X3 （0.2824 0.177 10 x0.394 10 X2）8）齿轮轴的最大挠度；max不大于许用值，得g14（x） =117.04x；.（X2X3X4） -0.003x4 乞 09）齿轮轴的弯曲应力；w不大于许用值w，得gi5（x）=厶（2.850 生）2 +2.4X101-5.5 兰0X5 VX2X3| 61 J 2.85 10 沧2 丄c12gi6（x）=p 1（） +6X10 5.5 兰 0X6 tX2X34优化方法的选择由于该问题有6个设计变量，16个约束条件的优化设计问题， 采用传统的优化设计方法比较繁琐，比较复杂，所以选用Matlab优化工具箱中的fmincon函数来求解此非线性优化问题，避免了 较为繁重的计算过程。5数学模型的求解将已知及数据代入上式，该优化设计的数学优化模型表示为:2 2 2 2min f （x 0.785398（4.75x1x2x385x1x2x3 -85x1x32220.92x1xx1x50.8x1x2x3x -1.6x1x3x6 x4x52 2 2X4X628x532X6）Subject to:g,x） =17 - X2 乞 0g2（x） =0.9 -为.（X2X3）乞 0g3（x）=为.（X2X3） T.4 乞0g4（x） =2 -X3 岂0g5（x）二 X2X3 -300 - 0g6（x） =100 -X5 乞 0g7（x） = X5 T50 - 0g8（x） =130 -X6 -0g9(x) =X6 200 _0gio (x) = Xi 0.5x6 - X4 - 40 0= 1468250 (x2x3 . xj -550 乞 0= 7098224 2 -400 _0/x1 x2x3 (0.169 + 0.6666 x 10 x2 0.854 X10 x；)二709822/ 2 -400 0/X1X2X3(0.2824 +0.177 x10 x；0.394x10 x；)=117.04x：(X2X3X4)-0.003x4 _0gii(x)gi2(x)gi3（x）gi4（x）12.85 106x4)2 2.41012 5.5 _0g15(x)62.85 10 x4X2X3)2 6 1012 5.5 _0运用Matlab优化工具箱对数学模型进行程序求解首先在Matlab优化工具箱中编写目标函数的 M文件myfun.m, 返回x处的函数值f:function f = myfun(x)f=0.785398*(4.75*x(1)*x (2) A2*x (3)A2+85*x(1)*x (2)*x (3)A2-85 *x(1)*x(3)A2+0.92*x(1)*x (6) A2-x(1)*x (5F2+0.8*x(1)*x (2) *x (3) *x(6 )-1.6*x(1)*x (3) *x (6)+x *x (5) A2+x *x (6) A2+28*x (5F2+32*x ( 6)a 2)由于约束条件中有非线性约束，故需要编写一个描述非线性约束条件的 M文件mycon.m：fun ctio n c,ceq=myobj(x) c=17-x(2);0.9-x(1)/(x(2)*x(3);x(1)/(x(2)*x(3)-1.4;2-x(3);x(2)* x(3)-300;100-x (5) ;x (5)-150;130-x (6);x (6)-200;x(1)+0.5*x (6)-x -40; 1486250/(x (2)*x (3)*sqrt(x(1)-550;7098/(x(1)*x (2) *x (3) A2*(0.169+0.006666*x (2)-0.0000854*x (2)八2)-4 00;7098/(x(1)*x (2)*x(3)A2*(0.2824+0.00177*x (2)-0.0000394*x (2)八2) )-400;117.04*x(4)A4/(x (2)*x (3) *x(5F4)-0.003*x (4) ;(1/(x (5) A3)*sqrt( (2850000*x /(x( 2)*x(3)A2+2.4*10A12)-5.5;(1/(x (6) A3)*sqrt(2850 000*x /(x( 2)*x(3)A2+6*10A13)-5.5;ceq=;最后在command window里输入：x0=230;21;8;420;120;160;% 给定初始值x,fval,exitflag,output=fmi ncon(myfun, xO,my obj,output)%调用优化过程最优解以及结果分析运行结果如下图所示:Command Windowx =123.356599. 8E171. 7561147. 3757150.4904129. 5006fval 二2.3163e+007exitflag =0output -iterationE： 42furcCeunt: 604stepsize;1algorithn:n jnedium-Etale: SQP, Quasi-Nawton, line-search?firstordetopt:L9950e+OO7cgiterations:nescase:1k79 char由图可知，优化后的最终结果为x=123.356599.85171.7561147.3157150.4904129.5096f(x)=2.36e*107由于齿轮模数应为标准值，齿数必须为整数，其它参数也要进行圆整 ,所以最优解不能直接采用，按设计规范，经标准化和圆 整后：x=124 100 2 148 150 130f(x)=6.16 *10 76. 结果对比分析若按初始值减速器的体积V大约为6.32 xi07mm3,而优化后 的体积V则为6.16X 107mm3,优化结果比初始值体积减少为： V = 1 (6.16X 107/6.32 X 107)X 100% = 2.5% 所以优化后的体积比未优化前减少了 2.5%，说明优化结果相 对比较成功。7. 学习心得体会 学习机械优化设计课程的心得体会 通过将近一学期的学习，对这门课有了初步的了解和认识， 学期伊始，浏览全书，发现全是纯理论知识，觉得这门课会很枯 燥，但是又回过头来想想，作为 21 世纪的大学生，要使自己适应 社会需求，首先在做任何事之前都应该有正确的态度看待问题， 把这些想法作为促使自己进步的动力，再去学习课本知识，效果 应该很不一样，有了想法就付诸行动，随着对课本内容的学习跟 老师的讲解，发现并不是像自己在学期初想的那样困难，特别是 在老师介绍了一些与机械优化设计相关的计算机语言和计算机软 件后，真正体会到科学优化设计的强大跟简洁明了，与传统优化 设计方法相比较，大大提高了设计效率和质量。传统设计方法常在调查分析的基础上，参照同类产品通过估 算，经验类比或试验来确定初始设计方案，如不能满足指标要求， 则进行反复分析计算性能检验参数修改，到满足设计指标要 求为止。整个传统设计过程就是人工凑试和定性分析比较的过程， 是被动地重复分析产品性能，不是主动设计产品参数。按照传统 设计方法做出的设计方案，有改进余地，但不是最佳设计方案。而现代化设计工作是借助电子计算机， ，应用一些精确度较高 的力学数值分析方法，优化软件进行分析计算，找最优设计方案， 实现理论设计代替经验设计，用精确计算代替近似计算，用优化 设计代替一般的安全寿命可行性设计。在进行程序求解的过程中，因为是初学 Matlab 软件，对很多 问题的关键点不能够掌握，非线性约束如何书写，上、下限如何 选择，函数格式如何书写，变量未定义等等或大或小的问题，但 是在一步步排除错误、重新编写程序的过程中，渐渐的对 Mtalab 熟悉起来，懂得了一些优化方法的简单计算过程和原理，省去了繁琐复杂的优化计算过程 在学完课程之后，反思自己在学习过程中的得失，深深体会 到，不论在人生的哪个阶段，都要对自己负责，做任何事都要耐 心，细致，“千里之行，始于足下” ，学会在物欲横流的社会大潮 中，坚持踏踏实实走好人生的每一步。8. 参考文献1 孙靖民 ,梁迎春 . 机械优化设计 . 北京 :机械工业出版 社,2006.2 濮良贵,纪名刚. 机械设计 . 8版. 北京:高等教育出版社,2006.3 孙桓,陈作模,葛文杰 . 机械原理. 7 版. 北京:高等教育出版 社,2006.4 李涛,贺勇军,刘志俭. MATLAB 工具箱应用指南应用数 学篇M.北京:电子工业出版社，2000.机械优化设计案例 2复杂刀具优化设计数学模型的建立及算法改进摘 要 : 目的 建立复杂刀具优化的数学模型 , 提高优化算法速度 . 方法 采用优化设计与 CAD 相结合的方法 . 结果与结论 解决了传统刀具 设计的缺点 , 改进后的算法速度大幅度提高 .关键词 : 数学模型；优化；算法在传统的刀具设计中 , 通过查表和经验公式来确定各种结构参数和几何 参数 ,然后,反复计算来得到相对较优的刀具参数 . 这种方法使设计过程复杂 费时,且得不到最优化的参数 ,设计出的刀具成本高 ,加工效率低 . 因而刀具 的计算机辅助设计应采用优化设计与 CAD相结合的方法，欲进行优化设计， 必需首先建立刀具优化设计的数学模型 , 由于复杂刀具的种类繁多 , 结构变 化多样 , 优化目标不同 , 因而需分门别类地建立模型 1 ，此篇仅以轮切式拉刀 为例 .1 拉刀优化设计的数学模型在拉刀参数设计过程中需要选择的主要参数有拉削余量 A,齿升量af, 齿距 t, 容屑槽形状和深度 h, 容屑系数 k, 同时工作齿数等 , 这些参数可分为 两类 , 一类是独立参数 , 如拉削余量和容屑槽形状等 , 这些参数基本不受其他 参数的影响 .另一类参数是非独立参数 ,如齿升量、 齿距、容屑槽深度、 容屑 系数等 , 这些参数既相互限制又相互依赖 , 第一类参数的选择比较容易 . 可以 用经验公式和数据库来解决 .第二类参数比较复杂 ,只有通过优化的方法才 能得到较好的结果 .粗切齿升量的选择是一个比较复杂的问题 增大af可使齿数减少，拉刀长度变短 , 但同时又要求容屑槽深度增加 . 另外齿升量的增加又会引起拉削 力的增大 ，受到拉床和拉刀拉应力的限制 .齿距是决定拉刀长度的一个重要因素 ，t 越大，拉刀越长 ，同时工作齿数 越少 .这样会在拉削过程中引起振动 ，生产效率低 ，降低刀具的使用寿命 ;t 过小 ， 又会使容屑空间变小 ， 从而限制了齿升量的增大 .其他参数如同时工作齿数 z i ， 容屑槽深度 h， 容屑系数 k 都是 a f 和 t 的函数，只有当af和t选择后才能确定.从上述参数分析可知，a f和t是拉刀设计的关键，在af和t之间应有一最佳组合值，使得af在拉床的额定应力和 拉刀的许用应力范围内达到最高 ， 即使拉刀的长度最小 .1.1 目标函数的建立确定以af和t为优化的自变量,A为切削余量拉刀长度是与拉削生产率、 成本及其工艺性能有关的参数 ，拉刀越短对使用和制造越有利 ，因而取粗切齿 部分长度 L 作为优化目标 4F= minL(a f ，t) =tA/(2a f ).(1) 1.2 约束条件的建立1) 容屑槽空间的限制h- 1.13(ka f L w ) 1/2 0(2)式中 h 是与 t 有关的参数 ;k 为容屑系数 ， 是与 t 和 af 有关的参数 ;Lw 为 拉削长度 .2) 拉床额定拉力的限制F e -p n D wz/zc 0.(3)式中 Fe 为拉床额定拉力 ;Dw 为拉削后孔直径 ;p 为单位切削力 ;zi 为同时工 作齿数 ，zi=INT(Lw/t)+1;zc 为组齿数 .3) 拉刀许用拉应力的限制(T 2pD wzi/zcdmin0.(4)式中 T 为拉刀许用拉应力 ;dmin 为拉刀最小直径 .4) 最大同时工作齿数的限制11 -z i 0.5) 最 小(6)同时工作齿数的限制zi-3 06)最大齿距的限制25-t0(7)7)最小齿距的限制t-40(8)8) 弧形槽能保证稳定的分屑要求的最大齿升量h-a f O；f(D,n z ,z c ) -a f 0.(9)9) 齿 距 应 为 0.5 的 整 数 倍t- Int(2t)/2 = 0.(10)1.3 优化模型自变量:af ,t;目标函数 :F=minL(a f ,t)=tA/2a f ；约束方程：g(1)=h-1.13(kafLw) 1/2 0;g(2)=Fe-p n D wzi/z角 0；g(3)= (T -2pD wzi/zcdmin 0Q(4)=11-z： 0;g(5)=zi -3 0;g(6)=25-t0;g(7)=t-40;g(8)=h-a f 0；g(9)=f(D,n z,z)af 0;g(10)=t-Int(2t)/2=0.2 优化算法2.1标准算法复合形法是一种采用直接搜索方式求解非线性规划问题的数值计算方法，这个方法可以在 N维非线性约束的空间中自动选择并改进设计点，该方法的一般步骤为2.3:1) 在可行域内生成 mn+1个点xj(i=1, 2, ,n,n+1, , ,m)构成初始复合形，这里需要注意两个问题：初始顶点的形成，可以人工选定， 也可随机产生；需要检验初始顶点是否满足约束条件，即检验其可行性；2) 计算各顶点的目标函数值，将其由小到大的顺序重新编号，f(x 1) f(X 2) ， f(X m)；3) 确定除去最坏点x m后复合形中其余 m-1个点的中心点，即1兰xc= xi ;m-1 m 14) 确定最坏点 xm对中心点xc的映射点xa=xc+ a (x c-x m)， a为映 射系数，一般取0.91.3;5) 检验映射点x a的可行性：如果违背了某个约束条件，则按(x c +xa )/2 xa ,把映射点xa向中心移动一半距离，反复直至映射点 x a是可行点；6) 计算新的可行点的函数值f(x a ),用它代替最坏点 x m ,完成一次迭代，回到第二步；7) 重复以上过程，直到满足f(x m )-f(x1) e ,则终止.以上称为复合形法的“标准算法”，由于该算法的概念简单、容易实现，且能有效灵活地处理 不等式约束问题，所以在结构化设计中得到广泛的应用 .2.2 存在的问题把上述标准算法应用于工程实际时 , 就会发现它还存在以下几个问题 :1) 过多的可行性检验限制了其在优化设计中的有效应用 . 初始顶点生成和 映射点的确定 , 都要进行可行性检验 , 在结构优化设计中 , 可行性检验其实质 上就是结构分析过程 ,其计算量通常要占总工作量的 80%以上, 因此结构分 析次数过多 , 必然会导致因计算时间过长而降低算法的效率 .2) 迭代过程中向极值点逼近的速度问题 .开始若干次迭代 ( 一次迭代是对于 选取一个既满足约束条件又使目标函数值有所改善的新点所需的计算 ), 目 标函数值下降得很快，各顶点迅速接近极值点，一般来说，最初的(510)次 迭代函数值下降得最快 .随着迭代次数的增加 ,函数值的变化却越来越缓慢 , 也就是说 ，这时要使目标函数值有微小的改善 ，都要付出宝贵的计算时间 .3) 局部最优点问题 .用上述算法得到的最优点有可能是局部最优点 ，虽然可 通过多取几个初始点 ， 经计算后得到几个最优点 ， 然后比较得到全局最优点 ， 但这样必然会导致计算工作量的成倍增加 .2.3 分层复合形法针对标准算法中存在的问题 ， 采用“分层复合形法” ，它是对标准复合形 法的改进 ， 其基本思想是 : 充分利用复合形法开始时目标函数值急剧下降的 特点 ， 以迭代次数为控制参数 ， 进行两层优化计算 ， 为避免产生局部最优点 ， 在第一层迭代中 ， 选取多组复合形分别地进行计算 ， 经过若干次有效地迭代 ， 各顶点迅速地逼近最优点 ， 分布在最优点附近 . 分层复合形法的基本步骤如 下:(ng)1) 选择ng组初始顶点x i(), x i, (i=1,m)构成讥个初始复合形 ，ng =Intn/2+1，n 为设计变量数 . 这里只要设计变量所取的值 不太小且相互间离得远些 ， 就可不对初始顶点作可行性检验 .2) 对各初始复合形标准算法第2) 6)步的计算是第一层迭代 ， 取映射率为a 1,迭代次数为n t1;3) 当各复合形都迭代nt1次以后，第一层迭代结束，取两个最好设计点组成新的复合形进入第二层迭代 ，取映射率为 a 2 (a 2a 1)，迭代次数为 n t2 ;第二层迭代得到的最优点可被认为全局最优点 .分层复合形法有以下几个优点 : 迭代次数大大减少 ; 以迭代次数为停 止准则 , 可根据需要人工控制计算工作量 ; 第二层迭代能有效地产生全局 最优点 .3 结论依据本文所述方法，已开发出具有高效率优化CAD系统,证明对传统算法的改进是有效的参考文献 :1 唐锡荣 . CAD/CAM 技术 M. 北京 : 北京航空航天大学出版社 , 1994.18-36.2 蔡锁章 . 计算方法 M. 北京 :中国科学技术出版社 , 1993. 54-60.3 徐灏.机械设计手册 . 第二卷 M. 北京: 机械工业出版社 , 1991. 40-41.4 吴伏家，刘兆华.圆孔拉刀 CAD系统研制J.华北工学院院报，1996, （ 增刊 ）: 74-78.机械优化设计案例3直齿圆柱齿轮传动的优化设计摘要：一、问题描述：现有一单级渐开线直齿圆柱齿轮减速器，其输入功率N=280kW输入转速ni=980r/min，传动比i=5。小齿轮为实体结构， 大齿轮为腹板式结构 （带 有四个减轻孔），两齿轮各部分尺寸的符号如图一所示：原用常规设计方法的设计结果为：齿宽B=B=13cm,小齿轮齿数zi=21，模数 m=0.8cm, I i=42cm, dsi=12cm, ds2=16cm。现要求在保证承载能力的条件 下，通过优选上述有关参数，使减速器的体积达到最小。二、建立优化设计目标函数：齿轮传动优化设计中，设计变量一般选为齿轮传动的基本几何参数或性 能参数，例如齿数、模数、齿宽系数、传动比、螺旋角、变位系数和中心距 分离系数等。齿轮传动的优化目标，较常见的是体积或质量最小，传动功率最大，工作寿命最长，振动最小，启动功率最小等。现在选体积最小为优化目标，而减速器的体积主要是取决于内部零件（齿轮和轴）的尺寸大小，在齿轮和轴的结构尺寸确定之后，箱体的尺寸将随之确定，因此将齿轮和轴的总体积达到最小作为优化目标。减速器内部有两个齿轮和两根轴，为了简化计算，将轴视为光轴，则有V = Vsi Vs2 Vgi Vg2:2 - -2二 dsi（li l3 ds2 （li I2）44 2 2- -2 2 : 2 2 : 2（di-dsi ）B （d2-ds2）B2 - （D2-Di ）（B2-C）-4（ d。C）4444式中：Vs1 , Vs2两轴体积，cm;3Vg1 , Vg2两齿轮体积，cmds1, ds2 两轴的直径，cmli , I2, I3轴的长度，cm；d1 , d2 两齿轮的分度圆直径，cm 4 =mz!, d2；m 两齿轮的模数，cm；B-i, B2两齿轮的宽度，近似取B = B2 = B ,cm。根据结构设计经验公式，齿轮各部分尺寸关系为：、二 5mD1-1.6ds2C =0.2BD2 =d2 -2、d0 =0.25（D2 -D；）并取：l2 =32cm l3 =28cm优化设计中的设计变量取为：X 二为，X2,X3,X4,X5,xJT 二B,Zi,m,li,dsi,ds2T将目标函数整理后得到：f （X） =0.78539815（4.752乂32(0.169 0.6666 10x2 -0.854 10，x22)-0 g13(X217098/ x1x2x32(0.2824 0.177 10x0.314 卩冷2)-。(9) 主动轴刚度条件PI34Jf式中：P 作用在小齿轮上的法向压力,N , P =2Mj / mz cosa ,其中为齿轮压力角，取=20 ；J 轴的惯性矩，对圆形剖面,J - -：dsi4 /64 ；E 轴材料的弹性模量，E =2 105MPa ；f 轴的许用挠度，取f= 0.003b所以可以得到:g/X) =0.003x4 -0.01233x43X2X3x/ 一0(10) 主动轴的弯曲强度条件:耳斗式中：T 轴上的扭矩， T M,;轴 上 的 弯 矩，N -cm ,M =pk 业 29050上；2 mz1 cos ：mzi:考虑扭矩和弯矩的作用性质差异的系数，取:=0.58 ；3si b轴的许用弯曲应力， 卜巾=55MPa ；W 轴的抗弯剖面系数，对实心轴W :“ 0.1d带入各参数，并整理得:g15(X) =55-x3” (29050 X4 )2 (0.58 27300)2 -0XX2X3(11) 仿照前面的处理方法可得从动轴弯曲强度条件:gi6(X) =55-X6”(29050 人)2(0.58 5 27300)2 _ 0VX2X3总结上述各式，可得到优化设计的数学模型为:min f (X)X 三 E6s.t gj(X) 一 0 j =1,2,.,16即一个具有十六个不等式约束的六维优化问题。四、优化方法选择及优化结果：1、采用MATLAB工具箱进行优化首先在当前 MATLAB的工作目录下建立目标函数文件myfun.m文件:fun ctio n f =myfun(x)f=0.78539815* (4.75*x(1)*x (2)A2*x(3)A2+85*x(1)*x (2) *x(3)A2-85*x(1) *x(3)A2+0.92*x(1)*x(6)A2-x(1)*x(5)A2+0.8*x(1)*x (2)*x (3)*x (6)-1. 6*x(1)*x(3)*x(6)+x(4)*x(5F2+x(4)*x(6)A2+28*x(5)A2+32*x(6)A2)然后建立约束条件程序confun 1.mfun ctio nc,ceq=con strai nt(x)c(1)=x( 2)*x (3)-33c(2)=16-x(1)/x(3)c(3)=x(1)/x(3)-35c(4)=44163/(x (2) *x (3) *sqrt(x(1)-855c(5)=-261+7098/(x(1)*x (2)*x(3)A2*(0.169+0.6666/100*x(2)-0.854/1 0000*x (2) A2)c(6)=-213+7098/(x(1)*x (2)*x(3)A2*(0.2824+0.177/100*x(2)-0.314/1 0000*x (2) A2)c=-0.03*x (4)+0.01233*x (4) A3/x(2)/x(3)/(x (5)人4)c(8)=-55+1/(x (5) A3)*sqrt(29050*x(4)/x(3)/x(2)A2+(0.58*27300)A2)c(9)=-55+1/(x (6)A3)*sqrt(29050*x(4)/x(3)/x(2)A2+(0.58*5*27300)A2)ceq=;在命令窗口键入：x(0)=13,21,0.8,42,12,16;A=；b=；LBn d=;constraint )x,f=fmincon( myfun ,x(0),A,b,LBnd, 优化结果为：2、使用复合型法进行优化：在命令窗口键入：f=0.78539815* (4.75*x(1)*x (2)A2*x(3)A2+85*x(1)*x (2) *x(3)2-85*x(1)*x(3)A2+0.92*x(1)*x(6F2-x(1)*x(5)A2+0.8*x(1)*x (2)*x (3)*x (6)-1. 6*x(1)*x(3)*x(6)+x(4)*x(5)A2+x(4)*x(6)A2+28*x(5)A2+32*x(6)A2);g=x(2)*x(3)-33,16-x(1)/x(3),x(1)/x(3)-35,44163/(x(2)*x(3)*sqrt (x(1)-855,-261+7098/(x(1)*x(2)*x(3)A2*(0.169+0.6666/100*x(2)-0. 854/10000*x(2)A2),-213+7098/(x(1)*x(2)*x(3)A2*(0.2824+0.177/100* x(2)-0.314/10000*x(2)A2),-0.03*x(4)+0.01233*x(4)A3/x(2)/x(3)/(x( 5)A4),-55+1/(x(5)A3)*sqrt(29050*x(4)/x(3)/x(2)A2+(0.58*27300)A2 ),-55+1/(x(6)A3)*sqrt(29050*x(4)/x(3)/x(2)A2+(0.58*5*27300)A2), 0.2-x(3),10-x(5),13-x(6),x(5)-15,x(6)-20,4-x(4)+x(1)+0.5*x(6);X=;x,minf = minconSimpSearch(f,g,x,1.2,0.5,8,0.3);所得结果为：3、使用内点罚函数法进行优化：