Abstract

Wood truss joist floors are increasingly used to replace traditional solid timber joist floors in low-rise timber houses. An understanding of the vibration performance of wood truss joist floors is critical for the design and serviceability of the floors. It is difficult to model wood truss joist floors accurately because of the complicated boundary conditions and numerous sophisticated flexible connections. This paper discusses three simplified modeling methods for the wood truss joist floor system. The modeling results were validated by a series of static deflection tests and vibration modes and frequencies tests of a full-size floor. And predictive analysis of human-induced vibration of the floor was also conducted. The vibration characteristics of the wood truss joist floor were investigated. The examination of the applicability of these modeling methods was provided. The results indicate that the point loading deflection more easily affects the deflection of the adjacent joist. However, the deflection influence on other joists that are three spaces away is minimal. Walking on the wood truss joist floor produces steep vibration acceleration fluctuations at the floor center for a relatively long time period. The sheathing-to-joist connections and the metal plate connections of the joists have significant influences on the vibration response of the wood truss joist floor. The modeling method, which considers the flexible metal plate connections and flexible sheathing-to-joist connections, performs best for predicting the vibration performance of the floor.

1. Introduction

Lightweight joist wooden floors are widely used in residential and low-rise commercial buildings and include solid lumber joist floors, wood I-joist floors (stiffening ribs), and wood truss joist floors, as shown in Figure 1. The timber floor systems consist of several parallel joist members supporting the wood flooring or sheathing, which is connected to the joists in a semirigid manner. Solid lumber joist floors and wood I-joist floors were popular in the past in traditional light-frame timber construction. However, the use of solid lumber and I joist for floors has the disadvantages of potential low quality and a limitation of the floor span. What is worse is the improper web holes and flange notch made to allow the passage of electrical cables and ventilation systems; it easily leads to crack growth and even catastrophic damages in the joist [1, 2]. Wood truss joist floors with truss joists connected by metal plates and sheathing connected to the joists by fasteners (nails or screws) are increasingly used to replace the traditional solid timber joist floors in low-rise residential houses because of the advantages of a large span, high strength, easy prefabrication, and convenient installation and access for service pipes and electrical cables. In addition, it is easier to ensure the high quality of wood truss joist floors, and more design freedom is provided for floors and roofs in domestic buildings.

An understanding of the human-induced vibration performance of wood truss joist floors is critical for the design and serviceability of the floors. The timber floor subjected to the footstep loads is supposed to generate a transient vibration response from each footstep declining to close to zero before the next footstep arrives [3]. Moreover, timber floors are more susceptible to vibration than conventional steel or concrete floors due to their low mass. Numerous experimental studies have been conducted on assessing the vibration performance of different kinds of timber floors. Different authors proposed various design criteria for ensuring vibration serviceability of wood floors, listed in Table 1. Ohlsson [4] investigated the vibration performance of wood I-joist floors and proposed design criteria by limiting the fundamental frequency to values exceeding 8 Hz, the point load deflection to 1.5 mm, and the impulse velocity response to. Smith and Chui [5] conducted experimental tests on solid lumber joist floors and suggested limiting the frequency-weighted root-mean-square acceleration to less than 0.45 m/s2 and the fundamental frequency exceeding 8 Hz. Johnson [6] and Dolan et al. [7] proposed a simple frequency-based design criterion without the consideration of other factors. The fundamental natural frequency of occupied and unoccupied floors is required to be greater than 14 and 15 Hz, respectively. Patrick [8] pointed out that the only limitation of fundamental frequency of 15 Hz criterion is equivalent to a criterion limiting the deflection of a single joist under dead load to 1.4 mm (0.055 inches). It was similar to the deflection control of a wood floor. Eurocode 5 [9] uses a frequency-based standard of 8 Hz limit with deflection limitation and unit impulse velocity limit. The vibration performances of a broad range of timber floors were studied by experimental tests. Khokhar [10] analyzed the influence of bridging (lateral element) between solid lumber joists and found that this method was economical and effective for reducing vibration. Weckendorf et al. [11, 12] investigated the vibration performances of wood I-joist floors and open metal-web joist timber floors using laboratory tests. Zhang et al. [13, 14] presented an experimental study on the dynamic performance of open metal-web joist floors with strongbacks and evaluated the effects of the joist spacing, strongbacks, and ceiling on the mode shapes, damping performance, and unit point load deflection. Meanwhile, Zhang et al. [15, 16] investigated EC5 and British National Annex on the design of timber floor joists. Bernard [17] studied the effects of gypsum ceilings, battens, furring channels, posttensioning, blocking, and strapping within the lightweight engineered timber floor in control of perceptible vibration. Jarnero et al. [18, 19] also investigated the impacts of various boundary conditions and the addition of floor layers on vibration response of the floor. Ding et al. [20] conducted a series of field tests on two-storey residential lightwood frame structures (including the fundamental frequency test, ball excitation dynamic vibration test, and pedestrian dynamic vibration test) and assessed the vibration comfort.

