Universiti Teknologi Mara: Shah Alam - Selangor - Malaysia
Universiti Teknologi Mara: Shah Alam - Selangor - Malaysia
CLASS :EMD5M9A
ASIGNMENT:
Mini-project
PREPARED BY
1
Group
9
QUESTION 3
2
%graph for varying damping ratio;
wdr=3;
wn=3.5;
fo=4;
tf=10;
t=0:tf/1000:tf;
for k=1:3
if (zeta(k)<0j|zeta(k)>=1)
end
end
for k=1:3
wd=wn*sqrt(1-zeta(k)^2);
3
Ao=fo/sqrt((wn^2-wdr^2)^2+(2*zeta(k)*wn*wdr)^2);
phi=atan2(2*zeta(k)*wn*wdr,(wn^2-wdr^2));
Z1=-zeta(k)*wn-wdr*tan(phi);
Z2=sqrt((zeta(k)*wn)^2+2*zeta(k)*wn*wdr*tan(phi)+ ...
(wdr*tan(phi))^2+wd^2); % l.
Z=(Z1+Z2)/wd;
Anum=Ao*((zeta(k)*wn*Z-wd)*cos(phi)+wdr*Z*sin(phi));
Aden=Z*wd;
A=Anum/Aden;
theta=2*atan(Z);
x(k,:)=A*exp(-zeta(k)*wn*t).*sin(wd*t+theta)+Ao*cos(wdr*t-phi);
end
for k=1:3
subplot(3,1,k)
plot(t,x(k,:))
ylabel('Response x')
grid
end
xlabel('Time, seconds')
4
5
Figure 1: varies in damping ratio.
wdr=3;
zeta=0.05;
fo=4;
tf=10;
t=0:tf/1000:tf;
for k=1:3
6
end
for k=1:3
wd=wn(k)*sqrt(1-zeta^2);
Ao=fo/sqrt((wn(k)^2-wdr^2)^2+(2*zeta*wn(k)*wdr)^2);
phi=atan2(2*zeta*wn(k)*wdr,(wn(k)^2-wdr^2));
Z1=-zeta*wn(k)-wdr*tan(phi);
Z2=sqrt((zeta*wn(k))^2+2*zeta*wn(k)*wdr*tan(phi)+ ...
(wdr*tan(phi))^2+wd^2);
Z=(Z1+Z2)/wd;
Anum=Ao*((zeta*wn(k)*Z-wd)*cos(phi)+wdr*Z*sin(phi));
Aden=Z*wd;
A=Anum/Aden;
theta=2*atan(Z);
x(k,:)=A*exp(-zeta*wn(k)*t).*sin(wd*t+theta)+Ao*cos(wdr*t-phi);
end
for k=1:3
subplot(3,1,k)
plot(t,x(k,:))
ylabel('Response x')
grid
end
xlabel('Time, seconds')
7
8
Figure 2: varies in natural frequency
9
DISCUSSION
Both graph show a transient period. Transient mean an event of a short-lived burst
of energy in a system caused by a sudden change of state. In figure 1, when the damping
ratio is varied the transient period of vibration varies inversely with the damping ratio.
The length of transient period varies from about 4.5 seconds for ζ=0.05 and about 2
second for ζ=0.5. While for figure 2 when driving frequency(3rad/s) is near the natural
frequency which wn= 2 and wn=4 the transient period of vibration quite long almost 10
seconds. However for large natural frequency(16rad/s) the transient period is less than 4
second.
The response of the system differ slightly as shown on figure 1 and 2 when the
damping is added. Both graph shows when ζ=0.05 the amplitude of the response is
almost 2. While when ζ=0.5 the amplitude of the response shows less than 1. thus we can
conclude that the damping ratio effect the amplitude in inversely relationship.
10