哈工大机械原理大作业 凸轮机构设计 题

合集下载

机械原理大作业凸轮结构设计

机械原理大作业凸轮结构设计

机械原理大作业(二) 作业名称:机械原理设计题目:凸轮机构设计院系: 机电工程学院班级:设计者:学号:指导教师:丁刚陈明设计时间:哈尔滨工业大学机械设计1、设计题目如图所示直动从动件盘形凸轮机构,根据其原始参数设计该凸轮。

表一:凸轮机构原始参数序号升程(mm) 升程运动角(º)升程运动规律升程许用压力角(º)回程运动角(º)回程运动规律回程许用压力角(º)远休止角(º)近休止角(º)12 80 150正弦加速度30 100 正弦加速度60 60 502、凸轮推杆运动规律(1)推杆升程运动方程S=h[φ/Φ0-sin(2πφ/Φ0)]V=hω1/Φ0[1-cos(2πφ/Φ0)]a=2πhω12sin(2πφ/Φ0)/Φ02式中:h=150,Φ0=5π/6,0<=φ<=Φ0,ω1=1(为方便计算)(2)推杆回程运动方程S=h[1-T/Φ1+sin(2πT/Φ1)/2π]V= -hω1/Φ1[1-cos(2πT/Φ1)]a=-2πhω12sin(2πT/Φ1)/Φ12式中:h=150,Φ1=5π/9,7π/6<=φ<=31π/18,T=φ-7π/63、运动线图及凸轮线图运动线图:用Matlab编程所得源程序如下:t=0:pi/500:2*pi;w1=1;h=150;leng=length(t);for m=1:leng;if t(m)<=5*pi/6S(m) = h*(t(m)/(5*pi/6)-sin(2*pi*t(m)/(5*pi/6))/(2*pi));v(m)=h*w1*(1-cos(2*pi*t(m)/(5*pi/6)))/(5*pi/6);a(m)=2*h*w1*w1*sin(2*pi*t(m)/(5*pi/6))/((5*pi/6)*(5*pi/6));% 求退程位移,速度,加速度elseift(m)<=7*pi/6S(m)=h;v(m)=0;a(m)=0;% 求远休止位移,速度,加速度elseif t(m)<=31*pi/18T(m)=t(m)-21*pi/18;S(m)=h*(1-T(m)/(5*pi/9)+sin(2*pi*T(m)/(5*pi/9))/(2*pi));v(m)=-h/(5*pi/9)*(1-cos(2*pi*T(m)/(5*pi/9)));a(m)=-2*pi*h/(5*pi/9)^2*sin(2*pi*T(m)/(5*pi/9));%求回程位移,速度,加速度elseS(m)=0;v(m)=0;a(m)=0;% 求近休止位移,速度,加速度endend推杆位移图推杆速度图推杆加速度图4、确定凸轮基圆半径与偏距在凸轮机构得ds/dφ-s线图里再作斜直线Dt dt与升程得[ds/dφ-s(φ)]曲线相切并使与纵坐标夹角为升程许用压力角[α],则D t d t线得右下方为选择凸轮轴心得许用区。

哈工大、机械原理大作业、凸轮机构设计20题

哈工大、机械原理大作业、凸轮机构设计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 绘制出位移、速度、加速度线图。

哈工大机械原理连杆和凸轮大作业24题

哈工大机械原理连杆和凸轮大作业24题

班级学号机械原理大作业说明书题目 1、连杆机构运动分析2、凸轮机构设计学生姓名1连杆机构运动分析1.设计题目:一、先建立如下坐标系:二、划分杆组如下,进行结构分析:该机构由I级杆组RR(如图1)、II级杆组RPR(如图2、3)和II级杆组RRP(如图4)组成。

(1)(2)(3)(4)三、运动分析数学模型:(1)同一构件上点的运动分析:如右图所示的原动件1,已知杆1的角速度=10/rad s ω,杆长1l =170mm,A y =0,A x =110mm 。

可求得下图中B 点的位置B x 、B y ,速度xB v 、yB v ,加速度xB a 、yB a 。

θcos 1l xB =,θsin 1l yB =θωυsin 1l xB -=,θωυcos 1l yB =,222B2==-cos =-BxB i d x a l x dt ωϕω2222==-sin =-B yB i B d y a l y dtωϕω。

(2)RPRII 级杆组的运动分析:a. 如右图所示是由2个回转副和1个移动副组成的II 级组。

已知两个外运动副C 、B 的位置(B x 、B y 、c x =110mm 、C y =0)、速度(xB υ,yB υ,xC υ=0,yC υ=0)和加速度(0,0,,==yC xC yB xB a a a a )。

可确定下图中D 点的位置、速度和加速度。

确定构件3的角位移1ϕ、角速度1ω、角加速度1α。

