Abstract

Based on classic vibrational bending theory on beams, this paper provides comprehensive analytical formulae for dynamic characteristics of two equal span continuous timber flooring systems, including frequency equations, modal frequencies, and modal shapes. Four practical boundary conditions are considered for end supports, including free, sliding, pinned, and fixed boundaries, and a total of sixteen combinations of flooring systems are created. The deductions of analytical formulae are also expanded to two unequal span continuous flooring systems with pinned end supports, and empirical equations for obtaining the fundamental frequency are proposed. The acquired analytical equations for vibrational characteristics can be applied for practical design of two-span continuous flooring systems. Two practical design examples are provided as well.

1. Introduction

Vibrational serviceability performance of timber floors has become an important issue in the world due to their resonance frequencies and low material masses. In Europe, Eurocode 5 [1] has been widely used for design of timber floors to satisfy serviceability limit state criteria [2], in particular vibrations, because they are often governing the design of timber floors. The fundamental frequency, unit point load deflection, and unit impulse velocity response as the key vibrational parameters need to be checked. The methods for obtaining these parameters and their limits are presented in EN 1995 Part 1-1 [2] and the National Annexes of individual European countries [3].

Human beings are considered as precarious sensors of vibrations, and their distress to timber floor vibrations concern many researchers, and human activities and machine-induced vibrations can cause distress. Human sensitivity and perception are basically related to structural vibrations.

2. Previous Research Studies on Timber Floor Vibration

Over past decades, extensive investigations have been conducted on evaluating the dynamic performance of timber floors and human vibrational perception in many European countries, Canada, Australia, and Japan. Ohlsson [4] investigated human-induced vibrations by testing a number of timber floors, and his work has been implemented in Eurocode 5 for timber floor design [1]. Chui [5] proposed use of the root mean square (r.m.s.) acceleration for design based on his field tests on timber floors and. Hu [6] dynamically tested I-joist floors and numerically simulated the vibrational behaviour of these floors. Eriksson [7] measured the dynamic forces of low frequencies caused by walking, running, and jumping. Smith [8] explored the vibrational issues for timber floors. Bahadori-Jahromi et al. [9, 10] innovated timber floors with multiwebbed joists and conducted static and dynamic tests on these floors. Homb [11] assessed human perception based on the results of impact-induced low frequency vibrations on timber floors. Ljunggren [12] assessed the dependence of human perception on the dynamic properties of lightweight steel framed floors. Toratti and Talja [13] established body perception scales to reflect the significant dependence of disturbance on different vibrational sources. Zhang et al. [14, 15] systematically assessed the effects of influencing parameters on the vibrational performance of timber floors constructed with various joists. Weckendorf [16] measured the vibrational characteristics of I-joist timber floors. Labonnote [17] categorised damping in timber components. Zhang et al. [18, 19] assessed the vibrational responses of timber floors constructed with metal web joists by altering structural configurations. Jarnerö et al. [20] measured the dynamic characteristics of a prefabricated timber floor with various boundary conditions at different construction stages. Two textbooks were also written on structural timber design to Eurocode 5 [21, 22].

Recently, considerable attention has been paid to the dynamic responses of timber floors constructed with LVL (laminated veneer lumber) beams, glulam beams and CLT (cross-laminated timber) panels, and TCC (timber-concrete composite) floors. Basaglia et al. [23] assessed the vibrational measures on timber floors and compared the design criteria between Australia and Japan by testing a 9 m span LVL ribbed-deck cassette floor. Ebadi et al. [24] explored the vibrational response of glulam floors. Ui Chulain and Harte [25] investigated the influence of modern timber fixing systems and extra mass on the serviceability behaviour of one-way or two-way CLT floors in laboratory. Koyama et al. [26] measured the vibrational characteristics of a CLT floor subjected to walking vibrations. Bui et al. [27] experimentally investigated the natural frequencies, damping ratios, and mode shapes of multilayered timber beams. Lanata et al. [28] attempted to establish correlations between the dynamic response of in situ timber and TCC floors and human comfort perception by varying building typology.

2.1. Design Codes in European Countries

