2007年6月15日 星期五

第十三次機動學作業

B94611034 張延瑋


1.試設計一組複式齒輪,使其轉速比為125(請說明思考步驟及結果)。
(1)決定組合數:
每組之轉速比以維持在10以內為佳,超過此值時則需考慮增加齒輪組數。
因此轉速比為125時必須用較多的組數才成達到目標。
首先,將125開方,其值為11.18034,仍然比10大,故使用兩組組合仍嫌不足
,如將125開立方,其值為5,低於10甚多,故三組以上齒輪的組合都可用。

(2)若能直接找到轉速比為5之齒輪組合時,應是最好的選擇,因為只
要將這種組合串聯三組即可。設驅動之小齒輪數最小為10齒,則
依序可以得到對應大齒輪之齒數:5*10=50

(3)大齒輪齒數必須為整數。若此複式齒輪作成之變速齒輪箱不作為計時之用,而僅
作一般吊車的變速系統,則總轉速比與理想比值差異是可接受。由於大部份製造機
械均需要正確的轉速比,以免產生脫序運動
。這時硬要三組之轉速比均相等並不一
定適當,有些時候,可以就總轉速比值之性質求得不同組合之轉速比。
可設計如下
齒數順序為
10:50 20:4016:408:4020:50 等五組。



2.請指出本學期中你自己最感得意的一次作業(請說明其原因,且該
作業必須在自己的部落格內)。
作業五,因為這次為了做圖,花了很多時間在嚐試點的座標該如何安排,
程式碼一大串,光要找錯誤修改,也不知看了多少行的程式碼才找到。
同時每一題也都有做出圖來,整體而言,是很用心也花很多時間製作的一次。

2007年6月14日 星期四

齒輪齒形曲線

*齒輪齒形曲線有下列兩種
1.漸開線齒輪
漸開線的形成這裡有圖示
係將繞於基圓外之直線展開,其端點所的弧線,即為漸開線。

漸開線齒輪的特性
(1) 基圓愈小,漸開線曲率就愈大,基圓愈大,漸開線愈
平直。
(2) 漸開線上任何一點作切線,再作法線,法線必定與基
圓相切。
(3) 漸開線齒輪之周節與壓力角相同即可嚙合。
(4) 兩齒輪之轉速比與基圓直徑成反比。


2.擺線齒輪
擺線齒輪形成這裡有圖示
一滾圓在節圓之外方滾動,滾圓上任意點之軌跡,即為外擺線常作為齒面屈曲線。
另一滾輪圓在節圓之內部滾動,滾圓上任意點之軌跡,即為內擺線,此內擺線
常作為齒腹曲線。

擺線齒輪之特性
(1) 擺線齒輪相互作用的齒面與齒腹,必須由同一滾圓所
產生的內擺線與外擺線。
(2) 嚙合齒輪組的軸心距離要保持一定。
(3) 產生同一齒輪的齒腹與齒面,所用之滾圓可以不同。
(4) 接觸路線永遠在母圓上。

2007年6月13日 星期三

直傘齒輪齒形

斜齒輪是傳遞相交軸與交錯軸傳動的主要零件,在傳遞交軸傳動的零件中,
它能以較高的效率、更平穩地傳遞轉矩。它不但能傳遞兩軸線相交成直角的轉動
,也能傳遞兩軸線相交成鈍角或銳角的轉動。

直齒傘齒輪是以最單純的形狀所形成的斜齒輪,若將它的齒轍向內側延長時
,就會具有軸心上會於一點筆直的齒形,這種直齒傘齒輪,一般都用於較低速之
下的的操作情形來使用,無論在汽車工業或其他產業上對直齒傘齒輪的需求是相
當大的,因直齒傘齒輪的優點在於安裝時儘篦有些許的誤差或由於負荷而產生了
變位等,也不會有應力集中到齒之端部的情形而對於破壞也是安全的。且直齒傘
齒輪的齒形非常適用於高精度鍛造齒輪的設計製造上。

傘齒輪的齒形工方法很多,必須根據零件的結構、材料、熱處理、齒輪精度
及生產批量等一系列因素,選擇最適合的加工方法。一般工業上傘齒輪的切削加
工方法有三大類:
(1)滾切法(generate),圖2為Gleason NO.710機台以雙刀盤銑齒法製成直傘齒。
(2)成形法(format),圖3 為Gleason NO.710機台就是以成形法加工直傘齒。
(3)仿形法
圖2Gleason NO.104機台 圖3Gleason NO.710機台 齒輪的主要切削製造程序為原料->鍛胚->粗切削->熱處理->精切削->研磨。
其程序複雜且費時,無論在粗切削、精切削、研磨時以成形法來製造必須一齒一
齒往複式銑削,齒數越多越費時,對於精度的要求高,成本當然高。鍛造為以高壓
一次成形,在時間上就相當節省;在材料的強度上,鍛造品在抗彎強度和衝擊強度
的實驗證明下都比切削品高,如以精密鍛造齒輪的技術來製成,可在製成上大量的
節省成本及時間,提升產業的競爭力。
圖4鍛造直齒傘齒輪應用於沙灘車 針對直傘齒輪之設計,本實驗室開發一套直傘齒輪設計軟體,可以模擬Gleason 公司的直傘齒創成機台,輸出齒輪製造數據及量測程式如圖5,並輸出3D模型 如圖6以供製造齒輪。 <-圖5鍛造直齒傘齒輪應用於沙灘車 <-圖6直傘齒輪3D模型

如何測定齒厚呢?

在課本上有提到齒厚之計算用以分析齒輪軸承受應力及其應變
書上內容是利用漸開線函數 或用MATLAB函數之ainv()來求得齒厚
關於齒厚測定之種類 是可分成三種1.弦齒厚2.跨齒厚3.跨銷值
1.弦齒厚(節圓弦齒厚) (量測弦齒厚應在節圓處)
圖示
在量測時是以齒頂為基準,量測左右兩側齒面與節圓相交處的齒厚。
關於將齒厚游標卡尺的用法及順序說明於下:
(1)先把本尺爪對準滑塊爪,使齒厚尺的 0 點確定。此時因卡尺兩爪部的
前端容易磨損,必須小心。

(2)把量齒高的榫舌對合卡尺的齒高,用進給螺絲固定。然後用進給輪來
移動齒高滑塊,使它對準游標卡尺,再用固定螺絲把滑塊固定。

(3)把榫舌的測定面接觸到齒輪的齒冠圓(外周),推動進給輪來進給齒厚
滑塊,使滑塊爪的測定面與齒面接觸時,由游標卡尺的刻度讀出其尺寸 。
齒厚滑塊的進給法與高滑塊相同。

2.跨齒厚 (跨齒數的選擇最好能使量測的接觸點置於節圓處)
圖示
測定部位在齒輪圓周上各差90度的齒寬中央處。用錶規與
卡尺組合而成的測定器,在小齒輪的測定時非常方便。

3.跨銷值 (量測的直線距離要能垂直齒輪軸線方向)
圖示
齒厚測定方面還有銷徑(Over Pin)法。這是在相對於正齒輪直徑上的齒溝(偶數溝)
中各插入銷,然後測定其外側尺寸(內齒輪時測其內側尺寸),而求出齒厚的方法。
在螺旋齒輪的情況,外齒是量測兩量銷或珠子外側的距離,內齒輪是量測兩個珠
子中間的距離。

*另外齒厚的量測位置在弦齒厚與跨齒厚的量測位置應取齒面寬一半處,每隔大約90度作
一選擇來量測。跨銷值法亦置於齒面寬的一半處,以大約垂直的 兩跨銷處作量測。

2007年6月8日 星期五

機動學作業十二

B94611034 張延瑋
本人5/31日有來上課
一組標準全齒輪齒輪之徑節為8,齒數分別為30T與48T,其工作壓力角為20度或14.5度。

