Skip to main content

Footfall Analysis

Footfall analysis (or in full, footfall induced vibration analysis) is used to calculate the elastic vertical nodal responses (acceleration, velocity, response factor etc.) of structures to human footfall loads (excitations). The human footfall loads f(t)f(t) are taken as periodic loads. Using to Fourier Series, the period footfall loads can be expressed as:

f(t)=G[1+h=1Hrhsin(2πhTt)]f(t) = G\left\lbrack 1 + \sum_{h = 1}^{H}{r_{h}\sin\left( \frac{2\pi h}{T}t \right)} \right\rbrack

where GG is the body weight of the individual, and rhr_{h} are the Fourier coefficients (or dynamic load factor), the actual values of dynamic load factors can be found from reference 24, 25 and 26 in the bibliography, TT is the period of the footfall (inverse of walking frequency) and HH the number of Fourier (harmonic) terms to be considered, 4 is used for walking on floor using CCIP-016 method, 3 is used for walking on floor using SCI method and 2 is used for walking on stairs.

After subtracting the static weight of the individual (since it does not vary with time and does not induce any dynamic response), the dynamic part of the footfall loads are the sum of a number of harmonic loads

f(t)=Gh=1Hrhsin(2πhTt)f(t) = G\sum_{h = 1}^{H}{r_{h}\sin\left( \frac{2\pi h}{T}t \right)}

There are two distinctive responses from the footfall excitation, the resonant (steady state) and transient. If the minimum natural frequency of a structure is higher than 4 times the highest walking frequency (see reference 24), the resonant response is normally not excited since the natural frequencies of the structure are so far from the walking (excitation) frequency, therefore the transient response is normally in control, otherwise, the resonant response is probably in control. Both resonant (steady state) and transient analyses are considered in GSA footfall analysis, so the maximum responses will always be captured.

Resonant response analysis

As footfall loads are composed of a number of harmonic loads (components), harmonic analysis is used to get the responses for each of the harmonic components of footfall loads and then to combined them to get the total responses. From one of the harmonic components (h)(h) of the footfall loads in equation above and the given walking frequency (1T)\left( \frac{1}{T} \right), the following dynamic equation of motion can be obtained

Mu¨h+Cu˙h+Kuh=δkGrhsin(2πhTt)\mathbf{M}{\ddot{\mathbf{u}}}_{h} + \mathbf{C}{\dot{\mathbf{u}}}_{h} + \mathbf{K}\mathbf{u}_{h} = \delta_{k}Gr_{h}\sin\left( \frac{2\pi h}{T}t \right)

where δk\delta_{k} is a unit vector used to define the location of the harmonic load. All the components in this vector are zero except the term that corresponds to the vertical direction of the node subjected footfall load.

Since the number of footfalls is limited and the full resonant response from the equation above may not always be achieved, a reduction factor ρh,m\rho_{h,m} for the dynamic magnification factors μ\mu is needed to account for this non-full resonant response. The reduction factor can be calculated from

ρh,m=1e2πζmN\rho_{h,m} = 1 - e^{- 2\pi\zeta_{m}N}

Where ζm\zeta_{m} is damping ratio of mode m and N=0.55hWN = 0.55hW with hh the harmonic load number and WW the number of footfalls.

Applying this reduction factor to the dynamic magnification factors (μ)(\mu) in Harmonic Analysis, this equation can be solved using the method described in Harmonic Analysis Theory section. Repeating this analysis, the responses from the other harmonic loads of the footfall can also be obtained. The interested results from this analysis are the total vertical acceleration and the response factor from all harmonic loads of the footfall. The total vertical acceleration is taken as the square root of the sum of squares of the accelerations from each of the harmonic analyses. The response factor for each of the harmonic loads is the ratio of the nodal acceleration to the base curve of the Root Mean Square acceleration given in reference 25 as shown below. This total response factor is then taken as the square root of the sum of squares of the response factors from each of the harmonic loads. According to this, the total acceleration and response factor can be calculated from

ai=h=1Hu¨i,h2fi=h=1Hfi,h2=h=1H(u¨i,h2wi0.0052)\begin{aligned}a_{i} &= \sqrt{\sum_{h = 1}^{H}{{\ddot{u}}_{i,h}}^{2}} \\ f_{i} &= \sqrt{\sum_{h = 1}^{H}{f_{i,h}}^{2}} = \sqrt{\sum_{h = 1}^{H}\left( \frac{{{\ddot{u}}_{i,h}}^{2}w_{i}}{0.005\sqrt{2}} \right)}\end{aligned}

where

aia_{i} is the maximum acceleration at node ii

u¨i,h{\ddot{u}}_{i,h} is the maximum acceleration at node i by the excitation of harmonic loadhh