EN 1995-1-1 [1] has adopted three parameters proposed by Ohlsson [4] for controlling the vibrational serviceability design of timber flooring systems. For residential floors which are simply supported along all four edges, it requires that the fundamental frequency or the lowest first-order modal frequency f1 in Hz should satisfy the following equation:where is the floor span in m, is the equivalent plate bending stiffness of the floor about an axis perpendicular to the beam direction in Nm2/m, and is the mass per unit area in kg/m2. Equation (1) is a simplified design equation which is actually applied for two-side supported floors, and the effect of the transverse stiffness is omitted because of small errors caused. EN 1995-1-1 does not clearly indicate how the participating mass should be calculated and whether the composite effect of floor joists and deck in the floor direction should be considered. The design equations and the corresponding limits for f1 proposed from EN 1995-1-1 and various National Annexes in European countries have been previously summarised [3].

The majority of European countries have directly adopted equation (1) and the limit of 8 Hz specified in EN 1995-1-1 [29]. However, Austria [30] only adopts equation (1) for two-side supported floors and provides a fairly accurate equation for four-side supported floors by including a quartic term about the floor span to width ratio, , to reflect the effect of the transverse stiffness as

Similarly, Finland [31] provides a more accurate equation for four-side supported floors by including both second- and fourth-order terms about for the effect of the transverse stiffness with a raised frequency limit of 9 Hz as

Both Austria and Finland specify that the floor mass m should be determined by using a quasipermanent combination of dead and imposed loads, as specified in EN 1990 [2]:where is the mass due to the characteristic dead load and is the mass due to the characteristic imposed load . is the factor for quasipermanent value of a variable action, and its values for different building categories can be taken from the relevant tables of EN 1990 [2] or the National Annexes to EN 1990. In general, can be taken as 0.3 for residential and office buildings, 0.6 for congregation areas, and 0.8 for storage areas.

Spain [32] limits the values for the fundamental frequency f1 for all construction materials including timber. More recently, Hu et al. [33] attempted to develop a baseline vibration design method for timber floors by combining all available databases in the world so as to eventually develop an ISO standard. The future work will be to harmonise the calculation methods and to rederive the baseline design criterion by using the collective floor database.

2.2. Significance of Current Research

With the development of engineered timber products, floor joist sizes can be manufactured much larger for longer floor spans or multispans. Some research work on structural dynamic characteristics of general continuous beams has been reported [34, 35]. The information for determining vibrational characteristics of continuous timber flooring systems, however, is still limited. In this paper, the characteristic equations for modal frequencies and shapes of continuous floors of two equal spans with various boundary conditions or floors of two unequal spans with pinned boundary conditions are deduced. Two design examples of two-span timber flooring systems are also presented.

3. Equations of Motion for Two-Span Continuous Floors

For two-span continuous timber flooring systems with equal joist spacing but various end supports, they can be treated as two-span beams and analysed on a single two-span continuous timber floor joist (see Figure 1 in which 0 ≤ x1 ≤ L1 and 0 ≤ x2 ≤ L2).

3.1. Equations of Motion

For a two-span continuous Bernoulli–Euler beam with a uniform cross-sectional area, mass density, and flexural stiffness, the equation of motion for each beam span, i.e., transverse displacement yi (xi, t) versus time t for span i, is given aswhere is the flexural stiffness, is the viscous damping coefficient, is the mass density per unit volume, is the cross-sectional area, is the time varying external load, and is the span number (i = 1, 2) [36]. Assume can be separated for and as

Substituting equation (6) in equation (5) and ignoring the harmonic part yields the homogeneous solutions as

The natural frequencies and coefficients Ci1 to Ci4 of the beam can be obtained by applying appreciate boundary conditions. The natural frequencies fn are determined aswhere is the circular modal frequency for nth mode in rad/sec, is the equivalent plate bending stiffness of the timber floor about an axis perpendicular to the beam direction in Nm2/m and is given as , is the flexural stiffness of the floor joist in Nm2, s is the floor joist spacing in m, and m is the mass per unit area in kg/m2.

3.2. Boundary Conditions

Four typical boundary conditions for two-span continuous beams are considered with various end supports, e.g., free, sliding, pinned, and fixed boundaries (Table 1), which forms sixteen combinations of flooring systems.

4. Vibrations for Continuous Floors with Two Equal Spans

For continuous timber flooring systems with two equal spans, L1 = L2 = L.

4.1. Frequency Equations

