Issue 
Acta Acust.
Volume 8, 2024



Article Number  11  
Number of page(s)  18  
Section  Building Acoustics  
DOI  https://doi.org/10.1051/aacus/2024001  
Published online  26 February 2024 
Technical & Applied Article
Vibration and impact sound properties of hybrid steeltimber floor structures
^{1}
Chair of Vibroacoustics of Vehicles and Machines, Department of Engineering Physics and Computation, TUM School of Engineering and Design, Technical University of Munich, Boltzmannstr. 15, 85748 Garching n. Munich, Germany
^{2}
Laboratory of Soft Matter and Biophysics, Katholieke Universiteit Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium
^{*} Corresponding author: bettina.chocholaty@tum.de
Received:
31
October
2023
Accepted:
19
January
2024
Lightweight floor structures, such as timber or hybrid timber floors, face challenges associated with excessive vibrations and elevated levels of lowfrequency impact sound. Especially here, accurate prediction of a floor’s vibration and acoustic behavior is essential. However, typical laboratory testing of building elements is costly and timeconsuming. To reduce costs, in this study, adapted simulations are carried out on two types of hybrid steeltimber floor structures to evaluate vibrations and impact sound. The hybrid elements are made of laminated veneer lumber as the top and bottom layers and a trapezoidal steel component as the web. Vibration measurements are used in combination with Bayesian optimization to efficiently calibrate Finite Element models, which are subsequently utilized to quantify and validate the floor structures regarding vibrations and impact sound. The two types of crosssections, i.e., closed and open, are investigated and compared. The impact sound pressure level computations reveal promising results in predicting the behavior of the hybrid structures. However, further countermeasures are required to fulfill vibration serviceability requirements.
Key words: Hybrid steeltimber floor / Vibrations / Impact sound / Laminated veneer lumber
© The Author(s), Published by EDP Sciences, 2024
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
With the growing significance of carbon neutrality and environmental sustainability, the construction industry is actively developing sustainable solutions. Research on sustainable construction has recently expanded beyond conventional materials to encompass hybrid structures that incorporate timber in conjunction with steel [1–3], and concrete [4, 5]. These hybrid building elements have led to a range of investigations involving aspects such as vibration characteristics [1, 2], economic viability, environmental sustainability [3], and safety performance [6]. One notable area of concern in these modern hybrid structures, particularly in the context of extended spans and slender designs for lightweight floors such as wood and hybrid wood floors, evolves around their susceptibility to low resonance frequencies and annoying vibrations, which can potentially cause discomfort among occupants [7, 8]. Consequently, a significant challenge in the context of largespan timber floors lies in the exceedance of allowable vibration limits, thus making vibration serviceability a paramount criterion for design [9].
As an example, research by Hassan et al. [10] illustrates that in terms of deformations, CrossLaminated Timber slabs with spans of 7 m and relatively modest thickness can be accommodated. However, issues arise as natural frequencies dip below 8 Hz for spans of approximately 4.5 m [10]. In such scenarios, composite or hybrid timber structures emerge as potential solutions. Perković et al. [11], for instance, determined a natural frequency of approximately 10 Hz for a hybrid concretetimber floor with a 7m span.
Furthermore, lightweight structures like timber or hybrid timber floors often encounter challenges with impact sound insulation [12]. Typically, building standards and regulations rely on laboratory testing for assessing the suitability of floor structures. However, these laboratory tests, especially those concerning impact sound insulation, are timeconsuming and costly [13, 14] and are conducted in certified testing facilities. In the context of vibrational characteristics, either calculations or certified measurements are imperative [15]. Therefore, when adequately executed, simulations have emerged as an appealing alternative to laboratory tests.
Simulationbased engineering design often employs numerical models such as those based on the Finite Element Method (FEM), a tool extensively used by engineers across various disciplines. In practice, there tends to be a discrepancy between FEM predictions and actual test results, attributed to model form errors and approximation errors, both considered types of epistemic uncertainty [16]. To minimize the discrepancy between simulations and measurements, model parameters are iteratively adjusted, a procedure referred to as model updating. This iterative process reduces approximation errors and enhances both the reliability and accuracy of the model. The procedure involves the identification of model parameters by minimizing a problemspecific objective function [17]. Various approaches have been proposed for this purpose, whereas the approach of Bayesian optimization, in particular, offers a means to minimize the need for extensive model evaluations [18]. Its efficiency results from employing Gaussian process regression and an acquisition function for an iterative search of optimal model parameters [19].
Furthermore, recognizing the challenges presented by vibroacoustic properties of hybrid steeltimber floors, recent research [2, 20, 21] has explored innovative floor elements. In [2], the analysis focuses on a structure comprising a wooden plate combined with an Hshaped steel beam, primarily evaluating its vibration serviceability. In [20], a modular hybrid steeltimber floor constructed by Cross Laminated Timber and a Ushaped steel beam is studied experimentally. Meanwhile, [21] delves into the vibrational behavior of a hybrid steeltimber structure featuring top and bottom laminated veneer lumber (LVL) plates and a trapezoidal steel web floor element in between, with a particular emphasis on the influence of joints on overall vibroacoustics.
However, despite these valuable contributions, there remains a need for a comprehensive assessment of the acoustic properties of hybrid timbersteel structures, which are relevant for the applicability of the building elements in practice. Moreover, the structure examined in [21] has yet to be subjected to an investigation regarding its vibration serviceability. This study addresses this gap by studying both vibration serviceability and impact sound levels of the hybrid floors presented in [21]. To achieve this objective, the applied approach involves the calibration of a Finite Element model (FE model) using Bayesian optimization techniques [22]. The Bayesian optimization framework is chosen since it requires fewer model evaluations for a minimization than commonly used approaches by applying Gaussian process regression [18]. The calibration process is based on data derived from experimental vibration analysis, which allows the identification of critical parameters associated with the joints within these hybrid structures. Subsequently, using the calibrated simulation model, vibration limits and impact sound criteria are evaluated, providing valuable insights into the applicability and performance of these innovative floor elements. As a result, this research offers a comprehensive numerical exploration of the vibroacoustic behavior of a novel LVLtrapezoidal steel floor design. These findings, in turn, offer initial assessments for the prospective integration of these building components into sustainable construction practices.
2 Hybrid steeltimber floor structures
This study numerically analyzes hybrid structures composed of LVL and trapezoidal steel frames. In detail, it is constructed by top and bottom LVL plates with a trapezoidal steel core in between. The structural components are connected by fasteners whose stiffness values are unknown in advance. The numerical analysis is conducted using FEM. First, a calibration of the unknown joint parameters p is performed utilizing numerically determined natural frequencies as well as mode shapes and results from experimental modal analyses on smallscale test specimens, i.e., S1* and S2*. Subsequently, the updated joint parameters are used in an FE simulation to analyze largescale floor specimens, i.e., S1 and S2, regarding vibration dose values (VDV) and impact sound pressure levels. The impact sound pressure level simulations use a numerically implemented tapping machine load. The overall workflow is visualized in Figure 1.
Figure 1 Overall workflow of the study: calibration of unknown model parameters p is performed utilizing numerically determined natural frequencies as well as mode shapes and results from experimental modal analyses on the smallscale test specimens, S1* and S2*. The updated parameters are used in a subsequent Finite Element simulation to analyze the largescale floor specimens, S1 and S2, regarding vibration dose values and impact sound pressure levels. The impact sound pressure level simulations use a numerically implemented tapping machine load. 
2.1 Description of investigated structures and simulation models
Two configurations are investigated: a closed (S1) and an open (S2) crosssection. Snippets of the test specimens are visualized in Figure 2 together with their geometric properties, which are also listed in Table 1.
Figure 2 Snippets of closed (left, S1) and open (right, S2) hybrid LVL and trapezoidal steel elements. The lower timber plate is built of a continuous lower wooden panel for closed crosssections and made of several separate panels for open configurations. The respective geometrical parameters are given in Table 1. 
According to ISO 101405 [23], the recommended size of a floor structure, which is to be analyzed concerning impact sound, is given as 10–20 m^{2}, and the shorter edge should be greater than 2.3 m. Hence, the dimensions of 3.6 m × 2.8 m are chosen for the hybrid floors (Tab. 1). However, due to the testing facility’s limited spatial capacities, it has been impossible to carry out vibration measurements on the entire floor structures as given in Table 1. Hence, smaller specimens are used for the model validation and calibration, i.e., a closed S1* and an open crosssection S2*.
The main difference between the open and closed crosssections lies in the design of the lower timber plate, which is either a continuous wooden panel for closed crosssections or made of several separate panels for open configurations. Based on preliminary structural investigations, it has been determined that the upper screw connection may be too weak in terms of loadcarrying capacity and that nails would be beneficial for the upper joint. However, as the nail connection can only be made from the steel side, the open variant has been designed. The respective mounting direction of the fasteners can also be seen in Figure 2. Obviously, the different geometry, fasteners, and materials used for both types of crosssections influence the dynamic behavior, e.g., the modal properties, of the investigated systems, which is why the two types are subsequently compared concerning VDV and impact sound insulation.
In the case of the closed form, both LVL plates are built of the material BauBuche^{1} Type Q (BBQ). For the open crosssection, the upper plate is also made of BBQ whereas the lower one is made of BauBuche Type S (BBS) [24]. The respective material properties provided by the manufacturer [24] or taken from literature [25] are given in Table 2.
The numerical studies use the commercial FE software ANSYS [26]. For the FE model, the structural components are represented by quadratic shell elements (SHELL281) with the provided orthotropic material characteristics for the timber parts and isotropic material properties for the steel sheets. Both specimen crosssection types use a trapezoidal core using steel with a nominal yield strength of f_{y,k} = 320 N/mm^{2}. However, the thickness of the steel core differs for the open, d_{S} = 1.5 mm, and closed, d_{S} = 0.75 mm, configuration. Furthermore, comparing the open and closed crosssections, the fastener type varies due to the manufacturing process of the building elements. For the open configuration, nails^{2} join the steel plates with the upper and the lower LVL plate. In the closed case, nails are only applied to join the lower LVL plate to the steel core. The upper LVL plate is attached through screws^{3}. The respective positions of the fasteners are visualized in Figure 2.
For all fasteners, the ANSYS element type “MPC184 general joint” is utilized. These elements allow specifying a stiffness matrix with 21 entries:
Here, the indices “1”, “2”, and “3” denote displacements in the x, y, and zdirection, respectively. Moreover, the indices “4”, “5”, and “6” refer to rotational degrees of freedom (DOFs) around the x, y and zaxis, respectively. The coordinate directions are visualized in Figure 1. A preliminary analysis conducted as part of this work revealed only three stiffness parameters, i.e., d_{11} = d_{xx}, d_{22} = d_{yy}, and d_{12} = d_{xy}, that significantly influence the vibrational behavior of the structure. Hence, all other entries of the fasteners’ stiffness matrix are set to zero.
The concept for the fastener model is depicted in Figure 3. In addition to the general joint elements, coupling constraints for the displacement in the zdirection are applied to closely located nodes at the interface of steel and timber surrounding the joints. The surrounding of the joints means within a radius of 0.022 m. This radius is determined through an empirical evaluation of the parameter space and leads to the best match of simulation and experimental results with the fixed value. In this surrounding area, nodes that are coincident within a certain tolerance behave as infinitely stiff coupled related to the displacement in zdirection. The tolerance is set just large enough to cover the distance between steel and timber components. This way, pairs of steel and timber component nodes are coupled concerning their displacement in zdirection.
Figure 3 Fastener model: an MPC184 general joint element (represented by a black line) using a stiffness matrix between nodes related to fasteners on the substructures (grids). Coincident nodes at the interface of steel and timber are coupled related to the zdisplacement (green lines) within an area surrounding the fastener (green circle). 
Since the open crosssections only use nails, denoted by the index “n”, the same joint stiffness parameters are assumed for the upper (u) and lower (l) joints d_{ij,l} = d_{ij,u} = d_{ij,n}. For the case of closed configurations, different joint parameters for upper and lower connections are set, i.e., d_{ij,u} = d_{ij,s} and d_{ij,l} = d_{ij,n}, where the index “s” refers to properties related to screws. As the fastener models’ parameters are initially unknown, they constitute the unknown model parameters p, identified using model calibration and vibration measurements in the first analysis step.
The FE model uses a mesh size of l_{ele,t} ≈ 0.015 m for the timber components and of l_{ele,s} ≈ 0.013 m for the steel core. In [27], six to eight quadratic elements per wavelength have led to an acceptable error in eigenfrequencies compared to a converged solution. Hence, the bending wavelengths λ [28] are computed using the respective material properties of Table 2. The comparison of the sixth of the timber’s bending wavelength λ/6 with the mesh size shows that the finite element size of the timber components is valid up to 10,000 Hz, which is sufficient for the current studies. In the case of the steel parts, the used mesh size accounts for six to eight quadratic elements per wavelength in the lowfrequency range. However, the mesh size is only valid up to 1300 Hz, for S1 with a steel thickness of d_{S} = 0.75 mm, and up to 2400 Hz, for S2 with a steel thickness of d_{S} = 1.5 mm. Due to the limited computational resources, it becomes impossible to reduce the mesh size further. The current model for the vibroacoustic evaluation already uses approximately 0.9 million DOFs. Hence, the simulations are performed until the specified maximum frequency for the respective cases.
2.2 Model calibration
For the model calibration, the closed and open samples S1* and S2* are used. The inplane dimensions of the test specimens are changed to 3.50 m × 0.84 m due to the limited spatial capacities of the testing facilities. All other properties remain as specified for the largescale samples in Table 1.
2.2.1 Numerical approach
The FE model of the respective sample S1* or S2* is calibrated by iteratively adapting the model parameters p, i.e., the joint stiffnesses values of the nails d_{ij,n} and the screws d_{ij,s}, and computing an error e using the numerical and experimental results. An overview of the model calibration workflow is visualized in Figure 4.
Figure 4 Flowchart of the model calibration: An FE model of the respective measured sample S1* or S2* is calibrated by iteratively adapting the unknown parameters p and computing an error e using the numerical and experimental results. 
Six parameters in the FE model calibration are unknown at the beginning:
The closed specimen uses all six of them, whereas the open specimen utilizes only the three parameters related to nails, i.e., d_{ij,n}. First, the calibration of the nails’ parameters is conducted using the data of the open specimen. Subsequently, those resulting nail parameters are used in the model of the closed specimen. Second, the screws’ parameters are calibrated utilizing the model and the measurement data of the closed specimen.
The free boundary conditions and the load position are adopted in the model in correspondence with measurement setups. Furthermore, a frequency resolution of Δf = 1.0 Hz is applied in the frequency range of interest f_{low} ∈ [10, 312.5] Hz, which covers the measured frequency range.
2.2.2 Experimental approach
Vibration measurements are performed on the two test samples S1* and S2* to identify modal properties through experimental modal analysis [29]. Therefore, a pseudorandom excitation is applied to the specimens by means of a B&K^{4} modal exciter type 4284 connected to the test samples utilizing a stinger at one point on the specimen. Heavyduty slings are used to suspend the test samples in the measurement setup (Fig. 5) to simulate free boundary conditions in the FE model.
Figure 5 Scan setup using heavyduty slings for suspension with Laser Doppler Vibrometer (LDV) (a) and shaker (b). 
A force transducer (B&K Force Transducer Deltatron Type 8230) measures the force applied to the specimens. Moreover, a scanning Laser Doppler Vibrometer PSV500^{5} (LDV) records the surface velocity on the specimen’s side opposite the force application position on a distributed grid of scan points. The measured force and velocity data are combined to compute the frequency response function. A Fast Fourier Transformation transforms the recorded time data to the frequency domain. The resulting spectra are averaged over 25 measurements for each scan point in the complex plane. Furthermore, the signals are adapted using a rectangular window. The sampling rate f_{smp} ≈ 781 Hz and the frequency resolution of Δf ≈ 390 mHz are applied. All data acquisition tasks utilize the Polytec software. Further details on the measurements are given in [21].
After performing experimental modal analyses, the experimentally determined modes visualized in Figure 6 for the closed specimen S1* and in Figure 7 for the open specimen S2* are used to calibrate the model.
Figure 6 Mode shapes for the sample S1*: experimentally determined (left) and numerically simulated (right) after model calibration. 
Figure 7 Mode shapes for the sample S2*: experimentally determined (left) and numerically simulated (right) after model calibration. 
2.2.3 Parameter optimization
The respective stiffness values of fasteners are calibrated using experimental modal analysis data and the FE modal analysis results of the samples S1* and S2*. A multiobjective model updating [30] is used, which considers a frequency error e_{1} and a mode shape error e_{2}:
where f_{exp} and f_{fem} are the experimentally and numerically determined natural frequencies, respectively, and MAC_{m} is the modal assurance criterion of the mode shapes computed, as described in [31]:
Here, N_{m} equals the number of eigenfrequencies and modes considered in the optimization and is assigned as N_{m} = 6 as six modes per specimen, i.e., S1* and S2*, are taken into account. A multiobjective optimization is chosen instead of a singleobjective optimization to circumvent the choice of weightings of the errors e_{1} and e_{2}. The optimization goal is to identify the best solution in the Pareto optimal front [32]. Hence, based on the evaluations of the objective functions of equations (4) and (5), the Pareto optimal front, e = (e_{1}(p), e_{2}(p)), is computed. For conflicting objectives, a singular optimal solution does not exist. Instead, there is a collection of alternative solutions referred to as Pareto optimal solutions. These solutions are optimal because no other solutions within the parameter space surpass them when considering all objectives [33]. Thus, a criterion is established to determine the most suitable solution at the Pareto optimal front by summing both errors for each point on the Pareto optimal front (Fig. 8) and taking the set of parameters, which leads to the smallest sum of errors as the best solution:
Figure 8 Pareto optimal front of the multiobjective optimization for the closed S1* (a) and open S2* (b) specimen. 
As an optimization method, Bayesian Optimization is utilized due to its benefits explained hereafter. First, regarding the number of function evaluations needed, Bayesian optimization is a highly efficient method. Furthermore, it performs well even when the objective function contains multiple local maxima. The method can integrate prior beliefs concerning the problem, which aids in guiding the sampling process and balancing the exploration and exploitation of the search space. As it relies on Bayes’ theorem, the prior knowledge embodies our initial beliefs concerning the range of potential objective functions. Despite the unknown nature of the cost function, it is justifiable to assume that there is prior knowledge about specific properties, such as smoothness, which renders certain objective functions more credible than others. As observations are accumulated, the prior distribution and the likelihood function are merged, leading to updated beliefs in the posterior distribution. For efficient sampling, Bayesian optimization employs an acquisition function to find appropriate sampling locations. This decision inherently involves a tradeoff between exploration (in areas where the objective function is highly uncertain) and exploitation (sampling values where the objective function is anticipated to be high) [34]. Due to its goal to minimize the number of objective function evaluations and, thus, to reduce the computational effort, Bayesian Optimization is applied in this study.
Since the parameters for the fasteners are initially unknown, a manual search is performed first by picking and testing reasonable values. Using the best choice for the joints’ stiffness values d_{ij,s} and d_{ij,n} of the manual search, an already acceptable match is achieved, as is visible in the columns denoted by “before” in Table 3. The label “before” is assigned since these results are computed utilizing manually found parameter values. The label “after” refers to results calculated with the parameters found by the Bayesian optimization.
Numerically f_{fem} and experimentally f_{exp} determined natural frequencies [Hz], relative frequency error and MAC before and after the automatic model calibration: m ∈ [1, 6] for S1* and S2*.
The initial joints’ parameters, as identified by the manual search, are given in Table 4 denoted by p_{0}. The properties in Table 4 refer to the values found by the automatic model calibration performed using Bayesian optimization.
Initial p_{0} and calibrated parameters for the closed S1* and open S2* specimen: u and l denote the upper and lower joint, respectively. All values in the table refer to the properties of nails except for the entries marked by †, which refers to screw parameters.
The optimization, i.e., the automatic model calibration, utilizes the manually found parameters as guiding values for the bounds specified in the calibration procedure. Narrow, i.e., ±25% p_{0}, and broader, i.e., −90%/+1000% p_{0}, bounds are applied to thoroughly search the parameter space. Subsequently, the pairs of errors at the Pareto front are found, and the pair of those data points with the least sum, i.e., e = min(e_{1}(p) + e_{2}(p)) is chosen as the best solution. The respective results of the parameter search are visualized in Figure 8, together with the Pareto optimal front for both specimens.
From the optimization computations, the parameters labeled as in Table 4 and, thus, the natural frequencies and MAC values denoted by “after” of Table 3 result.
Figure 8 shows that two groups of error combinations are formed for the test specimen S1*, while there is only one for S2*. For S1*, one group of model solutions results in a smaller frequency error with almost equal mode error compared to the other group. Hence, a better solution is found for the group with the smaller frequency error. This aspect can also be observed in Table 3. The relative frequency error before the automatic model calibration lies above 10% for modes 1, 3, and 6, which is strongly reduced to errors below 10%, i.e., e_{r,m} ∈ [0.02, 0.07], for all modes except the third, where it lies slightly above, i.e., e_{r,3} = 0.12. This improvement in frequency error comes with the price of reduced MAC values for modes 3, 5, and 6, i.e., MAC_{m} ∈ [0.74, 0.80]. However, the relative reduction of MAC values before and after the automatic updating lies below 0.03 for modes 3 and 5 and is only exceeded for mode 6 with a MAC value reduction of 0.08. Generally, the MAC values for specimen S1* are not as good as for the open specimen S2*. Only the MAC for the first mode lies above 0.90. For modes 2 and 5, the MAC values are still above 0.80, but the values are even lower for modes 3, 4, and 6. The reason for this behavior is presumably the shape of the modes. As visualized in Figure 6, only mode 1 shows a nice and clear behavior. Modes 2–6 are either challenging to identify, nonsymmetric, or both. This circumstance could stem from manufacturing inaccuracies, e.g., slight relative rotations of the components, which lead to a nonsymmetric builtup, prestress, and predeformation. Furthermore, sample S1* uses a steel sheet of thickness d_{s} = 0.75 mm as opposed to d_{s} = 1.5 mm for S2*, leading to individual parts of the steel frame vibrating locally, which influences the global structural behavior. Moreover, the modes are challenging to identify in the experiments, e.g., due to closely spaced modes. Consequently, due to the overall improved agreement between numerical and experimental results, the updated parameters are deemed acceptable.
For S2*, in Figure 8b, only one group of model solutions exists, and the Pareto front is discernible. The Pareto criterion assesses whether a state improves by changing one target value without deteriorating another. The Pareto front shows the values that represent the best compromise. Here, no solution with a smaller frequency error could be found without causing an increase in the mode error. Due to this behavior, only a slight improvement of frequency and mode error is achieved for S2* as can be seen in Table 3: only the natural frequencies and become slightly closer to the experimentally determined natural frequencies. However, simultaneously, the match for becomes slightly worse, and the MAC values do not increase. In particular, the fourth mode, which appears to be less clear than the others (Fig. 7), shows some discrepancy with the simulation results regarding natural frequency. This mode also shows a MAC of only 0.61, presumably due to the complexity of the mode shape. Moreover, the natural frequencies of the first and third modes, i.e., the bending modes in the lateral (y) direction, are somewhat apart from the numerical results. However, a better match is challenging to obtain since, for an advantageous adaptation of the joint parameters, the other natural frequencies deviate more from their experimental counterparts. Still, all MAC values lie above 0.90 except for mode 4 due to a complex mode shape. Furthermore, except for the second lateral bending mode 3 and the unclear mode 4, the relative frequency error lies below 10%, which is satisfying. Hence, also for the open crosssection, the calibrated parameters are assumed to be suitable for the simulations related to the vibroacoustic validation, which is described in the following section.
The updated values of the fasteners are calibrated for natural frequencies up to approximately 115 Hz and are thus adaptable for the VDV calculations, which consider modes around and below this frequency. However, it should be noted that frequencies of up to 1.3 kHz and 2.4 kHz are examined in the context of impact sound tests. Therefore, uncertainty is added to the model by assuming that the joint stiffnesses from the updating procedure are suitable for these impact sound computations.
2.3 Model for vibroacoustic validation
FE analyses are performed to validate the vibration serviceability and impact sound characteristics of the hybrid steeltimber structures, as measurements are often timeconsuming and costly [13, 14]. Here, the floor structures S1 and S2 (Tab. 1) without additional components, such as floating floors or walls located on the edges of the floors, are considered in the simulations. Even though these assumptions result in discrepancies between the predicted and presumably observed behavior of insitu floors, the FE models are assumed to provide a first estimate of the applicability of the proposed hybrid floors.
The parameters (Tab. 4), i.e., the joint parameters identified by calibration, are utilized for these FE models. Given the profound influence of boundary conditions on structural vibration and sound radiation in the lowfrequency regime, and in light of the unavailability of measurements for samples S1 and S2 in a builtin configuration, the investigations focus on the two extreme scenarios regarding boundary conditions: namely, clamped (“c”) and free (“f”). In the clamped configuration, all degrees of freedom are constrained along the outer edges of the floor structure, whereas in the free boundary condition scenario, none are restrained. Furthermore, to increase computational efficiency and manage system matrices with respect to memory resources and disk space, symmetry boundary conditions are employed in the FE model. This approach involves discretizing and solving only onequarter of the floor element, assuming that the remaining threequarters can be derived through symmetry principles. The respective boundary conditions are applied along the cutting planes of the quartered floor structure. The load for impact sound computations is located in the center of the structure, which is the point of symmetry in x and ydirections (see Fig. 1 for the definition of the coordinate system). This excitation force only excites symmetric modes in both the x and ydirections. Therefore, applying only symmetric boundary conditions is theoretically sufficient, and no additional antisymmetric boundaries need to be considered in the response analysis.
As the VDV calculation is based on the load case of a person crossing the floor, all kinds of modes are excited and thus need to be considered. Therefore, for the computation of the VDV, symmetric, antisymmetric, and all possible combinations of support conditions are applied in the simulation to incorporate both antisymmetric and symmetric mode shapes.
The implementation of the antisymmetric and symmetric boundary conditions is performed by imposing constraints on nodes, respectively DOFs, in the planes of symmetry or antisymmetry. Hence, the following DOFs are constrained:
• For symmetry conditions

