Analytical and Numerical Veriﬁcation of Vibration Design in Timber Concrete Composite Floors

: The TCC concept has been studied and developed over the past decades. The variety of solutions shows the meaningfulness and functionality of this system, as well as the continuous work of scientists over time. To beneﬁt from these advantages, the composite needs to provide sufﬁcient stiffness to meet the serviceability criteria and load capacity to resist loading at every stage of the building life. An example of connector types and load slip curves according to EN 1995 is given. This paper discusses possible limitations related to residential areas, and additionally, the possible solutions that EN 1995 does not discuss in the case of resonant response (f 1 < 8 Hz). The theoretical studies were accompanied by numerical analyses considering certain simpliﬁcations suitable for practical use.


Introduction
The composite system based on concrete and timber represents a very successful engineering solution in the process of optimization of construction. The timber-concrete composites (TCC) date back to the beginning of the 20th century. The first scientific approach was applied to the system of nails and steel braces as connectors between the concrete slab and timber beams, patented by Müller in 1922 [1]. However, it has been widely used only in the last 30 years [2], when it proved to be a good solution for the renovation of existing timber floors due to today's code requirements of load-bearing capacity, deformability, vibration, fire and earthquake resistance, and sound insulation [3]. Moreover, due to the conditions of overcoming larger spans without the use of walls (pillar-slab or frame-slab system), as well as the requirement for the floor structure to be a rigid diaphragm, it has become a must-have solution in the design of tall timber buildings [4][5][6][7][8][9]. The advantages of TCC construction in the building industry are achieved through the mass of concrete, compared with timber floors, as well as on the positive characteristics of the wood (high tensile strength) and concrete (high compressive strength). The disadvantages, on the other side, could be found in greater dead loads, which results in higher manufacturing and assembly demands. Analysis of the possible application of TCC is given in [10,11].
Economic considerations of the cost-effectiveness of TCC given by [12] tried to estimate market opportunities and obstacles when using TCC. Through the questioning of experts, it provided a prognosis and assessment of the market potentials for TCC in Germany based on a realistic scenario. The scenario presented in [12] results in a potential floor area for TCC for multistorey residential construction, single and two-family house construction, and nonresidential building. Altogether, a potential of approx. 1.6 million m 2 is the result (the market figures are based on Germany in 2014). The establishment of TCC for the wood-based sector makes significantly higher demands than what would normally occur with technical innovations where TCC has a considerable impact on the identity of building with wood. Such an approach to the market requires conditions that will bring the benefits of TCC closer to the user, with a sense of comfort being one of the most important factors. Therefore, the vibration conditions, or the way of the connection (between wood and concrete) principles which defined them, become the main accent in the design of floor structures.
The usage of optimization in product development is growing as both computational capabilities and newly developed numerical algorithms are increasing. At the same time, all design problems have more than one conflicting objective in design optimization, and these objectives need to be integrated to yield one final design. Thus, the design problem can be formulated as a multi-objective optimization problem [13]. In the serviceability design of the TCC floor, the evaluation index plays a critical role, which makes it essential to check the vibration behavior and occupant comfort. In this paper, two objectives of the cost and demanded function and comfort of the use of the TCC floor system have been chosen as two capital objectives. To achieve a well-engineered TCC floor design, a framework based on optimal design variables include the optimized thickness of the concrete and the optimized smeared stiffness of connectors, which are considered regarding the TCC floor structure vibrations.
FEA parametric analysis was performed to obtain the fundamental frequency of the composite floor with different types of connectors, different concrete slab thicknesses, and two types of span (short-and large-span system) for the usual type of boundary conditions (simply supported). Based on the numerical analysis, three existing analytical models were applied to estimate the correlation of fundamental frequency results. Conclusions are given on the optimal design of the elements.

