凸轮机构设计 (1),凸轮机构

("《凸轮机构设计》目录·······················姓名：学号：学院（一）、题目及原始数据···············（二）、推杆运动规律及凸轮廓线方程·········（三）、计算程序方框图···············（四）、计算源程序·················（五）、程序计算结果及分析·············（六）、凸轮机构图·················（七）、心得体会··················（八）、参考书···················（一）、题目及原始数据试用计算机辅助设计完成偏置直动滚子推杆盘形凸轮机构的设计，凸轮以1rad/s的角速度沿逆时针方向转动。要求：（1）、推程运动规律为等加速等减速运动，回程运动规律为五次多项式运动规律；（2）、打印出原始数据；（3）、打印出理论轮廓和实际轮廓的坐标值；（4）、打印出推程和回程的最大压力角，以及出现最大压力角时凸轮的相应转角；（5）、打印出凸轮实际轮廓曲线的最小曲率半径，以及相应的凸轮转角；（6）、打印出凸轮运动的位移；（7）、打印最后所确定的凸轮的基圆半径。原始数据如下：r0=0.015;初选的基圆半径r0Deltar0=0.0005;当许用压力角或许用最小曲率半径不满足时，r0以Δr0为步长增加重新计算rr=0.010;滚子半径rrh=0.028;推杆行程he=0.005;偏距eomega=1;原动件凸轮运动角速度，逆时针ωdelta1=pi/3;近休止角δ1delta2=2pi/3;推程运动角δ2delta3=pi/2;远休止角δ3delta4=pi/2;回程运动角δ4alpha1=pi/6;推程许用压力角[α1]alpha2=(70/180)pi;回程许用压力角[α2]rho0min=0.3rr;许用最小曲率半径ραmin（二）、推杆运动规律及凸轮廓线方程推杆运动规律：（1）近休阶段：0o≤δ<60os=0v=0a=0（2）推程阶段：60o≤δ<180o等加速运动规律：60o≤δ<120os=2h(δ-60o)2/(120o)2v=4hω(δ-60o)/(120o)2a=4hω2/(120o)2等减速运动规律：120o≤δ<180os=h-2h(120o-(δ-60o))2/(120o)2v=4hω(120o-(δ-60o))/(120o)2a=-4hω2/(120o)2（3）远休阶段：180o≤δ<270os=hv=0a=0（4）回程阶段：270o≤δ≤360o五次多项式运动规律：s=h-(10h(δ-270o)3/(90o)3-15h(δ-270o)4/(90o)4+6h(δ-270o)5/(90o)5)v=-(30hω(δ-270o)2/(90o)3-60hω(δ-270o)3/(90o)4+30hω(δ-270o)4/(90o)5)a=-(60hω2(δ-270o)/(90o)3-180hω2(δ-270o)2/(90o)4+120hω2(δ-270o)3/(90o)5)凸轮廓线方程：（1）理论廓线方程：s0=sqrt(r02-e2)x=(s0+s)sinδ+ecosδy=(s0+s)cosδ-esinδ(2)实际廓线方程先求x，y的一阶导数x’=(v/ω-e)sinδ+(s0+s)cosδy’=(v/ω-e)cosδ-(s0+s)sinδ再求sinθ，cosθsinθ＝x’/sqrt((x’)2+(y’)2)cosθ＝-y’/sqrt((x’)2+(y’)2)最后求实际廓线方程x1=x-rrcosθy1=y-rrsinθ（四）、计算源程序%凸轮机构大作业Matlab语言源程序%选题：偏置直动滚子推杆盘形凸轮机构5—A%推程运动规律：等加速等减速运动%回程运动规律：五次多项式运动%作者:罗永恒学号：2008301390076%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clearall;%closeall;clc;%赋初值r0=0.015;Deltar0=0.0005;rr=0.010;h=0.028;e=0.005;omega=1;%原动件凸轮运动角速度，逆时针delta1=pi/3;%近休止角delta2=2pi/3;%推程运动角delta3=pi/2;%远休止角delta4=pi/2;%回程运动角alpha1=pi/6;%推程许用压力角alpha2=(70/180)pi;%回程许用压力角rho0min=0.3rr;%许用最小曲率半径del1=delta1+delta2;del2=del1+delta3;temp=0;%判断是否执行r0＝r0＋Deltar0的变量while(temp==0)temp=1;s0=sqrt(r0^2-e^2);%求解s0alpha1max=0;delta1max=0;%定义alpha1的最大值以及对应的delta1值alpha2max=0;delta2max=0;%定义alpha2的最大值以及对应的delta2值rhoamin=r0-rr;deltamin=0;%定义rhoa的最小值以及对应的delta值forI=0:120;%圆周120等分delta=(I3/180)pi;ifdelta>=0&delta=delta1&delta<(delta2/2)+delta1%等加速推程s=2h(delta-delta1)^2/delta2^2;v=4homega(delta-delta1)/delta2^2;a=4homega^2/delta2^2;elseifdelta>=(delta2/2)+delta1&delta=del1&delta=del2&delta<=2pi%五次多项式运动规律回程s=h-(10h(delta-del2)^3/delta3^3-15h(delta-del2)^4/delta3^4+6h(delta-del2)^5/delta3^5);v=-(30homega(delta-del2)^2/delta4^3-60homega(delta-del2)^3/delta4^4+30homega(delta-del2)^4/delta4^5);a=-(60homega^2(delta-del2)/delta4^3-180homega^2(delta-del2)^2/delta4^4+120homega(delta-del2)^3/delta4^5);endx=(s0+s)sin(delta)+ecos(delta);%理论轮廓方程式y=(s0+s)cos(delta)-esin(delta);x_=(v/omega-e)sin(delta)+(s0+s)cos(delta);%理论轮廓对delta求一次导数y_=(v/omega-e)cos(delta)-(s0+s)sin(delta);x__=(a/omega^2-(s0+s))sin(delta)+(2v/omega-e)cos(delta);%理论轮廓对delta求二次导数y__=(a/omega^2-(s0+s))cos(delta)-(2v/omega--e)sin(delta);x1=x-rr(-y_/sqrt(x_^2+y_^2));%实际轮廓方程式y1=y-rr(x_/sqrt(x_^2+y_^2));alpha=atan((v-e)/(sqrt(r0^2-e^2)+s));%求压力角ifdelta>=del2&delta<=2pi%判断是否为回程ifabs(alpha)>alpha2%判断是否大于回程许用压力角r0=r0+Deltar0;temp=0;break;elseifabs(alpha)>alpha2max%满足许用压力角，则找出回程最大压力角alpha2max=abs(alpha);delta2max=delta;endendelseifabs(alpha)>alpha1%判断是否大于推程许用压力角r0=r0+Deltar0;%不满足许用压力角，则增大基圆半径重新计算temp=0;break;elseifabs(alpha)>alpha1max%满足许用压力角，则找出推程最大压力角alpha1max=abs(alpha);delta1max=delta;endendendrho=(x_^2+y_^2)^(3/2)/(x_y__-y_x__);%计算曲率半径ifrho<0rhoa=abs(rho)-rr;ifrhoa>=rho0min%满足最小曲率半径ifrhoa