Reference

NBS

Publications

NBSIR 75 - 951

Correlation of Floor Vibration

to Human Response

James R. Shaver

Center for Building Technology

Institute for Applied Technology

National Bureau of Standards

Washington, D. C. 20234

December 1975

Final Report

Prepared for Office of Policy Development and Research

Department of Housing and Urban Development

Washington, D. C. 20410

75-951

1975

NBSIR 75 - 951

National Bureau of Standards

Library, E-01 Admin. Bldg.

OCT '4 5 1982

CORRELATION OF FLOOR VIBRATION

James R. Shaver

Center for Building Technology

Institute for Applied Technology

National Bureau of Standards

Washington, D. C. 20234

December 1 975

Final Report

Prepared for Office of Policy Development and Research

Department of Housing and Urban Development

Washington, D. C. 20410

U.S. DEPARTMENT OF COMMERCE, Rogers C.B. Morton, Secretary

James A. Baker, III, Under Secretary

Dr. Betsy Ancker-Johnson, Assistant Secretary for Science and Technology

NATIONAL BUREAU OF STANDARDS, Ernest Ambler, Acting Director

TO HUMAN RESPONSE

Cic IOO

ABSTRACT

A new approach to the problem of perceptible floor vibrations

is presented based on the assumption that human activity, which

produces the vibration, and human response to vibration are random variables. Techniques for data reduction are discussed

and a detailed description of one approach is given along with the

associated computer programs. Data from floor vibrations is

compared with current criteria for human response.

Key words: Analysis; experimental; floor systems; human response;

random process; spectral analysis; vibration.
TABLE OF CONTENTS

Page

Abstract • ii

1 . Introduction 1

2 . Background 2

3. Approach to Problem 4

4. Data Reduction and Analysis 6

5. Damping 10

6. Human Response to Vibration 11

7. Summary 14

8. References 15

9. Appendix A: Field Test 16

10. Appendix B: Computer Program •. . . . 18

iii
CORRELATION OF FLOOR VIBRATION

TO HUMAN RESPONSE

BY

James R. Shaver

1. INTRODUCTION

Historically, the design of floor systems for serviceability has been

based on limiting the deflection that can occur for a given design load to

prevent structural damage. It was believed by most engineers that this design

procedure indirectly accounted for the effect of human activity on floor

systems which could produce quite perceptible and annoying vibrations for other

human occupants. Research and state-of-the-art studies [1,2,3]- conducted

during the past decade on the various aspects of vibrating floor systems have

increased the awareness of the engineering community to the fact that the traditional

serviceability design procedure does not necessarily lead to floor systems with

satisfactory dynamic characteristics which prevents objectionable vibrations

from the occupant view point when the floor is subjected to human activity.

Recent efforts [4.5] to study the problem of perceptible floor vibrations

have primarily employed experimental techniques. Results from these studies

have been inconclusive as to what constitutes objectionable floor vibration

from an occupant viewpoint. Questions have been raised about the proper

forcing function to be used in the test to simulate human activity [2,5] and

the lack of information on the human response to this type of vibration [4]

.

The experimental procedures used have been based on the assumption that both

human activity which produces vibration and occupant response to vibration

of a floor system are deterministic variables. However most scientific studies

which deal with human activity and human response as variables treat them as

a random variable rather than a deterministic one.

This report describes a new approach to the problem of perceptible floor

vibrations based on the assumption that human activity and human response to

2/ vibration are random variables. Presented is a random variable methodologybased on physical testing and subsequent analysis of the dynamic response of

floor systems subjected to human activity.

—^Numbers in brackets refer to literature references. —^This methodology was developed as part of a larger research effort into the

nature of structural deflection limitations sponsored by the Department of

Housing and Urban Development.

1
2 . BACKGROUND

In order to understand why both human activity and occupant response

should be characterized as random variables, it is necessary to understand

the difference between a deterministic variable and one that is random. A

deterministic variable is one whose value can be calculated or determined

from knowledge of its behavior at any previous instant in time, while a random variable is one whose instantaneous value cannot be predicted with

complete certainty for any given instant in time. This means that the

instantaneous value for a random variable can only be specified by a probability distribution function which gives the probable fraction of the total

time that the instantaneous value lies within a specified range.

If one considers all the possible parametric variations that can affect

the instantaneous value of the dynamic loading that occurs from human

activities such as walking, running, jumping, it becomes clear that each

activity is, in itself, a random process composed of many variables. For

example, considering footfall or walking since it is the most prevalent

activity that occurs on a floor system, variations in weight, gait, heel-toball of foot contact and foot wear of individuals will all affect the instantaneous value of the dynamic loading that occurs.

It appears, from an examination of the literature, that no definitive

study has been made to determine the statistical data needed for characterization of footfall as a random variable. However, a 1961 English floor

abrasion study [4] made some quantitative measurement of the forces applied

to the floor by walking and does give some insight into the nature of footfall from a static viewpoint. A typical plot of the vertical force component

for a footfall, obtained in this 1967 study, appears in Figure 1. Figure 1 shows two peaks and a trough. As pointed out by the authors, the first and

second peaks will vary in time of occurrence from initial contact and there

