基于MATLAB的凸輪設(shè)計(jì)

1、0 2/35 = 2。1/ = 0TI2/311320160S =- 140 H承 3 -不1T / 320 320V = - - + 61r /a =-320/tta25 MOrt中國(guó)地質(zhì)大學(xué)（武漢）1 .凸輪要求設(shè)計(jì)一對(duì)心直動(dòng)滾子推桿盤形凸輪機(jī)構(gòu)，滾子半徑r二10mm凸輪以等角速度逆時(shí)針回轉(zhuǎn)。凸輪轉(zhuǎn)角5=0120時(shí)，推桿等速上升20mm 5=120180 時(shí)，推桿遠(yuǎn)休止；書=180270時(shí)，推桿等加速等減速下降20mm耳=270360 時(shí)，推桿近休止。要求推程的最大壓力角4二30,試選取合適的基圓半徑，并 繪制凸輪的廓線。問此凸輪是否有缺陷，應(yīng)如何補(bǔ)救。2 .列出凸輪運(yùn)動(dòng)方程30S = +

2、 67T30V -L Q = 04801G0S = 360 -* 6 + 干* 8A2 兄nz4B0 320V = -* 635/4R 2i7r /a=- 320/7TA23 .由方程寫MATLAB源程序%i.已知參數(shù)clear;r0=50;%基圓半徑rr=10;%滾子半徑h=20;%亍程delta01=120;%推程運(yùn)動(dòng)角delta02=60;% 遠(yuǎn)休角delta03=90;沏程運(yùn)動(dòng)角hd=pi/180;du=180/pi;n1=delta01+delta02;n2=delta01+delta02+delta03;%2凸輪曲線設(shè)計(jì)n=360;for i=1:360%計(jì)算推桿運(yùn)動(dòng)規(guī)律if id

3、elta01 & in1 & i(n1+delta03/2) & in2 & i=ns(i)=0;ds=0;end%計(jì)算凸輪軌跡曲線xx(i)=(r0+s(i)*sin(i*hd);%計(jì)算理論輪廓曲線yy(i)=(r0+s(i)*cos(i*hd);dx(i)=ds*sin(i*hd)+(r0+s(i)*cos(i*hd);%計(jì)算導(dǎo)數(shù)dy(i)=ds*cos(i*hd)-(r0+s(i)*sin(i*hd);xp(i)=xx(i)+rr*dy(i)/sqrt(dx(i)A2+dy(i)A2);yp(i)=yy(i)-rr*dx(i)/sqrt(dx(i)A2+dy(i)A2);end%3. 輸

4、出凸輪輪廓曲線figure(1);holdon ;gridon ;axis equal ;axis(-(r0+h-30) (r0+h+10) -(r0+h+10) (r0+rr+10);text(r0+h+3,4,X );text(3,r0+rr+3,Y );text(-6,4,O );title( 對(duì)心直動(dòng)滾子推桿盤形凸輪設(shè)計(jì) );xlabel(x/mm);ylabel(y/mm);plot(-(r0+h-40) (r0+h),0 0,k );plot(0 0,-(r0+h) (r0+rr),k );plot(xx,yy,r-); %?繪凸輪實(shí)際輪廓曲線ct=linspace(0,2*pi)

5、;plot(r0*cos(ct),r0*sin(ct),plot(rr*cos(ct),r0+rr*sin(ct), plot(0,r0, o ); %滾子圓中心 plot(0 0,r0 r0+30,k );plot(xp,yp, b );%繪凸輪實(shí)際輪廓曲線%4. 凸輪機(jī)構(gòu)運(yùn)動(dòng)仿真 %計(jì)算凸輪滾子轉(zhuǎn)角 xp0=0;yp0=r0-rr;dss=sqrt(diff(xp).A2+diff(yp)-2);ss(1)=sqrt(xp-xp02+(xp-yp02);for i=1:359ss(i+1)=ss(i)+dss(i); endphi=ss/rr;%計(jì)算滾子轉(zhuǎn)角%運(yùn)動(dòng)仿真開始figure(2)

6、;m=moviein(20);j=0;for i=1:360j=j+1;delta(i)=i*hd;%凸輪轉(zhuǎn)角xy=xp,yp;%凸輪實(shí)際輪廓曲線坐標(biāo)A1=cos(delta(i),sin(delta(i);-sin(delta(i),cos(delta(i);xy=xy*A1;%旋轉(zhuǎn)后實(shí)際凸輪曲線坐標(biāo)clf;%繪凸輪plot(xy(:,1),xy(:,2);hold on ;axis equal ;axis(-(120) (470) -(100) (140);plot(-(r0+h-40) (r0+h),0, plot(0 0,-(r0+h) (r0+rr), plot(r0*cos(ct)

7、,r0*sin(ct), plot(rr*cos(ct),r0+s(i)+rr*sin(ct),%計(jì)算實(shí)際廓曲線長(zhǎng)度g ); %繪凸輪基圓k ); %繪滾子圓%對(duì)輪廓曲線進(jìn)行差分計(jì)算%輪廓曲線第一點(diǎn)長(zhǎng)度%凸輪坐標(biāo)旋轉(zhuǎn)矩陣k); %繪凸輪水平軸k); %繪凸輪垂直軸g ); %繪基圓k ); 繪滾子圓plot(0 rr*cos(-phi(i),r0+s(i) r0+s(i)+rr*sin(-phi(i),k );% 繪滾子圓標(biāo)線plot(0 0,r0+s(i) r0+s(i)+40,k ); %繪推桿%繪推桿曲線plot(1:360+r0+h,s+r0);plot(r0+h) (r0+h+360

8、),r0 r0,k );plot(r0+h) (r0+h),r0 r0+h,k );plot(i+r0+h,s(i)+r0,* );title( 對(duì)心直動(dòng)滾子推桿盤形凸輪設(shè)計(jì) );xlabel(x/mm);ylable(y/mm);m(j)=getframe;endmovie(m);4 . 運(yùn)動(dòng)仿真結(jié)果在MATLA沖可以看出輪廓曲線有一處缺口。應(yīng)用圓弧連接起來(lái)5 .計(jì)算結(jié)果由于數(shù)據(jù)太多，只等間隔取了三十六組數(shù)據(jù)Xx8.971818.241127.500036.424644.685951.961557.947762.371265.000065.653964.212360.621853.62314

9、4.995135.000023.941412.15540.0000-12.0696-23.2658-32.7778-39.9163-44.3549-47.1503-48.8408-49.7267-50.0000-49.2404-46.9846-43.3013-38.3022-32.1394-25.0000-17.1010-8.6824-0.0000Yy50.881750.116947.631443.409237.495930.000021.091210.99770.0000-11.5765- 23.3714-35.0000-44.9951-53.6231-60.6218-65.7785-68.

10、9365-70.0000-68.4502-63.9223-56.7728-47.5704-37.2182-27.2222-17.7766-8.7682-0.00008.682417.101025.000032.139438.302243.301346.984649.240450.0000Dx52.540053.383052.406049.547444.811138.269930.064620.40199.5493-2.1723-14.3980-26.7301-44.9951-53.6231-60.6218-65.7785-68.9365- 70.0000-67.4676-60.0514-48.2845-33.0207-19.8785-12.5201-7.1415-3.19530.00008.682417.101025.000032.139438.302243.301346.984649.240450.0000Dy0.4324-9.2677-19.2301-29.1094-38.5478-47.1869-54.6817-60.7129-65.0000-67.3121-67.4784-65.3964-53.