To establish the frequency equations of two equal span flooring systems and to determine the modal frequencies and mode shapes, the displacement equation (7) should be used in conjunction with the boundary conditions listed in Table 1. To demonstrate the procedure, a typical two equal span continuous beam with left end fixed and right end simply supported is used here (Figure 2).

For spans 1 and 2, equation (7) can be rewritten as

Differentiating equations (9) and (10) for three times yields

4.1.1. Boundary Conditions for Span 1

(a) at x1 = 0, so equation (9) degenerates to(b) at x1 = 0, so equation (11) degenerates toThus, equation (9) becomes(c) at x1 = L, so equation (19) degenerates to

4.1.2. Boundary Conditions for Span 2

(d) at x2 = 0, so equation (10) degenerates to(e) at x2 = 0, so equation (15) degenerates toSolving equations (22) and (23) simultaneously yieldsThus, equation (10) becomes(f) at x2 = L, so equation (25) degenerates to

4.1.3. Boundary Conditions at the Interior Support

Combining equations (11), (17), (18), and (21) and letting x1 = L yieldsCombining equations (14), (24), and (27) and letting x2 = L yieldsCombining equations (12), (17), (18), and (21) and letting x1 = L yieldsCombining equations (15), (24), and (27) and letting x2 = L yields(g) at x1 = x2 = L, so combining equations (28) and (29) yields(h) at x1 = x2 = L, so combining equations (30) and (31) yields

A nontrivial solution for equations (33) and (35) exists only if the determinant of the coefficient matrix vanishes, which yields

Thus, the frequency equation can be obtained as

Similarly, based on equation (7) and the boundary conditions for different two equal span flooring systems which are listed in Table 1, the frequency equations for other fifteen cases with two equal span floors and various end support conditions can be deduced, as illustrated in Table 2.

4.2. Modal Frequencies and Shapes

From the frequency equations, the modal frequencies for two equal span timber flooring systems can be obtained numerically. Here commercial software MathCAD is used for such purpose, and Table 3 illustrates the values of βiL for the first four modal frequencies of sixteen timber flooring systems with two equal spans and various boundary conditions. The symbol in Table 3 indicates the rigid mode for free-free boundary conditions.

The vibrational mode shapes only for the first modes of two equal span continuous floor beams with various boundary conditions are illustrated in Figures 36. Figure 3 illustrates the first vibrational modes for two equal span continuous timber floor beams with free left ends and various boundary conditions for right ends. All the mode shapes are not in scale. The left ends sustain free drops, and the right ends have various shapes which depend on predefined boundary conditions. At the interior support, the slopes of the mode shapes for both left and right spans are continuous. In general, the fundamental frequency parameter β1L increases in the following order for the right end boundary conditions: free, sliding, pinned, and fixed. It should be mentioned that the first vibrational mode for the free-free boundary conditions should be rigid with a zero modal frequency.

Figure 4 illustrates the first vibrational modes for two equal span continuous timber floor beams with sliding left ends and various boundary conditions for right ends. The left ends have zero slopes, and the right spans have various shapes which depend on predefined boundary conditions. The fundamental frequency parameter β1L increases in the same order for the right end boundary conditions: free, sliding, pinned, and fixed.

Figure 5 illustrates the first vibrational modes for two equal span continuous timber floor beams with pinned left ends and various boundary conditions for right ends. The left ends have zero displacements, and the right spans have various shapes which depend on predefined boundary conditions. The fundamental frequency parameter β1L increases in the same order for the right end boundary conditions: free, sliding, pinned, and fixed.

Finally, Figure 6 illustrates the first vibrational modes for two equal span continuous timber floor beams with fixed left ends and various boundary conditions for right ends. The left ends have zero displacement and slope, and the right spans have various shapes which depend on predefined boundary conditions. The fundamental frequency parameter β1L increases in the same order for the right end boundary conditions as mentioned above.

The shapes for higher modes are more complex and will not be discussed here further.

5. Vibrations for Continuous Floors with Two Unequal Spans

5.1. Frequency Equation

For continuous floors with two unequal spans (L2 = αL1) and pinned end supports (Figure 7), the frequency equation can be obtained based on equation (7) and the boundary conditions in Table 1 for pinned ends.

5.1.1. Boundary Conditions for Span 1