1.試求其接觸線長度,與接觸比。
本題利用MATLAB函式contact_ratio()計算,此函式之輸入值為徑節、兩齒輪之齒數及壓力角,可輸出項為接觸比和接觸長度、齒冠、周節及基周節、兩齒輪節圓直徑、兩齒輪之接近角、遠退角及作用角。
>> [c_ratio,c_length,ad,pc,pb,d2,d3,ag]=contact_ratio(8,30,48,[20 14.5]')

c_ratio =
1.7005 ...壓力角20度時的接觸比
2.0308 ...壓力角14.5度時的接觸比 <--接觸比不小於1.2 表示其各對牙所承受的負荷較輕, 力量傳遞較均勻

c_length =
0.6275 ...壓力角20度時的接觸線長度
0.7721 ...壓力角14.5度時的接觸線長度

ad = 0.1250 ...齒冠

pc = 0.3927 ...周節
pb =
0.3690 ...基周節
0.3802

d2 = 3.7500 ...齒輪二直徑
d3 = 6 ...齒輪三直徑


ag = 接近角 遠退角 作用角
10.4850 9.9211 20.4061 ...齒輪二(壓力角20度時)
12.6898 11.6797 24.3695 ...齒輪二(壓力角14.5度時)
6.5532 6.2007 12.7538 ...齒輪三(壓力角20度時)
7.9311 7.2998 15.2309 ...齒輪三(壓力角14.5度時)


2.兩齒輪之節圓、基圓直徑各為如何?請列式計算其結果。

節圓直徑*徑節=齒數 (徑節為8)
齒輪二的直徑d2=30/8=3.75
齒輪三的直徑d3=48/8=6 與函式contact_ratio()之結果 相符

基圓直徑=節圓直徑*cosθ (θ為工作壓力角)
齒輪二 d2*cos20°=3.523 d2*cos14.5°=3.631
齒輪三 d3*cos20°=5.638 d3*cos14.5°=5.809

3.此組齒輪是否會產生干涉現象?試列式證明之。
根據課本9.47式,N2、N3若滿足以下條件(N2²+2*N2*N3)*sin²(θ) >= 4(1+N3)
則不會發生干涉。

(i)N2為30,N3為48,壓力角為20度時
30(30+2*48)*sin²20°=442.176>4*(1+48)=196
>> isinterf(30,48,20)
ans =0 no interferenc
(ii)N2為30,N3為48,壓力角為14.5度時
30(30+2*48)*sin²14.5°=236.969>4*(1+48)=196
>> isinterf(30,48,14.5)
ans =0 no interferenc
兩種情形皆不會產生干涉

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

使用教學網頁上老師提供的函式draw_gear與函式move2_gear
動畫的部分

當工作壓力角=20度
>>move2_gear(8,30,48,20,10)
壓力角=20°的動畫

當工作壓力角=14.5度
>>move2_gear(8,30,48,14.5,10)
壓力角=14.5°的動畫

本來輸入齒輪二、齒輪三的數據時,兩齒輪並不會契合
為了使兩齒輪能契合,讓齒輪三在執行draw_gear.m
coord1=rotate2D(coord,-i+3.75,x0,y0) 做點修改 多轉3.75度,
若修改齒輪二亦可得到相同的圖。
壓力角20度
>> [coords]=draw_gear(8,30,20,360,0,0);
>> [coords]=draw_gear(8,48,20,360,9.75/2,0);
圖示
壓力角14.5度
>> [coords]=draw_gear(8,30,14.5,360,0,0);
>> [coords]=draw_gear(8,48,14.5,360,9.75/2,0);
圖示

2007年6月1日 星期五

機動學作業十一

作業十一 B94611034 張延瑋

1.我有上本週(5/24)的機動學課程
2.參考老師提供的parabol_cam
編寫出一函式test 可做出
位移,速度,加速度-角度關係圖

程式如下:
function test
% 計算返程對應點之資料:位移,速度及加速度
theta=0:10:100;
for i=1:length(theta)%共11項
ss(i)=0
vv(i)=0
aa(i)=0
end;
%(theta:對應角度; ss:位移, vv:速度, aa:加速度)
theta=0:10:200;
for i=11:[length(theta)] %共11
[ss(i),vv(i),aa(i)]=parabol_cam(theta(i),100,100,1,5,0);
end;
%(theta:對應角度; ss:位移, vv:速度, aa:加速度)
theta=0:10:260;
for i=22:[length(theta)] %共6項
ss(i)=5
vv(i)=0
aa(i)=0
end;
%(theta:對應角度; ss:位移, vv:速度, aa:加速度)
theta=0:10:360;
for i=27:[length(theta)] %共11項
[ss(i), vv(i), aa(i)]=parabol_cam(theta(i),260,100,-1,5,0);
end;
%(theta:對應角度; ss:位移, vv:速度, aa:加速度)
[theta' ss' vv' aa'] %列出來
axis([0 360 -40 20]); %固定座標軸
plot(theta,ss,' b-','LineWidth',3.5) %繪出角度與位移之關係圖(藍色)
hold on;
plot(theta,vv,'r:','LineWidth',2.5) %繪出角度與速度之關係圖(紅色)
hold on;
plot(theta,aa,'g+','LineWidth',1.5) %繪出角度與加速度之關係圖(綠色)
hold on;
axis([0 360 -40 20]);
title('B94611034位移,速度,加速度-角度');

function createlegend(axes1)
legend1 = legend(axes1,{'位移','速度','加速度'},'Position',[0.1312 0.1091 0.1607 0.1444]);

3.
使用程式pincam 可做出如下的圖
其將運動的型式分成5種 可以不同代號輸入二元素之列矩陣
而升程及返程之範圍,以三元素列矩陣表示
凸輪半徑為15公分,順時針方向旋轉,梢型從動件梢型,垂直接觸,長為10公分
參照第二題的運動形式
輸入[x y]=pincam([0:10:360],15,5,0,10,[100 200 260],[2 2],-1)
可得 凸輪之工作曲線

4.
將程式pincam 最後部份
加入rotate
顏色部分以三原色做調配
編改後程式如下
function [x,y]=pincam(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)
figure(1);
clf;
th=cth*pi/180;
for i=1:length(cth)
t=th(i)*cw;
A=[cos(t) -sin(t);sin(t) cos(t)]; %此為A點座標(s0,e)
s0=sqrt(r0*r0-e*e); %從動中心線與偏置圓之切點至A點之距離為s0
[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);
line(X(1,1:2),X(2,1:2)); %由基圓圓心至凸輪外廓邊緣的射線
line(X(1,2:3),X(2,2:3),'linewidth',3,'color',[0.5 0.3 0]) %梢型從動件垂直於工作曲線(棕色)
end
hold on;
ok=plot([0 x],[0 y],'ro',x,y,'color',[0.4 1 0.5]) %將各點連接在一起而成工作曲線(綠色)
axis equal
for n=1:360 %分成360等份旋轉
rotate(ok,[0 0 1],-1,[0 0 0]) %ok為物件工作曲線
pause(0.0034)
end;
如此一來 就能跑出凸輪迴轉動畫

2007年5月25日 星期五

機動學作業十




1.我本週機動學課程(5/17)有來上課
2.速度與加速度的問題,當一桿以某特定點M等角速度迴轉時,其端點P之速度方向如何?
其加速度方向如何?若該特定點M復以等速水平運動,則同一端點P之速度與加速度方向會變為如
何?若M點同時也有加速度,則點P會有何變化?若以此推理四連桿的運動,則點P與Q之速度與
加速度方向會與桿一(固定桿)之兩端點之關係如何?與我們前面的作業分析結果有無共通之處?
(參看第六章之四連桿機構之運動分析)
(1)若M,P兩點相距L(m) 則M以等角速度w rad/s 旋轉

使得P點有速度L*w(m/s),其方向與L為垂直方向

其向心加速度為L*(w^2) (m/s^2)

(2)M復以等速V(m/s)水平運動

P點速度為 V+L*w*sin(θ) i - L*w*cos(θ) j (m/s)

向心加速度仍為L*(w^2) (m/s^2)

(3)M點同時也有加速度a (m/s^2)時,

p點保持原速,加速度為V+L*w^2*cos(θ) i+L*w^2*sin(θ) j (m/s)

以此推理四連桿的運動,可知ORPQ形成一四連桿 如下圖
四連桿

P點繞著O點 以PO長為半徑轉動

Q點繞著R點 以QR長為半徑旋轉
速度分析
VQ=VP+VQ/P
VP=VQ+VP/Q

加速度分析
由P點觀察
AQ = AP+ tAQ/P+ rAQ/P
rAQ/P = w32 *RQ/P
其中tAQ/P 項之角加速度α3是未知,但其加速度方向與RQ/P 垂直

由D點觀察
AQ = AR+ tAQ/R+ rAQ/R
rAQ/R =w4*(w4*RQ/P)
tAQ/R項因為角加速度α4仍為未知,但其方向與RQ/R垂直
由上兩個方向可決定Q點的加速度AQ
同理若求P ,則可從O和Q觀察而得

3.設有一運動之曲柄滑塊連桿組合,設滑塊之偏置量為零,在水平方向移動
可寫出一動畫 桿長輸入20cm可得 動畫

程式如下
function sliderdrew4limt2(R,L,e)
clf;
grid on
hold on %保留住圖
a=cosd(linspace(R+L,90,100)); %連桿的x座標
b=sind(linspace(R+L,90,100)); %連桿的y座標
c=R*a; %桿上c點的位置
d=L*b; %桿上d點的位置
data=slider_solved(linspace(R+L,90,100),R,L,e,1); % 使用slider_sovled 資料以data表示

for n=1:99
hold on
axis ([-10 50 -10 60]); %座標範圍
line([50,-50],[-0,0],'color','r') %水平線以紅色表示
plot(c(n),d(n),'mo:') %粉紅色的連線
plot(data(n),e,'mo:') %將瞬心點以圓圈標示
g=[0,0] %A(和曲柄相接的原點)
k=[c(n),d(n)] %B點
j=[data(n),e] %C點(與滑塊連接的)
linkshape(g,k,2) %曲柄桿R linkshape以綠色呈現
linkshape(k,j,2) %聯結桿L
plot([0,0],[0,e-data(n)*(d(n)-e)/(c(n)-data(n))],'mo:') %從原點為固定點 瞬心點到連桿之間連線
plot([data(n),data(n)],[0,d(n)*data(n)/c(n)],'mo:') %將瞬心點以圓圈標示
plot([c(n),data(n)],[d(n),d(n)*data(n)/c(n)],'mo:') %將瞬心軌跡以圓圈標示
line([data(n)-3,data(n)+3,data(n)+3,data(n)-3,data(n)-3],[e-2,e-2,e+2,e+2,e-2]); %繪製滑塊的部分
plot([c(n),0],[d(n),e-data(n)*(d(n)-e)/(c(n)-data(n))],'mo:')
axis ([-20 50 -20 50]);
pause(0.12)
clf
end

若以曲桿長度及角度,以及連結桿之長度為輸入項 只需將for的部分作些微修改即可
圖示


六點瞬心與曲柄角間之關係
曲柄滑塊組合之瞬心
由圖可知此滑塊機構的基本瞬心有12 23 34 14 等四點

點14因為滑塊的關係 其方向應與地平線垂直 即點14存在於無窮遠的位置
與其他點配合 只要過該點作平行線於此線的線 即為所求兩點間的連線
利用瞬心點13之位置可求得桿2端的速度及桿2的迴轉速度
就桿2來看 點12和點23為其瞬心 前者固定成轉動軸 端點23速度為v
則桿2轉速w=v/r2

2007年5月20日 星期日

作業九

我有上本週5/3的機動學課程

若以第二桿驅動時,發現四種情況皆無法成立
因此無法移動

因為機構連結桿與曲桿之長度L-R=10+34-(44+5)差小於偏置量z
所以此曲桿無法進行完整的360度一圈的旋轉 僅於某一範圍活動
R=34+10=44
L=44+5=49
當以第三桿驅動時,符合第一類情形
0≦r4,且0≦r2-r4≦r3
所以帶入程式中
>> drawsldlimits([1 44 49 10],0,-1,1)

Qstart =-43.9378

Qstop =223.9378
得到第三桿所能運動的最小及最大角度,
範圍在 -43.9378 < θ < 223.9378
利用以下程式
for th=-43.9378:( Qstop- Qstart)/20: 223.9378

%(Qstop- Qstart)為總角度的變化量
pause(0.2);
drawsldlinks([1 44 49 10],0,th,-1,1) ;
axis([-110 80 -50 50]);
end;
做出動畫如下 動畫1

當以滑塊為驅動時,發現符合第二類的情形
|r2-r3|< r4

由程式drawsldlimits([1 44 49 10],0,-1,2)

Qstart =-92.4608

Qstop =92.4608

結果r1的範圍在 -92.4608到92.4608之間
根據公式
θmin=θ1 + sin-1(r4/( r2+r3))
θmax=θ1 +π- sin-1(r4/ (r2+r3))
在這範圍間第二桿角度的變化為
θmin=6.172
θmax=173.827
再利用老師提供的程式
move_sldpaths([10,44,49,10],0,0,3,0,10,0,1,0,4,100)
(r,r6,th6,nlink,th1,td2,tdd2,sigma,driver,ntimes,npts)
輸入參數:
r : 四連桿之長度,th1:第一桿的水平角度(通常為零)。
th2: 驅動桿或聯結桿的水平角度。
td2, tdd2:驅動桿之角速度及角加速度。
Sigma:連桿之組合模式選擇(需為+1或-1)
Driver:驅動桿之設定( 0表示f曲桿驅動,1聯結桿驅動,3滑塊驅動)。
Npts: 設定分割的點數或位置。
r6,rh6,nlink:桿上特定點之位置。
做出動畫如下 動畫2

2007年5月11日 星期五

作業八

B94611034 張延瑋
A:我(4/26)曾來上課。


B:
問題一:設桿2角度theta2=45度時,求各點之位置、速度與加速度為何?






虛線為加速度  


綠線為速度





第2、3節點之速度為   30/s 
第2、3節點之加速度為  300/s^2
第3、4節點之速度為   24.8387/s
第3、4節點之加速度為  1922/s^2







問題二:繪出此四連桿之相關位置及標明各點之速度方向及大小(以程式為之)。
能此函數能會製出四連桿和其加速度
function f4barline(r,theta1,theta2,td2,tdd2,mode,linkdrive)
if nargin<7,linkdrive=0;end mode="1;end" data="zeros(4,6);" linkdrive="=" r="[r(1)" rr="r.*r;d2g=" t1="theta(1);tx=" s1="sin(t1);c1=" sx="sin(tx);cx=" a="2*r(1)*r(4)*c1-2*r(2)*r(4)*cx;" c="rr(1)+rr(2)+rr(4)-rr(3)-2*r(1)*r(2)*(c1*cx+s1*sx);" b="2*r(1)*r(4)*s1-2*r(2)*r(4)*sx;" pos="B*B-C*C+A*A;">=0,
form=1;
if abs(C-A)>=1e-5
t4=2*atan((-B+mode*sqrt(pos))/(C-A));
s4=sin(t4);c4=cos(t4);
t3=atan2((r(1)*s1+r(4)*s4-r(2)*sx),(r(1)*c1+r(4)*c4-r(2)*cx));
s3=sin(t3);c3=cos(t3);
else
A=-2*r(1)*r(3)*c1+2*r(2)*r(3)*cx;
B=-2*r(1)*r(3)*s1+2*r(2)*r(3)*sx;
C=rr(1)+rr(2)+rr(3)-rr(4)-2*r(1)*r(2)*(c1*cx+s1*sx);
pos=B*B-C*C+A*A;
if pos>=0,
t3=2*atan((-B-mode*sqrt(pos))/(C-A));
s3=sin(t3); c3=cos(t3);
t4=atan2((-r(1)*s1+r(3)*s3+r(2)*sx),...
(-r(1)*c1+r(3)*c3+r(2)*cx));
s4=sin(t4);c4=cos(t4);
end
end
theta(3)=t3;theta(4)=t4;
td(2)=td2;
AM=[-r(3)*s3, r(4)*s4; -r(3)*c3, r(4)*c4];
BM=[r(2)*td(2)*sx;r(2)*td(2)*cx];
CM=AM\BM;
td(3)=CM(1);td(4)=CM(2);
tdd(2)=tdd2;
BM=[r(2)*tdd(2)*sx+r(2)*td(2)*td(2)*cx+r(3)*td(3)*td(3)*c3-...
r(4)*td(4)*td(4)*c4;r(2)*tdd(2)*cx-r(2)*td(2)*td(2)*sx-...
r(3)*td(3)*td(3)*s3+r(4)*td(4)*td(4)*s4];
CM=AM\BM;
tdd(3)=CM(1);tdd(4)=CM(2);
if linkdrive==1,
c2=c3;c3=cx;s2=s3;s3=sx;
r(2:3)=[r(3) r(2)];theta(2:3)=[theta(3) theta(2)];
td(2:3)=[td(3) td(2)];tdd(2:3)=[tdd(3) tdd(2)];
else
c2=cx;s2=sx;
end
for j=1:4,
data(j,1:4)=[r(j)*exp(i*theta(j)) theta(j)/d2g td(j) tdd(j)] ;
end % position vectors
data(1,5)=r(2)*td(2)*exp(i*theta(2));%velocity for point Q
data(2,5)=r(4)*td(4)*exp(i*theta(4));%velocity for point P
data(3,5)=r(2)*(i*tdd(2)-td(2)*td(2))*exp(i*theta(2));%acc of Q
data(4,5)=r(4)*(i*tdd(4)-td(4)*td(4))*exp(i*theta(4));%acc of P
data(1,6)=data(2,1);%position of Q, again
data(2,6)=data(1,1)+data(4,1);% position of P
rotate=abs(data(:,2));
fun=[0 r(2)*cosd(rotate(2)) r(4)*cosd(rotate(4)+r(1)) r(1) ,0;0 r(2)*sind(rotate(2)) r(4)*sind(rotate(4)+0) 0,0]'
clf
line(fun(:,1),fun(:,2),'color','r')
triangle1 = [[fun(2,1),fun(2,1)+real(data(1,5))/10]' [fun(2,2),fun(2,2)+imag(data(1,5))/10]']
line(triangle1(:,1),triangle1(:,2)','color','g')
triangle2 = [[fun(2,1),fun(2,1)+real(data(3,5))/100]' [fun(2,2),fun(2,2)+imag(data(3,5))/100]']
line(triangle2(:,1),triangle2(:,2)','color','b','linestyle',':')
triangle3 = [[fun(3,1),fun(3,1)+real(data(2,5))/10]' [fun(3,2),fun(3,2)+imag(data(2,5))/10]']
line(triangle3(:,1),triangle3(:,2)','color','g')
triangle4=[[fun(3,1),fun(3,1)+real(data(4,5))/100]',[fun(3,2),fun(3,2)+imag(data(4,5))/100]']
line(triangle4(:,1),triangle4(:,2),'color','b','linestyle',':')
else
form=0;
if linkdrive==1,
r=[r(1) r(3) r(2) r(4)];
for j=1:4, data(j,1)=r(j).*exp(i*theta(j));end % positions
end
end




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

使用網頁上之move_4paths函數 求得其限制角度為29度及331度















問題四:設theta2=[0:20:360],試繪出此組四連桿之重疊影像,解釋為何有些沒有值。

其為雙搖桿型連桿 任何桿都不可能產生完整的轉動 因此題目中的主動桿必存在死點且
位在第二桿為29度及331度時















問題五:若將問題三考慮在內 ,只在可迴轉的範圍內迴轉,請問你能讓此組四連桿作成動畫方式迴轉嗎?
可以
動畫1
動畫2

2007年4月26日 星期四

作業七

B94611034 張延瑋
延瑋的部落格

本人(4/19)有上課

rho=[14, 19, 9] a=14
theta=[0, 0, 0]
td=[0.2, 0.5, 0.3]
tdd=[0, 0.1, 0.2]

1.
t=1時端桿之對應方位

t=2時端桿之對應方位

t=3時端桿之對應方位

t=4時端桿之對應方位

t=5時端桿之對應方位
2.
三桿之結點之速度,加速度與時間t之關係
桿一之加速度與時間
桿一之加速度與時間
桿二之加速度與時間
桿二速度與時間
桿三速度與時間
桿三速度與時間
3.
t=[1 2 3 4 5]的動畫
應用老師網路上所提供的程式dyade draw 製作出與時間相關的函式 如下
function dyad_use(rho,theta,td,tdd,time)
for n=1:length(rho)
for m=1:length(rho)
a(m)=td(m)+tdd(m)*(n-1);
b(m)=theta(m)+td(m)*(n-1)+0.5*tdd(m)*(n-1)^2;
pause(0.1)
end
[A B]=dyad_draw(rho,b,a,tdd) ;
end
if n > 1
for k=1:length(rho)
k2=n-1
figure1(k,k2)=A(k);
figure2(k,k2)=B(k);
figure
plot(figure1(k,:))
figure
plot(figure2(k,:))
end
end

機動學作業六

B94611034 張延瑋
我有上本週(十二日)的課。

6.1.1
如圖所標之桿序 (P處之滑塊亦算一桿)
平面組合機構之可動圖

總桿數 N=12
M處為(3-1)=2結
N處為(3-1)=2結
Q處為(3-1)=2結
P處桿12與1間有一滑塊 連結度為2 有1個結
R有一個結 連結度為2
S處為(3-1)=2結
T處為槽中梢 2個結 連結度為3

總結數 J=2+1(P處)+1(R處)+2(M處)+2(N處)+2(Q處)+2(S處)+1(滑塊)+2(槽中梢)=15
M=3(N-J-1)+f
N=12 J=15

6.1.2
套入公式計算 4桿的f值為2 10桿的f值為2
f為12個旋轉結+1個滑動結+2個滑槽結
f=12*1+1*1+2*2=17
M=3(12-15-1)+17=5 df=5
故此機構的可動度應為5

6.1.3
>>df=gruebler(12,[12 1 2])
df=5

6.1.4
P處的滑塊會滑動 也將其視為一桿 也多一滑動的自由度 T槽中梢提供滑動與轉動的自由度=2
R滑槽提供滑動與轉動的自由度=2


6.2.1
如圖所示立體組合機構之可動圖

球結的自由度為3
旋轉結的自由度為1
圓柱結的自由度為2

6.2.2
M=6(N-J-1)+f=6(6-6-1)+13
 M=7 其自由度為7

6.2.3
>>df=gruebler(6,[2 0 0 3 1])
df=7

6.2.4
公式中所計算之自由度有時與實際之運作仍有出入
因為球結的存在 有些連桿會有自轉運動 經觀察得知
連桿2 連桿6可自轉 故具有2個惰性自由度
其實際之可動度減為5 比古魯伯公式所求得結果少
惰性自由度主要為桿的自轉所引起 因為當其自轉時不會牽動到機構中其他部分的運動
同時也不會影響系統外型 另外若在可自轉的該桿上多安插一桿,則可增加自由度


6.3
在一四連桿組中
g=最長桿之長度
s=最短桿之長度
p q=中間長度桿之長度
(1)當最短桿與最長桿之和小於其他兩桿之和時,則至少有一桿可為旋轉桿。稱為葛拉索型機構
s+g <> p+q

6.3.1
第一組:桿1-桿4分別為7,4,6,5cm
7+4=6+5,屬於葛拉索第三類桿,即是中立連桿組
行四邊形機構和平行四邊形反置機構,

第二組:桿1-桿4分別為8,3.6,5.1,4.1cm
8+3.6>5.1+4.1,屬於葛拉索第二類桿,
即是非葛拉索連桿 Crank-Rocker Linkage機構

第三組:桿1-桿4分別為5.4,3.1,6.6,4.7cm
6.6+3.1<5.4+4.7,屬於葛拉索型>> ans=grashof(1,[7 4 6 5])
ans=Neutral Linkage

第二組:
>> ans=grashof(1,[8 3.6 5.1 4.1])
ans=Non-Grashof Linkage


第三組:
>> ans=grashof(1,[5.4 3.1 6.6 4.7])
ans=Crank-Rocker Linkage

6.3.3
上述三組連桿 只有第二組為非葛拉索型 若要改成為葛拉索機構
把最長桿或最短桿減短,
或增加第二和第三長的連桿長度,
進而使之滿足葛拉索機構最長與最短之和小於另外兩桿之和。

2007年4月12日 星期四

機動學 作業五

5.1







5.1.1.





上臂的圖形 長度L1
arm1=[0 -L1/5;L1/5*cosd(-80) L1/5*sind(-80);L1/5*cosd(-70) L1/5*sind(-70);L1/5*cosd(-60) L1/5*sind(-60);L1/5*cosd(-50) L1/5*sind(-50);L1/5*cosd(-40) L1/5*sind(-40);L1/4*cosd(-45) L1+L1/4*sind(-45);L1/4*cosd(-30) L1+L1/4*sind(-30);L1/4*cosd(-10) L1+L1/4*sind(-10);L1/4*cosd(10) L1+L1/4*sind(10);L1/4*cosd(30) L1+L1/4*sind(30);L1/4*cosd(45) L1+L1/4*sind(45);L1/4*cosd(60) L1+L1/4*sind(60);L1/4*cosd(75) L1+L1/4*sind(75);L1/4*cosd(90) L1+L1/4*sind(90);L1/4*cosd(105) L1+L1/4*sind(105); L1/4*cosd(120) L1+L1/4*sind(120);L1/4*cosd(135) L1+L1/4*sind(135);L1/4*cosd(150) L1+L1/4*sind(150);L1/4*cosd(170) L1+L1/4*sind(170);L1/4*cosd(190) L1+L1/4*sind(190);L1/4*cosd(210) L1+L1/4*sind(210);L1/4*cosd(225) L1+L1/4*sind(225);L1/5*cosd(220) L1/5*sind(220);L1/5*cosd(230) L1/5*sind(230);L1/5*cosd(240) L1/5*sind(240) ;L1/5*cosd(250) L1/5*sind(250);L1/5*cosd(260) L1/5*sind(260);L1/5*cosd(270) L1/5*sind(270)]

下手臂的圖形 長度L2
arm2=[L2/5*cosd(80) L2/5*sind(80);L2/5*cosd(90) L2/5*sind(90);L2/5*cosd(100) L2/5*sind(100);L2/5*cosd(120) L2/5*sind(120);L2/5*cosd(140) L2/5*sind(140);L2/5*cosd(160) L2/5*sind(160);L2/5*cosd(180) L2/5*sind(180);L2/5*cosd(-160) L2/5*sind(-160);L2/5*cosd(-140) L2/5*sind(-140);L2/5*cosd(-120) L2/5*sind(-120);L2/5*cosd(-100) L2/5*sind(-100);L2/5*cosd(-90) L2/5*sind(-90);L2/5*cosd(-80) L2/5*sind(-80);L2+L2/7*cosd(-100) L2/7*sind(-100);L2+L2/7*cosd(-90) L2/7*sind(-90);L2+L2/7*cosd(-70) L2/7*sind(-70);L2+L2/7*cosd(-50) L2/7*sind(-50);L2+L2/7*cosd(-30) L2/7*sind(-30);L2+L2/7*cosd(-10) L2/7*sind(-10);L2+L2/7*cosd(10) L2/7*sind(10);L2+L2/7*cosd(30) L2/7*sind(30);L2+L2/7*cosd(50) L2/7*sind(50);L2+L2/7*cosd(70) L2/7*sind(70);L2+L2/7*cosd(90) L2/7*sind(90);L2+L2/7*cosd(100) L2/7*sind(100);L2/5*cosd(80) L2/5*sind(80)]

手掌的圖形 長度L3
palm=[L2+L3/4*cosd(-90) L3/4*sind(-90);L2+L3/4*cosd(-100) L3/4*sind(-100);L2+L3/4*cosd(-120) L3/4*sind(-120);L2+L3/4*cosd(-140) L3/4*sind(-140);L2+L3/4*cosd(-160) L3/4*sind(-160);L2+L3/4*cosd(-180) L3/4*sind(-180);L2+L3/4*cosd(-200) L3/4*sind(-220);L2+L3/4*cosd(-240) L3/4*sind(-240);L2+L3/4*cosd(-260) L3/4*sind(-260);L2+L3/4*cosd(-270) L3/4*sind(-270);L2+L3/4 L3/4;L2+L3/3+L3/3.2*cosd(30) L3/4+L3/3.2*sind(30);L2+L3/3+L3/3.2*cosd(26) L3/4+L3/3.2*sind(26);L2+L3/3+L3/8 L3/4;L2+L3/3+L3/4+L3/2*cosd(11) L3/8+L3/2*sind(11);L2+L3/3+L3/4+L3/2*cosd(9) L3/8+L3/2*sind(9);L2+L3/3+L3/4 L3/8;L2+L3/3+L3/4+L3/2*cosd(1) L3/2*sind(1);L2+L3/3+L3/4+L3/2*cosd(-1) L3/2*sind(-1);L2+L3/3+L3/4 0;L2+L3/3+L3/4+L3/2*cosd(0) -L3/8+L3/2*sind(0);L2+L3/3+L3/4+L3/2*cosd(-2) -L3/8+L3/2*sind(-2);L2+L3/3+L3/4 -L3/8;L2+L3/3+L3/4+L3/3*cosd(-2) -L3/4+L3/3*sind(-2);L2+L3/3+L3/4+L3/3*cosd(-4) -L3/4+L3/3*sind(-4);L2+L3/3+L3/4 -L3/4;L2+L3/4*cosd(-90) L3/4*sind(-90)]

line(arm1(:,1),arm1(:,2),'Color','r')
將arm1的點連起來成一圖 (紅)
line(arm2(:,1),arm2(:,2),'Color','m')
將arm2的點連起來成一圖 (粉紅)
line(palm(:,1),palm(:,2),'Color','b')
將palm的點連起來成一圖(綠)
axis equal

5.1.2.


% X1,X2,X3為輸入角度 
% L1 L2 L3分別為上手臂,下手臂及手掌的長度







function body(L1,L2,L3,X1,X2,X3) %程建立一程式函數body



theta1=-(X1-90)



theta2=-(180-X2+X1)



theta3=-(360-X3-X2+X1)



arm1=[0 -L1/5;L1/5*cosd(-80) L1/5*sind(-80);L1/5*cosd(-70) L1/5*sind(-70);L1/5*cosd(-60) L1/5*sind(-60);L1/5*cosd(-50) L1/5*sind(-50);L1/5*cosd(-40) L1/5*sind(-40);L1/4*cosd(-45) L1+L1/4*sind(-45);L1/4*cosd(-30) L1+L1/4*sind(-30);L1/4*cosd(-10) L1+L1/4*sind(-10);L1/4*cosd(10) L1+L1/4*sind(10);L1/4*cosd(30) L1+L1/4*sind(30);L1/4*cosd(45) L1+L1/4*sind(45);L1/4*cosd(60) L1+L1/4*sind(60);L1/4*cosd(75) L1+L1/4*sind(75);L1/4*cosd(90) L1+L1/4*sind(90);L1/4*cosd(105) L1+L1/4*sind(105); L1/4*cosd(120) L1+L1/4*sind(120);L1/4*cosd(135) L1+L1/4*sind(135);L1/4*cosd(150) L1+L1/4*sind(150);L1/4*cosd(170) L1+L1/4*sind(170);L1/4*cosd(190) L1+L1/4*sind(190);L1/4*cosd(210) L1+L1/4*sind(210);L1/4*cosd(225) L1+L1/4*sind(225);L1/5*cosd(220) L1/5*sind(220);L1/5*cosd(230) L1/5*sind(230);L1/5*cosd(240) L1/5*sind(240) ;L1/5*cosd(250) L1/5*sind(250);L1/5*cosd(260) L1/5*sind(260);L1/5*cosd(270) L1/5*sind(270)]



arm2=[L2/5*cosd(80) L2/5*sind(80);L2/5*cosd(90) L2/5*sind(90);L2/5*cosd(100) L2/5*sind(100);L2/5*cosd(120) L2/5*sind(120);L2/5*cosd(140) L2/5*sind(140);L2/5*cosd(160) L2/5*sind(160);L2/5*cosd(180) L2/5*sind(180);L2/5*cosd(-160) L2/5*sind(-160);L2/5*cosd(-140) L2/5*sind(-140);L2/5*cosd(-120) L2/5*sind(-120);L2/5*cosd(-100) L2/5*sind(-100);L2/5*cosd(-90) L2/5*sind(-90);L2/5*cosd(-80) L2/5*sind(-80);L2+L2/7*cosd(-100) L2/7*sind(-100);L2+L2/7*cosd(-90) L2/7*sind(-90);L2+L2/7*cosd(-70) L2/7*sind(-70);L2+L2/7*cosd(-50) L2/7*sind(-50);L2+L2/7*cosd(-30) L2/7*sind(-30);L2+L2/7*cosd(-10) L2/7*sind(-10);L2+L2/7*cosd(10) L2/7*sind(10);L2+L2/7*cosd(30) L2/7*sind(30);L2+L2/7*cosd(50) L2/7*sind(50);L2+L2/7*cosd(70) L2/7*sind(70);L2+L2/7*cosd(90) L2/7*sind(90);L2+L2/7*cosd(100) L2/7*sind(100);L2/5*cosd(80) L2/5*sind(80)]



palm=[L2+L3/4*cosd(-90) L3/4*sind(-90);L2+L3/4*cosd(-100) L3/4*sind(-100);L2+L3/4*cosd(-120) L3/4*sind(-120);L2+L3/4*cosd(-140) L3/4*sind(-140);L2+L3/4*cosd(-160) L3/4*sind(-160);L2+L3/4*cosd(-180) L3/4*sind(-180);L2+L3/4*cosd(-200) L3/4*sind(-220);L2+L3/4*cosd(-240) L3/4*sind(-240);L2+L3/4*cosd(-260) L3/4*sind(-260);L2+L3/4*cosd(-270) L3/4*sind(-270);L2+L3/4 L3/4;L2+L3/3+L3/3.2*cosd(30) L3/4+L3/3.2*sind(30);L2+L3/3+L3/3.2*cosd(26) L3/4+L3/3.2*sind(26);L2+L3/3+L3/8 L3/4;L2+L3/3+L3/4+L3/2*cosd(11) L3/8+L3/2*sind(11);L2+L3/3+L3/4+L3/2*cosd(9) L3/8+L3/2*sind(9);L2+L3/3+L3/4 L3/8;L2+L3/3+L3/4+L3/2*cosd(1) L3/2*sind(1);L2+L3/3+L3/4+L3/2*cosd(-1) L3/2*sind(-1);L2+L3/3+L3/4 0;L2+L3/3+L3/4+L3/2*cosd(0) -L3/8+L3/2*sind(0);L2+L3/3+L3/4+L3/2*cosd(-2) -L3/8+L3/2*sind(-2);L2+L3/3+L3/4 -L3/8;L2+L3/3+L3/4+L3/3*cosd(-2) -L3/4+L3/3*sind(-2);L2+L3/3+L3/4+L3/3*cosd(-4) -L3/4+L3/3*sind(-4);L2+L3/3+L3/4 -L3/4;L2+L3/4*cosd(-90) L3/4*sind(-90)]



B=[0 -L1] %此點為上手臂和下手臂的連結點 設此點當連結點便不會在關節處看到類似linshape的固定點 外型較美觀



m1=[cosd(theta1) sind(theta1);-sind(theta1) cosd(theta1)]; %為一旋轉theta1的旋轉矩陣



b=B*m1+[0 L1]%設B點轉動後為b點



Arm1=arm1(:,2)-L1



Ax1=arm1(:,1)*cosd(theta1)-Arm1(:,1)*sind(theta1)



%上手臂的旋轉



Ay1=arm1(:,1)*sind(theta1)+Arm1(:,1)*cosd(theta1)+L1patch(Ax1,Ay1,'r')



C=[L2 0]



%此點為下手臂和手掌的連結點 設此點當連結點便不會在關節處看到類似linshape的固定點



m2=[cosd(theta2) sind(theta2); -sind(theta2) cosd(theta2)];



%為一旋轉theta2的旋轉矩陣



c=C*m2 %設C點轉動後為c點



Ax2=arm2(:,1)*cosd(theta2)-arm2(:,2)*sind(theta2)+b(1)



%下手臂的旋轉



Ay2=arm2(:,1)*sind(theta2)+arm2(:,2)*cosd(theta2)+b(2)



patch(Ax2,Ay2,'m')



Palm=palm(:,1)-L2



Ax3=Palm(:,1)*cosd(theta3)-palm(:,2)*sind(theta3)+c(1)+b(1)



%手掌的旋轉



Ay3=Palm(:,1)*sind(theta3)+palm(:,2)*cosd(theta3)+c(2)+b(2)



patch(Ax3,Ay3,'g')



axis equal





5.1.3.
%theta1=90,theta2=-45,theta3=-30
%L1=30 L2=26 L3=18
輸入body(30,26,18,90,-45,-30)手臂尺寸輸入之圖

5.1.4.
模擬動畫
L1=30
L2=26
L3=18
function body2(L1,L2,L3,X1,X2,X3)
arm1=[0 -L1/5;L1/5*cosd(-80) L1/5*sind(-80);L1/5*cosd(-70) L1/5*sind(-70);L1/5*cosd(-60) L1/5*sind(-60);L1/5*cosd(-50) L1/5*sind(-50);L1/5*cosd(-40) L1/5*sind(-40);L1/4*cosd(-45) L1+L1/4*sind(-45);L1/4*cosd(-30) L1+L1/4*sind(-30);L1/4*cosd(-10) L1+L1/4*sind(-10);L1/4*cosd(10) L1+L1/4*sind(10);L1/4*cosd(30) L1+L1/4*sind(30);L1/4*cosd(45) L1+L1/4*sind(45);L1/4*cosd(60) L1+L1/4*sind(60);L1/4*cosd(75) L1+L1/4*sind(75);L1/4*cosd(90) L1+L1/4*sind(90);L1/4*cosd(105) L1+L1/4*sind(105); L1/4*cosd(120) L1+L1/4*sind(120);L1/4*cosd(135) L1+L1/4*sind(135);L1/4*cosd(150) L1+L1/4*sind(150);L1/4*cosd(170) L1+L1/4*sind(170);L1/4*cosd(190) L1+L1/4*sind(190);L1/4*cosd(210) L1+L1/4*sind(210);L1/4*cosd(225) L1+L1/4*sind(225);L1/5*cosd(220) L1/5*sind(220);L1/5*cosd(230) L1/5*sind(230);L1/5*cosd(240) L1/5*sind(240) ;L1/5*cosd(250) L1/5*sind(250);L1/5*cosd(260) L1/5*sind(260);L1/5*cosd(270) L1/5*sind(270)]

arm2=[L2/5*cosd(80) L2/5*sind(80);L2/5*cosd(90) L2/5*sind(90);L2/5*cosd(100) L2/5*sind(100);L2/5*cosd(120) L2/5*sind(120);L2/5*cosd(140) L2/5*sind(140);L2/5*cosd(160) L2/5*sind(160);L2/5*cosd(180) L2/5*sind(180);L2/5*cosd(-160) L2/5*sind(-160);L2/5*cosd(-140) L2/5*sind(-140);L2/5*cosd(-120) L2/5*sind(-120);L2/5*cosd(-100) L2/5*sind(-100);L2/5*cosd(-90) L2/5*sind(-90);L2/5*cosd(-80) L2/5*sind(-80);L2+L2/7*cosd(-100) L2/7*sind(-100);L2+L2/7*cosd(-90) L2/7*sind(-90);L2+L2/7*cosd(-70) L2/7*sind(-70);L2+L2/7*cosd(-50) L2/7*sind(-50);L2+L2/7*cosd(-30) L2/7*sind(-30);L2+L2/7*cosd(-10) L2/7*sind(-10);L2+L2/7*cosd(10) L2/7*sind(10);L2+L2/7*cosd(30) L2/7*sind(30);L2+L2/7*cosd(50) L2/7*sind(50);L2+L2/7*cosd(70) L2/7*sind(70);L2+L2/7*cosd(90) L2/7*sind(90);L2+L2/7*cosd(100) L2/7*sind(100);L2/5*cosd(80) L2/5*sind(80)]

palm=[L2+L3/4*cosd(-90) L3/4*sind(-90);L2+L3/4*cosd(-100) L3/4*sind(-100);L2+L3/4*cosd(-120) L3/4*sind(-120);L2+L3/4*cosd(-140) L3/4*sind(-140);L2+L3/4*cosd(-160) L3/4*sind(-160);L2+L3/4*cosd(-180) L3/4*sind(-180);L2+L3/4*cosd(-200) L3/4*sind(-220);L2+L3/4*cosd(-240) L3/4*sind(-240);L2+L3/4*cosd(-260) L3/4*sind(-260);L2+L3/4*cosd(-270) L3/4*sind(-270);L2+L3/4 L3/4;L2+L3/3+L3/3.2*cosd(30) L3/4+L3/3.2*sind(30);L2+L3/3+L3/3.2*cosd(26) L3/4+L3/3.2*sind(26);L2+L3/3+L3/8 L3/4;L2+L3/3+L3/4+L3/2*cosd(11) L3/8+L3/2*sind(11);L2+L3/3+L3/4+L3/2*cosd(9) L3/8+L3/2*sind(9);L2+L3/3+L3/4 L3/8;L2+L3/3+L3/4+L3/2*cosd(1) L3/2*sind(1);L2+L3/3+L3/4+L3/2*cosd(-1) L3/2*sind(-1);L2+L3/3+L3/4 0;L2+L3/3+L3/4+L3/2*cosd(0) -L3/8+L3/2*sind(0);L2+L3/3+L3/4+L3/2*cosd(-2) -L3/8+L3/2*sind(-2);L2+L3/3+L3/4 -L3/8;L2+L3/3+L3/4+L3/3*cosd(-2) -L3/4+L3/3*sind(-2);L2+L3/3+L3/4+L3/3*cosd(-4) -L3/4+L3/3*sind(-4);L2+L3/3+L3/4 -L3/4;L2+L3/4*cosd(-90) L3/4*sind(-90)]

for n=[1:10]
theta1=180-1.5*n
%上手臂旋轉-90度至-75度
theta2=-150+2.5*n
%下手臂旋轉-45度至-35度
theta3=-315-4.5*n
%手掌旋轉由-10度至-30度
B=[0 -L1]
%設此點用以當連結點(上手臂和下手臂)便不會在關節處看到類似linshape的固定點 外型較美觀
m1=[cosd(theta1) sind(theta1);-sind(theta1) cosd(theta1)];
b=B*m1+[0 L1]
Arm1=arm1(:,2)-L1
Ax1=arm1(:,1)*cosd(theta1)-Arm1(:,1)*sind(theta1)
Ay1=arm1(:,1)*sind(theta1)+Arm1(:,1)*cosd(theta1)+L1
patch(Ax1,Ay1,'r')
C=[L2 0]
%設此點用以當連結點(下手臂和手掌)便不會在關節處看到類似linshape的固定點 外型較美觀
m2=[cosd(theta2) sind(theta2); -sind(theta2) cosd(theta2)];
c=C*m2
Ax2=arm2(:,1)*cosd(theta2)-arm2(:,2)*sind(theta2)+b(1)
Ay2=arm2(:,1)*sind(theta2)+arm2(:,2)*cosd(theta2)+b(2)
patch(Ax2,Ay2,'m')
Palm=palm(:,1)-L2
Ax3=Palm(:,1)*cosd(theta3)-palm(:,2)*sind(theta3)+c(1)+b(1)
Ay3=Palm(:,1)*sind(theta3)+palm(:,2)*cosd(theta3)+c(2)+b(2)
patch(Ax3,Ay3,'g')
AXIS([-30 15 20 70])
%訂出範圍 則圖形在運轉時座標固定
pause(1.0)
end



5.2



5.2.1



手指關節在同一平面上移動時 有一運動連桿未與運動結相連可視為平面的三連桿開放型的運動鏈,有一運動連桿未與運動結相連在平面上手指上的三個結點均為自由度為1的旋轉對。手指可分三部分,將手指靠近掌的部分稱前端,接著中段,靠近指尖那段稱末端,前端從整個手掌攤平開始可向掌心部分轉動約90度(與手掌夾90度)。中段與前端夾角最大為90度,中段與末端夾角最大亦約90度。



使用古魯伯公式: m=3(N-J-1)+sum(fi)=3(4-3-1)+3=3,自由度為三。



5.2.2



function body(L1,L2,L4,X1,X2,X4) %將body function 的arm3修改後成arm4>



arm1=[0 -L1/5;L1/5*cosd(-80) L1/5*sind(-80);L1/5*cosd(-70) L1/5*sind(-70);L1/5*cosd(-60) L1/5*sind(-60);L1/5*cosd(-50) L1/5*sind(-50);L1/5*cosd(-40) L1/5*sind(-40);L1/4*cosd(-45) L1+L1/4*sind(-45);L1/4*cosd(-30) L1+L1/4*sind(-30);L1/4*cosd(-10) L1+L1/4*sind(-10);L1/4*cosd(10) L1+L1/4*sind(10);L1/4*cosd(30) L1+L1/4*sind(30);L1/4*cosd(45) L1+L1/4*sind(45);L1/4*cosd(60) L1+L1/4*sind(60);L1/4*cosd(75) L1+L1/4*sind(75);L1/4*cosd(90) L1+L1/4*sind(90);L1/4*cosd(105) L1+L1/4*sind(105); L1/4*cosd(120) L1+L1/4*sind(120);L1/4*cosd(135) L1+L1/4*sind(135);L1/4*cosd(150) L1+L1/4*sind(150);L1/4*cosd(170) L1+L1/4*sind(170);L1/4*cosd(190) L1+L1/4*sind(190);L1/4*cosd(210) L1+L1/4*sind(210);L1/4*cosd(225) L1+L1/4*sind(225);L1/5*cosd(220) L1/5*sind(220);L1/5*cosd(230) L1/5*sind(230);L1/5*cosd(240) L1/5*sind(240) ;L1/5*cosd(250) L1/5*sind(250);L1/5*cosd(260) L1/5*sind(260);L1/5*cosd(270) L1/5*sind(270)]



arm2=[L2/5*cosd(80) L2/5*sind(80);L2/5*cosd(90) L2/5*sind(90);L2/5*cosd(100) L2/5*sind(100);L2/5*cosd(120) L2/5*sind(120);L2/5*cosd(140) L2/5*sind(140);L2/5*cosd(160) L2/5*sind(160);L2/5*cosd(180) L2/5*sind(180);L2/5*cosd(-160) L2/5*sind(-160);L2/5*cosd(-140) L2/5*sind(-140);L2/5*cosd(-120) L2/5*sind(-120);L2/5*cosd(-100) L2/5*sind(-100);L2/5*cosd(-90) L2/5*sind(-90);L2/5*cosd(-80) L2/5*sind(-80);L2+L2/7*cosd(-100) L2/7*sind(-100);L2+L2/7*cosd(-90) L2/7*sind(-90);L2+L2/7*cosd(-70) L2/7*sind(-70);L2+L2/7*cosd(-50) L2/7*sind(-50);L2+L2/7*cosd(-30) L2/7*sind(-30);L2+L2/7*cosd(-10) L2/7*sind(-10);L2+L2/7*cosd(10) L2/7*sind(10);L2+L2/7*cosd(30) L2/7*sind(30);L2+L2/7*cosd(50) L2/7*sind(50);L2+L2/7*cosd(70) L2/7*sind(70);L2+L2/7*cosd(90) L2/7*sind(90);L2+L2/7*cosd(100) L2/7*sind(100);L2/5*cosd(80) L2/5*sind(80)]



arm4=[L2+L4/5*cosd(80) L4/5*sind(80);L2+L4/5*cosd(90) L4/5*sind(90);L2+L4/5*cosd(100) L4/5*sind(100);L2+L4/5*cosd(120) L4/5*sind(120);L2+L4/5*cosd(140) L4/5*sind(140);L2+L4/5*cosd(160) L4/5*sind(160);L2+L4/5*cosd(180) L4/5*sind(180);L2+L4/5*cosd(-160) L4/5*sind(-160);L2+L4/5*cosd(-140) L4/5*sind(-140);L2+L4/5*cosd(-120) L4/5*sind(-120);L2+L4/5*cosd(-100) L4/5*sind(-100);L2+L4/5*cosd(-90) L4/5*sind(-90);L2+L4/5*cosd(-80) L4/5*sind(-80);L2+L4+L4/7*cosd(-100) L4/7*sind(-100);L2+L4+L4/7*cosd(-90) L4/7*sind(-90);L2+L4+L4/7*cosd(-70) L4/7*sind(-70);L2+L4+L4/7*cosd(-50) L4/7*sind(-50);L2+L4+L4/7*cosd(-30) L4/7*sind(-30);L2+L4+L4/7*cosd(-10) L4/7*sind(-10);L2+L4+L4/7*cosd(10) L4/7*sind(10);L2+L4+L4/7*cosd(30) L4/7*sind(30);L2+L4+L4/7*cosd(50) L4/7*sind(50);L2+L4+L4/7*cosd(70) L4/7*sind(70);L2+L4+L4/7*cosd(90) L4/7*sind(90);L2+L4+L4/7*cosd(100) L4/7*sind(100);L2+L4/5*cosd(80) L4/5*sind(80)]



%L1,L2,L4分別為指頭三段靠近手掌 中段 末端指頭 的長度



%theta1,theta2,theta4 分別為各結點的角度



B=[0 -L1]



%設此點用以當連結點 便不會在關節處看到類似linshape的固定點 外型較美觀



m1=[cosd(theta1) sind(theta1);-sind(theta1) cosd(theta1)];







b=B*m1+[0 L1]Arm1=arm1(:,2)-L1



Ax1=arm1(:,1)*cosd(theta1)-Arm1(:,1)*sind(theta1)



%靠近手掌那段的旋轉



Ay1=arm1(:,1)*sind(theta1)+Arm1(:,1)*cosd(theta1)+L1



patch(Ax1,Ay1,'r')



C=[L2 0]



%設此點用以當連結點 便不會在關節處看到類似linshape的固定點 外型較美觀



m2=[cosd(theta2) sind(theta2); -sind(theta2) cosd(theta2)];



c=C*m2



Ax2=arm2(:,1)*cosd(theta2)-arm2(:,2)*sind(theta2)+b(1)



%中段的部分旋轉



Ay2=arm2(:,1)*sind(theta2)+arm2(:,2)*cosd(theta2)+b(2)patch(Ax2,Ay2,'m')



Arm4=arm4(:,1)-L2



Ax3=Arm4(:,1)*cosd(theta3)-arm4(:,2)*sind(theta3)+c(1)+b(1)



%指頭的旋轉



Ay3=Arm4(:,1)*sind(theta3)+arm4(:,2)*cosd(theta3)+c(2)+b(2)patch(Ax3,Ay3,'g')



axis equal



pause(1.0)



end



食指食指軌跡



L1=4 L2=2.5 L4=2



theta1=100-8.8*n theta2=20-9*n*2 theta3=25-9*n*3




中指中指軌跡



L1=4.5 L2=3 L4=2.5



theta1=100-9*n theta2=20-9.6*n*2 theta3=25-8.5*n*3




無名指無名指軌跡



L1=4 L2=2.7 L4=2.4



theta1=100-7.8*n theta2=20-9.5*n*2 theta3=25-9*n*3




小姆指小姆指軌跡



L1=3.5 L2=2 L4=1.5



theta1=100-9*n theta2=20-10*n*2 theta3=26-9.6*n*3




5.2.3.



設投出的球速為150km/hr由於投球動作的組成包括了跨步,扭腰,手臂,手腕等多項因素,在此僅做指頭的分析(假設只動到指頭)將食指(其餘指頭視為與食指相同)的運動情形,由右向左投擲,令球出手時的速度等於指尖速度。



指頭上各節之長度分別為[4.5 2.5 2]cm



角度(對原點逆時針為正)[90 125 130]度



由假設算出



角速度[.462 .458 .449]rad/s



角加速度 [.501 .486 .453]rad/s^2




dyad_draw之函式分析



紅色為各指節點之加速度,藍色為速度,綠色為合速度圖形如下



投球的分析圖

2007年3月30日 星期五

機動學作業四

4.1
程式碼
AXIS([-80 80 -80 80])
L=34+10
x=[L/sqrt(3)*cosd(30) 0 -L/sqrt(3)*cosd(30) L/sqrt(3)*cosd(30)]
y=[-L/sqrt(3)*sind(30) L/sqrt(3) -L/sqrt(3)*sind(30) -L/sqrt(3)*sind(30)]
TRI=line(x',y')
xlabel('X');
ylabel('Y');
grid on
for n=1:360
rotate(TRI,[0 0 1],-1,[x(3) y(3) 0])
pause(0.0034)
end;

for n=1:360
rotate(TRI,[0 0 1],-1,[x(2) y(2) 0])
pause(0.0034)
end;

for n=1:360
rotate(TRI,[0 0 1],-1,[x(1) y(1) 0])
pause(0.0034)
end;
三角形連桿
4.2
程式碼
x=[0 0]
y=[10 0]
for n=1:12
m=[cosd(30*n) sind(30*n); sind(-30*n) cosd(30*n)];
Y=y*m
linkshape([Y(1) Y(2)],[x(1) x(2)],4)
axis equal
line([Y(1) 15],[Y(2) 0],'Color','g','LineWidth',1.5)
end
連桿羽彈簧
藍色線為彈簧的變化

4.3
4.3.1
程式碼
A=[0 0]
B=[3 4]
C=[13 4]
D=[10 0]
linkshape(A,D,2)
linkshape(A,B,3)
linkshape(B,C,1.5)
linkshape(C,D,2)
四連桿

4.3.2
程式碼
A=[0 0]
B=[3 4]
C=[13 4]
D=[10 0]
linkshape(A,D,2)
for n=[1:12]
m=[cosd(30*n) sind(30*n); sind(-30*n) cosd(30*n)];
b=B*m
m1=[cosd(30*n) sind(30*n); sind(-30*n) cosd(30*n)];
c(1)=b(1)+10
c(2)=b(2)
linkshape([c(1) c(2)],D,2)
linkshape([b(1) b(2)],[c(1) c(2)],1.5)
linkshape(A,[b(1) b(2)],3)
axis equal
pause(2.0)
end
四連桿軌跡

2007年3月22日 星期四

機動作業三

P3.1
1.每上升5cm時之手臂、手肘與關節之對應位置與軌跡 。
位置與軌跡圖
x=[20:5:60]
y=atand(x./15)
c=x.^2+15^2
C=(c+31^2-37^2)./[2*sqrt(c).*31]
z=acosd(C)
X=31*cosd(180-y-z)
Y=31*sind(y+z)
m=[0 0 0 0 0 0 0 0 0]
l=[14 14 14 14 14 14 14 14 14]
line([X ;l],[Y; m])
line([X ;m],[Y; x])

x1=[20:5:60]
y1=atand(x1./15)
c1=x1.^2+15^2
C1=(c1+31^2-37^2)./[2*sqrt(c).*31]
z1=acosd(C1)
X2=-31*cosd(180-y1-z1)
Y2=31*sind(y1+z1)
m2=[0 0 0 0 0 0 0 0 0]
l2=[-14 -14 -14 -14 -14 -14 -14 -14 -14]
line([X2 ;l2],[Y2; m2])
line([X2 ;m2],[Y2; x1])
axis equal
GRID ON
2.求出拳頭最大移動範圍為40cm
3.有設定出肩膀的寬度 脖子和頭的高度 而頭部的存在 使得拳頭在往上移動時 會受到影響 因此兩隻手和肩膀所形成的並不是三角形 主要是手肘和手臂仍會有角度 無法完全打直

P3.2
1.日常生活應用的器材中
旋轉對 為低對的運動結 2D或3D 亦稱R型對 兩剛體共軸 僅能迴轉 無法平移 DOF=1
門閂

2.高對運動結
兩物件接觸於一點或一線 作用力集中
保齡球與跑道的接觸

低對運動結
物體間的接觸屬於面或多面 接觸面積大 應力較小
小朋友溜滑梯

3.型式閉合結 藉運動結形狀限制另一物的運動範圍
伸縮喇叭
外力閉合結 藉由外在之力,使兩桿維持接觸狀態
避震器

P3.3
假設人在圓內 手打平直 有軀幹的長度 共四肢外接著圓
此圓的半徑和人體的手臂或臂長增加而增大
人外接圓

以人的身高80cm輸入結果如上

2007年3月20日 星期二

歡迎光臨!!

感謝大家的支持~

~鐺鐺鐺鐺~

新開版 希望這小版 未來能共同充實你我!!

2007年3月16日 星期五

機動學作業二

P2.1
連桿的連接點,又稱為運動結,其連接形式亦稱為對偶,對偶可再分成高對和低對,高和低的差別在於兩物件接觸時應力的高低
Ball-and-Socket Joint
而左上圖所顯示的臀部的運動結為S-型對,連接時僅能約束在XYZ三軸方向旋轉,無法平移
三度空間中的六個自由度減去平移的三個自由度後,其自由度為3,DOF=3,拘束度亦為3
解剖學網站
此網站提供解剖學相關資料,內容豐富 包括有骨骼系統 肌肉 神經系統 消化系統 淋巴系統等等方面的圖示與介紹,對身體的部位點選後還會有進一步的說明

P2.2
科技始終來自於人性!!
汽車的演進來自於人類想輕鬆且快速地到達自己想到的地區--->以風力來帶動車子前進 看天氣行車--->蒸汽來帶動車子的發動 研發失敗--->構想被應用在火車上,而蒸氣火車成功了蒸氣火車圖
--->汽油來啟動車子的汽車,但可以行使的距離很短暫--->西元1885年德國人賓士製作了汽油引擎的汽車外觀是三輪形圖示 --->當時汽車製造緩慢,價錢昂貴,不普遍--->美國人福特想到成立生產線(名為魯殊的汽車生產廠房提供一個能自給自足,所有原料到汽車製成品,因而降低代工配件的成本)品質佳 生產量大大增加,日後的汽車生產更容易,大眾的接受度高--->環保的重視加上化石燃料趨减,價格上漲,開發出油電混合車(車子設計追求低耗油量)

p2.3
利用Matlab計算數值
M=34+10
M =44
>> X=[1:1:10]
X = 1 2 3 4 5 6 7 8 9 10
>> y=X.^(1/M)
y = 1.0000 1.0159 1.0253 1.0320 1.0373 1.0416 1.0452 1.0484

1.0512 1.0537
>> plot(X,y)
結果

2007年3月1日 星期四

機動學作業一

題目一
姓名:張延瑋
學號:B94611034
本學期聯絡電話:0919586797
e-mail網址:yw0403@gmail.com
第一次這樣使用部落格 還有點生疏 相信往後頻繁的使用部落格繳交作業 更能體會到其效率以及方便性
現今的機械基本原理多由槓桿 輪軸 滑輪 楔形錐與螺旋等單元機械衍生組合來的 看似簡單的機械具有著潛在的力量 帶給人生活的一大進步
作業二
機械設計之能否成功取決於是否選用了最適當的機械機構,而機構設計和分析能力的訓練乃是任何一位機械工程師養成教育中最重要的一環。http://mmrl.cgu.edu.tw/rehab/lab/mdesign.htm 此網頁建構一個結合機構設計及網路多媒體技術之「機動學輔助教學系統」,其主要設計理念在於將機動學相關之原理、設計案例,以文字、圖片、動畫方式透過網路多媒體技術呈現

機構與機械兩者,其所代表的意義,不完全相同;機構可說是機械中的一個組合,即機構是屬於機械中的一份子,而機械則是代表整體,泛指能使能量轉變為功者
http://mmrl.cgu.edu.tw/rehab/mme/rehab/organize/chap2.htm 此網頁對機電整合做基本的介紹