In the plane with the normal in xdirection: displacements in xdirection, rotations around the yaxis, rotations around the zaxis.
In the plane with the normal in ydirection: displacements in ydirection, rotations around the xaxis, rotations around the zaxis.
• For antisymmetric conditions
In the plane with the normal in xdirection: displacements in y and zdirection, rotations around the xaxis.
In the plane with the normal in ydirection: displacements in x and zdirection, rotations around the yaxis.
Furthermore, the global impedance for the impact load computation taken from the respective FE model is scaled due to symmetry considerations [35].
In contrast, the “c” and “f” boundary conditions are implemented along the outer edges as previously described. Hence, four cases are investigated:

Closed crosssection, clamped boundary condition: S1, c;

Closed crosssection, free boundary condition: S1,f;

Open crosssection, clamped boundary condition: S2, c;

Open crosssection, free boundary condition: S2,f.
The FE models of the analyzed onequarter of the floor structures are visualized in Figure 9.
Figure 9 FE model of closed sample S1 (a) and open sample S2 (b): bottom view and side views. Steel is colored purple, and timber parts are cyan. The element thickness is visually scaled to represent the thickness of the respective component. 
A structural modal damping ratio of 1% is chosen according to the recommendation of the European timber standard [36]. The damping identified for the test specimens S1* and S2* is deemed inappropriate for the FEmodel of the floor structures, i.e., S1 and S2, since boundary conditions significantly influence the damping of a structure [37] as well as the overall builtup. Furthermore, no structural vibration tests have been possible on largescale floor structures in a builtin situation. The chosen damping ratio is applied as a structural loss factor η_{s} = 0.01 × 2 in the FEmodel for the whole analyzed frequency ranges f_{impact,S1} ∈ [10, 1300] Hz and f_{impact,S2} ∈ [10, 2400] Hz. Here, uncertainty is introduced in the simulation model since the damping value proposed by Eurocode 5 [36] applies to low frequencies related to vibrations. A timber joist floor investigated by Wang et al. [38] showed an approximate structural loss factor of 0.02 up to 2500 Hz. Although this work deals with a pure timber instead of a hybrid steeltimber floor, the chosen loss factor of 0.02 is deemed acceptable as a first estimate, as currently, no investigations exist on the impact sound behavior of hybrid steeltimber floors such as the presently studied ones.
3 Vibroacoustic criteria
The updated and adapted FE model is used in the following sections to evaluate the floor structures, S1 and S2, concerning the vibration serviceability and impact sound insulation.
3.1 Vibration serviceability
A recent review on steeltimber composite floors [39] states AISC Design 11 [40] and CCIP016 [41] as currently the most adequate references concerning vibrations of hybrid steeltimber systems. According to [41], the regulations described in the document apply to structures made of plenty of construction materials with walking excitation. For floor structures with natural frequencies above 4 Hz, the impulsive response calculation procedure is adopted as shown in Algorithm 3.1.
The response factor can be calculated utilizing the acceleration time history from Algorithm 3.1, which is related to continuous vibration. Nevertheless, a continuous vibration is seldom the case. Consequently, the concept of VDV proposed in [40] is used to estimate the acceptability of intermittent vibration:
with a_{w,rms} the RMSacceleration response from Algorithm 3.1, n_{day} = 32, n_{night} = 16 and the time it takes to cross the floor
Consequently, FE modal analysis is adopted to identify the fundamental frequency f_{1}, the natural frequencies f_{m} of modes up to twice the fundamental frequency, and the corresponding modal masses of these modes. The modal damping ratios are taken as 1%, as recommended by [36].
3.2 Impact sound
The standard ISO 101405 [23] specifies procedures to measure impact sound properties. Here, the tapping machine is described as excitation for impact sound measurements. A numerical version of the tapping machine is used in this work. The procedure to compute the loading that is applied to the FE model is described in the following.
Brunskog et al. [43] provide a method to calculate the discrete frequency spectrum of the impact loading due to a standard tapping machine. The floor is described by a general frequencydependent drivingpoint impedance Z_{dp}. The hammer impact force f_{0}(t) = F_{0}δ(t) is applied to the structure. This yields the frequencydependent force f_{0}(ω) = F_{0} by a Fourier transform which varies between 2M_{0}v_{0} and M_{0}v_{0} for an elastic and plastic impact where v_{0} = 0.886 m/s is the initial hammer velocity and M_{0} = 0.5 kg equals the weight of a hammer. Here, the geometric mean is considered [44]. The spectrum of the continuing impact is given as
Then, an inverse Fourier transform computes the corresponding force in time . The time of the zero crossing is found and the force spectrum F_{1}(ω) is computed by a Fourier transform. Finally, the Fourier series components F_{n} = F_{1}(nf_{h})f_{h} with f_{h} = 1/T_{h} and T_{h} = 0.5 s for each of the five hammers are applied to the floor structure at the respective frequencies [43]. An exemplary plot of the tapping machine load over frequency is given in Figure 10.
Figure 10 Example plot of load F_{n} against frequency. 
Furthermore, the drivingpoint impedance of the floor is computed from the local Z_{L} and global impedance Z_{G} as follows [38]:
The global impedance
is found using the dynamic stiffness of the floor
taken from the FE model at the DOF j of the impact with M_{jj} and K_{jj} as mass and stiffness matrix entries and a material damping coefficient g = 0.0084 for the upper timber plate, which is derived from vibration measurements of the LVL plates.
The local impedance is analytically approximated [38]:
with r_{h} the hammer radius, the Poisson’s ratio , and the elasticity modulus of the respective material.
The computed force is used in the desired frequency range of f_{impact} as excitation on the structure in a timeharmonic analysis. The considered frequency range differs between open and closed crosssections due to requirements concerning finite element sizes. Hence, f_{impact,S1} ∈ [10, 1300] Hz and f_{impact,S2} ∈ [10, 2400] Hz are applied as frequency ranges in the simulations. The resolution of Δf = 2 Hz for the impact sound calculations results from the repetition frequency of the hammers of the tapping machine. To determine impact sound characteristics, the sound power P is computed from the surface velocity of the lower surfaces of the samples, i.e., timber panels and radiating steel surfaces, using the Rayleigh integral [45, 46]. The sound power level is calculated from the computed sound power P as
Moreover, thirdoctave band data is calculated from the equidistant frequency data.
Since the goal is to compare normalized sound pressure levels for the hybrid steeltimber floors, sound power levels L_{W} are converted to sound pressure levels L_{p} according to the following relation from [47]:
where the equivalent absorption surface of the room equals
A reverberation time of T_{60} = 0.45 s is considered, as specified in [48] for the respective room volume V = l · b · h_{room} = 24.3 m^{3} where the room height is assigned to be h_{room} = 2.4 m and length l and width b are taken from Table 1. Moreover, the reference surface A_{0} = 1 m^{2} and the sound wave velocity c = 343 m/s are specified. The overall surface of the room is computed as S = l · b · h_{room} · 2 + b · h_{room} · 2 = 51.1 m^{2}. Furthermore, the center frequency of thirdoctaves f [Hz] and the correction terms
as well as
are required for equation (16). The calculation of the correction terms in equations (18) and (19) uses the static air pressure p_{s} = 95 kPa for an altitude of 500 m, the reference static air pressure p_{s,0} = 101.325 kPa, the air temperature θ = 293.2 K (=20°), θ_{0} = 314 K and θ_{1} = 296 K.
Subsequently, the following relation from [23] is applied
to compute the normalized impact sound pressure.
4 Results
In this section, the findings related to vibration serviceability and impact sound properties of the samples denoted as S1 and S2 are presented. The boundary conditions play a significant role in influencing the structural vibration and sound radiation in the lowfrequency range. Unfortunately, measurements in a builtin situation have not been feasible. Thus, a comparative analysis is conducted between two extreme cases with different boundary conditions, namely, clamped (“c”) and free (“f”).
4.1 Vibration serviceability properties
The VDV is a critical metric for evaluating the potential for adverse comments in residential buildings, as prescribed by the British Standard BS 64721 [42]. The respective limit values for the VDV are given in Table 5, which are compared with the VDV calculated for all samples j ∈ [S1,f; S1,c; S2,f; S2,c] using Algorithm 3.1 and equation (8). This computation involves determining fundamental frequencies f_{1,j}, modes to be considered f_{m,j}, and modal masses through FE modeling. The results for the different sample configurations are presented in Table 6.
VDV ranges for the probability of adverse comments
Fundamental frequencies f_{1,j} [Hz], modes for the VDV evaluation f_{m,j} [Hz], and modal masses of the samples j ∈ [S1,f; S1,c; S2,f; S2,c].
Upon evaluating the VDV for the test specimens, as shown in Table 7, several observations emerge. Firstly, it is evident that S2,c exhibits a higher number of modes within the specified frequency range Δf_{S2,c} ∈ [0, 2 · f_{1,S2,c}], primarily due to its relatively high fundamental frequency f_{1,S2,c} and the presence of numerous localized modes originating within the steel frame, which is not as pronounced in S1. Secondly, S2,f has a relatively low fundamental natural frequency f_{1,S2,f}, necessitating consideration of only the fundamental mode within the observed range of Δf_{S2,f} ∈ [0,2 · f_{1,S2,f}] = [0, 10.58] Hz.
VDV during the day VDV_{d,j} [ms^{−1.75}] and VDV during the night VDV_{n,j} [ms^{−1.75}] of the samples j ∈ [S1, f; S1, c; S2, f; S2, c].
In terms of day and night VDV, samples S1,c, S1,f, and S2,c exhibit similar values, around 0.57. Notably, S2,f stands as an outlier, with a VDV reaching up to 1.67, indicating problematic vibrational behavior. The other samples fall within the “adverse comments possible” range, while S2,f exceeds this threshold and falls into the category of “adverse comments probable”. However, the behavior of a builtin floor structure, which presumably lies between the extreme boundary conditions of free and clamped, yields VDV values that also fall within a range between “possible” and “probable” adverse comments.
As specified in Algorithm 3.1, the structural response depends on the damping ratio, which has been fixed using normative specifications. Hence, uncertainty related to the damping remains and should be kept in mind when interpreting the results.
1: Compute modes and eigenfrequencies f_{m} up to twice the fundamental frequency f_{1}
2: Calculate effective footfall impulse with the maximum walking frequency f_{w} = 1.8 Hz and the considered eigenfrequencies f_{m}
3: Compute peak velocity in each mode m: with Φ_{e,m} the value of the mode shape m at the excitation position, Φ_{r,m} the value of the mode shape m at the receiving position, the modal mass of mode m
4: From this, calculate the velocity response in each mode over the period of one footfall T with ζ as the modal damping ratio: with the modal velocity magnitude
5: Compute total response to each footfall with a weighting w_{m} taken from [42] as suggested in [41]
6: For a timeharmonic analysis, the total response to each football can be converted to accelerations by computing the acceleration magnitude from the velocity magnitude as
7: From the resulting acceleration time history, rootmeansquare (RMS) response evaluated over one footfall
4.2 Impact sound properties
The results of impact sound pressure levels are examined from two perspectives: the influence of different boundary conditions, i.e., free “f” and clamped “c”, and the effect of the two distinct crosssection types, i.e., open and closed. The comparison of the impact sound results under the different boundary conditions “f” and “c” is depicted in Figure 11.
Figure 11 Normalized impact sound pressure of samples S1 (a) and S2 (b) for the cases “f” and “c” for thirdoctave bands. 
Two insights can be gained from these plots. Firstly, as anticipated, the impact of boundary conditions is most pronounced in the lower frequency range, up to 80 Hz for S1 and up to 500 Hz for S2. Secondly, the impact of boundary conditions is more substantial for the closed specimen S1 compared to the open one S2 up to 100 Hz. This discrepancy arises because, in the clamped case, all edges of the lower timber panel are fixed for S1, whereas for S2, only the edges at the longitudinal ends of the lower timber panel are fixed, resulting in a greater influence of boundary conditions in the lowfrequency range for S1.
Furthermore, a comparison of the different crosssection types, open and closed, under both boundary conditions is presented in Figure 12.
Figure 12 Comparison of the sound pressure levels of the samples S1 (closed) and S2 (open) for the cases “f” (a) and “c” (b) for thirdoctave bands. 
Both crosssection types exhibit similar quantitative behavior, with S1 displaying slightly higher sound pressure levels in some third octaves, specifically about 5 dB higher than S2. This observation suggests that S2 radiates slightly less sound than S1. Qualitatively, the two configurations share similar curve shapes, but slight shifts in characteristics are observed in frequency due to their different structural setups. For instance, under clamped support conditions, both specimens initially exhibit a low sound pressure level at lower frequencies, at approximately 30 Hz. Subsequently, a distinct peak is shown at roughly 50 Hz for S1 and approximately 63 Hz for S2. Notably, the peak for S2 shifts to the adjacent thirdoctave band due to the resonant behavior of the structure. Furthermore, qualitative disparities become evident in the clamped support scenario within the frequency range of 160–800 Hz, where S2 demonstrates more fluctuations in sound pressure levels than S1.
To provide context for the impact sound behavior of the hybrid steeltimber floors, the impact sound pressure level values from Table B.2 of ISO 123542 [49] for concrete floors (100 mm and 180 mm thick) with a bonded screed layer (20 mm and 50 mm thick) are included for comparison in Figure 13. The lower and upper bounds of sound pressure levels for the hybrid steeltimber floors are established, offering a range of insitu behavior. This means that in every thirdoctave band, the minimum value of the four sound pressure levels, i.e., S1,f, S1,c, S2,f and S2,c, is identified and taken as a value for the lower bound and vice versa for the upper bound.
Figure 13 Comparison of the upper and lower bounds of the sound pressure levels computed for the hybrid floor samples S1 and S2 with the values for a 100 mm thick concrete floor with a 20 mm thick bonded screed and a 180 mm thick concrete floor with a 50 mm thick bonded screed provided in [40]. Results are given in thirdoctave bands for the upper and lower bounds and in octaves for the concrete floors. 
In the lowfrequency range up to 200 Hz, the lower and upper bounds cover values similar to those of concrete floors. In the midfrequency range, up to around 1000 Hz, the hybrid floors’ sound pressure levels fall between the two concrete floors. Beyond 1000 Hz, the hybrid floors exhibit sound pressure levels higher than the concrete floors, with a decreasing trend at 2000 Hz. This indicates that the levels may approach those of concrete floors at higher frequencies. However, this could not be definitively verified due to the limited computational resources and, hence, limited observed frequency ranges. Additionally, it should be noted that the damping and joint stiffness values are chosen based on assumptions explained in previous sections. Thus, the interpretation of the results is subject to the respective modeling choices. In conclusion, the hybrid floors demonstrate behavior similar to concrete floors. However, influences of additional structural components of a builtin situation, e.g., floating floors, are not considered in this study. Hence, a more realistic setting in a building might give more insight into the comparison of the floor types.
5 Discussion
In the present study, the FE models of the hybrid floor structures rely on parameters derived from the calibration of FE models of smaller specimens. It is worth noting that this step assumes the transferability of these parameters, a presumption that should be validated through future largescale tests. Moreover, the construction of all models, both for small and largescale specimens, is based upon certain assumptions. These assumptions cover factors such as the appropriateness of the realization of free boundaries in the measurements and the suitability of the employed joint and structural model. The results indicate that the structural behavior has been reasonably approximated.
The assessment of vibration serviceability, as evaluated by the VDV, emphasizes concerns related to potential resident annoyance. Observably, none of the examined cases exhibit a low probability of adverse comments. One potential remedy for enhancing the performance of these floor structures involves increasing the distributed mass of the floor. However, such a mass increase necessitates accurate adjustment, as it concurrently affects the natural frequency of the floor – a factor of substantial significance for vibration serviceability as a toolow fundamental frequency may pose a challenge. Additionally, another potential strategy for mitigating the issue of high VDV involves increasing inherent damping. Investigating the damping characteristics in greater detail may provide valuable insights for effectively addressing this concern.
Regarding impact sound pressure levels, as expected, the influence of boundary conditions is most pronounced in the lower frequency range. Furthermore, the open and closed crosssections display quite similar behavior, thus making the open setup only slightly more convenient for practical applications. When the performance of the hybrid structures is compared with commonly used concrete floors with bonded screed, it becomes evident that a similar level of performance can be achieved. Still, it is to be noted that a builtin situation in buildings differs from the tested setting as usually floating floors or soft coverings are applied to further reduce impact sound noise. This step is commonly performed after the floor installation for both concrete and hybrid steeltimber floors. Compared to the concrete floors, the hybrid floors might be considered beneficial in terms of weight and sustainability.
6 Conclusions
In conclusion, this investigation explores novel hybrid steeltimber floor elements, primarily focusing on their vibration serviceability and impact sound performance. Leveraging Finite Element analysis techniques, this study provides a comprehensive examination of the behavior of these innovative building components in terms of vibroacoustics. First, the Finite Element model is validated by aligning it with experimental data on natural frequencies and modes and determining joint stiffness values. This calibrated model serves as the basis for the indepth analysis of vibration serviceability. This analysis expressed through vibration dose values, offers an initial insight into the realworld performance of these hybrid floor elements. Additionally, this study numerically explores the hybrid steeltimber floors’ impact sound insulation properties, an aspect that has received limited attention in previous research. The outcomes presented herein reveal that the proposed building elements exhibit comparable behavior to that of conventional concrete ceilings, thus making the application of hybrid steeltimber floors feasible for future construction projects. The assumptions related to the structural damping and the applicability of the joint stiffness parameters at higher frequencies must be kept in mind, and further investigations in this regard are still required. Moreover, analyses in a realistic setting utilizing additional structural components, e.g., floating floors, will lead to more detailed insights into the impact sound insulation. It is worth noting that while the resultant vibration amplitudes, as quantified by vibration dose values, mostly fall within the range of possible adverse comments, there remain possibilities to further investigate the damping behavior of largescale floor structures, which might support their broader adoption in construction practices. Laboratory tests on largescale floor structures might further be used to validate the impact force applied to the FE model using an ISO tapping machine.
Acknowledgments
The authors gratefully acknowledge the help of the chair of Timber Structures and Building Construction at Karlsruhe Institute of Technology for assistance with the measurements. Furthermore, the authors are thankful for the support of the manufacturers in providing the materials for the samples. The authors thank Semere Gebrekidan for the valuable comments and corrections that helped to improve the manuscript.
Funding
This study was financially supported by the German Federal Ministry of Food and Agriculture through the Fachagentur Nachwachsende Rohstoffe e.V. (FNR) (grant number 22009817) and the German Research Foundation (DFG) (grant number MA2395/152).
Conflict of interest
Author declared no conflict of interests.
Data availability statement
Data are available on request from the authors.
Author contribution statement
Bettina Chocholaty: Conceptualization, Methodology, Software, Data Curation, Validation, Formal Analysis, Visualization, Writing – original draft, Writing – review & editing. Nicolaas Bernardus Roozen: Conceptualization, Methodology, Writing – review, Supervision. KarlAlexander Hoppe: Methodology, Writing – review & editing. Marcus Maeder: Writing – review, Supervision. Steffen Marburg: Funding acquisition, Resources, Supervision, Writing – review.
References
 A. Chiniforush, M.M. Alamdari, U. Dackermann, H. Valipour, A. Akbarnezhad: Vibration behaviour of steeltimber composite floors, part (1): experimental & numerical investigation. Journal of Constructional Steel Research 161 (2019) 244–257. [CrossRef] [Google Scholar]
 A. Hassanieh, A.A. Chiniforush, H.R. Valipour, M.A. Bradford: Vibration behaviour of steeltimber composite floors, part (2): evaluation of humaninduced vibrations. Journal of Constructional Steel Research 158 (2019) 156–170. [CrossRef] [Google Scholar]
 M. Fujita, J. Sakai, H. Oda, M. Iwata: Building system for a composite steeltimber structure. Steel Construction 7 (2014) 183–187. [CrossRef] [Google Scholar]
 A.T. Balasbaneh, W. Sher, D. Yeoh, K. Koushfar: LCA & LCC analysis of hybrid glued laminated timber–concrete composite floor slab system. Journal of Building Engineering 49 (2022) 104005. [CrossRef] [Google Scholar]
 W. Zhu, H. Yang, W. Liu, B. Shi, Z. Ling, H. Tao: Experimental investigation on innovative connections for timber–concrete composite systems. Construction and Building Materials 207 (2019) 345–356. [CrossRef] [Google Scholar]
 A. Hassanieh, H. Valipour, M. Bradford: Experimental and numerical study of steeltimber composite (STC) beams. Journal of Constructional Steel Research 122 (2016) 367–378. [CrossRef] [Google Scholar]
 J. Negreira, A. Trollé, K. Jarnerö, L.G. Sjökvist, D. Bard: Psychovibratory evaluation of timber floors – towards the determination of design indicators of vibration acceptability and vibration annoyance. Journal of Sound and Vibration 340 (2015) 383–408. [CrossRef] [Google Scholar]
 J. Weckendorf, G. Hafeez, G. Doudak, I. Smith: Floor vibration serviceability problems in wood lightframe buildings. Journal of Performance of Constructed Facilities 28 (2014) A4014003. [CrossRef] [Google Scholar]
 A. Pavic: Results of IStructE 2015 survey of practitioners on vibration serviceability. In: Proceedings of the SECED 2019 Conference: Earthquake Risk and Engineering towards a Resilient Word, London, UK, 9–10 September, 2019. [Google Scholar]
 O.A. Hassan, F. Öberg, E. Gezelius: Crosslaminated timber flooring and concrete slab flooring: a comparative study of structural design, economic and environmental consequences. Journal of Building Engineering 26 (2019) 100881. [CrossRef] [Google Scholar]
 N. Perković, V. Rajčić, J. Barbalić: Analytical and numerical verification of vibration design in timber concrete composite floors, Forests 12 (2021) 707. [CrossRef] [Google Scholar]
 N.G. Vardaxis, D. Bard Hagberg, J. Dahlström: Evaluating laboratory measurements for sound insulation of crosslaminated timber (CLT) floors: configurations in lightweight buildings. Applied Sciences 12 (2022) 7642. [CrossRef] [Google Scholar]
 F.G. Branco, L. Godinho: On the use of lightweight mortars for the minimization of impact sound transmission. Construction and Building Materials 45 (2013) 184–191. [CrossRef] [Google Scholar]
 L. Godinho, R. Masgalos, A. Pereira, F. Branco: On the use of a smallsized acoustic chamber for the analysis of impact sound reduction by floor coverings. Noise Control Engineering Journal 58 (2010) 658–668. [CrossRef] [Google Scholar]
 Technical Committee ISO/TC 98: ISO 10137 – Bases for design of structures – Serviceability of buildings and walkways against vibrations, International Organization for Standardization, Geneva, Switzerland, 2007. [Google Scholar]
 B. Szabó, I. Babuška: Finite Element Analysis: Method, Verification and Validation. John Wiley & Sons, 2021. [CrossRef] [Google Scholar]
 J.E. Mottershead, M. Link, M.I. Friswell, C. Schedlinski: Model Updating. In: R. Allemang, P. Avitabile, Eds. Handbook of experimental structural dynamics. Springer, New York, NY, 2021, pp. 1–53. [Google Scholar]
 B. Shahriari, K. Swersky, Z. Wang, R.P. Adams, N. De Freitas: Taking the human out of the loop: A review of bayesian optimization. Proceedings of the IEEE 104 (2015) 148–175. [Google Scholar]
 S. Tao, A. Van Beek, D.W. Apley, W. Chen: Multimodel bayesian optimization for simulationbased design. Journal of Mechanical Design 143 (2021) 111701. [Google Scholar]
 D. Owolabi, C. Loss, J. Zhou: Vibration properties and serviceability performance of a modular crosslaminated timbersteel composite floor system. Journal of Structural Engineering 149 (2023) 04023171. [CrossRef] [Google Scholar]
 B. Chocholaty, N.B. Roozen, M. Maeder, S. Marburg: Vibroacoustic response of steel–timber composite elements. Engineering Structures 271 (2022) 114911. [CrossRef] [Google Scholar]
 F. Nogueira: Bayesian Optimization: Open source constrained global optimization tool for Python, 2014. Available at https://github.com/fmfn/BayesianOptimization. [Google Scholar]
 Technical Committee ISO/TC 43: DIN EN ISO 10140–5 – Acoustics – Laboratory measurement of sound insulation of building elements – Part 5: Requirements for test facilities and equipment, European Committee for Standardization, Brussels, Belgium, 2021. [Google Scholar]
 Pollmeier Massivholz GmbH & Co.KG: Leistungserklärung einer Platte in Baubuche.2018. Available at https://www.pollmeier.com/de/produkte/ueberbaubuche/baubucheplatte (accessed: 20220119). [Google Scholar]
 A. Andrej: Schneider: Bautabellen für Ingenieure, vol. 21. Bundesanzeiger Verlag, 2014. [Google Scholar]
 ANSYS, Inc.: Ansys engineering simulation software, 2019. Available at https://www.ansys.com/. [Google Scholar]
 P. Langer, M. Maeder, C. Guist, M. Krause, S. Marburg: More than six elements per wavelength: The practical use of structural finite element models and their accuracy in comparison with experimental results. Journal of Computational Acoustics 25 (2017) 1750025. [CrossRef] [Google Scholar]
 L. Cremer, M. Heckl: Abstrahlung von Körperschall. In: L. Cremer, M. Heckl, Körperschall. Springer, 1996, pp. 459–545. https://doi.org/10.1007/9783662081822. [CrossRef] [Google Scholar]
 D.J. Ewins: Modal testing – theory, practice and application. John Wiley & Sons, New York, 2009. [Google Scholar]
 S. Ereiz, I. Duvnjak, J.F. JiménezAlonso: Review of finite element model updating methods for structural applications. Structures 41 (2022) 684–723. [CrossRef] [Google Scholar]
 R.J. Allemang: The modal assurance criterion – twenty years of use and abuse. Sound and vibration 37 (2003) 14–23. [Google Scholar]
 O. Cuate, O. Schütze: Pareto explorer for finding the knee for many objective optimization problems. Mathematics 8 (2020) 1651. [CrossRef] [Google Scholar]
 K. Christodoulou, E. Ntotsios, C. Papadimitriou, P. Panetsos: Structural model updating and prediction variability using pareto optimal models. Computer Methods in Applied Mechanics and Engineering 198 (2008) 138–149. [CrossRef] [Google Scholar]
 E. Brochu, V.M. Cora, N. De Freitas: A tutorial on bayesian optimization of expensive cost functions, with application to active user modeling and hierarchical reinforcement learning. 2010. Preprint available at https://arxiv.org/abs/1012.2599. [Google Scholar]
 C.Q. Howard, B.S. Cazzolato: Acoustic Analyses using MATLAB^{®} and ANSYS^{®}. CRC Press, 2014. [Google Scholar]
 Technical Committee CEN/TC 250: DIN EN 1995 1–1 – Eurocode 5: Design of timber structures, Part 1–1: General – Common rules and rules for buildings, European Committee for Standardization, Brussels, Belgium, 2010. [Google Scholar]
 C. Geweth, S. Baydoun, F. Saati, K. Sepahvand, S. Marburg: Effect of boundary conditions in the experimental determination of structural damping. Mechanical Systems and Signal Processing 146 (2021) 107052. [CrossRef] [Google Scholar]
 P. Wang, C. Van Hoorickx, G. Lombaert, E. Reynders: Numerical prediction and experimental validation of impact sound radiation by timber joist floors. Applied Acoustics 162 (2020) 107182. [CrossRef] [Google Scholar]
 N. CheraghiShirazi, K. Crews, S. Malek: Review of vibration assessment methods for steeltimber composite floors. Buildings 12 (2022) 2061. [CrossRef] [Google Scholar]
 T. Murray, D. Allen, E. Ungar, D. Davis: Vibrations of steelframed structural systems due to human activity: AISC design guide 11. American Institute of Steel Construction, USA, 2016. [Google Scholar]
 M.R. Willford, P. Young: A design guide for footfall induced vibration of structures. Concrete Society for the Concrete Centre London, London, UK, 2006. [Google Scholar]
 Committee GME/21: BS 6472–1 – Guide to evaluation of human exposure to vibration in buildings – part 1: vibration sources other than blasting, 2008. https://doi.org/10.3403/19971044U. [Google Scholar]
 J. Brunskog, P. Hammer: The interaction between the ISO tapping machine and lightweight floors. Acta Acustica united with Acustica 89 (2003) 296–308. [Google Scholar]
 J. Lietzén, J. Miettinen, M. Kylliäinen, S. Pajunen: Impact force excitation generated by an ISO tapping machine on wooden floors. Applied Acoustics 175 (2021) 107821. [CrossRef] [Google Scholar]
 D. Fritze, S. Marburg, H.J. Hardtke: Estimation of radiated sound power: a case study on common approximation methods. Acta Acustica united with Acustica 95 (2009) 833–842. [CrossRef] [Google Scholar]
 F. Fahy: Sound and structural vibration: radiation, transmission and response, Academic Press, London and Orlando, FL, 1985. [Google Scholar]
 Technical Committee CEN/TC 211: DIN EN ISO 3741 – Acoustics – determination of sound power levels and sound energy levels of noise sources using sound pressure – Precision methods for reverberation test rooms, 2011. https://doi.org/10.31030/1642422. [Google Scholar]
 Subcommittee NA 001020303 UA: DIN 18041 – Acoustic quality in rooms – Specifications and instructions for the room acoustic design, 2016. https://doi.org/10.31030/2395845. [Google Scholar]
 Technical Committee CEN/TC 126: DIN EN ISO 123542 – Building acoustics – estimation of acoustic performance of buildings from the performance of elements – part 2: impact sound insulation between rooms, 2017. https://doi.org/10.31030/2625099. [Google Scholar]