(a) at x1 = 0, so equation (9) degenerates to(b) at x1 = 0, so equation (12) degenerates toSolving equations (38) and (39) simultaneously yieldsThus, equation (9) becomes(c) at x1 = L1, so equation (41) degenerates to

5.1.2. Boundary Conditions for Span 2

(d) at x2 = 0, so equation (10) degenerates to(e) at x2 = 0, so equation (15) degenerates toSolving equations (44) and (45) simultaneously yieldsThus, equation (10) becomes(f) at x2 = L2, so equation (47) degenerates to

5.1.3. Boundary Conditions at the Interior Support

Combining equations (11), (40), and (43) and letting x1 = L1 yieldsCombining equations (14), (46), and (49) and letting x2 = L2 yieldsCombining equations (12), (40), and (43) and letting x1 = L1 yieldsCombining equations (15), (46), and (49) and letting x2 = L2 yields(g) at x1 = L1 and x2 = L2, so combining equations (50) and (51) yields(h) at x1 = L1 and x2 = L2, so combining equations (52) and (53) yields

A nontrivial solution for equations (55) and (57) exists only if the determinant of the coefficient matrix vanishes, which yields

Thus, the frequency equation can be obtained as

5.2. Model Frequencies and Shapes

From the frequency equations (58), (59), and (60), the modal frequencies of timber flooring systems with two unequal spans and pinned end supports for various span ratios can be obtained numerically. Again, commercial software MathCAD is used for such purpose, and Table 4 illustrates the values of βiL for the first four modal frequencies of timber flooring systems with eleven typical floor span ratios between 0 and 1. Figure 8 shows the fundamental vibrational mode shapes for α = 0.00001, 0.2, 0.4, 0.6, 0.8, and 1.0, respectively.

When is very small, e.g., 0.00001, the mode shape for the left span approaches to the shape for a single span beam with pinned left end and fixed right end, and the fundamental modal frequency is the largest. The corresponding amplitude for the right span is very small. With the increase in α value, the fundamental vibrational frequency f1 decreases and the mode shape for the right span which always follows a half sine wave becomes more extended with higher amplitude. When α = 1.0, a full sine curve is obtained for the whole two equal span beam. This trend continues for α larger than 1.0. For this case, however, the right span can be treated as L1 and the left span treated as L2, so the same calculations can be conducted. At the interior support, the slopes of the mode shapes for both left and right spans are continuous.

5.3. Empirical Equations for Practical Design

Figure 9 illustrates that for α = 0.0 to 1.0, a cubic relationship between the frequency parameter β1L and the span ratio α can be established empirically with the linear regression coefficient R2 = 0.9999 as

This equation can be used to determine the fundamental vibrational frequency, f1, of continuous timber floors with two unequal spans for practical serviceability design as soon as the span ratio α is known.

Sometimes, a linear relationship between β1L1 and can be more convenient for practical engineers to quickly estimate the fundamental frequency of a timber floor. Here, only the values of β1L1 for α = 0 and 1.0 are used to form a simple linear empirical equation as follows without large errors from the accurate values (Figure 9):

6. Design Examples

6.1. Two Equal Span Floor with Solid Timber Joists for Fixed-Pinned Ends

A two-span timber floor is designed for a domestic timber frame building. It is constructed with continuous solid timber joists (Figure 10). The floor has a width B = 5.0 m and a length L1 = L2 = 4.8 m for each span. It is constructed with 75 mm × 220 mm C24 solid timber joists at a spacing s = 600 mm. The P5 particleboard with a thickness of 22 mm is chosen for the decking, and the Gyproc plasterboard with a thickness of 12.5 mm is chosen for the ceiling. The total self-weight of the flooring system including the timber joists is assumed to be 60 kg/m2, and Service Class 2 is assumed. The imposed load is Qk = 1.5 kN/m2 from EN 1991-1-1 [37].

From equation (8), the fundamental frequency f1 can be determined aswhere is the frequency parameter for the first mode and β1L = 3.3932 (quoted from Table 3), is the left span length of the floor and L1 = 4.8 m, is the right span length of the floor and L2 = 4.8 m, b is the width of the solid timber joist and b = 75 mm, is the depth of the solid timber joist and h = 220 mm, is the elastic modulus for C24 timber and E = 11 GPa [38], I is the second moment of area for a C24 timber joist and I = bh3/12 = 0.075 × 0.223/12 = 66.55 × 10−6 m4, (EI)L is the equivalent plate bending stiffness of the timber floor about an axis perpendicular to the joist direction and (EI)L = EI/s = 11 × 109 × 66.55 × 10−6/0.6 = 1.2201 × 106 Nm2/m, s is the floor joist spacing and s = 0.6 m, and m is the floor mass per unit area and m = 60 kg/m2.

