2/16/07

9.3 兩輪接觸滾動

兩輪接觸滾動,雖然在機動學中不常見,但其基本觀念則是齒輪的運動狀態。

這個範例所用到的程式函數與前面之單圓滾動之情形相同,但由於兩輪的直徑可能不同,兩輪的中心可以置於水平線上。表面上看來,這兩個輪子僅是作固定位置旋轉,不必作任何平移動作。實際則不然,因為座標的原點不可能放在兩個輪的軸心上,若然,亦僅能將其中一個置於原點,另一個還是要進行平移。

本題目係將座標置於兩輪之接觸點上,因此各輪均要平移與轉動。目前因為兩輪均是圓的,故轉動起來,即使採用座標轉動的變化也看不出來,當然投機取巧不讓它轉動也是可以的,但將來若表面是有齒輪面的話,也必須表現出來。

為使表現能夠明顯。因此在各輪上設有一個徑向把,可以看出其真正轉動情形。其初始動作如圖4所示。

這個函數有三個參數,前二者為二輪之半徑,後者為r1轉動的速度。

<3>程式內容
這個程式呼叫circle函數,也是繪圓,與前面的circ功能相同,也是屬不同的寫法。讀者可以印證其功能。兩主圓分別為h1及h2,分別使用藍色與紅色;且為突顯其轉動狀態,特別改用虛線。其次是利用h11與h21分別設定兩桿,注意各桿之兩端均有標記,其型式也是以圓圈。

控制兩圓迴轉之速度為s1與s2,其關係與兩圓之半徑有關,故可以利用下式計算:

theta1=theta1+s1;
theta2=theta2-s1*r1/r2;

程式內容


每個圓圈之迴轉則利用tans4這個副程式。這個程式內容在上一節中已經有討論過,可以參考。亦可自行撰寫簡單的程式應用。

function demo_mot7(R1,R2,v)
%demo_mot7.m
% Inputs:R1,R2:radius of two circles
% v: speed of rotation
%Two matched cicles in motion
% Example: demo_mot7(50,100,2)
%Author:DS Fon, Bime. NTU Date: Jan. 24, 2007
if nargin==0,R1=50;R2=100;v=100;end
clf;
r1=min(abs(R1),abs(R2));
r2=max(abs(R1),abs(R2));
height=max(r1,r2)*1.5;
axis([-2.5*r1 2.4*r2 -height height]);
axis image;axis off;
line([-2*r1 2*r2]',[0 0]','color','k',...
'linewidth',1,'linestyle','-.'); %The ground line
circle(r1/4,r2,0,10);circle(r1/4,-r1,0,10);
h1=circle(r1,0,0); %Draw the 1st circle of radius r1
h2=circle(r2,0,0); %Draw the 2nd circle of radius r2
h11=line([0 r1]',[0 0]','erasemode','xor','color','b',...
'linewidth',4,'marker','o','markersize',8); %Rod 1
h21=line([0 -r2]',[0 0]','erasemode','xor','color','r',...
'linewidth',8,'marker','o','markersize',16); %Rod 2
hold on;
%Get coordinates (x1,y1) & (x2,y2) of both circles
x1=get(h1,'xdata');y1=get(h1,'ydata');
x2=get(h2,'xdata');y2=get(h2,'ydata');
x11=get(h11,'xdata');y11=get(h11,'ydata');% rod 1
x21=get(h21,'xdata');y21=get(h21,'ydata');% rod 2
set(h1,'EraseMode','xor','color','b',...
'linewidth',1,'linestyle','-.')
set(h2,'EraseMode','xor','color','r',...
'linewidth',2,'linestyle','-.')
grid on;
title('Press Ctl-C to stop');
theta1=0;theta2=0;s1=v*pi/20;
while 1,
drawnow;
pause(0.1);
theta1=theta1+s1;
theta2=theta2-s1*r1/r2;
z1=trans4([x1',y1'],theta1,5);%rotate circle 1
z1=trans4(z1,[-r1 0],1);%translate circle 1
z2=trans4([x2',y2'],theta2,5);%rotate circle 2
z2=trans4(z2,[r2 0],1);%translate circle 2 to (r2,0)
z11=trans4([x11',y11'],theta1,5);%rotate rod 1
z11=trans4(z11,[-r1 0],1);%translate rod 1 to (-r1,0)
z21=trans4([x21',y21'],theta2,5);%rotate rod 2
z21=trans4(z21,[r2 0],1);%translate rod 2 to (r2,0)
set(h1,'xdata',z1(:,1),'ydata',z1(:,2)); % circle 1 data
set(h2,'xdata',z2(:,1),'ydata',z2(:,2)); % Circle 2 data
set(h11,'xdata',z11(:,1),'ydata',z11(:,2)); % Rod 1
set(h21,'xdata',z21(:,1),'ydata',z21(:,2)); % Rod 2
if theta1>2*pi,theta1=theta1-2*pi;s1=-s1;end;
end

function h = circle(r,x0,y0,nn)
% The inputs:
% r = radius of circle
% x0, y0= coordinates of the circular center
% Example: circle(10,0,0,10)
if nargin==3,nn=50;end;
t=0:2*pi/nn:2*pi;
h=line(x0+r*cos(t),y0+r*sin(t));


執行結果





看看這個例子,同學可以思考未來還可加一些什麼,因為可用的指令都有了。想一想較困難的吧,諸如:

  1. 把一個人字形的人圖像放在圓圈中央,讓它耍把戯。
  2. 試著在輪外圍加上齒輪的外形,其外形僅作一個就好,其餘用轉動拷貝的,如此就可以作成實際的齒輪動作了。
  3. 能否加上第三個輪,其變換方向又有不同。
  4. 能否作一個齒輪組。
  5. 不知將圓輪變成多角形時會怎樣?

No comments:

Post a Comment