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

("大作业（二）凸轮机构设计（题号：8）班级：姓名、学号：成绩：完成日期：目录1.凸轮机构大作业题目……………………………………22.推杆运动规律及凸轮廓线方程…………………………33.程序流程图………………………………………………34.源程序…………………………………………………^55.计算结果…………………………………………………146.凸轮机构图………………………………………………167.体会及建议………………………………………………198.参考资料…………………………………………………20一、凸轮机构大作业题目试用计算机辅助设计完成下列摆动滚子推杆盘形凸轮机构的设计，已知数据如下表所示，凸轮沿着逆时针方向做匀速转动。表1凸轮机构的推杆运动规律推程运动规律正弦加速度运动回程运动规律等加速等减速运动表2两种凸轮机构的推杆在近休、推程、远休及回程阶段的凸轮转角题号近休凸轮转角推程凸轮转角远休凸轮转角回程凸轮转角A0゜～60゜60゜～180゜180゜～270゜270゜～360゜B0゜～45゜45゜～210゜210゜～260゜260゜～360゜C0゜～30゜30゜～210゜210゜～280゜280゜～360゜表3摆动滚子推杆盘形凸轮机构的已知参数题号初选基圆半径r0/mm机架长度lOA/mm摆杆长度lAB/mm滚子半径rr/mm推杆摆角许用压力角许用最小曲率半径[ρamin][α1][α2]A1560551024゜35゜70゜0.3rrB2070651426゜40゜70゜0.3rrC227275681828゜45゜65゜0.35rr要求：每组（每三人为一组，每人一题）至少打印出一份源程序，每人打印出原始数据；凸轮理论轮廓和实际轮廓的坐标值；推程和回程的最大压力角，以及出现最大压力角时凸轮的相应转角；凸轮实际轮廓曲线的最小曲率半径，以及相应的凸轮转角；和最后说确定的基圆半径。计算点数N=72~120。绘出凸轮的理论轮廓和实际轮廓（可用计算机绘图）。二、推杆运动规律及凸轮廓线方程：推程（正弦加速度）：s=h[(δ/δ0)-sin(2πδ/δ0)/(2π)]回程（等加速段）：s=h-2hδ2/δ'02回程（等减速段）：s=2h(δ'0-δ)2/δ'02凸轮理论廓线方程：x=lOAsinδ-lABsin(δ+φ+φ0)y=lOAcosδ-lABcos(δ+φ+φ0)式中，φ0为推杆的初始位置角，其值为:φ0=arccosα>[α1]？开始读入：r0，Δr0，rt，h或（φ），e或（lAB、lOA）δ1，δ2，δ3，δ4，[α1]，[α2]，[ρamin]，N计算：s0I=1计算：s，x，y，ds/dδ，dx/dδ，dy/dδ，x′，y′计算：αr0=r0+Δr0r0=r0=Δr0是回程？α>[α2]？选出α1max及相应的凸轮转角δ1选出α2max及相应的凸轮转角δ2计算：ρρ<0?ρ-rt≥[ρamin]？计算ρa选出ρamin及相应的凸轮转角δaminI=I+1I≤N？三、程序流程图打印：x，y，x′，y′，ρamin，δamin，α1max，δ1max，α2max，δ2max，r0，δ，s结束四、源程序clear；r0=22;%初选的基圆半径dr0=0.05;a=72;%机架长度L=68;%摆杆长度rr=18;%滚子半径fai=28pi/180;%推杆摆角PI=3.141592653;alpha1=45;%许用压力角α1alpha2=65;%许用压力角阿尔法2lambda=6.3;%许用最小曲率半径N=120;%取用点的个数delta1=180pi/180;%推程凸轮最大转角delta2=70pi/180;%远休凸轮最大转角delta3=80pi/180;%回程凸轮最大转角delta4=30pi/180;%近休凸轮最大转角alphamax1=0;%推程最大压力角初值alphamax2=0;%回程最大压力角初值roumin=100;%凸轮最小曲率半径初值X=ones(1,121);Y=ones(1,121);XP=ones(1,121);YP=ones(1,121);FAI=ones(1,121);I=0;%初始化while(I<=N)F0=acos((aa+LL-r0r0)/(2aL));delta=3Ipi/180;if(delta<180pi/180&delta>=0)%正弦推程F0=acos((aa+LL-r0r0)/(2aL));F=fai((delta/delta1)-sin(2pidelta/delta1)/(2pi));x=asin(delta)-Lsin(delta+F+F0);y=acos(delta)-Lcos(delta+F+F0);dF=fai(1/delta1-cos(2pidelta/delta1)/delta1);ddF=fai2pisin(2pidelta/delta1)/(delta1delta1);dx=acos(delta)-Lcos(delta+F+F0)(1+dF);ddx=-asin(delta)+Lsin(delta+F+F0)(1+dF)^2-Lcos(delta+F+F0)ddF;dy=-asin(delta)+Lsin(delta+F+F0)(1+dF);ddy=-acos(delta)+Lcos(delta+F+F0)(1+dF)^2+Lsin(delta+F+F0)ddF;stheta=dx/sqrt(dxdx+dydy);ctheta=-dy/sqrt(dxdx+dydy);xp=x-rrctheta;yp=y-rrstheta;alpha=atan((Labs(dF)+(acos(F+F0)-L))/(asin(F+F0)))180/pi;alpha=abs(alpha);if(alpha>alphamax1)alphamax1=alpha;deltamax1=delta180/pi;endendif(delta>=180pi/180&delta<250pi/180)%远休F0=acos((aa+LL-r0r0)/(2aL));F=fai;x=asin(delta)-Lsin(delta+F+F0);y=acos(delta)-Lcos(delta+F+F0);dF=0;ddF=0;dx=acos(delta)-Lcos(delta+F+F0)(1+dF);ddx=-asin(delta)+Lsin(delta+F+F0)(1+dF)^2-Lcos(delta+F+F0)ddF;dy=-asin(delta)+Lsin(delta+F+F0)(1+dF);ddy=-acos(delta)+Lcos(delta+F+F0)(1+dF)^2+Lsin(delta+F+F0)ddF;stheta=dx/sqrt(dxdx+dydy);ctheta=-dy/sqrt(dxdx+dydy);xp=x-rrctheta;yp=y-rrstheta;endif(delta>=250pi/180&delta<290pi/180)%等加速回程F0=acos((aa+LL-r0r0)/(2aL));F=fai-2fai(delta-250pi/180)(delta-250pi/180)/(delta3delta3);x=asin(delta)-Lsin(delta+F+F0);y=acos(delta)-Lcos(delta+F+F0);dF=-4fai(delta-250pi/180)/(delta3delta3);ddF=-4fai/(delta3delta3);dx=acos(delta)-Lcos(delta+F+F0)(1+dF);ddx=-asin(delta)+Lsin(delta+F+F0)(1+dF)^2-Lcos(delta+F+F0)ddF;dy=-asin(delta)+Lsin(delta+F+F0)(1+dF);ddy=-acos(delta)+Lcos(delta+F+F0)(1+dF)^2+Lsin(delta+F+F0)ddF;stheta=dx/sqrt(dxdx+dydy);ctheta=-dy/sqrt(dxdx+dydy);xp=x-rrctheta;yp=y-rrstheta;alpha=atan((Labs(dF)-(acos(F+F0)-L))/(asin(F+F0)))180/pi;alpha=abs(alpha);if(alpha>alphamax2)alphamax2=alpha;deltamax2=delta180/pi;endendif(delta>=290pi/180&delta<330pi/180)%等减速回程F0=acos((aa+LL-r0r0)/(2aL));F=2fai(delta3-(delta-250pi/180))(delta3-(delta-250pi/180))/(delta3delta3);dF=-4fai(delta3-(delta-250pi/180))/(delta3delta3);ddF=4fai/(delta3delta3);x=asin(delta)-Lsin(delta+F+F0);y=acos(delta)-Lcos(delta+F+F0);dx=acos(delta)-Lcos(delta+F+F0)(1+dF);ddx=-asin(delta)+Lsin(delta+F+F0)(1+dF)^2-Lcos(delta+F+F0)ddF;dy=-asin(delta)+Lsin(delta+F+F0)(1+dF);ddy=-acos(delta)+Lcos(delta+F+F0)(1+dF)^2+Lsin(delta+F+F0)ddF;stheta=dx/sqrt(dxdx+dydy);ctheta=-dy/sqrt(dxdx+dydy);xp=x-rrctheta;yp=y-rrstheta;alpha=atan((Labs(dF)-(acos(F+F0)-L))/(asin(F+F0)))180/pi;alpha=abs(alpha);if(alpha>alphamax2)alphamax2=alpha;deltamax2=delta180/pi;endendif(delta>=330pi/180&delta<=360pi/180)%近休F0=acos((aa+LL-r0r0)/(2aL));F=0;x=asin(delta)-Lsin(delta+F+F0);y=acos(delta)-Lcos(delta+F+F0);dF=0;ddF=0;dx=acos(delta)-Lcos(delta+F+F0)(1+dF);ddx=-asin(delta)+Lsin(delta+F+F0)(1+dF)^2-Lcos(delta+F+F0)ddF;dy=-asin(delta)+Lsin(delta+F+F0)(1+dF);ddy=-acos(delta)+Lcos(delta+F+F0)(1+dF)^2+Lsin(delta+F+F0)ddF;stheta=dx/sqrt(dxdx+dydy);ctheta=-dy/sqrt(dxdx+dydy);xp=x-rrctheta;yp=y-rrstheta;endif(alphamax1>alpha1)(alphamax2>alpha2)%优化条件r0=r0+dr0;I=0;alphamax1=0;alphamax2=0;continueendrou=((dx^2+dy^2)^(3/2))/(dxddy-dyddx);%计算曲率半径if(rou<0)rou=-rou;if((rou-rr)<(0.35rr))%优化条件r0=r0+dr0;I=0;alphamax1=0;alphamax2=0;continueendif(rou