Owing to the difficulty of the analysis, present studies on the vibration of timber floors have largely relied on laboratory testing and there are fewer studies on numerical vibration analyses of lightweight wooden floors. Hu [21] conducted a numerical simulation of the dynamic behavior of ribbed plates with consideration of the rotary inertia and shear deformation. Weckendorf et al. [22] also performed a numerical study on the modal frequencies, mode shapes, and static point load deformation of an I-joist timber flooring system. Glisovic and Stevanovic [23] developed the numerical models of the solid lumber joist floor system to investigate the parameter effect on the vibration performances (joist spacing, depth, sheathing thickness, nail spacing, bridging, and support condition).

However, in the respect of wood truss joist floors, the static and dynamic performances related to the vibration serviceability have not been clarified to date. Numerical simulation combined with laboratory studies improves our understanding of the complexities of the vibration response of wood truss joist floor systems. The vibration performances include the fundamental frequency, vibration modes, point load deflection, and transient vibration response caused by the occupants’ activities. In this study, three modeling methods of the wood truss joist floor were developed and the numerical modeling results were verified with a series of static test results, vibration frequencies, and modes of a full-size floor. Subsequently, the models were also used to predict the floor vibration response induced by walking. An examination of the applicability of the modeling methods was discussed. This paper is intended to contribute to the simple and feasible numerical modeling methods of the wood truss joist floor in aspect of vibration performance for a future vibration serviceability study.

2. Overview of the Wood Truss Joist Floor

The wood truss joist floor system consisted of 15 truss joists with a spacing of 400 mm, oriented strand board (OSB) panels (15 mm) covering the joists, and laminated veneer lumber (LVL) plates as the rim board for the whole floor, which was connected to the ends of the joists, shown in Figure 2. The floor system (6.114 m long and 5.689 m wide) was attached to walls with 1.85 m height. The walls consisted of sawn lumber and were sheathed with OSB panels. The OSB subfloor was connected to the top chords of the wood trusses using 4 × 50 mm screws of 150 mm oc around the perimeter and 300 mm oc in the field, shown in Figure 3(a). The major direction of the OSB sheathing in-plane is placed perpendicular to the joist. The mean density of the OSB sheathing was measured as 650 kg/m3. The mean major elastic modulus of the OSB sheathing has been determined as 4080 MPa with the minor elastic modulus of 2080 MPa. Self-drilling screws with the length of 70 mm and the diameter of 4 mm were used to fasten the edge joists to the walls with a distance of 250 mm (Figure 3(b)). The LVL rim board with a thickness of 38 mm and a height of 440 mm was connected to each joist with the same self-drilling screws. The mean elastic modulus parallel to the grain of the LVL board has been determined as 13000 MPa and 1000 MPa for other directions. The mean density of the LVL board was measured as 600 kg/m3. Both the ends of the joists were fastened to the walls with the same two self-drilling screws (Figure 3(c)).

Dimensional lumbers of spruce-pine-fir (SPF) (38 mm × 89 mm) used for the flange and web members compose the herringbone truss joist with an overall joist depth of 440 mm, illustrated in Figure 4. The mean elastic modulus parallel to the grain EL for dimensional lumbers chords had been determined as 8700 N/mm2, with the elastic modulus for the radial direction ER of 600 N/mm2 and for the chord direction ET of 400 N/mm2. The mean density of dimensional lumbers, ρ, was determined as 560 kg/m3. The flange and web members were connected with metal plates produced by MiTek Ltd. The metal plates were made using thin steel plates of 1 mm thickness with zinc coat and teeth of 8.4 mm height. The mean yield strength of the metal plates has been measured as 271.95 MPa, and the tension strength has been measured as 341.32 MPa. The mean elastic modulus of the steel plate is assessed as 203 GPa.

