Matlab在四杆机构运动分析中的应用.doc

Matlab在四杆机构运动分析中的应用MATLAB软件由美国MathWorks公司于1982年推出，经过十几年的发展和竞争，现已成为国际公认的最优秀的科技应用软件之一。MATLAB提供了强大的矩阵处理和绘图功能。它主要包括两部分内容：核心函数和工具箱。 Matlab编程代码接近数学推导公式，简洁直观，与科技人员的思维方式和书写习惯相适应，操作简易，人机交互性能好，且可以方便迅速地用三维图形、图像、声音、动画等表达计算结果，拓展思路。编制相应的M函数文件。Pos.m用于求解位置、角度和角加速度。th1为初始角度，w1为杆1角速度，其余为杆长。将课本P35（i）、（ii）、（iii）式用MATLAB语言表述，即可编制为四杆机构求解函数文件pos.m。function f=pos(th1,w1,l1,l2,l3,l4)syms x21 x31 x22 x32x1=th1*pi/180;x11=cos(x1);x12=sin(x1);eq1=l1*x11+l2*x21-l3*x31-l4;eq2=l1*x12+l2*x22-l3*x32;eq3=x212+x222-1;eq4=x312+x322-1;s=solve(eq1,eq2,eq3,eq4,x21,x22,x31,x32);s1=double(s.x21); s2=double(s.x22);s3=double(s.x31); s4=double(s.x32);x2=(acos(s1(1,1)/pi*180;x3=(acos(s3(1,1)/pi*180;A=l2*s2(1,1),-l3*s4(1,1);l2*s1(1,1),-l3*s3(1,1);B=-l1*x12;-l1*x11;w=A(w1*B);w2=w(1,1);w3=w(2,1);C=-l2*w(1,1)*s1(1,1),l3*w(2,1)*s3(1,1);l2*w(1,1)*s2(1,1),-l3*w(2,1)*s4(1,1);D=w(1,1);w(2,1);E=-l1*w1*x11;l1*w1*x12;F=l2*s2(1,1),-l3*s4(1,1);l2*s1(1,1),-l3*s3(1,1);t=F(C*D+w1*E);a2=t(1,1);a3=t(2,1);l1=cat(1,th1,w1,0);l2=cat(1,x2,w2,a2);l3=cat(1,x3,w3,a3);f=(cat(2,l1,l2,l3);subplot(2,2,1);plot(th1,w2,r-);hold on;title(连杆2角速度分析);subplot(2,2,2);plot(th1,a2,b-);hold on;title(连杆2角加速度分析);subplot(2,2,3);plot(th1,w3,r-);hold on;title(连杆3角速度分析);subplot(2,2,4);plot(th1,a3,b-);hold on;title(连杆3角加速度分析);return在MA TLAB命令窗口输入命令： th1=60;w1=pi/3;l1=20;l2=30;l3=40;l4=45; function f=pos(th1,w1,l1,l2,l3,l4)即可得到 theta=60，=pi/3时的值采用实时动画的方法编写draw.m文件：function dw = draw(l1,l2,l3,l4,x1,x2,x3)figure(2);th1=x1*pi/180;th2=x2*pi/180;th3=x3*pi/180;x=0:0.001:l4;plot(x,0,r-);axis(-25,70,-25,60);hold on for i=0:0.1:l1; s=i*cos(th1); c=i*sin(th1);plot(s,c,b-);hold on;end;for x=0:0.1:l3;s2=x*cos(th3)+l4;c2=x*sin(th3);plot(s2,c2);hold on;endfor x=0:0.1:l2; s3=l1*cos(th1)+x*cos(th2); c3=l1*sin(th1)+x*sin(th2);plot(s3,c3);hold on;endreturn代入pos.m中所得的结果，输入命令：draw(l1,l2,l3,l4,y(1,1),y(2,1),y(3,1)即可得到当前位置下的四杆图形。实例计算假设已知各杆的尺寸和杆1的初始角度th1及角速度w1，其分别为th1=60，w1=pi/3，l1=25，l2=40，l3=50，l4=60。现求在图1所示位置时的杆2和3的角位移，角速度和角加速度及当前位置下的四杆机构图形。仅需输入： th1=60;w1=pi/3;l1=25;l2=40;l3=50;l4=60; pos(th1,w1,l1,l2,l3,l4)draw(l1,l2,l3,l4,y(1,1),y(2,1),y(3,1)若要得到四杆机构的运动图形，则输入：for m=0:1:360;%步距为1，可更改；y=pos(m,w1,l1,l2,l3,l4);draw(l1,l2,l3,l4,y(1,1),y(2,1),y(3,1);%不断重画；end最后可得到杆2和杆3的角速度、角加速度曲线，如图示：希望能在直角坐标系中用MATLAB绘图d2r = pi/180;%角度精度r1 =0.220;%如图示r2 =0.140;%如图示r3 =0.280;%如图示r4 =0.310;%如图示r5 =0.800;%如图示n=3000;%转速w=2*pi*n/60;%角速度T=2*pi/w;%周期dt=T/360;%每度转过需要的时间% for loop for input (crank) anglei = 0;for t2deg = 0:1:360% update counteri = i+1;% convert angles to radianst2 = t2deg*d2r;% geometric computation of other anglese = sqrt(r12 + r22 - 2*r1*r2*cos(t2);a = asin(sin(t2)*r2/e);g = acos(-e2 + r32 + r42)/(2*r3*r4);b = asin(sin(g)*r3/e);t4 = 180*d2r - a - b;t3 = t4-g;% compute and store position of point of interest% where 30 is the distance along link 3 from point A% STUDENTS: MODIFY HEREax(i)=r2*cos(t2);%A点横坐标ay(i)=r2*sin(t2);%A点纵坐标bx(i)=r4*cos(t4)+r1;%B点横坐标by(i)=r4*sin(t4);%B点纵坐标m=-(r3+r5)/r5;dx(i) = (ax(i)+m*bx(i)/(1+m);%D点横坐标dy(i) = (ay(i)+m*by(i)/(1+m);%D点纵坐标%cx(i) = r2*(cos(t2)-cos(t2-d2t) +(r3+r5)*(cos(t3)-cos(t3-d3t)%cy(i) = r2*(sin(t2)-sin(t2-d2t) +(r3+r5)*(sin(t3)-sin(t3-d3t)end% generate a new figurehold onxlabel(A点横坐标 (m)ylabel(A点纵坐标 (m)title(A点位置)plot(ax,ay);grid offsubplot(221)hold onxlabel(B点横坐标 (m)ylabel(B点纵坐标 (m)title(B点位置)plot(bx,by)subplot(222)hold onxlabel(C点横坐标 (m)ylabel(C点纵坐标 (m)title(C点位置)plot(dx,dy)grid offsubplot(223)hold onxlabel(D点X轴方向速度 (m/s)ylabel(D点Y轴方向速度(m/s)title(D点速度)vdx=diff(dx);vdy=diff(dy);plot(vdx,vdy);subplot(224)hold onxlabel(D点X轴加方向速度 (m/s2)ylabel(D点Y轴加方向速度(m/s2)title(D点加速度)adx=diff(vdx);ady=diff(vdy);plot(adx,ady);速度和加速度如何表示 我的有问题 请高手帮帮看下disp * * * * * * 平面四杆机构的运动分析 * * * * * * L1=input(输入连架杆长度 L1=); L2=input(输入曲柄长度(主动件) L2=); L3=input(输入连杆长度 L3=); L4=input(输入输入从动杆#4长度 L4=); rs=L1 L2 L3 L4; minimum=min(rs); maximum=max(rs); if (minimum+maximum)=(sum(rs)-minimum-maximum)