is often, although not always, a discontinuity in the curve as it rises to

the first peak. The study also indicated that if the vertical force component was considered in relation to the weight of the person, it appeared

that the trough in the curve is always less than body weight while the

peaks may be either above or below depending on the manner in which the

body trunk is accelerated. Additional information is needed on the duration

of footfall, its frequency content and probability of peak values in order

to describe it as a random dynamic forcing function.

2
The human response to vibration involves many variables: vibration

input where intensity, frequency, direction and duration must be considered;

psychological influences in the form of mental state, motivation and experience; and physical influences from sound and sight. Given this multiplicity

of parameters, it is reasonable to classify human response as a random

process

.

It should be pointed out that the structural floor systems are not

considered to be a random variable. Although there are several types of

floor systems, their dynamic properties can be described deterministically

from an engineering viewpoint. Rather it is the forcing function (human

activity) and the human response to this activity which must be considered

as random processes. This makes the overall problem of human response to

floor vibration induced by human activity a random vibration problem. This

can be stated in functional form as:

F(t) + G(t) -> Y(t) D(t)

where F(t) is the random forcing function, G(t) is the deterministic floor

system characteristics, Y(t) is the random response of the floor system and

D(t) is the random human response to the floor vibration, Y(t)

.

3
3. APPROACH TO PROBLEM

The dynamic response of a floor system subjected to human activity

represents a random physical phenomenon which cannot be described by an explicit mathematical relationship because any given observation will

represent only one of the many possible results which might have occurred.

A single time history representation or observation of the random phenomenon

is called a sample function (a sample record when observed over a finite

time interval) . Then the collection of all possible sample functions which

the random phenomenon might have produced is called a random process. From

a given sample function it is possible to estimate statistical parameters

necessary to describe the random process.

Four main types of statistical functions are used to describe the

basic properties of random data: (a) mean square values, (b) probability

density functions, (c) autocorrelation functions, and (d) power spectral

density functions. The mean square value furnishes a rudimentary description

of the intensity of the data. The probability density function furnishes

information concerning the properties of the data in the amplitude domain.

The autocorrelation function and the power spectral density function furnish

similar information in the time domain and frequency domain, respectively.

While it is possible to determine the required parameters directly

from the time history under certain simplifying assumptions, the majority

of the techniques used to date perform the analysis in the frequency domain.

Analysis in this domain is usually termed spectral decomposition.

There are three reasons for selecting spectral decomposition over time

history to represent one's data. First, spectral decomposition . presents

data in a form that is related to quantities familiar to engineers. Many

physical systems (including the human body) have been modeled with a fair

degree of success in this domain. Most models, for example, consist of a set of simple second order oscillators connected in a manner which simulates

the various components of the particular physical systems. Such a model is

shown in Figure 2. Each oscillator has its own oscillatory or resonant

frequency which, in turn, is dependent upon the size, shape and material

of the component it is simulating. Since each oscillator will be

excited primarily by energy at its resonant frequency, the knowledge of the

frequency content of the response provides the key to the determination of

the system response to impulse or shock, such as that produced by footfall.

Generally, the response of a floor system at resonant frequencies will

predominate in the total response. The second reason for utilizing spectral

4
/

decomposition is the independence of each frequency component from all

others. This independence in the frequency domain greatly simplifies

analysis. By comparison, the determination of a response in the time

domain at some instant in time is dependent not only on the value of the

excitation at that time but on all previous values of the excitation.

Third, this type of analysis can be used to describe input data, F(f) , response data, Y(f), or, if the input and response are measured simultaneously,

it can be used to describe the frequency response function, G(f), of the

system as shown in Eq. 3.1. Specifically, given any two of the above

items, the third can be determined.

Y(f ) = |G(f) | 2 F(f ) (3.1)

The primary approach to obtaining a meaningful dynamic response of a floor system subjected to human activity is to first obtain a time history

and then transform this time history into the frequency domain using Fourier

techniques. In addition, the mean square value and probability density

function should be obtained. It does not appear at this time that the

autocorrelation function or the power spectral density function will provide

any useful information.

5
4. DATA REDUCTION AND ANALYSIS

Spectral decomposition or analysis using Fourier techniques consists

of decomposing a time history in terms of trigonmetric functions. The

magnitude and phase of each trigonometric term needed to reconstruct the

time history are plotted as a function of frequency. This technique is

simply the evaluation of the finite, Fourier transform of the time history

as expressed by Equation 4.1 where the time history x(t) is multiplied by

the Fourier kernel, a complex exponential containing both the frequency and

time variables, and then integrating this product over the record length;

T/2 f

X(f) = / xUJe'^^cit (4.1)

-T/2

Using Euler's transformation, the kernel may be expanded so that Eq.

4.1 may be rewritten as,

T/2 T/2

X(f) = / x(t) cos 27Tftdt - jf X(t) sin 2nftdt (4.2)

-T/2 -T/2

This equation shows how spectral analysis can be considered as a decomposition

of the time history in terms of sine and cosine components.

In order to evaluate this integral by digital techniques, one must

solve the discrete, finite Fourier transform which is given as

