哈工大-机械原理大作业1-连杆-22
哈工大机械原理大作业1 连杆机构运动分析 牛头刨床 学号5号
大作业1 连杆机构运动分析运动分析题目图1-13 牛头刨床机构简图一、按比例画出机构运动简图见A3图纸二、机构的结构分析及基本杆组的划分（1）机构的结构分析该机构为平面连杆机构，活动构件数n=5；有A、B、C、E、F五个转动副和B、F两个移动副，没有高副，故低副数PL=7，高副数PH=0。
机构的自由度F=3n-2PL-PH=3*5-2*7=1。
机构有确定运动轨迹。
（2）基本杆组的划分图1-13中，原动件构件1为Ⅰ级杆组，移除之后按照拆杆组法对机构进行拆分，可得由滑块2和杆3组成的RPRⅡ级杆组，以及由杆4和滑块5组成的RRPⅡ级杆组。
拆分图如下：图1-13-1 Ⅰ级杆组图1-13-2 RPRⅡ级杆组图1-13-3 RRPⅡ级杆组三、各杆组的运动分析数学模型图1-13-1 Ⅰ级杆组 x B =x A +l AB cos ϕy B =y A +l AB sin ϕ图1-13-2 RPR Ⅱ级杆组 x E =x B +(l CE -s)cos ϕj y E =y B +(l CE -s)sin ϕjϕj =arctan(B 0/A 0)A 0=xB -x CB 0=y B -y CS =2020B A + 图1-13-3 RRP Ⅱ级杆组F 的位置 x F =x E -l EF cos ϕiy F =y E +l EF sin ϕi ϕi =arcsin(A 1/l EF )A 1=H-H 1-y EF 的速度F 的加速度四、建立坐标系以A 为原点建立坐标系，则A （0,0）五、编程使用MATLAB 编写程序如下lAB=108;lCD=620;lEF=300;H1=350;H=635;xC=0;yC=-350;syms t;f=(255*pi/30)*t;xB=lAB*cos(f);yB=lAB*sin(f);A0=xB-xC;B0=yB-yC;S=sqrt(A0^2+B0^2);f1=atan(B0/A0);xE=xB+(lAB-S)*cos(f1);yE=yB+(lAB-S)*sin(f1);A1=H-H1-yE;f2=asin(A1/lEF);xF=xE-lEF*cos(f2);yF=yE+lEF*sin(f2);a=0:0.001:0.5;xF=subs(xF,t,a);plot(a,xF)title('位移x随时间t变化图像') xlabel('t(s)'),ylabel(X')lAB=108;lCD=620;lEF=300;H1=350;H=635;xC=0;yC=-350;syms t;f=(255*pi/30)*t;xB=lAB*cos(f);yB=lAB*sin(f);A0=xB-xC;B0=yB-yC;S=sqrt(A0^2+B0^2);f1=atan(B0/A0);xE=xB+(lAB-S)*cos(f1);yE=yB+(lAB-S)*sin(f1);A1=H-H1-yE;f2=asin(A1/lEF);xF=xE-lEF*cos(f2);yF=yE+lEF*sin(f2);vF=diff(xF,t);a=0:0.001:0.5;vF=subs(vF,t,a);plot(a,vF)title('速度v随时间t变化图像') xlabel('t(s)'),ylabel(V')lAB=108;lCD=620;lEF=300;H1=350;H=635;xC=0;yC=-350;syms t;f=(255*pi/30)*t;xB=lAB*cos(f);yB=lAB*sin(f);A0=xB-xC;B0=yB-yC;S=sqrt(A0^2+B0^2);f1=atan(B0/A0);xE=xB+(lAB-S)*cos(f1);yE=yB+(lAB-S)*sin(f1);A1=H-H1-yE;f2=asin(A1/lEF);xF=xE-lEF*cos(f2);yF=yE+lEF*sin(f2);aF=diff(xF,t,2);a=0:0.001:0.5;aF=subs(aF,t,a);plot(a,aF)title('加速度a随时间t变化图像') xlabel('t(s)'),ylabel(A')六、位移、速度、加速度图像位移xE随时间t变化图像速度vE随时间t变化图像加速度aE随时间t变化图像七、运动分析由MATLAB编程后绘制所得图线分析可知，牛头刨床的刨头F行程约为500mm，在0.06s~0.17s之间速度和加速度波动较大，在0.17s~0.29s 之间速度和加速度波动较小，这就保证了刨头在空行程时有急回运动，而在工作行程是由较均匀的切削速度，运动连续，且运动性能良好稳定。
哈尔滨工业大学机械原理大作业-连杆机构运动分析
%打印图像
figure(1);
plot(fi1,sF);
title('位移变化曲线');
figure(2);
plot(fi1,vF);
title('速度变化曲线');
figure(3);
plot(fi1,aF);
title('加速度变化曲线');
六、计算结果
图8：推杆位移变化曲线
图9：推杆速度变化曲线
xD=400;yD=500;vDx=0;vDy=0;aDx=0;aDy=0;
xK=0;yK=600;vKx=0;vKy=0;aKx=0;aKy=0;
l1=150;l2=600;l3=500;lBE=480;l4=600;
fi5=pi;
n1=50;
w1=2*pi*n1/60;
fi1=linspace(0,2*pi,1000);
aEy(i)=aBy(i)-w2(i)^2*lBE*sin(fi2(i))+a2(i)*lBE*cos(fi2(i));
%求F点的运动参数
A1(i)=(yE(i)-yK)*cos(fi5)-(xE(i)-xK)*sin(fi5);
fi4(i)=fi5-asin(A1(i)/l4);
xF(i)=xE(i)+l4*cos(fi4(i));
（2）速度和加速度分析
B点的速度
B点的加速度
2、由“RRRⅡ级杆组”，已知B点和D点的运动参数，可求得构件2、构件3的运动参数
图5
D点的坐标方程
D点的速度
D点的加速度
（1）位置方程
构件2、构件3的长度
先求出 和
其中
则可求得
哈工大机械原理大作业
H a r b i n I n s t i t u t e o f T e c h n o l o g y机械原理大作业一课程名称：机械原理设计题目：连杆机构运动分析院系：机电学院班级： 1208105分析者：殷琪学号：指导教师：丁刚设计时间：哈尔滨工业大学设计说明书1 、题目如图所示机构，一只机构各构件的尺寸为AB=100mm，BC=4.28AB，CE=4.86AB，BE=8.4AB，CD=2.14AB，AD=4.55AB，AF=7AB，DF=3.32AB，∠BCE=139?。
构件1的角速度为ω1=10rad/s，试求构件2上点E的轨迹及构件5的角位移、角速度和角加速度，并对计算结果进行分析。
2、机构结构分析该机构由6个构件组成，4和5之间通过移动副连接，其他各构件之间通过转动副连接，主动件为杆1，杆2、3、4、5为从动件，2和3组成Ⅱ级RRR基本杆组，4和5组成Ⅱ级RPR 基本杆组。
如图建立坐标系3、各基本杆组的运动分析数学模型1) 位置分析2) 速度和加速度分析 将上式对时间t 求导，可得速度方程：将上式对时间t 求导，可得加速度方程:RRR Ⅱ级杆组的运动分析如下图所示 当已知RRR 杆组中两杆长L BC 、L CD 和两外副B 、D 的位置和运动时，求内副C的位置、两杆的角位置、角运动以及E 点的运动。
1) 位置方程由移项消去j ϕ后可求得i ϕ：式中，可求得j ϕ：E 点坐标方程：其中2) 速度方程两杆角速度方程为式中，点E 速度方程为3) 加速度方程两杆角加速度为式中，点E 加速度方程为RPR Ⅱ级杆组的运动分析（1） 位移方程（2）速度方程其中（3）加速度方程4、 计算编程利用MATLAB 软件进行编程，程序如下:% 点B 和AB 杆运动状态分析>>r=pi/180;w 1=10;e 1=0;l 1=100;Xa=0;Ya=0;Vax=0;Vay=0;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);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);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;for n=1:length(dX)if dX(n)>0&dY(n)>=0f3(n)=atan(dY(n)/dX(n));elseif dX(n)==0&dY(n)>0f3(n)=pi/2;elseif dX(n)<0&dY(n)>=0f3(n)=pi+atan(dY(n)/dX(n));elseif dX(n)<0&dY(n)<0f3(n)=pi+atan(dY(n)/dX(n));elseif dX(n)==0&dY(n)<0f3(n)=1.5*pi;elseif dX(n)>0&dY(n)<0f3(n)=2*pi+atan(dY(n)/dX(n));endendC2=l2*cos(f2);C3=l3*cos(f3);S2=l2*sin(f2);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;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);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=646.2912088;Yf=-268.9008617;l5=sqrt((Xe-Xf).^2+(Ye-Yf).^2);dX=Xe-Xf;dY=Ye-Yf;for n=1:length(dX)if dX(n)>0&dY(n)>=0f5(n)=atan(dY(n)/dX(n));elseif dX(n)==0&dY(n)>0f5(n)=pi/2;elseif dX(n)<0&dY(n)>=0f5(n)=pi+atan(dY(n)/dX(n));elseif dX(n)<0&dY(n)<0f5(n)=pi+atan(dY(n)/dX(n));elseif dX(n)==0&dY(n)<0f5(n)=1.5*pi;elseif dX(n)>0&dY(n)<0f5(n)=2*pi+atan(dY(n)/dX(n));endendw5=(-Vex.*sin(f5)+Vey.*cos(f5))./l5;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杆角位移');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');gridon;title('5杆角加速度');Warning: Unable to interpret TeX string "Xe/\mm"5、计算结果图一：E点的运动轨迹图二：5杆角位移图三：5杆角速度图四：5杆角加速度6、计算结果分析由E点位置图像可看出，构件4做周期往复运动，由图二、三、四可看出，构件5的角位移、角速度、角加速度均成周期性变化。
哈工大机械原理大作业连杆
Harbin Institute of Technology机械原理大作业一课程名称： 机械原理 设计题目： 连杆机构运动分析 院 系： 机电工程学院 班 级： 设 计 者： 学 号： 指导教师： 设计时间：1.运动分析题目(11)在图所示的六杆机构中，已知：AB l =150mm, AC l =550mm, BD l =80mm, DE l =500mm,曲柄以等角速度1w =10rad/s 沿逆时针方向回转，求构件3的角速度、角加速度和构件5的位移、速度、加速度。
2.机构的结构分析2.1建立以点A 为原点的固定平面直角坐标系A-x, y,如下图： 2.2机构结构分析该机构由Ⅰ级杆组RR （原动件1）、Ⅱ级杆组RPR （杆2及滑块3）和Ⅱ级杆组RRP （杆4及滑块5）组成。
3.建立组成机构的各基本杆组的运动分析数学模型3.1原动件1（Ⅰ级杆组RR ）由图所示，原动件杆1的转角a=0-360°，角速度1w =10rad/s ，角加速度1a =0,运动副A 的位置坐标A x =A y =0,速度(A ，A)，加速度（A ，A ），原动件1的长度AB l =150mm 。
求出运动副B 的位置坐标（B x , B y ）、速度（B ，B ）和加速度（B ，B ）。
3.2杆2、滑块3杆组（RPR Ⅱ级杆组）已出运动副B 的位置（B x , B y ）、速度（B ，B ）和加速度（B ，B ），已知运动副C 的位置坐标C x =0, C y =550mm,速度，加速度，杆长AC l =550mm 。
求出构件2的转角b,角速度2w 和角加速度2a . 3.3构件二上点D 的运动已知运动副B 的位置（B x , B y ）、速度（B ，B ）、加速度（B ，B ），已经求出构件2的转角b ，角速度2w 和角加速度2a ，杆BD 的长度BD l =80mm 。
根据Ⅰ级杆组RR 的运动分析数学模型求出点D 的位置坐标（D x ,D y ）、速度（D ，D ）和加速度（D ，D ）。
哈工大机械原理大作业连杆
哈工大机械原理大作业-连杆连杆是机械原理中常见的机构之一，也是机械工程中非常重要的部件。
它由两个旋转接头和一个连接两个旋转接头的杆件组成。
连杆广泛应用于各种机械设备中，如汽车发动机、泵、机床等。
本文将介绍连杆的工作原理、应用以及设计要点。
连杆的工作原理是将旋转运动转化为直线运动或将直线运动转化为旋转运动。
它通过两个旋转接头的运动将杆件上的一个点的运动转化为另一个点的运动。
连杆的运动有两种基本形式：一是曲柄连杆机构，二是摇杆连杆机构。
曲柄连杆机构中，一个旋转接头为曲柄，另一个旋转接头为连杆；摇杆连杆机构中，一个旋转接头为摇杆，另一个旋转接头为连杆。
连杆广泛应用于各种机械设备中。
在汽车发动机中，连杆将曲轴的旋转运动转化为活塞的直线运动，从而驱动汽缸的工作；在泵中，连杆将电机的旋转运动转化为柱塞的直线运动，从而产生压力；在机床中，连杆将电机的旋转运动转化为工作台的直线运动，从而实现加工。
设计连杆时需要考虑一些要点。
首先是连杆的材料选择和尺寸设计。
连杆需要承受较大的力和扭矩，因此需要选择具有较高强度和刚度的材料。
同时，根据应用需求和力学原理，设计连杆的尺寸，以确保其能够承受正常工作条件下的负荷。
其次是连杆的润滑和密封。
连杆在工作过程中需要润滑剂来减少摩擦和磨损，同时需要密封装置来防止润滑剂泄漏。
因此，设计连杆时需要考虑润滑剂的供给和密封装置的设计。
最后是连杆的制造和装配。
连杆的制造需要保证其精度和质量，以确保其运转平稳和可靠。
在装配过程中，需要按照设计要求进行装配，同时进行必要的调试和检测，以确保连杆的工作性能符合要求。
总之，连杆是机械工程中非常重要的部件，广泛应用于各种机械设备中。
设计和制造连杆需要考虑材料选择、尺寸设计、润滑和密封以及制造和装配等方面的要点。
通过合理的设计和制造，可以确保连杆的工作性能和可靠性，从而提高机械设备的工作效率和寿命。
哈工大-机械原理大作业-连杆机构运动分析
机械原理大作业（一）作业名称：连杆机构运动分析设计题目： 20院系：英才学院班级： XXXXXXX设计者：邵广斌学号： XXXXXXXXXX指导教师：林琳设计时间： 2013年05月19日哈尔滨工业大学机械设计1.运动分析题目如图所示机构，已知机构各构件的尺寸为150AB mm =，97β=︒，400BC mm =，300CD mm =，320AD mm =，100BE mm =，230EF mm =，400FG mm =，构件1的角速度为110/rad s ω=,试求构件2上点F 的轨迹及构件5上点G 的位移、速度和加速度，并对计算结果进行分析。
2. 机构分析该机构由原动件AB （Ⅰ级杆组）、BCD （RRR Ⅱ级杆组）和FG （RRP Ⅱ级杆组）组成。
3. 建立坐标系如图3，建立以定点A 为原点的平面直角坐标系A-xy 。
图1 运动机构结构图4. 运动分析数学模型4.1 原动件AB原动件AB 的转角： 10~2ψπ= 原动件AB 的角速度：110/rad s ω=原动件AB 的角加速度： 10α= 运动副A 的位置坐标： 0A x = 0A y =运动副A 的速度： 0xA v = 0yA v = 运动副A 的加速度： 0xA a = 0yA a =原动件AB 长度：150AB l mm =运动副B 的位置坐标： 1B A AB x x l cos ψ=+1B A AB y x l sin ψ=+运动副B 的速度： 11 xB xA AB v v l sin ωψ=-11 yB yA AB v v l cos ωψ=+运动副B 的加速度： 2 1111 xBxA AB AB a a l cos l sin ωψαψ=--21111yB yA AB AB a a l sin l cos ωψαψ=-+4.2 RRR Ⅱ级杆组BCD运动副D 的位置坐标： 320D x mm = 0D y = 运动副D 的速度： 0xD v = 0yD v = 运动副D 的加速度： 0xD a = 0yD a = 杆BC 长度： 400BC l mm = 杆CD 长度：300CD l mm =BC 相对于x 轴转角：200ψ=其中02BC D B A l x x =-（） 0 2 BC D B B l y y =-（）2220B B C C l C l D l D =+- 222())(BDD B D B l x x y y =-+- CD 相对于x 轴转角： 3C DC Dy y arctanx x ψ-=-求导可得BC 角速度2ω、角加速度2α以及CD 角速度3ω、角加速度3α。
机械原理大作业(一)--连杆(34题)
机械原理大作业（一）作业名称：连杆机构运动分析设计题目：（34）题院系：船舶学院班级： 1213101设计者：学号：哈尔滨工业大学机械设计一、运动分析题目如图所示机构，已知机构各构件的尺寸为BF=200mm,EF=1.25BF,DE=1.13BF,EH=0.85BF,HF=0.65BF,CH=0.81BF,GC=1.56BF,BD=0.58BF,BG=1.85BF,GD=1.6BF,构件1的角速度为W1=10rad/s，试求构件2上点H的轨迹及构件5的角位移，角速度，角加速度，并对计算结果进行分析。
机构结构分析：二、机构的结构分析及基本杆组划分机构各构件都在同一平面内运动，活动构件数n=5，=7，=0则机构的自由度为:F=3×n-2×-1×=3×5-2×7-0=12.基本杆组划分（1）去除虚约束和局部自由度本机构中无虚约束或局部自由度。
（2）拆杆组。
从远离原动件（即杆1）进行拆分，就可以得到由杆4,5 组成的RRRⅡ级杆组GCH，2,3 组成的RRRⅡ级杆组EDF，最后剩下Ⅰ级机构杆1。
（3）确定机构的级别 由（2）知，机构为Ⅱ级机构三、各基本杆组的运动分析数学模型1）Ⅰ级杆组BF （原动件）在Ⅰ级杆组BF 中，即已知构件上B 点的运动参数，求同一构件上F 点（回转副）的运动参数。
调用Ⅰ级机构子程序即可求解 ①位置分析 由图可得F 点的矢量方程F B BF r r l =+x,y 轴上的投影坐标方程为cos sin F B BF BF F B BF BF x x l y y l ϕϕ=+⋅⎫⎬=+⋅⎭（1）②速度和加速度分析 将式（1）对时间t 求导即可得出速度方程：cos sin BF BF BF BF FF B BF FF B BF dx x x l dtdy y x l dt ϕϕϕϕ⋅⋅⋅⋅⋅⋅⎫⎪==-⎪⎬⎪==+⎪⎭ （2）2222cos sin sin cos BF BF BF BF BF BF BF FF B BF BF BF F F B BFBF d x x x l l dtd y y y l l dt ϕϕϕϕϕϕϕϕ⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⎫⎪==--⎪⎬⎪==--⎪⎭（3）其中因为设B 为原点：B x =0；B y =0 ；B x ⋅=0 ；B y ⋅=0 ；B x ⋅⋅=0 ；B y ⋅⋅=0由上（1）（2）（3）方程可求出F 点的位移，速度，加速度2）RRR Ⅱ级杆组DEF 分析，求出F 点的角位移，角速度，角加速度上面1）中已求得F 点的位移，速度，加速度。
哈工大机械原理大作业一连杆-15题
Harbin Institute of Technology机械原理大作业一课程名称：机械原理设计题目：连杆传动设计院系：船舶与海洋工程学院班级：1513102班分析者：刘康哲学号：*********指导教师：***设计时间：2017.06.051.题目：如图所示机构，已知机构各构件的的尺寸为l AB=100mm,l BD=400mm,l AE=400mm,l BC=200mm,曲柄AB的角速度为w1=10rad/s,试求构件2上点D的轨迹，构件5的角位移、角速度和角加速度，并对结果进行分析。
（题中构件尺寸满足l BD-l AB<l AE<l BD+l AB）。
2．机构的结构分析及基本杆组划分机构各构件都在同一平面内运动，可拆分成一个原动件（构件1），一个RRPⅡ级基本杆组（构件2、3）和一个RPR Ⅱ级基本杆组（构件4、5）。
活动构件数n=5，P L，P H，因此，机构的自由度为F=3*n-2* P L- P H =3*5-2*7=1原动件：RRPⅡ级基本杆组：RPR Ⅱ级基本杆组:3．坐标系的建立以A点为坐标原点，AE方向为x轴正方向，建立平面直角坐标系。
4．各基本杆组的运动分析数学模型（1）原动件：位置分析： X B = X A + L AB * cos(wt) Y B = Y A + L AB * sin(wt) 其中，X A = Y A = 0 ，w = 10rad/s , L AB = 100 mm 速度分析：将上式分别对时间求导，可得速度方程,V XB = V XA – w*L AB *sin(wt)V YB = V YA + w*L AB *cos(wt)其中，V XA = V YA = 0 加速度分析:将上式分别对时间求导，可得速度方程,a XA = a XA – α* L AB * sin(wt) – w*w*L AB *cos(wt)a YA = a YA + α* L AB * cos(wt) – w*w*L AB *sin(wt)其中，a XA = a YA = 0 α= dw/dt（2）RRP Ⅱ级基本杆组:位置分析： X C = LAB*cos(wt)+wt Lab Lbc 222sin *Y C = 0 可求得，β= arcsin (Y B / L BC ) 速度分析：将上式分别对时间求导，可得速度方程,V XC =-w* LAB*cos(wt) – (wt Lab Lbc 222sin * )-1/2 * sin(wt)*wV YC = 0加速度分析：将上式分别对时间求导，可得加速度方程, （3）RPR Ⅱ级基本杆组位置分析： X D = X C + L CD * cos βY D = - L CD * sin β可求得， θ = arctan (-YD / (LAE -XD) ) 速度分析：将上式分别对时间求导，可得速度方程,5．计算编程在MATLAB 中编写如下程序：function [ XD,YD,Q,Q1,Q2] = calculate(t) %UNTITLED4 Summary of this function goes here % Detailed explanation goes here Lab=100; Lbc=200; Lcd=200; Lae=400; w=10;XB=Lab*cos(w.*t); YB=Lab*sin(w.*t); P=asin(YB./Lbc); XC=sqrt(Lbc.^2-YB.^2)+XB; XD=Lcd.*cos(P)+XC; YD=-Lcd.*sin(P); Q=atan(YD./(Lae-XD)); Q1=diff(Q); Q2=diff(Q1); end5．计算结果t 时间 XD D 点横坐标YD D 点纵坐标Q构件5角度Q1 构件5角速度 Q2构件5角加速度5000.10050.003008D点运动轨迹：构件5的角位移线图：构件5的角速度图：构件5的角加速度线图：6．计算结果分析原动件1可做整周回转，同时构件5也可做整周运动。
机械原理大作业连杆机构22号
构件5角位移随时间变化先变大后变小，当t=0.26s左右时构件5角位移达到最大，为4.4rad左右；当t=5.5左右时构件5角位移达到最小为2.3rad左右。
构件5角速度先增大后减小再增大，当t=0.1s左右时达最大角速度为10rad/s左右，当t=0.48s左右时，构件5角速度为-12rad/s左右。
fEGF=acos((lEF.^2+lFG.^2-lEG.^2)./(2*lEF.*lFG));
fEGX=atan2((yE-yG),(xE-xG));
fFGX=fEGF+fEGX;
xF=xG-lFG.*cos(fFGX);
yF=yG+lFG.*sin(fFGX);
subplot(2,2,1);
plot(xE,yE);
lBC=670;
lEF=390;
lFG=130;
lCD=350;
xA=0;
yA=0;
xG=465;
yG=-248;
xD=730;
yD=-138;
t=[0:0.000001:0.2*pi];
f1=10.*t;
xB=lAB.*cos(f1);
yB=lAB.*sin(f1);
lBD=sqrt((xB-xD).^2+(yB-yD).^2);
大作业1连杆机构运动分析
1.运动分析题目
（31）如图1-31所示机构，已知机构各构件的尺寸为AB=200mm，BC=670mm，CD=350mm，BE=335mm，EF=380mm，FG=130mm，GH=100mm，HK=486mm，β=124°，xD=730mm，yD=138mm，yG=138mm，h=80mm，构件1的角速度为ω1=10rad/s，试求构件2上点E的轨迹及构件5的角位移、角速度和角加速度，并对计算结果进行分析。
哈工大机械原理大作业连杆机构运动分析完美满分版哈尔滨工业大学
连杆机构运动分析说明书院（系）机电工程学院专业机械设计制造及其自动化姓名李乾学号1130810904班号1308109指导教师唐德威、赵永强日期2015年6月20日哈尔滨工业大学机电工程学院2015年6月一、题目如图1所示机构，已知机构各构件的尺寸为l AB=200mm，l BD=700mm，l AC=400mm，l AE=800mm，构件1的角速度为ω1=10rad/s，试求构件2上点D的轨迹及构件5的角位移、角速度和角加速度，并对计算结果进行分析。
（题中构件尺寸满足l BD-l AB<l AE<l BD+l AB)。
图 1 机构运动简图二、建立数学模型分析1.建立坐标系建立以点A为原点的平面直角坐标系A-x，y，如图2所示图 2 建立坐标系2.对机构进行结构分析该机构由Ⅰ级机构AB、两个RPRⅡ级基本杆组BCD、ED组成。
杆组拆分结果如图3、图4、图5所示。
图 3 Ⅰ级杆组AB图 4 RPRⅡ级基本杆组BCD图 5 RPRⅡ级基本组DE3.确定已知参数和求解流程（1）原动件AB（I级杆组）已知原动件1的转角φ=0~360°运动副A的运动参数x A=0y A=0原动件AB的长度l AB = 200mm代入I级杆组子程序，得到运动副B的位置坐标（x B，y B）根据《机械原理》第三版书中第36页的公式推导可知：A，B两点坐标在x轴，y轴上投影，得方程x B = x A+l AB*cosφy B = y A+l AB*sinφ（2）BCD（RPR II级杆组）已知运动副B的位置坐标（x B，y B）运动副C的坐标位置：x C=l AC=400mmy C=0代入RPR II级杆组子程序，求出构件2上D点的位置坐标（x D，y D）根据《机械原理》第三版书中第339页的公式推导可知：当杆件处于图所示位置,即x B>x D并且y B≥y D时,l j杆角位移：φj=arctan B0s+A0C0 A0s−B0C0式中：A0=x B-x DB0=y B-y DC0=l i+l ks=√A02+B02−C02而当x B<x D并且y B≥y D时，φj=arctan B0s+A0C0A0s−B0C0+180o 当x B<x D并且y B<y D时，φj=arctan B0s+A0C0A0s−B0C0+180o 当x B>x D并且y B<y D时，φj=arctan B0s+A0C0A0s−B0C0+360o图 6 RPR II级杆组分析内移动副C的位置：x C=x B-l i sinφjy C=y B-l i cosφj导杆上E点的位置：x E=x C+(l j-s)cosφjy E=y C+(l j-s)sinφj（3）DE（RPR II级杆组）已知运动副D的位置坐标（x D，y D），运动副E的坐标：x E=l AE=800mmy E=0代入RPR II级杆组子程序，求出构件5的转角φ5。
哈工大机械原理大作业1
[键入公司名称]机械原理课程设计[键入文档副标题]p[选取日期]设计题目连杆机构运动分析机电工程学院1008103班H100811109学号设计者王鹏[在此处键入文档的摘要。
摘要通常是对文档内容的简短总结。
在此处键入文档的摘要。
摘要通常是对文档内容的简短总结。
]1、运动分析题目(12)如图所示的六杆机构中，各构件的尺寸分别为：l AB=200mm ，l BC=500mm ，l CD=800mm ,x F= 400mm ，x D=350mm ，y D=350mm ，ω=100rad s⁄，求构件5上点F的位移，速度和加速度。
2、建立坐标系建立以点A为原点的固定平面直角系A−x,y3、对机构进行结构分析该机构由I级杆组RR（原动件1）、II级杆组RRR（杆2、杆3）和II级杆组PRP（滑块4及滑块5）组成。
I级杆组RR，如图2所示；II级杆组RRR，如图2所示；II级杆组PRP，如图3所示。
4、确定已知参数和求解流程图1所示，规定当φ=10° 时，F 点纵坐标为0 （1）如图2所示，已知原动件杆1的转角，φ=0~360° {x B =l AB ×cos φy B =l AB ×sin φ（2）如图3所示，已知B ，D 两点坐标分别为（x B ,y B ）（x D ,y D ）和 l BC l CD利用方程组{（x −x B ）2+(y −y B )2=l BC 2(x −x D )2+(y −y D )2=l CD2 可以求解出C 点坐标（3）如图4所示，已知C 点坐标、x F 、x D 、y D利用几何关系可以求解出E 点坐标同时，当φ=10° 时，可以求出杆EF 长，记为 l EF进而，可以求出F 点坐标，即F 点位移（4）利用导数的定义与其物理意义v n =s n+h −s n−h2ha n =s n+h −2s n +s n−hh 2利用上述公式，选取适当的步长h ，利用F 点位移就可以得出速度与加速度5、用VC 编程#include <stdio.h>#include <math.h>#define pi 3.14159265358979323846//定义全局变量double Lab,Lbc,Lcd,Xf,Xd,Yd;//定义已知位置量double Wab;//定义角速度量//定义被调用函数void Ccorner (double *a,double *b,double c);//声明C点坐标函数求解函数void Pcorner (double *a,double *b,double c,double d,double e);//声明P点坐标函数求解函数double Kcd (double a,double b);//声明CD直线倾斜角求解函数//主函数main (){Lab=0.200;Lbc=0.500;Lcd=0.800;Xf=0.400;Xd=0.350;Yd=0.350;//赋位置量值Wab=100;//赋角速度值//未知几何与位置参量double Xb,Yb;//定义B点坐标double Xc,Yc;//定义C点坐标double Ye;//定义E点纵坐标double Xp,Yp;//定义瞬心p点坐标double Lef;//定义bp,cp,ef,bd杆长double Yf[720];//定义F点纵坐标//未知速度参量double Vf[720];//定义EF杆速度//未知加速度参量double Af[720];//定义加速度//其余参量double o=10*pi/180,k,k1;//主动杆角度变量与CD杆倾斜角double t=1*pi/180/100;//时间参量，用定义法求速度与加速度int i;//循环控制变量//主函数主体//求位移量for (i=0;i<=361;i++){//准备几何量Xb=Lab*cos(o), Yb=Lab*sin(o);Ccorner (&Xc,&Yc,o);//求C点坐标Pcorner (&Xp,&Yp,Xc,Yc,o);//求瞬心P点坐标k=Kcd (Xp,Yp);//求CD杆倾斜角//求解位移量（规定主动杆10度为Yf零点）if (i==0){k1=k;}Lef=tan(k1)*Xd+Yd;//ef杆长Ye=tan(k)*Xd+Yd;Yf[i]=Ye-Lef;o=o+1*pi/180;}//用定义求速度for (i=1;i<=361;i++){Vf[i]=(Yf[i+1]-Yf[i-1])/(2*t);}//用定义求加速度for (i=1;i<=361;i++){Af[i]=(Yf[i+1]-2*Yf[i]+Yf[i-1])/pow(t,2);}//输出语句for (i=1;i<=180;i++){printf ("%d,%lf,%lf,%lf,\t",i+10,Yf[i]*1000,Vf[i],Af[i]);printf ("%d,%lf,%lf,%lf\n",i+190,Yf[i+180]*1000,Vf[i+180],Af[i+180]); }}//C点坐标值函数void Ccorner(double *a,double *b,double c)//&Xc,&Yc,o{double i,j,x,y,z;//中间参数i=(pow(Lcd,2)-pow(Lbc,2)+pow(Lab*sin(c),2)-pow(Yd,2)+pow(Lab*cos(c),2)-pow((Xd+Xf),2))/(2*(Lab*cos(c)-Xd-Xf));j=(Lab*sin(c)-Yd)/(Lab*cos(c)-Xf-Xd);x=pow(j,2)+1;y=2*(Xf+Xd)*j-2*i*j-2*Yd;z=pow(i,2)-2*(Xf+Xd)*i+pow((Xf+Xd),2)+pow(Yd,2)-pow(Lcd,2);*b=(-y+sqrt(pow(y,2)-4*x*z))/(2*x);*a=i-j*(*b);}//求瞬心P点坐标函数void Pcorner(double *a,double *b,double c,double d,double e)//&Xp,&Yp,Xc,Xp,o{double a1,b1,c1,a2,b2,c2,i,j,k;//中间参量a1=tan(e),a2=(d-Yd)/(c-Xd-Xf);b1=b2=-1;c1=0,c2=(c*(d-Yd)/(c-Xd-Xf))-d;i=a1*b2-b1*a2;j=c1*b2-b1*c2;k=a1*c2-c1*a2;*a=j/i;*b=k/i;}//CD杆倾斜角double Kcd(double a,double b)//Xp,Yp{double k1;//中间参量k1=-atan ((b-Yd)/(a-Xd-Xf));return k1;}6、计算结果（程序计算结果附在图像之后）6.1位移、速度、加速度的图像点F的位移线图如图5所示。
机械原理大作业1连杆机构27题
大作业1 连杆机构运动分析1、运动分析题目如图所示机构，已知机构各构件的尺寸为280mm AB =,350mm BC =,320mm CD =,160mm AD =，175mm BE = 220mm EF =，25mm G x =，80mm G y =，构件1的角速度为110rad/s ω=，试求构件2上点F 的轨迹及构件5的角位移、角速度和角加速度，并对计算结果进行分析。
2、建立坐标系建立以点A 为原点的固定平面直角坐标系图13、对机构进行结构分析该机构由I级杆组RR（原动件1）、II级杆组RRR（杆2、杆3）和II级杆组RPR（滑块4及杆5）组成。
I级杆组RR，如图2所示；II级杆组RRR，如图3所示；II级杆组RPR，如图4所示。
图2图 3图 44、各基本杆组运动分析的数学模型（1）同一构件上点的运动分析：图 5如图5所示的构件AB,,已知杆AB 的角速度=10/rad s ω，AB 杆长i l =280mm,可求得B 点的位置B x 、B y ，速度xB v 、yB v ，加速度xB a 、yB a 。
=cos =280cos B i x l ϕϕ;=sin =280sin B i y l ϕϕ;==-sin =-BxB i B dx v l y dt ωϕω； ==cos =;B yB i B dyv l x dt ωϕω222B 2==-cos =-BxB i d x a l x dt ωϕω；2222==-sin =-ByB i B d y a l y dtωϕω。
（2）RRRII 级杆组的运动分析：图 6如图6所示是由三个回转副和两个构件组成的II 级组。
已知两杆的杆长2l 、3l 和两个外运动副B 、D 的位置（B x 、B y 、D x 、D y ）、速度（xB yB xD yD v v v v 、、、）和加速度（xB yB xD yD a a a a 、、、）。
求内运动副C 的位置（C C x 、y ）、速度（xC yC v 、v ）、加速度（xC yC a 、a ）以及两杆的角位置（23ϕϕ、）、角速度（23ϕϕ、）和角加速度（23ϕϕ、）。
哈工大机械原理大作业——连杆机构运动分析报告162014
课程名称： 机械原理
院
系：
能源科学与工程学院
班
级：
完成
者：
学
号：
题
号：
16
任课教师： 丁刚
完成内容： 在完成题目计算要求的同时， 扩展了内容， 程序为该结构的通用程 序，可解决机构在不同条件下的运动情况，文本最末为几种情况的分析
哈尔滨工业大学
16、如图所示机构，已知机构各构件的尺寸为
77
1
120]
[ 117
78
1
123]
[ 119
79
1
127]
[ 120
80
1
129]
[ 121
80
1
131]
[ 121
81
1
132]
[ 123
82
1
132]
[ 124
83
1
132]
[ 125
84
1
132]
[ 125
85
1
131]
[ 127
86
1
130]
[ 128
87
1
128]
[ 129
88
0
126]
这是同一根杆， °，λ，λ
现在来看 ， ， ，由题目得， ， 和 是未知的，但不影响整 体，不然给 和 一个初值 °， ，当然，这是可以随意更改的。
基于以上的基本原理，
syms theta theta1 theta2 lamuda lamuda1 lamuda2 sigma sigma1 sigma2 beta beta1 beta2 l1 l11 l2 l21 t output i
哈工大机械原理大作业-连杆
Harbin Institute of Technology机械原理大作业一课程名称：机械原理设计题目：连杆运动分析院系：机电工程学院班级：设计者：学号：指导教师：明设计时间： 2013年6月25日1、运动分析题目在图1-10中所示的干草压缩机中，已知LAB=150mm,LBC=600mm,LCE=120mm,LCD=500mm,LEF=600mm,XD=400mm,YD=500mm,YF=600mm,曲柄1作等速转动，其转速n1=50r/min。
求在一个运动循环中活塞5的位移、速度和加速度的变化曲线。
图1-102、机构的结构分析（1）基本杆组的划分①AB即杆件1为原动件②DECB即杆件2、3为RRR型II级杆组，其中CE为同一构件上点。
③ EF 和滑块即4、5为RRP 型II 级杆组（2）、建立以点A 为原点的固定平面直角系3、确定已知参数和求解流程(1)原动件1（I 级杆组RR ）如图所示，已知原动件1的转角πϕ2~01=原动件杆1的角速度s rad /236.51=ω原动件1的角加速度01=α运动副A 的位置坐标0,0==A A y x运动副A 的速度0,0==yA xA v v运动副A 的加速度0,0==yA xA a a原动件杆I 的长度mm l 1501=可求出B 的位置B 的速度B 的加速度(2)构件2、3（II 级杆组RRR ）D 的位置500400==D D y xD 的速度 00==yD xD v vD 的加速度00==yD xD a a杆长mm l l CD j 500==，mm l l BC i 600==由关系j j D i i B C l x l x x ϕϕcos cos +=+=j j D i i B C l y l y y ϕϕsin sin +=+=其中)(20sin cos B D i i i x x l A C B A -==-+ ϕϕ)(2B D i y y l B -=222j BD i l l l C -+=222)()(B D B D BD y y X x l -+-=可解得C A C B A B i +-++=222arctan 2ϕ DC D C j x x y y --=arctan ϕ 由上面两个式子可以得到两杆的角速度 1)()(G y y s x x c w B D j B D j i i -+-==•ϕ132G s G c G ij i i +==ϕα其中i j j i s c s c G -=1，i i i l c ϕcos =，i i i l s ϕsin =，j j j l c ϕcos =，j j j l s ϕsin =可得E 的位置iB E iB E y y x x ϕϕsin 480cos 480+=+=E 的速度 i i yB E yE i i xB E xE v yv v x v ωϕωϕcos 480sin 480+==-==E 的加速度i i i i xB E xE a xa αϕωϕsin 480cos 4802--== i i i i yB E yE a ya ϕαϕωcos 480sin 4802+-== （3)、构件4、5杆组（II 级杆组RRP ）在建立的坐标系中取一参考点K 600,0==K K y x则速度0,0==yK xK v v加速度0,0==yK xK a a杆长mm l EF 600=,设F 位移为s由s x l x x K i EF E F +=+='cos ϕ K i EF E F y l y y =+='sin ϕ由上面两个式子可以得到 600600arcsin arcsin'E EF E k i y l y y -=-=ϕ所以：F 点位移 )600600arcsin(cos 600E E F y x x s -+== 速度F F xv = 加速度F F xa =四、编程计算并输出结果（VB 编程）主程序：Private Sub Command1_Click()Dim s5(3600) As DoubleDim v5(3600) As DoubleDim a5(3600) As DoubleDim pi As DoubleDim pa As Doublepi = 3.1415926pa = pi / 180Dim i As LongDim f1(3600) As Double Dim RR1 As RRDim RR2 As RRDim RRR1 As RRRDim RRP1 As RRPSet RR1 = New RRSet RR2 = New RRSet RRR1 = New RRRSet RRP1 = New RRPFor i = 0 To 3600 Step 1 f1(i) = i * pa / 10RR1.delt = 0RR1.f = f1(i)RR1.w = 5.24RR1.e = 0RR1.L = 150RR1.xA = 0RR1.yA = 0RR1.vxA = 0RR1.vyA = 0RR1.axA = 0RR1.ayA = 0RR1.calRRR1.Li = 600RRR1.Lj = 500RRR1.xB = RR1.xBRRR1.yB = RR1.yBRRR1.vxB = RR1.vxBRRR1.vyB = RR1.vyBRRR1.axB = RR1.axBRRR1.ayB = RR1.ayBRRR1.xD = 400RRR1.yD = 500RRR1.vxD = 0RRR1.vyD = 0RRR1.axD = 0RRR1.ayD = 0RRR1.M = 1RRR1.calRRRRR2.delt = 0RR2.f = RRR1.fi RR2.w = RRR1.wi RR2.e = RRR1.ei RR2.L = 480RR2.xA = RR1.xB RR2.yA = RR1.yB RR2.vxA = RR1.vxB RR2.vyA = RR1.vyB RR2.axA = RR1.axB RR2.ayA = RR1.ayB RR2.calRRP1.Li = 600RRP1.Lj = 0RRP1.fj = piRRP1.wj = 0RRP1.ej = 0RRP1.xB = RR2.xB RRP1.yB = RR2.yB RRP1.vxB = RR2.vxB RRP1.vyB = RR2.vyB RRP1.axB = RR2.axB RRP1.ayB = RR2.ayB RRP1.xK = 0RRP1.yK = 600RRP1.vxK = 0RRP1.vyK = 0RRP1.axK = 0RRP1.ayK = 0RRP1.M = 1RRP1.cals5(i) = RRP1.ssv5(i) = RRP1.vssa5(i) = RRP1.assNext iPicture1.Scale (-30, 700)-(360, 580)Picture1.Line (0, 0)-(360, 0) 'XPicture1.Line (0, 580)-(0, 700) 'YFor i = 0 To 360 Step 10 'X轴坐标Picture1.DrawStyle = 2Picture1.Line (i, 700)-(i, 580)Picture1.CurrentX = i - 10: Picture1.CurrentY = 0 Picture1.Print iNext iFor i = 580 To 700 Step 10 'Y轴坐标Picture1.DrawStyle = 2Picture1.Line (0, i)-(360, i)Picture1.CurrentX = -10: Picture1.CurrentY = iPicture1.Print iNext iFor i = 0 To 3600 Step 1Picture1.PSet (i / 10, s5(i))Next iEnd SubRR:Public L As DoublePublic f As DoublePublic delt As DoublePublic w As DoublePublic e As DoublePublic xA As DoublePublic yA As DoublePublic vxA As DoublePublic vyA As DoublePublic axA As DoublePublic ayA As DoublePublic xB As DoublePublic yB As DoublePublic vxB As DoublePublic vyB As DoublePublic axB As DoublePublic ayB As DoublePublic Sub cal()xB = xA + L * Cos(f + delt)yB = yA + L * Sin(f + delt)vxB = vxA - w * L * Sin(f + delt)vyB = vyA + w * L * Cos(f + delt)axB = axA - w ^ 2 * L * Cos(f + delt) - e * L * Sin(f + delt) ayB = ayA - w ^ 2 * L * Sin(f + delt) + e * L * Cos(f + delt) End SubRRR:Public Li As DoublePublic Lj As DoublePublic fi As DoublePublic fj As DoublePublic wi As DoublePublic wj As DoublePublic ei As DoublePublic ej As DoublePublic xB As DoublePublic yB As DoublePublic vxB As DoublePublic vyB As DoublePublic axB As DoublePublic ayB As DoublePublic xC As DoublePublic yC As DoublePublic vxC As DoublePublic vyC As DoublePublic axC As DoublePublic ayC As DoublePublic xD As DoublePublic yD As DoublePublic vxD As DoublePublic vyD As DoublePublic axD As DoublePublic ayD As DoublePublic M As DoublePublic Sub calRRR()Dim fDB As DoubleDim Ci As DoubleDim Cj As DoubleDim Si As DoubleDim Sj As DoubleDim G1 As DoubleDim G2 As DoubleDim G3 As DoubleDim LBD As DoubleDim JCBD As DoubleDim val As Doublepi = 3.1415926LBD = Sqr((xB - xD) ^ 2 + (yD - yB) ^ 2)If LBD < Li + Lj And LBD > Abs(Li - Lj) Thenval = (Li ^ 2 + LBD ^ 2 - Lj ^ 2) / (2 * Li * LBD) JCBD = Atn(-val / Sqr(-val * val + 1)) + 2 * Atn(1) End IfRRP:Public Li As DoublePublic Lj As DoublePublic fi As DoublePublic fj As DoublePublic wi As DoublePublic wj As DoublePublic ei As DoublePublic ej As DoublePublic xB As DoublePublic yB As DoublePublic vxB As DoublePublic vyB As DoublePublic axB As DoublePublic ayB As DoublePublic xK As DoublePublic yK As DoublePublic vxK As DoublePublic vyK As DoublePublic axK As DoublePublic ayK As DoublePublic xC As DoublePublic yC As DoublePublic vxC As DoublePublic vyC As DoublePublic axC As DoublePublic ayC As DoublePublic xD As DoublePublic yD As DoublePublic vxD As DoublePublic vyD As DoublePublic axD As DoublePublic ayD As DoublePublic M As SinglePublic ss As DoublePublic vss As DoublePublic ass As DoublePublic Sub cal()Dim A0 As DoubleDim Q1 As DoubleDim Q2 As DoubleDim Q3 As DoubleDim Q4 As DoubleDim Q5 As DoubleDim val As DoubleDim pi As Doublepi = 3.14159216A0 = Lj + ((yK - yB) * Cos(fj) - (xK - xB) * Sin(fj)) val = A0 / Lifi = M * Atn(val / Sqr(-val * val + 1)) + fjxC = xB + Li * Cos(fi)yC = yB + Li * Sin(fi)ss = (xC - xK) * Cos(fj) + (yC - yK) * Sin(fj)xD = xK + ss * Cos(fj)yD = yK + ss * Sin(fj)Q1 = vxK - vxB - wj * (ss * Sin(fj) + Lj * Cos(fj))Q2 = vyK - vyB + wj * (ss * Cos(fj) - Lj * Sin(fj))Q3 = Li * Sin(fi) * Sin(fj) + Li * Cos(fi) * Cos(fj) wi = (-Q1 * Sin(fj) + Q2 * Cos(fj)) / Q3vss = -(Q1 * Li * Cos(fi) + Q2 * Li * Sin(fi)) / Q3vxC = vxB - wi * Li * Sin(fi)vyC = vyB + wi * Li * Cos(fi)vxD = vxK + vss * Cos(fj) - ss * wj * Sin(fj)vyD = vyK + vss * Sin(fj) + ss * wj * Cos(fj)Q4 = axK - axB + wi ^ 2 * Li * Cos(fi) - ej * (ss * Sin(fj) + Lj * Cos(fj)) - wj ^ 2 * (ss * Cos(fj) - Lj * Sin(fj)) - 2 * vss * wj * Sin(fj)Q5 = ayK - ayB + wi ^ 2 * Li * Sin(fi) + ej * (ss * Cos(fj) - Lj * Sin(fj)) - wj ^ 2 * (ss * Sin(fj) + Lj * Cos(fj)) + 2 * vss * wj * Cos(fj)ei = (-Q4 * Sin(fj) + Q5 * Cos(fj)) / Q3ass = (-Q4 * Li * Cos(fi) - Q5 * Li * Sin(fi)) / Q3axC = axB - ei * Li * Sin(fi) - wi ^ 2 * Li * Cos(fi)ayC = ayB + ei * Li * Cos(fi) - wi ^ 2 * Li * Sin(fi)axD = axK + ass * Cos(fj) - ss * ej * Sin(fj) - ss * wj ^ 2 * Cos(fj) - 2 * vss * wj * Sin(fj) ayD = ayK + ass * Sin(fj) + ss * ej * Cos(fj) - ss * wj ^ 2 * Sin(fj) + 2 * vss * wj * Cos(fj) End Sub五、计算结果数据如图：位移曲线：速度曲线：加速度曲线：六、计算结果分析主动件转角为0时,滑块的位移为628mm,随着转角的匀速增加,滑块位移先上升,速度为负快速下降,加速度为负且开始值较小并逐渐下降,到达最低点-2200左右,此时速度为0, 然后位移开始继续下降,速度继续下降,然后到达最小值-100左右,此时加速度为0,又开始上升(向左运动),速度正向增大,在后面一段时期速度继续增大,加速度也正向增大,然后速度下降,加速度下降.滑块分别在4.2s,11.4s速度达到正向最大，1.7s,7.7s速度达到负向最大，0.8s,2.7s,5.8s,9.6s加速度达到极值，滑块就是这样周期性的左右运动。
哈尔滨工业大学机械原理大作业_连杆26题
1.运动分析题目如图所示机构，已知机构各构件的尺寸为AB=108mm，EF=320mm，BC=CE=CD=200mm，FG=162mm，AD=258mm，AG=514mm，DG=384mm，β=80º，构件1的角速度为ω1=10rad/s，试求构件2上点E的轨迹及构件5的角位移、角速度和角加速度，并对计算结果进行分析。
2.机构的结构分析，组成机构的基本杆组划分3.各基本杆组的运动分析数学模型（1）RR基本杆组：delt=0xB = xA + AB * Cos(f + delt)yB = yA + AB* Sin(f + delt)vxB = vxA - w * AB * Sin(f + delt)vyB = vyA + w * AB* Cos(f + delt)axB = axA - w ^ 2 * AB * Cos(f + delt):ayB = ayA - w ^ 2 * AB* Sin(f + delt)(2)RRR基本杆组Ci = lbc * Cos(fbc)Si = lbc * Sin(fbc)Cj = lcd * Cos(fcd)Sj = lcd * Sin(fcd)G1 = Ci * Sj - Cj * Siwbc = (Cj * (vxD - vxB) + Sj * (vyD - vyB)) / G1wcd = (Ci * (vxD - vxB) + Si * (vyD - vyB)) / G1vxC = vxB - wbc * lbc * Sin(fbc)vyC = vyB + wbc * lbc * Cos(fbc)G2 = axD - axB + wbc ^ 2 * Ci - wcd ^ 2 * CjG3 = ayD - ayB + wbc ^ 2 * Si - wcd ^ 2 * Sjebc = (G2 * Cj + G3 * Sj) / G1ecd = (G2 * Ci + G3 * Si) / G1axC = axB - ebc * lbc * Sin(fbc) - wbc ^ 2 * lbc * Cos(fbc)ayC = ayB + ebc * lbc * Cos(fbc) - wbc ^ 2 * lbc * Sin(fbc)EB = 2 * BC * Cos(febc)xE = xB + EB* Cos(fbc + febc)yE = yB + EB * Sin(fbc + febc)vxE = vxB – wbc * EB* Sin(fbc + febc)vyE = vyB + wbc * EB* Cos(fbc + febc)axE = axB - wbc ^ 2 * EB * Cos(fbc + delt) - ebc * EB * Sin(fbc + febc)ayE = ayB - wbc ^ 2 * leb * Sin(fbc + delt) + ebc * leb * Sin(fbc + febc) Ci = lef * Cos(fef)Si = lef * Sin(fef)Cj = lfg * Cos(ffg)Sj = lfg * Sin(ffg)G1 = Ci * Sj - Cj * Siwef = (Cj * (vxG - vxE) + Sj * (vyG - vyE)) / G1wfg = (Ci * (vxG - vxE) + Si * (vyG - vyE)) / G1vxF = vxE - wef * lef * Sin(fef)vyF = vyE + wef * lef * Cos(fef)G2 = axG - axE + wef ^ 2 * Ci - wfg ^ 2 * CjG3 = ayG - ayE + wef ^ 2 * Si - wfg ^ 2 * Sjeef = (G2 * Cj + G3 * Sj) / G1efg = (G2 * Ci + G3 * Si) / G1axF = axE - eef * lef * Sin(fef) - wef ^ 2 * lef * Cos(fef)ayF = ayE + eef * lef * Cos(fef) - wef ^ 2 * lef * Sin(fef)4.计算编程Dim xA As DoubleDim yA As DoubleDim vxA As DoubleDim vyA As DoubleDim axA As Double 'A '点加速度x轴分量Dim ayA As Double 'A '点加速度y轴分量Dim xB As Double 'B'点'x轴坐标Dim yB As Double 'B点y轴坐标Dim vxB As Double 'B点速度x轴分量Dim vyB As Double 'B点速度y轴分量Dim axB As Double 'B点加速度x轴分量Dim ayB As Double 'B点加速度y轴分量Dim xC As Double 'C点x轴坐标Dim yC As Double C'点y轴坐标Dim vxC As Double 'C点速度x轴分量Dim vyC As Double 'C点速度y轴分量Dim axC As Double 'C点加速度x轴分量Dim ayC As Double 'C点加速度y轴分量Dim xD As Double 'D点x轴坐标Dim yD As Double 'D点y轴坐标Dim vxD As Double 'D点速度x轴分量Dim vyD As Double 'D点速度y轴分量Dim axD As Double 'D点加速度x轴分量Dim ayD As Double 'D点加速度y轴分量Dim xE As Double 'E点x轴坐标Dim yE As Double 'E点y轴坐标Dim vxE As Double 'E点速度x轴分量Dim vyE As Double 'E点速度y轴分量Dim axE As Double 'E点加速度x轴分量Dim ayE As Double 'E点加速度y轴分量Dim xF As Double 'F点x轴坐标Dim yF As Double 'F点y轴坐标Dim vxF As Double 'F点速度x轴分量Dim vyF As Double 'F点速度y轴分量Dim axF As Double 'F点加速度x轴分量Dim ayF As Double 'F点加速度y轴分量Dim xG As Double 'G点x轴坐标Dim yG As Double 'G点y轴坐标Dim vxG As Double 'G点速度x轴分量Dim vyG As Double 'G点速度y轴分量Dim axG As Double 'G点加速度x轴分量Dim ayG As Double 'G点加速度y轴分量Dim delt As Double ' AB杆初始转角Dim lab As Double 'AB杆长Dim lbc As Double 'BC杆长Dim lcd As Double ' CD杆长Dim lce As Double 'CE杆长Dim lef As Double 'EF杆长Dim lfg As Double 'FG杆长Dim leb As Double 'ED杆长Dim f As Double 'AB杆转角Dim fbc As Double 'BC杆转角Dim fcd As Double 'CD杆转角Dim fce As Double 'CE杆转角Dim fef As Double 'EF杆转角Dim ffg As Double 'FG杆转角Dim fge As Double 'ge杆转角Dim w As Double 'AB杆角速度Dim wbc As Double ' BC角速度Dim wcd As Double 'CD角速度Dim wce As Double 'CE角速度Dim wef As Double 'EF角速度Dim wfg As Double 'FG角速度Dim e As Double 'AB杆角加速度Dim ebc As Double ' BC杆角加速度Dim ecd As Double 'CD杆角加速度Dim ece As Double 'CE杆角加速度Dim eef As Double 'EF杆角加速度Dim efg As Double 'FG杆角加速度Dim LBD As Double 'BD距离Dim leg As Double 'EG距离Dim JCBD As Double '角CBDDim jfeg As Double '角FEGDim fBD As Double 'BD转角Dim feg As Double 'EG转角Dim Ci As DoubleDim Cj As DoubleDim Si As DoubleDim Sj As DoubleDim G1 As DoubleDim G2 As DoubleDim G3 As DoubleDim val As DoubleDim pi As DoubleDim pa As DoubleDim febc As Double '角EBCDim i As DoubleDim fj1 As DoublePrivate Sub Command1_Click() '求点E的轨迹Picture1.Scale (-300, 400)-(10, -15)Picture1.Line (-300, 0)-(10, 0) 'XPicture1.Line (0, 400)-(0, -15) 'YFor i = -300 To 0 Step 50 'X轴坐标Picture1.DrawStyle = 2Picture1.Line (i, 400)-(i, 0)Picture1.CurrentX = i - 10: Picture1.CurrentY = 0 Picture1.Print iNext iFor i = 0 To 350 Step 50 'Y轴坐标Picture1.DrawStyle = 2Picture1.Line (0, i)-(-400, i)Picture1.CurrentX = -20: Picture1.CurrentY = i + 7 Picture1.Print iNext iFor fj1 = 0 To 360 Step 0.01f = fj1 * paCall RR1Call RRR1Call RR2Picture1.PSet (xE, yE)Next fj1End SubPrivate Sub Command2_Click() '求构件5的角位移Picture2.Scale (-20, 5)-(380, -0.5)Picture2.Line (-20, 0)-(380, 0) 'XPicture2.Line (0, 3)-(0, -0.5) 'YFor i = 0 To 360 Step 30 'X轴坐标Picture2.DrawStyle = 2Picture2.Line (i, 3)-(i, 0)Picture2.CurrentX = i - 10: Picture2.CurrentY = 0 Picture2.Print iNext iFor i = -0.5 To 3 Step 0.5 'Y轴坐标Picture2.Line (0, i)-(380, i)Picture2.CurrentX = -25: Picture2.CurrentY = i Picture2.Print iNext iFor fj1 = 0 To 360 Step 0.01f = fj1 * paCall RR1Call RRR1Call RR2Call RRR2Picture2.PSet (fj1, ffg)Next fj1End SubPrivate Sub Command3_Click() '求构件5的角速度Picture3.Scale (-20, 10)-(380, -10)Picture3.Line (-20, 0)-(380, 0) 'XPicture3.Line (0, 10)-(0, -10) 'YFor i = 0 To 360 Step 30 'X轴坐标Picture3.DrawStyle = 2Picture3.Line (i, 10)-(i, -10)Picture3.CurrentX = i - 10: Picture3.CurrentY = 0 Picture3.Print iNext iFor i = -8 To 8 Step 2 'Y轴坐标Picture3.Line (0, i)-(380, i)Picture3.CurrentX = -20: Picture3.CurrentY = i Picture3.Print iNext iFor fj1 = 0 To 360 Step 0.01f = fj1 * paCall RR1Call RRR1Call RR2Call RRR2Picture3.PSet (fj1, wfg)Next fj1End SubPrivate Sub Command4_Click() '求构件5的角加速度Picture4.Scale (-20, 300)-(380, -200)Picture4.Line (-20, 0)-(380, 0) 'XPicture4.Line (0, 300)-(0, -200) 'YFor i = 0 To 360 Step 30 'X轴坐标Picture4.DrawStyle = 2Picture4.Line (i, 300)-(i, -200)Picture4.CurrentX = i - 10: Picture4.CurrentY = 0 Picture4.Print iNext iFor i = -200 To 300 Step 50 'Y轴坐标Picture4.Line (0, i)-(380, i)Picture4.CurrentX = -25: Picture4.CurrentY = i + 5 Picture4.Print iNext iFor fj1 = 0 To 360 Step 0.01f = fj1 * paCall RR1Call RRR1Call RR2Call RRR2Picture4.PSet (fj1, efgNext fj1End SubPrivate Sub Form_Load() '赋初值lab = 108lce = 200lbc = 200lcd = 200lef = 320lfg = 162w = 10e = 0delt = 0xA = 0yA = 0vyA = 0axA = 0ayA = 0xD = -178.311284yD = 186.464704vxD = 0vyD = 0axD = 0ayD = 0xG = -514yG = 0vxG = 0vyG = 0axG = 0ayG = 0pi = 3.1415926pa = pi / 180febc = pa * 50End SubPrivate Sub RR1() 'RR基本杆组xB = xA + lab * Cos(f + delt)yB = yA + lab * Sin(f + delt)vxB = vxA - w * lab * Sin(f + delt)vyB = vyA + w * lab * Cos(f + delt)axB = axA - w ^ 2 * lab * Cos(f + delt) - e * lab * Sin(f + delt)ayB = ayA - w ^ 2 * lab * Sin(f + delt) + e * lab * Sin(f + delt)End SubPrivate Sub RR2() 'RR基本杆组leb = 2 * lbc * Cos(febc)xE = xB + leb * Cos(fbc + febc)yE = yB + leb * Sin(fbc + febc)vxE = vxB - wbc * leb * Sin(fbc + febc)vyE = vyB + wbc * leb * Cos(fbc + febc)axE = axB - wbc ^ 2 * leb * Cos(fbc + delt) - ebc * leb * Sin(fbc + febc) ayE = ayB - wbc ^ 2 * leb * Sin(fbc + delt) + ebc * leb * Sin(fbc + febc) End SubPrivate Sub RRR1() 'RRR基本杆组LBD = Sqr((xD - xB) ^ 2 + (yD - yB) ^ 2)If LBD > lbc + lcd And LBD < Abs(lbc - lcd) ThenIf MsgBox("RRR杆组杆长不符合要求", vbOKOnly, "提示") = 1 Then EndEnd IfElseEnd IfIf LBD < lbc + lcd And LBD > Abs(lbc - lcd) Then val = (lbc ^ 2 + LBD ^ 2 - lcd ^ 2) / (2 * lbc * LBD) JCBD = Atn(-val / Sqr(-val * val + 1)) + 2 * Atn(1) ElseEnd IfIf LBD = lbc + lcd ThenJCBD = 0ElseEnd IfIf LBD = Abs(lbc - lcd) ThenIf lbc > lcd ThenJCBD = 0ElseEnd IfIf lbc < lcd ThenJCBD = piElseEnd IfElseEnd IfIf xD > xB And yD >= yB Then '第一象限fBD = Atn((yD - yB) / (xD - xB))ElseEnd IfIf xD = xB And yD > yB ThenfBD = pi / 2ElseEnd IfIf xD < xB And yD >= yB Then '第二象限fBD = pi + Atn((yD - yB) / (xD - xB))ElseEnd IfIf xD < xB And yD < yB Then '第三象限fBD = pi + Atn((yD - yB) / (xD - xB))ElseEnd IfIf xD = xB And yD < yB ThenfBD = 3 * pi / 2ElseEnd IfIf xD > xB And yD <= yB Then '第四象限fBD = 2 * pi + Atn((yD - yB) / (xD - xB))ElseEnd Iffbc = fBD - JCBDxC = xB + lbc * Cos(fbc)yC = yB + lbc * Sin(fbc)If xC > xD And yC >= yD Then '第一象限fcd = Atn((yC - yD) / (xC - xD))ElseEnd IfIf xC = xD And yC >= yD Thenfcd = pi / 2ElseEnd IfIf xC < xD And yC >= yD Then '第二象限fcd = pi + Atn((yC - yD) / (xC - xD))ElseEnd IfIf xC < xD And yC < yD Then '第三象限fcd = pi + Atn((yC - yD) / (xC - xD))ElseEnd IfIf xC = xD And yC < yD Thenfcd = 3 * pi / 2ElseEnd IfIf xC > xD And yC <= yD Then '第四象限fcd = 2 * pi + Atn((yC - yD) / (xC - xD))ElseEnd IfCi = lbc * Cos(fbc)Si = lbc * Sin(fbc)Cj = lcd * Cos(fcd)Sj = lcd * Sin(fcd)G1 = Ci * Sj - Cj * Siwbc = (Cj * (vxD - vxB) + Sj * (vyD - vyB)) / G1 wcd = (Ci * (vxD - vxB) + Si * (vyD - vyB)) / G1 vxC = vxB - wbc * lbc * Sin(fbc)vyC = vyB + wbc * lbc * Cos(fbc)G2 = axD - axB + wbc ^ 2 * Ci - wcd ^ 2 * CjG3 = ayD - ayB + wbc ^ 2 * Si - wcd ^ 2 * Sj ebc = (G2 * Cj + G3 * Sj) / G1ecd = (G2 * Ci + G3 * Si) / G1axC = axB - ebc * lbc * Sin(fbc) - wbc ^ 2 * lbc * Cos(fbc)ayC = ayB + ebc * lbc * Cos(fbc) - wbc ^ 2 * lbc * Sin(fbc)End SubPrivate Sub RRR2() 'RRR基本杆组leg = Sqr((xG - xE) ^ 2 + (yG - yE) ^ 2)If leg > lef + lfg And leg < Abs(lef - lfg) ThenIf MsgBox("RRR杆组杆长不符合要求", vbOKOnly, "提示") = 1 Then EndElseEnd IfElseEnd IfIf leg < lef + lfg And leg > Abs(lef - lfg) Thenval = (lef ^ 2 + leg ^ 2 - lfg ^ 2) / (2 * lef * leg)jfeg = Atn(-val / Sqr(-val * val + 1)) + 2 * Atn(1)ElseEnd IfIf leg = lef + lfg Thenjfeg = 0ElseEnd IfIf leg = Abs(lef - lfg) ThenIf lef > lfg Thenjfeg = 0ElseEnd IfIf lef < lfg Thenjfeg = piElseEnd IfElseEnd IfIf xG > xE And yG >= yE Then '第一象限feg = Atn((yG - yE) / (xG - xE))ElseEnd IfIf xG = xE And yG > yE Thenfeg = pi / 2ElseEnd IfIf xG < xE And yG >= yE Then '第二象限feg = pi + Atn((yG - yE) / (xG - xE))ElseEnd IfIf xG < xE And yG < yE Then '第三象限feg = pi + Atn((yG - yE) / (xG - xE)) ElseEnd IfIf xG = xE And yG < yE Thenfeg = 3 * pi / 2ElseEnd IfIf xG > xE And yG <= yE Then '第四象限feg = 2 * pi + Atn((yG - yE) / (xG - xE)) ElseEnd Iffef = feg - jfegxF = xE + lef * Cos(fef)yF = yE + lef * Sin(fef)If xF > xG And yF >= yG Then '第一象限ffg = Atn((yF - yG) / (xF - xG))ElseEnd IfIf xF = xG And yF >= yG Thenffg = pi / 2ElseEnd IfIf xF < xG And yF >= yG Then '第二象限ffg = pi + Atn((yF - yG) / (xF - xG)) ElseEnd IfIf xF < xG And yF < yG Then '第三象限ffg = pi + Atn((yF - yG) / (xF - xG)) ElseEnd IfIf xF = xG And yF < yG Thenffg = 3 * pi / 2ElseEnd IfIf xF > xG And yF <= yG Then '第四象限ffg = 2 * pi + Atn((yF - yG) / (xF - xG)) ElseEnd IfCi = lef * Cos(fef)Si = lef * Sin(fef)Cj = lfg * Cos(ffg)Sj = lfg * Sin(ffg)G1 = Ci * Sj - Cj * Siwef = (Cj * (vxG - vxE) + Sj * (vyG - vyE)) / G1wfg = (Ci * (vxG - vxE) + Si * (vyG - vyE)) / G1vxF = vxE - wef * lef * Sin(fef)vyF = vyE + wef * lef * Cos(fef)G2 = axG - axE + wef ^ 2 * Ci - wfg ^ 2 * CjG3 = ayG - ayE + wef ^ 2 * Si - wfg ^ 2 * Sjeef = (G2 * Cj + G3 * Sj) / G1efg = (G2 * Ci + G3 * Si) / G1axF = axE - eef * lef * Sin(fef) - wef ^ 2 * lef * Cos(fef)ayF = ayE + eef * lef * Cos(fef) - wef ^ 2 * lef * Sin(fef)End Sub5．结果及分析图1 E点的运动轨迹（1）由图1所示，E点的运动轨迹呈稍倾斜“8字形”。
哈工大机械原理大作业(连杆机构)
建立坐标系：以C为原点，水平方向为X轴，CA所在直线为Y轴建立直角坐标系（如图4）。
取曲柄1水平且位于A点右侧为初始时刻，设曲柄1角速度为w，由题意知w= =8.5π rad/s………………（1）
设曲柄1转角为θ，则B点坐标：
xB=ιABcosθ=ιABcoswt
yB=H1+ιABsinθ=H1+ιABsinwt………………（2）
form=1:length(t)-1
ddxF(m)=(dxF(m+1)-dxF(m))/0.0001;
end
ddxF(length(t))=ddxF(length(t)-1);
figure
plot(t,ddxF)
title('¼ÓËÙ¶ÈͼÏñ');
xlabel('t /s'),ylabel('v /(m/s^2)');
输出图像：
xE(m)=yE(m)/k(m);
xF(m)=xE(m)-(-H^2+lEF^2-yE(m)^2+2*yE(m)*H)^(1/2)+0.1142;
end
form=1:length(t)-1
dxF(m)=(xF(m+1)-xF(m))/0.0001;
end
dxF(length(t))=dxF(length(t)-1);
∵ιEF+ιCE>H且ιCE<H
∴E点始终在F点的右下方
∴xF<xE,所以x2舍去，只取xF=x1……………（8）
∴点F坐标为（xF，H）
当t=0时，可得F点初始位置坐标，不妨设为（xo,H）。
则F点位移（通过计算，t=0时，得xo=-0.1142）