机械原理课程设计偏置直动滚子推杆盘形凸轮

凸轮大作业选题：凸轮5-C一：题目及原始数据：运用计算机辅助设计完毕偏置直动滚子推杆盘形凸轮机构旳设计，已知数据如下所示，凸轮沿逆时针方向做匀速运动。具体规定如下：1.推程运动规律为等加速等减速运动，回程运动规律为五次多项式运动；2.近休凸轮转角为0°-30°；推程凸轮转角30°-210°；远休凸轮转角210°-280°；回程凸轮转角280°-360。°3.初选基圆半径为22mm；4.偏距为+14mm5.滚子半径为18mm6.推杆行程为35mm7.许用压力角为α1=35°，α2=65°。8.最小曲率半径为0.35rr9.计算点数取120.二：推杆运动规律及凸轮轮廓线方程1.推程加速阶段：s1=70.*a1.*a1/pi/pi;x1=(s0+s1).*sin(a1)+e*cos(a1);y1=(s0+s1).*cos(a1)-e*sin(a1);k1=140*a1/pi^2;i1=[(k1-e).*sin(a1)+(s0+s1).*cos(a1)].*[(k1-e).*(k1-e)+(s0+s1).*(s0+s1)].^(-1/2);j1=[-(k1-e).*cos(a1)+(s0+s1).*sin(a1)].*[(k1-e).*(k1-e)+(s0+s1).*(s0+s1)].^(-1/2);x10=x1-18*j1;y10=y1-18*i1;2.推程减速阶段：s11=35-70.*(pi-a11).*(pi-a11)/pi/pi;x11=(s0+s11).*sin(a11)+e*cos(a11);y11=(s0+s11).*cos(a11)-e*sin(a11);k11=140.*(pi-a1)/pi^2;i11=[(k11-e).*sin(a11)+(s0+s11).*cos(a11)].*[(k11-e).*(k11-e)+(s0+s11).*(s0+s11)].^(-1/2);j11=[-(k11-e).*cos(a11)+(s0+s11).*sin(a11)].*[(k11-e).*(k11-e)+(s0+s11).*(s0+s11)].^(-1/2);x101=x11-18*j11;y101=y11-18*i11;3.远休阶段：x2=(s0+s2).*sin(a2)+e*cos(a2);y2=(s0+s2).*cos(a2)-e*sin(a2);k2=0;i2=[(k2-e).*sin(a2)+(s0+s2).*cos(a2)].*[(k2-e).*(k2-e)+(s0+s2).*(s0+s2)].^(-1/2);j2=[-(k2-e).*cos(a2)+(s0+s2).*sin(a2)].*[(k2-e).*(k2-e)+(s0+s2).*(s0+s2)].^(-1/2);x20=x2-18*j2;y20=y2-18*i2;4.回程阶段：a30=33*pi/18-a3;s3=(350.*a30.*a30.*a30/(4*pi/9)^3-525.*a30.*a30.*a30.*a30/(4*pi/9)^4+210.*a30.*a30.*a30.*a30.*a30/(4*pi/9)^5);x3=(s0+s3).*sin(a3)+e*cos(a3);y3=(s0+s3).*cos(a3)-e*sin(a3);k3=-1050.*a30.*a30/(4*pi/9)^3+2100.*a30.*a30.*a30/(4*pi/9)^4-1050.*a30.*a30.*a30.*a30/(4*pi/9)^5;i3=[(k3-e).*sin(a3)+(s0+s3).*cos(a3)].*[(k3-e).*(k3-e)+(s0+s3).*(s0+s3)].^(-1/2);j3=[-(k3-e).*cos(a3)+(s0+s3).*sin(a3)].*[(k3-e).*(k3-e)+(s0+s3).*(s0+s3)].^(-1/2);x30=x3-18*j3;y30=y3-18*i3;5.近休阶段：s4=0;x4=(s0+s4).*sin(a4)+e*cos(a4);y4=(s0+s4).*cos(a4)-e*sin(a4);k4=0;i4=[(k4-e).*sin(a4)+(s0+s4).*cos(a4)].*[(k4-e).*(k4-e)+(s0+s4).*(s0+s4)].^(-1/2);j4=[-(k4-e).*cos(a4)+(s0+s4).*sin(a4)].*[(k4-e).*(k4-e)+(s0+s4).*(s0+s4)].^(-1/2);x40=x4-18*j4;y40=y4-18*i4;三：matlab计算程序e=14;r0=22s0=sqrt(r0^2-e^2);c111=1;c3=1;for(i=1:0.1:200)%由压力角条件循环求合适基圆半径if(c111<35/180)&&(c3<65/180)%判断条件elsea1=0:pi/60:pi/2;%推程加速阶段s1=70.*a1.*a1/pi/pi;x1=(s0+s1).*sin(a1)+e*cos(a1);y1=(s0+s1).*cos(a1)-e*sin(a1);k1=140*a1/pi^2;%对s1求导i1=[(k1-e).*sin(a1)+(s0+s1).*cos(a1)].*[(k1-e).*(k1-e)+(s0+s1).*(s0+s1)].^(-1/2);j1=[-(k1-e).*cos(a1)+(s0+s1).*sin(a1)].*[(k1-e).*(k1-e)+(s0+s1).*(s0+s1)].^(-1/2);x10=x1-18*j1;y10=y1-18*i1;a11=pi/2:pi/60:pi;%推程减速阶段s11=35-70.*(pi-a11).*(pi-a11)/pi/pi;x11=(s0+s11).*sin(a11)+e*cos(a11);y11=(s0+s11).*cos(a11)-e*sin(a11);k11=140.*(pi-a1)/pi^2;i11=[(k11-e).*sin(a11)+(s0+s11).*cos(a11)].*[(k11-e).*(k11-e)+(s0+s11).*(s0+s11)].^(-1/2);j11=[-(k11-e).*cos(a11)+(s0+s11).*sin(a11)].*[(k11-e).*(k11-e)+(s0+s11).*(s0+s11)].^(-1/2);x101=x11-18*j11;y101=y11-18*i11;a2=pi:pi/60:25*pi/18;%凸轮远休阶段s2=35;%推杆行程x2=(s0+s2).*sin(a2)+e*cos(a2);y2=(s0+s2).*cos(a2)-e*sin(a2);k2=0;i2=[(k2-e).*sin(a2)+(s0+s2).*cos(a2)].*[(k2-e).*(k2-e)+(s0+s2).*(s0+s2)].^(-1/2);j2=[-(k2-e).*cos(a2)+(s0+s2).*sin(a2)].*[(k2-e).*(k2-e)+(s0+s2).*(s0+s2)].^(-1/2);x20=x2-18*j2;y20=y2-18*i2;a3=25*pi/18:pi/60:33*pi/18;%推杆回程阶段a30=33*pi/18-a3;s3=(350.*a30.*a30.*a30/(4*pi/9)^3-525.*a30.*a30.*a30.*a30/(4*pi/9)^4+210.*a30.*a30.*a30.*a30.*a30/(4*pi/9)^5);x3=(s0+s3).*sin(a3)+e*cos(a3);y3=(s0+s3).*cos(a3)-e*sin(a3);k3=-1050.*a30.*a30/(4*pi/9)^3+2100.*a30.*a30.*a30/(4*pi/9)^4-1050.*a30.*a30.*a30.*a30/(4*pi/9)^5;i3=[(k3-e).*sin(a3)+(s0+s3).*cos(a3)].*[(k3-e).*(k3-e)+(s0+s3).*(s0+s3)].^(-1/2);j3=[-(k3-e).*cos(a3)+(s0+s3).*sin(a3)].*[(k3-e).*(k3-e)+(s0+s3).*(s0+s3)].^(-1/2);x30=x3-18*j3;y30=y3-18*i3;a4=33*pi/18:pi/60:2*pi;%推杆近休阶段s4=0;x4=(s0+s4).*sin(a4)+e*cos(a4);y4=(s0+s4).*cos(a4)-e*sin(a4);k4=0;i4=[(k4-e).*sin(a4)+(s0+s4).*cos(a4)].*[(k4-e).*(k4-e)+(s0+s4).*(s0+s4)].^(-1/2);j4=[-(k4-e).*cos(a4)+(s0+s4).*sin(a4)].*[(k4-e).*(k4-e)+(s0+s4).*(s0+s4)].^(-1/2);x40=x4-18*j4;y40=y4-18*i4;plot(x10,y10,'-g*',x101,y101,'-g+',x20,y20,'-r*',x30,y30,'-b*',x40,y40,'-k*',x1,y1,'-g*',x11,y11,'-g+',x2,y2,'-r*',x3,y3,'-b*',x4,y4,'-k*')%凸轮轮廓曲线绘制title(‘凸轮轮廓曲线绘制');xlabel('VariableX');%X轴ylabel('VariableY');%Y轴text(-250,-200,'工作廓线')%文字标注text(100,-100,'理论廓线');gridon%加网格axisequal%坐标相等a1=0:pi/60:pi/2;%压力角计算force1=abs(atan((k1).*(s0+s1).^-1));c1=max(force1);a11=pi/2:pi/60:pi;force11=abs(atan((k11).*(s0+s11).^-1));c11=max(force11);c111=max(c1,c11);a2=pi:pi/60:25*pi/18;force2=abs(atan((k2).*(s0+s2).^-1));c2=max(force2);a3=25*pi/18:pi/60:33*pi/18;force3=abs(atan((k3).*(s0+s3).^-1));c3=max(force3);a4=33*pi/18:pi/60:2*pi;force4=abs(atan((k4).*(s0+s4).^-1));c4=max(force4);r0=r0+1;%每循环一次基圆半径+1s0=sqrt(r0^2-e^2);endend%求最大压力角位置c111c3[m1,n1]=sort(force1);bend1=n1(end-1+1:end)jiaodu1=(bend1(end)-1)*3[m11,n11]=sort(force11);bend11=n11(end-1+1:end)jiaodu11=(bend11(end)-1)*3+90[m3,n3]=sort(force3);bend3=n3(end-1+1:end)jiaodu=(bend3(end)-1)*3+250vv=i1./j1;vv1=diff(vv);vv2=diff(vv,2);vv22=[0,vv2];p=(1+vv1.^2).^(3/2)./vv22;g=min(p)vv1=i11./j11;vv11=diff(vv1);vv21=diff(vv1,2);vv221=[0,vv21];p1=(1+vv11.^2).^(3/2)./vv221;g1=min(p1)vv2=i2./j2;vv222=diff(vv2);vv223=diff(vv2,2);vv2211=[0,vv223];p2=(1+vv222.^2).^(3/2)./vv2211;g2=min(p2)vv=i3./j3;vv1=diff(vv);vv2=diff(vv,2);vv22=[0,vv2];p=(1+vv1.^2).^(3/2)./vv22;g3=min(p)vv=i4./j4;vv1=diff(vv);vv2=diff(vv,2);vv22=[0,vv2];p=(1+vv1.^2).^(3/2)./vv22;g4=min(p))ro=sqrt(s0^2+e^2)%求合适基圆半径x1y1x11y11x2y2x3y3x4y4x10y10x101y101x20y20x30y30x40y40force1force11force2force3force4c111c3jiaodu11jiaodu3四：计算成果及分析：（1）工作廓线坐标：1;推程加速阶段X坐标14 24.8 35.8 46.3 57.6 67.5 78.0 88.50 98.5 108.3 117. 127.14255136 136.6 144. 153.5 161.0 168. 176.0 182. 189. 195.6 200. 205. 210. 214. 217. 220. 223.4 225. 226.30872647 227.2推程加速阶段Y坐标209. 208. 206. 204. 202. 199. 195.6 191. 186.86029371 181. 176.7 170. 163. 156. 149.6 141. 133.3 124.5 115. 106.24932735 96.1 86.9 76.7 65.6 54.7 43.6 32.3 21.67 9.54 -2.948 -14.03推程减速阶段X坐标227. 227.6 226. 225. 223. 220. 217. 213. 209. 203. 198. 191. 184. 176. 168. 159.9 150. 140. 130. 120.3 109.6 97.7 86.64 74.5 62.0 49.8 37.0 24.5 11.1 -1.846 -14.04推程减速阶段Y坐标-14.0 -25.9 -37.3 -49.1 -61.3 -73.6 -85.43 -97.07 -108. -119. -130. -141. -151. -161. -170. -179.7 -188. -196. -203. -210. -217. -222. -227. -232.4 -236. -239. -241. -243. -244. -244. -244.5远休阶段X坐标-14.0 -26.2 -39.2 -52.3 -64.8 -76.2 -88.7 -100. -112. -123. -134. -144. -155. -164. -174. -182. -191. -198. -206. -212.70743182 -218. -224. -229. -233.6远休阶段Y坐标-244. -243. -241. -239. -236. -232. -228. -223. -217. -211. -204. -197. -189. -181. -172. -163. -153. -143. -132. -121. -110.8 -98.54 -86.2 -74.27减速阶段X坐标-234. -237. -240. -242. -243.4 -243. -242. -240.3 -237. -233. -228. -223. -217. -210. -203. -196.2 -188. -180. -171.6 -163. -154. -145. -136. -127. -118.9 -108. -99.688减速阶段Y坐标-70.0 -58.9 -45.54 -32.0 -20.3 -7.250 5.72 18.00 30.7 42.6 54.3 65.5 76.2 87.47 97.49 106. 115. 124. 132. 140. 147. 154. 161. 167. 173. 179. 185.19近休X坐标-92.1 -82.65 -72.0 -62.05 -51.1 -40.0 -29.9 -18.5 -7.64 3.076 13.910近休Y坐标1188. 193. 197. 200. 203. 206.0 207. 209.8 209. 209. 209.（2）理论廓线坐标;推程加速阶段X坐标12.0 22.4 32.6 42.6 52.66 62.2 71. 81.4 90.2 99.1 108. 116. 125. 133.13203958 140. 148. 155. 161. 168. 174.3 179. 184.7 189. 193. 197. 200. 203. 205. 207.8 208. 209.推程加速阶段Y坐标191. 190. 189. 187. 184. 182. 178. 174. 170. 165. 160. 155. 149. 142. 136. 128. 121.2 113. 104. 96.1 87.0 78.0 68.1 59.6 49.5 38.5 28.44 17.5 7.054 -3.10 -14.1推程减速阶段X坐标209.0 209.9 208. 207. 205. 203. 199. 196. 191. 187. 181. 175. 169. 161. 154.3 146. 137. 128. 119.4 109. 99. 88.6 77.2 66.85 55.1 44.69 32.6 20. 9.87 -2.37 -14.5推程减速阶段Y坐标-16.40 -27.33 -38.5 -49.8 -60.04 -71.6 -81.8 -92.1 -102. -113.08170624 -123. -132. -142.3 -151.3 -159.8 -168. -175. -183.9 -190.05865654 -196. -202.1 -207. -211. -215. -219.4 -222.2 -224. -225. -226. -226. -226.远休阶段X坐标-12.4 -24.4 -36.8 -48.4 -59.6 -71.91 -82.9 -93.1 -104. -114.4 -124.5 -134. -143. -152. -161. -169. -177. -184. -190.9 -197. -202. -207. -212. -216.4远休阶段Y坐标-226. -225. -223. -221. -218.9 -215. -211. -206. -201. -195. -189. -182. -175. -167. -159. -151. -141. -132. -122. -112.5 -102. -91.0 -80.30 -69.7回程阶段X坐标-217. -220. -222. -224. -225.5 -225.5 -224. -222. -219. -216. -212. -207.5 -202. -195. -189. -182. -174. -167. -159. -151. -142. -134. -125. -117. -108. -99.7 -90.0回程阶段Y坐标-65.2 -53.6 -42.6 -30.3 -19.8 -7.05 3.39 14.4 25.94 36. 47.8 57.2 67.1 77.5 86.0 95.7 103. 111. 119. 126. 133. 140. 146. 153. 158. 164.2 16