Universiti Malaysia Pahang

Faculty of Manufacturing & Mechatronic Engineering Technology

MATLAB Assignment
Assessment Task
Course Code BFM3103 Course Name Vibration
Week/Session/Year Week 12 / Sec 01 / Lecturer En. ZUBAIR BIN KHALIL
20222023

Student Number Name Signature

FB20065 Wong Jun Yin
FB20037  Zaqwan Hazim bin Abu Zarin 
FA20017  Muhammad Amirul Husni Bin Shaharuzaman 

Mark/grade given (For official use only)

Marker’s comments

Since this work counts toward your formal assessment for this course, please write your name and student
number were indicated above, and sign. Attach this cover sheet to the front of your submission, so that
your name and student number can be seen without any cover needing to be opened. We declare that this
assessment item is our own work, except where acknowledged, and has not been submitted for academic
credit elsewhere, and acknowledge that the assessor of this item may, for the purpose of assessing this
item:
• Reproduce this assessment item and provide a copy to another member of the University; and/or,
• Communicate a copy of this assessment item to a plagiarism checking service (which may then retain a
copy of the assessment item on its database for the purpose of future plagiarism checking).
We certify that we have read and understood the University Rules in respect of Student Academic
Misconduct.
Date:
 M atLab Assignment

TABLE OF CONTENT

Index Pages
1. Introduction
2. Methodology (coding)
3. Result (graph and answer)
4. Analysis (manual calculation)
5. Conclusion
 M atLab Assignment
INTRODUCTION

METHODOLOGY (CODING)

1. For a vibrating system, m = 10 kg, k = 3500 N/m, and c = 45 Ns/m. A

harmonic force of amplitude 180 N and frequency 3.5 Hz acts on the
mass. If the initial displacement and velocity of the mass are 0 mm and
10 m/s, plot the complete solution representing the motion of the mass
using MATLAB for at least 3 cycles.
m=10
k = 3500
c = 45
F0= 180
f=3.5
v0=10

w= 2.*pi.*f
wn=sqrt(k./m)
dst=F0./k
r=w./wn
cc=2.*m.*wn
z=c./cc
wd=wn.*(sqrt(1-z.^2))
x=dst./(sqrt(((1-r.^2).^2)+(2.*z.*r.^2)))
p=atand((2.*z.*r)./(1-r.^2))

%total solution

%x0 = X0 cos(p0) + xcos(wt-p), x0= 0

% 0=X0 cos(p0) + x cos(p), x=0.0744
k = x.*(cos(p)) % X0=X0 cos(p0)
% v = -z*wn*X0*cos(p0) + wd*X0*sin(p0) + w*x*sin(p)
h = (5-((-z).*wn.*(-k))-(w.*x.*sin(p)))./wd
%(X0*cos(p0))^2 + (X0*sin(p0))^2, 1=cos2(p0)+sin2(p0)
X0 = sqrt(k+h)

% X0sin(p0)/X0cos(p0)
p0 = atand(h./k)
%x(t)=X0*e^(-z*wn*t)*cos(wd*t-p0)+xcos(w*t-p)
t = [0:pi./180:2.*pi];
xt=X0.*exp(-z.*wn.*t).*cos(wd.*t-p0)+x.*cos(w.*t+p)
plot(t,xt)
 M atLab Assignment
2. A four-cylinder automobile engine is to be supported in parallel using on
four shock mounts. The engine block assembly has a mass of 225 kg. if the
unbalanced force generated by the engine is given by 900 sin100t N,
design the shock mounts so that the amplitude of the vibration is less than
30 mm. Then plot the total solution of the system using MATLAB.

3. A suspension system has a spring constant of 40,000 N/m and damping ratio
of 0.1. The car are moving on a road profile with an amplitude of 50 mm
and about 6 m to complete one period of cycle sinusoidally. Write a
MATLAB program to find the displacement amplitude of the car for the
following conditions;
a. m = 600 kg (empty), m = 1000 kg (loaded)
b. velocity is 10 km/h, 50 km/h, and 100 km/h

k=40000; %spring constant

z=0.1; %zeta
x=0.05; %X amplitude 50mm
p=6; %period of 6 meter

