MATLAB大作业讲解.docx

偏置直动滚子推杆盘形凸轮机构设计推程运动规律：余弦加速度运动回程运动规律：正弦加速度运动凸轮机构的推杆在近休，推程，远休及回程段的凸轮转角近休凸轮转角推程凸轮转角远休凸轮转角回程凸轮转角0°~45°45°~210°210°~260°260°~360°偏置直动滚子推杆盘形凸轮机构已知参数初选基圆半?径ro/mm偏距e/mm滚子半径rr/mm推杆行程?h/mm许用压力角许用最小曲率半径[Pmina][a1][a2] 2070°0.3rr推杆运动规律推程运动规律：余弦加速度运动推程运动方程式：s =?h［1 - cos(πδ/δ0)］/ 2,v?= πhωsin(πδ/δ0)/( 2δ0),a?= π2hω2cos(πδ/δ0)/( 2δ02)?回程运动规律：正弦加速度运动回程运动方程式：s = h [1-(δ/δ0’)+sin(2πδ/δ0’)/(2π) v = hω [cos(2πδ/δ0’-1)]/δ0’a = -2πhω2sin(2πδ/δ0’)/(δ0’2)滚子中心坐标：x=(s+s0)sinδ+ecosδy=(s+s0)cosδ-ecosδ其中实际轮廓曲线坐标：其中计算压力角：计算曲率半径：计算程序符号的表示含义：r0 初选基圆半径x 理论轮廓曲线的x坐标y 理论轮廓曲线的y坐标X实际轮廓曲线的x坐标Y实际轮廓曲线的x坐标p曲率半径s推杆位移rr滚子半径a压力角paamin许用最小曲率半径pamin实际曲率半径的最小值a1max回程前压力角的最大值a2max回程后压力角的最大值e偏距h推杆行程a1推程许用压力角a2回程程许用压力角dar0每次修正r0所加的长度N计算点数k1近休角k2推程角k3远休角k4回程角kamin实际曲率半径最小值所对应的角度k1max回程前压力角的最大值所对应的角度k2max回程后压力角的最大值所对应的角度子程序1计算最大压力角：function [amax,kmax]=yalijiao1(a,N1,b)amax=0;kmax=0; amax=max(a);for i=1:N1 if a(i)==amaxkmax=i*b; endendend子程序2计算最小曲率半径：function [pamin,kamin]=qulv(p,b,N)pamin=0;kamin=0;for i=1:N-1 p(i)=abs(p(i));endpamin=min(p);for i=1:N-1pamin=min(p) if p(i)==paminkamin=i*b; endendend子程序3计算回程前的压力角及修正r0：function [A,r0]=haha(r0,rr,e,h,a1,dar0,N,k1,k2,k3,k4) N1=(k1+k2)/360*N;eta=pi/180;%化为弧度while(1)s0=(r0^2-e^2)^(1/2);syms k s dk ds d2k d2s n1=k1/360*N; k11=linspace(0,k1,n1); s1=k11./k11.*0; x1=e.*cos(k11.*pi/180)+(s0+s1).*sin(k11.*pi/180); y1=(s0+s1).*cos(k11.*pi/180)-e.*sin(k11.*pi/180); n2=k2/360*N; k22=linspace(k1,k1+k2,n2); s2=(h/2).*(1-cos(pi.*(k22-k1)/k2)); x2=e.*cos(k22.*pi/180)+(s0+s2).*sin(k22.*pi/180); y2=(s0+s2).*cos(k22.*pi/180)-e.*sin(k22.*pi/180); x=[x1,x2]; y=[y1,y2]; k=[k11,k22]; s=[s1,s2]; dx=diff(x)./diff(k);dk=k(1:end-1);dy=diff(y)./diff(k);dk=k(1:end-1); ds=diff(s)./diff(k);dk=k(1:end-1); ds(N1)=0; for i=1:N1 f(i)=(ds(i)-e)/(sqrt(r0^2-e^2)+s(i));a(i)=atan(f(i));%计算a压力角 A(i)=abs(a(i))/eta;if A(i)a1 r0=r0+dar0;break;end