3. Numerical Simulation Methods

Park et al. [24] presented a finite element model for lightweight steel floors (C or Z light-gauge steel joists were welded to the beams); a rigid rink was used to deal with the differences in the centroid of the beam, joist, and flooring material. Zhang et al. [25] proposed a damped plate-oscillator model to obtain the dynamic properties of coupled floor-occupant systems, and the model was verified using laboratory tests of the full-size lightweight cold-formed steel floor. Setareh [26] conducted a comparison of modal testing and computer analysis and sensitivity studies on the dynamic responses of a long steel truss cantilevered floor. However, flexible connections in wood truss joist floors differ significantly from the rigid welded connections in lightweight steel floors or steel truss floors. Barron and Kim [27] added nonlinear springs to characterize the rotation of the heel joint of the metal plate-connected wood roof trusses in the Sap2000 software. Xu [28] determined that each metal plate connection exhibited axial deformation, shear deformation, and rotation deformation and suggested the use of three spring elements to quantify the axial, shear, and rotation behaviors. However, for a wood truss joist floor system, the trade-off between model complexity and computational effort should be considered. Detailed models that consider the behavior of the connection between each web member and the chord member with different angles in three directions are computationally complex, especially for multiangle and multimember intersections in a wood truss floor system.

3.1. Model Description

In this study, the finite element software Ansys was used for the numerical analysis of the wood truss joist floor. Three simplified methods were proposed, in which the connections were simplified as follows: (i) pinned connections for metal plate connections and rigid connections for sheathing-to-joist connections (method I); (ii) orthogonal springs to characterize the in-plane shear and tension or compression performances of the metal plate connections and rigid connection for sheathing-to-joist connections (method II); and (iii) orthogonal springs to characterize the in-plane shear and tension or compression performances of the metal plate connections and springs to quantify the in-plane shear deformation for orthogonal direction between the fasteners and the sheathing (method III). The simplification of the fixed support between the edge joists and the walls and the simply supported connection between the other joists and the walls were assumed in the numerical analysis.

It should be mentioned that the connections of the wood floor commonly exhibit recoverable linear-elastic deformation during floor vibration. The orthogonal springs were used to quantify the in-plane elastic shear deformation for orthogonal direction between the fasteners (4 × 50 mm screws) and the OSB sheathing, illustrated in Figure 5. Based on different definitions of the yield point from load-displacement curves used in different countries [2934] and numerous actual curves for timber connections, the initial periods defined between 0% and 40% of the maximum load of the load-displacement curves of timber connections are conservatively accepted as linear-elastic periods. The curve from 40% of the maximum load to peak load obviously experiences nonlinear plastic deformation. The average envelope curves for the load-slip test of a single-shear screwed 15 mm OSB panel joint were quantified by multilinear curves, shown in Figure 6. The multilinear curves are determined using the point at 40% of the maximum load and the points of obvious slope variation and the peak load of the curves. The numerical vibration analysis related to floor serviceability indicated that the deformations of the timber connections are in the elastic phase between 0% and 40% of the maximum load of the load-displacement curves. Thus, the linear-elastic stiffness parameters of screwed connections (Disp 1 and Load 1) are defined as the displacement and load corresponding to 40% of the maximum load of the curves, listed in Table 2.

For metal plate connections in the model, the connection stiffness of each web member with various angles to chords cannot be quantified, respectively. In order to simply consider the in-plane shear and tension or compression performances of metal plate connections, it was assumed that the multiweb members produce a composite force at a reference node (composite force node). No matter what kind of joint, the elastic stiffness of the connection between composite force node for multiweb members and the continuous chord is assumed to the same. The composite force node for multiweb members is connected to the continuous chords with orthogonal zero-length nonlinear springs, which enables us to quantify the relationship between the component force and the deformation in the directions of the major axis (shear) and minor axis (tension) of the connections while neglecting the effect of the rotation stiffness, taking three web members connecting to the chord as an example, shown in Figure 7.

