机械原理课程设计--内燃机机构设计及其运动分析.doc

机械原理课程设计说 明 书设计题目 内燃机机构设计及其运动分析 xx 学院xx班学号xx 设计者: xxx指导教师: xxx2010 年06 月 16 日 目录一、设计过程 （1）齿轮机构传动设计（2）连杆机构设计和运动分析（3）凸轮机构设计 二、连杆机构与凸轮机构运动相关图像 （1）连杆机构图像 （2）凸轮机构从动件运动图像 （3）凸轮理论廓线与实际廓线图像 三、相关电算数据（主要利用matlab求得数据） （1）连杆机构程序运算结果（2）凸轮机构轮廓线电算结果 四、相关程序及源代码（1）连杆机构的matlab源代码（2）凸轮轮廓曲线matlab源代码（3）凸轮轮廓曲线c+源代码（4）凸轮机构从动件源代码五、主要参考文献 二、连杆机构与凸轮机构运动相关图像（1）连杆机构图像 滑块位移角度图像滑块速度角度图像滑块加速度角度图像（2）凸轮机构从动件运动图像从动件位移图像从动件类速度图像（3）凸轮理论廓线与实际廓线图像凸轮理论廓线凸轮实际廓线内位理论廓线，外为实际廓线三、相关电算数据（主要利用matlab求得数据）（1）连杆机构程序运算结果角度i/deg位移y mm/s速度v mm/s加速度a mm/s*si=0.0y=-78.899v=0.0a=214286.675i=5.0y=-78.553v=-3170.8a=213168.982i=10.0y=-77.519v=-6243.1a=214595.729i=15.0y=-75.804v=-9120.5a=218646.334i=20.0y=-73.421v=-11710.9a=225285.858i=25.0y=-70.386v=-13928.2a=234352.851i=30.0y=-66.724v=-15694.7a=245548.646i=35.0y=-62.459v=-16943.4a=258428.172i=40.0y=-57.623v=-17619.8a=272392.523i=45.0y=-52.252v=-17684.3a=286684.115i=50.0y=-46.385v=-17114.2a=300386.278i=55.0y=-40.065v=-15906.0a=312430.567i=60.0y=-33.337v=-14076.5a=321616.659i=65.0y=-26.250v=-11664.4a=326650.645i=70.0y=-18.858v=-8730.8a=326206.964i=75.0y=-11.213v=-5358.0a=319015.974i=80.0y=-3.373v=-1648.3a=303972.999i=85.0y=4.604v=2279.9a=280256.639i=90.0y=12.658v=6296.3a=247436.938i=95.0y=20.728v=10264.8a=205551.382i=100.0y=28.752v=14049.5a=155131.037i=105.0y=36.668v=17521.5a=97170.083i=110.0y=44.416v=20563.9a=33045.527i=115.0y=51.934v=23077.1a=-35595.203i=120.0y=59.163v=24981.9a=-106956.691i=125.0y=66.047v=26221.2a=-179208.582i=130.0y=72.530v=26760.7a=-250577.000i=135.0y=78.562v=26588.6a=-319410.128i=140.0y=84.095v=25714.1a=-384219.943i=145.0y=89.084v=24166.2a=-443705.298i=150.0y=93.491v=21991.0a=-496762.170i=155.0y=97.280v=19250.0a=-542485.877i=160.0y=100.422v=16017.8a=-580168.569i=165.0y=102.892v=12379.7a=-609293.799i=170.0y=104.670v=8429.7a=-629528.985i=175.0y=105.743v=4268.3a=-640716.014i=180.0y=106.101v=-0.0a=-642860.024i=185.0y=105.743v=-4268.3a=-636116.505i=190.0y=104.670v=-8429.7a=-620777.067i=195.0y=102.892v=-12379.7a=-597254.513i=200.0y=100.422v=-16017.8a=-566068.143i=205.0y=97.280v=-19250.0a=-527830.496i=210.0y=93.491v=-21991.0a=-483236.831i=215.0y=89.084v=-24166.2a=-433058.608i=220.0y=84.095v=-25714.1a=-378141.816i=225.0y=78.562v=-26588.6a=-319410.128i=230.0y=72.530v=-26760.7a=-257871.445i=235.0y=66.047v=-26221.2a=-194624.448i=240.0y=59.163v=-24981.9a=-130859.640i=245.0y=51.934v=-23077.1a=-67847.792i=250.0y=44.416v=-20563.9a=-6908.675i=255.0y=36.668v=-17521.5a=50644.315i=260.0y=28.752v=-14049.5a=103582.212i=265.0y=20.728v=-10264.8a=150851.855i=270.0y=12.658v=-6296.3a=191663.829i=275.0y=4.604v=-2279.9a=225557.113i=280.0y=-3.373v=1648.3a=252424.174i=285.0y=-11.213v=5358.0a=272490.204i=290.0y=-18.858v=8730.8a=286252.761i=295.0y=-26.250v=11664.4a=294398.055i=300.0y=-33.337v=14076.5a=297713.708i=305.0y=-40.065v=15906.0a=297014.700i=310.0y=-46.385v=17114.2a=293091.832i=315.0y=-52.252v=17684.3a=286684.115i=320.0y=-57.623v=17619.8a=278470.649i=325.0y=-62.459v=16943.4a=269074.862i=330.0y=-66.724v=15694.7a=259073.984i=335.0y=-70.386v=13928.2a=249008.231i=340.0y=-73.421v=11710.9a=239386.284i=345.0y=-75.804v=9120.5a=230685.620i=350.0y=-77.519v=6243.1a=223347.647i=355.0y=-78.553v=3170.8a=217768.492i=360.0y=-78.899v=-0.0a=214286.675（2）凸轮机构廓线电算结果凸轮理论廓线角度/degx/mmy/mm0x=0.000y=35.0005x=3.061y=34.98610x=6.159y=34.93015x=9.324y=34.79820x=12.569y=34.53425x=15.888y=34.07230x=19.250y=33.34235x=22.602y=32.27940x=25.872y=30.83345x=28.974y=28.97450x=31.815y=26.69655x=34.307y=24.02260x=36.373y=21.00065x=38.065y=17.75070x=39.355y=14.32475x=40.116y=10.74980x=40.352y=7.11585x=40.097y=3.50890x=39.406y=0.00095x=38.353y=-3.355100x=37.023y=-6.528105x=35.498y=-9.512110x=33.853y=-12.321115x=32.146y=-14.990120x=30.414y=-17.560125x=28.670y=-20.075130x=26.812y=-22.498135x=24.749y=-24.749140x=22.498y=-26.812145x=20.075y=-28.670150x=17.500y=-30.311155x=14.792y=-31.721160x=11.971y=-32.889165x=9.059y=-33.807170x=6.078y=-34.468175x=3.050y=-34.867180x=0.000y=-35.000185x=-3.050y=-34.867190x=-6.078y=-34.468195x=-9.059y=-33.807200x=-11.971y=-32.889205x=-14.792y=-31.721210x=-17.500y=-30.311215x=-20.075y=-28.670220x=-22.498y=-26.812225x=-24.749y=-24.749230x=-26.812y=-22.498235x=-28.670y=-20.075240x=-30.311y=-17.500245x=-31.721y=-14.792250x=-32.889y=-11.971255x=-33.807y=-9.059260x=-34.468y=-6.078265x=-34.867y=-3.050270x=-35.000y=-0.000275x=-34.867y=3.050280x=-34.468y=6.078285x=-33.807y=9.059290x=-32.889y=11.971295x=-31.721y=14.792300x=-30.311y=17.500305x=-28.670y=20.075310x=-26.812y=22.498315x=-24.749y=24.749320x=-22.498y=26.812325x=-20.075y=28.670330x=-17.500y=30.311335x=-14.792y=31.721340x=-11.971y=32.889345x=-9.059y=33.807350x=-6.078y=34.468355x=-3.050y=34.867360x=-0.000y=35.000凸轮实际廓线角度/degx1/mmy1/mm0x1=0.000y1=25.0005x1=2.960y1=24.98610x1=5.883y1=24.93415x1=8.739y1=24.81520x1=11.506y1=24.59025x1=14.168y1=24.22130x1=16.705y1=23.67135x1=19.089y1=22.91740x1=21.290y1=21.94545x1=23.277y1=20.75550x1=25.021y1=19.35855x1=26.504y1=17.76860x1=27.713y1=16.00065x1=29.002y1=13.52470x1=29.756y1=11.52075x1=30.208y1=9.39380x1=30.352y1=7.16285x1=30.188y1=4.85390x1=29.721y1=2.49395x1=28.972y1=0.107100x1=27.967y1=-2.286105x1=26.743y1=-4.679110x1=25.339y1=-7.075115x1=23.799y1=-9.482120x1=22.165y1=-11.906125x1=20.479y1=-14.339130x1=19.151y1=-16.070135x1=17.678y1=-17.678140x1=16.070y1=-19.151145x1=14.339y1=-20.479150x1=12.500y1=-21.651155x1=10.565y1=-22.658160x1=8.551y1=-23.492165x1=6.470y1=-24.148170x1=4.341y1=-24.620175x1=2.179y1=-24.905180x1=0.000y1=-25.000185x1=-2.179y1=-24.905190x1=-4.341y1=-24.620195x1=-6.470y1=-24.148200x1=-8.551y1=-23.492205x1=-10.565y1=-22.658210x1=-12.500y1=-21.651215x1=-14.339y1=-20.479220x1=-16.070y1=-19.151225x1=-17.678y1=-17.678230x1=-19.151y1=-16.070235x1=-20.479y1=-14.339240x1=-21.651y1=-12.500245x1=-22.658y1=-10.565250x1=-23.492y1=-8.551255x1=-24.148y1=-6.470260x1=-24.620y1=-4.341265x1=-24.905y1=-2.179270x1=-25.000y1=-0.000275x1=-24.905y1=2.179280x1=-24.620y1=4.341285x1=-24.148y1=6.470290x1=-23.492y1=8.551295x1=-22.658y1=10.565300x1=-21.651y1=12.500305x1=-20.479y1=14.339310x1=-19.151y1=16.070315x1=-17.678y1=17.678320x1=-16.070y1=19.151325x1=-14.339y1=20.479330x1=-12.500y1=21.651335x1=-10.565y1=22.658340x1=-8.551y1=23.492345x1=-6.470y1=24.148350x1=-4.341y1=24.620355x1=-2.179y1=24.905360x1=-0.000y1=25.000四、相关程序及源代码（1）连杆机构的matlab源代码for i=0:360;w=68.0678;w0=2*3.141592654/(360*w);t=i*w0;c=w*t;m= sqrt(185*185-92.5*92.5*sin(c)*sin(c);y=-92.5*cos(c)+sqrt(m);v=92.5*w*sin(c)-92.5*92.5*w*sin(2*c)/(2*sqrt(m);a=w*w*92.5*cos(c)-(92.5*w)*(92.5*w)*cos(2*c)/sqrt(185*185-92.5*92.5*sin(c)-0.25*(92.5*92.5*w)*(92.5*92.5*w)*sin(2*c)*sin(2*c)/(m*m*m);plot(c,v);fprintf(i=%.1f,y=%.3f,v=%.1f,a=%.3fn,i,y,v,a);hold on;end（2）凸轮轮廓曲线matlab源代码1、理论廓线源代码for t=0:0.1:60;x=(38.5-3.5*cos(3*(t*pi/180)*sin(t*pi/180);y=(38.5-3.5*cos(3*(t*pi/180)*cos(t*pi/180);a=10.5*(sin(3*(t*pi/180)*sin(t*pi/180)+(38.5-3.5*(cos(3*(t*pi/180)*cos(t*pi/180);b=10.5*(sin(3*(t*pi/180)*cos(t*pi/180)-(38.5-3.5*(cos(3*(t*pi/180)*sin(t*pi/180);x1=x+10*b/sqrt(a*a+b*b);y1=y-10*a/sqrt(a*a+b*b);plot(x,y);fprintf(x=%.3f,y=%.3fn,x,y);hold on;endfor t=60:0.1:125;x=(38.5+3.5*cos(3*(t*pi/180-65*pi/180)*sin(t*pi/180);y=(38.5+3.5*cos(3*(t*pi/180-65*pi/180)*cos(t*pi/180);a=-10.5*(sin(3*(t*pi/180-65*pi/180)*sin(t*pi/180)+(38.5+3.5*(cos(3*(t*pi/180-65*pi/180)*cos(t*pi/180);b=-10.5*(sin(3*(t*pi/180-65*pi/180)*cos(t*pi/180)-(38.5+3.5*(cos(3*(t*pi/180-65*pi/180)*sin(t*pi/180);x1=x+10*b/sqrt(a*a+b*b);y1=y-10*a/sqrt(a*a+b*b);plot(x,y);fprintf(x=%.3f,y=%.3fn,x,y);hold on;endfor t=125:0.1:360;x=35*sin(t*pi/180);y=35*cos(t*pi/180);a=35*cos(t*pi/180);b=(-35)* sin(t*pi/180);x1=x+10*b/sqrt(a*a+b*b);y1=y-10*a/sqrt(a*a+b*b);plot(x,y);fprintf(x=%.3f,y=%.3fn,x,y);hold on;end2、实际廓线源代码for t=0:0.1:60;x=(38.5-3.5*cos(3*(t*pi/180)*sin(t*pi/180);y=(38.5-3.5*cos(3*(t*pi/180)*cos(t*pi/180);a=10.5*(sin(3*(t*pi/180)*sin(t*pi/180)+(38.5-3.5*(cos(3*(t*pi/180)*cos(t*pi/180);b=10.5*(sin(3*(t*pi/180)*cos(t*pi/180)-(38.5-3.5*(cos(3*(t*pi/180)*sin(t*pi/180);x1=x+10*b/sqrt(a*a+b*b);y1=y-10*a/sqrt(a*a+b*b);plot(x1,y1);fprintf(x1=%.3f,y1=%.3fn,x1,y1);hold on;endfor t=60:0.1:125;x=(38.5+3.5*cos(3*(t*pi/180-65*pi/180)*sin(t*pi/180);y=(38.5+3.5*cos(3*(t*pi/180-65*pi/180)*cos(t*pi/180);a=-10.5*(sin(3*(t*pi/180-65*pi/180)*sin(t*pi/180)+(38.5+3.5*(cos(3*(t*pi/180-65*pi/180)*cos(t*pi/180);b=-10.5*(sin(3*(t*pi/180-65*pi/180)*cos(t*pi/180)-(38.5+3.5*(cos(3*(t*pi/180-65*pi/180)*sin(t*pi/180);x1=x+10*b/sqrt(a*a+b*b);y1=y-10*a/sqrt(a*a+b*b);plot(x1,y1);fprintf(x1=%.3f,y1=%.3fn,x1,y1);hold on;endfor t=125:0.1:360;x=35*sin(t*pi/180);y=35*cos(t*pi/180);a=35*cos(t*pi/180);b=(-35)* sin(t*pi/180);x1=x+10*b/sqrt(a*a+b*b);y1=y-10*a/sqrt(a*a+b*b);plot(x1,y1);fprintf(x1=%.3f,y1=%.3fn,x1,y1);hold on;end（3）凸轮轮廓曲线c+源代码#include using namespace std;#include #include#includeint main() fstream datafile; double pi,s0,h,e,p0,w,rr; s0=35; pi=3.141592654; h=7.0; rr=7.0; p0=60.0; w=325*2*pi/60; int p=0,p3=0; double x=0.0,y=0.0,x1=0.0,y1=0.0,dx=0.0,dy=0.0; double s=0.0,v=0.0,a=0.0,v1=0,a1=0; for(p=0;p=0.0)&(p=60.0)&(p=65.0)&(p=125.0)&(p=360.0) s=0;v=0;a=0; v1=v/w; a1=a/(w*w); x=-(s0+s)*sin(p*2*pi/360.0); y=(s0+s)*cos(p*2*pi/360.0); dx=-(v/w-e)*sin(p*2*pi/360.0)+(s0+s)*cos(p*2*pi/360.0); dy=(v/w-e)*cos(p*2*pi/360.0)+(s0+s)*sin(p*2*pi/360.0); x1=x+rr*dy/sqrt(dx*dx)+(dy*dy); y1=y-rr*dx/sqrt(dx*dx)+(dy*dy); datafile.precision(2); datafile.setf(ios:fixed|ios:showpoint); datafile.open(x.txt,ios:out|ios:app); datafilexendl; datafile.close( ); datafile.open(y.txt,ios:out|ios:app); datafileyendl; datafile.close( ); datafile.open(x1.txt,ios:out|ios:app); datafilex1endl; datafile.close(); datafile.open(y1.txt,ios:out|ios:app);