Learning Module 9 -Continuous Structures
Learning Module 9 -Continuous Structures
𝑑𝑥 m (unit length), EI
L
𝑃𝑑𝑥
𝜕𝑀
𝑀 V 𝑀+ 𝑑𝑥
𝜕𝑥
𝜕𝑉
𝑑𝑥 V+ 𝑑𝑥
𝜕𝑥
𝜕 2𝑣
(𝑚𝑑𝑥) ∙ 2
𝜕𝑡
Flexural vibration of uniform beams (cont…)
From the above figure, and assuming Bernoulli – Euler theory,
𝜕𝑉 𝜕2 𝑣
∴ V− 𝑉 + 𝑑𝑥 + 𝑃𝑑𝑥 − 𝑚𝑑𝑥 2 =0
𝜕𝑥 𝜕𝑡
𝜕𝑉 𝜕2 𝑣
⇒ + 𝑚 2 = 𝑃 (1) 𝑃𝑑𝑥
𝜕𝑥 𝜕𝑡
𝜕2 𝑣 𝑀
𝜕𝑀
But 𝑀 = 𝐸𝐼 2 V 𝑀+ 𝑑𝑥
𝜕𝑥 𝜕𝑥
𝜕𝑀
and 𝑉 =
𝜕𝑥
3
𝜕 𝑣
∴ 𝑉 = 𝐸𝐼 𝜕𝑉
𝑑𝑥 V+ 𝑑𝑥
𝜕𝑥 3 𝜕𝑥
𝜕4 𝑣 𝜕2 𝑣 𝜕 2𝑣
⇒ 𝐸𝐼 +𝑚 =𝑃 (2) (𝑚𝑑𝑥) ∙ 2
𝜕𝑥 4 𝜕𝑡 2 𝜕𝑡
Flexural vibration of uniform beams (cont…)
Note:
Equation (2) does not take into account shear deflections and rotary inertia.
4 𝑚𝜔2
where 𝑎 =
𝐸𝐼
𝐸𝐼𝑎 4 𝐸𝐼 2
∴𝜔= =𝐺 where G = 𝑎𝐿
𝑚 𝑚𝐿4
𝑣 𝐿, 𝑡 = 0, 𝑀 𝐿, 𝑡 = 0
∴ At 𝑥=0 ; 𝜙 0 =0 , 𝜙 ′′ 0 = 0
𝑥=𝐿 ; 𝜙 𝐿 =0 , 𝜙 ′′ 𝐿 = 0
Example 1: Simply Supported Beam (cont…)
Thus, from Equation (10),
𝜙 0 =0 ⇒ 𝐵+𝐷=0 𝜙 𝑥 = 𝐴 sin 𝑎𝑥 + 𝐵 cos 𝑎𝑥 + 𝐶 sinh 𝑎𝑥 + 𝐷 cosh 𝑎𝑥
𝜙 ′′ 0 = 0 ⇒ −𝐵 + 𝐷 = 0
∴𝐵=𝐷=𝑂
𝜙 𝐿 = 0 ⟹ 𝐴 sin 𝑎𝐿 + 𝐶 sinh 𝑎𝐿 = 0
𝜙 ′′ 𝐿 = 0 ⟹ 𝑎2 −𝐴 sin 𝑎𝐿 + 𝐶 sinh 𝑎𝐿 = 0 (a)
∴ 2𝐶 sinh 𝑎𝐿 = 0.
⟹𝐶=0
𝐸𝐼 𝜔=𝐺
𝐸𝐼
where G = 𝑎𝐿 2
∴𝜔=𝐺 𝑚𝐿4
𝑚𝐿4
∴ 𝐺 = 𝑛2 𝜋 2
𝐸𝐼
i.e. 𝜔 = 𝑛 2 𝜋 2
𝑚𝐿4
and since 𝐵 = 𝐶 = 𝐷 = 0
𝑛𝜋𝑥
∴ 𝜙𝑛 𝑥 = 𝐴 sin 𝜙 𝑥 = 𝐴 sin 𝑎𝑥 + 𝐵 cos 𝑎𝑥 + 𝐶 sinh 𝑎𝑥 + 𝐷 cosh 𝑎𝑥
𝐿
𝑛𝜋𝑥
i.e. 𝑣𝑛 𝑥, 𝑡 = sin 𝐴𝑛 cos 𝜔𝑛 𝑡 + 𝐵𝑛 sin 𝜔𝑛 𝑡
𝐿
𝑑4 𝜙𝑛 (𝓍) 2𝑚
𝜔𝑛
⇒
𝑑𝓍 4
=
𝐸𝐼
𝜙𝑛 (𝑥) (2a)
𝐸𝐼 𝜙 𝐼𝑉 𝑥 𝑓ሷ 𝑡
=− = 𝜔2
𝑚 𝜙 𝑥 𝑓 𝑡
𝑑4 𝜙𝑚 (𝑥) 2 𝑚
𝜔𝑚
𝑑𝑥 4
=
𝐸𝐼
𝜙𝑚 (𝑥) (2b)
Orthogonality of Modes of Vibration
Multiply Equation (2a) by 𝜙𝑚 (𝑥) and Equation (2b) by 𝜙𝑛 (𝑥) and subtract one from the other
and integrate from o to L.
2 −𝜔2
𝜔𝑛 𝐿 𝐿 𝑑4 𝜙𝑛 𝑑𝜙𝑚
𝑚
𝑚𝜙 𝑜 𝑥 𝜙𝑛 𝑥 𝑑𝑥 = 𝑜 𝜙𝑚 − 𝜙𝑛 𝑑𝑥
𝐸𝐼 𝑑𝑥 4 𝑑𝑥 4
𝐿
𝐿 𝑑4 𝜙𝑛 𝑑3 𝜙𝑛𝐿 𝑑𝜙𝑚 𝑑3 𝜙𝑛
𝑥𝑑 𝑚𝜙 𝑜4 𝑑𝑥 = 𝜙𝑚 3 − 𝑜 . 3 dx
𝑑𝑥 𝑜 𝑑𝑥 𝑑𝑥
𝐿 𝐿
𝑑3 𝜙𝑛 𝑑𝜙𝑚 𝑑2 𝜙𝑛 𝐿 𝑑2 𝜙𝑚 𝑑2 𝜙𝑛
= 𝜙𝑚 3 − ∙ 2 + 𝑜 ∙ 𝑑𝑥
𝑑𝑥 𝑜 𝑑𝑥 𝑑𝑥 𝑜 𝑑𝑥 2 𝑑𝑥 2
𝐿 𝑑2 𝜙𝑚 𝑑2 𝜙𝑛
= 𝑥𝑑 𝑜2 ∙ 𝑑𝑥
𝑑𝑥 2
Similarly,
𝐿 𝑑4 𝜙𝑚 𝐿 𝑑2 𝜙 2
𝑛 𝑑 𝜙𝑚
න 𝜙𝑛 4
𝑑𝑥 = න 2
∙ 2
𝑑𝑥
𝑜 𝑑𝑥 𝑜 𝑑𝑥 𝑑𝑥
Orthogonality of Modes of Vibration (Cont…)
𝜔𝑛2 − 𝜔𝑚2 𝐿
∴ න 𝜙𝑚 (𝑥)𝜙𝑛 (𝑥)𝑑𝑥 = 0
𝐸𝐼 𝑜
Thus for 𝜔𝑛 ≠ 𝜔𝑛
𝐿
⇒ = 𝑥𝑑 𝑛𝜙)𝑥( 𝑚𝜙 𝑜0 (3)
𝐿 𝑑2 𝜙𝑚 𝑑2 𝜙𝑛
𝑥𝑑 𝑜2 ∙ 𝐸𝐼 𝑥 𝑑𝑥 = 0 (5)
𝑑𝑥 2
MODAL ANALYSIS OF STRUCTURES
Forced Vibration of Beams
𝜕2 𝒱 𝜕4 𝒱
m 2 + 𝐸𝐼 4 = 𝑃(𝑥, 𝑡) (6)
𝜕𝑡 𝜕𝑥
• 𝑃(𝑥, 𝑡) excites several natural frequencies and modes of the
structure
• Every mode is a possible solution
• Therefore, dynamic deflection under arbitrary loading is a
linear combination of all possible modes.
∴ 𝑣(𝑥, 𝑡) = σ∞ 𝑛=1 𝜙𝑛 (𝑥) ∙ 𝑌𝑛 (𝑡) (7)
𝑌𝑛 (𝑡) are time-varying modal amplitudes
Forced vibration of beams
Equation (7) into Equation (6) 𝜕2𝒱 𝜕4𝒱
m 𝜕𝑡 2 + 𝐸𝐼 𝜕𝑥 4 = 𝑃 𝑥, 𝑡 (6)
𝜕 2 𝑌𝑛 (𝑡) 𝐸𝐼 𝜕 4𝜙𝑛 (𝑥) 𝑃(𝑥,𝑡)
σ∞
𝑛=1 𝜙𝑛 (𝑥) + σ∞
𝑛=1 ∙ 𝑌𝑛 (𝑡) = (8)
𝜕𝑡 2 𝑚 𝜕𝑥 4 𝑚 𝑣(𝑥, 𝑡) = σ∞
𝑛=1 𝜙𝑛 (𝑥) ∙ 𝑌𝑛 (𝑡) (7)
But;
𝐸𝐼 𝜕 4𝜙𝑛 (𝑥)
∙ 𝜕𝑥 4 = 𝜔𝑛2 (9) 𝐸𝐼 𝜙 𝐼𝑉 𝑥 𝑓ሷ 𝑡
𝑚𝜙𝑛 (𝑥)
=− = 𝜔2
𝑚 𝜙 𝑥 𝑓 𝑡
Where;
𝐿
𝑄𝑛 (𝑡) = 𝑥 𝑝 𝑥 𝑛𝜙 𝑜, 𝑡 𝑑𝑥 (13)
𝐿
𝑀𝑛 = 𝑚 𝑛𝜙 𝑜2 𝑥 𝑑𝑥 (14)
Forced vibration of beams (Cont…)
𝑄𝑛 (𝑡) : Generalized Force
𝐿
𝑄𝑛 (𝑡) = 𝑥 𝑝 𝑥 𝑛𝜙 𝑜, 𝑡 𝑑𝑥
𝑀𝑛 : Generalized Mass
𝐿
𝑀𝑛 = 𝑚 න 𝜙𝑛2 𝑥 𝑑𝑥
𝑜
𝑛2 𝜋 2 𝐸𝐼
𝜔𝑛 = 2 = 𝑛2 𝜔1
𝐿 𝑚
∴ 𝜔𝑛 = 30.84𝑛2
𝑛𝜋𝑥
𝜙𝑛 (𝑥) = 𝐶1 sin ; 𝐶1 = 1 (arbitrary)
𝐿
𝑛𝜋𝑥
∴ 𝜙𝑛 (𝑥) = sin
𝐿
Example
Beam subjected to Impulsive Load
Generalized mass (𝑀𝑛 )
𝐿 𝐿 𝑛𝜋𝑥
𝑀𝑛 = 𝑚 න 𝜙𝑛2 𝑥 𝑑𝑥 = 𝑚 න sin 2
𝑑𝑥
𝑜 𝑜 𝐿
𝑚𝐿
∴ 𝑀𝑛 = = 3 000 𝑘𝑔
2
Note: Generalized mass is the same for every mode of vibration and is exactly one half of the
mass of the beam.
Example (Cont…)
F(t) Consider a moving constant force (F(t)
𝑥 = 𝑉𝑡
𝑃(𝑥, 𝑡)
𝑃 𝑉𝑡 = 𝑃(𝑥, 𝑡)
𝑥 = 𝑉𝑡 𝑑𝑥 𝑥
Example (Cont…)
𝐿 𝐿
𝑛𝜋𝑥 𝑛𝜋𝑥
𝑄𝑛 𝑡 = න 𝜙𝑛 𝑥 ∙ 𝑝 𝑥, 𝑡 𝑑𝑥 = න sin ∙ 𝑝 𝑥, 𝑡 𝑑𝑥 𝜙𝑛(𝑥) = sin
𝐿
𝑜 0 𝐿
But 𝑥 = 𝑉 ∙ 𝑡 ;
𝐿 𝑛𝜋𝑉𝑡 𝑛𝜋𝑉𝑡 𝐿
∴ 𝑄𝑛 (𝑡) = 𝑜sin 𝐿 ∙ 𝑝 𝑥, 𝑡 𝑑𝑥 = sin 𝑥(𝑝 𝑜, 𝑡) ∙ 𝑑𝑥
𝐿
𝐿
But 0 𝑝 𝑥, 𝑡 𝑑𝑥 = 𝐹 𝑡
𝑛𝜋𝑉𝑡 (15)
∴ 𝑄𝑛 (𝑡) = 𝐹(𝑡) ∙ sin
𝐿
Example (Cont…)
𝐿
For 𝑥 = 𝑉𝑡 = (Load applied at mid-span)
2
𝐹 𝑡 ; 𝑛 = 1,5,9, ⋯ 𝑛𝜋𝑉𝑡
𝑛𝜋
𝑄𝑛 (t)= 𝐹 𝑡 sin =ቐ 0; 𝑛 𝑒𝑣𝑒𝑛 𝑄𝑛 (𝑡) = 𝐹(𝑡) ∙ sin
2 𝐿
−𝐹 𝑡 ; 𝑛 = 3; 7; 11, ⋯
Response of a structure subject to an impact 𝐼 where 𝐼 = 𝐹(𝑡) ∙ 𝑑𝑡
𝐼
𝑢= ∙ sin 𝜔𝑡
𝑚𝜔
𝑄𝑛 𝑡
𝑌𝑛ሷ 𝑡 + 𝜔𝑛2 𝑌𝑛 𝑡 = ; 𝑛 = 1,2, ⋯ ∞
𝑀𝑛
𝐼𝑛
∴ 𝑌𝑛 (𝑡) = ∙ sin 𝑛 2 𝜔1 𝑡
𝑀𝑛 𝜔𝑛
Example (Cont…)
∞
𝑣(𝑥, 𝑡) = 𝜙𝑛 𝑥 𝑌𝑛 (𝑡)
𝑛=1
∞
𝐼𝑛 sin 𝑛𝜋ൗ2
∴ 𝑣 𝐿ൗ2 , 𝑡 = ∙ sin(𝑛 2
𝜔1 𝑡)
𝑀𝑛 𝑛2 𝜔1
𝑛=1
𝐿 70 1 1
∴ 𝑣 ൗ2 , 𝑡 = sin 𝜔1 𝑡 + sin 9𝜔1 𝑡 + sin 25𝜔1 𝑡 + ⋯
3000 × 30.84 9 25
where 𝜔1 = 30.84 𝑟𝑎𝑑/𝑠
General Solutions to Moving Loads
General Solutions to Moving Loads
ሷ𝑌𝑛 𝑡 + 𝜔𝑛2 𝑌𝑛 𝑡 = 𝑄𝑛 𝑡
; 𝑛 = 1,2, ⋯ ∞ ; 𝑄𝑛 (𝑡) = 𝐹(𝑡) ∙ sin
𝑛𝜋𝑉𝑡
𝑀𝑛 𝐿
𝐿 𝐿 𝑛𝜋𝑥 𝑚𝐿
𝑀𝑛 = 𝑚 𝑛𝜙 𝑜2 𝑥 𝑑𝑥 = 𝑚 𝑜sin2 𝑑𝑥 ∴ 𝑀𝑛 =
𝐿 2
2𝐹(𝑡) 𝑛𝜋𝑉𝑡
ሷ 2
∴ 𝑌𝑛 𝑡 + 𝜔𝑛 𝑌𝑛 𝑡 = ∙ sin ; 𝑛 = 1,2, ⋯ ∞ (16)
𝑚𝐿 𝐿
𝑛𝜋𝑉
Periodic function with 𝜔
ഥ=
𝐿
2𝐹ൗ 𝑛𝜋𝑉𝑡ൗ
𝑚𝐿 ∙sin 𝐿
∴ Steady state solution is: 𝑌𝑛 (𝑡) = 2 − 𝑛𝜋𝑉ൗ 2 (17)
𝜔𝑛 𝐿
when: 𝑉𝑡 = 𝐿Τ2
𝐹𝐿3
𝛿𝑐 (𝑡) =
48 𝐸𝐼
Vibrating String
𝑣 T
𝜕𝜃
𝜃+ 𝑑𝑥
𝜕𝑥
𝜃 𝑣
T
𝑥
Vibrating String (cont…)
Let the mass per unit length =𝜌
𝜕𝜃 𝜌 𝜕2𝑣
= ∙ (1)
𝜕𝑥 𝑇 𝜕𝑡 2
𝜕𝑣
𝜃=
𝜕𝑥
𝜕2𝑣 1 𝜕2𝑣 𝑇
∴ = ; where 𝐶 =
𝜕𝑥 2 𝐶 2 𝜕𝑡 2 𝜌
𝑣 0, 𝑡 = 𝑣 𝑙, 𝑡 = 0
𝑣 0, 𝑡 = 0 ⇒ 𝐵 = 0