For the determination of shear and tension stiffness of metal plate connections in the model, the experimental results of single web member connected to the chord with different angles were only taken as the reference. However, the curves of specimens demonstrated discreteness with stiffness and strength [3538]. The interaction of multiweb members connecting to the chord differs the situation of the single web member connecting to the chord. We use the experimental results of the third-point bending loading test for single wood truss to calibrate the elastic stiffness of metal plate connections in the model. As long as the load deformation of a wood truss joist was obtained from an experiment, the calibration of the shear and tension stiffness of the metal plate connections could be easily performed by modeling the static response of the single wood truss test. This ensures the accuracy of numerical modeling in the aspect of the component level and reduces the dispersion of numerical results. When vertical compression behavior occurs, the installment gap of the connection is eliminated under the assumption of infinite compression stiffness of the connection; at this time, the vertical compression was only affected by the stiffness of the members.

3.2. Model Verification Based on the Wood Truss Joist Test

When designing wood floors, a commonly used limit of the maximum deflection produced by a uniformly distributed live load of 1.9 kPa is less than 1/360 of the span for floors with sawn lumber joists [39]. Thus, the midspan deflection of a single wood truss under a total vertical load of 3.34 kN (equivalent to a uniformly distributed load of 1.9 kPa) was evaluated in a third-point bending loading test. Based on the standards JGJ/T265-2012 [40] and ASTM E73-13 [41], the timber truss joist was loaded synchronously at one-third points of the upper chord using the special H steel beam and hydraulic synchronous loading system, illustrated in Figure 8. Both the ends of the wood truss were supported by steel blocks to create simple supports, and an out-of-plane supporting system was used for both the sides near the three equal points of the truss to prevent out-of-plane overturning, as shown in Figure 9. There was a small gap with no contact between the out-of-plane support and the truss to prevent interference with the vertical deflection. The load increase is controlled according to each loading level of 0.25 kN with the maximum load of up to 3 kN and the loading time interval of 1 minute. For floor serviceability issue, the load of 3 kN for the wood truss joist is enough. The displacement sensors (CDP-50) are placed to measure the displacement of both ends of the joist (P1 and P2) and the displacement of midspan of the joist (P3). The average experimental deformations at the midspan of three trusses are shown in Figure 10. For the test data, the deformation of the midspan at the bottom of the truss exhibited a nearly linear increase up to less than 6 mm with an increase in the load from 0 to 3 kN. The maximum deflection of the designed wood truss joist produced by a uniformly distributed live load of 1.9 kPa is about 1/1105 of the span, which is far less than the design limit (1/360), to keep in the elastic status.

In numerical modeling, the material parameters of the wood-engineered products were tested and are provided in Table 3. The shear and tension stiffness of metal plate connections is calibrated to fit the experimental load deformation of single wood truss joist. Table 4 presents the calibrated linear-elastic parameters for metal plate connections. The numerical deformation for modeling the static bending test of single wood truss fits the experimental results well, shown in Figure 10. The numerical increasing trend of midspan deformation is close to the tested results. That is to say, the stiffness of the numerical model of the wood truss joist is close to that of the tested actual wood truss joist. Meanwhile, the modeling peak displacement is also close to the measured data. The modeling result considering flexible metal plate connections (5.56 mm) was slightly larger compared to the tested value (5.43 mm) subjected to 3 kN. The reason for the error may be the connection gap existing of the fabricated wood truss joist floor and discreteness and variability of the natural wood materials including finger joints of flange members and the scarring defect of the dimension wood.

4. Vibration Performance Analysis of the Floor

4.1. Experimental and Numerical Analysis of Static Deflection

The approach of limiting the maximum deflection of floors under a concentrated live load has been widely used because such a limitation relates more directly to floors subjected to the force exerted by the footfall of people. The measurement results were used to assess the performance of the three numerical modeling methods.

4.1.1. Methods for Static Deflection Tests of the Floor

The floor was tested to measure the deformation at the midspan of the joists under a concentrated load of 1 kN. We placed 100 kg iron rollers at different points on the sheathing, and different places at the bottom of the joists were selected to monitor the deflection response and determine the vibration mechanics and integration performance of the floor. The loading points were located on every third joist at the center of the joists and were numbered from W1 at the bottom to W5 at the top (Figure 11(a)). Meanwhile, other loading points were distributed at the same distance in the middle of the joist and were numbered from L1 on the right to L5 on the left (Figure 11(a)). Displacement transducers with a resolution of 0.001 mm (ID-C ABSOLUTE Digimatic543-490B, Mitutoyo, Japan) were installed by attaching to the steel hangers at the monitoring points at the bottom of the trusses. The monitoring points were distributed as follows: at the center of the second joist (P1), the fifth joist (P2), the seventh joist (P3), the middle joist (P4), the ninth joist (P5), the eleventh joist (P6), and the fourteenth joist (P7), as illustrated in Figure 11(b).

