Module 3 - Fall 2023
Module 3 - Fall 2023
Engineering
MODULE 3
2. Duhamel’s Integral
4. Numerical methods
This is the solution to equation of motion for undamped free vibration in the form
of Simple Harmonic Motion with an angular velocity, ωn
4
Initial Conditions
The constants A and B in the solution can be determined by evaluating the solution at two
different times or more commonly from the velocity and displacement at time t = 0
……… 2
At 𝑡=0 , u ( 𝑡 )= u ( 0 )= 𝐴𝐶𝑜𝑠 0+ 𝐵𝑆𝑖𝑛0= 𝐴
⇒ 𝐴 =𝑢( 0)
Similarly 𝑢˙ ( 𝑡 ) =− 𝐴 𝜔𝑛 𝑆𝑖𝑛(𝜔𝑛 𝑡)+𝐵 𝜔 𝑛 𝐶𝑜𝑠 (𝜔 𝑛 𝑡 )
At 𝑡=0 , 𝑢˙ ( 𝑡 )= 𝑢˙ ( 0 )=− 𝐴 𝜔 𝑛 𝑆𝑖𝑛0 +𝐵 𝜔𝑛 𝐶𝑜𝑠 0=𝐵 𝜔𝑛
𝑢(0
˙ )
⇒ 𝐵=
𝜔𝑛
Using the value of A & B in equation 2 results in:
𝑢˙ ( 0 )
𝑢 ( 𝑡 ) =𝑢 ( 0 ) 𝐶𝑜𝑠 ( 𝜔 𝑛 𝑡 ) + 𝑆𝑖𝑛 ( 𝜔𝑛 𝑡 ) … … … 𝟑
𝜔𝑛 5
Displacement Response History: Undamped free vibration
𝑢(0)
˙
𝑢(𝑡 )=𝑢(0)𝐶𝑜𝑠(𝜔𝑛 𝑡)+ 𝑆𝑖𝑛(𝜔 𝑛 𝑡)
𝜔𝑛
u(t)
t 6
Natural frequencies of free vibration
= √𝑘/𝑚
The Natural* Circular Frequency, can be calculated by using relation, 𝜔 𝑛
𝑘 canalso be determined using deflection of a systemunder gravity load as follows
√ [ ]
√ [ ]
𝑊 𝑊
𝜔𝑛 =
𝑢 𝑠𝑡
𝑚
=
𝑢 𝑠𝑡
(𝑊 /𝑔 )
=
√ 𝑔
𝑢 𝑠𝑡
* The term natural comes from the fact that these properties are intrinsic to the system. These
properties depend solely on the mass and the stiffness of the system and are constant. these do not
depend on the initial conditions for a linear elastic system.
7
Natural Period of Free vibration
In earthquake engineering the preferred measure of the dynamic characteristic of structures
is the Natural Period of free vibration, Tn , measured in units of seconds.
1 2𝜋 𝑚 𝑢𝑠𝑡
𝑇𝑛= = =2𝜋 =2𝜋
𝑓𝑛 𝜔𝑛 𝑘 𝑔 √√
8
Natural period of free vibration
Torsion
E-
W
Steel structure, 26 stories.
Periods of vibration:
9
Natural period of free vibration
Periods of vibration:
18.2
Transverse: 18.2 s s
3.81 s
Vertical : 10.9 s 4.43 s
Longitudinal: 3.81 s
Torsional: 4.43 s
10
Natural period of free vibration: Rule of thumb
This rule is not very accurate, specially when building height exceeds 11 story height,
but is good for understanding the general ballpark.
11
Problem 3.1
A heavy table is supported by flat steel legs. Its natural
period in lateral vibration is 0.5 sec. When a 50-lb
plate is clamped to its surface, the natural period in
lateral vibration is lengthened to 0.75 sec. What are
the weight and the lateral stiffness of the table?
Solution Initial time period ,𝑇 𝑛 1= 2 𝜋
√ 𝑚
𝑘
∨ 0.5= 2 𝜋
√ 𝑚
𝑘
…………𝟏
𝑇 𝑛 2= 2 𝜋
√ 𝑀
𝑘
∨0.75 =2 𝜋
√ (𝑚+ 1.55 )
𝑘
…………𝟐
Dividing Eqn 1 by 2,
0.5
0.75
=
√ 𝑚
𝑚+1.55
∨0.444=
𝑚
𝑚+1.55
or
or
𝐄𝐪𝐧 𝟏 ⇒ 0.5=2 𝜋
√ 𝑚
𝑘
𝑜𝑟 0.0063 𝑘=𝑚 𝑜𝑟 𝑘 =
1.24
0.0063
=196 lb / ft 12
Alternate mathematical form of expressing solution to EOM for
free vibration
+
A simpler and more instructive form of equation 2 can be obtained in replacing the
arbitrary constants A and B by two new constants.
A= CCos ϕ and B = -CSin ϕ
√ [ ]
Eqn 6 ⇒𝑢0= √ 𝐴 2+𝐵 2= [𝑢 ( 0 ) ] 2+
𝑢˙ ( 0 )
𝜔𝑛
2
14
Equivalent Static force
One of the most appealing approach for analyzing a structural system subjected to
dynamic loading is to use Equivalent Static force approach.
According to this approach, at any instant of time, t , the equivalent static force, fs ,
is the external force that will produce the deformation u at the same t in the stiffness
component of structure (i.e., the system without mass and damping).
𝑓 𝑠 (𝑡)=𝑘.𝑢(𝑡)
15
Equivalent Static force
A mass of 5 kg is dropped onto the end of a cantilever beam with a velocity of 0.5 m/s, as shown
in Figure. The impact causes vibrations of the mass, which sticks to the beam. The beam is made
of steel (E= 210*109 N/m2), is 3 m long, and has a moment of inertia I= 3*10–6 m4.
Determine:
1. Natural time period and displacement response of beam
2. Peak dynamic displacement
3. Max shear force and Bending Moment in the beam
17
Solution: Problem 3.2
−6
3 𝐸𝐼 3 (210 ∗ 10 9 N / m 2)( 3 ∗10 m 4) N
𝑘= = = 70,000
𝐿3 (3 m )3 m
𝜔𝑛 =
𝑘
𝑚 √ √
=
70000
5
=118 rad/ s
⇒
Let u(t) represent the displacement of the mass, which is measured positive downward from the
equilibrium position of the mass after it is attached to the beam.
The beam is in equilibrium at t = 0 when the particle hits, therefore, u is measured from the
equilibrium position of the system with the particle attached, thus
˙ (0 )
𝑢
𝑢 (0 )
19
Solution: Problem 3.2
√ [ ]
𝑢˙ ( 0 ) ⇒
𝑢 0= [ 𝑢 ( 0 ) ] 2+ 2
𝜔𝑛
⇒ W+fso = 5*9.81+574
=623.1 N
3
-623.1 N SFD
21
Exercise 3.1
3. An electromagnet weighing 3,000 lb is at rest while holding an automobile of weight 2,000 lb in
a junkyard. The electric current is turned off, and the automobile is dropped. Assuming that the
crane and the supporting cable have an equivalent spring constant of 10,000 lb/in.,find the
following: (a) the natural Time period of vibration of the electromagnet, (b) the resulting motion of
the electromagnet, and (c) the maximum tension developed in the cable during the motion
Ans: Tn= 0.175 s, u(t)=0.5 Cos (35.92t) -in, Tmax= 5000 lb
Ans: 1656 lb
23
Free vibrations Response of
Viscously Damped systems
24
Viscously Damped Free Vibration
The equation of free vibration for damped free vibration has the form 𝑚 𝑢+𝑐
¨ 𝑢+𝑘𝑢=0
˙
The solution to this equation will be taken in the same form as for the undamped form
st
i.e 𝑢 ( 𝑡 ) =G e ˙ ( 𝑡 )=G s e st
⇒ 𝑢
2
¨ ( 𝑡 )=G 𝑠 e
⇒𝑢
st
2 𝑠𝑡 𝑠𝑡 𝑠𝑡
Substituting this value in equation of motion result in: 𝑚(𝑠 𝐺𝑒 )+𝑐(𝑠𝐺𝑒 )+𝑘𝐺𝑒 =0
2 𝑠𝑡
or( 𝑠 𝑚+ 𝑠𝑐+ 𝑘)( 𝐺 𝑒 )=0
2
G e ≠ 0 ⇒ 𝑠 𝑚+ 𝑠𝑐 +𝑘= 0
st
Since
−𝑐 ± √ 𝑐 − 4 𝑘𝑚
√
2 2
𝑐 𝑐 4 𝑘𝑚
𝑜𝑟 𝑠= =− ± −
2𝑚 2𝑚 4𝑚
2
4𝑚
2
√ √[ ]
2 2
𝑐 𝑐 2 𝑐 𝑐
k = m ω n 2 ⇒ 𝑠=− ± − 𝜔𝑛 ¿ 𝑠=− 𝜔 ± 𝜔 1−
2𝑚 4 𝑚 2
2𝑚𝜔𝑛 𝑛 𝑛
2𝑚𝜔 𝑛
25
Viscously Damped Free
Vibration
√[ ] √[ ]
2 2
𝑐 𝑐 𝑐 𝑐
𝑠=− 𝜔𝑛 ± 𝜔𝑛 1 − 𝑜𝑟 𝑠=− 𝜔 𝑛 ± 𝜔𝑛 1 −
2𝑚 𝜔𝑛 2 𝑚 𝜔𝑛 𝑐𝑐𝑟 𝑐 𝑐𝑟
2 𝑚 𝜔 𝑛=𝑐 𝑐𝑟 is known as Critical damping coefficient
Also Where (Greek alphabet for Zeta) is known as Damping ratio or fraction of critical damping
𝑠=− 𝜁 𝜔𝑛 ± 𝜔𝑛 √ 1 −𝜁 2 ¿ 𝑠=− 𝜁 𝜔 𝑛 ± 𝜔 𝐷
Where ω n √ 1− 𝜁 2
=𝜔 𝐷is known as Damped Natural Frequency
26
Viscously Damped Free Vibration
There are 3 forms of solutions available depending on the magnitude of the
Damping Coefficient, c
27
Critically damped Over damped systems
systems
𝑐=2𝑚 𝜔 𝑛 𝑐 >2 𝑚 𝜔𝑛
There is no vibration in the response. The Like critically damped systems, there is no
vibration. However, the structure returns to its
structure returns to its initial position without
initial position slowly as compared to critically
vibrating about the zero position but in the
damped system.
shortest time.
29
Free vibration response of under damped, critically damped , and
over damped systems videos 1,2,3
31
Free vibration Response of
Underdamped systems
32
Displacement Response: Underdamped systems
The solution for the under-damped system have the form
𝑢(𝑡 )= e
− 𝜁 𝜔𝑛 𝑡
[ A . 𝐶 os(𝜔 𝐷 𝑡)+ 𝐵 . 𝑆𝑖𝑛(𝜔 𝐷 𝑡)where
] 𝜔 𝐷=𝜔𝑛 √ 1 −𝜁 2
This is similar to the case for undamped free vibration except that the frequency is slightly smaller and
there is a decay of the response with time.
Again the constants A and B may be found from the solution at two different times or from the initial
[ ]
conditions at time t = 0
−𝜁 𝜔𝑛 𝑡 𝑢(0)+𝑢(0)𝜁
˙ 𝜔𝑛
𝑢(𝑡 )=𝑒 𝑢(0)𝐶𝑜𝑠(𝜔 𝐷 𝑡)+ 𝑆𝑖𝑛(𝜔 𝐷 𝑡)
𝜔𝐷
Peak dynamic displacement ,uo , and equivalent static force, fso
√ [ ] √ [ ]
2 2
𝑢˙ ( 0 ) +𝜁 𝜔 𝑛 𝑢(0) 𝑢˙ ( 0 ) +𝜁 𝜔𝑛 𝑢 (0)
𝑢𝑜 = [ 𝑢 ( 0 ) ] 2+ ⇒ 𝑓 𝑠𝑜=𝑘𝑢0=𝑘 [ 𝑢 ( 0 ) ] 2+
𝜔𝐷
𝜔𝐷
33
Underdamped Free Vibration
Problem 3.3
For the cantilever beam in Problem 3.2
Develop and solve the EOM and determine:
1. Damped Natural time period
2. Amplitudes of dynamic displacement and Equivalent static force
𝑢(𝑡 )=𝑒
−𝜁 𝜔𝑛 𝑡
[ 𝑢 (0)𝐶𝑜𝑠(𝜔 𝐷 𝑡)+
𝑢˙ (0)+𝑢 (0)𝜁 𝜔𝑛
𝜔𝐷
𝑆𝑖𝑛 (𝜔 𝐷 𝑡)
]
˙ ( 0 )=0.5 m / s
From Problem 3.2, u(0) = 7*10-3 m and 𝑢
2𝜋 2𝜋
𝑇 𝐷 = = =0.00532 s
𝜔𝐷 117.96
8 uo
6
Displacement -mm
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
-2
-4
-6
-8
-10
Time-s
Solution: Problem 3.3 contd…
√ [ ]
2
𝑢˙ ( 0 )+𝜁 𝜔𝑛 𝑢(0)
P eak dynamic disp . ,𝑢0= [ 𝑢 ( 0 ) ] 2+
𝜔𝐷
𝑁
⇒ 𝑓 𝑠𝑜=𝑘𝑢0=70000 ∗0.00827 m=579 N
𝑚
Effect of damping ratio on response of civil engineering structures
Damping has the effect of lowering the natural frequency from ωn to ωD and
lengthening the natural period from Tn to TD.
These effects are negligible for damping ratios below 20%, a range that includes most
structures, as shown in Figure. For most structures the damped properties ωD and TD
are approximately equal to the undamped properties ωn and Tn, respectively.
𝜔 𝐷=𝜔𝑛 √ 1− 𝜁 2
𝜁 <1 ⇒ 𝜔 𝐷< 𝜔 𝑛
Decay of Response Due To Damping
The decay observed in the response of a structure to some initial disturbance can be used
to obtain a measure of the amount of viscous damping ,c, present in the structure.
Consider two successive positive peaks ui and ui+1 during viscously damped free vibration
response (shown in the figure) occurring at time ti and ti+1, respectively.
ti
ti Following relation can be derived that:
uui =
i
uui+1
i+1
u(t)
tti+1
i+1
t 39
Decay of Response Due To Damping
Taking the natural logarithm of both sides, we get the so-called logarithmic decrement of damping,
δ, defined by the following equation.
𝛿=𝑙𝑛
𝑢𝑖
𝑢𝑖 +1( )
[ ( )] 2 𝜋𝜁
2 𝜋𝜁
√ 1− ζ 2 ¿ 𝛿=
⇒ 𝛿=𝑙𝑛 𝑒 √ 1− 𝜁 2
≈
Therefore, for Civil engineering systems
40
Decay of Response after j no. of cycles
Decay of response after specific no of j cycles can be determined by considering peak
displacements over a range of j+1 peaks. ti
u1
. . ….. uui
i
uuj+1
i+1
u(t)
. ….. =
j tcycles
𝑙𝑛
( 𝑢1
𝑢 𝑗 +1 ) 𝑗δ
=ln 𝑒 = 𝑗 δ
i+1
t- s
( )
𝑢1 √
( )
𝑢1
2
2 𝜋𝜁 1− 𝜁
⇒ 𝑙𝑛 = 𝑗. ¿ 𝑗= 𝑙𝑛
𝑢 𝑗 +1 √1 − 𝜁 2
2 𝜋𝜁 𝑢 𝑗+ 1
41
Problem 3.4 :A free vibration test was conducted on an empty water tank
shown in figure. A force of 60 kips, applied through a cable attached to the tank,
displace the tank by 2 " in Horizontal direction.
The cable is suddenly cut and the resulting vibration is recorded. At the end of 5 cycles, which
complete in 2.55 sec., the amplitude of displacement is 0.9 ".
Ignore the vertical vibration of tank and compute the following: 2"
1. Damping ratio
2. Natural period of undamped vibration
3. Stiffness of structure
4. Weight of tank 60 k
5. Damping coefficient
6. Number of cycles to reduce the displacement
60o
amplitude to 0.5 ″
42
Solution: problem 3.4
u1 = 2 in.
⇒ ζ =0.0254 or 2.54%
2. Tn =?
5 cycles of damped vibration takes 2.55 s
or 1 cycle of damped vibration takes 2.55/5= 0.51 s ⇒TD =0.51 s
= 0.5098 s
3. k =?
43
Solution: problem 3.4 (contd…)
4. Weight of empty tank, W =?
𝑘𝑔 𝑘𝑔 2
𝑇
√ √
¿𝑊 = = 𝑛
𝑘 𝑘𝑔 ¿ ∗ 𝑘𝑔
( )
2 2
𝜔𝑛 2𝜋
𝜔𝑛= = 4𝜋 𝑇𝑛
2
𝑚2 𝑊
0.51 lb ft
𝑊= 2
∗ 18000 ∗ 32.2 2 =74,875 lb=74.9 k
4𝜋 ft s
5. c=?
𝑐=ζ ∗2 𝑚 𝜔 𝑛=ζ ∗2 𝑚 ∗
74,875
2𝜋
=
4 ∗ 0.0254 ∗ ( 32.2 )
¿ 𝑐 = 1455 lb .
s
ft
𝑇𝑛 0.51
6. cycles required for displ. amplitude to reduce from 2 in. to 0.5 in. , j=?
√ 𝑢1
( ) 𝑢1
( ) ( )
2
1− 𝜁 1 1 2
𝑗= 𝑙𝑛 = 𝑙𝑛 = 𝑙𝑛 =8.69∨9𝑐𝑦𝑐𝑙𝑒𝑠
2𝜋𝜁 𝑢 𝑗+1 2𝜋𝜁 𝑢 𝑗+ 1 2 𝜋 ∗ 0.0254 0.5
44
Exercise 3.2
1. Consider the R.C. slender tower shown in Figure, which vibrates in the transverse
direction shown in the figure. An estimate for the first natural frequency of this system is
0.15 Hz. The logarithmic decrement values measured for the tower with uncracked
reinforced concrete material and cracked reinforced concrete material are 0.04 and 0.10,
respectively. If a wind gust induces an initial displacement of 0.5 m and an initial
velocity of 0.2 m/s, determine the peak displacement amplitudes in the cases with
uncracked concrete material and cracked concrete material
Ans: uo=0.5429 m and 0.5426 m for uncracked and cracked concrete, respectively ;
therefore, the maximum displacement remains virtually unchanged.