1
Tokyo Institute of Technology | Takeuchi Lab
Footfall Induced Vibration Arup Unified Method for Floors, - - PowerPoint PPT Presentation
Tokyo Institute of Technology | Takeuchi Lab Footfall Induced Vibration Arup Unified Method for Floors, Footbridges, Stairs and Other Structures Ben Sitler, PE 1 Tokyo Institute of Technology | Takeuchi Lab 2 Outline Introduction to
1
Tokyo Institute of Technology | Takeuchi Lab
2
Tokyo Institute of Technology | Takeuchi Lab 2
3
Tokyo Institute of Technology | Takeuchi Lab 3
4
Tokyo Institute of Technology | Takeuchi Lab 4
Introduction to Footfall Induced Vibration
5
Tokyo Institute of Technology | Takeuchi Lab 5
Introduction to Footfall Induced Vibration :: Loading :: Modal Properties :: Response
Modal (structural) properties: Loading:
Response
6
Tokyo Institute of Technology | Takeuchi Lab 6
Introduction to Footfall Induced Vibration :: Loading :: Modal Properties :: Response
7
Tokyo Institute of Technology | Takeuchi Lab 7
Introduction to Footfall Induced Vibration :: Loading :: Modal Properties :: Response
8
Tokyo Institute of Technology | Takeuchi Lab 8
Introduction to Footfall Induced Vibration :: Loading :: Modal Properties :: Response
9
Tokyo Institute of Technology | Takeuchi Lab 9
Introduction to Footfall Induced Vibration :: Loading :: Modal Properties :: Response
5 10 15 20 25 0.5 1 1.5 2 2.5 3
Amplification Factor (excitation frequency)/(natural frequency)
Steady State Resonant Amplification Factors for SDOF systems
0.02 damping 0.05 damping 0.10 damping 0.20 damping
Damping :: Frequency :: Modal Mass :: Mode Shape
Be suspicious of >2% Footbridges 0.5~1.5% Stairs 0.5% Floors 1~3%
10
Tokyo Institute of Technology | Takeuchi Lab 10
Introduction to Footfall Induced Vibration :: Loading :: Modal Properties :: Response
Damping :: Frequency :: Modal Mass :: Mode Shape
11
Tokyo Institute of Technology | Takeuchi Lab 11
Introduction to Footfall Induced Vibration :: Loading :: Modelling :: Response
Situation R RMS Acc @ 5Hz Low vibration 1 0.005 m/s2 Residential (night) 1.4 0.007 m/s2 Residential (day) 2-4 0.01~0.02 m/s2 Office (high grade) 4 0.02 m/s2 Office (normal) 8 0.04 m/s2 Footbridge (heavy) 24 0.12 m/s2 Footbridge (inside) 32 0.16 m/s2 Footbridge (outside) 64 0.32 m/s2 Stair (high use) 24 0.12 m/s2 Stair (light use) 32 0.16 m/s2 Stair (very light use) 64 0.32 m/s2
0.001 0.01 0.1 1 10 100 Frequency (Hz) rms acceleration (m/s
2)
Approximate threshold of human perception to vertical vibration
T
2 t t
RMS RMS Acc Peak,Harmonic RMS T
R=Acc Acc Acc Acc = = 2
12
Tokyo Institute of Technology | Takeuchi Lab 12
13
Tokyo Institute of Technology | Takeuchi Lab 13
Computing the Structural Response
14
Tokyo Institute of Technology | Takeuchi Lab 14
Force Time
Computing the Structural Response
Force Time
15
Tokyo Institute of Technology | Takeuchi Lab 15
( ) ( ) ( ) ( )
i t
my t cy t ky t P t P e
Computing the Structural Response
, F W DLF f harmonic
2 2 2 2 2 2 2 2 2 2 2
1 1 1 2 1 2 1 2 1 2
h h m m h h h h m m m m h h h h m m m m
disp f f f f f f f f f f f f f f f f f f f f
k FRF k m ic i i
SDOF Harmonic Steady State Response
excition reponse acc
F Acc FRF m
…Concrete Centre eq 4.4 ƒh = harmonic forcing freq. ƒm = modal (structural) freq.
16
Tokyo Institute of Technology | Takeuchi Lab 16
Computing the Structural Response
Force Time
SDOF Impulse – RMS Velocity Response
( ) ( ) ( ) ( )
i t
my t cy t ky t P t P e
1.43 1.3
54 [Ns]
w eff n
f I f
eff excition reponse Peak
I Vel m ( ) sin
t Peak
Vel t Vel e t
( )
RMS m
Vel RMS Vel t
ƒw = walking forcing freq. ƒm = modal (structural) freq. …Concrete Centre eq 4.10~4.13
17
Tokyo Institute of Technology | Takeuchi Lab 17
Computing the Structural Response
1 2 2 2
( ) 1
s a s s a a a disp a a a a a
k k m i c c k i c FRF K M iC k i c k m i c
s s s a p a p s s a p a p s a a a a a a a a p p p p p p p p
m y c c c c c y k k k k k y P m y c c y k k y P m y c c y k k y
2 2 2
1 (1,1) (1,2)
a disp a a a s a s a
m DMF FRF FRF FRF f f k D k D D f f
2
1 2
a a a a
f f D i f f
2
1 2
s s s s
f f D i f f
2 disp
Acc P DMF ( ) ( ) ( ) ( )
i t
My t Cy t Ky t P t P e
ƒ ,
a
P W GLF harmonic scenario
MDOF Harmonic Response – crowd interacting with structure …IStructE Route 2 a = active crowd p = passive crowd s = structure
18
Tokyo Institute of Technology | Takeuchi Lab 18
Computing the Structural Response
foreach node // Get response of governing excitation frequency RNode = MAX(RNode(ƒ)) foreach excitation frequency // Get R Factor RNode(ƒ) = AccRMS,Node(ƒ) / AccRMS,0(ƒ) // Get SRSS acceleration AccRMS,Node(ƒ) = √(∑Mode∑HarmonicAccRMS(n,m,h,ƒ*h)2) foreach mode // Get modal properties
// modal frequency (ƒ), damping (ξ) & mass (m) // displacement at excitation(φne) & response(φnr) nodes // <participation factor(ρm)>,<static mass (W)>
ƒm,ξm,mm,φne,φnr,Wm,ρm = ... foreach harmonic // Get harmonic properties
// dynamic amplification factor(DLF), harmonic freq (ƒ)
DLFƒ,h,ƒm,h = ... // Get RMS acceleration AccRMS(n,m,h,ƒ*h)=RMS(F(ƒm,h,ξm,mm,φne,φnr,DLFƒ,h,Wm,ρm))
2000 nodes x 20 modes x 4 harmonics x 1.5Hz frequency range = 2.4E7 calculations …not practical back in ’90s, hence the simplified methods in AISC DG11, AS 5100-2, etc 6000 nodes x 40 modes x 2 harmonics x 6Hz frequency range (running) = 2.8E8 calculations …Excel VBA is single threaded so no parallel processing …Suggest building with parallel libraries in C#, VB, Python, etc
Footbridge GSA model Library floor GSA model
19
Tokyo Institute of Technology | Takeuchi Lab 19
20
Tokyo Institute of Technology | Takeuchi Lab 20
Examples
21
Tokyo Institute of Technology | Takeuchi Lab 21
Questions?