m=600; m2= 1000; %masses

v= [10, 50, 100] %velocity
f=(v*1000/3600)*(1/p)

%natural freq. for both masses

wn1=sqrt(k/m)
wn2=sqrt(k/m2)
%frequency for all velocities
w=2*pi*f
%ratio of all frequencies
rwn1= w/wn1
rwn2= w/wn2
r=[rwn1,rwn2]

%amplitude ratio
xy=sqrt((1+(2*z*r).^2)./((1-r.^2).^2+(2*z*r).^2))
y=0.05./xy %displacement amplitude Y
 M atLab Assignment

RESULT (GRAPH AND ANSWER)

Q1
m = 10
k = 3500
c = 45
F0 = 180
f = 3.5000
v0 = 10
w = 21.9911
wn = 18.783
dst = 0.0514
r = 1.1755
cc = 374.1657
z = 0.1203
wd = 18.5725
x = 0.07
p = -36.5259
k = 0.0288
h = 0.1845
X0 = 0.4619
p0 = 81.1293

xt = 1×361
0.4220 0.2853 0.1277 -0.0319 -0.1758 -0.2892 -0.3620 -
0.3890 -0.3709 -0.3134 -0.2260 -0.1210 -0.0119 0.0887 0.1701
0.2248 0.2490 0.2427 0.2098 0.1566 0.0914 0.0231 -0.0400
-0.0910 -0.1254 -0.1411 -0.1386 -0.1205 -0.0913 -0.0559 -0.0198
0.0125 0.0376 0.0535 0.0600 0.0580 0.0497 0.0377 0.0248
0.0133 0.0049 -0.0000 -0.0015 -0.0007 0.0010 0.0019 0.0008
-0.0032 -0.0099 -0.0187

Question 3 results
 M atLab Assignment
The result below is the answer using MATLAB to calculate,
v = 1×3
10 50 100
f = 1×3
0.4630 2.3148 4.6296
wn1 = 8.1650
wn2 = 6.3246
w = 1×3
2.9089 14.5444 29.0888
rwn1 = 1×3
0.3563 1.7813 3.5626
rwn2 = 1×3
0.4599 2.2997 4.5993
r = 1×6
0.3563 1.7813 3.5626 0.4599 2.2997 4.5993
xy = 1×6
1.1445 0.4821 0.1048 1.2651 0.2552 0.0673
y = 1×6
0.0437 0.1037 0.4770 0.0395 0.1959 0.7424

Given the masses are 600 kg and 1000kg, and the velocities of 10 m/s, 50m/s and 100m/s. Hence,
the displacement amplitude, Y is 0.0437, 0.1037, 0.4770, 0.0395, 0.1959 and 0.7424.
 M atLab Assignment
ANALYSIS (MANUAL CALCULATION)
QUESTION 1

QUESTION 2
 M atLab Assignment
QUESTION 3

Figure 3.1
Given spring constant k = 40000 N/m, zeta = 0.1, road profile Amplitude = 0.05m and
period = 6m

Figure 3.2
The empty mass is 600 kg and 1000 kg, the natural frequency for each mass is 8.16 and
6.32 respectively.

Figure 3.3
Given the velocities of the car are 10 km/h, 50 km/h and 100km/h, the frequencies will be
8.16 rad/s, 14.544 rad/s and 29.089 rad/s. Then, the frequencies are divided by the natural
frequency of the car with empty load (600kg), to find the ratios which then used to find the
displacement amplitude, Y. The same step is done for the car with loaded weight (1000 kg) as
shown in the next figure.
 M atLab Assignment

Figure 3.4
The results of the handwritten calculation are tabulated as follows:
Mass, Natural Velocities, Frequencies Ratios Amplitude Displacement
kg frequencies, km/h
(Wn)

600 8.16 10 2.9089 0.3565 0.0436

50 14.544 1.7824 0.1039
100 29.089 3.565 0.1047

1000 6.32 10 2.9089 0.4603 0.0395

50 14.544 2.3013 4.6027
100 29.089 4.6027 0.7433
 M atLab Assignment
CONCLUSION

From the result generated by MATLAB and hand calculated while solving the
vibration question, it is found out that the results have a slight difference.