matlab四连杆带代码

用Matlab对四连杆运动模拟引言四连杆机构因其结构敏捷、能够传递动力并有效地实现预定动作，在许多领域得到了广泛应用。进行连杆机构运动分析，传统方法主要是图解法或分析法，无论设计精度还是设计效率都相对低下，无法满意现代机械高速高精度的 要求 。随着计算机技术的飞速发展 ，特殊是以MATLAB为代表的数值计算软件的出现， 为进行机构分析供应了有力的工具。1、四连杆介绍1.1、四连杆介绍及分类全部运动副均为转动副的四杆机构称为铰链四杆机构，它是平面四杆机构的基本形式，其他四杆机构都可以看成是在它的基础上演化而来的。选定其中一个构件作为机架之后，干脆及机架链接的构件称为连架杆，不干脆及机架连接的构件称为连杆，能够做整周回转的构件被称作曲柄，只能在某一角度范围内往复摇摆的构件称为摇杆。假如以转动副连接的两个构件可以做整周相对转动，则称之为整转副，反之称之为摆转副。铰链四杆机构中，根据连架杆是否可以做整圆周转动，可以将其分为三种基本形式，即曲柄摇杆机构，双曲柄机构和双摇杆机构。 曲柄摇杆机构，两连架杆中一个为曲柄一个为摇杆的铰链四杆机构。 双曲柄机构，具有两个曲柄的铰链四杆机构称为双曲柄机构。其特点是当主动曲柄连续等速转动时，从动曲柄一般做不等速转动。在双曲柄机构中，假如两对边构件长度相等且平行，则成为平行四边形机构。这种机构的传动特点是主动曲柄和从动曲柄均以相同的角速度转动，而连杆做平动。 双摇杆机构。双摇杆机构是两连架杆均为摇杆的铰链四杆机构。1.2、格拉霍夫定理杆长之和条件：平面四杆机构的最短杆和最长杆的长度之和小于或者等于其余两杆长度之和。在铰链四杆机构中，假如某个转动副能够成为整转副，则它所连接的两个构件中，必有一个为最短杆，并且四个构件的长度关系满意杆长之和条件。在有整装副存在的铰链四杆机构中，最短杆两端的转动副均为整转副。此时，假如取最短杆为机架，则得到双曲柄机构；若取最短杆的任何一个相连构件为机架，则得到曲柄摇杆机构；假如取最短杆对面构件为机架，则得到双摇杆机构。假如四杆机构不满意杆长之和条件，则不论选取哪个构件为机架，所得到机构均为双摇杆机构。上述系列结论称为格拉霍夫定理。 运用条件分析： 1.Lmax+Lmin其余两杆之和,-此铰链四杆机构为双摇杆机构; 2.Lmax+Lmin其余两杆之和,要详细分析: (1)Lmin为机架-为双曲柄机构 (2)Lmin为连架杆-为曲柄摇杆机构(Lmin为曲柄) (3)Lmin为连杆-为双摇杆机构1.3、理论分析对图1的四连杆机构列出方程组。其中L1、L2、L3、L4为四根杆，L4为机架。 L1cos1+L2cos2=L3cos1+L4 L1sin1+L2sin2=L3sin3 由此两个方程消去2，便可得到一个关于1、3的函数F(1, 3)=（L3cos3+L4-L1cos1）2+（L3sin3-L1sin1）2-L22然后用隐函数求出 3=3(1)另 解出角速度：1=w1tw2=d(3)/d(t)2、Matlab的实现(一) 数据分析给出了3组数据：a) L1=8;L2=25;L3=20;L4=20; b) L1=7;L2=9;L3=11;L4=4; c) L1=7;L2=6.2;L3=10;L4=6; 用格拉霍夫定理进行分析得到：a为曲柄摇杆机构;b为双曲柄机构；c双摇杆机构；(二) Matlab求解并进行动画模拟a) L1=8;L2=25;L3=20;L4=20; （曲柄摇杆机构）3=3(1) 图像实时动画模拟（部分截屏）W2的图像说明：其中w2(1)为先符号求导再赋值计算，w2(2)为干脆数值求导。（下同）b) L1=7;L2=9;L3=11;L4=4; （双曲柄机构） 3=3(1) 图像实时动画模拟（部分截屏）W2的图像c) L1=7;L2=6.2;L3=10;L4=6; （双摇杆机构）3=3(1) 图像实时动画模拟（部分截屏）W2的图像3、结果分析对3的角速度w2的求值得两种方法中，分析图像知：符号求导再代入数值更简单出现大的偏差（w2的图像中的红线）。干脆用数值求导得到的解更平稳。4、程序1、函数文件 dd.mfunction F=dd(y)global L1 L2 L3 L4 p x=p;F=(L3*cos(y)+L4-L1*cos(x)2)+(L3*sin(y)-L1*sin(x)2)-(L22);2、计算3的文件 y_figure.mglobal L1 L2 L3 L4 p %传递参数% L1=8;L2=25;L3=20;L4=20; %曲柄摇杆机构% L1=7;L2=9;L3=11;L4=4; %双曲柄机构L1=7;L2=6.2;L3=10;L4=6; %双摇杆机构% L1=9;L2=5;L3=10;L4=7; %双摇杆机构 Ltotal=L1+L2+L3+L4; Lmax=max(L1 L2 L3 L4); Lmin,I=min(L1 L2 L3 L4); m=100; b=0;while (Lmax+Lmin)=(Ltotal/2) str=双摇杆机构; b=1; breakendwhile (Lmax+Lmin)m; k=1; end xb=L1*cos(x(k)+xa; yb=L1*sin(x(k)+ya; xc=L3*cos(y(k)+L4;yc=L3*sin(y(k); set(stick_1,xdata,xa;xb,ydata,ya;yb); set(stick_2,xdata,xc;xb,ydata,yc;yb); set(stick_3,xdata,xd;xc,ydata,yd;yc); set(dot_b,xdata,xb,ydata,yb); set(dot_c,xdata,xc,ydata,yc); set(gcf,doublebuffer,on); drawnow;pause k=k+1;end