Cite this article as: Chocholaty B. Roozen NB. Hoppe KA. Maeder M. & Marburg S. 2024. Vibration and impact sound properties of hybrid steeltimber floor structures. Acta Acustica, 8, 11.
All Tables
Numerically f_{fem} and experimentally f_{exp} determined natural frequencies [Hz], relative frequency error and MAC before and after the automatic model calibration: m ∈ [1, 6] for S1* and S2*.
Initial p_{0} and calibrated parameters for the closed S1* and open S2* specimen: u and l denote the upper and lower joint, respectively. All values in the table refer to the properties of nails except for the entries marked by †, which refers to screw parameters.
Fundamental frequencies f_{1,j} [Hz], modes for the VDV evaluation f_{m,j} [Hz], and modal masses of the samples j ∈ [S1,f; S1,c; S2,f; S2,c].
VDV during the day VDV_{d,j} [ms^{−1.75}] and VDV during the night VDV_{n,j} [ms^{−1.75}] of the samples j ∈ [S1, f; S1, c; S2, f; S2, c].
All Figures
Figure 1 Overall workflow of the study: calibration of unknown model parameters p is performed utilizing numerically determined natural frequencies as well as mode shapes and results from experimental modal analyses on the smallscale test specimens, S1* and S2*. The updated parameters are used in a subsequent Finite Element simulation to analyze the largescale floor specimens, S1 and S2, regarding vibration dose values and impact sound pressure levels. The impact sound pressure level simulations use a numerically implemented tapping machine load. 