6.2. Two Unequal Span Floor with I-Joists for Pinned-Pinned Ends

A two-span timber floor is designed for an office timber frame building. It has a width B = 6.0 m, together with lengths L1 = 7.3 m and L2 = 2.92 m for left and right spans, which gives a span ratio α = L2/L1 = 2.92/7.3 = 0.4. It is constructed with the engineered I-joists (JJI joists) produced by James Jones & Sons Ltd in the UK [39] (Figure 11). The top and bottom flanges are manufactured from C24 solid timber with the width b ranging from 47 mm to 97 mm (A to D) and a constant height of hf = 45 mm. The web is manufactured from 9 mm OSB3 which is embedded into the flanges by 12 mm. The 22 mm P5 particleboard is chosen for the decking, and the Gyproc plasterboard with a thickness of 12.5 mm is chosen for the ceiling. The total self-weight of the flooring system including the engineered I-joists is assumed to be 75 kg/m2, and also Service Class 2 is assumed. The imposed load is taken as Qk = 2.5 kN/m2 [37]. The adopted JJI joists for this study are JJI 400D joists at a spacing s = 300 mm and with an overall joist depth h = 400 mm.

From equation (8), the fundamental frequency f1 can be determined aswhere is the frequency parameter for the first mode and β1L = 3.6070 (quoted from Table 4), L1 is the left span length of the floor and L1 = 7.5 m, L2 is the right span length of the floor and L2 = 3.0 m, α is the span ratio and a = L2/L1 = 3.0/7.5 = 0.4, b is the width of the solid timber flanges of JJI 400D joists and b = 97 mm, hf is the depth of the solid timber flanges of JJI 400D joists and h = 45 mm, tw is the thickness of the OSB3 web of JJI 400D joists and tw = 9 mm, h is the overall depth of JJI 400D joists and h = 400 mm, E is the elastic modulus for C24 timber and E = 11 GPa [38], EP5 is the mean elastic modulus for 9 mm OSB3 web and EP5 = 3000 MPa, I is the equivalent second moment of area with respect to C24 solid timber for a JJI 400D joist and I = 2.8345 × 10−4 m4, (EI)L is the same as above and (EI)L = EI/s = 11 × 109 × 2.8345 × 10−4/0.3 = 10.3932 × 106 Nm2/m, s is the floor joist spacing and s = 0.3 m, and m is the floor mass per unit area and m = 75 kg/m2.

7. Conclusions

Based on classic vibrational bending theory on beams, comprehensive analytical formulae for dynamic characteristics of two-span continuous timber flooring systems have been established, including frequency equations, modal frequencies, and mode shapes. Four practical boundary conditions are considered for end supports, including free, sliding, pinned, and fixed supports, and a total of sixteen combinations of flooring systems are created.

The characteristic equations for modal frequencies of two equal span continuous beams with various boundary conditions have been deduced, and four sets of mode shapes have been illustrated. A rigid mode exists for free-free end boundary conditions. A full example for fixed-pinned boundary conditions has been presented to show the deduction procedure.

The characteristic equations for modal frequencies of two unequal span continuous beams with pinned-pinned end conditions have been deduced, and one set of the corresponding mode shapes for various span ratios has been illustrated. A full example for pinned-pinned boundary conditions with a general span ratio has been presented to show the deduction procedure. Also, two empirical equations, cubic and linear, for determining the modal frequency parameters with respect to varying span ratios have been proposed and can be used directly for practical design.

Finally, two practical design examples for determining fundamental modal frequencies have been presented, one for a two equal span continuous floor constructed with solid timber joists and fixed-pinned ends and the other for a two unequal span continuous floor constructed with JJI 400D I-joists pinned at both ends.

Data Availability

The data used in this study will be deposited in a repository.

Disclosure

Part of the work was presented at World Conference on Timber Engineering (WCTE 2016), 22–25 August 2016 in Vienna, Austria.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The assistance of Dr. Martin Cullen toward this paper is highly appreciated.