2007年6月16日 星期六

作業十三

1試設計一組複式齒輪,使其轉速比為125(請說明思考步驟及結果)。

由於每組之轉速比最好維持在10以內,超過此值時則需考慮增加齒輪組數。因此轉速比為125時必須用較多的組數才成達到目標。首先,將125開方,其值為11.1803,仍然比10大,故使用兩組組合仍嫌不足,如將125開立方,其值為5,低於10甚多,且為整數,故使用三組齒輪之組合,只要將這種組合串聯三組即可。設驅動之小齒輪數最小為12齒(N1>=12)*,則依序可以得到對應大齒輪之齒數如下:

N2 = N1 * 5 = 60

N3 = N2 * 5 = 300

所以取12齒,60齒以及300齒的齒輪作為串聯,可得到題目的要求
再配合matlab程式即可得到所設計之動畫。

for wa=0:0.1:2
clf
[coords]=draw_gear(5,25,20,360,0,0,wa)
wb=-5*wa
[coords]=draw_gear(5,5,20,360,3,0,wb)
[coords]=draw_gear(5,25,20,360,3,0,wb)
wc=-5*wb
[coords]=draw_gear(5,5,20,360,6,0,wc)
[coords]=draw_gear(5,25,20,360,6,0,wc)
wd=-5*wc
[coords]=draw_gear(5,5,20,360,9,0,wd)
pause(0.1)
end

2請指出本學期中你自己最感得意的一次作業(請說明其原因,且該作業必須在自己的部落格內)。
我想,在這學期最感得意的應該是就是第五次作業吧。那時剛接觸MATLAB沒多久,對MATLAB還有些害怕。剛開始看到陌生題目的時候,不知所措,那時候覺得只有小畫家可以畫出這樣的東西,之前學的C++都只有運算之類的,沒有畫圖,而AUTOCAD也不用寫程式。於是我再只上先自己畫出一之手臂,然後分析每點之間的關係,在輸入視窗。其實這步驟沒有想像中困難,就只是普通的算數。當我用line把手連起來時,我覺得自己學會了一個很不錯的軟體。
接下來,要寫一個可以輸入角度與每節長度的function。這花了我一些時間思索。後來把期間關係式寫出來後,發現怎麼樣都跑不出來。非常沮喪,好在同學提醒我需檢查”.”,原來要加了才可以被當作數字算。
還要讓這手臂轉動,天啊,那時我想我絕對做不出來。後來利用老師上課講的for迴圈,還有pause,發現沒幾個式子,就可以得到想要的動畫。我發現,只要將所想呈現的動畫,稍作分析,就可以畫出來。讓我實在佩服製造MATLAB軟體的人。之前因為要研究MATLAB,有到書局翻閱相關書籍,我發現,Matlab的運用範圍十分廣,有工程數學、商業、地質分析等等玲瑯滿目,我想這真是一個很好機會讓我們接觸這軟體。

2007年6月7日 星期四

作業十二




1. 本人5/31有來上課

2.一組標準全齒輪齒輪之徑節為8(亦可使用自設值),齒數分別為30T與48T,其工作壓力角為20度(可為14.5或25度,自選)。

•試求其接觸線長度,與接觸比。

在老師的課程網頁中有contact_ratio程式,
輸入參數:
Pd:徑節=8
n2, n3:兩齒輪之齒數=30,48
phi壓力角=20

輸出參數:
cr_ratio接觸比
cr_length:接觸長度
ad:齒冠
pc, pb:周節及基周節
d2, d3:兩齒輪節圓直徑。
ag:兩齒輪之接近角、遠退角及作用角

輸入[c_ratio,c_length,ad,pc,pb,d2,d3,ag]=contact_ratio(8,30,48,20 )
得到
c_ratio =

1.7005


c_length =

0.6275


ad =

0.1250


pc =

0.3927


pb =

0.3690


d2 =

3.7500


d3 =

6


ag =

10.4850 9.9211 20.4061
6.5532 6.2007 12.7538
(注意單位)
也就是
接觸比=1.7005
 接觸長度=0.6275吋
 齒冠=0.1250吋
 周節=0.3927
 基周節pb=0.3690
 齒輪節圓直徑齒輪一=3.75吋
       齒輪二=6吋
 齒輪一之接近角=10.4850度
     遠退角=9.9211度
     作用角=20.4061度
 齒輪二之接近角=6.5532度
     遠退角=6.2007度
     作用角=12.7538度

