AML 543 EQ DYNAMICS

ASSIGNMENT II
Panaiyappan pl
MT17SDE014
 PANAIYAPPAN PL
MT17SDE014

1) A SDOF with 2% damping, and natural frequency of 0.85Hz is subjected to base acceleration
of 𝑥̈ g (t) = R/(1+0.5t). Here R is your roll number. Using program of Duhamel integral, obtain
the displacement time history of SDOF for 10 seconds?

Given data:

η=0.02

Natural frequency = 0.85Hz

Tmax = 10 seconds

Solution:

Natural frequency = 2*3.14*0.85 = 5.34 rad/seconds

Tn = (2*3.14)/5.34 =1.176 seconds

dτ = Tn/20 = 0.06 seconds

Tmax = 10 seconds (given)

R = 14

t varies from 0 to 10 seconds in an interval of 0.06.

Make an excel sheet for 𝑥̈ g(t) and save as .txt file format in the same folder that contains
MATLAB program. Then run the program and get the output results.

1
 PANAIYAPPAN PL
MT17SDE014

Time history graph

0.2
Displacement in metres

-0.2

-0.4 MATLAB output

-0.6

-0.8

-1
0 2 4 6 8 10
Time in seconds

Fig 1 DISPLACEMENT TIME HISTORY

2) For the 2DoF, M1=M2=5t & K1=10kN/m & K2=8kN/m. Obtain natural frequencies and mode
shapes. Also write the modal equations if base acceleration of 𝑥̈ g (t) =12sin (0.5t) m/sec2 is
applied?

Given data:

Two degree of freedom system

M1 = 5 tonnes

M2 = 5 tonnes

K1 =10KN/m; K2 = 8 KN/m

2
 PANAIYAPPAN PL
MT17SDE014

Solution:

Step 1

For that we determine  by |𝑘 − (𝜔2 ) ∗ 𝑚| = 0

1 = 0.844 ; 2 2.11

Step 2

Now we find the mode shapes,

Where i is the mode number, Assuming ∅1𝑖 = 1, we find ∅2𝑖.

1 1
∅1 = { } And ∅2 = { }.
1.804 −0.553

Step 3

Now we write the modal equations,

𝑎j̈ + (j2) = -𝑥𝑔̈*Pj

∑𝑛 𝑚𝑖∗∅𝑖𝑗
pj = ∑𝑛1
1 𝑚𝑖∗∅𝑖𝑗 2

so we get modal equations as,

𝑎1̈ + (0.844)2 = -7.908 sin(0.5t) –(For mode1)

𝑎2̈ + (2.11)2 = -4.1076 sin(0.5t) -(For mode 2)

The graphs are given below in fig (a) and (b) respectively from MATLAB.

(Here we have 3 value of  so we take minimum  to calculate Tmax and maximum  to calculate dtau.
We get Tmax = 251.2s and dtau = 0.1s.

Step 4

Now we use

xi = ∑𝑛1 ∅𝑖j*ai

3
 PANAIYAPPAN PL
MT17SDE014

x1 = (∅11 ∗ 𝑎1) + (∅12 ∗ 𝑎2).

x2 = ( ∅21 ∗ 𝑎1) + (∅22 ∗ 𝑎2).

The graphs are given in Fig2 (c) and Fig 2(d) respectively.

a1(t)
30

20
Displacement in m

10

0
a1
-10

-20

-30
0 50 100 150 200 250 300
Time(s)

Fig 2(A)

a2(t)
1.5

1
Displacement in m

0.5

0
a2
-0.5

-1

-1.5
0 50 100 150 200 250 300
Time(s)

Fig 2(b)

X1(t) = M1* ϕ11 + M2* ϕ12

4
 PANAIYAPPAN PL
MT17SDE014

X1(t)
40

20
Displacement (m)

0
x1
-20

-40
0 50 100 150 200 250 300
Time(s)

Fig 2(c)

X2(t) = M1* ϕ21 + M2* ϕ22

X2(t)
60
Displacement(m)

30

0
x2
-30

-60
0 50 100 150 200 250 300
Time(s)

Fig 2(d)

DISPLACEMENT TIME HISTORY

5
 PANAIYAPPAN PL
MT17SDE014