哈工大(威海)机械原理_连杆14题

1、Harbin Institute of Technology课程名称： 机械原理大作业 设计题目： 连杆运动分析（14） 院 系： 班 级： 设 计 者： 学 号： 指导教师： 设计时间： 2015年5月 哈尔滨工业大学（威海）平面连杆机构运动分析1题目及原始数据题目中图上yG=210mm所标记的位置错误，应该为G到A点间纵坐标间的距离。2机构结构分析、组成机构的基本杆组划分拆分基本杆组： RR RRR RRR构件1与机架构成一级机构构件2和3构成二级杆组构件4与构件5构成二级杆组。所以整个连杆机构为二级机构3各基本杆组的运动分析及数学模型以A点为原点建立坐标系，可根据题目条件得到个点坐标，如

2、下：A（0，0） D（300,-500） G（-430，210）设AB以的角速度逆时针旋转，B点坐标随时间的变化公式：C点坐标：xc=xd+CD*cos(a2+a3)yc=yd+CD*sin(a2+a3)E点坐标： F点坐标： 令=a7=30° EGF=a104用 matlab编程（程序）4.1 构件5的角位移xa=0;ya=0;xd=300;yd=-500;xg=-430;yg=210;fe=265;ab=100;gf=670;be=460;bc=460;cd=250;a7=pi/6;ec=2*bc*sin(a7/2);a1=0:0.001:2*pi;t=a1./10;xb=ab.

3、*cos(a1);yb=ab.*sin(a1);a2=atan(yd-yb)./(xd-xb)+pi;bd=sqrt(yb-yd).*(yb-yd)+(xb-xd).*(xb-xd);a3=acos(bd.*bd+cd.*cd-bc.*bc)./(2*bd.*cd);xc=xd+cd.*cos(a2+a3);yc=yd+cd.*sin(a2+a3);for i=1:length(a1) if atan(yc(i)-yb(i)./(xc(i)-xb(i)>0 a4(i)=atan(yc(i)-yb(i)./(xc(i)-xb(i); else a4(i)=atan(yc(i)-yb(i).

4、/(xc(i)-xb(i)+pi; endendfor i=1:length(a1) if a4(i)+5*pi/12<pixe=xc-ec.*sin(a4-1/12*pi);ye=yc+ec.*cos(a4-1/12*pi);else xe(i)=xc(i)-ec.*cos(a4(i)-7*pi/12);ye(i)=yc(i)-ec.*sin(a4(i)-7*pi/12);endendfor i=1:length(a1) if atan(ye(i)-yg)./(xe(i)-xg)>0 a5(i)=atan(ye(i)-yg)./(xe(i)-xg); else a5(i)=ata

5、n(ye(i)-yg)./(xe(i)-xg)+pi; endendge=sqrt(ye-yg).*(ye-yg)+(xe-xg).*(xe-xg);a6=acos(ge.*ge+fe.*fe-gf.*gf)./(2*ge.*fe);xf=xe+fe.*cos(a5+a6);yf=ye+fe.*sin(a5+a6);a10=acos(ge.*ge+gf*gf-fe*fe)./(2*gf.*ge);plot(t,a5-a10);>> title('构件5的角位移');xlabel('t/s');ylabel('角度rad');grid

6、on;4.2构件5的角速度xa=0;ya=0;xd=300;yd=-500;xg=-430;yg=210;fe=265;ab=100;gf=670;be=460;bc=460;cd=250;a7=pi/6;ec=2*bc*sin(a7/2);a1=0:0.001:2*pi;t=a1./10;xb=ab.*cos(a1);yb=ab.*sin(a1);a2=atan(yd-yb)./(xd-xb)+pi;bd=sqrt(yb-yd).*(yb-yd)+(xb-xd).*(xb-xd);a3=acos(bd.*bd+cd.*cd-bc.*bc)./(2*bd.*cd);xc=xd+cd.*cos(

7、a2+a3);yc=yd+cd.*sin(a2+a3);for i=1:length(a1) if atan(yc(i)-yb(i)./(xc(i)-xb(i)>0 a4(i)=atan(yc(i)-yb(i)./(xc(i)-xb(i); else a4(i)=atan(yc(i)-yb(i)./(xc(i)-xb(i)+pi; endendfor i=1:length(a1) if a4(i)+5*pi/12<pixe=xc-ec.*sin(a4-1/12*pi);ye=yc+ec.*cos(a4-1/12*pi);else xe(i)=xc(i)-ec.*cos(a4(i)-

8、7*pi/12);ye(i)=yc(i)-ec.*sin(a4(i)-7*pi/12);endendfor i=1:length(a1) if atan(ye(i)-yg)./(xe(i)-xg)>0 a5(i)=atan(ye(i)-yg)./(xe(i)-xg); else a5(i)=atan(ye(i)-yg)./(xe(i)-xg)+pi; endendge=sqrt(ye-yg).*(ye-yg)+(xe-xg).*(xe-xg);a6=acos(ge.*ge+fe.*fe-gf.*gf)./(2*ge.*fe);xf=xe+fe.*cos(a5+a6);yf=ye+fe.*

9、sin(a5+a6);a10=acos(ge.*ge+gf*gf-fe*fe)./(2*gf.*ge);plot(t(1:6283),diff(a5-a10)./diff(t);title('构件5的角速度');xlabel('t/s');ylabel('角速度rad/s');grid on;4.3 构件5的角加速度xa=0;ya=0;xd=300;yd=-500;xg=-430;yg=210;fe=265;ab=100;gf=670;be=460;bc=460;cd=250;a7=pi/6;ec=2*bc*sin(a7/2);a1=0:0.00

10、1:2*pi;t=a1./10;xb=ab.*cos(a1);yb=ab.*sin(a1);a2=atan(yd-yb)./(xd-xb)+pi;bd=sqrt(yb-yd).*(yb-yd)+(xb-xd).*(xb-xd);a3=acos(bd.*bd+cd.*cd-bc.*bc)./(2*bd.*cd);xc=xd+cd.*cos(a2+a3);yc=yd+cd.*sin(a2+a3);for i=1:length(a1) if atan(yc(i)-yb(i)./(xc(i)-xb(i)>0 a4(i)=atan(yc(i)-yb(i)./(xc(i)-xb(i); else a

11、4(i)=atan(yc(i)-yb(i)./(xc(i)-xb(i)+pi; endendfor i=1:length(a1) if a4(i)+5*pi/12<pixe=xc-ec.*sin(a4-1/12*pi);ye=yc+ec.*cos(a4-1/12*pi);else xe(i)=xc(i)-ec.*cos(a4(i)-7*pi/12);ye(i)=yc(i)-ec.*sin(a4(i)-7*pi/12);endendfor i=1:length(a1) if atan(ye(i)-yg)./(xe(i)-xg)>0 a5(i)=atan(ye(i)-yg)./(xe(

12、i)-xg); else a5(i)=atan(ye(i)-yg)./(xe(i)-xg)+pi; endendge=sqrt(ye-yg).*(ye-yg)+(xe-xg).*(xe-xg);a6=acos(ge.*ge+fe.*fe-gf.*gf)./(2*ge.*fe);xf=xe+fe.*cos(a5+a6);yf=ye+fe.*sin(a5+a6);a10=acos(ge.*ge+gf*gf-fe*fe)./(2*gf.*ge);yv=diff(a5-a10)./diff(t);plot(t(1:length(diff(yv)./diff(t(1:6283),diff(yv)./diff(t(1:6283);title('构件5的角加速度');xlabel('t/s');ylabel(