1sin 31..ϕϕl x dt dx C B -= 1s i n 131c o s 13.....2ϕϕϕϕl l x dt x d C B --= 1cos 31..ϕϕl y dt dy C B += 1c o s 131s i n 13.....2ϕϕϕϕl l y dty d C B +-= 根据关系:1111d 122..11.αϕϕωϕϕ====dtd dt , 故可得出: D x =)1cos(4βϕ++l x CD y =)1sin(4βϕ++l y Cb. 如右图所示是由2个回转副和1个移动副组成的II级组。

哈工大机械原理大作业凸轮结构设计3

哈工大机械原理大作业凸轮结构设计3

仅供个人参考Harbin Institute of Technology机械原理大作业二课程名称:机械原理设计题目:凸轮结构设计院系:机电工程学院班级:设计者:学号:指导教师:设计时间:哈尔滨工业大学1、设计题目2、凸轮机构推杆升程、回程运动方程,推杆位移、速度、加速度线图。

(1)推杆各行程运动方程(设定角速度为s s rad /2/1πω==)①从动件推程运动方程(650πϕ≤≤) ⎥⎦⎤⎢⎣⎡⎪⎪⎭⎫ ⎝⎛-=ϕφπφω0012cos 1h v ; 代入数据,可得:②从动件远休程运动方程(πϕπ≤≤65) ③从动件回程运动方程(914πϕπ≤≤) 代入数据,可得:④从动件进休程运动方程(πϕπ2914≤≤) (2)、推杆位移、速度、加速度线图①推杆的位移线图如下②推杆的速度线图如下③推杆的加速度线图如下 3、凸轮机构的dss d ϕ-线图,并依次确定凸轮的基圆半径和偏距凸轮机构的s d ds -ϕ线图如下图所示 由图中范围选定点(-10,-50)为凸轮转轴O 点,则mm r 99.501050220=+=取基圆半径为r 0 =51mm ,偏距e = 10mm 。

4、滚子半径的确定及凸轮理论廓线和实际廓线的绘制由程序计算得凸轮理论轮廓线最小曲率半径mm r 51min = .由滚子半径选择范围∆-<min ρr r ,mm 5~3=∆得到滚子半径mm r r 46≤.又因为凸轮整体尺寸较小,此范围明显过大,故适当减小滚子半径,这里取半径为mm r r 12= .得到图线为:附录1.求位移、速度、加速度的程序(matlab )function f = tulunh=50;x1=150;t1=30;x2=100;t2=80;w=2*pi;x1=x1*pi/180;x2=x2*pi/180;t1=t1*pi/180;t2=t2*pi/180;%升程x=0:0.001:x1;s = h*(x/x1-sin(2*pi*x/x1)/(2*pi));v = h*w*(1-cos(2*pi*x/x1))/x1;a = 2*pi*h*w*w*sin(2*pi*x/x1)/(x1*x1);subplot(3,1,1),plot(x,s),hold onsubplot(3,1,2),plot(x,v),hold onsubplot(3,1,3),plot(x,a),hold on%远休x = x1:0.001:x1+t1;s = h;v=0;a=0;subplot(3,1,1),plot(x,s),hold onsubplot(3,1,2),plot(x,v),hold onsubplot(3,1,3),plot(x,a),hold on%回程x= x1+t1:0.001:x1+t1+x2;s = h*(1+cos(pi*(x-(x1+t1))/x2))/2;v = -pi*h*w*sin(pi*(x-(x1+t1))/x2)/(2*x2);a = -pi*pi*h*w*w*cos(pi*(x-(x1+t1))/x2)/(2*x2*x2);subplot(3,1,1),plot(x,s),hold onsubplot(3,1,2),plot(x,v),hold onsubplot(3,1,3),plot(x,a),hold on%近休x=x1+t1+x2:0.001:x1+x2+t1+t2;s = 0;v = 0;a = 0;subplot(3,1,1),plot(x,s),xlabel('φ/rad'),ylabel('S/mm'),title('位移-转角图线'),hold onsubplot(3,1,2),plot(x,v),xlabel('φ/rad'),ylabel('v/(mm/s)'),title('速度-转角图线'),hold onsubplot(3,1,3),plot(x,a),xlabel(φ/rad'),ylabel('a/(mm/s^2)'),title('加速度-转角图线'),hold on2.绘制凸轮机构s d ds -ϕ线图 function f= jiyuan;x1=150;t1=30;x2=100;t2=80;h=50;x1=x1*pi/180;x2=x2*pi/180;t1=t1*pi/180;t2=t2*pi/180;x= 0:0.001:150*pi/180;%升程 v/ws = h*(x/x1-sin(2*pi*x/x1)/(2*pi));k =-h*(1-cos(2*pi*x/x1))/x1;plot(k,s,'r'),hold on ;x=180*pi/180:0.001:280*pi/180;%回程 v/ws = h*(1+cos(pi*(x-(x1+t1))/x2))/2;k = pi*h*sin(pi*(x-(x1+t1))/x2)/(2*x2);plot(k,s,'g'),hold on ;%回程切线for i=-11:1:-11;f=@(k)k*tan(pi/6)+i;k =-40:0.1:50;s=f(k);plot(k,s),hold on ;end%升程切线for i=-45:0.2:-45;f=@(k)-k*tan(60*pi/180)+i;k =-40:0.1:50;s=f(k);plot(k,s),hold on ;endgrid onf=@(k)k*tan(50*pi/180);k=-50:0.1:0;s=f(k);plot(k,s),hold onxlabel('ds/d φ');ylabel('s(φ)');title('类速度-位移图线 ');plot(-10,-50,’o ’);3.绘制凸轮轮廓曲线function f= lunkuo;h=50;x1=150;t1=30;x2=100;t2=80;x1=x1.*pi./180;x2=x2.*pi./180;t1=t1.*pi./180;t2=t2.*pi./180;s0=51;e=10;rr=12;%升程x=0:pi/200:150.*pi/180;s = h.*(x./x1-sin(2.*pi.*x./x1)./(2.*pi));X1=(s0+s).*cos(x)-e.*sin(x);Y1=(s0+s).*sin(x)+e.*cos(x);%实际轮廓X11=X1-(rr.*(cos(x).*(s + s0) - e.*sin(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);Y11=Y1-(rr.*(sin(x).*(s + s0) + e.*cos(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);plot(X1,Y1,'r',X11,Y11,'r'),hold on;%远休x=150.*pi/180:pi/180:180.*pi/180;s=50;X2=(s0+s).*cos(x)-e.*sin(x);Y2=(s0+s).*sin(x)+e.*cos(x);X22=X2-(rr.*(cos(x).*(s + s0) - e.*sin(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);Y22=Y2-(rr.*(sin(x).*(s + s0) + e.*cos(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);plot(X2,Y2,'g',X22,Y22,'g'),hold on;%回程x=180.*pi/180:pi/180:280.*pi/180;s = h.*(1+cos(pi.*(x-(x1+t1))./x2))./2;X3=(s0+s).*cos(x)-e.*sin(x);Y3=(s0+s).*sin(x)+e.*cos(x);X33=X3-(rr.*(cos(x).*(s + s0) - e.*sin(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);Y33=Y3-(rr.*(sin(x).*(s + s0) + e.*cos(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);plot(X3,Y3,'k',X33,Y33,'k'),hold on;%近休x=280*pi/180:pi/180:2*pi;s=0;X4=(s0+s).*cos(x)-e.*sin(x);Y4=(s0+s).*sin(x)+e.*cos(x);X44=X4-(rr.*(cos(x).*(s + s0) - e.*sin(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);Y44=Y4-(rr.*(sin(x).*(s + s0) + e.*cos(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);plot(X4,Y4,'b',X44,Y44,'b'),hold on;x=0:pi/200:2*pi;X4=(s0+s).*cos(x)-e.*sin(x); Y4=(s0+s).*sin(x)+e.*cos(x); plot(X4,Y4,'b');text(-40,90,'理论轮廓线');text(-40,72,'实际轮廓线');text(-5,55,'基圆');grid on;axis equal仅供个人用于学习、研究;不得用于商业用途。

哈工大机械原理大作业——凸轮——5号

哈工大机械原理大作业——凸轮——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号

哈工大机械原理大作业——凸轮——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ϕϕ⨯=⨯通过编程求其角度。

哈工大机械原理大作业-凸轮结构设计 -05

哈工大机械原理大作业-凸轮结构设计 -05

机械原理大作业二课程名称:机械原理设计题目:凸轮机构设计院系:机电工程学院班级:完成者:学号:指导教师:设计时间:设计说明书1 设计题目:如图所示直动从动件盘形凸轮机构,其原始参数见表2-1。

从表2-1中选择一组凸轮机构的原始参数,据此设计该凸轮机构。

注: 5题参数有误,我将远休止角改为45 ,将近休止角改为105 .2、推杆升程、回程运动方程及位移、速度、加速度线图:2.1凸轮运动分析: 推程运动方程; 01cos 2h s πϕ⎡⎤⎛⎫=-⎢⎥ ⎪Φ⎝⎭⎣⎦100sin 2h v πωπϕ⎛⎫=⎪ΦΦ⎝⎭序号 升程 (mm ) 升程运 动角( ) 升程运 动规律 升程 许用压力角 () 回程运动角 () 回程运 动规律 回程许用 压力角 ( ) 远休 止角 ( ) 近休 止角 ( ) 5(改)60120余弦加速度3590等速6545105ω221200cos 2h a πωπϕ⎛⎫=⎪ΦΦ⎝⎭回程运动方程: ()0'01s s h ϕ-Φ+Φ⎡⎤=-⎢⎥Φ⎣⎦1'0h v ω=-Φ 0a =2.2求位移、速度、加速度线图MATLAB 源程序: pi=3.1415926; c=pi/180; h=60; f0=120; fs=45; f01=90; fs1=105;f=0:1:360; %升程阶段 for n=0:f0s(n+1)=h/2*(1-cos(pi/f0*f(n+1)));v(n+1)=pi*h/(2*f0*c)*sin(pi/f0*f(n+1));a(n+1)=pi^2*h/(2*f0^2*c^2)*cos(pi/f0*f(n+1)); end%远休止阶段 for n=f0:f0+fs s(n+1)=60; v(n+1)=0; a(n+1)=0; endfor n=f0+fs:f0+fs+f01 %回程阶段 s(n+1)=h*(1-(f(n+1)-(f0+fs))/f01); v(n+1)=-h/(f01*c); a(n+1)=0; end;for n=f0+fs+f01:360; %近休止阶段 s(n+1)=0; v(n+1)=0; a(n+1)=0; endfigure(1);plot(f,s,'r');xlabel('\phi/\circ');ylabel('s/mm');grid on;title('推杆位移线图')figure(2);plot(f,v,'r');xlabel('\phi/\circ');ylabel('v/\ommiga');grid on;title('推杆速度线图')figure(3);plot(f,a,'r');xlabel('\phi/\circ');ylabel('a/\ommiga');grid on;title('推杆加速度线图')2.3位移、速度、加速度线图:3 凸轮机构的dss d ϕ-线图,确定基圆半径和偏心距: 3.1理论分析:机构压力角α应按下式计算: 0-tan =+ds d es sϕα作出-ds d s ϕ曲线,因凸轮机构压力角0-tan =+ds d e s sϕα,式中,2200=-s r e ,右侧为升程,作与s 轴夹角等于升程许用压力角35的切界线D d t t ,则在直线上或其右下方取凸轮轴心时,可使[]0-=arctan+ds d es sϕαα≤,同理左侧回程,作与s 轴夹角等于回程许用压力角65的切界线D 'd 't t ,则在直线上或其左下方取凸轮轴心时,使[]0-=arctan+ds d es sϕαα≤。

哈工大—机械原理凸轮大作业

哈工大—机械原理凸轮大作业

一、题目要求及机构运动简图如图1所示直动从动件盘形凸轮机构。

其原始参数见表1。

图一凸轮运动简图表一凸轮原始参数二、计算流程框图凸轮机构分析建立数学模型位移方程速度方程加速度方程速度线图位移线图加速线图ds/dΨ-s曲线升程压力角回程压力角确定轴向及基圆半径压力角图确定滚子半径实际轮廓理论轮廓轮廓图结束三、建立数学模型1.从动件运动规律方程首先,由于设计凸轮轮廓与凸轮角速度无关,所以不妨设凸轮运动角速度为w = 1rad/s。

(1)推程运动规律(0 < φ < 90°)s=φ)v=φφ)a=φφ)式中:h=65mm,Φ0=π/2(2)远休程运动规律(90°< φ < 190°)s = 65mmv = 0a = 0(3)回程运动规律(190°< φ < 240°)(190°< φ < 196.25°)(196.25°< φ < 233.75°)(233.75°< φ < 240°) 回程运动中的速度和加速度为位移对时间t的倒数:(4)近休程运动规律(240°< φ < 360°)s = 0v = 0a = 02.从动件位移、速度、加速度线图(1)位移线图(2)速度线图(3)加速度线图(4)位移、速度、加速度线图MATLAB源程序%% 已知条件h = 65; %mmphi_0 = 90./180*pi; %radalpha_up_al = 35./180*pi; %升程许用压力角phi_00 = 50./180*pi;alpha_down_al = 70./180*pi; %回程许用压力角phi_s = 100./180*pi;phi_ss = 120./180*pi;w = 1;%% 绘制从动件位移、速度、加速度线图% 推程阶段t_up = 0 : 0.5 : 90;t_up1 = t_up./180*pi;syms t_up1phi_up s_up v_up a_upphi_up = w.*t_up1;s_up = h./2.*(1 - cos(pi.*phi_up./phi_0));v_up = diff(s_up,t_up1);a_up = diff(v_up,t_up1);s_up1 = double(subs(s_up,t_up./180*pi));v_up1 = double(subs(v_up,t_up./180*pi));a_up1 = double(subs(a_up,t_up./180*pi));% 远休程t_s = 90 : 0.5 : (90+100);t_s1 = t_up./180*pi;s_s(1:201) = h;v_s(1:201) = 0;a_s(1:201) = 0;% 回程阶段1t_down1 = (90+100) : 0.5 : (90+100+50/8);t_down11 = t_down1./180*pi;syms t_down11phi_down1s_down1v_down1a_down1phi_down1 = w.*t_down11;s_down1 = h - h./(4+pi).*(pi.*(phi_down1 - phi_0 - phi_s)./phi_00 - ...sin(4.*pi.*(phi_down1 - phi_0 - phi_s)./phi_00)./4);v_down1 = diff(s_down1,t_down11);a_down1 = diff(v_down1,t_down11);s_down11 = double(subs(s_down1,t_down1./180*pi));v_down11 = double(subs(v_down1,t_down1./180*pi));a_down11 = double(subs(a_down1,t_down1./180*pi));% 回程阶段2t_down2 = (90+100+50/8) : 0.5 : (90+100+7*50/8);t_down22 = t_down2./180*pi;syms t_down22phi_down2s_down2v_down2a_down2phi_down2 = w.*t_down22;s_down2 = h - h./(4+pi).*(2+pi.*(phi_down2 - phi_0 - phi_s)./phi_00 - 9.*sin(pi./3 + 4.*pi.*(phi_down2 - phi_0 - phi_s)./(3.*phi_00))./4);v_down2 = diff(s_down2,t_down22);a_down2 = diff(v_down2,t_down22);s_down22 = double(subs(s_down2,t_down2./180*pi));v_down22 = double(subs(v_down2,t_down2./180*pi));a_down22 = double(subs(a_down2,t_down2./180*pi));% 回程阶段3t_down3 = (90+100+7*50/8) : 0.5 : (90+100+50);t_down33 = t_down3./180*pi;syms t_down33phi_down3s_down3v_down3a_down3phi_down3 = w.*t_down33;s_down3 = h - h./(4+pi).*(4+pi.*(phi_down3 - phi_0 - phi_s)./phi_00 - …sin(4.*pi.*(phi_down3 - phi_0 - phi_s)./phi_00)./4);v_down3 = diff(s_down3,t_down33);a_down3 = diff(v_down3,t_down33);s_down33 = double(subs(s_down3,t_down3./180*pi));v_down33 = double(subs(v_down3,t_down3./180*pi));a_down33 = double(subs(a_down3,t_down3./180*pi));% 近休程t_ss = (90+100+50) : 0.5 : 360;s_ss(1:241) = 0;v_ss(1:241) = 0;a_ss(1:241) = 0;% 绘图位移t = [t_up t_s t_down1 t_down2 t_down3 t_ss];phi = w .* t ./ 180 .*pi;s = [s_up1 s_s s_down11 s_down22 s_down33 s_ss];v = [v_up1 v_s v_down11 v_down22 v_down33 v_ss];a = [a_up1 a_s a_down11 a_down22 a_down33 a_ss];figure('Name','从动件位移-时间线图');plot(t,s,'k','linewidth',1.0);grid on;title('从动件位移-时间线图');xlabel('转角\phi / 度');ylabel('位移h/mm');% 绘图速度figure('Name','从动件速度-时间线图');plot(t,v,'k','linewidth',1.0);grid on;title('从动件速度-时间线图');xlabel('转角\phi / 度');ylabel('速度v/mm*s^{-1}');% 绘图加速度figure('Name','从动件加速度-时间线图');plot(t,a,'k','linewidth',1.0);grid on;title('从动件加速度-时间线图');xlabel('转角\phi / 度');ylabel('加速度a/mm*s^{-2}');3.绘制ds/dΦ线图并确定基圆半径和偏距(1)绘制ds/dΦ线图及源程序①MATLAB源程序:%% 绘制ds/dphi-s线图,确定基圆半径和偏距ds_dphi = v ./ w;figure('Name','凸轮ds/dphi - s线图');plot(ds_dphi,s,'k','linewidth',1.5);hold on;axis([-150 150 -70 70]);grid on;title('凸轮ds/dphi - s线图');xlabel('ds/dphi / (mm*s^{-2})');ylabel('s/mm');% 三条临界线x = linspace(-150,150,301);k_up = tan(pi/2 - alpha_up_al);y_up = k_up.*x - 66;plot(x,y_up,'linewidth',1.5);k_down = - tan(pi/2 - alpha_down_al);y_down = k_down.*x - 24.7;plot(x,y_down,'linewidth',1.5);x0 = linspace(0,150,151);k0 = - tan(alpha_up_al);y0 = k0.*x0;plot(x0,y0,'--');% 由图像选取凸轮基圆半径为r0 = sqrt(23^2 + 34^2) = 41 mm,偏距e = 23mm plot(23,-34,'or');r0 = 41;e = 23;plot(linspace(0,23,10),linspace(0,-34,10),'r',linspace(0,23,10),linspace(-34,-34,10),'r',li nspace(23,23,10),linspace(0,-34,10),'r','linewidth',1.0);(2)确定基圆半径和偏距在凸轮机构的ds/d φ-s 线图里再作斜直线Dt-dt 与升程的[ds/d φ-s]曲线相切并使与纵坐标夹角为升程许用压力角[α],则Dt-dt 线的右下方为选择凸轮轴心的许用区。

