哈工大机械原理大作业连杆14-2015版
哈工大机械原理大作业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 之间速度和加速度波动较小，这就保证了刨头在空行程时有急回运动，而在工作行程是由较均匀的切削速度，运动连续，且运动性能良好稳定。
哈工大机械设计大作业平面连杆机构设计说明书
机械设计基础大作业计算说明书题目：平面连杆机构设计学院：材料学院班号：学号：姓名：日期： 2014年9月30日哈尔滨工业大学机械设计基础大作业任务书题目：平面连杆机构设计设计原始数据及要求：()目录1设计题目2设计原始数据3设计计算说明书3.1计算极位夹角3.2设计制图3.3验算最小传动角4参考文献1 设计题目平面连杆机构的图解法设计2 设计原始数据设计一曲柄摇杆机构。
已知摇杆长度，摆角，摇杆的行程速比系数，要求摇杆靠近曲柄回转中心一侧的极限位置与机架间的夹角为，试用图解法设其余三杆的长度，并检验（测量或计算）机构的最小传动角。
()3 设计计算说明书3.1 计算极位夹角极位夹角代入数值3.2 设计制图3.2.1 在图纸上取一点作为点，从点垂直向上引出一条长为的线段，终点为；3.2.2 从点在左侧引出一条与夹角为的射线；3.2.3 以点为圆心，以为半径画圆，与射线交于点；3.2.4 分别从、两点向下引两条射线，射线与夹角为，两射线交于点，点即为曲柄的回转中心；3.2.5 以点为圆心以为半径画圆；3.2.6 过点向左侧引出一条射线，射线与夹角，与圆交于点；3.2.7 连接，并量取其长度，以为半径画圆，直线，与圆的交点分别为，；3.2.8 在图中量取，，3.3 验算最小传动角3.3.1 在处根据余弦定理3.3.2 在处根据余弦定理所以最小传动角4 参考文献[1]宋宝玉，王瑜，张锋主编.机械设计基础.哈尔滨：哈尔滨工业大学出版，2010.[2]王瑜主编.机械设计基础大作业指导书.哈尔滨：哈尔滨工业大学出版社，2014.。
哈工大机械原理大作业
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大作业设计说明书课程名称：机械原理设计题目：凸轮机构院系：机电学院班级：姓名：学号：指导教师：丁刚设计时间：2014.5.29哈尔滨工业大学1.设计题目第31题：升程/mm 升程运动角/。
升程运动规律升程许用压力角/。
回程运动角/。
回程运动规律回程许用压力角/。
远休止角/。
近休止角/。
150 90 等加等减速40 80 余弦加速度70 40 1502.运动方程式及运动线图由题目要求凸轮逆时针旋转（1）确定凸轮机构推杆升程、回程运动方程，并绘制推杆位移、速度、加速度线图。
升程第一段：（0 <φ< pi /4）φ0=pi/2;s1 = 73*φ^2;v1=146*w*φ;a1 = 146*w^2;升程第二段：（pi/4 <φ< pi /2）s2 =90-73*(pi/2-φ)^2;v2=146*w*( pi/2-φ);a2 =-146*w.^2;远休止程：（pi/2 <φ< 10*pi/9)s3 = 90;v3 = 0;a3 = 0;回程：(10*pi/9)< φ< ( 14*pi/9)s4 =45*(1+cos(9/4*(φ-10*pi/9)));v4 =-101.25*w*sin(9/4*(φ-10*pi/9)) ;a4 =-227.8*w^2* cos(9/4*(φ-10*pi/9)); 近休止程：(14*pi/9)< φ < ( 2*pi);s5 =0;v5 =0;a5 =0;1.由上述公式通过编程得到位移、速度、加速度曲线如下：（编程见附录）.基圆半径为r0 = (50^2+100^2)0.5=112mm，偏距e = 50mm。
3.凸轮实际轮廓，理论轮廓，基圆，偏距圆绘制4.整体图像附录1.求位移、速度、加速度的程序（matlab）w = input('请输入W=');x = 0:(pi/1000):(pi/4);s1 = 73*x.^2;v1=146*w*x;a1 = 146*w.^2;y = (pi/4):(pi/1000):(pi/2);s2 =90-73*(pi/2-y).^2;v2=146*w*( pi/2-y);a2 =-146*w.^2;z = (pi/2 ):(pi/1000):(10*pi/9);s3 = 90;v3 = 0;a3 = 0;c = (10*pi/9):(pi/1000):( 14*pi/9);s4 =45*(1+cos(9/4*(c-10*pi/9)));v4 =-101.25*w*sin(9/4*( c-10*pi/9)) ;a4 =-227.8*w.^2* cos(9/4*(c-10*pi/9));d=(14*pi/9):(pi/1000):( 2*pi);s5 =0; v5 =0; a5 =0;subplot (2,2,1)plot(x,s1,'b',y,s2,'b',z,s3,'b',c,s4,'b', d,s5,'b');xlabel('转角/rad')ylabel('位移/(mm/s)')title('位移与转角曲线')grid onsubplot (2,2,2)plot(x,v1,'g',y,v2,'g',z,v3,'g ',c,v4,'g', d,v5,'g')ds4 =45*9/4*sin(9/4*(c-10*pi/9));d=(14*pi/9):(pi/1000):( 2*pi);s5 =0;ds5 =0;plot(ds1,s1,'b',ds2,s2,'b',ds3,s3,'b',ds4,s4,'b',ds5,s5,'b'); xlabel('ds/dp');ylabel('(位移s/mm)')title('ds/dp 与位移s曲线')grid onhold onx3=-150:0.001:0;y3 = 0.577*x3;plot (x3,y3,'r');hold onx1=-150:0.001:150;for i=1:1:250;k1=(s1(i+1)-s1(i))/ (ds1(i+1)-ds1(i));if(k1>=-1.733 && k1<=-1.731)y1=k1*(x1-ds1(i))+s1(i);plot (x1,y1,'r');end3.确定滚子半径（1）先求凸轮理论轮廓曲线，程序如下：Clc；clear;w = input('请输入w=');s0 = 100;s = 90;e = 50; x = 0:(pi/100):(pi/4);x1 = (s + s0)*cos(x)-e*sin(x);y1 = (s0 + s)*sin(x) - e*cos(x);y = (pi/4):(pi/100):(pi/2);x2 = (s + s0)*cos(y)-e*sin(y);y2 = (s0 + s)*sin(y) - e*cos(y);z = (pi/2 ):(pi/100):(10*pi/9);x3 = (s + s0)*cos(z)-e*sin(z);y3 = (s0 + s)*sin(z) - e*cos(z);c = (10*pi/9):(pi/1000):( 14*pi/9);x4 = (s + s0)*cos(c)-e*sin(c);y4 = (s0 + s)*sin(c) - e*cos(c);d=(14*pi/9):(pi/1000):( 2*pi);x5 = (s + s0)*cos(d)-e*sin(d);y5 = (s0 + s)*sin(d) - e*cos(d);plot(x1,y1,'b',x2,y2,'b',x3,y3,'b',x4,y4,'b',x5,y5,'b');xlabel('x/mm')ylabel('y/mm')title('理论轮廓曲线')grid on（2）理论轮廓线最小曲率半径编程代码：由下程序结果可知pmin =81.6667这里滚子半径为 r r < pmin-ΔΔ=3~5mm；取r r=10mm;clear;clc;v=[];syms x1 x2 x3 x4 x5s0 = 100;e = 50;s1 = 73*x1.^2;t1 = (s1 + s0)*cos(x1)-e*sin(x1);y1 = (s0 + s1)*sin(x1) - e*cos(x1);tx1=diff(t1,x1);txx1=diff(t1,x1,2);yx1=diff(y1,x1);yxx1=diff(y1,x1,2);for xx1= 0:(pi/100):(pi/4);k1=subs(abs((tx1*yxx1-txx1*yx1)/(tx1^2+yx1^2)^1.5),{x1},{xx1}); v=[v,1/k1];ends2 =90-73*(pi/2-x2).^2;t2 = (s2 + s0)*cos(x2)-e*sin(x2);y2 = (s0 + s2)*sin(x2) - e*cos(x2);tx2=diff(t2,x2);txx2=diff(t2,x2,2);yx2=diff(y2,x2);yxx2=diff(y2,x2,2);for xx2=(pi/4):(pi/100):(pi/2);k2=subs(abs((tx2*yxx2-txx2*yx2)/(tx2^2+yx2^2)^1.5),{x2},{xx2});k4=subs(abs((tx4*yxx4-txx4*yx4)/(tx4^2+yx4^2)^1.5),{x4},{xx4}); v=[v,1/k4];ends5 =0;t5 = (s5 + s0)*cos(x5)-e*sin(x5);y5 = (s0 + s5)*sin(x5) - e*cos(x5);tx5=diff(t5,x5);txx5=diff(t5,x5,2);yx5=diff(y5,x5);yxx5=diff(y5,x5,2);for xx5=(10*pi/9):(pi/100):( 4*pi/3);k5=subs(abs((tx5*yxx5-txx5*yx5)/(tx5^2+yx5^2)^1.5),{x5},{xx5}); v=[v,1/k5];endpmin=min(v)4.绘制凸轮轮廓曲线clear ;clc;syms x y z c dw= input('请输入w=');n3 = diff(x3);m3 = diff(y3);xt3= subs(x3 + (r*m3)./sqrt(m3.^2+n3.^2),z,zz);yt3 = subs(y3 - (r*n3)./sqrt(m3.^2+n3.^2),z,zz);cc= (10*pi/9):(pi/1000):( 14*pi/9);s4 =45*(1+cos(9/4*(c-10*pi/9)));x4 = (s4 + s0).*cos(c)-e*sin(c);y4 = (s0 +s4).*sin(c) - e*cos(c);n4 = diff(x4);m4 = diff(y4);xt4= subs(x4 + (r*m4)./sqrt(m4.^2+n4.^2),c,cc);yt4 =subs( y4 - (r*n4)./sqrt(m4.^2+n4.^2),c,cc);dd=(14*pi/9):(pi/1000):( 2*pi);s5 =0;x5 = (s5 + s0).*cos(d)-e*sin(d);y5 = (s0 +s5).*sin(d) - e*cos(d);n5 = diff(x5);m5 = diff(y5);xt5= subs(x5 + (r*m5)./sqrt(m5.^2+n5.^2),d,dd);yt5 =subs( y5 - (r*n5)./sqrt(m5.^2+n5.^2),d,dd);plot(xt1,yt1,'b',xt2,yt2,'b',xt3,yt3,'b',xt4,yt4,'b',xt5,yt5,'b')for i=1:3601if yy(1,i)<=y0/2s(1,i)=2*h*(yy(i)./y0).^2;v(1,i)=4*h*w*yy(i)./(y0.^2);a(1,i)=4*h*w.^2./(y0.^2);elseif yy(1,i)>y0/2 && yy(1,i)<y0s(1,i)=h-2*h./y0.^2.*(y0-yy(i)).^2;v(1,i)=4*h*w*(y0-yy(i))./(y0.^2);a(1,i)=-4*h*w.^2./(y0.^2);elseif yy(1,i)>=y0 && yy(1,i)<y0+yss(1,i)=h;v(1,i)=0;a(1,i)=0;elseif yy(1,i)>=y0+ys && yy(1,i)<y0+ys+y01s(1,i)=h/2*(1+cos(pi/y01*(yy(1,i)-y0-ys)));v(1,i)=-pi*h*w/2/y01*sin(pi/y01*(yy(1,i)-y0-ys));a(1,i)=-pi^2*h*w^2/2/y01^2*cos(pi/y01*(yy(1,i)-y0-ys)); elseif yy(1,i)>=y0+ys+y01 && yy(1,i)<=360s(1,i)=0;v(1,i)=0;a(1,i)=0;subplot(2,3,5) ;plot(x,y,'r',xt,yt,eex,eey,'g',r0x,r0y,'k')%画图title('实际理论轮廓线')axis equal%使坐标轴比例相等grid on%画网格线。
哈工大机械原理大作业——齿轮——1号
Harbin Institute of Technology机械原理大作业3课程名称：机械原理设计题目：齿轮传动设计哈尔滨工业大学一、设计题目：如下图一个机械传动系统，运动运动由电动机1输入，经过机械传动系变速后由圆锥齿轮16输出三种不同转速。
选择一组传动系统的原始参数，据此设计该传动系统。
序号电机转速〔r/min〕输出轴转速〔r/min〕带传动最大传动比滑移齿轮传动定轴齿轮传动最大传动比模数圆柱齿轮圆锥齿轮一对齿轮最大传动比模数一对齿轮最大传动比模数7 1450 17 23 30 ≤2.8 ≤4.5 2 ≤4.5 3 ≤4 3二、传动比的分配计算：电动机转速n=1450r/min，输出转速n1=17r/min，n2=23 r/min，n3=30 r/min，带传动的最大传动比=2.8,滑移齿轮传动的最大传动比=4.5，圆柱齿轮传动的最大传动比=4.5，圆锥齿轮最大传动比=4。
根据传动系统的原始参数可知，传动系统的总传动比为：i1=1450/30=48.333i2=1450/23=63.043i3=1450/17=85.294传动系统的总传动比由带传动、滑移齿轮传动和定轴齿轮传动三局部实现。
设带传动的传动比为ipmax =2.8，滑移齿轮的传动比为iv1,iv2和iv3，令iv3=ivmax=4.5,那么定轴的传动比为if =85.294/(4.5*2.8)=6.769，从而iv1=48.333/〔6.769*2.8〕=2.550,iv2=3.326。
定轴齿轮每对的传动比为id==1.89。
三、滑移齿轮变速传动中每对齿轮的几何尺寸及重合度:经过计算、比拟，确定出三对滑移齿轮的齿数，其分别为：z5=17，z6=44，z 7=14,z8=47，z9=11，z10=50。
变位系数确实定：x5=x6=0； x7≥ha*(17-14)/17=0.176，取x7=0.18，x8=-0.18；x9≥ha*(17-11)/17=0.353，取x9=0.36；x10=-0.36。
哈工大机械原理大作业齿轮传动设计(word文档良心出品)
机械原理大作业（三）课程名称：设计题目：院系：姓名：学号：指导教师：设计时间：哈尔滨工业大学（威海）设计说明书1.设计题目 (2)2.传动比的分配计算 (3)3. 计算滑移齿轮变速传动中每对齿轮的基本几何尺寸 (3)（1）齿轮5、齿轮6 (4)（2）齿轮7、齿轮8 (4)（3）齿轮9、齿轮10 (5)4.计算定轴齿轮传动中每对齿轮的基本几何尺寸。
(5)（1）齿轮11、齿轮12 (5)（2）齿轮13、齿轮14 (6)（3）齿轮15、齿轮16 (6)5.每对齿轮的几何尺寸及重合度。
(6)7.实际设计参数 (14)1.设计题目如图所示一个机械传动系统，运动由电动机1输入，经过机械传动系统变速后由圆锥齿轮16输出三种不同的转速。
根据表中的传动系统原始参数设计该传动系统。
1.15，16.圆锥齿轮表机械传动系统原始参数2.传动比的分配计算电动机转速ni=1450r/min，输出转速n1=12r/min，n2=17r/min，n3=23r/min，带传动的最大传动比idmax=2.5,滑移齿轮传动的最大传动比ihmax=4，定轴齿轮传动的最大传动比ifmax=4。
根据传动系统的原始参数可知，传动系统的总传动比为i1=nin1=1450÷12=120.833i2=nin2=1450÷17=85.294i3=nin3=1450÷23=63.043传动系统的总传动比由带传动、滑移齿轮传动和定轴齿轮传动三部分实现。
设带传动的传动比为id，滑移齿轮的传动比为ih1、ih2和ih3，定轴齿轮传动的传动比为if，则总传动比i1=id*ih1*ifi2=id*ih2*ifi3=id*ih3*if令=ih1=ihmax=4则可得：定轴齿轮传动部分的传动比为if=i1/(id*ih1)=120.833/(2.5*4)=12.083滑移齿轮传动的传动比ih2=i2/(id*if)=85.294/(2.5*12.083)=2.824Ih3=i3/(id*if)=63.043/(2.5*12.083)=2.087定轴齿轮传动由3对齿轮传动组成，则每对齿轮的传动比为iv*iv*iv=if=12.083,iv=2.2953.计算滑移齿轮变速传动中每对齿轮的基本几何尺寸根据滑移齿轮变速传动系统中对齿轮齿数的要求，选择齿轮5、6、7、8、9和10为角度变位齿轮。
哈工大机械原理大作业连杆机构运动分析29
机械原理大作业一课程名称：机械原理设计题目：连杆机构运动分析-29院系：机电工程学院班级：设计者：学号：指导教师：设计时间：哈尔滨工业大学一、设计题目(29)二、VB编程得到动画、轨迹、位移、速度、加速度线图1.用VB编写程序得到的机构动画2.构件2上的E点轨迹4.构件7的速度线图三、计算说明书(一)建立机构运动分析数学模型对机构进行结构分析，找出基本杆组1.基本杆组划分a)杆件1为原动件b)构件2、3为RRR型II级杆组c)构件4、5为RRR型II级杆组d)构件6、7为RRP型II级杆组2.结构分析与自由度计算机构为平面连杆机构，活动构件数n=7,L P =10,转动副为9个，移动副为1个，高副为0个，则机构自由度为F=3n-2P L -P H =3x7-2x10=1。
各基本杆组的运动分析数学模型1.原动件（I 级构件）图1（1）位置分析B 点位置方程cos sin B A i iBA i i x x l y y l ϕϕ=+⎧⎨=+⎩ （1-1）（2）速度分析把（1-1）对时间t 求导得B 速度方程sin cos B B A i i i BB A i i i dx x x l dt dy y y l dt ϕϕϕϕ⎫==-⎪⎪⎬⎪==+⎪⎭g g g g g g （1-2）（3）加速度分析把（2-2）对时间他求导得B 加速度方程222222cos sin sin cos B B A i i i i i i BB A i i i i i i d x x x l l dt d y y y l l dt ϕϕϕϕϕϕϕϕ⎫==--⎪⎪⎬⎪==-+⎪⎭g g g g g g g g g g g g g g （1-3） 上式中 i i d dt ϕϕω==g22i i d dtϕϕα==g g 2. RRR 杆组的运动分析3.RRP杆组的运动分析图2（1） 位置分析内运动副C 的位置方程cos cos sin sin sin cos C B i i K j j jCB i i K j j j x x l x s l y y l y s l ϕϕϕϕϕϕ=+=+-⎧⎪⎨=+=++⎪⎩ （2-1）i l 的转角 0arcsin()ji j iA l l ϕϕ+=+ （2-2）式中0()sin ()cos B k j B K j A x x y y ϕϕ=---满足装配条件0j i A l l +≤ 滑块D 相对参考点K 的位移为sin cos cos sin C k j jC K j jjjx x l y y l s ϕϕϕϕ-+--==（2-3）滑块D 的位置方程 cos sin D K iD K ix x s y y s ϕϕ=+⎧⎨=+⎩ （2-4）（2） 速度分析杆i l 的角速度123sin cos j ji i Q Q Q ϕϕωϕ-+==g（2-5）滑块D 沿导路的移动速度 123cos sin i i i iQ l Q l s Q ϕϕ-+=g（2-6）式中1(sin cos )K B j j j j Q x x s l ϕϕϕ=--+ggg2(cos sin )K B j j j j Q y y s l ϕϕϕ=-+-g g g3sin sin cos cos i i j i i j Q l l ϕϕϕϕ=+内运动副C 的速度为sin cos C B i i iC B i i ix x l y y l ϕϕϕϕ⎧=-⎪⎨⎪=+⎩gg gg g g （2-7）外移动副D 的速度为cos sin sin cos D K i j jD K i j j x x s s y y s s ϕϕϕϕϕϕ⎧=+-⎪⎨⎪=++⎩gg g g g g g g（2-8） （3） 加速度分析 杆i l 的角加速度453sin s j ji i Q Q co Q ϕϕεϕ-+==g g（2-9）滑块D 沿导路移动的加速度为453cos sin i i i i Q l Q l s Q ϕϕ--=g g （2-10） 式中224cos (sin cos )(cos sin )2sin K B i i i j j j j j j j j j j Q x x l s l s l s ϕϕϕϕϕϕϕϕϕϕ=-+-+---g g g g g g g g g g225sin (cos sin )(sin cos )2cos K B i i i j j j j j j j j j j Q y y l s l s l s ϕϕϕϕϕϕϕϕϕϕ=-++--++g g g g g g g g g g内运动副C 的加速度22cos sin sin cos C B i i i j i j C B i i i j i jx x l l y y l l ϕϕϕϕϕϕϕϕ⎧=--⎪⎨⎪=-+⎩g gg g g g g g gg g g g g（2-11） 滑块D 的加速度22cos sin cos 2sin sin cos sin 2cos D K j j j j j j jD K j j j j j j jx x s s s s y y s s s s ϕϕϕϕϕϕϕϕϕϕϕϕϕϕ⎧=+---⎪⎨⎪=++-+⎩g g g g g g g g g g gg g g g g g g gg g g（2-12） 以上为基本杆组的通用分析，根据已知条件进行赋值均在程序中体现(二) 计算流程框图（见下页）(三)VB编写程序清单'定义变量Dim xA As Double 'A点的坐标、速度、加速度Dim yA As DoubleDim vxA As DoubleDim vyA As DoubleDim axA As DoubleDim ayA As DoubleDim xB As Double 'B点Dim yB As DoubleDim vxB As DoubleDim vyB As DoubleDim axB As DoubleDim ayB As DoubleDim xD As Double 'D点Dim yD As DoubleDim vxD As DoubleDim vyD As DoubleDim axD As DoubleDim ayD As DoubleDim xE As Double 'E点Dim yE As DoubleDim vxE As DoubleDim vyE As DoubleDim axE As DoubleDim ayE As DoubleDim xG As Double 'G点Dim yG As DoubleDim vxG As DoubleDim vyG As DoubleDim axG As DoubleDim ayG As DoubleDim xH As Double 'H点Dim yH As DoubleDim vxH As DoubleDim vyH As DoubleDim axH As DoubleDim ayH As DoubleDim xM As Double 'M点Dim yM As DoubleDim vxM As DoubleDim vyM As DoubleDim axM As DoubleDim xN As Double '参考点NDim yN As DoubleDim vxN As DoubleDim vyN As DoubleDim axN As DoubleDim ayN As DoubleDim xK As Double 'K点Dim yK As DoubleDim vxK As DoubleDim vyK As DoubleDim axK As DoubleDim ayK As DoubleDim delt1 As Double 'AB杆初始转角Dim L1 As Double 'AB杆长Dim f1 As Double 'AB杆转角Dim w1 As Double 'AB杆角速度Dim e1 As Double 'AB杆角加速度Dim L2 As Double 'BC杆Dim f2 As DoubleDim w2 As DoubleDim e2 As DoubleDim L3 As Double 'CD杆Dim f3 As DoubleDim w3 As DoubleDim e3 As DoubleDim LBE As Double 'BE杆Dim delt2 As DoubleDim L4 As Double 'EF杆Dim f4 As DoubleDim w4 As DoubleDim e4 As DoubleDim LFG As Double 'FG杆Dim f5 As DoubleDim w5 As DoubleDim e5 As DoubleDim deltGH As Double 'GH杆Dim LGH As DoubleDim L6 As Double 'HM（HK）杆Dim f6 As DoubleDim w6 As DoubleDim e6 As DoubleDim Lj As Double 'MK杆Dim wj As DoubleDim ej As DoubleDim ss As Double '移动副K的位移Dim vss As Double '移动副K的速度Dim ass As Double '移动副K的加速度Dim pi As DoubleDim pa As DoubleDim i As DoubleDim fj1 As DoublePrivate Sub Command1_Click() '点E的轨迹Picture1.Scale (-50, 250)-(600, -100)Picture1.Line (-50, 0)-(600, 0) 'XPicture1.Line (0, 250)-(0, -100) 'YFor i = -50 To 600 Step 50 'X轴坐标Picture1.DrawStyle = 2Picture1.Line (i, 250)-(i, -100)Picture1.CurrentX = i - 0: Picture1.CurrentY = 0 Picture1.Print iNext iFor i = -100 To 250 Step 50 'Y轴坐标Picture1.DrawStyle = 2Picture1.Line (-50, i)-(600, i)Picture1.CurrentX = -20: Picture1.CurrentY = i + 7 Picture1.Print iNext iFor fj1 = 0 To 360 Step 0.01f1 = fj1 * paCall RR1Call RRR1Call RR2Picture1.PSet (xE, yE)Next fj1End SubPrivate Sub Command2_Click() '求点G的位移Picture2.Scale (-20, 300)-(380, 200)Picture2.Line (-20, 0)-(380, 0) 'XPicture2.Line (0, 300)-(0, 200) 'YFor i = 0 To 360 Step 30 'X轴坐标Picture2.DrawStyle = 2Picture2.Line (i, 300)-(i, 0)Picture2.CurrentX = i - 10: Picture2.CurrentY = 0Picture2.Print iNext iFor i = 200 To 300 Step 10 'Y轴坐标Picture2.Line (0, i)-(380, i)Picture2.CurrentX = -25: Picture2.CurrentY = i Picture2.Print iNext iFor fj1 = 0 To 360 Step 0.01f1 = fj1 * paCall RR1Call RRR1Call RR2Call RRR2Call RR3Call RRP1Picture2.PSet (fj1, ss)Next fj1End SubPrivate Sub Command3_Click() '求点G的速度Picture3.Scale (-30, 600)-(380, -600)Picture3.Line (-30, 0)-(380, 0) 'XPicture3.Line (0, 600)-(0, -600) 'YFor i = 0 To 360 Step 30 'X轴坐标Picture3.DrawStyle = 2Picture3.Line (i, 600)-(i, -600)Picture3.CurrentX = i - 10: Picture3.CurrentY = 0 Picture3.Print iNext iFor i = -600 To 600 Step 50 'Y轴坐标Picture3.Line (0, i)-(380, i)Picture3.CurrentX = -28: Picture3.CurrentY = i Picture3.Print iNext iFor fj1 = 0 To 360 Step 0.01f1 = fj1 * paCall RR1Call RRR1Call RR2Call RRR2Call RR3Call RRP1Picture3.PSet (fj1, vss)Next fj1End SubPrivate Sub Command4_Click() '求点G的加速度Picture4.Scale (-20, 20000)-(380, -10000)Picture4.Line (-20, 0)-(380, 0) 'XPicture4.Line (0, 20000)-(0, -10000) 'YFor i = 0 To 360 Step 30 'X轴坐标Picture4.DrawStyle = 2Picture4.Line (i, 20000)-(i, -10000)Picture4.CurrentX = i - 10: Picture4.CurrentY = 0 Picture4.Print iNext iFor i = -10000 To 20000 Step 1000 'Y轴坐标Picture4.Line (0, i)-(380, i)Picture4.CurrentX = -25: Picture4.CurrentY = i + 5 Picture4.Print iNext iFor fj1 = 0 To 360 Step 0.01f1 = fj1 * paCall RR1Call RRR1Call RR2Call RRR2Call RR3Call RRP1Picture4.PSet (fj1, ass)Next fj1End SubPrivate Sub Form_Load() '赋初值xA = 0yA = 0vxA = 0vyA = 0axA = 0ayA = 0xD = 730yD = -138vxD = 0vyD = 0axD = 0ayD = 0xG = 465yG = -248vxG = 0vyG = 0axG = 0ayG = 0xN = 545yN = 0vxN = 0vyN = 0axN = 0ayN = 0delt1 = 0L1 = 200w1 = 10e1 = 0L2 = 670L3 = 350delt2 = 0LBE = 335L4 = 380LFG = 130deltGH = 56 * paLGH = 100L6 = 486Lj = 0fj = pi / 2wj = 0ej = 0pi = 3.14pa = pi / 180End SubPublic Sub RR1() '杆ABxB = xA + L1 * Cos(f1 + delt1)yB = yA + L1 * Sin(f1 + delt1)vxB = vxA - w1 * L1 * Sin(f1 + delt1)vyB = vyA + w1 * L1 * Cos(f1 + delt1)axB = axA - w1 ^ 2 * L1 * Cos(f1 + delt1) - e1 * L1 * Sin(f1 + delt1) ayB = ayA - w1 ^ 2 * L1 * Sin(f1 + delt1) + e1 * L1 * Cos(f1 + delt1) End SubPublic Sub RRR1() '杆BCDDim xC As DoubleDim yC As DoubleDim vxC As DoubleDim vyC As DoubleDim axC As DoubleDim ayC As DoubleDim LBD As DoubleDim JCBD As DoubleDim val1 As DoubleDim fDB As DoubleDim C2 As DoubleDim C3 As DoubleDim S2 As DoubleDim S3 As DoubleDim G11 As DoubleDim G12 As DoubleDim G13 As DoubleLBD = Sqr((xD - xB) ^ 2 + (yD - yB) ^ 2)If LBD > L2 + L3 And LBD < Abs(L2 - L3) ThenIf MsgBox("RRR杆组杆长不符合要求", vbOKOnly, "提示") = 1 Then EndElseEnd IfElseEnd IfIf LBD < L2 + L3 And LBD > Abs(L2 - L3) Thenval1 = (L2 ^ 2 + LBD ^ 2 - L3 ^ 2) / (2 * L2 * LBD)JCBD = Atn(-val1 / Sqr(-val1 * val1 + 1)) + 2 * Atn(1)ElseEnd IfIf LBD = L2 + L3 ThenJCBD = 0ElseEnd IfIf LBD = Abs(L2 - L3) ThenIf L2 > L3 ThenJCBD = 0ElseEnd IfIf L2 < L3 ThenJCBD = piElseEnd IfElseEnd IfIf xD > xB And yD >= yB Then '第一象限fDB = Atn((yD - yB) / (xD - xB))ElseEnd IfIf xD = xB And yD > yB ThenfDB = pi * 2ElseEnd IfIf xD < xB And yD >= yB Then '第二象限fDB = pi + Atn((yD - yB) / (xD - xB)) ElseEnd IfIf xD < xB And yD < yB Then '第三象限fDB = pi + Atn((yD - yB) / (xD - xB)) ElseEnd IfIf xB = xD And yD < yB ThenfDB = 3 * pi / 2ElseEnd IfIf xD > xB And yD <= yB Then '第四象限fDB = 2 * pi + Atn((yD - yB) / (xD - xB)) ElseEnd Iff2 = fDB + JCBDxC = xB + L2 * Cos(f2)yC = yB + L2 * Sin(f2)If xC > xD And yC >= yD Then '第一象限f3 = Atn((yC - yD) / (xC - xD))ElseEnd IfIf xC = xD And yC > yD Thenf3 = pi / 2ElseEnd IfIf xC < xD And xC >= xD Then '第二象限f3 = pi + Atn((yC - yD) / (xC - xD))ElseEnd IfIf xC < xD And xC < xD Then '第三象限f3 = pi + Atn((yC - yD) / (xC - xD))ElseEnd IfIf xC = xD And xC < xD Thenf2 = 3 * pi / 2ElseEnd IfIf xC > xD And xC <= xD Then '第四象限f2 = 2 * pi + Atn((yC - yD) / (xC - xD))ElseEnd IfC2 = L2 * Cos(f2)S2 = L2 * Sin(f2)C3 = L3 * Cos(f3)S3 = L3 * Sin(f3)G11 = C2 * S3 - C3 * S2w2 = (C3 * (vxD - vxB) + S3 * (vyD - vyB)) / G11w3 = (C2 * (vxD - vxB) + S2 * (vyD - vyB)) / G11vxC = vxB - w2 * L2 * Sin(f2)vyC = vyB + w2 * L2 * Cos(f2)G12 = axD - axB + w2 ^ 2 * C2 - w3 ^ 2 * C3G13 = ayD - ayB + w2 ^ 2 * S2 - w3 ^ 2 * S3e2 = (G12 * C3 + G13 * S3) / G11e3 = (G12 * C2 + G13 * S2) / G11axC = axB - e2 * L2 * Sin(f2) - w2 ^ 2 * L2 * Cos(f2)ayC = ayB - e2 * L2 * Cos(f2) - w2 ^ 2 * L2 * Sin(f2)End SubPublic Sub RR2() '杆BExE = xB + LBE * Cos(f2 + delt2)yE = yB + LBE * Sin(f2 + delt2)vxE = vxB - w2 * LBE * Sin(f2 + delt2)vyE = vyB + w2 * LBE * Cos(f2 + delt2)axE = axB - w2 ^ 2 * LBE * Cos(f2 + delt2) - e2 * LBE * Sin(f2 + delt2) ayE = ayB - w2 ^ 2 * LBE * Sin(f2 + delt2) + e2 * LBE * Cos(f2 + delt2) End SubPublic Sub RRR2() '杆EFGDim xF As DoubleDim yF As DoubleDim vxF As DoubleDim vyF As DoubleDim axF As DoubleDim ayF As DoubleDim LEG As DoubleDim JFEG As DoubleDim val2 As DoubleDim fGE As DoubleDim C4 As DoubleDim C5 As DoubleDim S4 As DoubleDim S5 As DoubleDim G21 As DoubleDim G22 As DoubleDim G23 As DoubleLEG = Sqr((xG - xE) ^ 2 + (yG - yE) ^ 2)If LEG > L4 + LFG And LEG < Abs(L4 - LFG) ThenIf MsgBox("RRR杆组杆长不符合要求", vbOKOnly, "提示") = 1 Then EndElseEnd IfElseEnd IfIf LEG < L4 + LFG And LEG > Abs(L4 - LFG) Thenval2 = (L4 ^ 2 + LEG ^ 2 - LFG ^ 2) / (2 * L4 * LEG)JFEG = Atn(-val2 / Sqr(-val2 * val2 + 1)) + 2 * Atn(1)ElseEnd IfIf LEG = L4 + LFG ThenJFEG = 0ElseEnd IfIf LEG = Abs(L4 - LFG) ThenIf L4 > LFG ThenJFEG = 0ElseEnd IfIf L4 < LFG ThenJFEG = piElseEnd IfElseEnd IfIf xG > xE And yG >= yE Then '第一象限fGE = Atn((yG - yE) / (xG - xE))ElseEnd IfIf xG = xE And yG > yE ThenfGE = pi * 2ElseEnd IfIf xG < xE And yG >= yE Then '第二象限fGE = pi + Atn((yG - yE) / (xG - xE)) ElseEnd IfIf xG < xE And yG < yE Then '第三象限fGE = pi + Atn((yG - yE) / (xG - xE)) ElseEnd IfIf xE = xG And yG < yE ThenfGE = 3 * pi / 2ElseEnd IfIf xG > xE And yG <= yE Then '第四象限fGE = 2 * pi + Atn((yG - yE) / (xG - xE)) ElseEnd Iff4 = fGE - JFEGxF = xE + L4 * Cos(f4)yF = yE + L4 * Sin(f4)If xF > xG And yF >= yG Then '第一象限f5 = Atn((yF - yG) / (xF - xG))ElseEnd IfIf xF = xG And yF > yG Thenf5 = pi / 2ElseEnd IfIf xF < xG And xF >= xG Then '第二象限f5 = pi + Atn((yF - yG) / (xF - xG))ElseEnd IfIf xF < xG And xF < xG Then '第三象限f5 = pi + Atn((yF - yG) / (xF - xG))ElseEnd IfIf xF = xG And xF < xG Thenf4 = 3 * pi / 2ElseEnd IfIf xF > xG And xF <= xG Then '第四象限f4 = 2 * pi + Atn((yF - yG) / (xF - xG)) ElseEnd IfC4 = L4 * Cos(f4)S4 = L4 * Sin(f4)C5 = LFG * Cos(f5)S5 = LFG * Sin(f5)G21 = C4 * S5 - C5 * S4w4 = (C5 * (vxG - vxE) + S5 * (vyG - vyE)) / G21w5 = (C4 * (vxG - vxE) + S4 * (vyG - vyE)) / G21vxF = vxE - w4 * L4 * Sin(f4)vyF = vyE + w4 * L4 * Cos(f4)G22 = axG - axE + w4 ^ 2 * C4 - w5 ^ 2 * C5G23 = ayG - ayE + w4 ^ 2 * S4 - w5 ^ 2 * S5e4 = (G22 * C5 + G23 * S5) / G21e5 = (G22 * C4 + G23 * S4) / G21axF = axE - e4 * L4 * Sin(f4) - w4 ^ 2 * L4 * Cos(f4)ayF = ayE - e4 * L4 * Cos(f4) - w4 ^ 2 * L4 * Sin(f4)End SubPublic Sub RR3() '杆GHxH = xG + LGH * Cos(f5 + deltGH)yH = yG + LGH * Sin(f5 + deltGH)vxH = vxG - w5 * LGH * Sin(f5 + deltGH)vyH = vyG + w5 * LGH * Cos(f5 + deltGH)axH = axG - w5 ^ 2 * LGH * Cos(f5 + deltGH) - e5 * LGH * Sin(f5 + deltGH) ayH = ayG - w5 ^ 2 * LGH * Sin(f5 + deltGH) + e5 * LGH * Cos(f5 + deltGH) End SubPublic Sub RRP1() '杆组H(M)KDim A0 As DoubleDim Q1 As DoubleDim Q2 As DoubleDim Q3 As DoubleDim Q4 As DoubleDim Q5 As DoubleDim val3 As DoubleA0 = Lj + (yN - yH) * Cos(fj) - (xN - xH) * Sin(fj)val3 = A0 / L6fi = Atn(val3 / Sqr(-val3 * val3 + 1)) + fjxM = xH + L6 * Cos(fi)yM = yH + L6 * Sin(fi)ss = (xM - xN) * Cos(fj) + (yM - yN) * Sin(fj)xK = xN + ss * Cos(fj)yK = yN + ss * Sin(fj)Q1 = vxN - vxH - wj * (ss * Sin(fj) + Lj * Cos(fj))Q2 = vyN - vyH + wj * (ss * Cos(fj) - Lj * Sin(fj))Q3 = L6 * Sin(fi) * Sin(fj) + L6 * Cos(fi) * Cos(fj)w6 = (-Q1 * Sin(fj) + Q2 * Cos(fj)) / Q3vss = -(Q1 * L6 * Cos(fi) + Q2 * L6 * Sin(fi)) / Q3vxM = vxH - w6 * L6 * Sin(fi)vyM = vyH + w6 * L6 * Cos(fi)vxK = vxN + vss * Cos(fj) - ss * wj * Sin(fj)vyK = vyN + vss * Sin(fj) + ss * wj * Cos(fj)Q4 = axN - axH + w6 ^ 2 * L6 * Cos(fi) - ej * (ss * Sin(fj) + Lj * Cos(fj)) - wj ^ 2 * (ss * Cos(fj) - Lj * Sin(fj)) - 2 * vss * wj * Sin(fj)Q5 = ayN - ayH + w6 ^ 2 * L6 * Sin(fi) + ej * (ss * Cos(fj) - Lj * Sin(fj)) - wj ^ 2 * (ss * Sin(fj) + Lj * Cos(fj)) + 2 * vss * wj * Cos(fj)e6 = (-Q4 * Sin(fj) + Q5 * Cos(fj)) / Q3ass = (-Q4 * L6 * Cos(fi) - Q5 * L6 * Sin(fi)) / Q3axM = axH - e6 * L6 * Sin(fi) - w6 ^ 2 * L6 * Cos(fi)ayM = ayH + e6 * L6 * Cos(fi) - w6 ^ 2 * L6 * Sin(fi)axK = axN + ass * Cos(fj) - ss * ej * Sin(fj) - ss * wj ^ 2 * Cos(fj) - 2 * vss * wj * Sin(fj)ayK = ayN + ass * Sin(fj) + ss * ej * Cos(fj) - ss * wj ^ 2 * Sin(fj) + 2 * vss * wj * Cos(fj)End Sub(四)计算结果分析用VB编程可得构件6的位移、速度、加速度图像，由图像可得机构运动一个周期内每个时刻的运动参数值，较为直观。
机械原理大作业一连杆机构参考例子
机械原理大作业课程名称：机械原理设计题目：机械原理大作业院系：汽车工程学院车辆工程班级：1101201姓名：。
学号：。
指导教师：游斌弟大作业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的角位移、角速度和角加速度，并对计算结果进行分析。
图 12、对机构进行结构分析该机构由I 级杆组RR （原动件1）、II 级杆组RRR （杆2、杆3）和II 级杆组RPR （滑块4及杆5）组成。
I 级杆组RR ，如图2所示；II 级杆组RRR ，如图3所示；II 级杆组RPR ，如图4所示。
图2 图 2图4 3、建立坐标系建立以点A为原点的固定平面直角坐标系4、各基本杆组运动分析的数学模型（1）同一构件上点的运动分析：如图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 ϕϕ;图 3==-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 Bd y a l y dtωϕω。
图 4（2）RRRII 级杆组的运动分析如图6所示是由三个回转副和两个构件组成的II 级组。
已知两杆的杆长2l 、3l 和两个外运动副B 、D 的位置（B x 、B y 、D x 、D y ）、速度（ xB yB xD yD v v v v 、、、 ） 图6和加速度（xB yB xD yD a a a a 、、、）。
机械原理大作业-连杆机构
设计内容
结论
1. 机构结构分析 1）计算机构自由度，确定机构是否有确定运动。
机构自度 F=1
本机构中，n= 5 ， pL= 7
，pH = 0
则有：F=3n-2PL –PH= 3*5-2*7=1
是否有确定
机构确定运动判断： 因 F= 1 ，原动件个数= 1 ，可知： 自由度等于原动件数，机构有确定 的运动。
➢ 求解速度 vB 、角速度 2 ： 矢量方程：
大小 ？ √ ？ 方向 ⊥BD ⊥OA ⊥AB
逆时针方向 ➢ 求解速度 vC 、2 杆质心 S2 的速度 vs2 ：
vB =1.57m/s 2 =0.36rad/s， 逆时针方向
vC =1.56m/s vs2 =1.53m/s
-3-
➢ 求解速度 vE 、角速度 4 ： 矢量方程：
动态静力学参数：m2=20 kg ， JS2 = 1.1 kg m2 ，m5=50 kg
表 2 阻抗力参数表
班级序号 1
2
3
4
5
6
7
8
9
10
Fr /N
500 600 700
800
900
1000 1100 1200 1300 1400
-1-
机械原理模块训练一
二、 训练要求
通过对干草压缩机六杆机构进行结构分析、运动分析和力分析，对该设备的运动性能做出 定量的计算，为新设备的设计与评价提供依据。 机构结构分析：了解机构组成，学会对机构工作原理简图表达，判断机构的结构组成是否可行。 机构运动分析：已知原动件运动参数情况下，学会求解机构输出端和机构中关键点的运动参数。 机构动态静力分析：当已知工作阻力时，需要给机器配置动力，通过平衡力求解可以获得；当 机器安装时需要知道支座的反力，可通过动态静力分析获得。
哈工大机械原理大作业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所示。
哈工大机械原理大作业齿轮传动设计
Harbin Institute of Technology机械原理说明书课程名称：机械原理设计题目：齿轮传动设计院系：机电学院班级：设计者：学号：指导教师：**设计时间：2013/6/16哈尔滨工业大学一、设计题目如图所示一个机械传动系统，运动由电动机1输入，经过机械传动系统变速后由圆锥齿轮16输出三种不同的转速。
二、机构运动简图1.电动机 2，4.皮带轮 3.皮带5,6,7,8,9,10,11,12,13,14.圆柱齿轮15,16.圆锥齿轮 序号 电机转速 （r/min ） 输出轴转速 （r/min ）带传动最大传动比 滑移齿轮传动 定轴齿轮传动最大传动比 模数 圆柱齿轮圆锥齿轮 一对齿轮最大传动比模数 一对齿轮最大传动比模数 34 1450 50 57 65 ≤2.5≤42≤43≤43四、传动比分配计算电动机转速n i =1450r/min,输出转速为n 1o =50 r/min;n 2o =57 r/min 和n 3o =65 r/min ，带传动的最大传动比为i m ax p =2.5,滑移齿轮传动的最大传动比为i m ax v =4,定轴齿轮传动的最大传动比为i m ax d =4。
根据传动系统的原始参数可知，传动系统的总传动比为i 1=1o i n n =501450=29i 2=2o i n n =571450=25.438i 3=3o i n n =651450=22.308传动系统的总传动比由带传动、滑移齿轮传动和定轴齿轮传动三部分实现。
设带传动的传动比为i m ax p ，滑移齿轮的传动比为i 1v ,i 2v 和i 3v ,定轴齿轮传动的传动比为i f ,则传动比 i 1= i m ax p i 1v i f i 2= i m ax p i 2v i fi 3= i m ax p i 3v i f令i 1v = i m ax v =4则可得定轴齿轮传动部分的传动比为i f =max max 1v p i i i =4*5.229=2.9则可得滑移齿轮的传动比i 2v =fp i i i max 2=9.2*5.2438.25=3.51i 3v =fp i i i max 3=9.2*5.2308.22=3.08定轴齿轮传动由3对齿轮传动组成，则每对齿轮的传动比为i d =3f i =39.2=1.324五、滑移齿轮传动齿数的确定（m=2）56z z =2268=3.091 78z z =2070=3.5 910z z =1872=4 齿轮9/齿轮10：角度变位正传动。