matlab法畫凸輪輪廓線

1、班級： 姓名： 學(xué)號：基于matlab的凸輪輪廓設(shè)計一、 設(shè)計凸輪機構(gòu)的意義在工業(yè)生產(chǎn)中，經(jīng)常要求機器的某些部件按照規(guī)定的準確路線運動，僅應(yīng)用連桿機構(gòu)已難以滿足這個要求，所以需要利用工作表面具有一定形狀的凸輪。凸輪在所有基本運動鏈中，具有易于設(shè)計和能準確預(yù)測所產(chǎn)生的運動的優(yōu)點。如果設(shè)計其他機構(gòu)來產(chǎn)生給定的運功、速度、和加速度，其設(shè)計工作是很復(fù)雜的，但是設(shè)計凸輪機構(gòu)則比較容易，而且運動準確、有效。所以在許多機器中，如紡織機、包裝機、自動機床、自動化專用機床、數(shù)控機床、印刷機、內(nèi)燃機、建筑機械、礦山機械、計算機的輔助裝備及農(nóng)業(yè)機具等，都可以找到凸輪機構(gòu)。在進行研究時，先設(shè)計一個簡單的凸輪，在給定的

2、旋轉(zhuǎn)角度內(nèi)有一定的總升距。設(shè)計凸輪輪廓的基本方法是把凸輪固定，使從動件以其與凸輪的相關(guān)位置繞凸輪回轉(zhuǎn)而形成凸輪輪廓。因此設(shè)計凸輪時，必須畫出足夠多的點，使凸輪輪廓平滑可靠。Matlab軟件提供了強大的矩陣處理和繪圖功能，具有核心函數(shù)工具箱。其編程代碼接近數(shù)學(xué)推導(dǎo)公式，簡潔直觀，操作簡易，人機交互性能好。因此，基于matlab軟件進行凸輪機構(gòu)的設(shè)計可以解決設(shè)計工作量大的問題。運用解析法進行設(shè)計，matlab可以精確的計算出輪廓上每一點的坐標，然后更為精確的繪制出凸輪輪廓曲線。二、設(shè)計凸輪機構(gòu)的已知條件凸輪做逆時針方向轉(zhuǎn)動，從動件偏置在凸輪軸心右邊。從動件在推程做等加/減速運動，在回程做余弦加速運

3、動?；鶊A半徑rb=50mm，滾子半徑rt=10mm，推桿偏距e=10mm，推程升程h=50mm，推程運動角ft=100º，遠休止角fs=60º，回程運動角fh=90º。三、 分析計算1、建立坐標系 以凸輪軸心為坐標原點建立平面直角坐標系XOY，取桿件上升方向為Y軸正方向。2、推桿運動規(guī)律計算凸輪運動一周可分為5個階段：推程加速階段、推程減速階段、遠休止階段、回程階段、進休止階段。根據(jù)已知條件，推程階段為等加/減速，故推程階段的運動方程為：推程加速階段（0ft/2）sf=2*h*f2ft2vf=dsf=4*h*fft2推程減速階段（ft/2ft）sf=h-

4、2*h*ft-f2ft2vf=dsf=4*h*ft-fft2遠休止階段（ftft+fs）推桿運動方程為sf=0vf=dsf=0根據(jù)已知條件，在回程做余弦加速運動，因此回程階段（ft+fsft+fs+fh）的運動方程為sf=h*1+cos*f-ft-fsfh2vf=dsf=-*h*sin*f-ft-fsfh2*fh近休止階段（ft+fs+fh360°）的運動方程為sf=0vf=dsf=03、凸輪理論輪廓線計算xf=se+sf*sinf+e*cosfyf=se+sf*cosf-e*sinf式中se=rb2-e2，為推桿滾子中心到X軸的垂直距離。4、實際輪廓線計算根據(jù)3的計算結(jié)果有dxf=

5、dsf-e*sinf+se+sf*cosfdyf=dsf-e*cosf-se+sf*sinf可得sink=dxfdx2+dy2cosk=-dyfdx2+dy2凸輪實際輪廓線為xpf=xf-r*coskypf=yf-r*sink四、程序代碼 rb = 50; rt = 10; e = 10; h = 50; ft = 100; fs = 60; fh = 90; hd= pi / 180; du = 180 / pi; se=sqrt( rb2 - e2 ); d1 = ft + fs; d2 = ft + fs + fh; n = 360; s = zeros(n); ds = zeros(n

6、); x = zeros(n); y = zeros(n); dx = zeros(n); dy = zeros(n); xx = zeros(n); yy = zeros(n); xp = zeros(n); yp = zeros(n); for f = 1 : n if f <= ft/2 s(f) = 2 * h * f 2 / ft 2; s = s(f); ds(f) = 4 * h * f * hd / (ft * hd) 2; ds = ds(f); elseif f > ft/2 & f <= ft s(f) = h - 2 * h * (ft - f

7、) 2 / ft 2; s = s(f); ds(f) = 4 * h * (ft - f) * hd / (ft * hd) 2; ds = ds(f); elseif f > ft & f <= d1 s = h;ds = 0; elseif f > d1 & f <= d2 k = f - d1; s(f) = .5 * h * (1 + cos(pi * k / fh); s = s(f); ds(f)= -.5 * pi * h * sin(pi * k / fh) / (fh * hd); ds = ds(f); elseif f >

8、d2 & f <= n s = 0;ds = 0; end xx(f) = (se + s) * sin(f * hd) + e * cos(f * hd); x = xx(f); yy(f) = (se + s) * cos(f * hd) - e * sin(f * hd); y = yy(f); dx(f) = (ds - e) * sin(f * hd) + (se + s) * cos(f * hd); dx = dx(f); dy(f) = (ds - e) * cos(f * hd) - (se + s) * sin(f * hd); dy = dy(f); xp(

9、f) = x + rt * dy / sqrt(dx 2 + dy 2);xxp = xp(f); yp(f) = y - rt * dx / sqrt(dx 2 + dy 2);yyp = yp(f); end disp ' 凸輪轉(zhuǎn)角 理論x 理論y 實際x 實際y' for f = 10 : 10 :ft nu = f xx(f) yy(f) xp(f) yp(f); disp(nu) end disp ' 凸輪轉(zhuǎn)角 理論x 理論y 實際x 實際y' for f = d1 : 10 : d2 nu = f xx(f) yy(f) xp(f) yp(f); d

10、isp(nu) end plot(xx,yy,'r-.') axis (-(rb+h-10) (rb+h+10) -(rb+h+10) (rb+rt+10) axis equal text(rb+h+3,0,'X') text(0,rb+rt+3,'Y') text(-5,5,'O') title('偏置移動從動件盤形凸輪設(shè)計') hold on; plot(-(rb+h) (rb+h),0 0,'k') plot(0 0,-(rb+h) (rb+rt),'k') plot(e e,0 (rb+rt),'k-') ct = linspace(0,2*pi);