哈工大机械原理大作业2凸轮机构设计

哈工大机械原理大作业2凸轮机构设计

机械原理大作业(二)作业名称:凸轮机构设计设计题目:23题院系:班级:设计者:学号:指导教师:设计时刻:哈尔滨工业大学机械设计1.运动分析题目:设计直动从动件盘形凸轮机构,其原始参数见下表2.确信凸轮机构推杆升程、回程运动方程(设定角速度为ω=10 rad/s)升程:0°< Φ < 120°由公式可得:s=60-60*cos(3*Φ/2);v=90*ω*sin(3*Φ/2);a=135*ω2 *cos(3*Φ/2);远停止:120°< Φ < 200°由公式可得:s=120;v=0;a=0;回程:200°< Φ < 290°由公式可得:s=h[1-(10T23-15T24+6T25)]v=(-30hω1/Φ0')T22(1–2T2+T22)a=(-60hω12/Φ0'2)T2(1–3T2+2T22)式中:T2=(Φ-Φ0-Φs)/ Φ0'近停止: 290°< Φ < 360°由公式可得:s=0;v=0;a=0;3.绘制推杆位移、速度、加速度线图(设ω=10rad/s)1) 推拉位移曲线代码:%推杆位移曲线;x=0:(pi/1000):(2*pi/3);s1=60-60*cos(1.5*x);y=(2*pi/3):(pi/1000):(10*pi/9);s2=120;z=(10*pi/9):(pi/1000):(29*pi/18);T2=(z-10*pi/9)*2/pi;s3=120*(1-(10*T2.^3-15*T2.^4+6*T2.^5));m=(29*pi/18):(pi/1000):(2*pi);s4=0;plot(x,s1,'b',y,s2,'b',z,s3,'b',m,s4,'b'); xlabel('角度(rad)');ylabel('行程(mm)');title('推杆位移曲线');grid;2)推杆速度曲线代码:%推杆速度曲线;w=10;x=0:(pi/1000):(2*pi/3);v1=90*w*sin(1.5*x);y=(2*pi/3):(pi/1000):(10*pi/9);v2=0;z=(10*pi/9):(pi/1000):(29*pi/18);T2=(z-10*pi/9)*2/pi;v3=(-30*120*w/(pi/2))*T2.^2.*(1-2*T2.^2+T2.^2); % v3=-120*w*sin(2*z-20*pi/9);m=(29*pi/18):(pi/1000):(2*pi);v4=0;plot(x,v1,'r',y,v2,'r',z,v3,'r',m,v4,'r'); xlabel('角度(rad)');ylabel('速度(mm/s)');title('推杆速度曲线(w=10rad/s)');grid;3)凸轮推杆加速度曲线代码:%凸轮推杆加速度曲线;w=10;x=0:(pi/1000):(2*pi/3);a1=135*w^2*cos(3*x/2);y=(2*pi/3):(pi/1000):(10*pi/9);a2=0;z=(10*pi/9):(pi/1000):(29*pi/18);T2=(z-10*pi/9)*2/pi;a3=(-60*120*w^2/(pi/2)^2)*T2.*(1-3*T2.^2+2*T2.^2); m=(29*pi/18):(pi/1000):(2*pi);a4=0;plot(x,a1,'m',y,a2,'m',z,a3,'m',m,a4,'m');xlabel('角度(rad)');ylabel('加速度(mm/s^2)');title('凸轮推杆加速度曲线(w=10rad/s)');grid;4)绘制凸轮机构的dd/dd−d线图,并依次确信凸轮的基圆半径和偏距代码:%dd/dd−d线图,确信e,s0;x=0:(pi/1000):(2*pi/3);s1=60-60*cos(1.5*x);ns1=90*sin(1.5*x);y=(2*pi/3):(pi/1000):(10*pi/9);s2=120;ns2=0;z=(10*pi/9):(pi/1000):(29*pi/18);T2=(z-10*pi/9)*2/pi;s3n=120*(1-(10*T2.^3-15*T2.^4+6*T2.^5));ns3=-120*10*3*T2.^2+120*15*4*T2.^3-120*6*5*T2.^4 ;m=(29*pi/18):(pi/1000):(2*pi);s4=0;ns4=0;x1=0:pi/36000:pi/2;s1n=60-60*cos(1.5*x1);v1=90*sin(1.5*x1);m1=diff(s1n);%求切线1n1=diff(v1);z=m1./n1;for i=1:length(z);if abs(z(i)+tan(-55*pi/180))<0.001;breakendendb11=s1n(i)-z(i)*v1(i);x1=-300:200;y01=z(i)*x1+b11;%切线1k1=z(i);plot(x1,y01)x3=10*pi/9:pi/36000:14*pi/9;%求切线2s3n=120*(1-(10*T2.^3-15*T2.^4+6*T2.^5));v3=-120*10*3*T2.^2+120*15*4*T2.^3-120*6*5*T2.^4 ;m3=diff(s3n);n3=diff(v3);p=m3./n3;for o=1:length(p);if abs(p(o)-tan(-25*pi/180))<0.01;breakendendo;b33=s3n(o)-p(o)*v3(o);x3=-300:700;y03=p(o)*x3+b33;%切线2plot(x3,y03);sym uv[u,v]=solve('u= 1.4281*v-81.7665','u=-0.4663*v-59.6715');%v=11.66332347972972972972972972973 x%u=-65.110107738597972972972972972973 yplot(ns1,s1,'m',ns2,s2,'b',ns3,s3n,'b',ns4,s4,'b',x1,y01,'g',x3,y03,'g',v,u,'*'); xlabel('ds/d¦µ');ylabel('S');axis([-300,200,-300,300]);title('s0,e 的确信');grid;确信凸轮基圆半径与偏距:偏距e=90mm,d020mm;基圆半径为d0=150mm。