4.1.2. Comparison of the Measured and Modeled Static Deformations

Figure 12 shows the deflection of the floor at the monitoring points (P1 to P7) when a concentrated load of 1 kN was applied at W1, W2, W3, W4, and W5 successively in the experiment; the comparison between modeling results and experimental results is illustrated in Figure 13. The red lines indicate the results of method I, the blue lines show the results of method II, and the magenta lines show the results of method III (the same applies to subsequent figures). It was found that almost symmetrical deflection corresponds to the locations of the loading points on the middle joist. Since the elastic stiffness of the metal plate connections was not considered in method I, the deformation was significantly underestimated in this method. In contrast, the modeling results of methods II/III were in good agreement with the experimental results and exhibited a similar trend.

When the load was applied to W1, the maximum experimental deflection (0.9 mm) occurred at P7, which was located under the loading point; the minimal deflection near 0 mm was observed at P6, which was located three joists away from W1 (Figure 12). The modeling results showed that method I provided the lowest value of 0.62 mm in contrast to the experimental result of 0.9 mm at P7, and significant differences in the values were also observed at the other monitoring points (P1 to P6). Method II provided a slightly lower value at P7 (0.88 mm) than the experimental result (0.97 mm), and method III yielded a slightly higher value (1.1 mm). When the load was symmetrically applied to W5, a similar trend of the model results compared to the experimental results was found.

When the load was applied to W2, the maximum experimental deflection (1–1.12 mm) occurred at P6, which was under the loading point W2. The deflection was small (0.1 mm) at P5, which was two joists away from the loading point W2; the influence on the other monitoring points was negligible. Method II provided a slight underestimation of 0.98 mm at P6 and P2. Method III provided a value of 1.15 mm. Method I resulted in a much lower value (0.7 mm) than the experimental result. A similar trend also occurs at the symmetrically loading point of W4.

When the load was applied to W3, the maximum experimental deflection occurred at P4 (1.05 mm), which was under the loading point W3; the deflection was between 0.43 mm and 0.56 mm at P3 and P5, which were one joist away from the loading point. The influence on the points P2/P6 and P1/P7, which were more than three joists away, was negligible (Figure 12). For the modeling aspect, a similar trend was observed for methods II/III; the largest deflection obtained from method II was 0.98 mm (1.14 mm for method III) at P4 (Figure 13). However, the assumption of pin connection in method I results in a significant underestimation of the deformation of 0.7 mm. At the adjacent joists, method III (0.42 mm) provided results that were closest to those of the experiment at P3 and P5 (0.43 mm-0.56 mm), followed by method II (0.37 mm) and method I (0.27 mm). In general, under a static concentrated load of 1 kN, the deflection at the midspan of the joist reached the maximum value of about 1 mm, which affected the deflection of the adjacent joist. The influence on other joists that were three spaces away was minimal at about 0.1 mm, and the effect on the joists at the further distance was negligible.

Figure 14 shows the experimental and model results of the deflection at P1 to P7 when a concentrated load of 1 kN was applied to L1 and L5 successively. Figure 15 shows the respective results for applying a load to L2 and L4 successively. The results showed a good agreement between the experimental and model results. When the load was applied to L1 or L5, the measured deflection at P4 in the middle of the joist was about 0.3 mm and the deflection at the adjacent joists was less (0.15–0.25 mm). Method III provided a value that was closest to the experimental results at P4 (0.27 mm), followed by method II (0.23 mm) and method I (0.19 mm) (Figure 14). When the load was applied to L2 or L4, the measured deflection at P4 was 0.65–0.7 mm; the model results show close values of 0.79 mm for method III and 0.67 mm for method II at P4. However, method I underestimated the values, with relatively high errors (Figure 15).

