2007年5月30日 星期三

作業十一

第一題
本週(5/24)有來上課。


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

1.設ψ凸輪之迴轉角度,ψ=100-200°為升程,轉折點為150°,起點100°,
令θ = ψ-100,β = 100及h = 5代入公式8.9,可得升程之第一段(θ=100-200°)

y(θ)=2h(θ/β)2=2(5)[(ψ-100)/100]2

2. 而升程之第二段(θ=150-200°)則仍以ψ-100=θ,β=100及h=5代入公式8.10,
得:
y(θ)=h[1-2(1-θ/β)2]=(5)[1-2(1-{ψ-100}/100]2

3. 返程之區間為260°至360°,轉折點為(260+360)/2=310°,而β=360-260=100°。設ψ-210=θ,代入返程一段的位移公式8.11,得第一區段(θ=210-285°):

y(θ)=h[1-2(θ/β)2]=(5)[1-2({ψ-260}/100]2

4. 返程第二段之區間為310至360度,其轉折點與β均與第3項同,
代入返程二段的位移公式,可以得到:
y(θ)=2h[1-(θ/β)2]=2(5)[1-({ψ-260}/100]2


%(theta:對應角度; ss:位移, vv:速度, aa:加速度)
% 計算升程對應點之資料:位移,速度及加速度
theta=100:10:200;
for i=1:length(theta)
[ss(i), vv(i), aa(i)]=parabol_cam(theta(i),100,100,1,5,0);
end;
[theta' ss' vv' aa']
figure
plot(theta,ss,'b+:')
figure
plot(theta,vv,'r+:')
figure
plot(theta,aa,'g+:')
% Demo8_2
% 計算返程對應點之資料:位移,速度及加速度
theta=260:10:360;
for i=1:length(theta)
[ss(i), vv(i), aa(i)]=parabol_cam(theta(i),260,100,-1,5,0);
end;
[theta' ss' vv' aa']
figure
plot(theta,ss,'b+:')
figure
plot(theta,vv,'r+:')
figure
plot(theta,aa,'g+:')


位移



速度



加速度



第三題
依照題義,使用到教授的pincam function
設偏置量為零

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


第四題



for theta=1:100
r(theta)=15; %0~99度的時候r不變
end
for theta=101:201
r(theta)=15+5*(theta-100)/100; %在這個範圍,衝程5被分成100份
end
for theta=201:260
r(theta)=20; %r不變
end
for theta=261:360
r(theta)=20-5*(theta-260)/100; %將衝程再次分成100段扣回
end
theta2=linspace(0,2*pi,360)
for p=1:360;s=r(1);
for q=1:359;r(q)=r(q+1);
end
r(360)=s;
polar(theta2,r)
pause(0.01)
end

2007年5月16日 星期三

作業九

B94611009 楊子羲

本週(5/3)有來上課



若以曲桿作為驅動,則令整個系統的座標如下圖

由c點座標可知,系統的限制為L^2-(10-R*sinθ)^2>=0

算出來的結果得知-47.47度到227.46度為可移動範圍

以下是我用來畫出整個移動範圍的程式



R=19 %曲桿
L=24 %連桿
X=10 %偏置量

for theta1=-47.47:1:227.46
if L^2-(X-R*sind(theta1))^2>=0 %把還是虛數的部分排除
clf
line([0,R*cosd(theta1)], [0,R*sind(theta1)]) %曲桿的程式
line([R*cosd(theta1),(L^2-(10-R*sind(theta1))^2)^(1/2)+R*cosd(theta1)],[R*sind(theta1),X]) %連桿的程式
axis equal
axis ([-50 50 -50 50])
pause(0.01)
else
end
end

2007年5月9日 星期三

作業八

B94611009 楊子羲
我本週四4/26)曾來上課。
第一題

>> [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

第二題

我不清楚此題的題義
因此我畫出 theta2=45 度的圖
再加上第一題的數據
圖片


第三題

>> move_4paths([4 3 3 5],0,0,3,0,10,0,1,0,4,2)

Qstart =
28.9550 度
圖片

Qstop =
331.0450 度
圖片

第四題

檢查連桿類型 :
>> grashof(2,[4 3 3 5])
ans =
Non-Grashof Linkage

由此可知此四連桿在theta2=[0:20:360]的時候有部分會沒有值
圖片

第五題