机械原理课程设计程序.docx

附录：机械原理课设相关程序Part1.符号说明机构名称符号符号意义 曲 柄 滑 块 机 构s滑块位移v滑块速度a滑块加速度x曲柄转角x1连杆转角w1连杆角速度偏置直动滚子从动件盘形凸轮机构s从动件位移v从动件类速度a从动件类加速度x凸轮转角n0,m0凸轮基圆横纵坐标n,m凸轮理论廓线横纵坐标n_,m_凸轮实际廓线横纵坐标Part2.程序1. 曲柄滑块：（1） s-x图：clear ;x= linspace(0,360,2000);x0=180-atan(15/26.104)*180/pi;s=28.647*cos(x+x0)/180*pi)+58.754*cos(asin(-0.448)*sin(x+x0)/180*pi)-0.255)-26.7623;plot(x,s);anis(0 360 0 65)(2) v-x图clear ;x= linspace(0,360,2000);x0=180-atan(15/26.104)*180/pi;x1=asin(-0.488)*sin(x+x0)/180*pi)-0.255)+2*pi*(asin(-0.488)*sin(x+x0)/180*pi)-0.255)0);w1=(-173.995)*cos(x+x0)/180*pi-x1)/58.754;v=(-173.995)*sin(x+x0)/180*pi)-58.754*(w1).*sin(x1);plot(x,v)（3） a-x图clear ;x= linspace(0,360,2000);x0=180-atan(15/26.104)*180/pi;x1=asin(-0.488)*sin(x+x0)/180*pi)-0.255)+2*pi*(asin(-0.488)*sin(x+x0)/180*pi)-0.255)0);w1=(-173.995)*cos(x+x0)/180*pi)/28.647./cos(x1);alfa=1056.799*sin(x+x0)/180*pi)+58.754*(w1).2.*sin(x1)/58.754./cos(x1);a=(-1056.799)*cos(x+x0)/180*pi)-58.754*(alfa).*sin(x1)-58.754*(w1).2.*cos(x1);plot(x,a)2.凸轮机构（1）s-x图clear ;x1=0:0.1:120;x2=120:0.1:160;x3=160:0.1:280;x4=280:0.1:360;s1=18.5-18.5*cos(1.5*(x1)/180*pi);s2=37;s3=86.3439-17.67*(x3)*pi/180+5.89*sin(3*(x3)/180*pi-2/3*pi);s4=0;plot(x1,s1);hold on;plot(x2,s2);hold on;plot(x3,s3);hold on;plot(x4,s4)（2） v-x图clear ;x1=0:0.1:120;x2=120:0.1:160;x3=160:0.1:280;x4=280:0.1:360;v1=27.75*sin(1.5*(x1)/180*pi);v2=0;v3=17.67*cos(3*(x3)/180*pi-2/3*pi)-17.67;v4=0;plot(x1,v1);hold on;plot(x2,v2);hold on;plot(x3,v3);hold on;plot(x4,v4)(3) a-x图clear ;x1=0:0.1:120;x2=120:0.1:160;x3=160:0.1:280;x4=280:0.1:360;a1=41.625*cos(1.5*(x1)/180*pi);a2=0;a3=(-53)*sin(3*(x3)/180*pi-2/3*pi);a4=0;plot(x1,a1);hold on;plot(x2,a2);hold on;plot(x3,a3);hold on;plot(x4,a4)(4) 凸轮基圆半径和最大压力角的确定#include#include#include#define PI 3.1415926void ydgl();float d,am,r0,e,h,p,dt0,d0,s,ds,s0,al;int f; int main () /*输入凸轮基础尺寸和初始基圆半径*/printf(Determination the prime circle of camn); printf(input r0=);scanf(%f,&r0);printf(input h=);scanf(%f,&h);printf(input e=);scanf(%f,&e);printf(input dt0=);scanf(%f,&dt0);am=0;p=PI/180;printf(The initial data: r0=%fn,r0);printf(h=%f e=%f dt0=%fn,h,e,dt0);/*检验初始基圆半径，并求出符合要求的基圆半径*/dor0=r0+5;for(f=0;f30*p);printf(The intermediate results: am=%fn,am/p);printf(r0=%fn,r0);/*缩小基圆半径，求出最佳基圆半径和最大压力角*/dor0=r0-1;if(r0=50*/break;for(f=0;fdt0;f=f+1)d=PI*f/dt0;ydgl();while(amam)am=al;（求出r0=50.0mm,最大压力角11.537度）(5) 凸轮廓线（理论和实际）clear ;x=0:0.01:360;%画基圆n0=50*cos(x/180*pi);m0=50*sin(x/180*pi);plot(n0,m0,-b);hold on;x1=0:0.1:120;x2=120:0.1:160;x3=160:0.1:280;x4=280:0.1:360;s0=sqrt(50*50-10*10);s1=18.5-18.5*cos(1.5*(x1)/180*pi);%画理论廓线s2=37;s3=86.3439-17.67*(x3)*pi/180+5.89*sin(3*(x3)/180*pi-2/3*pi);s4=0;n1=(s0+s1).*sin(x1)/180*pi)+10*cos(x1)/180*pi);n2=(s0+s2).*sin(x2)/180*pi)+10*cos(x2)/180*pi);n3=(s0+s3).*sin(x3)/180*pi)+10*cos(x3)/180*pi);n4=(s0+s4).*sin(x4)/180*pi)+10*cos(x4)/180*pi);m1=(s0+s1).*cos(x1)/180*pi)-10*sin(x1)/180*pi); m2=(s0+s2).*cos(x2)/180*pi)-10*sin(x2)/180*pi); m3=(s0+s3).*cos(x3)/180*pi)-10*sin(x3)/180*pi); m4=(s0+s4).*cos(x4)/180*pi)-10*sin(x4)/180*pi);plot(n1,m1,-r,n2,m2,-r,n3,m3,-r,n4,m4,-r);hold on;v1=27.75*sin(1.5*(x1)/180*pi);%画实际廓线v2=0;v3=17.67*cos(3*(x3)/180*pi-2/3*pi)-17.67;v4=0;dn1=(v1-10).*sin(x1)/180*pi)+(s0+s1).*cos(x1)/180*pi); dn2=(v2-10).*sin(x2)/180*pi)+(s0+s2).*cos(x2)/180*pi);dn3=(v3-10).*sin(x3)/180*pi)+(s0+s3).*cos(x3)/180*pi);dn4=(v4-10).*sin(x4)/180*pi)+(s0+s4).*cos(x4)/180*pi);dm1=(v1-10).*cos(x1)/180*pi)-(s0+s1).*sin(x1)/180*pi);dm2=(v2-10).*cos(x2)/180*pi)-(s0+s2).*sin(x2)/180*pi);dm3=(v3-10).*cos(x3)/180*pi)-(s0+s3).*sin(x3)/180*pi);dm4=(v4-10).*cos(x4)/180*pi)-(s0+s4).*sin(x4)/180*pi);y1=atan(dn1)./(-dm1)+pi*(atan(dn1)./(-dm1).*(dn1)0);y2=atan(dn2)./(-dm2)+pi*(atan(dn2)./(-dm2).*(dn2)0);y3=atan(dn3)./(-dm3)+pi*(atan(dn3)./(-dm3).*(dn3)0);y4=atan(dn4)./(-dm4)+pi*(atan(dn4)./(-dm4).*(dn4)0);n1_=n1-5*cos(y1);n2_=n2-5*cos(y2);n3_=n3-5*c