HH is the number of harmonic components of the footfall loads considered in the analysis

fif_{i} is the response factor at node ii

fi,hf_{i,h} is the response factor at node ii by the excitation of harmonic loadhh

wiw_{i} is the frequency weighting factor and it is a function of frequency

For standard weighting factors see Table 3 of BS6841.

Transient response analysis

The transient response of structures to footfall forces is characterised by an initial peak velocity followed by a decaying vibration at the natural frequency of the structure until the next footfall. As the natural frequencies of the structure considered in this analysis is much higher the highest walking frequency, there is no tendency for the response to build up over time as it does in resonant response analysis. The maximum response will be at the beginning of each footfall. Each footfall is considered as an impulse to the structure, according to references 35 & 29, the design impulse can be calculated from

When walking on floor (Concrete Centre/Arup method)

Ides,m=54f1.43fm1.3I_{des,m} = 54\frac{f^{1.43}}{{f_{m}}^{1.3}}

When walking on floor (SCI P354 method)

Ides,m=60f1.43fm1.3Q700I_{des,m} = 60\frac{f^{1.43}}{{f_{m}}^{1.3}}\frac{Q}{700}

When walking on stairs (Concrete Centre/Arup method)

Ides,m=54150fmI_{des,m} = 54\frac{150}{f_{m}}

When walking on stairs (SCI P354 method)

Ides,m=0I_{des,m} = 0

where

Ides,mI_{des,m} is the design impulse for mode mmin NS

ffis the walking frequency in Hz

fmf_{m}is the natural frequency of the structure in mode mmin Hz

QQis the weight of the walker in N

For this impulse, the peak velocity in each mode is given by

v^m=ue,mur,mIdes,mm^m{\widehat{v}}_{m} = u_{e,m}u_{r,m}\frac{I_{des,m}}{{\widehat{m}}_{m}}

and the peak acceleration in each mode is given by

a^m=2πfmv^m=2πfmue,mur,mIdes,mm^m{\widehat{a}}_{m} = 2\pi f_{m}{\widehat{v}}_{m} = 2\pi f_{m}u_{e,m}u_{r,m}\frac{I_{des,m}}{{\widehat{m}}_{m}}

where

v^m{\widehat{v}}_{m} is the peak velocity in mode mm by the footfall impulse

a^m{\widehat{a}}_{m} is the peak acceleration in mode mm by the footfall impulse

ue,m,ur,mu_{e,m},u_{r,m} are the vertical displacements at the excitation and response nodes respectively in mode mm

m^m{\widehat{m}}_{m} is the modal mass in mode mm

The variation of the velocity with time of each mode is given by

vm(t)=v^me2πξfmtsin(2πfmt)v_{m}(t) = {\widehat{v}}_{m}e^{- 2\pi\xi f_{m}t}\sin\left( 2\pi f_{m}t \right)

and the variation of the acceleration with time of each mode is given:

am(t)=a^me2πξfmtsin(2πfmt)a_{m}(t) = {\widehat{a}}_{m}e^{- 2\pi\xi f_{m}t}\sin\left( 2\pi f_{m}t \right)

where ξ\xi is the damping ratio associated with mode mm

The final velocity and acceleration at the response node are the sum of the velocities and accelerations of all the modes (M)\left(M\right) that are considered

v(t)=m=1Mvm(t)a(t)=m=1Mam(t)\begin{aligned}v(t) &= \sum_{m = 1}^{M}{v_{m}(t)} \\ a(t) &= \sum_{m = 1}^{M}{a_{m}(t)}\end{aligned}

This gives the peak velocity and peak acceleration. The root mean square velocity and root mean square acceleration can be calculated from the period of the footfall

vRMS=1T0Tv2(t)dtaRMS=1T0Ta2(t)dt\begin{aligned}v_{RMS} = \sqrt{\frac{1}{T}\int_{0}^{T}{v^{2}(t)dt}} \\ a_{RMS} = \sqrt{\frac{1}{T}\int_{0}^{T}{a^{2}(t)dt}}\end{aligned}

The response factor at time tt (tt is from 00 to TT and TT is the period of the footfall loads) can be calculated from

fR(t)=10.005m=1Mam(t)wmf_{R}(t) = \frac{1}{0.005}\sum_{m = 1}^{M}{a_{m}(t)w_{m}}

where

wmw_{m} is the frequency weighting factor corresponding to the frequency of mode mm

The final transient response factor, based on the root mean square principle, is given by

fR=1T0TfR2(t)dtf_{R} = \sqrt{\frac{1}{T}\int_{0}^{T}{{f_{R}}^{2}(t)dt}}