TCC Floor System Subjected to Footfall Induced Vibration
Analysis of TCC floors in the area of serviceability, besides vibration, is developing progressively due to the increasing application of larger spans [14]. As vibrations performance depends on the floor system mass, the longitudinal and transverse bending stiffness, the damping of the floor system, and the shear stiffness of the joints between timber and concrete, changing each of these parameters will affect the dynamic response. The influence of the cracks developing in the concrete slab and the deflections of the supporting beams are also considered as one of the influencing parameters. At present, plenty of studies focus on such dynamic characteristics as natural frequencies, damping ratios, and mode shapes of TCC floor systems; in among the others mentioned here, there are more studies [15][16][17][18][19][20][21][22][23][24][25][26][27]. If a fixed-length system is considered, it can be concluded that the thickness of the concrete slab contributes to stiffness as well as mass, but the increased mass can outweigh the improvement to stiffness, thus causing the fundamental frequency to decrease. Furthermore, short-span systems (up to 6 m) determine acceptable dynamic responses to footfall-induced vibration, since such floors usually have sufficient transversal bending stiffness with natural eigenfrequency above 8 Hz. However, with the large-span length (above 6 m), the floors' natural eigenfrequency usually does not meet the requirements above 8 Hz [28].
According to [29] for floor performance levels I-VI, no further investigations are necessary if the requirements in respect to fundamental frequency, acceleration or velocity, and stiffness from Table 1 are satisfied. The root mean square acceleration or velocity responses are compared to the vibration perception base curve in [30]. The acceleration criterion is expressed as a multiple of the base curve value. This multiple is termed as the response factor R (equal to mean square velocity response v rms /0.0001) given in Table 1. For resonant vibration design situations, (when f 1 < 8 Hz), the minimum fundamental frequency, acceleration, and stiffness criteria of Table 1 should be fulfilled. For transient vibration design situations, (when f 1 ≥ 8 Hz), the velocity and stiffness criteria of Table 1 should be fulfilled. A structure with timber-concrete composite parts should be designed and constructed such that all relevant serviceability limit states are satisfied according to the principles of [31], where criteria are stated in Subclause 3.4 of [31]. The dynamic properties of floor beams, then, should satisfy the criteria of Subclause 1.4.4 of [31]. Lightweight (timber/wood) structures with a span length below approximately 6 to 8 m, and heavy (concrete) floors "with short span" are related to the high-frequency, while the low-frequency floor includes heavyweight floors with a longer span [32]. The vibration level, according to [33], should be estimated by measurement or by calculation, taking into consideration the expected stiffness of the member, component, or structure, and the modal damping ratio. Residential floor structures have often been treated as either low-or high-frequency floors, where suggestions for limits of vibration are given in Subclause 7.3 of [34]. People most often feel discomfort at floor frequencies between 4 and 8 Hz, which corresponds to the natural frequencies of some organs of the human body [35,36]. Therefore, the authors of [34] recognized that, in case of a fundamental frequency that is less than 8 Hz, a special investigation should be made, and, if a fundamental frequency is greater than 8 Hz, the following requirements should be satisfied: where w is the maximum instantaneous vertical deflection caused by a vertical concentrated static force F applied at any point on the floor, taking account of load distribution, v is the unit impulse velocity response, i.e., the maximum initial value of the vertical floor vibration velocity (in m/s) caused by an ideal unit impulse (1 Ns) applied at the point of the floor giving the maximum response (according to the [34], components above 40 Hz may be disregarded), and ς is the modal damping ratio. These calculations should be made under the assumption that the floor is unloaded, i.e., only the mass corresponding to the self-weight of the floor and other permanent actions. For a rectangular floor with an overall length dimension of l and width b, simply supported along all four edges and with timber beams having a span e, the natural undamped eigenfrequency f 1 may approximately be calculated as: where m is the mass per unit area in kg/m 2 , l is the floor span in m, and (El) l is the equivalent plate bending stiffness of the floor about an axis perpendicular to the beam direction in Nm 2 /m. The concrete slab in timber-concrete composite systems is relatively narrow compared to pure concrete slabs. Therefore, the reinforcement is often installed near the centroid of the cross-section. The cracking of concrete will lead to a significant drop in the bending stiffness. The effective width is caused by the distribution of the normal force in the concrete cross-section as a shell and the distribution of the bending moment as a plate. Due to the decrease of the bending stiffness by cracking, the load distribution in the concrete slabs tends to be only in the shell model. However, the effective width given in [37] is comparable to the shell mode, whereas in [38], the plate mode is also considered [33,39,40]. For that reason, it is appropriate to use the effective width according to [37] instead of [38].
Considering the need to simplify the calculation, the composite floor can be considered as a parallel arrangement of composite beams. The equation above can be represented as: where b c is the spacing between the timber beams [32]. The natural frequencies of the composite beams and the spacing of the composite beams are the same as those of the composite floor [41]. For the same type of floor, the value v may, as an approximation, be taken as: v ≤ 4(0.4 + 0.6n 40 )/(mbl + 200), in m/Ns 2 , where v is the unit impulse velocity response in m/(Ns 2 ), m is the mass in kg/m 2 , b is the floor width in m, l is the floor span in m, and n 40 is the number of first-order modes with natural frequencies up to 40 Hz which may be calculated from: where (EI) b is the equivalent plate bending stiffness, in Nm 2 /m, of the floor about an axis parallel to the beams, where (E1) b < (EI) l . Series of studies and research, of which the results are given in [42][43][44][45][46][47][48][49], were a background for choosing values for the appropriate modal damping ratio. According to [33], for floors, unless other values are proven to be more appropriate, a modal damping ratio of ς = 0.025 (i.e., 2.5%) should be used for timber-concrete composite slabs alone and ς = 0.035 slabs with a floating screed. The instantaneous elastic bending stiffness of the composite structure should be used in vibration analysis. Alternatively, damping ratio values can be calculated based on actual loads (i.e., precomposite dead load, postcomposite dead load, and sustained live load) [50]. More accurate damping values can alternatively be obtained by samples testing, applying [51]. At present, some other proposed models for the natural frequencies occurred. The authors of [52] proposed Equation: where y w is the weighted average of the static deflection of the beam due to the selfweight of the floor or the dead weight of the beam-slab system, and g is the gravitational acceleration. A similar method is proposed in [53], where the fundamental frequency of simply supported beams or floor systems under uniformly distributed additional loads, as shown in the equation: where E is Young's modulus of the beam, I is the moment of inertia of the beam, w is the uniform load per unit length, and l is the span. For the bending stiffness of composite beams with pinned boundary conditions, the γ-method determined in [33] should be used. The effective stiffness of the partially composite beam is calculated as follows: EI ef = E c (I c + γ c A c a c 2 ) + E t (I t + γ t A t a t 2 ), in Nmm 2 , where subscripts c and t refer to concrete and timber elements, respectively, E is Young's modulus of the material, A and I are the area and the second moment of area of the element cross-section, a is the distance from the centroid of the element to the neutral axis of the composite section, and γ is the shear connection reduction factor according to: γ c = 1/(1 + (π 2 E c A c s)/(K ser l 2 )), where s is the spacing of the connectors assumed as smeared along the span, l is the span, and K ser is the slip modulus of the connector. The authors of [54] proposed a method to calculate the effective stiffness of composite beam with an arbitrary boundary as: where µ is the buckling length factor, r is the distance between the centroid of timber section to the centroid of the concrete section, α is the nondimensional composite action (shear connector) parameter, EI 0 the bending stiffness of the noncomposite section (K→0), and EI ∞ is the bending stiffness of the fully composite section. EA 0 and EA p are the sum and product of axial stiffness of the sub-elements, respectively. Some design methods, such as the γ-method, assume that the connectors may be smeared along the beam axis. In this case, no discrete connector is considered, but the stiffness per unit length is taken into account. For short distances, this simplification works quite well. However, with increasing distance, the discrete connector leads to local stresses. This conclusion is reported in [55,56]. The authors of [55] propose a maximum distance of 5%, whereas [56] proposes a maximum distance of 3%. For notches as the connection between timber and concrete, [57] proposes an effective distance in order to enable smearing of these connectors. Providing that the effective spacing of the connections is less than or equal to 5% of the distance between the points of contra flexure, even distribution of the connection stiffness along the beam axis (smearing) may be used. If the spacing of the connections is greater than 5% of the distance between the points of contra flexure, then, when smearing, the connections should be distributed in proportion to the shear force. Comparing the sinusoidal course and the constant course of the normal force, the different deformation can be evaluated. In the hyperstatic composite cross-section, the stiffness influences the distribution of the forces. The stiffness of the cross-section connected with discrete connectors should be reduced by the factor 2/π to consider the same stiffness as in the smeared case. For the reason of simplification, this factor is rounded to 0.7, since the case of only one single connector at each of the beams is hardly used. However, up to now, no other values are given, since the variability of possible positions of the connectors is quite large. So, in the absence of a more precise model, only 70% of the axial stiffness of the attached cross-section should be considered for the calculation of stresses and deformation. For the calculation of the shear force in the connection, 100% of the axial stiffness of the cross-section should be considered [33]. If the spacing and/or stiffness of the connections are varied in proportion to the shear force, the effective spacing may be determined as: where s ef is the effective spacing of the connections, s min is the minimum spacing of the connections, K ref is the reference stiffness of the connection used with the corresponding value of s ef , K max is the maximum stiffness of the connection, K min is the minimum stiffness of the connection, and s max is the maximum spacing of the connections or the maximum distance between the connection and the point of zero shear force. The slip modulus for serviceability limit states that K ser of connections made with dowel-type fasteners inserted perpendicular to the shear plane should be determined using: Forests 2021, 12, 707 6 of 20 for dowels, bolts, screws, and nails (with predrilling), and nails (without predrilling): where K ser is the slip modulus for serviceability limit states, ρ m is the mean value of the timber member density in kg/m 3 , and d is the fastener diameter in mm. Research studies [58][59][60][61] show that this can be a good estimation when more detailed data is not available. For regular interlayers, with stiffness perpendicular to the shear plane similar to that of timber and with thickness up to 30 mm, the slip modulus of connections with dowel-type fasteners may be taken as that for a similar configuration without an interlayer, with a reduction of 30% [11]. In other cases where there is an intermediate nonstructural layer between the timber and the concrete, the slip modulus should be determined by tests or special analysis. For connections made with steel rebar glued into timber perpendicular to the shear plane using epoxy resin, the slip modulus for serviceability limit states K ser should be determined according to: where E tim is the mean modulus of elasticity of timber parallel to the grain in N/mm 2 , and d is the nominal diameter of the rebar in mm [33,58,60,62]. In many cases, when the manufacturer provides new, modern types of connectors, such as [63], slip modulus should be determined based on verified tests. For notched connections, the slip modulus for serviceability limit states that K ser should be determined according to: where h n is the depth of the notch. Linear interpolation may be used for h n between 20 mm and 30 mm [33,64]. By using the effective bending stiffness expressed above, three different calculation models have been combined as a function of the two design variables: concrete thickness hc and overspread connector stiffness K ser /s, where s is spacing between connectors, and K ser /s is stiffness per unit length of the element expressed in (N/mm)/mm. By applying the deflection criteria, the minimum bending stiffness fulfilling the eigenfrequency criteria can be predicted.

