平面四杆机构的设计与运动分析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=A\B

m= P(1)

n=-m/ P(2)

l=sqrt(m^2+n^2+1-2*n*P(3))

运行设计结果显示:

平面四杆机构的运动分析M代码----

%参数赋值

clc,clear

l0=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=l1^2+l2^2+l0^2-l3^2-2*l1*l0*cos(Theta1);

E=2*l1*l3*sin(Theta1);

F=2*l3*(l1*cos(Theta1)-l0);

G=l2^2-l1^2-l3^2-l0^2+2*l1*l0*cos(Theta1);

Theta3=2*atan((E+M*sqrt(E.^2+ F.^2- G.^2))./(F-G));

Theta2=2*atan((A+M*sqrt(A.^2+B.^2-C.^2))./(B-C));

Omiga2=Omiga1*l1*sin(Theta1-Theta3)./(l2*sin(Theta3-Theta2));

Omiga3=Omiga1*l1*sin(Theta1-Theta2)./(l3*sin(Theta3-Theta2));

Alfa3=(Omiga1^2*l1*cos(Theta1-Theta2)+Omiga2.^2*l2-Omiga3.^2*l3.*cos(Theta3-Theta2))./ (l3*sin(Theta3-Theta2));

Alfa2=(-Omiga1^2*l1*cos(Theta1-Theta3)+Omiga3.^2*l3-Omiga2.^2*l2.*cos(Theta2-Theta3))./ (l2*sin(Theta2-Theta3));

%绘图

Theta1=Theta1*180/pi;

Subplot(3,1,1)

plot(Theta1,Theta3*180/pi),grid on

xlabel('曲柄转角(^。) ');ylabel('CD杆角位移(^。) ');

Subplot(3,1,2)

plot(Theta1,Omiga3),grid on

xlabel('曲柄转角(^。) ');ylabel('CD杆角速度(rad/s) ');

Subplot(3,1,3)

plot(Theta1,Alfa3),grid on

xlabel('曲柄转角(^。) ');ylabel('CD杆角加速度(rad/s^2) ');

by Xu jianping

CD杆运动规律