平面连杆机构优化设计

1、平面连杆机构优化设计一、问题描述平面连杆机构是由所有构件均由低副连接而成的机构，四杆机构是最常用的平面连杆机构。一般情况下，四杆机构只能近似实现给定的运动规律或运动轨迹，精确设计较为复杂。在四杆机构中，若两连架杆中的一个是曲柄，另一个是摇杆，则该机构为曲柄摇杆机构。曲柄摇杆机构可将曲柄的连续转动转变为摇杆的往复摆动。设计一曲柄摇杆机构（如图1所示）。已知曲柄长度li=100mm机架长度14=500mm摇杆处于右极限位置时，曲柄与机架的夹角为（1）0,摇杆与机架的夹角为巾0。在曲柄转角小从口匀速增至小0+90°的过程中，要求摇一一2杆转角中=中。+已（9_中。2。为防止从动件卡死，连杆

2、与摇杆的夹角丫只3冗允许在45°135°范围内变化。图1机构运动简图、基本思路四杆机构的设计要求可归纳为三类，即满足预定的连杆位置要求、满足预定的运动规律要求、满足预定的轨迹要求。本案例中，要求曲柄作等2速转动时，摇杆的转角酒足预定运动规律中E=4%+2（邛中0）。优化设计3冗时，通常无精确解，一般采用数值方法得到近似解。本案例将机构预定的运动规律与实际运动规律观测量之间的偏差最小设为目标，由此建立优化设计数学模型，并运用MATLA就化工具箱的相关函数进行求解。三、要点分析优化设计数学模型的三要素包括设计变量、目标函数和约束条件。依次确定三要素后，编写程序进行计算。1,设计

3、变量的确定通常将机构中的各杆长度，以及摇杆按预定运动规律运动时，曲柄所处的初始位置角？列为设计变量，即X-(xiX2X3X4X5)=。1l2l3l4.0)考虑到机构各杆长按比例变化时，不会改变其运动规律，因此在计算可取11为单位长度，而其他杆长则按比例取为11的倍数。若曲柄的初始位置对应摇杆的右极限位置，则小0及巾。均为杆长的函数，即(|)2|2一|2o=arccos(l1l2)l4l32。1。4=arccos(li1)2-J32m因此，设计变量缩减为3个独立变量，即X=(xiX2X3)=(12l3l4)(4)2.目标函数的建立以机构预定的运动规律观测量巾Ei与实际运动规律观测量巾i之间的偏差

4、平方和最小为指标来建立目标函数，即mf(X)='('-Ei-1-i)2，min(5)iJ式中，m为输入角的等分数；巾E为预期输出角,巾Ei=巾E（小i）；巾i为实际输出角。由图2可知:冗一叫一Pi（0<%<力Tt->i'（九一i：二24=arccos：i2l3-lf2：h_,?222：i二arccos4-(8)i2,)1=.l；I"2IiLcos；(9)(a)00小i<兀(b)九0di<2几图2曲柄摇杆机构的运动学关系3.约束条件的确定(1)曲柄摇杆机构应满足曲柄存在条件，可得gi(x)=ii-%工0(10)g2(X)=ii-13

5、三0(ii)g3(X)=ii-14三0(i2)g4(X)=iiI4-I2-13-0(i3)gs(X)=iiI2-I3-I4-0(i4)ge(X)=iii3-i2-i4三0(i5)(2)连杆与摇杆的夹角应在Ymin和Ymax之|可，即1g7(X)=arccosijiZ)2_23max-0(16)g8(X)=minl|l1-(l4-li)2八-arccos0想3(17)四、具体步骤1 .选择设计变量已知li=100mm14=500mm且小0和巾0不是独立参数，它们可由下式(2)、式(3)求出，即(10012)2250000-1；:n'arccos(10012)2-250000-1；0=ar

6、ccos2310003所以该问题只有两个独立参数12和13,故设计向量为X=(XiX2)T=(1213)T2 .建立目标函数330,将输入角分成30等分，并依次取30个观测点巾1,巾2,.,得目标函数30f(X)八L日i)2i=1式中：i=Tt-C(i-Pi222_riX2-X12底222二iri13-12:arccos2rh:arccos/J1；=arccos/2400002r141000riri-,.j112142-21114cos；=.260000-100000cosi'-Ei-'-0AL023冗3 .确定约束条件约束函数按曲柄存在条件及对传动角的限制来建立，得g(X)=

7、100<0g2(X)=100-x2M0g3(X)=600-Xi-X2M0g4(X)-Xi-X2-400<0g5(X)=X2f-400三022g6(X)=XiX2-1.414XiX2-160000<022g7(X)=360000-x；-x|-1.414X1X2_04.MATLAB程序及优化结果这是一个具有2个设计变量、7个不等式约束条件的优化设计问题。应用MATLAB件的优化工具箱的fmincon函数对上述优化问题求解。(1)编写m文件Objfun.m定义目标函数。functionf=objfun(X)11=100；14=500;th0=acos(100+X(1)A2-X(2)

8、A2+250000)/(1000*(100+X(1);ps0=acos(100+X(1)A2-X(2)A2-250000)/(1000*X(2);f=0;forth=th0:pi/2/30:th0+pi/2r=(10000+250000-2*100*500*cos(th)A0.5a=acos(rA2+x(2)A2-x(1)A2)/(2*r*x(2)b=acos(rA2+240000)/(1000*r)；ps=pi-a-b；pse=ps0+2/(3*pi)*(th-th0)A2；f=f+(ps-pse)A2；end编写加件confun.m定义约束。functionc,ceq=confun(x)c

9、(1)=100-x(1)；c(2)=100-x(2);c(3)=600-x(1)-x(2);c(4)=x(1)-x(2)-400;c(5)=x(2)-x(1)-400;c(6)=x(1)A2+x(2)A2-1.414*x(1)*x(2)-160000c(7)=360000-x(1)A2-x(2)A2-1.414*x(1)*x(2)ceq=;编写mmC件run.m求解计算。x0=400400;options=optimset('LargeScale','off)；x,fval=fmincon(objfun,x0,口，口，口，口，口，口，confun)(4)运行m文件run.m,得最优解X*=(412.8926mm,232.2417mm),f(X*)=0.0076mm2。五、问题拓展满足