FEM Analysis of TCC Floor Slab
Numerical modeling is a powerful tool aimed at expanding knowledge and saving time and, ultimately, finance. Efficient data input and intuitive handling facilitate the modeling of simple and large structures. The numerical analyses aimed to extend the knowledge of the behavior of the tested system. Furthermore, the numerical simulations served to confirm and complement the experimental results. In this paper, the RFEM program is used, a powerful 3D FEA program helping structural engineers meet requirements in modern civil engineering. One such case is certainly composite systems, specifically, timber-concrete systems. An additional aggravating circumstance is defining and calculating a semi-rigid connection between the different elements and materials. There are several options to calculate a semi-rigid composite beam or floor. The main difference is in the modeling method itself. Some methods ensure simple modeling (such as the gamma method). However, there are other more complex methods (shear analogy). Another option for modeling a composite system is shown in this paper. Since the definition and analysis of connecting elements are time-consuming, it is recommended to connect the surfaces of the elements to the other surfaces directly. Although there are several options in this software, the coupling member surface with the line release option will be shown below.
The structural system and the cross-section dimensions are shown in Figure 1. The model consists of timber elements that are defined as a member, with a given eccentricity. At the contact lines of wooden beams and concrete slabs, the line hinge is created. Thus, it is possible to release one of the components and to control the relation to each other by a line release type. The flexibility is considered by defining the spring constant. This method has proven to be very accurate in determining the internal forces and the effect of shrinkage, but the disadvantage is that the shear flow cannot be read directly, and the solution is shown below. This paper looks into the HBV-shear connector for a TCC system. Recently, the HBV-shear connector has gained popularity because it is the only connector to create a wide area connection that corresponds to slip modulus, which describes the efficiency of a connector, which is high compared to other connectors [65]. It is important to note that this only applies to serviceability classes 1 and 2, and concrete minimum class C20/25. It is recommended that the thickness of the concrete slab is 70 mm minimum. The effect of cracks on the bearing was taken into account, as well as the effect of creep and shrinkage. Finally, vibration analysis was performed.