哈尔滨工业大学机械原理大作业凸轮设计

哈尔滨工业大学机械原理大作业凸轮设计

大作业1连杆机构运动分析1.题目(8)如图所示机构,已知机构各构件的尺寸为==100AC CE l l mm ,==200BC CD l l mm ,90BCD ∠=︒,构件1的角速度为1=10/w rad s ,试求构件5的角位移、角速度和角加速度,并对计算结果进行分析。

1.1机构的运动分析AB 为原动件,AB 转动通过转动导杆机构带动杆BCD 转动,BCD 转动通过转动导杆机构带动杆DE 摆动。

1.2 机构的结构分析杆组可以划分为一个RR I级杆组(杆1)、RRPII级杆组(滑块2,杆3)、RPRII 级杆组(滑块4,杆5)(1)RRI级杆组1:(2)RRPII级杆组2,3:(3)RPRII级杆组4,5:2.分析过程 2.1 建立坐标系建立以点E 为原点的固定平面直角坐标系x-E-y ,如图所示。

2.2 建立数学模型(1)构件1、2、3的分析原动件杆1的转角:1θ=0--360。

原动件杆1的角速度:1ω=.1θ=10/rad s 原动件杆1的角加速度:..1αθ==0运动副A 的坐标:0200A A x y mm =⎫⎬=⎭运动副A 的速度及加速度都为零。

构件1为BC (RRP Ⅱ级杆组)上滑块B 的导路 滑块B 的位置为:132cos cos B A C x x s x l θθ=+=+132sin sin B A C y y s x l θθ=+=+消去s,得:0212arcsinA l θθ=+式中:011()sin ()cos C A C A A x x y y θθ=---构件3的角速度i ω和滑块B 沿导路的移动速度D υ:.211213(Q sin Q cos )/Q ωϕθθ==-+ 1322323(Q cos Q sin )/Q D s l l υθθ⋅==-+式中:..11111211321212Q sin ;Q cos ;Q sin sin cos sin l l l θθθθθθθθ=-==+构件3的角加速度和滑块B 沿导路移动的加速度:..241513(Q sin Q cos )/Q αθθθ==-+..4325323(Q cos Q sin )/Q B s l l υθθ==-+式中:122......21142211111Q cos sin cos 2sin l l l s θθθθθθθθ=---- 122......21152211111Q sin cos sin 2cos l l l s θθθθθθθθ=+-+(2)构件3,4,5的分析构件3,4,5,由1个Ⅰ级基本杆组和一个RRP Ⅱ级杆组组成,与构件1,2,3结构相同,只运动分析过程与其相反。

哈工大机械原理大作业凸轮

哈工大机械原理大作业凸轮

1 / 17机械原理大作业课程名称:机械原理 设计题目:直动从动件盘形凸轮机构 院 系: 机电学院 班 级: 完 成 者: 学 号: 指导教师: 设计时间:哈尔滨工业大学题目:如图所示直动从动件盘形凸轮机构,其原始参数见表,据此设计该凸轮机构。

