matlab(四连杆优化设计)

1、机械优化设计在matlab中的应用东南大学机械工程学院*一优化设计目的：在生活和工作中，人们对于同一个问题往往会提出多个解决方案，并通过各方面的论证从中提取最佳方案。最优化方法就是专门研究如何从多个方案中科学合理地提取出最佳方案的科学。由于优化问题无所不在，目前最优化方法的应用和研究已经深入到了生产和科研的各个领域，如土木工程、机械工程、化学工程、运输调度、生产控制、经济规划、经济管理等，并取得了显著的经济效益和社会效益。二优化设计步骤：机械优化设计的全过程一般可以分为如下几个步骤：1）建立优化设计的数学模型；2）选择适当的优化方法；3）编写计算机程序；4）准备必要的初始数据并伤及计算；5）对

2、计算机求得的结果进行必要的分析。其中建立优化设计数学模型是首要的和关键的一步，它是取得正确结果的前提。优化方法的选取取决于数学模型的特点，例如优化问题规模的大小，目标函数和约束函数的性态以及计算精度等。在比较各种可供选用的优化方法时，需要考虑的一个重要因素是计算机执行这些程序所花费的时间和费用，也即计算效率。建立数学模型的基本原则与步骤设计变量的确定；设计变量是指在优化设计的过程中，不断进行修改，调整，一直处于变化的参数称为设计变量。设计变量的全体实际上是一组变量，可用一个列向量表示：x二，.，-.:.II。目标函数的建立；选择目标函数是整个优化设计过程中最重要的决策之一。当对某以设计性能有特

3、定的要求，而这个要求有很难满足时，则针对这一性能进行优化会得到满意的效果。目标函数是设计变量的函数，是一项设计所追求的指标的数学反映，因此它能够用来评价设计的优劣。目标函数的一般表达式为：f（x）=，要根据实际的设计要求来设计目标函数。约束条件的确定。一个可行性设计必须满足某些设计限制条件，这些限制条件称为约束条件，简称约束。由若干个约束条件构成目标函数的可行域，而可行域内的所有设计点都是满足设计要求的，一般情况下，其设计可行域可表示为在可行域中，任意设计点满足全部约束条件，称为可行解，但不是最优解，而优化设计就是要求出目标函数在可行域的最优解。三实例分析（机械优化设计P241页例8-5）设计

4、一曲柄摇杆机构如图，要求：曲柄.从且已范围内变化。分析：1）设计变量的确定决定机构尺寸的各杆长度，以及当摇杆按已知运动规律开始运行时，曲柄所载的位置角应列为设计变量，即：X=考虑到机构的杆长按比例变化时，不会改变其运动规律，因此在计算时常取，而其他杆长则按比例取为的倍数。若取曲柄的初始位置角为极位角，则及相应摇杆位置角均为杆长的函数，几何图形关系如右图，其关系式为：TOC o 1-5 h z=arcos（1）=arcos（2）将的长度代入上式（1）,（2）得到：=arcos=arcos因此，只有为独立变量，设计变量减少，故最后的设计变量为：:；7=目标函数的建立目标函数可根据已知的运动规律与机

5、构实际运动规律之间的偏差最小为指标来建立，即f(x)=min式中一期望输出角，m输入角等分数；实际输出角，由下图得:a)0W式中b)(3)将输入角分成30等分，并用近似公式计算，可得目标函数的表达式：f(x)=由题意知，传动角的变化范围是，则上式中变量的最后形式可以写成:arcos将带入（4）（5）得:为当时的理想输出角，其值在题目中已经给出:3）约束条件的确定曲柄摇杆机构应满足曲柄存在条件，可得曲柄摇杆机构的传动角可得把约束条件简化（=5x=g7(x)二云十云1.414阳16youhua4.15742.2909%最优解fval=5.1899e-004%目标函数最优点的值exitflag=%标

6、志值，5表示重要方向导数小于规定的容许范围并且约束违背小于options.ToLConoutput=iterations:12%迭代次数funcCount:40%函数的评价次数lssteplength:1stepsize:7.6955e-005algorithm:medium-scale:SQP,Quasi-Newton,line-search%采用的中型算法firstorderopt:1.0832e-006constrviolation:-1.0852e-006%一阶最优性条件message:1x780char%跳出信息lambda=lower:2x1doubleupper:2x1doubl

7、eeqlin:0 x1doubleeqnonlin:0 x1doubleineqlin:5x1doubleineqnonlin:2x1doublegrad=1.0e-003*0.48880.4445%函数在最优点处梯度信息hessian=%函数在最优点处海塞矩阵0.00160.00750.00750.0468结果分析采用fmincon求解的最优值:：.=4.1574;2.2909;采用算法：中型算法(mediun-scale)。这与课本给出的最优解：=4.1286;2.3325,=0.0156相比，计算精度更高，最优解的数值更精确，故计算准确度高。用matlab绘制输入输出曲线关系图上图中(单

8、位为“度”蓝色的线代表曲柄摇杆机构的实际输出角与输入角的关系，红色的线代表理想输出角与输入角的关系。可以看出：实际输出和理论输出曲线之间存在线性误差，其最大线性误差为，误差在允许的范围之内，故结果的可信度也较大，运用matlab优化工具箱计算所得结果正确。小结通过结合实际问题的分析，计算，求解，更加深入地了解和掌握机械优化设计的过程和步骤，比较重要的步骤是数学模型的建立，以及设计变量的选取，以及数学模型的尺度变换，根据机构实际工作需要，建立目标函数的约束条件等等，当数学模型建好以后，剩下的工作可以再matlab里面完成，而matlab里面的优化工具箱，给用户提供了多种优化函数，使用者只需要将数

9、学模型按要求编写成子程序嵌入已有的优化程序即可。在设计过程中也遇到一些困难，比如说在在用matlab计算时，计算机已知处于busy状态，得不到函数的最优解，最后反复的检查，终于找的了其原因，是由于初始点选择不恰当引起的，如果初始点选择得好，可以节省计算时间和计算空间，故初始点的选取比较重要。附录1.编写目标函数M文件myfun.m：functionf=myfun(x)f=0;aO=acos(1+x(1)入2-x(2)%函数f赋初值入2+25)/(10*(1+x(1);%初始计算点曲柄和摇杆的角度bO=acos(1+x(1)入2-x(2)i=2;入225)/(10*x(2);while(i=31

10、)a(i)=a0+(pi/2)*(i/30);b(i)=b0+2*(a(i)a0)入2/(3*pi);r=sqrt(2610*cos(a(i);c(i)=acos(r2+x(2)A2x(1)人2)/(2*x(2)*r);d(i)=acos(r入2+24)/(10*r);ifa(i)=pie(i)=pic(i)d(i);elseifa(i)=2*pie(i)=pic(i)+d(i);end%设置迭代次数为30次%计算曲柄各分度的角度值%计算摇杆各分度的角度值%计算摇杆输出的实际值enda(1)=a0;f=f+(b(i)e(i)入2)*(a(i)a(i1);i=i+1;%目标函数的计算end2.编写非线性不等式约束M文件constrain.m：functioncceq=constrain(x)c=36x(1)八2x(2)八21.414*x(1)*x(2);x(1)入2+x(2)入2-1.414*x(1)*x(2)-16;ceq=;3.调用fmincon优化函数，建立youhua.m文件：lb=1;1;%非线性不等式约束x0=4;2;A=-1,0;0,-1;-1,-1;1,-1;-1,1;b=-1;-1;-6;4;4;%设计变量的下界%迭代初始点%线性不等