X(kAf) = At^V X. e -j 2 ^ kAfiAt

(4>3)

i=o

where X^ is the digital representation of the time history obtained by

sampling the continuous record at At intervals and the natural frequency

spacing for the complete transform is

Af = — (4.4)

NAt

Then kAf is defined by

kAf = k = 0, 1, 2..., N/2 (4.5)

Equation 4.3 may then be rewritten as N-l -j2irik

X = At Z X. e N k = 0, 1, 2, . ., N/2 (4.6)

K 1 i=o

6
Employing Euler ' s formula again and rewriting, Eq. 4.6 becomes

X

k

= A

k

- jBk k = 0, 1, 2, N/2 (4.7)

where N-l

: I i=oA

R

= At I X

i cos (4.7a

N-l

: Z i=oB

k

= At I X

± sin i^ik (4.7b)

A computer program can be written which will evaluate Eq. 4.7 directly.

However, this program would require a rather large amount of running time

particularly when large volumes of data are being analyzed. Although the

most efficient techniques known to date for evaluating Eq. 4.6 are the Fast

Fourier Transform (FFT) methods, they were not used because these methods

are not straight forward in their application and do require some special knowledge and experience for proper utilization.

The evaluation procedure employed to evaluate Eq . 4.7, was first

derived by Goertzel [5] and is believed to be the most efficient of the

pre-FFT methods. - The method requires the generation of an auxiliary variable

Q as follows:

Q

Q

= 0

Q

l = Xn-1

Q

i = (2 cos ^ ) Q

i_1 -Q i _ 2 + Xn -i, i = 2, 2,. . . , N-l (4.8)

Then, A

k and B

k

in Eqs. 4.7a and 4.7b become

A

k

= At [(COS TT} QN-1

"QN-2

+ X

o

] k =

°' l ' • "

N/2

B

k

= At (sin Q

n

_ 1 (4.9)

This technique was utilized in the spectral analysis of the 40 second walking

record obtained using the test procedure given in Appendix A and required

less than 3 seconds for 12,000 data points on a UNIVAC 1108 with an EXEC 8 system. The computer program is given in Appendix B. It should be remembered that the evaluation of Xk by Eqs. 4.8 and 4.9 would require

7
generally unacceptable amounts of computer running time when extremely

large volumes of data are to be analyzed.

The most useful output from the Fourier spectrum is the modulus. This

modulus is obtained by realizing that the Fourier spectrum can be represented

by a complex expression.

X(f) = |x(f)| e-^ 0(f) (4 ' 10 >

where [ X ( f ) ] is the modulus and e 3®^' i s the phase.

From Eq. 4.7, the expression for the finite Fourier transform, the

modulus is,

|Xk

| = (A

2 + b£) 1/2 (4.11)

The modulus dimensionally represents a pseudovelocity and it can be

shown [6] that this is identical to the relative velocity of the residual

response for a linear undamped single degree-of -freedom system when the

system response is a pure sinusoid. In general, the two will be in close

agreement for any undamped system with the pseudovelocity lower at low

frequencies and higher at high frequencies than the relative velocity

because of its frequency dependence.

The phase presentation of the spectrum may prove to be useful in

determining damping in the system. Again from Eq. 4.7, the phase angle is

given by

-1 B

k

0 = tan 1 ~ (4.12) A

k

Figure 3 is a plot of the modulus of the spectral decomposition of

the 40 second walking record. The values of this modulus are in mq-s which

is a pseudovelocity in feet per second. This graph is also called an energy

density curve in random data analysis. Due to filters in the acceleration

signal conditioning equipment all input above 25 Hz has been removed thus

accounting for the low energy content above this value.

8
The acceleration record was sampled every 160 ms which is

equivalent, because of a difference in record and playback speed of the

analog data tape, to a 10 ms sampling rate in real time. At the 10 ms real

time sampling rate, 1191 samples were taken thus producing a natural frequency

spacing of 0.084 Hz in the Fourier transform. The modulus of the spectral

decomposition has been plotted at every 0.5 Hz.

An examination of Figure 3 indicates that significant amounts of

energy occurs in the range of 15-25 Hz with a pronounced peak at 18 Hz

which is the fundamental frequency for the floor tested. The existence of

energy in this range indicates that a vibrating floor system is capable of

inducing a human response to vibration at frequencies other than the

lowest natural frequency of the system. Thus, test procedures which determine

the vibration level only at the fundamental frequency for the floor system

may exclude significant energy levels at other frequencies.

The missing link at this time is a means of transforming the energy

level at these frequencies to a corresponding root-mean-square (rms) acceleration at the same frequency thereby permitting a comparison with the current

human response criteria. However, it is possible to determine the rms acceleration for the record and associate it with the lowest natural frequency

for comparison purposes. This is done in section 6.

9
5. DAMPING

Damping in floor systems is a very complex quantity yet its measurement

is important since it effects the level and duration of vibration. In

addition, it has been suggested that in the case of a single pulse input

the occupants will perceive motion at a lower threshold if the vibration of

the floor does not decay rapidly. Although the scope of the study did not

permit an investigation into more effective means of determining damping

and relating its effect to human response a brief discussion will be given

because of the significant role damping plays in any dynamic problem.

Structural damping that occurs in floor systems is composed of three

primary components: (1) material damping which is the ability of the

materials themselves to dissipate energy, (2) interlayer slip which is

energy dissipation through friction between the subfloor and joist and (3)

end conditions which may damp motion by inducing support vibration or dissipate energy through friction. Generally, the energy loss from material

damping is small as compared to the amount dissipated by slip and support

effects.

There are six generally accepted mathematical means of representing

damping: (1) logarithmic decrement, (2) amplification factor, (3) quality

factor, (4) equivalent dashpot constant, (5) complex modulus and, (6)

bandwidth. For the simple case of a linear single degree of freedom

system, a relationship can be established between all six definitions, but

for other cases this is no longer necessarily true. When structural damping

is of primary interest any one of the six representations could be used

with some degree of success depending upon the end result desired. Logarithmic decrement or decay rate is the simplest to measure experimentally but

is difficult to use in design or analysis. The most useful form of presentation for the designer is specific damping energy, in inch- pounds per

cycle, while for purposes of analysis the equivalent dashpot constant is

most useful.

However, from the viewpoint of the present approach taken in this

methodology and its emphasis on spectral techniques, the use of either the

bandwidth method or quality factor in characterizing the damping of structural

floor systems should lead to a better means of describing the damping

quality of floors as it relates to human response to vibration. While it

appears feasible at this time to obtain these damping characteristics from

a spectral decomposition, this derivation is beyond the scope of this

report.

10
6. HUMAN RESPONSE TO VIBRATION

Human body vibration is usually treated as an externally-applied

condition, although it is often self-induced such as during walking or running. The spectrum of human reaction to vibration is extremely wide

encompassing both "pleasant" effects which range from relaxiation to stimulation and "unpleasant" effects which span from discomfort to injury. Much

of the information and literature concerning the effects of vibration on man is of a provisional, specialized or contradictory nature. This is

particularly true in the area of human response to floor vibrations where

only in the past decade have attempts been made to quantify this response.

Thus, there does not exist, at this time, a consensus standard for determining whether or not floor vibration induced by human activity is acceptable

from the occupant viewpoint.

This section will present the current state of knowledge on human

reaction to externally applied vertical transient vibration against which

the results of the proposed floor vibration methodology already given might

be compared.

The current state of knowledge on human reaction to externally-applied

vertical transient vibration is based on a limited number of studies and

represents a wide range of thresholds. However, as basic as it is, this

knowledge is currently the only information on which to base a judgement of

whether or not a floor is acceptable from a human response viewpoint.

The ISO Guide for the Evaluation of Human Exposure to Whole-Body

Vibration ISO/DIS 2631 [7] provides allowable rms accelerations at different

frequencies for various levels of exposure time. It is primarily intended for

periodic (sinusoidal) vibrations in the frequency range from 1 Hz to 80 Hz but

does cover narrow-band or wide-band vibration within this frequency range.

In the case of either wide- or narrow-band random vibration, the allowable

accelerations are based on the limit given for the center frequency in each

third octave band.

Splittgerber [8] has proposed a standard based on the ISO Guide which

provides vibration and shock limits for occupants of buildings. This

proposed standard differentiates between continuous or intermittent vibration

and impulsive shock excitation during either a day or night environment.

It also contains vibration limits for different types of occupancy such as hospital, residential, office, etc.

11
A standard acceleration level is specified for each direction of

vibration and then permissable acceleration levels are obtained for the

different types of vibration input, occupant classification, and time

conditions by applying an appropriate weighting factor to the standard.

Figure 4 shows the ISO Guide's standard acceleration curve for vertical

vibrations and Splittgerber 1 s permissable acceleration level for intermittent vibration in residences during the daytime hours. This daytime

level is based on minimum complaints from the occupants and was obtained by

applying a weighting factor of two to the standard curve.

Plotted on Figure 4 for purposes of comparison is Lenzen's [9] modification of the original Reiher and Meister work [10] on steady state vibrations

for transient vibration. Lenzen's modification considers transient vibration

as a vibration which may last only a few seconds but is characterized by a buildup to a level which is maintained for several cycles. Realistically,

Lenzen's definition fits the vibration produced by a single footfall or someone jumping with Splittgerber 1 s definition fitting the vibration produced

by walking or running. The dashed lines on Figure 4 denote the slightly

perceptible range for Lenzen's modification which is again believed to

represent a minimum complaint level.

Also plotted in Figure 4 is Splittgerber ' s rms acceleration limit for

daytime impulsive shock in residences which he characterizes as a rapid

buildup to a peak followed by decay. It is comparable to Lenzen's definition

for transient vibration. This acceleration limit was obtained by weighting

the standard curve by sixteen. As can be seen from Figure 4 there is good

agreement between these two criteria for f > 8 Hz with a divergence of the

two below this frequency. No definite conclusion can be drawn at this time

on the reason for this divergence in the low frequency range.

Plotted as a triangle on Figure 4 is the peak rms acceleration of 0.95

mg as determined from the 40 second continuous walkinq test record. This value

was obtained using analog techniques with the time constant of the rms meter

and scan rate of the data adjusted so that the uncertainty is less

than 2% and is the peak rms value associated with maximum acceleration.

Since the major energy component in the frequency spectrum occurs at 18 Hz,

in accordance with the ISO Guide, this rms value is plotted at a frequency

of 16 Hz as this is the center frequency for the third octave band in which

18 Hz falls.

12
A reinterpretation of this rms acceleration with regard to frequency

and type of excitation was also made in order to compare the results with

Lenzen's criteria. Lenzen's criteria would consider it to be a single

pulse at the predominate floor frequency which in this case is 18 Hz. The

maximum rms acceleration is found from the expression:

rms peak

and is 1.2 mg. This value is plotted as a circle on the figure.

From the figure it can be seen that the treatment of the walking

record as an intermittent vibration gives a resulting rms acceleration

slightly below the standard curve and well below the minimum complaint

level. Treating the data as being comparable to single footfall or an impulsive shock which is the case for the plot of the peak rms acceleration

it is obvious that the results are not realistic.

An interview with the occuDants of the rhv«=l 1 i na where the test data were

obtained indicated that the vibrations produced by walking were not considered

annoying although at times they were perceptible. It is interesting to

note that a deaf couple living in this housing project used floor vibrations

(which are, in this case, quite low) as a means of sensing the presence of someone in the room.

13
7 . SUMMARY

In this report a new approach to the treatment of perceptible floor

vibrations is presented based on the concept that human activity and human

response to this activity are random variables. The reasons for choosing

the random process approach to this floor vibration problem are given along

with an overview of the method in functional form. A description of the

random process and the statistical functions used to treat random data is

presented

.

Spectral analysis is chosen over time, history, as the method for

analyzing floor vibration data. A computationally efficient procedure for

obtaining the Fourier spectrum is given along with the associated computer

software

.

The current state of knowledge on human reaction to externally applied

vertical transient vibration is discussed. One record of human activity

obtained in a field test is analyzed using the techniques presented in this

report and the results compared with a proposed standard for vibration and

shock limits for occupants of buildings.

14
8 . REFERENCES

1. Lenzen, K.H. , "Vibration of Steel Joist-Concrete Slabs Floor System - Final Report," Studies in Engineering Mechanics , No. 16, Center for

Research in Engineering Science, the Univ. of Kansas, Lawrence, Kansas,

August, 1962.

2. Onysko, D.M. , "Performance of Wood-Joist Floor Systems - A Literature

Review," Forest Products Laboratory, Information Report OP-X-24,

Ottawa, Ontario, Canada, January, 197 0.

3. Glamabos, T.V., Gould, P.L., Ravindra, M.R., Suryoutomo, H., Crist,

R.A., "Structural Deflections - A Literature and State-of-the-Art

Survey," Building Science Series 47, National Bureau of Standards,

Washington, D.C., October, 1973.

4. Harper, F.C., Warlow, W.J., Clarke, B.L., "The Forces Applied to the

Floor by the Foot in Walking," National Building Studies, Research

Paper 32, Building Research Station, Watford, England, January, 1961.

5. Geortzel, G., "An Algorithm for the Evaluation of Finite Trigonmetric

Series," American Mathematical Monthly , Vol. 65, No. 1, January, 1958,

34-35.

6. Rubin, S., "Concepts in Shock Data Analysis," Shock and Vibration Handbook

(Harris, CM., Crede, C.E., Eds.) Vol 2, Ch. 23, McGraw-Hill Book Co.,

New York, 1961.

7. International Organization for Standardization, "Guide for the Evaluation

of Human Exposure to Whole-Body Vibration," Draft International Standard

ISO/DIS 2631, May, 1972.

8. Splittgerber , H., "Vibration and Shock Limits for Occupants of Buildings,"

ISO/TC 108/SC4/WG2, Draft document, June, 1973.

9. Lenzen, K.H., Dorsett, L.P., "Effect of the Variation of Structural

Parameters on the Vibrational Characteristics of Steel Joist-Concrete

Slab Floor Systems and Suggested Designs," Studies in Eng. Mech. No.

32 , Univ. of Kansas, Lawrence, Kansas, Aug., 1968.

10. Reihler, H., Meister, F.J., "The Sensitivity of Humans to Vibrations,"

Forsc auf den Gebiere des Ingenieureuesens , Vol. 2, No. 11, November,

1931.

15
9. APPENDIX A Field Test Procedure

The physical test procedure to be used must be designed for field

application. Because of this requirement, the instrumentation must be

simple and require a minimum of set-up time with little inconvenience to

the occupants where the test is being made. A transducer which has the

required sensitivity and response time and measures acceleration is preferred

since it directly measures the quantity generally considered by most to be

related to human response. Velocity transducers may also satisfy the above

requirements but the subsequent differentiation needed to produce acceleration

is often difficult to achieve without loss of information. Displacement

transducers are difficult to satisfactorily employ in the field for this

type of test, due to physical restraints of installation and reference

datum. In addition, these transducers generally do not have the proper

time response. Transformation of the displacement data to acceleration

would also result in a considerable amount of lost information.

A displacement transducer was used in addition to an accelerometer in

the field test.

Description of Test Site

The site selected for application of this methodology was a single

family attached house (townhouse) with a wood-joist floor system located in

Gaithersburg , Maryland. The house was occupied at the time of the test and

furnished as shown in Figure 5. The floor system was typical of that found

in many houses and consists of 2 x 8 wood joists spaced 16 in on center

overlaid with 1/2-in plywood. The joists have a clear span of 12 feet 10

inches with diagonal bridging at mid-span. A synthetic wood floor tile is

attached to the plywood subfloor by an adhesive. The floor construction is

shown in Figures 6 and 7.

Floor response measurements were obtained from both an accelerometer

and linear variable differential transformer (LVDT) . The range of the

accelerometer is from 10 mg to 0.001 mg. Figure 6 shows the location of

this instrumentation. The LVDT was positively attached to the joist at

mid-point of the clear span and the accelerometer set on the floor directly

above the LVDT. Output from the LVDT and accelerometer were recorded as analog signals. While the frequency response of the electronic equipment

was adequate to reproduce the time histories, the signal conditioning

equipment for the accelerometer has a cut-off frequency of 25 Hz. This is

generally considered satisfactory for most wood-joist floor systems.

16
Test Procedure

A series of tests were made using two procedures: (1) human activity

represented by a 170-lb man walking in a normal traffic pattern, and (2)

dropping a 25-lb bag from several heights. The first procedure is typical

of a minimum human activity on a floor system while the second one has been

used in the past as a method to assess the adequacy of floor systems from a vibrational viewpoint.

During the activity test, the walker utilized a path which was considered

typical of the traffic patterns in that room. This path was followed for

approximately 4 0 seconds.

The bag was dropped from heights of 7 , 6, 5 and 4 inches. Several

tests were made for drop heights greater than 7 in and, in each instance,

the accelerometer raised from the floor in a rebound action.

Test Results

The walking test is the most significant data obtained since it is

representative of the human activity that most frequencly occurs on floor

systems. In addition, it is a minimum type of activity as compared to

running, jumping, etc. The accelerometer is a representative measure of

the human response as its output is that which is preceived by a human

standing or sitting in an unpadded chair at the location of the transducer.

A typical segment of this record is given in Figure 8

.

17
10. APPENDIX B Computer Program

The computer program developed to perform the spectral decomposition

of a digitized time history is presented. The main program acquires a channel of digitized data, called data counts, from tape and then converts

this to data in engineering units. The program assumes that the data is

going to be either a record of displacement measurements or acceleration and

makes the conversion upon being given the appropriate conversion factor and

displacement or acceleration identification. Following the conversion, the

data is made available to a subroutine which, using the Fourier techniques

described in section 4.0, performs the spectral decomposition. The program

is written in Fortran V and is operational on a UNIVAC 1108 system equipped

with Exec 8. Only minor modification needs to be made to the main program

for operation on any computer with a Fortran V compiler assuming the input

data tape is compatible with the one for which the program was written.

However, the subroutine can be used, without any modification, on a system

with a Fortran V compiler.

18
COMPUTER PROGRAM FOR FLOOR VIBRATION ANALYSIS

OIMrMs'ON JwEAOfl?), !tI m E(«)i 1 m ijXtU?), IGA|N(|2I, CNVFc(l2),

I r«UF{ 1 ?OoO) i Ml2"nO>, Tt*lF( 1 2nobj , 10COI I?) • IOP112)

Dl MpMS ' ON »EL VEC I AO"" j , F M r. v E C < 600" ) , FRFQ<6000)

EQUT VA' ENCF (IHF»n ( |),I'N

i 1 <THEani?>,IPM»g),(!HEA0(3l,NSCPPC),

1 t IhEAD ( «» ) ,-NdCpfL I % « I WFAO f «VJ. ,M5*-T.N > , ( I HE AO ( 6 ) .NSSEC i ,

1 f ThEAO { 7 ) ,"«MeFt J , f IMF AO f n I , T T I me ( | ) ) , ( IHFAO ( 9 ) , IT I ME (21 » ) ,

| ( I HE AD ( I 0 )• » t T t M E ( 3 ) ) i (IMF A D ( 11 ) » I T I M E. ( H ) ) . ( I HE Ad ( 1 2 ) , NG.HSN ) C»»»»»REAn poor, PA" Ann r u

A

hN F L Ca^ps • REAP<5,60n) NOCH, tPsOFC

600 FORMA T ' T 5 ,E 1 2 . t , ! S ) DO S I e | , NDCH

5 RFAOlS^OO) 10rP(I|, CMVFC ( 1 ) , T OP ( 1 ) N O C M R 0 x 0 C«»«»«RE A n H«"j DING "FC^P", 4 Nn f H 4

*' ME L 'OS, WITH G A I N •••••••••••••••••••

CALL N T BAN(7,2,!?,tHFAP,l ,2?l

IF ( l*i .

L

T . n i ro_To on CALL N

T RAN ( 7 , 2 , 'iCVK'l* 1 'UIK T ,i. • 2? 1 IF r L*l . L T t 0 ) r,O.Tp oo CALL M

T RA M ( 7 , ? , nrn^M ,

I G * ' N .1 » 2? > IF ( L*i . L T , n | /;0 TO RO

c

• • • • n R I t e "EAnlMf, amo f h»nmfi nos. »ith r a i n» «••••••••••••••••••••••• • *R!TE(*,60I ) f'*. Ip>'»'. MSCPPC, MPCPFI., N5MIN, NSSEC, NSMSEC.

1 I ITITt 1 1 1 1 1 . t ) . hChsh

601 FORMAT' IHI., It*, 'SITE ID *. ' S / I 5 X , 'RUN N0,«, IS / 15X,

I 'S^AMS / PFCOOn 15/ 15t, tRECPPDS / FILE IS /

1 ISXj 'SCAN INTEoVaL mn = « , I M , ? X , i Sf m •

, 1 « , 3X , • Ms- • , I < /

i tsx; »time , i « . • : ' . ' 2 , • : • ,

t 7 , • : •

, i 2 / isx,

1 iN". channels sca n '

*• // > *R I TF I «• ,602 ) ' 1 "IH' ( J > t ' r'A J M ( J ) , J-1.NCH5N 1 60? FORMAT t | 5 x, •CHAHnFl. •

, tj, 2X, • G A I N •

, T3 ) ^•••••RFAP A"0 m R I T F CHAnNrl. OF o»TA C 0 U N T «:•••••••••«••»•••••••«»••••• • • 500 CALL Mtran(7,2, |?OnoitBUF tl_,22)

\r ( L*l .LT. 0 I r,0 Tp 00 NDCHPD « NOChPO 1 NPTPTH MSCPRC • flRr p FL WR I TF ( * , 603 ) ! MUX T , NnCHR" 1 , I iniiFlM, I-l.NDTPCH ) 603 FORMAT' I H 1 • t»i VH^•IMFL ,

. 13, ISI7 / ( ISX, 1SI7 ) ) C»«»»»iF oaT« Is OiSPLACrMfNT ( c">OE = 1 ) PEMPVE OFF SET •••••••••••••••

IF (I "CO ( N^CHPO ) .Ep. 0 ) GO TO | 00

DO 10 ' 1 .NHTPCH

10 murin « 1 p of ( 1 » _ t«uF(i)

C ••• • •COMVFR'' I On FROM OATA COUMT TO ENfilNFFPINC UM I TS»»» ••••••••••••••••

100 DP 20 t I ,WDTPCH

20 XII) » ( J 0 • IBUFrt) > / I ?OMfl • 1 GA I N 1 NDCHRO ) • C N V F C ( NO C H R D ) ) C»»»»«FILL T'ME a»RAV •••,•«•••••••••••••••••••»•••••••••••••••••• •••••••

OFL T A T a ( NS"t N • 60 MSSFC MS^SFf • 0.0P1 ) • TPSDFP

T 1 mf ( M n, DO >o '^.'iDTPCH

30 T1MFU» TtMEJI-ii 0Fi.T a t C»»»»««RITF <"nKvFPTFO n» T A IF 1 |*rD<rjhcMPr)| '.F«. 0 ) GO to 20"

WR I TF I 4 , 6n«t 1 I muxt ( finCHPT 1 A0<» FORMAT! IHI, <», •rHA M *jFl ,

1 13, • IS P I S P L A C E ME N T IN INC W ES« > 60 TO >nl) 200 «P 1TF (

<• ,«riS » rMIHTjNriCHPPI

AOS F0R|| A T« ln|, •»t, VH,NHFL». 13. • 1 S ACCELERATION IN MG • ) 300 WR I TF ( * , 606 ) I T

I

M E (III XII). I«|,HDTPCH)

f,nt, fO B H 4 T' 1n( F>>.3, r 1 ? , >t 1

) c »«»».SPFC Tf'

,

l. ANALYSIS I rOiE I !«•*•••••• ••••••••••••••••••

IF ( l"p(MnfHRn) ,r 0 . p ) G n To ng-i

CALL F"li9TR » X t

*" n TPcH" n t LTAt , EMp«'FC , PFLVCF . FREO ) MOp IF t M» c n u n - "r>CH ) SCO

CALL "'RAH ( 7 ,

1

n )

r; o T 0 ""^

o 0 V R I T F 1 * , A 0 7 i 607 FORH*T'« ! =•*. 13, ) 205 JTOK

E

l .'P

19
SUBROUTINE FOR SPECTRAL DECOMPOSITION

SUBROUTINE ^OU« Tf?{x.»l»0f"lT4T,E M r;vEC,PFLV*"C.FREQ)

DIMENSION RECF0R(3>, RELVEc ( AOnn ) , E M r, v E C ( 6 0 0 0 ) , FREQ<6000>,

t X (

t 2000 ) L N - 1 M = N ' 2 00 I no K = n

,

m 0 M E G A * t> . 2 8 3 t 8 5 3 • < / H

C S 0 M E G = C0S(0mega)

SNOMEG = S I M ( 0

M E G A 1 C»»»«» EVALUATE PECDRRSInN FORMULA ••••••• • ••••••••••

R E C F 0 R I

\ ) = X ( N ) R E C F 0 R ' 2 1 = x( ) R E C F n R < 3 ) = x( N - 2 ) + 2. • C S 0 M E G • P E C F o R ( 2 > - RECFOR(I)

DO 2 n-0 t = 3,L

RECFOR Ml ? PFcFnp,2i

RECFOR I 2 ) = PfCFnR , 3

,

200 RECFOR ' 3 ) = x( M - 1 ) 2. . CsOMEG • RECTOR ( 2 J - RECFORH)

C • • • • • CO»"Pll T

E REAL A NO IMAGINARY VECTORS ••••••••••••••••••••••••••

F R E Q < K ' = K / ( N « nEl_TAT ) EMGVEC K > = R E f F o p ( 3 ) . «; N 0 M E G • DFlTat

100 R E L V E r ' k > = ( X ( n J + R c C F 0 w ( 3 ) « C S

0

M E G - RECF0R(2) ) • 0 E L T A T £••««« COMPUTE ABSOLUTE vAt.UE •••• ••«• ••«•••»•••*•

DO 3 n P <=0. M 300 X ( K ) = 5QPT( RELVEr<<'»«2 E M G V E C t K ) • * 2 ) WRlTF(A,2n!) ( ERErj(<>. RElV'EC ( < ) , EMGVEC ( K ) , X ( K ) , K=0,M ) 20 1 FORMAT* | H|, MX, nF\b,9 / ( , t*F\ 6.8 / ) ) RETURN

END

i IS( l "MMNHS-OC

20
Figure 1 VERTICAL LOAD FROM FOOTFALL [4]
F(t)

ure 2 ENGINEERING VIBRATION MODEL

22
s-6u»'sn-| naow U3iyn0d

23
f,Hz

Figure 4 HUMAN RESPONSE THRESHOLDS

TRANSIENT VERTICAL VIBRATION

24
25
I CM

I

I

I I

I I

I

I 8 JOISTS

6 in O.C 2V

i

i

LVDT I I ACCELEROMETER

i '

i

1

I

PARTITION

:r:r:r:c:r zX~.lzii '

SUPPORT BEAM

DOWN

F= I M «=

Figure 6 FLOOR PLAN OF TEST FLOOR

26
FIGURE 7 TEST FLOOR UNDERSIDE

27
ri v 1 . H

a.

->c

t-< Q>

-1

T

in 1

in

a.

E-1 OO

(Xuu< T1

i-l

Piwc

U:a<

Cmxw

Gwu

i-i aoH

dl

or, o-o

h4

l-l <uI—I

Cm

>l

E->

f\l oo Wu<

Cma:!—I a

OOI £

en£ Oi—i <

C<WWUU

CPE o!—I

En<w

\-A WUU< E Do
NBS-1 <A (REV. 7-73)

U.S. DEPT. 0 F COMM. BIBLIOGRAPHIC DATA

SHEET

1. PUB1.K A 'IT ON OK RKPOR 1 NO. 2. Gov't Accession

No.

3. Recipient's Accession No.

4. rn U-: and subtitli-:

Correlation of Floor Vibration to Human Response

5. Publication Date

6. Performing Organization ( ode 461 . 01

7. AUTHOR(S)

James R. Shaver

8. Performing Organ. Report No.

9. PERFORMING ORGANIZATION NAME AND ADDRESS

NATIONAL BUREAU OF STANDARDS

DEPARTMENT OF COMMERCE

WASHINGTON, D.C. 20234

10. Project/Task /'Work Unit No. 4618381

11. Gontract /Grant No.

12. Sponsoring Organization Name and C omplete Address (Street, City, State, ZIP)

Office of Policy Development and Research

Department of Housing and Urban Development

Washington, DC 20410

13. Type of Report & Period Covered

Final

14. Sponsoring Agency Code

Samp as 1?

15. SUPPLEMENTARY NOTES

16. ABSTRACT (A 200-word or less factual summary of most significant information. If document includes a significant

bibliography or literature survey, mention it here.) A new approach to the problem of perceptible floor vibrations is presented predicated

on the realization that human activity and human response to this activity are random variables. Techniques for data reduction are discussed and a detailed

description of one approach is given along with the associated computer program.

Data from floor vibrations is compared with current criteria for human response

to vibration.

17. KE.Y WORDS (six to twelve entries; alphabetical order; capitalize only the first letter of the first key word unless a proper name; separated by semicolons

)

Analysis; experimental; floor systems; human response;

random process; spectral analysis; vibration.

18. AVAILABILITY [XXUnl imited

[ lor Official Distribution. Do Not Release to N I IS

Order Krom Sup. of Doc, U.S. Government Printing Office Washington, D.C. 20-102, SD Cat. No. C H

[^J Order I'rom National Technical Information Service (NTIS)

Springfield, Virginia 22151

19. SECURITY CLASS

(THIS REPORT)

UNGL ASSIEIED

20. SEC URITY C LASS

(THIS PAGE)

UNCLASSIFIED

21. NO. OF PAGES

32

22. Price

$4 .00

USCOMM-DC 29042-P74