Although the distances between the loading points L1, W1, L5, and W5 and the monitoring point P4 at the center are the same, loading on L1 and L5 resulted in a larger deflection at P4 than loading on W1 and W5. The same trend was observed for L2, L4, W2, and W4. The reason is that, under a static concentrated load of 1 kN, the loading load is directly transferred to the same loading joist and causes the deflection, whereas in perpendicular to the joist direction, the loading force applied to the adjacent joists significantly decreases with the increasing distance from the loading point, and the influence on the joists that were three spaces away was minimal. The metal plate connection of the joists provides linear-elastic stiffness, resulting in a large recoverable deflection because the connection is not rigid. In addition, the nailed connection between the sheathing and the joists is flexible and results in larger deformation than that of the model II which only considers the flexible metal plate connections. Method III provides the best overall accuracy, followed by method II and method I.

4.2. Experimental and Numerical Analysis of Vibration Frequencies and Modes
4.2.1. Vibration Mode and Frequency Testing

The modal frequencies and shapes of the floor system were obtained by the dynamic test according to ISO-18324-2016 [42]. Equally distributed node points 7 × 15 = 105 forming a grid were plotted on the floor surface to serve as excitation points, and five accelerometers (INV9828) were used to measure the vibration acceleration of the floor (to measure the high vibration mode, accelerometers were not installed at the center of the floor system), as illustrated in Figure 16. Figure 17 shows several kinds of equipment used for the vibration performance testing, including the ICP hammer (INV9314), data acquisition instrument, and accelerometer (INV9828). A person sitting on a movable wood frame moved from one side to the other side of the floor and used an instrumented ICP hammer to hammer the floor and created impact excitation. The distribution of the excitation points was displayed in the DASP-V10 Software on the computer screen. The vertical vibration movements of the measured points were transformed into electrical signals and continually recorded by the data acquisition, and the electrical signals were processed based on fast Fourier transform (FFT) to obtain the frequency response function (FRF). The modal frequencies, shapes, and damping ratios of the tested wood truss joist floor were extracted from the FRF.

4.2.2. Comparison of the Measured and Modeled Vibration Mode and Frequencies

Laboratory model testing permits the unambiguous description of the floor boundary conditions. It is very difficult to model the actual support conditions of timber floor structures because the floor system is connected to the wall by screws or nails. In this study, the numerical results indicated that the simplification of the fixed support between the edge joists and the walls (AB/CD) and the simply supported connection between the other joists and the walls (AD/BC) provided closest vibration modes to those of the experiments. The modeling results of three modeling methods are consistent with the tested results for the first three-order mode shapes, illustrated in Figure 18. Vertical vibration throughout the floor is observed in the first-order mode. Opposite vibrations in the two parts in the width direction (perpendicular to the joist) were exhibited in the second-order mode, whereas fluctuations for three parts in the width direction were shown in the third-order mode.

Table 5 lists the details of the first three-order frequencies obtained from the experiment and the different numerical models. As to serviceability design, the fundamental frequency of a timberfloor in residential buildings is required to be larger than 8 Hz. The hammer impact test indicated that the first vertical mode of the vibration had a frequency of 15 Hz.

Methods I/II/III provided different values of the fundamental frequency of the floor. Method I, which ignores the effects of the flexible connections, provided the highest fundamental natural frequency (20.12 Hz) with great deviation up to 34.3%. Method II only considers the flexible characteristics of the metal plate connections, which resulted in a slight overestimation of the vibration stiffness of the floor and a high value of the fundamental frequency (17.73 Hz) with the difference of 18.2%. Method III (16.9 Hz) with small difference of 12.7% was closest to the results of the experiment because the method considers the flexible characteristics of the metal plate connections of the joists and the nailed/screwed connections of the sheath to the joist; and the linear-elastic stiffness of the floor was simulated. However, the modeling results of the frequency were slightly higher than the experimental results. The reason for the error may be that the fabricated wood truss joist floor had a gap in the connection and there was some variability of the natural wood materials and simplification boundary of numerical models.

Method III also provided good agreement with the second-order frequency (17.7 Hz) with the difference of less than 1% versus the tested frequency (17.8 Hz). As to the fluctuations for three parts in the width direction, method III predicted that the vibration mode occurred at the third-order frequency of 18 Hz and the sixteenth-order frequency of 20.85 Hz. There are many local mode occurrences at the frequency intervals between 18 Hz and 20.85 Hz. In the laboratory test, the third-order frequency of 21 Hz was provided, which is close to the modeling sixteenth-order frequency (20.85 Hz). There are some differences between the modeled and the tested third-order frequency. The test may miss the vibration modes at tiny frequency intervals between 18 Hz and 19 Hz.

