aby=aay-l1*w1.^2.*sin(r*f1);
% RRR2级杆组运动分析
%输入D点参数
l2=428;
l3=214;
Xd=455;
Yd=0;
Vdx=0;
Vdy=0;
adx=0;
ady=0;
%计算E点、2杆、3杆运动参数
lbe=840;
lce=486;
a0=2*l2*(Xd-Xb);
b0=2*l2*(Yd-Yb);
aex=abx-lbe*(e2.*sin(f2-w)+w2.^2.*cos(f2-w));
aey=aby+lbe*(e2.*cos(f2-w)-w2.^2.*sin(f2-w));
%计算杆5运动参数
Xf=;
Yf=;
l5=sqrt((Xe-Xf).^2+(Ye-Yf).^2);
dX=Xe-Xf;
dY=Ye-Yf;
Warning: Unable to interpret TeX string "Xe/\mm"
5、计算结果
图一：E点的运动轨迹
图二：5杆角位移
图三：5杆角速度
图四：5杆角加速度
6、计算结果分析
由E点位置图像可看出，构件4做周期往复运动，由图二、三、四可看出，构件5的角位移、角速度、角加速度均成周期性变化。构件5的角位移曲线为类正弦曲线，从与之对应的角速度曲线中可看出5杆做往复摆动运动，推程时速度较快，回程时速度较慢。
3)加速度方程
两杆角加速度为
式中，
点E加速度方程为
RPRⅡ级杆组的运动分析
（1）位移方程
（2）速度方程
其中
（3）加速度方程
4、 计算编程
利用MATLAB软件进行编程，程序如下:
%点B和AB杆运动状态分析
>>r=pi/180;
w1=10;
e1=0;
l1=100;
Xa=0;
Ya=0;
Vax=0;
Vay=0;
f3(n)=pi+atan(dY(n)/dX(n));
elseif dX(n)==0&dY(n)<0
f3(n)=*pi;
elseif dX(n)>0&dY(n)<0
f3(n)=2*pi+atan(dY(n)/dX(n));
end
end
C2=l2*cos(f2);
C3=l3*cos(f3);
S2=l2*sin(f2);
figure(3);plot(f1,w5,'k');xlabel('f/\circ');ylabel('w5/rad/s');grid on;title('5杆角速度');
figure(4);plot(f1,a5,'k');xlabel('f/\circ');ylabel('a5/rad/s2');grid on;title('5杆角加速度');
Harbin Institute of Technology
机械原理大作业一
课程名称：机械原理
设计题目：连杆机构运动分析
院系：机电学院
班级：1208105
分析者：殷琪
学号：
指导教师：丁刚
设计时间：
哈尔滨工业大学
设计说明书
1 、题目
如图所示机构，一只机构各构件的尺寸为AB=100mm，BC=，CE=，BE=，CD=，AD=，AF=7AB，DF=，∠BCE=139?。构件1的角速度为ω1=10rad/s，试求构件2上点E的轨迹及构件5的角位移、角速度和角加速度，并对计算结果进行分析。
aax=0;
aay=0;
f1=0:1: 360;
% B点位置
Xb=Xa+l1*cos(r*f1);
Yb=Ya+l1*sin(r*f1);
% B点速度
Vbx=Vax-w1*l1*sin(r*f1);
Vby=Vay+w1*l1*cos(r*f1);
% B点加速度
abx=aax-l1*w1.^2.*cos(r*f1);
for n=1:length(dX)
if dX(n)>0&dY(n)>=0
f3(n)=atan(dY(n)/dX(n));
elseif dX(n)==0&dY(n)>0
f3(n)=pi/2;
elseif dX(n)<0&dY(n)>=0
f3(n)=pi+atan(dY(n)/dX(n));
elseif dX(n)<0&dY(n)<0
S3=l3*sin(f3);
G1=C2.*S3-C3.*S2;
w2=(C3.*(Vdx-Vbx)+S3.*(Vdy-Vby))./G1;
w3=(C2.*(Vdx-Vbx)+S2.*(Vdy-Vby))./G1;
G2=adx-abx+(w2.^2).*C2-(w3.^2).*C3;
G3=ady-aby+(w2.^2).*S2-(w3.^2).*S3;
for n=1:length(dX)
if dX(n)>0&dY(n)>=0
f5(n)=atan(dY(n)/dX(n));
elseif dX(n)==0&dY(n)>0
f5(n)=pi/2;
elseif dX(n)<0&dY(n)>=0
f5(n)=pi+atan(dY(n)/dX(n)ห้องสมุดไป่ตู้;
elseif dX(n)<0&dY(n)<0
c0=l2^2+(Xb-Xd).^2+(Yb-Yd).^2-l3^2;
f2=2*atan((b0+sqrt(a0.^2+b0.^2-c0.^2))./(a0+c0));
% C点位置
Xc=Xb+l2*cos(f2);
Yc=Yb+l2*sin(f2);
% 2杆、3杆运动参数计算
dX=Xc-Xd;
dY=Yc-Yd;
将上式对时间t求导，可得加速度方程:
RRRⅡ级杆组的运动分析
如下图所示当已知RRR杆组中两杆长LBC、LCD和两外副B、D的位置和运动时，求内副C的位置、两杆的角位置、角运动以及E点的运动。
1)位置方程
由移项消去 后可求得 ：
式中，
可求得 ：
E点坐标方程：
其中
2)速度方程
两杆角速度方程为
式中，
点E速度方程为
a5=(-aex.*sin(f5)+aey.*cos(f5))./l5;
%画出各参数曲线
figure(1);plot(Xe,Ye,'k');xlabel('Xe/\mm');ylabel('Ye/mm');grid on;title('E点位置');
figure(2);plot(f1,f5,'k');xlabel('f/\circ');ylabel('f5/\circ');grid on;title('5杆角位移');
2、机构结构分析
该机构由6个构件组成，4和5之间通过移动副连接，其他各构件之间通过转动副连接，主动件为杆1，杆2、3、4、5为从动件，2和3组成Ⅱ级RRR基本杆组，4和5组成Ⅱ级RPR基本杆组。
如图建立坐标系
3、各基本杆组的运动分析数学模型
1)位置分析
2)速度和加速度分析
将上式对时间t求导，可得速度方程：
e2=(G2.*C3+G3.*S3)./G1;
% E点位置
w=acos((l2^2+lbe^2-lce^2)/(2*l2*lbe));
Xe=Xb+lbe*cos(f2-w);
Ye=Yb+lbe*sin(f2-w);
Vex=Vbx-lbe*w2.*sin(f2-w);
Vey=Vby+lbe*w2.*cos(f2-w);
f5(n)=pi+atan(dY(n)/dX(n));
elseif dX(n)==0&dY(n)<0
f5(n)=*pi;
elseif dX(n)>0&dY(n)<0
f5(n)=2*pi+atan(dY(n)/dX(n));
end
end
w5=(-Vex.*sin(f5)+Vey.*cos(f5))./l5;