Predicting Footfall Induced Vibration Pa
Predicting Footfall Induced Vibration Pa
Predicting Footfall Induced Vibration Pa
Civil Engineers
Structures & Buildings 160
April 2007 Issue SB2
Pages 65–72
Paper 14648
Received 10/02/2006
Accepted 24/11/2006
Keywords: Michael Willford Peter Young Caroline Field
bridges/design methods & aims/ Director, Arup, London, UK Associate, Arup, London, UK Associate, Arup, San
dynamics Francisco, USA
The current paper reviews different methods that are 2. HISTORICAL BACKGROUND AND CURRENT
presently available for the prediction of footfall-induced INDUSTRY PRACTICE
vibration in structures. Five such methods are The issue of vibration serviceability for structures has been
considered in detail and it will be shown that these vary acknowledged for many years now. Early work by Smith 8 and
considerably in formulation, applicability and predicted others related to footbridges and this led to the performance
response. Where possible, the theoretical basis, criterion that still forms the basis of BS 5400. 4 Little
complexity, robustness, rigour, reliability and limitations consideration was given to the dynamic performance of floors
of each are considered. This paper is presented in two in buildings until steel composite design became more
parts. Part 1 deals with bridges and floors having vertical prevalent in the 1960s and 1970s. Since then, architectural and
natural frequencies below 10 Hz, part 2 deals with floors economic considerations have led to floors having longer spans
whose first vertical natural frequency is over 10 Hz. and lower weight, and the issue of vibration serviceability has
become much more significant.
Structures & Buildings 160 Issue SB2 Predicting footfall-induced vibration: Part 1 Willford et al. 65
greater than 4 Hz’) or the impulse response under a heel drop. root sum of the squares’ calculation. The force data used are
Neither approach represents a performance-based assessment of harmonic loads derived statistically from a large number of
the vibration likely to arise in normal service. The weight and measured footfall force–time histories.
stiffness of modern steel composite and post-tensioned
concrete floors are such that the need for a more robust and Reliance on fundamental dynamic principles makes the method
comprehensive assessment technique became clear. robust and applicable to any structure of arbitrary complexity
whose modal properties may be measured or calculated. This
In 1989 the Steel Construction Institute (SCI) published their rigour does, however, make the method more involved than the
Guide P076,1 which had been prepared independently by Dr others described here. For a structure with several modes that
Tom Wyatt of Imperial College London with the support of a might respond it is impractical to do the full calculation by
seven-person steering committee representing various sectors hand and a simple spreadsheet is required. CCIP-016 5 provides
of the industry. This presented the first widely recognised flowcharts to illustrate how this might be done.
performance-based approach for floors in buildings, including
calculation procedures aimed at predicting the actual 3. TECHNICAL BACKGROUND
vibration that would arise in practice as a result of walking. It
made the important distinction between lower-frequency 3.1. Quantifying response
floors that are susceptible to resonant response under repeated Absolute levels of footfall-induced vibration are normally
footfalls, and higher-frequency floors, which are not. quantified in terms of velocity or acceleration. It is, however,
Although applicable only to steel structures with a regular well documented (for example reference 13) that human
framing arrangement, its low-frequency floor prediction perception of vibration is influenced not just by amplitude but
methodology has proved to be reasonably robust and reliable also by frequency and duration. Amplitude alone is therefore
over 15 years. Its principal shortcomings are the very limited an inadequate descriptor of footfall-induced vibration.
experimental data upon which the pedestrian dynamic forces ‘Response factor’ (R) and ‘vibration dose value’ (VDV) are two
are based and the difficulty in applying the method to metrics that attempt to incorporate these additional factors.
anything other than the specific steel floor framing layouts
covered in the guide. BS 6472 13 contains curves which nominally define the
threshold of human perception to vibration. Specifically, these
In 2004 the SCI published P331,2 which purports to ‘enable define the perception threshold as a level of harmonic
designers to calculate the vibration response of sensitive floors vibration of constant amplitude that is a function of frequency
more accurately’. Unlike its 1989 predecessor, these guidelines and the orientation of the vibration to the individual. These
have not been independently reviewed and no data or curves are used to define a response factor of 1 (R ¼ 1), and the
theoretical basis that supports the methodology have been response factor of any harmonic response of known frequency
made generally available to facilitate such a review. Although and constant amplitude is found by dividing the amplitude by
it claims to be based in part on the work of the current authors, the corresponding base value at that frequency. The response
this paper shows that the P331 method can give very different factor is therefore a direct measure of perceptibility, R ¼ 2
(and often much less conservative) predictions than the Arup being twice what is perceptible and so on.
method. Further concerns relating to the theoretical basis of the
method have been discussed in recent correspondence between The VDV is also defined in BS 6472. It incorporates the
the current authors and the SCI in The Structural Engineer.9 duration of exposure to vibration and also weights higher
responses more heavily. To calculate a VDV it is necessary to
In the USA the most commonly used predictive technique for know how often and for how long a given level of vibration
low-frequency floors was published in 1997 by AISC. 3 Again, occurs. If vibration is attributed to a nearby railway and the
this guide is applicable to steel structures only and is similar in frequency of train pass-bys is known then this information
basis to SCI P076.1 It does provide a means of assessing the may well be available. In the context of footfall-induced
response in edge and corner bays (which is generally higher vibration, however, the required additional data (speed, number
than in central bays) and also considers the response of and route of pedestrians) are seldom available to the designer
footbridges. Otherwise its limitations are generally similar to and the VDV is less commonly used. In the present paper all
those of SCI P076. calculated responses will be expressed as response factors.
As a reaction to the inconsistencies and limitations of the 3.2. Resonant response of floors and bridges
published methods, Arup has developed its own methodology, It is commonly accepted by the methods reviewed here that
which is now used worldwide within the firm. It was first made ‘low-frequency’ structures are susceptible to resonant response
available outside of Arup in 2001. 10 The method has been under repeated footfall loading. Each methodology therefore
independently peer reviewed 5 and has recently been adopted seeks to predict the extent of this resonant response when
by other organisations and industrial bodies. 5,11 It is most assessing the maximum acceleration that is likely to occur in
completely described in the recent CCIP-016. 5 service. They differ, however, in the following respects
For low-frequency floors, standard modal response theory is (a) the threshold between ‘low’- and ‘high’-frequency floors
used to calculate the response to each harmonic load in exactly (i.e. the highest frequency at which a floor is susceptible to
the way that a finite element (FE) package such as Nastran 12 resonant build-up of vibration from successive footfalls)
calculates modal forced response. Total response is then (b) the magnitude of the harmonic components of footfall
calculated from the four harmonic responses using a ‘square forces
66 Structures & Buildings 160 Issue SB2 Predicting footfall-induced vibration: Part 1 Willford et al.
(c) how the resonant response is calculated for real structural they are well above or below a reasonable representative value.
systems with several modes and a limited number of Over most of the range there is a factor of at least two between
footfalls per span. the lowest and highest values selected, and at some frequencies
this difference is more than a factor of four. Even if the
The significance of these factors is discussed in turn below. dynamic properties of a floor and its response to a particular
force are known exactly, therefore, the choice of prediction
3.2.1. Threshold frequency for low-frequency floors. Most method will very significantly affect the predicted performance
modern floor vibration prediction methods consider there to be of the structure.
a cut-off frequency above which resonant build-up of vibration
response to repeated footfalls is not significant compared with As well as being very different in absolute terms, the trends in
the transient vibration generated by a single footfall. In harmonic force with frequency also differ greatly from one
justifying this, at least three approaches are possible. method to the next. Depending on the method adopted, an
engineer could be led to believe that by increasing the
First, measured response data can be analysed using a series of frequency of his structure from 5.5 Hz to 7.5 Hz there is either
band pass filters to isolate response in individual modes. The a 50% decrease (AISC) or a 20% increase (Arup) in force to
filtered time histories can then be inspected for evidence of consider. The implication in the AISC guide that someone
resonant build-up in each mode. Modes with frequencies walking very slowly at 1.5 steps per second applies a 40%
higher than the fundamental are always present in a floor. greater force (in the first harmonic) than someone who is
Therefore, if resonant excitation of the higher modes by higher walking very fast at 2.5 steps per second is counterintuitive
harmonics were possible, then this response should be and at variance with the measured data.
identifiable. In practice response in these modes is
characterised more by impulses followed by decays, and Ellis6 is alone in recommending that resonance be considered
resonant build-up over a number of cycles is not seen. for modes over 10 Hz. The conservatism associated with this
assumed loading becomes clear when comparing the predicted
Second, measured or simulated footfall force–time histories resonant response owing to a dynamic load factor (DLF) of 0.1
can be applied to a series of single mode systems of varying (as proposed by Ellis) with the measured data and impulsive
natural frequency so that the sensitivity of the maximum responses calculated for the worked examples in part 2 of the
dynamic response to exact synchronisation of a harmonic of current paper, 7 which deals with high-frequency structures.
footfall rate to a modal frequency can be examined. CCIP-016 5
presents the results of such an exercise. 3.2.3. Dynamic properties of floors and response calculation
method. The calculation of dynamic response requires some
Third, a sensible limit on the highest likely harmonic can be estimate of the natural frequency, effective mass and damping
estimated by considering the natural variability of the footfalls of a floor. The SCI P076, SCI P331 and AISC methods are based
of an individual. Steady state harmonic response analysis upon ‘hand calculation’ techniques that attempt to approximate
assumes a periodically repeating force input where every cycle the true situation (in which there are usually several or many
is identical. The natural variation between left and right foot modes of lowish natural frequency) by a single equivalent
and perhaps a 10% (0.04 s) variation in exact pacing rate mode. This inevitably introduces a degree of approximation,
between steps can be used to show that effective harmonic which may result in significant error, as is illustrated in the
excitation beyond the fourth harmonic (10 Hz, 0.1 s period) worked example presented later in the current paper.
will be very low.
A more rigorous approach is to derive modal properties of a
The assumed cut-off frequency for each of the methods is shown structure by finite element analysis (FEA), which is the basis of
in Fig. 2 (see below). The lowest recommendation is 5 Hz in BS the Arup method. This approach enables a much more reliable
5400, 4 the highest, proposed by Ellis,6 is at least 19.2 Hz. It is assessment to be made, including the effect of simultaneous
the current authors’ experience that the practical significance of response in the many closely spaced modes that will be present
resonance above the fourth harmonic is very small. in large floor plates. It also enables the response to be
calculated at points distant from the walking route and the
3.2.2. Harmonic forces. The harmonic components of footfall effects of partitions to be included explicitly. Full-height
forces from over 800 measurements are plotted in Fig. 1. It can partitions can significantly increase the stiffness and damping
be seen that there is a very substantial scatter in these values, of floors, but this benefit must be evaluated on a case by case
and it is therefore important to realise that there is no single basis. It should also only be included in a response calculation
‘correct’ value, but rather that values can be selected within the if the associated compromise in space planning flexibility is
range to be more or less conservative. Ideally, therefore, the agreed with the client.
values selected should have a probability of exceedance
assigned to them. Most of the published methods select values As mentioned in section 3.2 above, the generally accepted
based on a much smaller data set and with no statistical basis. design process for low-frequency floors should consider the
The Arup method selects values with a 25% probability of maximum level of resonant response that can be induced by a
exceedance based on this large data set. person walking on the structure. There is, however,
considerable variation between the methods in how the
The harmonic force components that are assumed by each magnitude of this response should be calculated. The following
method are plotted in Fig. 2. It can be seen that while for the factors need to be considered in any reliable form of
most part they fall within the measured data set, sometimes assessment
Structures & Buildings 160 Issue SB2 Predicting footfall-induced vibration: Part 1 Willford et al. 67
1st harmonic 2nd harmonic
1·0 0·25
0·9
0·8 0·20
0·7
0·6 0·15
DLF
DLF
0·5
0·4 0·10
0·3
0·2 0·05
0·1
0·0 0·00
0·0 0·5 1·0 1·5 2·0 2·5 3·0 3·5 0·0 1·0 2·0 3·0 4·0 5·0 6·0 7·0
Frequency: Hz Frequency: Hz
0·16 0·14
0·14
0·12
0·12
0·10
0·10
DLF
DLF 0·08
0·08
0·06
0·06
0·04
0·04
0·02 0·02
0 0
0·0 2·0 4·0 6·0 8·0 10·0 0·0 2·0 4·0 6·0 8·0 10·0 12·0 14·0
Frequency: Hz Frequency: Hz
Kerr Galbraith and Barton Wheeler Ohlssen Rainer Ellis Mean Value Design Value
0·3 Ellis
BS 5400 (bridges)
0·2 4. REVIEW OF RESPONSE CALCULATION
19·2 Hz PROCEDURES: WORKED EXAMPLES
0·1
The inconsistency of the various methods is illustrated by
0 example below. Unless otherwise specified, a 750 N person
0 2 4 6 8 10 12 walking at between 90 and 150 paces per minute (1.5 to
Forcing frequency: Hz
2.5 Hz) with a third harmonic force of 75 N is assumed.
Fig. 2. DLFs for various footfall prediction methodologies Damping of 3% of critical damping has been used throughout.
68 Structures & Buildings 160 Issue SB2 Predicting footfall-induced vibration: Part 1 Willford et al.
Grid 7.5 m 3 7.5 m
Beam layout: 300ASB(FE)153 (asymmetric beam) secondaries spanning in one direction between columns
Decking Comdek SD225, 1.25 mm thick spanning 7.5 m between secondaries
Concrete Normal weight with overall slab depth of 316 mm
Floor size 16 bays in the concrete spanning direction, 2 to 16 bays in secondary spanning direction (the SCI P331
assumes 16 bays in its worked example)
(a) The ‘FE’ solution is a forced response calculation assuming Slimdek floor of varying size,
a 75 N harmonic force applied in the middle of a bay using 16 bays in beam spanning direction SCI P331
8 SCI P076
7 Arup 750N walker
B AISC
300 ASB153 6
Response factor
FE
A A 5
22·5 m
3 bays 4
3
B 2
120 m, 16 bays
1
Section A–A 91 mm
0
0 2 4 6 8 10 12 14 16 18
316 mm Number of bays in beam spanning direction
300
mm 7·5 mm Fig. 5. Predicted response factors for a Slimdek floor of
varying number of bays
91 mm
Section B–B
225 mm all modes up to 15 Hz. For a given set of modes this will
calculate the theoretically correct, full resonant response
400 200
mm mm
for this load using standard modal response theory.
Resonant response to a harmonic load of 75 N is the basis
of the SCI P331 method and so the FE and SCI P331
Fig. 3. Slimdek floor. Three bays in the spanning direction are
calculations would be expected to give very similar results.
shown
(b) The Arup method uses the mode-shapes from the FE model
Structures & Buildings 160 Issue SB2 Predicting footfall-induced vibration: Part 1 Willford et al. 69
(modes up to 15 Hz) to calculate the response using the different construction types. This is a considerable weakness of
method described in reference 5. both the SCI and AISC techniques, which makes like for like
(c) All calculations are for mid-span in a bay in the middle of comparisons very difficult.
the floor. (The AISC and Arup methods allow calculation of
response in edge and corner bays, where the response will A comparable concrete design to the Slimdek floor of example
be higher. The FE method could similarly calculate 1 is a 220 mm post-tensioned flat slab, and the first three
response away from the middle of the floor.) modes of the 16 33 bay floor predicted by FEA are shown in
Fig. 6. Each mode has a modal mass of approximately 430 t.
The following points arise from this graph. Data on the floor summary are presented in Table 2.
70 Structures & Buildings 160 Issue SB2 Predicting footfall-induced vibration: Part 1 Willford et al.
5·0 considered here only the Arup method statistically
FE–75 N harmonic force
4·5 quantifies this variation or indicates the statistical
Arup 750 N Walker
4·0
significance of the chosen loads.
(d ) The AISC method incorporates forces that are clearly
3·5
Response factor
Structures & Buildings 160 Issue SB2 Predicting footfall-induced vibration: Part 1 Willford et al. 71
methodologies. Arup Vibration Seminar, Institute of 13. BRITISH STANDARDS INSTITUTION. Guide to Evaluation of
Mechanical Engineers, London, 2001. Human Exposure to Vibration in Buildings (1 Hz to 80 Hz).
11. THE CONCRETE SOCIETY. Post-tensioned Concrete Floors – BSI, London, 1992, BS 6472.
Design Handbook. The Concrete Society Camberley, 2005, 14. NHS ESTATES. Acoustics, Design Considerations. Health
Technical report 43. Technical Memorandum 2045. Her Majesty’s Stationery
12. MSC Nastran. See www.mscsoftware.com. Office, Norwich, 1996.
72 Structures & Buildings 160 Issue SB2 Predicting footfall-induced vibration: Part 1 Willford et al.