凸轮运动分为五个阶段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'204s h v ωϕ=-Φ+Φ+Φ-Φ21'204h v ω=-Φ5.近休止 240~360其中,0Φ 推程运动角 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/fio1).^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 onplot(ds4,s4);fi5=fio1+fis+fio2:0.001:2*pi; s5=0+0*fi5; ds5=0*fi5; hold onplot(ds5,s5);title('ds/d φ-s');xlabel('ds/d φ(mm/rad)'),ylabel('s(mm)'); grid on hold off5.凸轮轴心位置的确定凸轮压力角的正切值s s ed ds +-=0/tan ϕα,右侧为升程,作与s 轴夹6π角等于升程许用压力角的切界线t t dD ,则在直线上或其下方取凸轮轴心时,可使[]αα≤,同理右侧回程,作与s 轴夹角等于回程许用压力角3π的切界线''t t d D ,则在直线上或其下方取凸轮轴心时,可使[]αα≤。

哈工大机械原理大作业二凸轮机构

哈工大机械原理大作业二凸轮机构
#define PI 3.14159265
void main()
{
float fi=0,fi0=90,h=100,S=0,v=0,a=0,fis=110,fi01=80,fis1=80,T2=0,w1=10,A=0,B=0,v0=0,v01=0,v02=0,f11=30,f31=70,A11=0,A12=0,B11=0,B12=0,A21=0,A22=0,B21=0,C11=0,C21=0,W=0;
{
W=fi*2*PI/360.0;
if(fi<=90)
{
S=h*(W/fi0-(1/(2*PI))*sin((2*PI/fi0)*W));
v=(h*w1/fi0)*(1-cos((2*PI/fi0)*W));
a=(2*PI*h*w1*w1/(fi0*fi0))*sin(2*PI*W/fi0);
}
elseif(fi>=90&&fi<200)
{
A=90;
A=A*2*PI/360.0;
S=h*(A/fi0-(1/(2*PI))*sin((2*PI/fi0)*A));
v=(h*w1/fi0)*(1-cos((2*PI/fi0)*A));
a=(2*PI*h*w1*w1/(fi0*fi0))*sin(2*PI*A/fi0);
}
elseif(fi>=200&&fi<=280)
}
printf("%f\n",a);
}
}
此程序计算推杆位移,速度,加速度线图。
程序二:
#include "stdio.h"
#include "math.h"

(完整word版)哈工大机械原理大作业凸轮DOC

(完整word版)哈工大机械原理大作业凸轮DOC

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机械原理大作业二课程名称:机械原理设计题目: 凸轮机构设计院系:班级:设计者:学号:指导教师:哈尔滨工业大学一、设计题目如右图所示直动从动件盘形凸轮机构,选择一组凸轮机构的原始参数,据此设计该凸轮机构。