•兩齒輪之節圓、基圓直徑各為如何?請列式計算其結果。
我們可以由定義知道,
節徑*徑節=齒數
 所以齒輪一之節徑為30/8=3.75吋
   齒輪二之節徑為48/8=6吋

 而基圓直徑=節徑*cos(th)
(其中th為壓力角)
 所以齒輪一之基圓直徑為3.75*cos(20)=3.523吋
   齒輪二之基圓直徑為 6*cos(20)=5.638吋

•此組齒輪是否會產生干涉現象?試列式證明之。
在課本9-39中,有詳細介紹干涉的檢驗並且逐步推導。 
根據公式我們知道,
N2及N3要滿足
 (N2^2+2N2*N3)sin^2(壓力角) >= 4(1+N3)
才不發生干涉

  依照題意N2為30,N3為48,壓力角為20 則
 (20^2+2*20*48)sin^2(20) >= 4*(1+48)
271.1 > 196
  故本組數據將不會發生干涉
可以由
function [x]=isinterf(phi,N1,N2) 來測試

(若輸出為1代表有干涉現象,0則否)
所以我們輸入isinterf(20,30,48)
Ans=0


•可否利用draw_gear.m繪出其接合情形,並繪出其動畫效果。

function [coords]=draw_gear(Dp,N,phi,range,x0,y0)



其中,輸入參數分別定義如下:

Dp: 節矩

N: 齒數

phi: 壓力角

range: 繪出之部份

x0,y0: 齒輪中心座標

draw_gear(8,48,20,360,2.5,0)


可得到接合情形

  



另外,我們可以使用
function move2_gear(Dpitch,nn1,nn2,phi,omega1)


輸入參數之定義如下:


Dpitch:節矩
nn1,nn2:兩齒輪之齒數
phi:壓力角
omega1: 齒輪1之角速度,rad/s

在螢幕上輸入move2_gear(8,30,48,20,5)
即可得到下面動畫


2007年5月30日 星期三

作業十一

1.本人5/24有來上課,謝謝老師。

2.某凸輪開始時先在0-100∘區間滯留,然後提升後在200至260∘區間滯留,其高度(衝程)為5公分,其餘l由260∘至360∘則為返程。升程採用等加速度運動,返程之運動型式自定。設刻度區間為10∘,試繪出其高度、速度及加速度與凸輪迴轉角度間之關係。


在老師課程網站上第十章凸輪運動有對凸輪詳細的剖析以及介紹,這題我們使用老師的給的function plot_dwell

程式碼如下
function plot_dwell(ctheta,s,pattern,range)
figure(1);clf;
[y,yy,yyy]=dwell(ctheta,range,pattern)
%[y,yy,yyy]=dwell(ctheta,range,pattern)亦取自老師的課程網頁

h1=plot(ctheta,y*s,'b-',ctheta,yy*s,'k-',ctheta,yyy*s,'r-')
legend('Displacement','Velocity','Acceleration',3)
xlabel('Elapsed Angle, degrees')
grid

而所需輸入的參數如下:
*ctheta = 需要計算之凸輪角度,單位為度數。(可以使用矩陣輸入之型式)
*s=衝程,極為所上升之高度。
*pattern = 運動的型式,二元素之列矩陣,其代碼如下:
1:等速運動uniform
2:抛物線parabolic
3:簡諧simple harmonic
4:擺線cycloidal 
5:多項式polynomial motion
*range =升程及返程之範圍,以矩陣表示

我們分析此程式中[y,yy,yyy]=dwell(ctheta,range,pattern) function中,我們知道所輸出的値表示下列意思:
*y:位移
*yy:對於凸輪角ctheta之第一導數;
* yyy :對於凸輪角ctheta之第二導數。

因此,有了這些function,我們只要在螢幕上輸入

plot_dwell(0:10:360,5,[2,5],[100 200 260 360])

*0:10:360 凸輪從0到360度每隔十度來計算
*5 衝程
* [2 5] 2表示升程維等加速度拋物線運動
5表示回程自定的多項式運動
* [100 200 260 360]
標示升程與回程的角度

