万能拉拨机凸轮设计.docx

例(万能拉拨机凸轮设计) 在万能拉拨机中有一个圆柱形凸轮， 其底圆半径 R=300mm,凸轮的上端面不在同一平面上 (如图)，而要 根据从动杆位移变化的需要进行设计制造。 根据设计要求，将底圆周 18等分，旋转一周。 图9.5万能拉拨机第i个分点对应柱高 yi, (i=0,1,2,… 图9.5万能拉拨机 i yi 0和18 502.8 1 525.0 2 514.3 3 451.0 4 326.5 5 188.6 i 6 7 8 9 10 11 yi 92.2 59.6 62.2 102.7 147.1 191.6 i 12 13 14 15 16 17 yi 236.0 280.5 324.9 369.4 413.8 458.3 我们将圆周展开，画出对应的柱高曲线 (图9.6)。 clear;close; x=li nspace(0,2*pi*300,19); y=[502.8 525.0 514.3 451.0 326.5 188.6 92.2 59.6 62.2 102.7… 147.1 191.6 236.0 280.5 324.9 369.4 413.8 458.3 502.8]; plot(x,y,o);axis([0,2000,0,550]); 可见柱高形成一条 V形曲线。现在的问题是，怎样给出分点之外的柱高呢？ 图9.6凸轮柱高数据 图9.6凸轮柱高数据 500 450 400 350 300 250 200 150 100 50 0 0 由于没有拟合经验公式，又要求严格按设计数据要求(而这些数据应认为是足够精确 的),所以应该用插值方法。又由于问题是一条 封闭曲线，所以考虑使用周期性端点条件样 条插值，它具有较好的光滑性。 clear;close; x=li nspace(0,2*pi*300,19); y=[502.8 525.0 514.3 451.0 326.5 188.6 92.2 59.6 62.2 102.7... 147.1 191.6 236.0 280.5 324.9 369.4 413.8 458.3 502.8]; pp=csape(x,y,periodic);fnplt(pp);axis([0,2000,0,550]); % 样条插值曲线图 9.7 pp.breaks,pp.coefs % 分段多项式 pp=mkpp(pp.breaks,pp.coefs); xi=0:2*pi*300;yi=ppval(pp,xi) % 插值数值结果 注意：封闭曲线要使用(1)角度作为自变量，如椭圆 (2)周期样条 图9.7凸轮样条插值 x=acosn y=bsinv