平面四杆机构的设计与运动分析Matlab代码

平面四杆机构的设计与运动分析 M 代码平面四杆机构的设计 M 代码A=cos(50*pi/180),cos(50-35)*pi/180),1;cos(75*pi/180),cos(75-80)*pi/180),1;cos(105*pi/180), cos(105-125)*pi/180),1;B=cos(35*pi/180);cos(80*pi/180);cos(125*pi/180);P=ABm= P(1)n=-m/ P(2)l=sqrt(mA2+nA2+1-2*n*P(3)运行设计结果显示:Command Window A= Leos 0*pi/180), cos(50-35)#pi/180)1;cos r75*pi/180)cos (75- B=cos f35*pi/180):cos (80*pi/180);cos f125*pi/180);P=ABP (1)n=-m/ P(2)l=sqrt1m 2+n 2+l-2*n* P (3)P =1.5815-1.26401.0235m 二1.5815n =1.25121.5829%参数赋值clc,clearl0=1.2512;l1=1.0;l2=1.5829;l3=1.5815;M=-1;Omiga1=10;Theta1=0:0.01:360;Theta1=Theta1*pi/180;%求解各个构件位移、速度、加速度A=2*l1*l2*sin(Theta1);B=2*l2*(l1*cos(Theta1)-l0);C=ll人2+12人2+10人2-13人2-2*ll*10*cos(Thetal);E=2*l1*l3*sin(Theta1);F=2*13*(1l*cos(Thetal)-10);G=12A2-11A2-13A2-10A2+2*11*10*cos(Theta1);Theta3=2*atan(E+M*sqrt(E.A2+ F42- GA2)./(F-G);Theta2=2*atan(A+M*sqrt(A.A2+B.A2-C.A2)./(B-C);Omiga2=Omiga1*11*sin(Theta1-Theta3)./(12*sin(Theta3-Theta2);Omiga3=Omiga1*11*sin(Theta1-Theta2)./(13*sin(Theta3-Theta2);A1fa3=(Omiga1A2*11*cos(Theta1-Theta2)+Omiga2.A2*12-Omiga3.A2*13.*cos(Theta3-Theta2)./ (13*sin(Theta3-Theta2);A1fa2=(-Omiga1A2*11*cos(Theta1-Theta3)+Omiga3.A2*13-Omiga2.A2*12.*cos(Theta2-Theta3)./ (12*sin(Theta2-Theta3);%绘图Theta1=Theta1*180/pi;Subp1ot(3,1,1)p1ot(Theta1,Theta3*180/pi),grid onxlabel(曲柄转角(a。) );y1abe1(CD 杆角位移(人。);Subp1ot(3,1,2)p1ot(Theta1,Omiga3),grid onx1abe1(曲柄转角(人。)；y1abe1(CD 杆角速度(rad/s);Subp1ot(3,1,3)p1ot(Theta1,A1fa3),grid onx1abe1(曲柄转角(人。)；y1abe1(CD 杆角加速度(rad/sA2);by Xu jianping.CD杆角加速M(rad/s2)CD杆角速度(日d虫)CD杆角位移* ) o CDC3o5CDF-4 Ed_* -i InweTt Hoo-M- DfDGktop w-ndow Hfp-P回cddfiEI卑潜