Damping is an intrinsic structural property and represents the ability of the floor system to absorb and dissipate vibrational energy. Damping cannot be calculated but is typically determined by a test. Based on the average measured damping ratios for the first three-order modes, a damping ratio of 6.6% or simply 5% is used for the vibration analysis of the wood truss joist floors in the subsequent analysis.

4.3. Numerical Prediction of Vibration Induced by Walking

Excessive vibration of a lightweight floor structure may occur when a person walks on the floor, resulting in a vibration serviceability problem. In this part, we discuss vibration performance of the floor obtained from different modeling methods. The responses of different floor models to a walking load were compared to optimize the model for a future vibration serviceability study.

4.3.1. The Human Walking Model and Paths

We used the single-step load function proposed by Chen and Li [43, 44]; a person with a mass of 90 kg and 2 Hz walking frequency was considered, and a stride length of 600 mm and 5% damping were assumed, in Figure 19.where Fs (t) represents the single footfall force, An stands for the load factors of each order, Ts is the period of a single footfall, fs represents the walking frequency, and m is the pedestrian mass.

The cases of a person walking perpendicular to the joist direction (Figure 20) and parallel to the joist direction (Figure 21) were simulated. Owing to the stiffness inhomogeneity and differences of the wood joisted floor, the deformation and acceleration responses were predicted and investigated at two points near the floor center. The monitoring point A is located on the sheathing at the center of the middle joist, and the monitoring point B is located at 100 mm in the vertical direction from A, as shown in Figures 2021. The use of different measuring points was required to obtain an accurate assessment of the dynamic response of the wood truss joist floor because there may be installment errors of the acceleration transducer in the laboratory test.

4.3.2. Comparison of the Vibration Responses of the Different Floor Models

The vertical deformation-time curves at points A and B for the simulation of a person walking perpendicular to the joists (W line) are presented in Figure 22. It was found that method III provided larger displacement values than methods II and I, and the differences between the methods were greater at point B than those at point A. The maximum vibration deflection occurred at the center of the floor, and the walking perpendicular to the joists at the far place produced little influence on vibration deformation at the floor center, where the curves showed almost horizontal line for a relatively long time period. There were significant differences in the modeling results for the three methods at the peak; low stiffness was observed at point B, and the peak for method III was 2.3 mm, whereas that of method II was 1.8 mm and method I was about 1.6 mm. The reason is that the sheathing-to-joist screwed connection has linear-elastic stiffness, resulting in larger vibration deflection than the model that only considers the metal plate connection, especially for the vibration between the joists. Without consideration of elastic stiffness of flexible connections, method I greatly underestimated the vibration amplitude and may lead to high errors.

The vibration acceleration-time curves at points A and B for the simulation of a person walking perpendicular to the joists (W line) are shown in Figure 23. The acceleration was larger at B than at A (about 7.58 m/s−2 at B and 1.45 m/s−2 at A for method II; 9 m/s−2 at B and 1.62 m/s−2 at A for method III; 3.71 m/s−2 at B and 1.34 m/s−2 at A for method I). That is because the acceleration is higher for the slender place (between the joists). Three models resulted in the peak acceleration at 2.6 s at A and 2 s at B.

Table 6 lists the peak acceleration and displacement at A and B for methods I, II, and III. Method III was used as the reference method for calculating the difference. At point A, method II, which only considers the linear-elastic characteristics of the metal plate connections, resulted in an underestimation of the peak acceleration (90%) and the peak displacement (86%). The simplification of pinned connections for metal plate connections and rigid connections for sheathing-to-joist connections from method I directly results in large underestimation (83% at the peak acceleration and 81% at the peak displacement). At point B, method II provided an underestimation of the peak acceleration (84%) and the peak displacement (76%); method I lead to large deviation at the peak acceleration (41%) and the peak displacement (70%). It is concluded that the linear-elastic characteristics of the metal plate connections directly affect the dynamic response of wood truss joist floors; meanwhile, the linear-elastic features of the nailed/screwed connections to the sheathing cannot be neglected in a dynamic analysis; the effect is greater between the joists than on the joist. The flexible characteristics of connections are more susceptible to great vibration response.

