哈工大机械原理大作业凸轮机构
哈工大、机械原理大作业、凸轮机构设计20题
Harbin Institute of Technology机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：能源科学与工程学院班级：1102301设计者：刘平成学号：1110200724指导教师：唐德威设计时间：2013年6月7日凸轮机构设计1．设计题目(1) 凸轮机构运动简图：(2)凸轮机构的原始参数表2-1.凸轮机构原始参数 序号 升程（mm ）升程运动角 升程运动规律升程许用压力角20 110 120° 正弦加速度35°回程运动角回程运动规律 回程许用压力角 远休止角近休止角 90°正弦加速度 65°90°60°（二）凸轮运动方程及相关图像、程序凸轮推杆升程、回程运动方程及推杆位移、速度、加速度线图： ○1 凸轮推杆升程、回程方程 πϕπϕϕs)650(πϕ≤≤140)(2=ϕs 511()69πφπ≤≤pi))*5708)/(23.2289)/1.-(sin(2+57083.2289)/1.-(-140(1)(3ϕπϕϕ=s1116()99πφπ≤≤)2914(πϕπ≤≤ 0)(4=ϕs ○2速度方程/2.0944;/2.09440))cos(2-140(1)(1πϕϕ=v 16(2)9πφπ≤≤ 0)(2=ϕv 511()69πφπ≤≤ 708;5708))/1.53.2289)/1.-(cos(2-140(1)(3ϕπϕ=v 1116()99πφπ≤≤ 0)(4=ϕv 16(2)9πφπ≤≤○3加速度方程 .0944^2;/2.0944)/2sin(2280)(1πϕπϕ=a )650(πϕ≤≤0)(2=ϕa 511()69πφπ≤≤08^25708)/1.573.2289)/1.-(sin(2280)(3ϕππϕ=a 1116()99πφπ≤≤ 0)(4=ϕa 16(2)9πφπ≤≤推杆位移、速度、加速度线图matlab编程clear,clcpu=0*pi/180:0.0001:120*pi/180; %升程运动角范围pf=120*pi/180:0.0001:210*pi/180; %远休止角范围pd=210*pi/180:0.0001:300*pi/180; %回程运动角范围pn=300*pi/180:0.0001:2*pi; %近休止角范围h=110e-3; %升程w=10; %凸轮角速度p0=120*pi/180; %升程运动角p01=90*pi/180; %回程运动角ps=90*pi/180; %远休止角%----------推程-----------------------------------------su=h.*(pu./p0-sin(2.*pi.*pu./p0)/(2*pi)); %推杆位移vu=h*w/p0*(1-cos(2*pi*pu./p0)); %推程速度au=2*pi*h*w^2/p0^2*sin(2*pi*pu./p0); %推程加速度%------------远休止----------------------------nf=size(pf);sf=h*ones(nf); %推杆位移vf=zeros(nf); %推程速度af=zeros(nf); %推程加速度%---------------回程------------------------------T=pd-(p0+ps);sd=h/2*(1+cos(pi/p01*T)); %回程位移vd=-pi*h*w/(2*p01)*sin(pi/p01*T); %回程速度ad=-pi^2*h*w^2/(2*p01^2)*cos(pi/p01*T); %回程加速度%--------------------近休止---------------------------------nn=size(pn);sn=zeros(nn); %推杆位移vn=zeros(nn); %推程速度an=zeros(nn); % 推程加速度%------画出推杆位移、速度、加速度线图---------------p=[pu,pf,pd,pn];s=[su,sf,sd,sn];subplot(2,3,1),hold onplot(p,s*1e3,'linewidth',2),xlabel('\phi/rad'),ylabel('s/mm'),grid on,title('推杆位移'),axis([0,2*pi,1.1*min(s)*1e3,1.1*max(s)*1e3]) subplot(2,3,2) v=[vu,vf,vd,vn];plot(p,v,'linewidth',2),xlabel('\phi/rad'),ylabel('v/m/s'),grid on,title('推杆速度'),axis([0,2*pi,1.1*min(v),1.1*max(v)]) subplot(2,3,3) a=[au,af,ad,an];plot(p,a,'linewidth',2),xlabel('\phi/rad'),ylabel('a/m/s^2'),grid on,title('推杆加速度'),axis([0,2*pi,1.1*min(a),1.1*max(a)]) hold off（三）凸轮机构s d ds-ϕ图像及程序代码 %--------------------求ds/d_phi-------------------subplot(2,3,4),plot(v/w*1e3,s*1e3,'linewidth',2),xlabel('ds/d\phi/mm'),ylabel('s/mm'),axis equal,grid on,title('ds/d\phi —s')%---------------------凸轮轴心许用区域--------------------------- alpha_up=35*pi/180; %升程许用压力角 alpha_down=65*pi/180; %回程许用压力角 p1=pi/2-alpha_up; %推程斜率角 p2=alpha_down-pi/2; %回程斜率角 ku=tan(p1); %推程切线斜率 kd=tan(p2); %回程切线斜率 R2=[cos(-p2),-sin(-p2);sin(-p2),cos(-p2)];%推程旋转矩阵 R1=[cos(-p1),-sin(-p1);sin(-p1),cos(-p1)];%推程旋转矩阵 nu=size(pu); for i=1:nu(2)Temp=R1*[vu(i)/w;su(i)];vut(i)=Temp(1); %旋转推程ds/dp-s 曲线 sut(i)=Temp(2); endnd=size(pd); for i=1:nd(2)Temp=R2*[vd(i)/w;sd(i)];vdt(i)=Temp(1); %旋转回程ds/dp-s 曲线 sdt(i)=Temp(2); endfor j=1:nu(2)if sut(j)==min(sut)temu=j; %旋转推程ds/dp-s 曲线后求最低点 end endfor j=1:nd(2)if sdt(j)==min(sdt)temd=j; %旋转回程ds/dp-s曲线后求最低点endendt1=1.2*min(vd/w):0.01:1.2*max(vu/w); %切线定义域t2=min(vd/w)/6:0.01:1.2*max(vu/w);t3=0:0.01:1.2*max(vu/w);s1=ku*(t2-vu(temu)/w)+su(temu); %推程切线s2=kd*(t1-vd(temd)/w)+sd(temd); %回程切线s3=tan(-p1)*t3; %推程起点压力角限制线subplot(2,3,5) %画图hold on,axis equal,grid onplot(v/w*1e3,s*1e3,'linewidth',2)plot(t2*1e3,s1*1e3,'linewidth',1,'color','r')plot(t1*1e3,s2*1e3,'linewidth',1,'color','r')plot(t3*1e3,s3*1e3,'linewidth',1,'color','r')xlabel('ds/d\phi/mm'),ylabel('s/mm'),hold off,title('ds/d\phi—s，轴向许用范围')（四）确定凸轮的基圆半径和偏距、绘制凸轮机圆、偏距圆、理论轮廓曲线---------------画理论廓线图-------------------------e=36e-3;s0=52e-3;r0=sqrt(s0^2+e^2);x=(s0+s).*cos(p)-e.*sin(p);y=(s0+s).*sin(p)+e.*cos(p);x1=r0*cos(p);y1=r0*sin(p);subplot(2,3,6)plot(x*1e3,y*1e3,'linewidth',1),axis equal,grid on,hold on,title('廓线图')plot(x1*1e3,y1*1e3,'linewidth',1,'color','r')%-------------求最小曲率半径-----------------------nx=size(x);nx1=nx(2)-2;dydp=diff(y)./diff(p);%求微分dxdp=diff(x)./diff(p);d2ydp2=diff(dydp)./diff(p(1:nx1+1));d2xdp2=diff(dxdp)./diff(p(1:nx1+1));rho=(dxdp(1:nx1).^2+dydp(1:nx1).^2).^1.5./abs((dxdp(1:nx1).*d2ydp2(1:nx1)-dydp(1:nx1).*d2x dp2(1:nx1)));%理论廓线曲率半径rhomin=min(rho);%最小曲率半径rr=rhomin-3e-3;%----------------实际廓线图----------------X=x(1:nx(2)-1)-rr*dydp./(dxdp.^2+dydp.^2).^0.5;%求实际廓线坐标Y=y(1:nx(2)-1)+rr*dxdp./(dxdp.^2+dydp.^2).^0.5;plot(X*1e3,Y*1e3,'linewidth',2,'color','k')%画实际廓线图Legend('理论廓线','基圆','实际廓线'),axis([1.1*min(x)*1e3,1.1*max(x)*1e3,1.1*min(y)*1e3,1.1*max(y)*1e3])得到基圆半径311mm、偏距36mm。
哈工大机械原理大作业凸轮机构设计第题
哈工大机械原理大作业-凸轮机构设计(第题)————————————————————————————————作者：————————————————————————————————日期：机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：机电学院班级：1208103完成者：xxxxxxx学号：11208103xx指导教师：林琳设计时间：2014.5.2哈尔滨工业大学凸轮机构设计一、设计题目如图所示直动从动件盘形凸轮机构，其原始参数见表，据此设计该凸轮机构。
序号 升程（mm ） 升程运动角（°） 升程运动规律 升程许用压力角（°） 回程运动角（°） 回程运动规律 回程许用压力角（°）远休止角（°） 近休止角（°）3 50 150 正弦加速度 30 100 余弦加速度60 30 80二、凸轮推杆升程、回程运动方程及其线图1 、凸轮推杆升程运动方程(650πϕ≤≤) 升程采用正弦加速度运动规律，故将已知条件mm h 50=，650π=Φ带入正弦加速度运动规律的升程段方程式中得：⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-=512sin 215650ϕππϕS ； ⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-=512cos 1601ππωv ； ω⎪⎭⎫ ⎝⎛=512sin 14421ϕπωa ； 2、凸轮推杆推程远休止角运动方程（πϕπ≤≤65） mm h s 50==；0==a v ；3、凸轮推杆回程运动方程（914πϕπ≤≤） 回程采用余弦加速度运动规律，故将已知条件mm h 50=，95'0π=Φ，6s π=Φ带入余弦加速度运动规律的回程段方程式中得：⎥⎦⎤⎢⎣⎡-+=)(59cos 125πϕs ； ()πϕω--=59sin451v ； ()πϕω-=59cos 81-a 21；4、凸轮推杆回程近休止角运动方程（πϕπ2914≤≤） 0===a v s ；5、凸轮推杆位移、速度、加速度线图根据以上所列的运动方程，利用matlab 绘制出位移、速度、加速度线图。
哈工大机械原理大作业直动从动件盘形凸轮机构满分完美版哈尔滨工业大学
end
%曲率半径
dx=diff(x);
dx(36001)=0;
dy=diff(y);
dy(36001)=0;
dydx=dy./dx;
ddy=diff(dydx);
ddy(36001)=0;
ddy=ddy./dx;
~7~
xlabel('凸轮转角φ/°');
ylabel('从动件加速度 a/(mm/s^2)');
title('加速度');
subplot(2,3,4);
plot(dsdPhi,s);
axis equal
axis([-30 50 -50 30]);
hold on
plot(x1,f1,'r');
hold on
axis equal
hold on
plot(x,y,'k')
legend('理论廓线',-1);
hold on
%轮廓图
%实际廓线
%理论廓线
~8~
哈尔滨工业大学
直动从动件盘型凸轮机构设计说明书
plot(x0,y0,':')
hold on
%基圆
plot(xe,ye,'k:')
%偏距圆
legend('实际廓线','理论廓线','基圆','偏距圆',-1);
subplot(2,3,1);
plot(Phi,s);
grid on
axis([0 360 0 27]);
哈工大机械原理大作业——凸轮——5号
Harbin Institute of Technology机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：机电工程学院班级：设计者：学号：指导教师：设计时间：凸轮机构设计一． 设计题目表1 凸轮机构原始参数二． 凸轮推杆升程、回程运动方程及推杆位移、速度、加速度线图1. 推杆升程运动方程（ϕ为转角）（203ϕπ≤≤） 令1/rad s ω= 位移：001212040[sin()]sin(3)2s h ϕπϕϕϕφπφππ=-=- 速度：1002120120[1cos()]cos(3)h v ωπϕϕφφππ=-=- 加速度：2120022360[sin()]sin(3)h a πωπϕϕφφπ==2. 推杆回程运动方程（ϕ为转角）7563πϕπ≤≤ 位移：''00128001604014[1sin()]sin(4)]233T s h T πϕϕπφπφππ=-+=-+- 速度：1''00216016014[1cos()]cos(4)]3h v T ωπϕϕπφφππ=--=-+- 加速度：21'2'002264014sin()sin(4)]3h a T πωπϕπφφπ=-=-- 式中：027()()326s T πππϕϕϕϕϕ=-+=-+=-由MATLAB编程得到线位移图像：线速度图像：线加速度图像：三．绘制机构的ds sdϕ-线图由dsds vdtdd wtϕϕ==可知线图即vsw-线图，由MATLAB编程后得图四．确定基圆半径和偏距1.画切线图如上页图，得其切线图2.求基圆半径和偏距从图中读取e=50mmr0=111.80mm五．画理论及实际轮廓附：MATLAB编程1.推杆位移线图代码%升程阶段t=0:0.001:2*pi/3;s=120.*t./pi-40./pi.*sin(3.*t);hold onplot(t,s);%远休止阶段t=2*pi/3:0.001:7*pi/6;s=80;hold onplot(t,s);%回程阶段t=7*pi/6:0.001:5*pi/3;s=800/3-160.*t./pi+40./pi.*sin(4.*t-14*pi/3) hold onplot(t,s);%近休止阶段t=5*pi/3:0.001:2*pi;s=0;hold onplot(t,s);grid onhold off2.推杆速度线图代码%升程阶段t=0:0.001:2*pi/3;v=120/pi-120/pi*cos(3.*t);hold onplot(t,s);%远休止阶段t=2*pi/3:0.001:7*pi/6;v=0;hold onplot(t,s);%回程阶段t=7*pi/6:0.001:5*pi/3;v=-160/pi+160/pi*cos(4.*t-14/3*pi) hold onplot(t,s);%近休止阶段t=5*pi/3:0.001:2*pi;v=0;hold onplot(t,s);grid onhold off3.推杆加速度线图代码%升程阶段t=0:0.001:2*pi/3;a=360/pi.*sin(3.*t)hold onplot(t,a);%远休止阶段t=2*pi/3:0.001:7*pi/6;a=0;hold onplot(t,a);%回程阶段t=7*pi/6:0.001:5*pi/3;a=-640/pi.*sin(4.*t-14/3*pi) hold onplot(t,a);%近休止阶段t=5*pi/3:0.001:2*pi;a=0;hold onplot(t,a);grid onhold off4.dssdϕ-线图代码%升程阶段t=0:0.001:2*pi/3;s=120.*t./pi-40./pi.*sin(3.*t);v=120/pi-120/pi*cos(3.*t);hold onplot(v,s);%远休止阶段t=2*pi/3:0.001:7*pi/6;s=80;v=0;hold onplot(v,s);%回程阶段t=7*pi/6:0.001:5*pi/3;s=800/3-160.*t./pi+40./pi.*sin(4.*t-14*pi/3);v=-160/pi+160/pi*cos(4.*t-14/3*pi)hold onplot(v,s);5.最终轮廓线图代码h=80;w=1;e=50;rr=20;s0=100;m=120*pi/180;ms=(120+90)*pi/180;m1=(120+90+90)*pi/180; for i=1:1:120mm(i)=i*pi/180.0;s1=h.*( mm(i)./m -(1/(2*pi)).*sin(2.*pi.*mm(i)./m));v1=(w*h./m)*(1-cos(2*pi*mm(i)/m));x(i)=(s0+s1)*sin(mm(i))+e*cos(mm(i));y(i)=(s0+s1)*cos(mm(i))-e*sin(mm(i));a(i)=(s0+s1)*cos(mm(i))-e*sin(mm(i))+v1/w*sin(mm(i)); b(i)=-(s0+s1)*sin(mm(i))-e*cos(mm(i))+v1/w*cos(mm(i));xx(i)=x(i)+rr*b(i)/sqrt(a(i)*a(i)+b(i)*b(i));yy(i)=y(i)-rr*a(i)/sqrt(a(i)*a(i)+b(i)*b(i));endfor i=121:1:210mm(i)=i*pi/180;s2=h;v2=0;x(i)=(s0+s2)*sin(mm(i))+e*cos(mm(i));y(i)=(s0+s2)*cos(mm(i))-e*sin(mm(i));a(i)=(s0+s2)*cos(mm(i))-e*sin(mm(i))+v2/w*sin(mm(i));b(i)=-(s0+s2)*sin(mm(i))-e*cos(mm(i))+v2/w*cos(mm(i));xx(i)=x(i)+rr*b(i)/sqrt(a(i)*a(i)+b(i)*b(i));yy(i)=y(i)-rr*a(i)/sqrt(a(i)*a(i)+b(i)*b(i));endfor i=211:1:300mm(i)=i*pi/180;s3=800/3-160/pi.*mm(i)+40/pi.*sin(4.*mm(i)-14/3*pi);v3=-160/pi+160/pi.*cos(4*mm(i)-14/3*pi)x(i)=(s0+s3)*sin(mm(i))+e*cos(mm(i));y(i)=(s0+s3)*cos(mm(i))-e*sin(mm(i));a(i)=(s0+s3)*cos(mm(i))-e*sin(mm(i))+v3/w*sin(mm(i));b(i)=-(s0+s3)*sin(mm(i))-e*cos(mm(i))+v3/w*cos(mm(i));xx(i)=x(i)+rr*b(i)/sqrt(a(i)*a(i)+b(i)*b(i));yy(i)=y(i)-rr*a(i)/sqrt(a(i)*a(i)+b(i)*b(i));endfor i=301:1:360mm(i)=i*pi/180;x(i)=(s0+0)*sin(mm(i))+e*cos(mm(i));y(i)=(s0+0)*cos(mm(i))-e*sin(mm(i));a(i)=(s0+0)*cos(mm(i))-e*sin(mm(i))+v3/w*sin(mm(i));b(i)=-(s0+0)*sin(mm(i))-e*cos(mm(i))+v3/w*cos(mm(i));xx(i)=x(i)+rr*b(i)/sqrt(a(i)*a(i)+b(i)*b(i));yy(i)=y(i)-rr*a(i)/sqrt(a(i)*a(i)+b(i)*b(i));endplot(x,y,'r',xx,yy,'g')text(0,20,'实际轮廓线')text(120,100,'理论轮廓线')hold on。
哈工大机械原理大作业-凸轮机构设计
哈工大机械原理大作业-凸轮机构设计(第3题)(共15页)-本页仅作为预览文档封面，使用时请删除本页-机械原理大作业二课程名称：机械原理设计题目：凸轮设计院系：机电学院班级： 1208103完成者： xxxxxxx学号： xx指导教师：林琳设计时间：工业大学凸轮设计一、设计题目如图所示直动从动件盘形凸轮，其原始参数见表，据此设计该凸轮。
二、凸轮推杆升程、回程运动方程及其线图1 、凸轮推杆升程运动方程(650πϕ≤≤) 升程采用正弦加速度运动规律，故将已知条件mm h 50=，650π=Φ带入正弦加速度运动规律的升程段方程式中得：⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-=512sin 215650ϕππϕS ；⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-=512cos 1601ππωv ； ⎪⎭⎫ ⎝⎛=512sin 14421ϕπωa ； 2、凸轮推杆推程远休止角运动方程（πϕπ≤≤65） mm h s 50==； 0==a v ；3、凸轮推杆回程运动方程（914πϕπ≤≤）回程采用余弦加速度运动规律，故将已知条件mm h 50=，95'0π=Φ，6s π=Φ带入余弦加速度运动规律的回程段方程式中得：⎥⎦⎤⎢⎣⎡-+=)(59cos 125πϕs ；()πϕω--=59sin451v ； ()πϕω-=59cos 81-a 21；4、凸轮推杆回程近休止角运动方程（πϕπ2914≤≤） 0===a v s ；5、凸轮推杆位移、速度、加速度线图根据以上所列的运动方程，利用matlab 绘制出位移、速度、加速度线图。
①位移线图 编程如下： %用t 代替转角 t=0::5*pi/6;s=50*((6*t)/(5*pi)-1/(2*pi)*sin(12*t/5)); hold on plot(t,s); t=5*pi/6::pi; s=50; hold on plot(t,s); t=pi::14*pi/9;s=25*(1+cos(9*(t-pi)/5));hold onplot(t,s);t=14*pi/9::2*pi;s=0;hold onplot(t,s),xlabel('φ/rad'),ylabel('s/mm'); grid onhold off所得图像为：②速度线图编程如下：%用t代替转角,设凸轮转动角速度为1t=0::5*pi/6;v=60/pi*(1-cos((12*t)/5));hold onplot(t,v);t=5*pi/6::pi;v=0;hold onplot(t,v);t=pi::14*pi/9;v=-45*sin(9*(t-pi)/5);hold onplot(t,v);t=14*pi/9::2*pi;v=0;hold onplot(t,v),xlabel('φ(rad)'),ylabel('v(mm/s)'); grid onhold off所得图像为：③加速度线图利用matlab编程如下：%用t代替转角,设凸轮转动角速度为1t=0::5*pi/6;a=144/pi*sin(12*t/5);hold onplot(t,a);t=5*pi/6::pi;a=0;hold onplot(t,a);t=pi::14*pi/9;a=-81*cos(9*(t-pi)/5);hold onplot(t,a);t=14*pi/9::2*pi; a=0; hold onplot(t,a),xlabel('φ(rad)'),ylabel('a(mm/s^2)'); grid on hold off所得图形：三、绘制s d ds -ϕ线图根据运动方程求得：()⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧≤≤≤≤-≤≤≤≤--=πϕππϕππϕπϕππϕπππϕ2914.0914,59sin 4565,0650),512cos 6060(d ds 利用matlab 编程：%用t 代替φ,a 代替ds/d φ, t=0::5*pi/6;a=-(60/pi-60/pi*cos(12*t/5));s=50*((6*t)/(5*pi)-1/(2*pi)*sin(12*t/5)); hold on plot(a,s); t=5*pi/6::pi; a=0; s=50; hold on plot(a,s); t=pi::14*pi/9;a=45*sin(9*(t-pi)/5); s=25*(1+cos(9*(t-pi)/5)); hold on plot(a,s);t=14*pi/9::2*pi; a=0; s=0; hold onplot(a,s),title('ds/d φ-s'),xlabel('ds/d φ(mm/rad)'),ylabel('s(mm)'); grid on hold off 得s d ds-ϕ图：凸轮压力角的正切值s s e d ds +-=0/tan ϕα，左侧为升程，作与s 轴夹6π角等于升程许用压力角的切界线t t d D ，则在直线上或其左下方取凸轮轴心时，可使[]αα≤，同理右侧回程，作与s 轴夹角等于回程许用压力角3π的切界线''t t d D ，则在直线上或其右下方取凸轮轴心时，可使[]αα≤。
哈工大机械原理大作业——凸轮——2号
哈工大机械原理大作业——凸轮——2号————————————————————————————————作者：————————————————————————————————日期：Harbin Institute of Technology机械原理大作业课程名称：机械原理设计题目：凸轮机构设计一、设计题目(1)凸轮机构运动简图：(2)凸轮机构的原始参数序号升程升程运动角升程运动规律升程许用压力角回程运动角回程运动规律回程许用压力角远休止角近休止角14 90°120°余弦加速度35°90°3-4-5多项式65°80°70°(1) 推杆升程、回程运动方程如下：A.推杆升程方程：设为1rad sω=升程位移为：()()1cos451cos1.52hsπψψψ⎡⎤⎛⎫=-=-⎢⎥⎪Φ⎝⎭⎣⎦23ψπ≤≤升程速度为：()()1100sin67.5sin1.52hvπωπψψωψ⎛⎫==⎪ΦΦ⎝⎭23ψπ≤≤升程加速度为：()()2221100cos101.25cos1.52haπωπψψωψ⎛⎫==⎪ΦΦ⎝⎭23ψπ≤≤B.推杆回程方程：回程位移为：()()345111110156s h T T T ψ⎡⎤=--+⎣⎦1029918ψπ≤≤ 回程速度为：()()22111103012h v T T T ωψ=--+'Φ 1029918ψπ≤≤ 回程加速度为：()()221111260132h a T T T ωψ=--+'Φ 1029918ψπ≤≤其中：()010s T ψ-Φ+Φ='Φ1029918ψπ≤≤ (2) 利用Matlab 绘制推杆位移、速度、加速度线图 A. 推杆位移线图clcclearx1=linspace(0,2*pi/3,300);x2=linspace(2*pi/3,10*pi/9,300); x3=linspace(10*pi/9,29*pi/18,300); x4=linspace(29*pi/18,2*pi,300); T1=(x3-10*pi/9)/(pi/2); s1=45*(1-cos(1.5*x1)) s2=90;s3=90*(1-(10*T1.^3-15*T1.^4+6*T1.^5)); s4=0;plot(x1,s1,'r',x2,s2,'r',x3,s3,'r',x4,s4,'r') xlabel('角度ψ/rad'); ylabel('位移s/mm') title('推杆位移线图') gridaxis([0,7,-10,100]) 得到推杆位移线图：B.推杆速度线图clcclearx1=linspace(0,2*pi/3,300);x2=linspace(2*pi/3,10*pi/9,300);x3=linspace(10*pi/9,29*pi/18,300);x4=linspace(29*pi/18,2*pi,300);T1=(x3-10*pi/9)/(pi/2);v1=67.5*1*sin(1.5*x1);v2=0;v3=-30*90*1*T1.^2/(pi/2).*(1-2*T1+T1.^2);v4=0;plot(x1,v1,'r',x2,v2,'r',x3,v3,'r',x4,v4,'r') xlabel('角度ψ/rad');ylabel('速度v/(mm/s)')title('推杆速度线图')Grid得到推杆速度线图：C.推杆加速度线图clcclearx1=linspace(0,2*pi/3,300);x2=linspace(2*pi/3,10*pi/9,300);x3=linspace(10*pi/9,29*pi/18,300);x4=linspace(29*pi/18,2*pi,300);T1=(x3-10*pi/9)/(pi/2);a1=101.25*1.^2.*cos(1.5*x1);a2=0;a3=-60.*90.*T1./((pi/2).^2).*(1-3*T1+2*T1.^2); a4=0;plot(x1,a1,'r',x2,a2,'r',x3,a3,'r',x4,a4,'r') xlabel('角度ψ/rad');ylabel('加速度a/')title('推杆加速度线图')Grid得到推杆加速度线图：三、凸轮机构的ds/dψ-s线图，并依次确定凸轮的基圆半径和偏距.1、凸轮机构的ds/dψ--s线图：x1=linspace(0,2*pi/3,300);x2=linspace(2*pi/3,10*pi/9,300);x3=linspace(10*pi/9,29*pi/18,300);x4=linspace(29*pi/18,2*pi,300);T1=(x3-10*pi/9)/(pi/2);s1=45*(1-cos(1.5*x1))s2=90;s3=90*(1-(10*T1.^3-15*T1.^4+6*T1.^5));s4=0;v1=67.5*1*sin(1.5*x1);v2=0;v3=-30*90*1*T1.^2/(pi/2).*(1-2*T1+T1.^2);v4=0;plot(v1,s1,'r',v2,s2,'r',v3,s3,'r',v4,s4,'r')xlabel('ds/dψ');ylabel('(位移s/mm)')title('ds/dψ—s曲线')gridaxis([-120,80,-10,100])得到ds/dψ—s曲线：2、确定凸轮的基圆半径和偏距：在dssdϕ-线图中，右侧曲线为升程阶段的类速度-位移图，作直线D t d t与其相切，且与位移轴正方向呈夹角[α1]=350, 故该直线斜率：32sin2=tan5533cos2okϕϕ⨯=⨯通过编程求其角度。
哈工大机械原理大作业凸轮设计
哈工大机械原理大作业凸轮设计Harb inIn stituteofTech no logy大作业设计说明书课程名称：设计题目：院班学级：机械原理凸轮机构设计1208103系：机械设计制造及其自动化设计指导教师：设计时间：林琳2019425哈尔滨工业大学一、运动分析题目如图所示直动从动件盘形凸轮机构，其原始数据参数见表2-1,。
从表2-1中选择一组凸轮机构原始参数，据此设计该凸轮机构。
二、凸轮运动规律升程运动角（°）90升程运动规律生程许回程运用压力动角角（°）等加等4080减速回程运动规律回程许远休用压力止角角（°）余弦加7040速度近休止角（°）150升程（mm）1501 、升程运动规律（0 /4）位移s=2h（速度v2 /2 4*150*w( /2)A24*150*w A2( /2)人2加速度a2 、升程运动规律(/4 /2)位移s 1502*150( /2 )A2(/2)A2速度v4*60*w( /2 )( /2)A24*60*wA2( /2)A2加速度a3 、回程运动规律(/2 2/2 2 ) 93 位移s 75*{1 cos[ (/2 2 )]}949速度vhw*sin*[ ( /2 2 )]92*4 499加速度aA2hw A2cos[ ( /2 2 )]94 2*(4 )A299根据运动规律做出的曲线以及源代码如图所示位移线图速度线图加速度线图位移线图源代码fl=pi/180;x0=0:fl:pi/4;x1=pi/4:fl:pi/2;x2=pi/2:fl:13*pi/18;x3=13*pi/18:fl:7*pi/6;x4=7*pi/6:fl:2*pi;s0=300*(2*x0/pi)A 2；s1=150-1200*(pi/2-x1).*(pi/2-x1)/(pi.*pi);s2=150+x2*0;s3=75*(1+cos(9/4*(x3-13*pi/18)));s4=x4*0;Plot(x0,s0,x1,s1,'b',x2,s2,'b',x3,s3,'b',x4,s4,'b')axis([070200])title('杆位移线图')xlabel(' 0 (rad)')ylabel('V(mm⑸')gridon速度源代码fl=pi/180;x0=0:fl:pi/4;x1=pi/4:fl:pi/2;x2=pi/2:fl:13*pi/18;x3=13*pi/18: fl:7*pi/6;x4=7*pi/6:fl:2*pi;w=30;v0=600.*w.*x0/(pi/2)A2;v1=600.*w.*(pi/2-x1)/(pi/2)A2;v2=0*x2;v3=-150*30*pi/(2*4*pi/9).*si n(9/4*(x3-13*pi/18));v4=0*x4;Plot(x0,v0,'b',x1,v1,'b',x2,v2,'b',x3,v3,'b',x4,v4,'b')title('推杆速度')xlabel(' 0 (rad)')ylabel('v(mm/s')gridon加速度源代码fl=pi/180;x0=0:fl:pi/4;x1=pi/4:fl:pi/2;x2=pi/2:fl:13*pi/18;x3=13*pi/18: fl:7*pi/6;x4=7*pi/6:fl:2*pi; w=30;a0=600*w.A2/(pi/2).A2+xO*0;a 仁-600*w.A2/(pi/2)A2+x1*0;a2=x2*0; a3=-pi*pi*150*30*30/(2*4*pi/9)A2.*cos(9/4*(x3-13*pi/18));a4=x4*0;Plot(x0,a0,'b',x1,a1,'b',x2,a2,'b',x3,a3,'b',x4,a4,'b')title('推杆加速度')xlabel(' 0 (rad)')ylabel('a(mm/sA2')gridon三、凸轮机构的dss曲线绘制d由凸轮机构位移公式可知4h(/2)A2(0 /4)ds 4*60 ( /2 )( /4 /2)d ( /2)A275*4*s in 9( 13 )( /2 2 /2 2 4 ) 9418999 则其曲线如图所示其源代码如下clcfl=pi/180;x0=0:fl:pi/4;x1=pi/4:fl:pi/2;x2=pi/2:fl:13*pi/18;x3=13*pi/18: fl:7*pi/6;x4=7*pi/6:fl:2*pi;w=30;d0=-600.*x0/(pi/2)A2;d1=-600.*(pi/2-x1)/(pi/2)A2;d2=0*x2;d3=75*4/9.*si n(9/4*(x3-13*pi/18));d4=0*x4;s0=300*(2*x0/pi)A2;s1=150-1200*(pi/2-x1).*(pi/2-x1)/(pi.*pi);s2=150+x2*0;s3=75*(1+cos(9/4*(x3-13*pi/18)));s4=x4*0;Plot(d0,s0,'b',d1,s1,'b',d2,s2,'b',d3,s3,'b',d4,s4,'b')title(' 类速度-位移曲线')xlabel(' 类速度(mm/rad)')ylabel(' 位移(mmm)')gridon四、确定凸轮的基圆半径和偏距以ds/df-s图为基础，可分别作出二条限制线，以这二条线可确定最小基圆半径及所对应的偏距e,在其下方选择一合适点，即可满足压力角的限制条件图像如图所示由图像可知，设置点(50，-100 )为凸轮轴心位置。
哈工大机械原理大作业 凸轮资料
机械原理大作业课程名称：机械原理设计题目：直动从动件盘形凸轮机构院系：机电学院班级：完成者：学号：指导教师：设计时间：哈尔滨工业大学题目：如图所示直动从动件盘形凸轮机构，其原始参数见表，据此设计该凸轮机构。
凸轮运动分为五个阶段 1.升程阶段0~5034500010156s h ϕϕϕ⎡⎤⎛⎫⎛⎫⎛⎫⎢⎥=⨯-⨯+⨯ ⎪ ⎪ ⎪ΦΦΦ⎢⎥⎝⎭⎝⎭⎝⎭⎣⎦ 23410000306030h v ωϕϕϕ⎡⎤⎛⎫⎛⎫⎛⎫⎢⎥=⨯-⨯+⨯ ⎪ ⎪ ⎪ΦΦΦΦ⎢⎥⎝⎭⎝⎭⎝⎭⎣⎦ 23212000060180120h a ωϕϕϕ⎡⎤⎛⎫⎛⎫⎛⎫⎢⎥=⨯-⨯+⨯ ⎪ ⎪ ⎪ΦΦΦΦ⎢⎥⎝⎭⎝⎭⎝⎭⎣⎦2.远休止50~1503.回程等加速150~195()20'202sh s h ϕ=--Φ-ΦΦ()10'204s h s ωϕ=--Φ-ΦΦ21'204h v ω=Φ4.回程等减速195~240()2'00'202s h s ϕ=Φ+Φ+Φ-Φ()'100'24s h v ωϕ=-Φ+Φ+Φ-Φ21'204h v ω=-Φ5.近休止240~360其中，Φ 推程运动角sΦ 远休止角 '0Φ 回程运动角使用Matlab实现1.从动件位移%用fi1,fi2,fi3,fi4，fi5代替转角h=20;fio1=5*pi/18;fis=5*pi/9;fio2=pi/2;fi1=0:0.01:fio1;dfi=fi1/fio1;s1=h*(10*(fi1/fio1).^3-15*(fi1/fio1).^4+6*(fi1/fio1).^5); hold onplot(180*fi1/pi,s1);fi2=fio1:0.01:fio1+fis;s2=h;hold onplot(180*fi2/pi,s2);fi3=fio1+fis:0.01:fio1+fis+fio2/2;s3=h-2*h/fio2^2*(fi3-fio1-fis).^2;hold onplot(180*fi3/pi,s3);fi4=fio1+fis+fio2/2:0.001:fio1+fis+fio2;s4=2*h/fio2^2*(fio1+fis+fio2-fi4).^2;hold onplot(180*fi4/pi,s4);fi5=fio1+fis+fio2:0.001:2*pi;s5=0;hold onplot(180*fi5/pi,s5);title('位移');xlabel('φ/度'),ylabel('s/mm');grid onhold off2.从动件速度%用fi1,fi2,fi3,fi4，fi5代替转角%设角速度为1rad/sw=1;h=20;fio1=5*pi/18;fis=5*pi/9;fio2=pi/2;fi1=0:0.01:fio1;dfi=fi1/fio1;v1=h*w/fio1*(30*(fi1/fio1).^2-60*(fi1/fio1).^3+30*(fi1/fio1).^4); hold onplot(180*fi1/pi,v1);fi2=fio1:0.01:fio1+fis;v2=0;hold onplot(180*fi2/pi,v2);fi3=fio1+fis:0.01:fio1+fis+fio2/2;v3=-4*h*w/fio2^2*(fi3-fio1-fis);hold onplot(180*fi3/pi,v3);fi4=fio1+fis+fio2/2:0.001:fio1+fis+fio2; v4=-4*h*w/fio2^2*(fio1+fis+fio2-fi4); hold onplot(180*fi4/pi,v4);fi5=fio1+fis+fio2:0.001:2*pi;v5=0;hold onplot(180*fi5/pi,v5);title('速度');xlabel('φ/度'),ylabel('v/(mm/s)');grid onhold off3.从动件加速度%用fi1,fi2,fi3,fi4，fi5代替转角%设角速度为1rad/sw=1;h=20;fio1=5*pi/18;fis=5*pi/9;fio2=pi/2;fi1=0:0.01:fio1;dfi=fi1/fio1;a1=h*w^2/fio1*(60*(fi1/fio1)-180*(fi1/fio1).^2+120*(fi1/fio1).^3) ;hold onplot(180*fi1/pi,a1);fi2=fio1:0.01:fio1+fis;a2=0;hold onplot(180*fi2/pi,a2);fi3=fio1+fis:0.01:fio1+fis+fio2/2;a3=-4*h*w/fio2^2;hold onplot(180*fi3/pi,a3);fi4=fio1+fis+fio2/2:0.001:fio1+fis+fio2;a4=4*h*w/fio2^2;hold onplot(180*fi4/pi,a4);fi5=fio1+fis+fio2:0.001:2*pi;a5=0;hold onplot(180*fi5/pi,a5);title('加速度');xlabel('φ/度'),ylabel('v/(mm/s^2)');grid onhold off4.sdsdϕ-线图%用fi1,fi2,fi3,fi4，fi5代替转角h=20;fio1=5*pi/18;fis=5*pi/9;fio2=pi/2;fi1=0:0.01:fio1;s1=h*(10*(fi1/fio1).^3-15*(fi1/fio1).^4+6*(fi1/fio1).^5);ds1=h*(3*10*(fi1/fio1).^2/fio1-4*15*(fi1/fio1).^3/fio1+5*6*(fi1/f io1).^4/fio1);hold onplot(ds1,s1);fi2=fio1:0.01:fio1+fis;s2=h+0*fi2;ds2=0*fi2;hold onplot(ds2,s2);fi3=fio1+fis:0.01:fio1+fis+fio2/2;s3=h-2*h/fio2^2*(fi3-fio1-fis).^2;ds3=-2*h/fio2^2*2*(fi3-fio1-fis);hold onplot(ds3,s3);fi4=fio1+fis+fio2/2:0.001:fio1+fis+fio2; s4=2*h/fio2^2*(fio1+fis+fio2-fi4).^2; ds4=-2*h/fio2^2*2*(fio1+fis+fio2-fi4); hold on plot(ds4,s4);fi5=fio1+fis+fio2:0.001:2*pi; s5=0+0*fi5; ds5=0*fi5; hold on plot(ds5,s5); title('ds/d φ-s');xlabel('ds/d φ(mm/rad)'),ylabel('s(mm)'); grid on hold off5.凸轮轴心位置的确定凸轮压力角的正切值s s e d ds +-=0/tan ϕα，右侧为升程，作与s 轴夹6π角等于升程许用压力角的切界线t t d D ，则在直线上或其下方取凸轮轴心时，可使[]αα≤，同理右侧回程，作与s 轴夹角等于回程许用压力角3π的切界线''t t d D ，则在直线上或其下方取凸轮轴心时，可使[]αα≤。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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凸轮机构设计课程名称：机械原理设计题目：凸轮机构院系：班级：设计者：学号：指导教师：设计时间：哈尔滨工业大学1. 设计题目1) 凸轮机构运动简图：(2)凸轮机构的原始参数序号升程升程运动角升程运动规律 升程许用压力角 回程运动角 回程运动规律 回程许用压力角 远休止角 近休止角 30140mm 150°正弦加速度30°100°等减等加速60°35°75°2. 凸轮推杆升程，回程运动方程及推杆位移,速度,加速度线图 2) 推杆升程，回程运动方程A.推程运动方程（正弦加速度）： (506ϕπ≤≤)()()116870sin 2.4s ϕϕϕππ=-()[]111681cos(2.4)v ωϕϕπ=⨯-()()211403.2sin 2.4a ωϕϕπ⨯=⨯远休止方程：(537636πϕπ≤≤)()2140s ϕ= ()20v ϕ= ()20a ϕ=B.回程运动方程（等减等加速度）： 等减： (37473636πϕπ≤≤)()2312907.21864.89443s piϕϕϕ-π=-+()113121814.41864.8v ωωϕϕ+ππ=-()213121814.4a ωϕπ=-加速： (47193612πϕπ≤≤)()2322907.22872.8s piϕϕϕ+2274.3π=-()113221814.42872.8v ωωϕϕ-ππ=()213221814.4a ωϕπ=近休止方程： (19212πϕπ≤≤)()4140s ϕ= ()40v ϕ=()40a ϕ=3) 推杆位移,速度,加速度线图A.推杆位移线图B.推杆速度线图C.推杆加速度线图4)位移、速度、加速度曲线源代码：x1=0:pi/180:5*pi/6;x2=5*pi/6:pi/180:37*pi/36;x31=37*pi/36:pi/180:47*pi/36;x32=47*pi/36:pi/180:19*pi/12;x4=19*pi/12:pi/180:2*pi;w1=1;s1=168*x1/pi-70*sin(2.4*x1)/pi;v1=168*w1*(1-cos(2.4*x1))/pi;a1=403.2*w1^2*sin(2.4*x1)/pi;s2=140+x2*0;v2=0*x2;a2=0*x2;s31=-907.2*x31.^2/pi^2+1864.8*x31/pi-818.3;v31=-1814.4*w1*x31/pi^2+1864.8*w1/pi;a31=-1814.4*w1^2/pi^2+x31*0;s32=907.2*x32.^2/pi^2-2872.8*x32/pi+2274.3;v32=1814.4*w1*x32/pi^2-2872.8*w1/pi;a32=1814.4*w1^2/pi^2+x32*0;s4=x4*0;v4=0*x4;a4=0*x4;subplot(3,1,1)plot(x1,s1,x2,s2,x31,s31,x32,s32,x4,s4)title('推杆线位移图 ')xlabel('φ(rad)')ylabel('S(mm)')gridsubplot(3,1,2)plot(x1,v1,x2,v2,x31,v31,x32,v32,x4,v4)title('推杆速度线图')xlabel('φ(rad)')ylabel('V(mm/s)')gridsubplot(3,1,3)plot(x1,a1,x2,a2,x31,a31,x32,a32,x4,a4)title('推杆加速度线图')xlabel('φ(rad)')ylabel('a(mm/s^2)')grid3.凸轮机构ds/dψ-s线图，并确定凸轮的基圆半径和偏距.1)凸轮机构的ds/dψ-s线图：使用作图法：推程许用压力角为30度，做淡蓝色直线，其与左端曲线相切并与纵轴夹角为30度；做深绿色直线与纵轴夹角60度且过原点。
回程许用压力角为60度，做深蓝色直线，与右端曲线相切并与纵轴夹角60度；做深绿色直线与纵轴夹角30度且过原点。
凸轮轴心必须位于四直线下方交集，可知凸轮轴心可以选在图中阴影部分。
2)确定凸轮的基圆半径和偏距：由图知:可选取轴心坐标（50，-150）即:基圆半径=158.1mm偏距e=50mm3)凸轮机构ds/dψ-s图源代码：x1=0:pi/180:5*pi/6;x2=5*pi/6:pi/180:37*pi/36;x31=37*pi/36:pi/180:47*pi/36;x32=47*pi/36:pi/180:19*pi/12;x4=19*pi/12:pi/180:2*pi;w1=1;s1=168*x1/pi-70*sin(2.4*x1)/pi;v1=168*w1*(1-cos(2.4*x1))/pi;s31=-907.2*x31.^2/pi^2+1864.8*x31/pi-818.3; v31=-1814.4*w1*x31/pi^2+1864.8*w1/pi; s32=907.2*x32.^2/pi^2-2872.8*x32/pi+2274.3; v32=1814.4*w1*x32/pi^2-2872.8*w1/pi; plot(v1,s1,0,140,v31,s31,v32,s32,0,0); axis([-200 150 -100 200]); title('凸轮机构ds/d ψ-s 线图') xlabeL('ds/d¦×/(mm/s^2)') ylabeL('s/mm') grid4. 滚子半径的确定及凸轮理论轮廓和实际轮廓的绘制。
4) 最小曲率半径数学模型 确定最小曲率半径数学模型如下：)/)(/()/)(/(])/()/[(22222/322ϕϕϕϕϕϕρd x d d dy d y d d dx d dy d dx -+=其中：ϕϕϕϕcos )(sin ])/[(/0s s e d ds d dx ++-= ϕϕϕϕsin )(cos ])/[(/0s s e d ds d dy +--=ϕϕϕϕϕsin ])/[(cos ])/(2[/02222s s d s d e d ds d x d --+-= ϕϕϕϕϕcos ])/[(sin ])/(2[/02222s s d s d e d ds d y d --+--=利用上式可求的最小曲率半斤，而后可确定实际廓线。
5) 凸轮理论轮廓线数学模型ϕϕϕϕsin cos )(cos sin )(00e s s y e s s x -+=++=6) 凸轮实际廓线数学模型22'22')/()/()/()/()/()/(ϕϕϕϕϕϕd dy d dx d dy r y y d dy d dx d dx r x x tt+-=++=(其中r t为确定的滚子半径)7)最小曲率半径计算求滚子半径，需要确定最小曲率半径，以防止凸轮工作轮廓出现尖点或出现相交包络线.8)最小曲率半径计算源代码%最小曲率半径计算函数function pm=qulv(psi,s,f,F,e,r0)e=50;r0=158.1;s0=sqrt(r0^2-e^2);Q1=(s0+s).*cos(psi)+(f-e).*sin(psi);Q2=-(s0+s).*sin(psi)+(f-e).*cos(psi);A0=sqrt(Q1.^2+Q2.^2);A=A0.^3;S1=(2*f-e).*cos(psi)+(F-s0-s).*sin(psi);S2=(F-s0-s).*cos(psi)-(2*f-e).*sin(psi);B=Q1.*S2-Q2.*S1;p=A./B;pm=300;for i=1:length(psi*180/pi)if abs(p(i))<pmpm=abs(p(i));endend%调用主程序w1=1;e=50;r0=158.1;x1=0:pi/180:5*pi/6;x2=5*pi/6:pi/180:37*pi/36;x31=37*pi/36:pi/180:47*pi/36;x32=47*pi/36:pi/180:19*pi/12;x4=19*pi/12:pi/180:2*pi;s1=168*x1/pi-70*sin(2.4*x1)/pi;v1=168*w1*(1-cos(2.4*x1))/pi;a1=403.2*w1^2*sin(2.4*x1)/pi;s2=140+x2*0;v2=0*x2;a2=0*x2;s31=-907.2*x31.^2/pi^2+1864.8*x31/pi-818.3;v31=-1814.4*w1*x31/pi^2+1864.8*w1/pi;a31=-1814.4*w1^2/pi^2+x31*0;s32=907.2*x32.^2/pi^2-2872.8*x32/pi+2274.3;v32=1814.4*w1*x32/pi^2-2872.8*w1/pi;a32=1814.4*w1^2/pi^2+x32*0;s4=x4*0;v4=0*x4;a4=0*x4;rt1=qulv(x1,s1,v1,a1,e,r0);rt2=qulv(x2,s2,v2,a2,e,r0);rt31=qulv(x31,s31,v31,a31,e,r0);rt32=qulv(x32,s32,v32,a32,e,r0);rt4=qulv(x4,s4,v4,a4,e,r0);比较函数四个输出结果rt1,rt2,rt31,rt32,rt4可知最小曲率半径为rt1=158.1mm,为基圆的半径。
9)凸轮轮廓曲线绘制综上,可以取滚子半径为10mm，根据上述数学模型，编程实现凸轮理论轮廓与实际轮廓的绘制。
10)凸轮轮廓曲线源代码%凸轮轮廓绘制函数function picture(psi,s,f,e,r0,rt)s0=sqrt(r0^2-e^2);x1=(s0+s).*cos(psi)-e*sin(psi);y1=(s0+s).*sin(psi)+e*cos(psi);Q1=(s0+s).*cos(psi)+(f-e).*sin(psi);Q2=-(s0+s).*sin(psi)+(f-e).*cos(psi);A0=sqrt(Q1.^2+Q2.^2);x2=x1-rt*Q1./A0;y2=y1+rt*Q2./A0;plot(r0.*cos(psi),r0.*sin(psi),'-',x1,y1,'--',x2,y2),grid on legend('基圆','凸轮理论轮廓','凸轮实际轮廓 ')axis equal%调用主函数x1=0:pi/180:5*pi/6;w1=1;x2=5*pi/6:pi/180:37*pi/36;x31=37*pi/36:pi/180:47*pi/36;x32=47*pi/36:pi/180:19*pi/12;x4=19*pi/12:pi/180:2*pi;s1=168*x1/pi-70*sin(2.4*x1)/pi;v1=168*w1*(1-cos(2.4*x1))/pi;s2=140;v2=0;s31=-907.2*x31.^2/pi^2+1864.8*x31/pi-818.3;v31=-1814.4*w1*x31/pi^2+1864.8*w1/pi;s32=907.2*x32.^2/pi^2-2872.8*x32/pi+2274.3;v32=1814.4*w1*x32/pi^2-2872.8*w1/pi;s4=0;v4=0;e=50;r0=158.1;rt=10;picture(x1,s1,v1,e,r0,rt)hold onpicture(x2,s2,v2,e,r0,rt)hold onpicture(x31,s31,v31,e,r0,rt)hold onpicture(x32,s32,v32,e,r0,rt)hold onpicture(x4,s4,v4,e,r0,rt)axis([-400 300 -300 300])xlabel('x坐标/mm');ylabel('y坐标/mm');title('凸轮理论轮廓与实际轮廓')。