还剩4页未读,继续阅读
文本内容:
用MatIab对四连杆运动模拟引言四连杆机构因其结构灵活、能够传递动力并有效地实现预定动作,在很多领域得到了广泛应用进行连杆机构运动分析,传统方法主要是图解法或分析法,无论设计精度还是设计效率都相对低下,无法满足现代机械高速高精度的要求随着计算机技术的飞速发展,特别是以MATLAB为代表的数值计算软件的出现,为进行机构分析提供了有力的工具
1、四连杆介绍四连杆介绍与分类所有运动副均为转动副的四杆机构称为钱链四杆机构,它是平面四杆机构的基本形式,其他四杆机构都可以看成是在它的基础上演化而来的选定其中一个构件作为机架之后,直接与机架链接的构件称为连架杆,不直接与机架连接的构件称为连杆,能够做整周回转的构件被称作曲柄,只能在某一角度范围内往复摆动的构件称为摇杆如果以转动副连接的两个构件可以做整周相对转动,则称之为整转副,反之称之为摆转副钱链四杆机构中,按照连架杆是否可以做整圆周转动,可以将其分为三种基本形式,即曲柄摇杆机构,双曲柄机构和双摇杆机构曲柄摇杆机构,两连架杆中一个为曲柄一个为摇杆的钱链四杆机构双曲柄机构,具有两个曲柄的钱链四杆机构称为双曲柄机构其特点是当主动曲柄连续等速转动时,从动曲柄一般做不等速转动在双曲柄机构中,如果两对边构件长度相等且平行,则成为平行四边形机构这种机构的传动特点是主动曲柄和从动曲柄均以相同的角速度转动,而连杆做平动双摇杆机构双摇杆机构是两连架杆均为摇杆的钱链四杆机构格拉霍夫定理杆长之和条件平面四杆机构的最短杆和最长杆的长度之和小于或者等于其余两杆长度之和在校链四杆机构中,如果某个转动副能够成为整转副,则它所连接的两个构件中,必有一个为最短杆,并且四个构件的长度关系满足杆长之和条件在有整装副存在的钱链四杆机构中,最短杆两端的转动副均为整转副此时,如果取最短杆为机架,则得到双曲柄机构;若取最短杆的任何一个相连构件为机架,则得到曲柄摇杆机构;如果取最短杆对面构件为机架,则得到双摇杆机构如果四杆机构不满足杆长之和条件,则不论选取哪个构件为机架,所得到机构均为双摇杆机构上述系列结论称为格拉霍夫定理运用条件分析l.Lmax+Lmin>其余两杆之和此钱链四杆机构为双摇杆机构;
2.Lmax+LminW其余两杆之和,要具体分析lLmin为机架----为双曲柄机构2Lmin为连架杆为曲柄摇杆机构Lmin为曲柄3Lmin为连杆.一-为双摇杆机构
1.
3、理论分析对图1的四连杆机构列出方程组其中LI、L
2、L
3、L4为四根杆,L4为机架LIxcos%4-L2xcos02=L3xcos%+L4LIxsin%4-L2xsin%=i3sin03由此两个方程消去2便可得到一个关于
1、3的函数F4032o=L3xcos03+L4-LIxcosOr4-L3xsin03-LIxsin0AJ-L22然后用隐函数求出另解出角速度:%=wlXtw2=d
(03)/d(t)
2、Matlab的实现(-)数据分析给出了3组数据Ll=8;L2=25;L3=20;L4=20;L1=7;L2=9;L3=11;L4=4;Ll=7;L2=
6.2;L3=10;L4=6;用格拉霍夫定理进行分析得到a为曲柄摇杆机构;b为双曲柄机构;c双摇杆机构;
(二)Matlab求解并进行动画模拟Ll=8;L2=25;L3=20;L4=20;(曲柄摇杆机构)3=3
(1)图像实时动画模拟(部分截屏)W2的图像说明其中w2(l)为先符号求导再赋值计算,w2
(2)为直接数值求导(下同)L1=7;L2=9;L3=11;L4=4;(双曲柄机构)3=%(%)图像实时动画模拟(部分截屏)W2的图像Ll=7;L2=
6.2;L3=10;L4=6;(双摇杆机构)3=3(%)图像实时动画模拟(部分截屏)W2的图像
3、结果分析对%的角速度w2的求值得两种方法中,分析图像知符号求导再代入数值更容易出现大的偏差w2的图像中的红线直接用数值求导得到的解更平稳
4、程序
1、函数文件dd.mfunctionF=ddyglobalLIL2L3L4px=p;F=L3*cosy+L4-Ll*cosxA2+L3*siny-Ll*sinxA2-L2A2;
2、计算出的文件yjigure.mglobalLIL2L3L4p%传递参数%Ll=8;L2=25;L3=20;L4=20;%曲柄摇杆机构%L1=7;L2=9;L3=11;L4=4;%双曲柄机构Ll=7;L2=
6.2;L3=10;L4=6;%双摇杆机构%Ll=9;L2=5;L3=10;L4=7;%双摇杆机构Ltotal=Ll+L2+L3+L4;Lmax=max[LlL2L3L4];[LminJ]=min[LlL2L3L4];m=100;b=0;whileLmax+Lmin=Ltotal/2str±双摇杆机构,;b=l;breakendwhileLmax+LminLtotal/2ifl==4sg双曲柄机构,;elseifl==l111==3str=曲柄摇杆机构;elseifl==2strd双摇杆机构,;b=l;endendendbreakendswitchbcase0m=100;x=linspace©2*pim;%确定隐函数自变量的范围y0=
0.7341;%第一个方程的初值case1a=
0.58;m=80;x=linspacea2*pi-am;y0=
1.2465;endy=[];f=H;fork=l:mP=xk;[ylJvaLexitflagoutput]=fzeroddyO;yO=yl;endfigureplotxl:myl:m匕%绘制隐函数图形gridon
3、计算w2的文件w_figure.mglobalLIL2L3L4wlwl=100;symsxxyy%xx为输入杆转角,yy为输出角转角F=Ll*cosxx+L2*l-L3*sinyy-Ll*sinxx/L2A2Al/2-L3*cosyy-L4;w=-l*diffFxx/diffFzyy*wl;w21;subsw{xxyy}{xy};%用隐函数求导t=x/wl;w22=diffy./difft;%直接数值求导plotw21/rholdonplotw22/k
4、实时动画文件siliangan.mpic=figurenamestr;setfpic/color/white;axisequalgridonaxisd-1030-1020];axisonglobalLIL3L4xa=0;ya=0;%A点坐标xd=L4;yd=0;%D点坐标xb=Ll*cosxl;yb=Ll*sinxl;%B点坐标xc=L3*cosyl+L4;yc=L3*siny⑴;%C点坐标stick_l=line[xa;xb]/[ya;yb]/color7red7linewidthz5/linestyle7-;stick_2=line[xb;xc][yb;yc]/color/g/linewidth/5/linestyle7-;stick_3=line[xc;xd]/[yc;yd]zcolor/b7linewidth/5zlinestyle7-;stick_0=line[xd;xa][yd;ya]/color7ylinewidth/5/linestyle7-;dot_a=linexaya/color7r7linestyle
7./markersize30;dot_b=linexbzyb/colorblackzlinestyle
7./,nnarkersize730;dot_c=linexcyccolor「「color[
0.
10.
70.3Llinestyle「.「markersize30;dot_d=linexdyd;colorr;linestylez.;markersize/30;dt=2*pi/m;k=l;while1ifkm;k=l;endxb=Ll*cosxk+xa;yb=Ll*sinxk+ya;xc=L3*cosyk+L4;yc=L3*sinyk;setstick_lxdata[xa;xb]ydata[ya;yb];setstick_2/xdata/[xc;xb]ydata/[yc;yb];setfstick-S/xdata[xdjxcl/ydata[yd;yc];setfdot-b/xdata\xb/ydata^yb;setdot_c「xdataxcydatayc;setfgcf/doublebuffer/on;drawnow;pausek=k+l;end。