The vibration responses of the floor at A and B when simulating a person walking parallel to the joist (at a vertical distance of 100 mm between the L line and the middle joist) are shown in Figure 24. This case is considered because a person may not always walk on the middle joist. Although the force of the footfall is on the sheathing and not on the joist, the load is transferred to the nearest joist (the middle joist) and causes vibration of the entire joist. It was found that method III provided larger fluctuations in the displacement than method II and method I produced the lowest fluctuation. Since A was located on the sheathing and just on the joist, there were slight differences between method I, method II, and method III; the peak displacement was 0.96 mm for method II, 1.1 mm for method III, and 0.91 mm for method I. At B, the peak displacement was 2.02 mm for method III, 1.5 mm for method II, and 1.34 mm for method I. Because point B was located on the sheathing between the joists, the displacement was larger because of the linear-elastic nailed connection.

Figure 25 shows the acceleration response at A and B when simulating a person walking along the L line with eight footfalls. The load was transferred to the nearest joist (which is the center joist) and excited the vibration of the center of the joist, with a peak acceleration of more than 1 m/s−2 at the start. At A, method I provided the lowest peak acceleration (1.73 m/s−2) and method II resulted in lower peak acceleration (1.87 m/s−2) than method III (2.13 m/s−2). Since B was on the sheathing and not on the joist, three models exhibited relatively large differences in the acceleration, with a peak acceleration of 4.7 m/s−2 for method II, 6.44 m/s−2 for method III, and 2.42 m/s−2 for method I.

Table 7 summarizes the peak acceleration and displacement at A and B for methods I, II, and III; method III was used as the reference method for calculating the difference. Method II resulted in an underestimation of the peak acceleration (88%) and peak displacement (87%) at A; method I led to greater underestimations of peak acceleration (81%) and peak displacement (82%). At B for peak acceleration, the underestimations of 73% for method II and 38% for method I were demonstrated; for peak displacement, the underestimations of 74% for method II and 38% for method I were given. Thus, in human-induced vibration analysis, the linear-elastic stiffness of the sheathing-to-joist nailed/screwed connections and metal plate connections cannot be ignored. Method III is the preferred method for floor vibration analysis because all factors are considered, including the linear-elastic stiffness of the metal plate connections of the joists and the sheathing-to-joist connections.

5. Conclusion

Numerical simulations and laboratory tests were conducted to investigate vibration performances of wood truss joist floors. The three floor models were developed to predict the vibration modes/frequencies, static response under a concentrated load, and human-induced vibration response of a full-size floor. The response of the floor models was compared with the tested results to optimize the model for a future vibration serviceability study.

Method I neglecting the flexible connections will result in the obvious underestimations of predictions of static deformation and vibration acceleration and deformation of the wood truss joist floor. Neglecting the flexible characteristics of the sheathing-to-joist fasteners connections (method II) may result in distinct underestimation of prediction of vibration response, especially at the slender place of the wood truss joist floor. The flexible metal plate connections have a significant influence on the vibration performance of wood truss joist floors, while the flexible characteristics of the sheathing-to-joist nailed/screwed connections contribute nonnegligible amount on the vibration performance.

For method III, metal plate connection is simplified by orthogonal springs between the composite force node for multiweb members and the chord node; sheathing-to-joist connection is also simplified by orthogonal springs to quantify the elastic stiffness. It is valid to predict a series of static response under concentrated loading, vibration frequencies, and modes of the full-size floor tests. It also considers all factors and provides the best performance for the vibration prediction of the wood truss joist floor.

For wood truss joist floors, the vibration induced by walking more easily affects the deflection of the adjacent joist; however, the deformation influence on other joists that are three spaces away is minimal. However, different walking ways produced steep vibration acceleration fluctuations at the floor center for a relatively long time period. The deformation and acceleration are much higher for the more slender place (between the joists) compared to the sheathing on the joists. For the future study of vibration performance of the wood truss joist floor, comprehensive assessments are needed for the place on the joist and the place between the joists.

In general, this study presents a feasible and accurate numerical method to predict the vibration performance of wood truss joist floor systems. Further research will be conducted on the effect of strong backs, thickness of subfloor, and other geometric configurations of floors on the vibration performance of the wood truss joist floor system so as to improve the predicting method. The results will be reported when available.

Data Availability

The data in this paper are given in the tables and figures within the manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

This project was supported by the National Natural Science Foundation of China (31770603 and 52008011) and Science and Technology Project of Beijing Municipal Education Commission (KM201910005019).