Compliance of the Connection-Kser
An additional line is generated and it is associated with the concrete surface, while the lower line is associated with the timber member. The spring is defined accordingly to the degree of freedom ( Figure 2). The model consists of timber elements that are defined as a member, with a given eccentricity. At the contact lines of wooden beams and concrete slabs, the line hinge is created. Thus, it is possible to release one of the components and to control the relation to each other by a line release type. The flexibility is considered by defining the spring constant. This method has proven to be very accurate in determining the internal forces and the effect of shrinkage, but the disadvantage is that the shear flow cannot be read directly, and the solution is shown below. This paper looks into the HBV-shear connector for a TCC system. Recently, the HBV-shear connector has gained popularity because it is the only connector to create a wide area connection that corresponds to slip modulus, which describes the efficiency of a connector, which is high compared to other connectors [65]. It is important to note that this only applies to serviceability classes 1 and 2, and concrete minimum class C20/25. It is recommended that the thickness of the concrete slab is 70 mm minimum. The effect of cracks on the bearing was taken into account, as well as the effect of creep and shrinkage. Finally, vibration analysis was performed.

Compliance of the Connection-K ser
An additional line is generated and it is associated with the concrete surface, while the lower line is associated with the timber member. The spring is defined accordingly to the degree of freedom ( Figure 2).