凸轮机构原始参数序号升程(mm)升程运动角升程运动规律升程许用压力角27130150正弦加速度30°回程运动角回程运动规律回程许用压力角远休止角近休止角100°余弦加速度60°30°80°二. 凸轮推杆升程、回程运动方程及推杆位移、速度、加速度线图凸轮推杆升程运动方程:)]512sin(2156[130s ϕππϕ-= )512sin(4.374)]512cos(1[156v 211ϕπϕπωω=-=a% t 表示转角,s 表示位移t=0:0.01:5*pi/6;%升程阶段s= [(6*t)/(5*pi )- 1/(2*pi )*sin(12*t/5)]*130; hold on plot(t ,s ); t= 5*pi/6:0。

01:pi; %远休止阶段s=130; hold on plot(t,s );t=pi :0.01:14*pi/9;%回程阶段s=65*[1+cos(9*(t-pi )/5)]; hold on plot(t ,s );t=14*pi/9:0.01:2*pi ;s=0;hold onplot(t,s);grid onhold off%t表示转角,令ω1=1t=0:0。

01:5*pi/6;%升程阶段v=156*1*[1-cos(12*t/5)]/pi hold onplot(t,v);t= 5*pi/6:0。

01:pi;v=0hold onplot(t,v);t=pi:0.01:14*pi/9;%回程阶段v=—117*1*sin(9*(t—pi)/5) hold onplot(t,v);t=14*pi/9:0。

哈工大机械原理大作业——凸轮——14号

哈工大机械原理大作业——凸轮——14号

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机械原理大作业二课程名称:机械原理设计题目:凸轮机构设计院系:机电工程学院班级:设计者:学号:指导教师:设计时间:一.设计题目设计直动从动件盘形凸轮机构,凸轮机构原始参数 序号 升程(mm) 升程运动角(º) 升程运 动规律升程许用压力角(º) 回程运动角(º)回程运 动规律回程许用压力角(º)远休止角 (º)近休止角 (º) 22 120 90余弦加速度 3080正弦加速度 70 60130二、凸轮推杆升程、回程运动方程及推杆位移,速度、加速度图1.运动规律A .从动件升程运动方程: 升程阶段采用余弦加速度运动规律当02πϕ≤≤时s 60[1cos2]ϕ=-1v 120sin2ωϕ=21a 240cos2ωϕ=B.从动件远休止运动方程:当526ππϕ≤≤时S=h=120v=0a=0C.从动件回程运动方程:回程阶段采用正弦加速度运动规律当523618ππϕ≤≤时91120 2.875sin(4.5 3.75)42sϕϕπππ⎡⎤=⨯-+-⎢⎥⎣⎦[]12701cos(4.5 3.75)vωϕππ=---211215sin(4.5 3.75)aωϕππ=--D.近休止运动阶段方程:当29π/18≤Φ≤2π时s=0v=0a=0推杆位移、速度、加速度线图采用matlab编程,为作图方便,取ω1=1rad/s:其源程序及图像如下:1.位移:编程源代码:t1=0:0.000:pi/2;s1=60.*(1-cos(2.*t1));t2=pi/2:0.0001:5*pi/6;s2=120;t3=5*pi/6:0.0001:23.*pi/18;s3=120.*(2.875-9.*t3./4./pi+1./(2.*pi).*sin(4.5.*t3-3.75.*pi));t4=23.*pi/18:0.0001:2.*pi;s4=0;plot(t1,s1,t2,s2,t3,s3,t4,s4)线图2.速度编程源代码:t1=0:0.000:pi/2;v1=120.*sin(2.*t1);t2=pi/2:0.0001:5*pi/6;v2=0;t3=5*pi/6:0.0001:23.*pi/18;v3=-270./pi.*(1-cos(4.5.*t3-3.75*pi)); t4=23.*pi/18:0.0001:2.*pi;v4=0;plot(t1,v1,t2,v2,t3,v3,t4,v4);线图:3.加速度编程源代码:4.t1=0:0.000:pi/2;5.a1=240.*sin(2.*t1);6.t2=pi/2:0.0001:5*pi/6;7.a2=0;8.t3=5*pi/6:0.0001:23.*pi/18;9.a3=-81.*15./pi.*sin(4.5.*t3-3.75.*pi);10.t4=23.*pi/18:0.0001:2.*pi;11.a4=0;12.plot(t1,a1,t2,a2,t3,a3,t4,a4);线图:三、凸轮机构ds/dψ和s的关系图与基圆半径的确定和偏距的确定:源代码:i=1;for t=0:0.01:2*pi;i=i+1;if t<pi/2ds=120.*sin(2.*t);s=60.*(1-cos(2.*t));elseif t<5*pi/6ds=0;s=0;elseif t<23*pi/18ds=-270./pi.*(1-cos(4.5.*t-3.75*pi));s=120.*(2.875-9.*t./4./pi+1./(2.*pi).*sin(4.5.*t-3.75.*pi));else t<2*pi; ds=0; s=0; end Ds(i)=ds; S(i)=s; end plot(Ds,S)凸轮机构的s d ds-φ线图凸轮的基圆半径和偏距确定 以s d ds-φ图为基础,可分别作出三条限制线(推程许用压力角的切界限D t d t ,回程许用压力角的限制线D t 'd t ',起始点压力角许用线B 0d ''),以这三条线可确定最小基圆半径及所对应的偏距e,在其下方选择一合适点,即可满足压力角的限制条件。

哈工大机械原理大作业第25题凸轮设计

哈工大机械原理大作业第25题凸轮设计

机械原理大作业二课程名称:机械原理设计题目:凸轮机构设计院系:机电工程学院指导教师:林琳福利设计时间:2014.06工业大学凸轮机构设计1.设计题目(1)凸轮机构运动简图:(2)凸轮机构的原始参数序号升程升程运动角升程运动规律升程许用压力角回程运动角回程运动规律回程许用压力角远休止角近休止角25 130mm 90°等加等减速40°80°等速70°60°130°2.确定凸轮推杆升程、回程运动方程并绘制推杆位移,速度,加速度线图:(1) 凸轮推杆升程,回程运动方程如下:A.推杆升程方程B.推杆回程方程(2)推杆位移,速度,加速度线图如下:A.推杆位移线图(使用matlab画图,程序详见附录1)B.推杆速度线图(使用matlab画图,程序详见附录2)C.推杆加速度线图(使用matlab画图,程序详见附录3)3.凸轮机构的-s线图,并依次确定凸轮的基圆半径和偏距.(1) 凸轮机构的-s线图:(使用matlab画图,程序详见附录4)(2)确定凸轮的基圆半径和偏距:以ds/df-s(f)图为基础,可分别作出三条限制线(推程许用压力角的切界限Dt d t ,回程许用压力角的限制线Dt 'dt',起始点压力角许用线Bd''),以这三条线可确定最小基圆半径及所对应的偏距e,在其下方选择一合适点,即可满足压力角的限制条件。

得图如下:(使用matlab画图,程序详见附录6)得最小基圆对应的坐标位置大约为(55.28,-65.88)经计算取偏距e=55mm,r0=90mm.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 --+--=利用上式可求的最小曲率半径后可确定实际廓线。

哈工大机械原理大作业2-凸轮

哈工大机械原理大作业2-凸轮

Harbin Institute of Technology机械原理大作业二课程名称:机械原理设计题目:凸轮结构设计院系:能源学院指导教师:唐德威赵永强设计时间:2013年6月27日哈尔滨工业大学一、设计题目1、凸轮机构运动简图:2、凸轮机构的原始参数序号升程升程运动角升程运动规律升程许用压力角回程运动角回程运动规律回程许用压力角远休止角近休止角21 110 150°3-4-5多项式40°100°3-4-5多项式60°45°65°二、凸轮推杆升程、回程运动方程及推杆位移,速度,加速度线图1、推杆升程,回程运动方程如下:A.推杆升程方程:设11/ rad sω=,由3-4-5多项式可知:当56ϕπ≤≤时,有:3451112211112211112(10156)30(12)60(132)s h T T Th Tv T Tha T T Tωφωφ=-+=-+=-+式中10/ϕφT=H=110 ,056φπ=B.推杆回程方程:当13π/12 ≤φ≤59π/36 时,有:2345222221222'0221222'0[1(10156)]30(12)60(132)s h T T T h v T T T h a T T T ωφωφ=--+=--+=--+ 式中 ()02'0s T ϕφφφ-+=h=110 056φπ= '059φπ= ϕs =5π/92、推杆位移,速度,加速度线图如下(用matlab 编程得):A 、推杆位移线图clear clcx1=linspace(0,5*pi/6,300);x2=linspace(5*pi/6,13*pi/12,300); x3=linspace(13*pi/12,59*pi/36,300); x4=linspace(59*pi/36,2*pi,300); t1=x1/(5*pi/6)s1=110*(10*t1.^3-15*t1.^4+6*t1.^5); s2=110;t2=9*x3/(5*pi)-39/20;s3=110*(1-(10*t2.^3-15*t2.^4+6*t2.^5)); s4=0;plot(x1,s1,'k',x2,s2,'k',x3,s3,'k',x4,s4,'k') ; xlabel('角度/rad'); ylabel('位移s/mm'); title('推杆位移线图'); grid;B、推杆速度线图clearclcx1=linspace(0,5*pi/6,300);x2=linspace(5*pi/6,13*pi/12,300);x3=linspace(13*pi/12,59*pi/36,300);x4=linspace(59*pi/36,2*pi,300);f1=5*pi/6;t1=x1/f1;f2=5*pi/9;t2=9*x3/(5*pi)-39/20;v1=(t1.^2-2*t1.^3+t1.^4)*3300/f1;v2=0;v3=-30*110*(t2.^2-2*t2.^3+t2.^4)/f2;v4=0;plot(x1,v1,'k',x2,v2,'k',x3,v3,'k',x4,v4,'k') xlabel('角度/rad ');ylabel('速度v/(mm/s)');title('推杆速度线图');grid;C、推杆加速度线图clearclcx1=linspace(0,5*pi/6,300);x2=linspace(5*pi/6,13*pi/12,300);x3=linspace(13*pi/12,59*pi/36,300);x4=linspace(59*pi/36,2*pi,300);f1=5*pi/6;t1=x1/f1;f2=5*pi/9;t2=9*x3/(5*pi)-39/20;a1=60*110*(t1-3*t1.^2+2*t1.^3)/f1^2;a2=0;a3=-60*110*(t2-3*t2.^2+2*t2.^3)/f2^2;a4=0;plot(x1,a1,'k',x2,a2,'k',x3,a3,'k',x4,a4,'k') xlabel('角度/rad');ylabel('加速度a/ ');title('推杆加速度线图');grid;三、凸轮机构的ds/dψ---s线图,并依次确定凸轮的基圆半径和偏距1、凸轮机构的ds/dψ--s线图:clearclcx1=linspace(0,5*pi/6,300);x2=linspace(5*pi/6,13*pi/12,300);x3=linspace(13*pi/12,59*pi/36,300);x4=linspace(59*pi/36,2*pi,300);f2=5*pi/9;f1=5*pi/6;t1=x1/(5*pi/6)s1=110*(10*t1.^3-15*t1.^4+6*t1.^5);s2=110;t2=9*x3/(5*pi)-39/20;s3=110*(1-(10*t2.^3-15*t2.^4+6*t2.^5));s4=0;v1=(t1.^2-2*t1.^3+t1.^4)*3300/f1;v2=0;v3=-30*110*(t2.^2-2*t2.^3+t2.^4)/f2;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曲线');grid;2、确定凸轮的基圆半径和偏距:clearclcx1=linspace(0,5*pi/6,300);x2=linspace(5*pi/6,13*pi/12,300);x3=linspace(13*pi/12,59*pi/36,300);x4=linspace(59*pi/36,2*pi,300);f2=5*pi/9;f1=5*pi/6;t1=x1/(5*pi/6)s1=110*(10*t1.^3-15*t1.^4+6*t1.^5);s2=110;t2=9*x3/(5*pi)-39/20;s3=110*(1-(10*t2.^3-15*t2.^4+6*t2.^5));s4=0;v1=(t1.^2-2*t1.^3+t1.^4)*3300/f1;v2=0;v3=-30*110*(t2.^2-2*t2.^3+t2.^4)/f2;v4=0;k1=tan(pi/2-40*pi/180);k2=-tan(pi/6);f=sym('-k1*(2*k/f1^3-6*k^2/f1^4+4*k^3/f1^5)+k^2/f1^3-2*k^3/f1^4+k ^4/f1^5=0');k=solve(f);t01=k/f1;s01=110*(10*t01.^3-15*t01.^4+6*t01.^5);v01=(t01.^2-2*t01.^3+t01.^4)*3300/f1;c=80.5056;d=41.7790;%求出推程切点坐标x=-200:1:200;y5=k1*(x-c)+d;f2=5*pi/9;k2=-tan(pi/6);f=sym('-k2*(-2*(k*9/(5*pi)-39/20)*9/(5*pi)+6*(k*9/(5*pi)-39/20)^2 *9/(5*pi)-4*(k*9/(5*pi)-39/20)^3*9/(5*pi))-(k*9/(5*pi)-39/20)^2+2 *(k*9/(5*pi)-39/20)^3-(k*9/(5*pi)-39/20)^4=0');k=solve(f);t02=k*9/(5*pi)-39/20;s02=110*(1-(10*t02.^3-15*t02.^4+6*t02.^5));v02=-30*110*(t02.^2-2*t02.^3+t02.^4)/f2;o=32.1715;p= -112.4712;%求出回程切点坐标y6=k2*(x-p)+o;y7=x*-k1;plot(v1,s1,v2,s2,v3,s3,v4,s4,x,y5,x,y6,x,y7);xlabel('ds/dψ');ylabel('位移s/mm');title(' ds/dψ—s曲线');grid;所以,由图就可以确定回转中心所在的区域,所以,可取偏距 e=20mm , 080s = mm , 所以 222200802082.46r s e =+=+= mm 。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 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.运动规律(等加速等减速运动) 推程 0450≤≤ϕ 推程 009045≤≤ϕ2.运动规律(等加速等减速运动) 回程 00200160≤≤ϕ 回程 00240200≤≤ϕ三.推杆位移、速度、加速度线图及凸轮s d ds-φ线图采用VB 编程,其源程序及图像如下: 1.位移:Private Sub Command1_Click() = True '开启计时器 End SubPrivate Sub Timer1_Timer() Static i As SingleDim s As Single, q As Single 'i作为静态变量,控制流程;s代表位移;q代表角度= 0= 0i = i +If i <= 45 Thenq = is = 240 * (q / 90) ^ 2Step(q, -s), vbRedElseIf i >= 45 And i <= 90 Thenq = is = 120 - 240 * ((90 - q) ^ 2) / (90 ^ 2)Step(q, -s), vbGreenElseIf i >= 90 And i <= 150 Thenq = is = 120Step(q, -s), vbBlackElseIf i >= 150 And i <= 190 Thenq = is = 120 - 240 * (q - 150) ^ 2 / 6400Step(q, -s), vbBlueElseIf i >= 190 And i <= 230 Thens = 240 * (230 - q) ^ 2 / 6400Step(q, -s), vbRedElseIf i >= 230 And i <= 360 Thenq = is = 0Step(q, -s), vbBlackElseEnd IfEnd Sub2.速度Private Sub Command2_Click()= True '开启计时器End SubPrivate Sub Timer2_Timer()Static i As SingleDim v As Single, q As Single, w As Single 'i为静态变量,控制流程;q代表角度;w代表角速度,此处被赋予50= 0= 0w = 50i = i +If i <= 45 Thenv = 480 * w * q / 8100Step(q, -v), vbRedElseIf i >= 45 And i <= 90 Thenq = iv = 480 * w * (90 - q) / 8100Step(q, -v), vbBlackElseIf i >= 90 And i <= 150 Then q = iv = 0Step(q, -v), vbGreenElseIf i >= 150 And i <= 190 Then q = iv = -480 * w * (q - 150) / 6400 Step(q, -v), vbBlueElseIf i >= 190 And i <= 230 Then q = iv = -480 * w * (230 - q) / 6400 Step(q, -v), vbRedElseIf i >= 230 And i <= 360 Then q = iv = 0ElseEnd IfEnd Sub3.加速度Private Sub Command3_Click()= True '开启计时器End SubPrivate Sub Timer3_Timer()Static i As SingleDim a As Single, w As Single, q As Single 'i为静态变量,控制流程;a代表加速度;q代表角度;w代表角速度w = 50= 0= 0i = i +If i <= 45 Thenq = ia = 480 * w ^ 2 / 8100Step(q, -a), vbRedElseIf i >= 45 And i <= 90 Thenq = ia = -480 * w ^ 2 / 8100ElseIf i >= 90 And i <= 150 Then q = ia = 0Step(q, -a), vbGreenElseIf i >= 150 And i <= 190 Then q = ia = -480 * w ^ 2 / 6400Step(q, -a), vbBlueElseIf i >= 190 And i <= 230 Then q = ia = 480 * w ^ 2 / 6400Step(q, -a), vbRedElseIf i >= 230 And i <= 360 Then q = ia = 0Step(q, -a), vbBlackElseEnd IfEnd Subdq---dsPrivate Sub Command4_Click()= True '开启计时器;建立坐标系(-400, -400)-(400, 400)End SubPrivate Sub Timer4_Timer()Static i As SingleDim x As Single, s As Single, q As Single, scaley As Single, t As Single 'i为静态变量,控制流程;x代表位移;s代表纵坐标ds/dq;q代表角度= 0= 0scaley = 1t = / 180i = i +If i <= 45 Thenq = i * tx = * qs = 240 * (2 * q / ^ 2Step(x, -s), vbRedElseIf i >= 45 And i <= 90 Thenq = i * tx = * / 2 - q)s = 120 - * / 2 - q) ^ 2Step(x, -s), vbRedElseIf i >= 90 And i <= 150 Thenq = i * tx = 0s = 120 * scaleyStep(x, -s), vbRedElseIf i >= 150 And i <= 190 Then q = i * tx = * (q - 5 * / 6)s = 120 - * (q - 5 * / 6) ^ 2 Step(x, -s), vbRedElseIf i >= 190 And i <= 230 Then q = i * tx = * (23 * / 18 - q)s = * (23 * / 18 - q) ^ 2Step(x, -s), vbRedElseIf i >= 230 And i <= 360 Then q = i * tx = 0s = 0Step(x, -s), vbRedElseEnd IfEnd Sub四.确定凸轮基圆半径和偏距1. 求切点转角在图中,右侧曲线为升程阶段的类速度-位移图,作直线Dt dt与其相切,且位移轴正方向呈夹角[ 1]=300,则切点处的斜率与直线D t d t的斜率相等,因为k Dtdt=tan600,右侧曲线斜率可以表示为:q;q=tan600继而求出切点坐标(,)。

同理可求出另一切点(,)2.确定直线方程=*()=*(x+)3.绘图确定基圆半径和偏距程序如下:Private Sub Command5_Click() '开启计时器;画出ds/dq-s图的切线= True(-400, -400)-(400, 400)(-50, -(200, , vbRed(-200, -(40, , vbBlue(-200, -(200, , vbBlack(-200, -(200, , vbGreenEnd SubPrivate Sub Timer5_Timer() 'i为静态变量,控制流程;x代表位移;s代表纵坐标ds/dq;q代表角度;scaley是图形比例系数;t为角度,弧度转化系数Static i As SingleDim x As Single, s As Single, q As Single, scaley As Single, t As Single = 0= 0scaley = 1t = / 180i = i +If i <= 45 Thenq = i * tx = * qs = 240 * (2 * q / ^ 2Step(x, -s), vbRedElseIf i >= 45 And i <= 90 Thenq = i * tx = * / 2 - q)s = 120 - * / 2 - q) ^ 2Step(x, -s), vbRedElseIf i >= 90 And i <= 150 Thenq = i * tx = 0s = 120 * scaleyStep(x, -s), vbRedElseIf i >= 150 And i <= 190 Then q = i * tx = * (q - 5 * / 6)s = 120 - * (q - 5 * / 6) ^ 2 Step(x, -s), vbRedElseIf i >= 190 And i <= 230 Then q = i * tx = * (23 * / 18 - q)s = * (23 * / 18 - q) ^ 2Step(x, -s), vbRedElseIf i >= 230 And i <= 360 Thenq = i * tx = 0s = 0Step(x, -s), vbRedElseEnd IfEnd Sub在轴心公共许用区内取轴心位置,能够满足压力角要求,现取直线红线与直线绿线的= 偏距e= s0=交点为轴心位置,通过解方程组:x= y= r五.绘制凸轮理论廓线和实际廓线1.绘制凸轮理论轮廓Private Sub Command6_Click() '开启计时器;建立坐标系(-500, -500)-(500, 500)= TrueEnd SubPrivate Sub Timer6_Timer() 'i为静态变量,控制流程;e为偏距;s为基圆半径平方减去偏距平方再开方所得;q为角度;x为横坐标;y为纵坐标Static i As SingleDim e As Single, q As Single, s As Single, x As Single, y As Singlee =s ===i = i +If i <= / 4 Thenq = ix = (s + 240 * (2 * q / ^ 2) * Cos(q) - e * Sin(q)y = (s + 240 * (2 * q / ^ 2) * Sin(q) + e * Cos(q)Step(x, -y), vbRedElseIf i >= / 4 And i <= / 2 Thenq = ix = (s + 120 - * / 2 - q) ^ 2) * Cos(q) - e * Sin(q)y = (s + 120 - * / 2 - q) ^ 2) * Sin(q) + e * Cos(q)Step(x, -y), vbRedElseIf i >= / 2 And i <= 5 * / 6 Thenq = ix = (s + 120) * Cos(q) - e * Sin(q)y = (s + 120) * Sin(q) + e * Cos(q)Step(x, -y), vbRedElseIf i >= 5 * / 6 And i <= 19 * / 18 Thenq = ix = (s + 120 - * (q - 5 / 6 * ^ 2) * Cos(q) - e * Sin(q) y = (s + 120 - * (q - 5 / 6 * ^ 2) * Sin(q) + e * Cos(q)Step(x, -y), vbRedElseIf i >= 19 / 18 * And i <= 23 / 18 * Thenq = ix = (s + * (23 / 18 * - q) ^ 2) * Cos(q) - e * Sin(q)y = (s + * (23 / 18 * - q) ^ 2) * Sin(q) + e * Cos(q)Step(x, -y), vbRedElseIf i >= 23 / 18 * And i <= 2 * Thenq = ix = s * Cos(q) - e * Sin(q)y = s * Sin(q) + e * Cos(q)Step(x, -y), vbRedElseEnd IfEnd Sub2.确定滚子半径并绘制曲率半径图Private Sub Command7_Click() '开启计时器,建立坐标系= True(-10, -500)-(10, 10)End SubPrivate Sub Timer7_Timer() 'i为静态变量,控制流程;x代表dx/dq;y代表dy/dq;p代表曲率;e为偏距;s为基圆半径平方减去偏距平方再开方所得Static i As SingleDim s As Single, q As Single, x As Single, y As Single, p As Singles =e == 0= -Sqr(e ^ 2 + s ^ 2)i = i +If i <= / 4 Thenq = ix = (480 * 4 * q / ^ 2) * Cos(q) + (s + 240 * (2 * q / ^ 2) * (-Sin(q)) - e * Cos(q)y = (480 * 4 * q / ^ 2) * Sin(q) - e * Sin(q) + (s + 240 * (2 * q / ^ 2) * Cos(q)p = Sqr(x ^ 2 + y ^ 2)Step(q, -p), vbRedElseIf i >= / 4 And i <= / 2 Thenq = ix = (2 * * / 2 - q)) * Cos(q) + (s + 120 - * / 2 - q) ^ 2) * (-Sin(q)) - e * Cos(q)y = (2 * * / 2 - q)) * Sin(q) + (s + 120 - * / 2 - q) ^ 2) * Cos(q) - e * Sin(q)p = Sqr(x ^ 2 + y ^ 2)Step(q, -p), vbRedElseIf i >= / 2 And i <= 5 / 6 * Thenx = (s + 120) * (-Sin(q)) - e * Cos(q)y = (s + 120) * Cos(q) - e * Sin(q)p = Sqr(x ^ 2 + y ^ 2)Step(q, -p), vbRedElseIf i >= 5 / 6 * And i <= 19 / 18 * Thenq = ix = (-2 * * (q - 5 / 6 * ) * Cos(q) + (s + 120 - * (q - 5 / 6 * ^ 2) * (-Sin(q)) - e * Cos(q)y = (-2 * * (q - 5 / 6 * ) * Sin(q) + (s + 120 - * (q - 5 / 6 * ^ 2) * Cos(q) - e * Sin(q)p = Sqr(x ^ 2 + y ^ 2)Step(q, -p), vbRedElseIf i >= 19 / 18 * And i <= 23 / 18 * Thenq = ix = (-2 * * (23 / 18 * - q)) * Cos(q) + (s + * (23 / 18 * - q) ^ 2) * (-Sin(q)) - e * Cos(q)y = (-2 * * (23 / 18 * - q)) * Sin(q) + (s + * (23 / 18 * - q) ^ 2) * Cos(q) - e * Sin(q)p = Sqr(x ^ 2 + y ^ 2)Step(q, -p), vbRedElseIf i >= 23 / 18 * And i <= 2 * Thenx = s * (-Sin(q)) - e * Cos(q)y = s * Cos(q) - e * Sin(q)p = Sqr(x ^ 2 + y ^ 2)Step(q, -p), vbRedElseEnd IfEnd Sub其中曲率最小值为可以得到滚子半径为3.绘制实际轮廓线Private Sub Command8_Click()(-500, -500)-(500, 500) '建立坐标系;开启计时器= TrueEnd SubPrivate Sub Timer8_Timer() 'i为静态变量,控制流程;q代表角度;x代表dx/dy;y代表dy/dq;a代表理论轮廓的横坐标;b代表理论轮廓的纵坐标;m代表横坐标;n代表纵坐标Static i As SingleDim s As Single, q As Single, x As Single, y As Single, p As Single, rr As Single, a As Single, b As Single, m As Single, n As Singlerr =s =e ===i = i +If i <= / 4 Thenq = ix = (480 * 4 * q / ^ 2) * Cos(q) + (s + 240 * (2 * q / ^ 2) * (-Sin(q)) - e * Cos(q)y = (480 * 4 * q / ^ 2) * Sin(q) - e * Sin(q) + (s + 240 * (2 * q / ^2) * Cos(q)a = (s + 240 * (2 * q / ^ 2) * Cos(q) - e * Sin(q)b = (s + 240 * (2 * q / ^ 2) * Sin(q) + e * Cos(q)p = Sqr(x ^ 2 + y ^ 2)m = a - rr * y / pn = b + rr * x / pStep(m, -n), vbBlackElseIf i >= / 4 And i <= / 2 Thenq = ix = (2 * * / 2 - q)) * Cos(q) + (s + 120 - * / 2 - q) ^ 2) * (-Sin(q)) - e * Cos(q)y = (2 * * / 2 - q)) * Sin(q) + (s + 120 - * / 2 - q) ^ 2) * Cos(q) - e * Sin(q)a = (s + 120 - * / 2 - q) ^ 2) * Cos(q) - e * Sin(q)b = (s + 120 - * / 2 - q) ^ 2) * Sin(q) + e * Cos(q)p = Sqr(x ^ 2 + y ^ 2)m = a - rr * y / pn = b + rr * x / pStep(m, -n), vbBlackElseIf i >= / 2 And i <= 5 / 6 * Thenq = ix = (s + 120) * (-Sin(q)) - e * Cos(q)y = (s + 120) * Cos(q) - e * Sin(q)a = (s + 120) * Cos(q) - e * Sin(q)b = (s + 120) * Sin(q) + e * Cos(q)p = Sqr(x ^ 2 + y ^ 2)m = a - rr * y / pn = b + rr * x / pStep(m, -n), vbBlackElseIf i >= 5 / 6 * And i <= 19 / 18 * Thenq = ix = (-2 * * (q - 5 / 6 * ) * Cos(q) + (s + 120 - * (q - 5 / 6 * ^ 2) * (-Sin(q)) - e * Cos(q)y = (-2 * * (q - 5 / 6 * ) * Sin(q) + (s + 120 - * (q - 5 / 6 * ^ 2) * Cos(q) - e * Sin(q)a = (s + 120 - * (q - 5 / 6 * ^ 2) * Cos(q) - e * Sin(q)b = (s + 120 - * (q - 5 / 6 * ^ 2) * Sin(q) + e * Cos(q)p = Sqr(x ^ 2 + y ^ 2)m = a - rr * y / pn = b + rr * x / pStep(m, -n), vbBlackElseIf i >= 19 / 18 * And i <= 23 / 18 * Thenq = ix = (-2 * * (23 / 18 * - q)) * Cos(q) + (s + * (23 / 18 * - q) ^ 2) * (-Sin(q)) - e * Cos(q)y = (-2 * * (23 / 18 * - q)) * Sin(q) + (s + * (23 / 18 * - q) ^ 2) * Cos(q) - e * Sin(q)a = (s + * (23 / 18 * - q) ^ 2) * Cos(q) - e * Sin(q)b = (s + * (23 / 18 * - q) ^ 2) * Sin(q) + e * Cos(q)p = Sqr(x ^ 2 + y ^ 2)m = a - rr * y / pn = b + rr * x / pStep(m, -n), vbBlackElseIf i >= 23 / 18 * And i <= 2 * Thenq = ix = s * (-Sin(q)) - e * Cos(q)y = s * Cos(q) - e * Sin(q)a = s * Cos(q) - e * Sin(q)b = s * Sin(q) + e * Cos(q)p = Sqr(x ^ 2 + y ^ 2)m = a - rr * y / pn = b + rr * x / pStep(m, -n), vbBlackElseEnd IfEnd Sub内部轮廓为实际轮廓,外部轮廓为理论轮廓。

相关文档
最新文档