In the text 
Figure 2 Snippets of closed (left, S1) and open (right, S2) hybrid LVL and trapezoidal steel elements. The lower timber plate is built of a continuous lower wooden panel for closed crosssections and made of several separate panels for open configurations. The respective geometrical parameters are given in Table 1. 

In the text 
Figure 3 Fastener model: an MPC184 general joint element (represented by a black line) using a stiffness matrix between nodes related to fasteners on the substructures (grids). Coincident nodes at the interface of steel and timber are coupled related to the zdisplacement (green lines) within an area surrounding the fastener (green circle). 

In the text 
Figure 4 Flowchart of the model calibration: An FE model of the respective measured sample S1* or S2* is calibrated by iteratively adapting the unknown parameters p and computing an error e using the numerical and experimental results. 

In the text 
Figure 5 Scan setup using heavyduty slings for suspension with Laser Doppler Vibrometer (LDV) (a) and shaker (b). 

In the text 
Figure 6 Mode shapes for the sample S1*: experimentally determined (left) and numerically simulated (right) after model calibration. 

In the text 
Figure 7 Mode shapes for the sample S2*: experimentally determined (left) and numerically simulated (right) after model calibration. 

In the text 
Figure 8 Pareto optimal front of the multiobjective optimization for the closed S1* (a) and open S2* (b) specimen. 

In the text 
Figure 9 FE model of closed sample S1 (a) and open sample S2 (b): bottom view and side views. Steel is colored purple, and timber parts are cyan. The element thickness is visually scaled to represent the thickness of the respective component. 

In the text 
Figure 10 Example plot of load F_{n} against frequency. 

In the text 
Figure 11 Normalized impact sound pressure of samples S1 (a) and S2 (b) for the cases “f” and “c” for thirdoctave bands. 

In the text 
Figure 12 Comparison of the sound pressure levels of the samples S1 (closed) and S2 (open) for the cases “f” (a) and “c” (b) for thirdoctave bands. 

In the text 
Figure 13 Comparison of the upper and lower bounds of the sound pressure levels computed for the hybrid floor samples S1 and S2 with the values for a 100 mm thick concrete floor with a 20 mm thick bonded screed and a 180 mm thick concrete floor with a 50 mm thick bonded screed provided in [40]. Results are given in thirdoctave bands for the upper and lower bounds and in octaves for the concrete floors. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.