Results
The process of designing and checking vibration and other construction requirements is shown in Figure 5. As it is a residential space, it is important to determine whether the condition of the allowed frequencies of the TCC floor slab are fulfilled. This part was done in a special module DYNAMPro (Figure 6), and the result is given in the tabular form (Table 4), where it is seen that the dominant first mode reaches 9.967 Hz, which is higher than the required 8 Hz so that the requirement is fulfilled. As it is a residential space, it is important to determine whether the condition of the allowed frequencies of the TCC floor slab are fulfilled. This part was done in a special module DYNAMPro (Figure 6), and the result is given in the tabular form (Table 4), where it is seen that the dominant first mode reaches 9.967 Hz, which is higher than the required 8 Hz so that the requirement is fulfilled. For practical use, the simplified natural frequency calculation rule is very useful. In this analysis, the deflection under permanent and quasi-permanent action at the ideal onespan beam may not exceed a limit value (6 mm according to DIN 1052). If the relationship shown in Figure 7 between natural frequency and deflection is considered for a hinged one-span beam that is subjected to a constant uniform load, then the 6 mm of deflection corresponds to a minimum natural frequency of about 7.2 Hz. Taking into account the fact that, in most National Annexes of EC 5, a minimum natural frequency of 8Hz is to be considered, a maximum deflection of about 5 mm is allowable. According to Figure 8 and the simplified method, frequency is higher than 8 Hz, so the postulate is fulfilled.    For practical use, the simplified natural frequency calculation rule is very useful. In this analysis, the deflection under permanent and quasi-permanent action at the ideal onespan beam may not exceed a limit value (6 mm according to DIN 1052). If the relationship shown in Figure 7 between natural frequency and deflection is considered for a hinged one-span beam that is subjected to a constant uniform load, then the 6 mm of deflection corresponds to a minimum natural frequency of about 7.2 Hz. Taking into account the fact that, in most National Annexes of EC 5, a minimum natural frequency of 8Hz is to be considered, a maximum deflection of about 5 mm is allowable. According to Figure 8 and the simplified method, frequency is higher than 8 Hz, so the postulate is fulfilled.  For practical use, the simplified natural frequency calcu this analysis, the deflection under permanent and quasi-perma span beam may not exceed a limit value (6 mm according to D shown in Figure 7 between natural frequency and deflection one-span beam that is subjected to a constant uniform load, corresponds to a minimum natural frequency of about 7.2 Hz. that, in most National Annexes of EC 5, a minimum natura considered, a maximum deflection of about 5 mm is allowabl the simplified method, frequency is higher than 8 Hz, so the p fe

Stiffness Design
According to [19], the stiffness criterion is very important to check, because the frequency criterion alone is not sufficient. In this paper, the stiffness under a static load of 2 kN is checked. The deflection due to the load needs to be less than the deflection limit value. If the vibrations are noticeable only when you concentrate on them and are not annoying, then w limit = 0.5 mm. This is usually the case for floors in residential or office buildings. On the other hand, if the vibrations are quite noticeable, but not annoying, then w limit = 1.0 mm (e.g., family house). The analytical approach of this design is to create a replacement simple beam system, with an appropriate span corresponding to the largest range of the TCC floor. According to [19,67,68], the deflection may be calculated as follows: w (2kN) = (2·l)/(48·EI l ·b w(2kN) ) ≤ w limit (31) b w(2kN) = min{(b ef @ b (width of the floor)) According to [19], the stiffness criterion is very important to check, because the frequency criterion alone is not sufficient. In this paper, the stiffness under a static load of 2 kN is checked. The deflection due to the load needs to be less than the deflection limit value. If the vibrations are noticeable only when you concentrate on them and are not annoying, then wlimit = 0.5 mm. This is usually the case for floors in residential or office buildings. On the other hand, if the vibrations are quite noticeable, but not annoying, then wlimit = 1.0 mm (e.g., family house). The analytical approach of this design is to create a replacement simple beam system, with an appropriate span corresponding to the largest range of the TCC floor. According to [19,67,68], the deflection may be calculated as follows: Although in [34], for stiffness design, a static load of 1 kN was used and a substitutional static system, based on research and work [10]. In this paper, a load of 2 kN and the real TCC system were examined.
It is shown (Figure 9) that the deflection is obtained from 0.3 mm, which satisfies the requirements following [19].
Although in [34], for stiffness design, a static load of 1 kN was used and a substitutional static system, based on research and work [10]. In this paper, a load of 2 kN and the real TCC system were examined.
It is shown (Figure 9) that the deflection is obtained from 0.3 mm, which satisfies the requirements following [19]. 0.03 < 0.5 mm (34)

Connector Design
Lastly, special attention needs to be paid to the forces that occur in the connector, which are transferred from the member to the surface and conversely. To determine the shear flow, sections were made exactly at the intersection of the wooden elements and the concrete slab, ie., at the joint location. bef = l/1.1•∜((EIb/(EIl) = b/(1.1•α)) (33) Although in [34], for stiffness design, a static load of 1 kN was used and a substitutional static system, based on research and work [10]. In this paper, a load of 2 kN and the real TCC system were examined.
It is shown (Figure 9) that the deflection is obtained from 0.3 mm, which satisfies the requirements following [19].

Connector Design
Lastly, special attention needs to be paid to the forces that occur in the connector, which are transferred from the member to the surface and conversely. To determine the shear flow, sections were made exactly at the intersection of the wooden elements and the concrete slab, ie., at the joint location.
Observing the middle beam and superimposing the effect on both sides, in Figure 10, the total force is obtained. The average internal force amounts to 5.01 kN/m', while the total force is 9.21 kN, which represents the average value multiplied by the length ( Figure  11). The total force on the connector is 9.21 kN multiplied by 2 (both sides). After applying the same principle in the area of every connector, it can be concluded that the maximum force that occurs is less than the allowable previously calculated load capacity (21). Observing the middle beam and superimposing the effect on both sides, in Figure 10, the total force is obtained. The average internal force amounts to 5.01 kN/m', while the total force is 9.21 kN, which represents the average value multiplied by the length (Figure 11). The total force on the connector is 9.21 kN multiplied by 2 (both sides). After applying the same principle in the area of every connector, it can be concluded that the maximum force that occurs is less than the allowable previously calculated load capacity (21).

Concrete Check
Stress analysis in concrete was performed. The tensile stress in the lower zone of the concrete slab was primarily checked according to EN 1992-1-1, NDP (7.2). The limitation of concrete compressive stress is σc, max = 11.25 N/mm 2 , and as Figure 12 shows, maximal   Figure 11. Sum and average internal force on the connector.

Concrete Check
Stress analysis in concrete was performed. The tensile stress in the lower zone of the concrete slab was primarily checked according to EN 1992-1-1, NDP (7.2). The limitation of concrete compressive stress is σ c, max = 11.25 N/mm 2 , and as Figure 12 shows, maximal stress in the TCC floor slab is 2.4 N/mm 2 .

Acceleration Design
The additional examination of acceleration is successful only in the case of timberconcrete composite systems, or other heavy systems with wide spans. If the frequency is smaller than 8 Hz, an acceleration check is required. Empirical investigations in the 80 s offered the results that the vibration with a frequency lower than 8 Hz is dependent on the acceleration, and higher than 8 Hz is dependent on the velocity. The velocity check is not often considerable, because a typical timber floor is heavy enough (acoustic reasons), therefore, the velocity check is not relevant. The surface-to-surface model is shown in Figure 13.

Acceleration Design
The additional examination of acceleration is successful only in the case of timberconcrete composite systems, or other heavy systems with wide spans. If the frequency is smaller than 8 Hz, an acceleration check is required. Empirical investigations in the 80 s offered the results that the vibration with a frequency lower than 8 Hz is dependent on the acceleration, and higher than 8 Hz is dependent on the velocity. The velocity check is not often considerable, because a typical timber floor is heavy enough (acoustic reasons), therefore, the velocity check is not relevant. The surface-to-surface model is shown in Figure 13. offered the results that the vibration with a frequency lower than 8 Hz is dependent on the acceleration, and higher than 8 Hz is dependent on the velocity. The velocity check is not often considerable, because a typical timber floor is heavy enough (acoustic reasons), therefore, the velocity check is not relevant. The surface-to-surface model is shown in Figure 13. After the dynamic calculation, the results show that the natural frequency, in this case, is 6.228 Hz, which is less than 8 Hz (see Table 5). According to the analytical equation, acceleration is calculated as follows: a = F dyn /(M·2D) = (0.4·1·700 N)/(m·0.5·l·0.5·b·2·D) = 56/(m·l·b·D) (35) where: After the analysis, the results shown in Figure 14 indicates that the maximum acceleration at the relevant point occurs before the 10th second, and it is lower than the limit value of 0.05 m/s 2 . Besides, the RMS value is displayed. For the floor to be deemed acceptable this value must be less than a limit based on the floor's fundamental natural frequency. However, RMS acceleration can only be used for more resonant excitation.

Constructive Design Requirements
Since the vibration criterion is closely related to mass and stiffness, it is recommended to use heavy screeds and generally heavier layers, not less than 60 kg/m 2 . In the illustrative example, Table 6 shows that the sum of the loads in the Z direction is equal to 119.27 kN. If this is divided by the total area of 7.0 × 4.8 m, it obtains a satisfactory weight of 354 kg/m 2 , thus fulfilling this condition.
• b-the width of the floor (<1.5 l) • D-damping (following [70,71]) After the analysis, the results shown in Figure 14 indicates that the maximum acceleration at the relevant point occurs before the 10th second, and it is lower than the limit value of 0.05 m/s 2 . Besides, the RMS value is displayed. For the floor to be deemed acceptable this value must be less than a limit based on the floor's fundamental natural frequency. However, RMS acceleration can only be used for more resonant excitation.

Constructive Design Requirements
Since the vibration criterion is closely related to mass and stiffness, it is recommended to use heavy screeds and generally heavier layers, not less than 60 kg/m 2 . In the illustrative example, table 6 shows that the sum of the loads in the Z direction is equal to 119.27 kN.

Comparison Study
When reviewing the existing literature, the difference in detailing steps for the calculation according to NDS 2018 and Eurocode 5 is noticeable. According to the results shown in [72], the main distinction between the two codes is where the deformation/creep factor is applied in the design process. Following the US calculation method, this factor is added at the end of the calculation and, consequently, affects only the deflections. On the other hand, EC5 applies the creep factor right at the beginning of the calculation to the modulus of elasticity. The long-term deflection shows the difference between the European and US method of calculations. The European method for long-term total load deflection is greater and more conservative than the US method [72]; therein lies the cause of the minor scattering of FEM results of this case study compared to certain results available in the literature.
Furthermore, the parametric analysis (software SAP 2000) shown in [73] was focused on the structural components of the TCC floor and the geometry of the floor. It can be concluded that timber joists were found to be the most influential on the natural frequency of the floor, and increasing the depth of the timber joists of TCC floors is the most effective means of increasing the floor's frequency.
Finally, the influence of the connector types on the behavior of the TCC system is highlighted. The reason for using the HBV connector for this case study is that HBV connectors have been shown to have about 10% better efficiency compared to mechanical fasteners (STS shear connectors) [74]. Finite element analysis and parametric study in [74] show that the thickness of the steel connector is closely related to the shear stiffness and strength of the TCC system. Another important conclusion useful in practice is that the gap between the concrete and timber, up to a certain width, does not affect the slip modulus of the connecting system, which leaves the space for the insulation material that can improve the vibration and sound performance of the TCC floor.

Discussion
Special attention should be paid to areas where there are larger slab spans or where there are higher loads. In this case, frequencies of less than 8 Hz occur, and then the acceleration requirement should be checked. However, EN 1995-1-1 does not prescribe special measures and conditions to be met. As a recommendation and plan for further research, it is proposed to see the influence of walls and boundary conditions of the TCC floor, and especially nonload-bearing walls influence of vibrations. Furthermore, the bending stiffness of a floating floor should be taken into account. There are a lot of technical topics and issues which need to be addressed and taken into account in the design process of TCC, which makes TCC a very exciting building component to design. Reviewing the differences between the Eurocode and North American codes, such as NDS 2018, the discrepancies are observed, and they need to be addressed and researched. Finally, there is a great variety of information and research on TCC connectors; therefore, a lot of further research can emerge from this paper, specifically with different connectors, finite element software analysis, and experimental testing.

Conclusions
This paper was intended to research TCC systems and their applications while focusing on FEM analysis and vibration performance of standard TCC systems in residential and office buildings, and potential problems that may occur related to vibration design. In addition to the proposed models for the calculation of vibration, velocity, and acceleration, analytical models for calculating different types of slip modulus are given. This is an important detail because the degree of coupling conditions the stiffness of the system and, thus, the vibration conditions themselves. The composite timber-concrete system is an efficient system that is applicable during new construction or restoration, and consists of a monolithic concrete slab connected to timber beams. Comparing the load-bearing capacity of this composite system with the load-bearing capacity of the constituent elements, it can be concluded that the system has up to four times higher stiffness and up to twice the load-bearing capacity. For wide-span ceilings, the vibration design is often governing. The advantage of the lighter material of timber over concrete becomes a disadvantage, because a high mass material is advantageous for a low natural frequency. Having long stiff support beams and shorter floor spans can ensure improved vibration performance when compared to systems with longer floor spans and shorter beam spans. An assessment and calculation of the TCC system was made, which was satisfactory for the frequency, additional acceleration, and stiffness limit criterion. Additionally, constructive design requirements are fulfilled. Finally, FE analysis showed that stresses occur locally within the connector area in the concrete slab on both sides of the joint. The advantages of TCC (vs. sole timber) that can be emphasized are increased stiffness through composite action, increased floor mass at decking level, improved sound insulation (airborne sound), and reduced sensitivity concerning vibrations. If, on the other hand, the TCC is compared with pure concrete, it can be concluded that the weight is reduced, the CO 2 emissions are reduced, the building process is faster, and reduced effort is needed for the props and formwork. Data Availability Statement: Data available on requzest due to restrictions eg privacy or ethical. The data presented in this stury are available on request from the corresponding author.