基于MATLAB的优化设计

基于MATLAB的曲柄摇杆机构优化设计1.问题的提出根据机械的用途和性能要求的不同，对连杆机构设计的要求是多种多样的，但这些设计要求可归纳为以下三种问题：（1）满足预定的运动规律要求；（2）满足预定的连杆位置要求；（3）满足预定的轨迹要求。在在第一个问题里按照期望函数设计的思想，要求曲柄摇杆机构的曲柄与摇杆转角之间按照f（称为期望函数）的关系实现运动，由于机构的待定参数较少，故一般不能准确实现该期望函数，设实际的函数为F（称为再现函数），而再现函数一般是与期望函数不一致的，因此在设计时应使机构再现函数F尽可能逼近所要求的期望函数f。这时需按机械优化设计方法来设计曲柄连杆，建立优化数学模型，研究并提出其优化求解算法，并应用于优化模型的求解，求解得到更优的设计参数。2.曲柄摇杆机构的设计在图1所示的曲柄摇杆机构中，1l、2l、3l、4l分别是曲柄AB、连杆BC、摇杆CD和机架AD的长度。这里规定0为摇杆在右极限位置0时的曲柄起始位置角，它们由1l、2l、3l和4l确定。图1曲柄摇杆机构简图设计时，可在给定最大和最小传动角的前提下，当曲柄从0转到090时，要求摇杆的输出角最优地实现一个给定的运动规律f。这里假设要求：20023Ef（1）对于这样的设计问题，可以取机构的期望输出角Ef和实际输出角F的平方误差之和作为目标函数，使得它的值达到最小。2.1设计变量的确定决定机构尺寸的各杆长度1l、2l、3l和4l，以及当摇杆按已知运动规律开始运行时，曲柄所处的位置角0应列为设计变量，即:12340Txllll(2)考虑到机构的杆长按比例变化时，不会改变其运动规律，通常设定曲柄长度1l=1.0，在这里可给定4l=5.0，其他杆长则按比例取为1l的倍数。若取曲柄的初始位置角为极位角，则及相应的摇杆l位置角均为杆长的函数，其关系式为:222221243230124225arccos210llllllllllll（3）22222124323034325arccos210llllllllll（4）因此，只有2l、3l为独立变量，则设计变量为2312TTxllxx。2.2目标函数的建立目标函数可根据已知的运动规律与机构实际运动规律之间的偏差最小为指标来建立，即：21minmEiiifx（5）式中，Ei-期望输出角；m-输出角的等分数；i-实际输出角，由图1可知：02iiiiiii(6)式中，222222322132arccosarccos22iiiiirllrxxrlrx(7)222241424arccosarccos210iiiiirllrrlr(8)2214142cos2610cosiiirllll(9)2.3约束条件曲柄存在条件:12131423;,llllllll24133412,llllllll曲柄与机架共线位置时的传动角（连杆BC和摇杆CD之间的夹角）:最小传动角minmin45rBCD最大传动角maxmax135rBCD由上面的分析可以算出：222222234112min231216arccos4522llllxxrllxx（10）222222234112max231236arccos13522llllxxrllxx（11）3.用MATLAB工具箱优化计算结果通过上面的分析后，将输入角分成30等分（m=30），经过转化为标准形式得到曲柄摇杆机构优化设计标准数学模型为:21minmEiiifx2312TTxllxx112231241252122612122271212101060..40401.414360361.4140gxxgxxgxxxstgxxxgxxxgxxxxxgxxxxx(12)机械优化设计中的问题，大多数属于约束优化问题，此为非线性约束优化问题，运用MATLAB优化工具箱的命令函数fmincon来处理有约束的非线性多元函数最小化优化问题。3.1编写程序求解（1）首先编写目标函数M文件optimfun.mfunctionf=optimfun（x）；s=30；qb=1；jj=5；fx=0；fa0=acos（（（qb+x（1））^2-x（2）^2+jj^2）/（2*（qb+x（1））*jj））；%曲柄初始角pu0=acos（（（qb+x（1））^2-x（2）^2-jj^2）/（2*x（2）*jj））；%摇杆初始角fori=1:sfai=fa0+0.5*pi*i/s；pui=pu0+2*（fai-fa0）^2?（3*pi）；ri=sqrt（qb^2+jj^2-2*qb*jj*cos（fai））；alfi=acos（（ri^2+x（2）^2-x（1）^2）/（2*ri*x（2）））；bati=acos（（ri^2+jj^2-qb^2）（/2*ri*jj））；iffai0&fai=pipsi=pi-alfi-bati;elseiffaipi&fai=2*pipsi=pi-alfi+bati；endfx=fx+（pui-psi）^2；endf=fx；（2）编写非线性约束函数M文件confun.mfunction［c，ceq］=confun（x）；qb=1；jj=5；m=45*pi/180；n=135*pi/180；c（1）=x（1）^2+x（2）^2-（jj-qb）^2-2*x（1）*x（2）*cos（m）；%最小传动角约束c（2）=-x（1）^2-x（2）^2+（jj+qb）^2+2*x（1）*x（2）*cos（n）；%最大传动角约束ceq=［］；（3）在MATLAB命令窗口调用优化程序x0=［6；4］；lb=［1；1］；ub=［］；%线性不等式约束a=［-1-1；1-1；-11］；b=［-6；4；4］；［x，fn］=fmincon（@optimfun，x0，a，b，［］，［］，lb，ub，@confun）；（4）运行结果5.结论MATLAB优化工具箱具有强大的优化工具，应用它求解优化问题时工作量小，操作简单，计算结果精确，大大地提高了设计的时效性和准确性。利用MATLAB优化工具箱对曲柄摇杆机构设计，达到了设计的预期目的。参考文献［1］孙桓，陈作模，葛文杰.机械原理［M］.北京:高等教育出版社，2006（5）：125-126.［2］龚水明，詹小刚.基于MATLAB工具箱的机械优化设计［J］.机械工程师，2008（10）：10-11.［3］孙靖民.现代机械设计方法［M］.哈尔滨：哈尔滨工业大学出版社，2003.［4］褚洪生，杜增吉，阎金华，等.MATLAB72优化设计实例指导教程［M］.北京：机械工业出版社，2006.