即可以得到位置速度加速度與角度之關係圖。





3.設凸輪之半徑為15公分,以順時針方向旋轉,其從動件為梢型,垂直接觸,長為10公分,從動件之運動係依照第二項之運動型式。試繪出此凸輪之工作曲線。

若要描繪凸輪的工作曲線,可以從老師的網頁中參考
function [x,y]=pincam(cth,r0,s,e,L,range,pattern,cw)
其中
*cth:凸輪角度,度數
*r0:凸輪基圓半徑
*e:偏置量
*s:從動件衝程
*L:從動件長度
*cw:凸輪轉動方向(反時鐘為正,順時鐘為負)
*pattern 運動的型式motion
*range 升程及返程之範圍

因此,依照題意,
在螢幕上輸入

[x,y]=pincam([0:10:360],15,5,0,10,[100 200 260 360],[2 3],-1)

可得到凸輪之工作曲線。



4. 你能讓此凸輪迴轉嗎?

要讓凸輪旋轉其實不難,只要將for迴圈妥善利用即可。我寫ㄧ個function pincam2,運用第三題老師所給的pincamrotate funtion,再加一些旋轉的程式碼分析其轉動情形。
我的程式碼如下:

function [x,y]=pincamrotate(cth,r0,s,e,L,range,pattern,cw)
%Find the pin type cam with an offsect e
%Inputs:
% cth:angle of cam, degrees
% r0:radius of base circle
% e:offset
% s:stroke
% L:length of pin
% cw:rotation direction of cam(-counterclockwise,+clockwise
%pattern = denote the type of motion used(a 3 element-row matrix)
% 1:uniform 2:parabolic 3:simple harmonic 4: cycloidal
% 5:polynomial motion
% example [4 3]
%range =the degrees the specific motion starts, eg.[90 180 240]
% Example: [x y]=pincam([10 60],5,2,1,10,[90 180 240],[4 3],-1)
clf;
th=cth*pi/180;
s0=sqrt(r0*r0-e*e);
for i=1:length(cth)
t=th(i)*cw;
A=[cos(t) -sin(t);sin(t) cos(t)];
[ym,yy,yyy]=dwell(cth(i),range,pattern);
x0=s0+ym*s;
Sx=[0 x0 x0+L;e e e];
X=A\Sx;
x(i)=X(1,2);y(i)=X(2,2);
end

%以下是我多加的,和第三題程式不同的地方

%寫一個for迴圈,從第一次跑到第36次
for n=1:36
%以旋轉公式寫出各點經旋轉後的座標

a=x+r0*cosd(10*n);
b=y+r0*sind(10*n);

%描畫各點,再將它們的軌跡以紅色o表達
plot([0 a],[0 b],'ro',a,b,'k-');

%設定座標範圍
axis ([-50 50 -50 50]);

%每次動作停留0.1秒
pause(0.1);
clf
end

依照題意,在螢幕上輸入

pincam2([0:10:360],15,5,0,10,[100 200 260],[2 1],-1)

可得動畫如下

2007年5月26日 星期六

作業十

1.本人5/17有來上課,謝謝老師。
2.
*請思考速度與加速度的問題,當一桿以某特定點M等角速度迴轉時,其端點P之速度方向如何?
其加速度方向如何?




設角速度為w,桿長r
P點的速度則為 w x r
方向為桿運動之切線方向

其加速度為w x w x r
其方向由P指回M點

*若該特定點M復以等速水平運動,則同一端點P之速度與加速度方向會變為如何?

此時的P點速度為 VM + w x r (其中V代表M以V速度在移動)
而加速度仍是為w x w x r(方向由P指回M點)
注意:這些加減並非純向量加減,必須考慮方向

*若M點同時也有加速度,則點P會有何變化?



M點有加速度,所以P點不再只有法線加速度而已,還有一個切線加速度。和M有水平速度不一樣的地方在於,M有加速度只有影響速度和加速度的大小而已。
AP=AM+ωx(ωxr)


*若以此推理四連桿的運動,則點P與Q之速度與加速度方向會與桿一(固定桿)之兩端點之關係如何?與我們前面的作業分析結果有無共通之處?



固定桿兩端為O點和R點著地,所以其速度與加速度皆為0,
若以第二桿為驅動桿,
則P點速度為
VP=Wop x Rp/o  方向與桿二垂直

AP=[αx Rp/o]+[Wop x (Wop x Rp/o)] 
前項為切線加速度與桿2垂直,後項為法線加速度指向O點


Q點的速度與加速度
由於P點有速度
所以Q點的速度為
Vq = Vp + Wpq x Rq/p
但是由於Vq 是硍R點連接 而R點為接地點
所以Q點的速度也為 Wq r x Rq/r  方向與桿四垂直

同樣加速度也是
Aq = Ap+α x Rq/p + Wpq x (Wpq x Rq/p)
也等於 α x Rq/r + Wqr x (Wqr x Rq/r) 
切線加速度與桿四垂直,法線加速度指向R點



3.設有一運動之曲柄滑塊連桿組合,設滑塊之偏置量為零,且在水平方向移動,試以此機構之曲桿長度及角度,以及連結桿之長度為輸入項,利用matlab寫出一程式計算在不同曲柄角度時,六點瞬心之對應位置。可順便探討六點瞬心與曲柄角間之關係。
關於這題,我想連桿部份就和第九次的作業大致相同,其重點部份在於瞬心的運用。我們知道有4桿,所以依照公式我們能得到有6個瞬心,(N*(N-1)/2),瞬心的定位法在課本5-20都有詳細敘述,為了標示這些點,我寫了個function slider_4,程式碼如下:


function slider_4(R,L,e)
%先假設從水平方向開始
th1=0;
if R>L
%算出th2
th2=asind(L/R);
else
%th2極限值到90度就不能再畫瞬心
th2=90;
end

th=linspace(th1,th2,200 )
%設定角度間隔,角度範圍分成200等份
[d,th3]=slider_solve(th,R,L,e,1)
%呼叫slider_solve 以求出d

x=R*cosd(th),
y=R*sind(th)
%利用th求出B點座標(x,y)

for n=1:200
%建立一個for迴圈

hold on;
%設定座標範圍
line([0,x(n),d(n)],[0,y(n),e]);
%將A B C三點聯成直線
line([d(n)-3,d(n)+3,d(n)+3,d(n)-3,d(n)-3],[e-2,e-2,e+2,e+2,e-2]);
%畫出滑塊
ground(0,0,2,0)
anchor(-2,0,0,0)
%畫底座

%用plot描繪出運動的各點
plot(0,0,'ro');
plot(x(n),y(n),'ro');
plot(d(n),e,'ro');
plot([0,0],[0,e-d(n)*(y(n)-e)/(x(n)-d(n))],'ro:');
plot([x(n),0],[y(n),e-d(n)*(y(n)-e)/(x(n)-d(n))],'ro:');
plot([x(n),d(n)],[y(n),y(n)*d(n)/x(n)],'ro:');
plot([d(n),d(n)],[0,y(n)*d(n)/x(n)],'ro:');

axis equal;
%設定座標範圍
axis ([-10 30 -10 100]);

pause(0.001);
%每個0.5秒停格一次
clf;
%清除動畫前一次畫面

end
在視窗中輸入slider_4(15,8,0)
其可得到下列動畫
動畫

2007年5月18日 星期五

作業九

A:本人5月3日有來上課,謝謝老師。

B:本人於5月10日下午兩點鐘有參加meeting

C:

請將偏置機構作另類分析,分析過程可採你所知的方式(包括講義中所列的方法)。運動中分以曲桿驅動及滑塊驅動的方式,並說明運動的界限或範圍。設此機構之曲桿長Rcm , 連桿Lcm,滑塊之偏置量為10cm等數據作分析。其中,R=10+(學號末二碼),L=R+5 。

ANS

所謂滑塊偏置機構,即是滑塊滑動方向和水平方向距離e(偏置量)。




我們若要分析此機構,我們首先必須求出偏移量對角度限制的影響。上圖中,我們若要讓此偏置機構動作,必須藉由其中的關係式,求出B點與C點的標。而在課本第四章中,我們參考其function[s,theta21,theta22]=slider_limit(R,L,e)的程式,可求出B點座標。但是我們可以將其函數稍微整理一下,用sind(),cosd()代替多出來角度的轉換,可使程式更為簡易。以下是我的function:

function [s,th1,th2]=slider_limit(R,L,e)
th1=asind(e./(R+L));
th2=180+asind((L-e)/R);
s=(R+L).*cosd(th1)-abs(R-L).*cosd(th2);

(其中s是滑塊的衝量(移動量),th1,th2是AB曲桿和水平夾角度的極限值)

在這裡我們必須注意的是th2的運算,課本所給的是R大於L的時候,所以輸出的値不會出現複數,但是老師給的値是L>R,故我們必須修改一下th2=180+asind((L-e)/R)。

爲了求出滑塊(C點)的座標,我們可以課本上的

function [d,theta3]=slider_solve(theta2,R,L,e,mode)
其程式如下:

function [d,theta3]=slider_solve(theta2,R,L,e,mode)
if nargin<5, mode="0;end" theta="theta2*d2g;" cc="(e-R.*sind(theta))./L;">=0,
theta3=asind(cc);
else
theta3=asind(-cc)+pi;
end
d=L.*cosd(theta3)+R.*cosd(theta);


*驅桿驅動

建立一個function slider_analysis(R,L,e)的函數,其程式碼如下:

function slider_analysis(R,L,e)
[s,th1,th2]=slider_limit(R,L,e)
%先呼叫slider_limit,以求出角度極限值
th=linspace(th1,th2,200 )
%設定角度間隔,角度範圍分成200等份
[d,th3]=slider_solve(th,R,L,e,1)
%呼叫slider_solve 以求出d

x=R*cosd(th),
y=R*sind(th)
%利用th求出B點座標(x,y)

for n=1:200
%建立一個for迴圈
axis ([-40 50 -30 40]);
axis equal;
%設定座標範圍
line([0,x(n),d(n)],[0,y(n),e]);
%將A B C三點聯成直線
line([d(n)-3,d(n)+3,d(n)+3,d(n)-3,d(n)-3],[e-2,e-2,e+2,e+2,e-2]);
%畫出滑塊
ground(0,0,2,0)
anchor(-2,0,0,0)
%畫底座
pause(0.05);
%每個0.5秒停格一次
clf;
%清除動畫前一次畫面

end


*滑塊驅動

若是滑塊驅動,則會比區柄驅動多一部份程式。
我建立了slider_analysis2,其程式碼如下:
function slider_analysis2(R,L,e)

[s,th1,th2]=slider_limit(R,L,e)
th=linspace( th1, th2,200 )
[d,th3]=slider_solve(th,R,L,e,1)
x=R*cosd(th),
y=R*sind(th)

for n=1:200
line([0,x(n),d(n)],[0,y(n),e]);
line([d(n)-3,d(n)+3,d(n)+3,d(n)-3,d(n)-3],[e-2,e-2,e+2,e+2,e-2]);
ground(0,0,2,0)
anchor(-2,0,0,0)
pause(0.01);
clf;
axis ([-100 100 -100 100]);
end
%多加的部份
th=linspace( th2,180-th1,200 )
[d,theta3]=slider_solve(th,R,L,e,-1)
%若mode輸入為-1,則輸出角度會和mode值為1時互補,這是滑塊驅動必須考慮的部份。
%這次角度的範圍要從th2到180-th1

x=R*cosd(th),
y=R*sind(th)

for n=1:200


line([0,x(n),d(n)],[0,y(n),e]);
line([d(n)-4,d(n)+4,d(n)+4,d(n)-4,d(n)-4],[e-3,e-3,e+3,e+3,e-3]);
ground(0,0,2,0)
anchor(-2,0,0,0)
pause(0.01);
clf;
axis ([-100 100 -100 100]);
end


在螢幕上輸入
slider_analysis(20,25,10)

slider_analysis2(20,25,10)

即可得到曲桿驅動以及滑塊驅動的圖形

曲桿驅動


滑塊驅動

2007年5月11日 星期五

作業八

A:本人4月26日曾來上課,謝謝老師。
B:
有一組四連桿,其桿長分別為r=[4 3 3 5],由桿2驅動,設第一固定桿角度theta1=0度; 角速度 td2=10rad/s; 角加速度tdd2=0 rad/s^2。

1. 設桿2角度theta2=45度時,求各點之位置、速度與加速度為何?

我們可以參考老師網頁的第六章,使用function f4bar,其中,function的輸入原則為f4bar(r,theta1,theta2,td2,tdd2,mode,linkdrive) r是各桿長度,theta1:第一桿之水平角,theta2:驅動桿之水平夾角,td2:驅動桿之角速度(rad/sec),tdd2:驅動桿之角加速度(rad/sec^2),mode:+1 or -1. 組合模數,負值表示閉合型,正值為分支型,linkdrive :0 (驅動桿為第二桿); 1 (驅動桿為第三桿),所以依照題目的意思,我們輸入
[val,form]=f4bar([4 3 3 5],0,45,10,0,-1,0)
得到
val =
1.0e+003 *
0.0040 0 0 0 0.0212 + 0.0212i 0.0021 + 0.0021i
0.0021 + 0.0021i 0.0450 0.0100 0 0.0041 - 0.0245i 0.0032 + 0.0049i
0.0011 + 0.0028i 0.0695 -0.0163 0.4914 -0.2121 - 0.2121i 0
-0.0008 + 0.0049i 0.0995 -0.0050 0.3836 -1.8712 - 0.4391i 0
form =
1
此處form=1表示可連成四連桿

*各點位置:
O(0,0)
P(2.12,2.12)
Q(2.12+1.05,2.12+2.81)=(3.17,4.93)
(即為value中第一欄第三行的向量值加上P點座標)
R(4,0)

*各桿跟水平線的夾角
在螢幕上輸入
>> abs(val(:,2))'
ans =

0 45.0000 69.4856 99.5246
即為各桿跟水平線的夾角
*各桿的角速度
>> abs(val(:,3))'

ans =

0 10.0000 16.2681 4.9677
即為各桿的角速度
*各桿的角加速度
>> abs(val(:,4))'

ans =

0 0 491.4428 383.6120
即為各桿的角加速度

2. 繪出此四連桿之相關位置及標明各點之速度方向及大小(以程式為之)。

首先,我們可以在老師第六張講義中找到function drawlinks,存取後在螢幕上輸入drawlinks([4 3 3 5],0,45,-1,0),及可得到下圖。



3. 當桿2迴轉時,求出此組四連桿之限制角度,並繪出其位置(以程式為之)。

我們可以清楚地從老師第六張講義了解function drawlimits在角度限制上的運用。另外,這程式亦結合了fb_angle_limits,因此我們必須先存取此檔案。
在螢附上輸入:
drawlimits([4 3 3 5],0,1,0)
可得
Qstart =

28.9550


Qstop =

331.0450

此類型角度限制在右邊,如圖所示。



4. 設theta2=[0:20:360],試繪出此組四連桿之重疊影像,解釋為何有些沒有值。
依照題意,我們輸入
clf;
for i=0:20:360
drawlinks([4 3 3 5],0,i,1,0);
end
可以得到一個圖,但是會出現:
Combination of links fail at degrees 0.0
Combination of links fail at degrees 20.0
Combination of links fail at degrees 340.0
Combination of links fail at degrees 360.0

這是因為我們剛剛已經藉由function drawlimits得知此四連桿由第二桿驅動時的角度限制(28.95度~331.04度),故在此範圍外的角度沒有值,也就是說,當旋轉到這些角度時,無法構成四連桿。

5. 若將問題三考慮在內,只在可迴轉的範圍內迴轉,請問你能讓此組四連桿作成動畫方式迴轉嗎?
我們首先必須找出角度的範圍,所以設一個迴圈從29~331每隔五度跑一次,其次是要注意標示pause的時間,為了讓畫面美觀,更可以設定座標範圍。
AXIS([-5 5 -5 5])

for i=29:5:331
drawlinks([4 3 3 5],0,i,1,0);
pause(0.2)
clf;

end
或者使用move_4paths的function,並輸入move_4paths([4 3 3 5 ],2,-30,3,0,10,0,1,0,4,100)及可得到動畫及其軌跡的移動。
動畫

2007年4月29日 星期日

作業六

第一部分
某一平面組合機構如下圖,其中包括兩滑塊元件一與地固定,另一分於固定於兩桿。青色者則為滑槽。試

6.1.1 標出桿號及結數,並計算共計有多少連桿及結數。
如圖所示



6.1.2利用古魯伯公式,計算此機構之可動度,請列出其計算方法。
由公式得知
M=3*(N-J-1)+FN=12 J=15 
F為12個旋轉結+1個滑動結+2個滑槽結F=12*1+1*1+2*2=17 M=-12+17=5 
所以自由度為5

6.1.3請利用function[df]=gruebler()函數計算其對應之可動度。

使用function gruebler我們可以在裡面輸入旋轉結12、滑動結1、滑槽結2可得gruebler(12,[12 1 2])自由度為5

6.1.4討論此機構中滑塊及滑槽對可動度之影響

滑槽的自由度為2,可旋轉並可延滑槽方向移動。
第2桿處的滑塊,自由度為2,可旋轉並可延x方向移動F處滑塊,自由度為2,可旋轉並可延滑塊方向移動








第二部份
下面為一個立體機構,分別由兩個旋轉結,一個筒結及三個球結組成。試說明:
· 利用古魯伯公式如何計算整個機構之自由度,可以動嗎?
· 請利用function[df]=gruebler()函數計算其對應之可動度,並相互印證。
· 這裡有所謂楕性自由度嗎?其對整個機構之影響如何?

6.2.1各結之自由度如何?
如圖


1.4為旋轉結,自由度1,為能滑動的圓柱體。
6為圓柱結,自由度2 ,為不能滑動的圓柱體。
2,3,5為球結,自由度3 。

6.2.2利用古魯伯公式如何計算整個機構之自由度,可以動嗎?

M=6(N-J-1)+f=6(6-6-1)+13
M=7 
可算出其自由度為7在減掉惰性自由度可得
M=7-2=5


6.2.3請利用function[df]=gruebler()函數計算其對應之可動度,並相互印證。

gruebler(6,[2 0 0 3 1])
可得自由度為7
減掉惰性自由度可得M=5

6.2.4這裡有所謂楕性自由度嗎?其對整個機構之影響如何?
M=6(7-7-1)+Σf=8
故機構可動
以gruebler函數算之,
gruebler(7,[3 0 0 3 1])
自由度為8
惰性自由度為4,
所以整體自由度為8-4=4。

第三部份
假設有三組四連桿,設第一桿為固定桿,各桿長度分別如下:
1. 第一組:桿1-桿4分別為7,4,6,5cm
2. 第二組:桿1-桿4分別為8,3.6,5.1,4.1cm
3. 第三組:桿1-桿4分別為5.4,3.1,6.6,4.7cm

6.3.1何謂葛拉索機構及非葛拉索機構?

在一四連桿當中,
令四桿的桿長為g(最長桿的長度)
s(最短桿的長度)
p,q(中間長度之兩桿的長度)
則當s+g<p+q時,為葛拉索第一類型至少有一桿可為旋轉桿。
s+g>p+q,為葛拉索第二類型為三搖桿運動。
s+g=p+q,為葛拉索第三類型連桿有兩種形式可繞一圈旋轉。
(一三型為葛索拉機構,第二行為非葛索拉)
6.3.2試問各組應屬何種機構?其迴轉情況會如何?
一:7+4=6+5,故屬於葛拉索第三類型桿,即中立連桿組二:8+3.6>5.1+4.1,故屬於葛拉索第二類桿,即非葛拉索連桿三:6.6+3.1<5.4+4.7,其屬於葛拉索第一類桿



6.3.3試用grashof()函數檢驗上述三組的連桿組合。


一:grashof(1,[7 4 6 5])
為Neutral Linkage
二:grashof(1,[8 3.6 5.1 4.1])
為Non-Grashof Linkage
三:grashof(1,[5.4 3.1 6.6 4.7])
為Crank-Rocker Linkage

6.3.4上述三組連桿若要成為葛拉索機構,則應如何改善?

從上述的分析,我們知道一跟三是葛拉索機構,第二組則為非葛拉索型,若要將其改成葛拉索型機構,只要將葛拉索機構中最長與最短之和小於另外兩桿之和即可變成一個葛拉索機構。