Matlab凸轮位移分析设计.docVIP

MATLAB机械工程 凸轮位移分析的实例 班 级: 姓 名: 学 号: 指导教师: 一．实验目的 通过MATLAB上机编程，掌握MATLAB语言变量、函数、MATLAB常用程序结构和流程控制语句、交互控制语句、程序调试运行等。?? 二．实验内容? 凸轮是把一种运动转化为另一种运动的装置，通常凸轮做旋转运动，并转化为从动件的振动、直线运动或者二者的结合。对基圆为 QUOTE 的凸轮， QUOTE 为最小圆周半径，与凸轮表面相切，并和凸轮的转轴同心。则从动件的运动方程为：? QUOTE ??? 对于旋转运动，假定有：? QUOTE QUOTE QUOTE QUOTE QUOTE QUOTE ?? 其中，h为从动件的最大位移， QUOTE 为凸轮转角， QUOTE 为从动件摆角，并且 QUOTE 。如果凸轮的角速度 QUOTE 为常量，试求出当 QUOTE 时从动件的位移s、速度v、加速度a、瞬时加速度j的变化曲线。? ? 三、建立M文件如下： beat=60*pi/180; phi=linspace(0,beat,40); phi2=[beat+phi]; ph=[phi phi2]*180/pi; arg=2*pi*phi/beat; arg2=2*pi*(phi2-beat)/beat; s=[phi/beat-sin(arg)/2/pil-(arg2-sin(arg2))/2/pi]; v=[1-cos(arg)/beat-(1-cos(arg2))/beat]; a=[2*pi/beat^2*sin(arg)2*pi/beat^2*sin(arg2)]; j=[4*pi^2/beat^3*cos(arg)4*pi^2/beat^3*cos(arg2)]; subplot(2,2,1); plot(ph,s,′k′); xlabel(′凸轮转角（度）′)； ylabel(′位移(s)′); g=axis;g(2)=120;axis(g); subplot(2,2,2); plot(ph,v,′k′,[0 120],[0 0],′k--′) xlabel(′凸轮转角（度）′)； ylabel(′速度（v）′)； g=axis;g(2)=120;axis(g); subplot(2,2,3); polt(ph,a,′k′,[0 120],[0 0],′k--′) xlabel(′凸轮转角（度）′)； ylabel(′加速度（v）′)； g=axis;g(2)=120;axis(g); subplot(2,2,4); polt(ph,j,′k′,[0 120],[0 0],′k--′) xlabel(′凸轮转角（度）′)； ylabel(′瞬时加速度（j）′)； g=axis;g(2)=120;axis(g); 四，输出图像