Next Article in Journal
Simplified Analytical Methods for Prefabricated Concrete Wall Panel Building System with Alveolar-Type Joints
Previous Article in Journal
Applying the Tier II Construction Management Strategy to Measure the Competency Level among Single and Multiskilled Craft Professionals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulation of Load–Slip Capacity of Timber–Concrete Connections with Dowel-Type Fasteners

by
Dragan Manojlović
1,
Andrija Rašeta
1,*,
Vladimir Vukobratović
1,*,
Arpad Čeh
2,
Ljiljana Kozarić
2,
Đorđe Jovanović
1 and
Anka Starčev-Ćurčin
1
1
Faculty of Technical Sciences, University of Novi Sad, Trg Dositeja Obradovica 6, 21102 Novi Sad, Serbia
2
Faculty of Civil Engineering Subotica, University of Novi Sad, Kozaracka 2a, 24000 Subotica, Serbia
*
Authors to whom correspondence should be addressed.
Buildings 2023, 13(5), 1171; https://doi.org/10.3390/buildings13051171
Submission received: 21 March 2023 / Revised: 25 April 2023 / Accepted: 25 April 2023 / Published: 28 April 2023
(This article belongs to the Section Building Structures)

Abstract

:
Quality assessment of stiffness and load-carrying capacity of composite connections is of great importance when it comes to designing timber–concrete composite structures. The new European regulation intended explicitly for timber–concrete structures has made a significant contribution to this field, considering that until today there was no adequate design standard. Due to the proposed general expressions for determining the stiffness and load-carrying capacity of composite connections made with dowel-type fasteners, which are incapable of describing most of the commonly applied fasteners, engineering, and scientific practice remained deprived of a quality assessment of the essential mechanical properties of the connection. In order to overcome this problem, this paper proposes a numerical model of the connection suitable for determining the whole load–slip curve, allowing it to estimate the stiffness and load-carrying capacity of the connection. The model was developed by considering the non-linear behavior of timber and fasteners, which is determined through simple experimental tests. For the numerical model validation, experimental tests were carried out at the level of the applied materials and on the models of the composite connection. Through numerical simulations, analysis of obtained results, and comparison with experimental values, it can be confirmed that it is possible to simulate the pronounced non-linear behavior of the timber–concrete connection using the proposed model. The estimated values of stiffness and load-carrying capacity are in agreement with the conducted experimental testing. At the same time, the deviations are much less than the ones obtained from recommendations given by the new regulation. Additionally, apart from evaluating the value and the simulation of the complete curve, it is possible to determine local effects, such as the crushing depth in timber and concrete, the fastener’s rotation, and the participation of forces in the final capacity of the connection.

1. Introduction

In the last few decades, the worldwide interest in timber–concrete composite (TCC) structures has increased because of their structural benefits, i.e., more rational construction is obtained, where a high degree of utilization of both materials can be achieved [1]. It is an adequate solution for renovating old timber floors and is suitable for new structures, such as floors and walls in buildings and timber bridge decks [2]. TCC elements usually comprise thin concrete plates connected to timber beams with shear connectors (dowel-type fasteners, steel-punched metal plates, split rings and toothed plates, notches, glued-in plates, notches with dowel-type fasteners) [3]. The function of a shear connection is to transfer shear forces between timber and concrete and prevent transverse separation. Under bending action, when mechanical fasteners (dowel-type fasteners) are used, adjacent layers of composite systems slide relative to each other in a partially constrained manner, and only partial composite action is developed. The behavior of TCC structures is quite complex because their analysis requires the consideration of the slip between timber and concrete [4].
The design of timber–concrete composite structures is based on the assumptions of the elastic coupling theory, which takes into account the slip between two materials [5]. Displacement/slip in the connection, i.e., the flexibility of the fasteners (stiffness of the connection), is introduced into the calculation using the value of the slip modulus. The efficiency of TCC structures is related to the quality of the connection system (through the stiffness of the connection) because their behavior directly affects the distribution of internal forces and structural deformations. Additionally, apart from stiffness, the load-carrying capacity of the fasteners, which ensures the transfer of shear forces in the shear plane, plays an essential role in the composite connection.
Ideally, the analysis of composite beams should be based on experimental load–slip curves that describe the complete behavior of the composite connection. However, the behavior of the connection is usually represented by a limited number of properties based on experimental curves, which are expected to represent the actual behavior of the connection. These properties are stiffness and load-carrying capacity. Therefore, the behavior of connections is most often described by a bilinear load–slip relationship, defined by the slip modulus (initial part) and the load-carrying capacity of the fasteners (final part) [3].
Due to the lack and inadequacy of the provisions in the existing regulations used in this field, Eurocode 5 (EN 1995-1 [6] and EN 1995-2 [7]) and Eurocode 4 (EN 1994-1 [8]), there was a necessity for an adequate standard, which resulted in CEN/TS 19103 [9], explicitly intended for TCC structures. It made a significant contribution to this field, considering that, so far, there was no specific standard for TCC structures. Due to the proposed general assessments (for stiffness and load-carrying capacity), which are incapable of describing most of the frequently used dowel-type fasteners (nails, screws, bolts, dowels, and rods), only some observations that were not considered in the proposed expressions are given below. According to CEN/TS 19103, the evaluation of the connection stiffness (slip modulus Kser) for all considered dowel-type fasteners is still determined as a double value of the slip modulus given by empirical expressions for steel-timber connections, neglecting deformations in concrete [1]. It should also be noted that CEN/TS 19103 expressions were developed for the ductile type of failure, i.e., only for the case where two plastic hinges are formed in the fastener, with an empirically based approximation of the instantaneous deformation at 40% of the load-carrying capacity of the nailed connection [10]. According to [11], the slip modulus is a function of the fastener embedment depth in timber. Different failure modes result in different values of the slip modulus for different embedment depths. When considering the different types of concrete, such as normal weight concrete (NC) and lightweight aggregate concrete (LWAC), in accordance with CEN/TS 19103, the proposed expressions for LWAC concrete may lead to an inaccurate estimation of the slip modulus value [12,13] since local failure effects in concrete are not considered. On the other hand, in expressions for the load-carrying capacity of connections (Fv,R), different types of mode failure are taken into account by using expressions for timber–timber connections, with the approximation of the deformation in concrete based on the local strength of the concrete through a coefficient β [3]. In this way, the influence of local failure effects on the concrete side on the load-carrying capacity of the connection (local crushing of concrete) was only partially assessed. In contrast, the effects of pull-out [13] and cone expulsion [14] for the fasteners with heads were not considered, which indicates that LWAC concretes are not fully covered. When determining the connection load-carrying capacity, the rope effect contribution is considered in increasing the connection load-carrying capacity through the timber withdrawal capacity, provided that the fasteners are anchored adequately in concrete and ensured respecting the provisions of EN 1992 [15]. The general assessment proposed through empirical expressions in CEN/TS 19103 cannot comprehensively describe the behavior of a wide range of timber–concrete connections when different types of dowel fasteners and types of concrete are used. Due to the observation mentioned above, it is reasonable to expect that, in the next generation of CEN/TS 19103, a division will be made for different types of fasteners and consideration of the influence of local failure effects in concrete, which will affect the innovation and expansion of expressions for evaluating stiffness and load-carrying capacity. In the lack of appropriate proposals, engineering practice is always left with the possibility of evaluating the properties of composite connections via experimental tests or consulting the available applicable literature.
In a linear analysis of composite beams, stresses and deformations are usually calculated based on the linear elastic behavior of materials and connections. Based on these assumptions, for most practical problems, linear models offer sufficient accuracy and reliability, decreasing the need for more complex non-linear models. When timber–concrete connections are assembled with dowel-type fasteners, due to their pronounced non-linear behavior, assumptions about the linear elastic behavior of the connections and materials do not provide the possibility of including non-linear phenomena (the elastic non-linearity and the plastic deformation capacity of the joint) and can lead to non-conservative results. In order to overcome this problem, non-linear analyses have been increasingly used, taking into account the non-linear behavior of the connections and materials [16,17,18,19]. Therefore, the behavior of composite connections with dowel-type fasteners cannot be represented by a simple relation as an ideal elastoplastic behavior that is described only by the fastener’s slip modulus and load-carrying capacity. This problem is precisely the weakness of numerical models for non-linear analysis. It is essential to use more sophisticated models that could describe the non-linear behavior of the timber–concrete connection and define the whole load–slip curve (e.g., finite element (FE) models [20,21,22,23]).
Due to the necessity for numerical simulations of composite beams and considering the pronounced non-linear behavior of connections in relation to complex FE models, most research was based on the development of analytical models for determining complete load–slip curves [24]. These models are based on exponential functions, which were mainly developed to describe the behavior of timber [25], fasteners [26], and timber connections. However, due to their wide application in other fields, the possibility of their application in timber–concrete connections as descriptive models were also proposed [27,28,29,30,31]. The limitation of the models is the finite number of parameters (stiffness and load-carrying capacity), previously adjusted to experimental values, which describe the mechanical behavior of the composite connection. To overcome this problem and to reduce the extent of experimental tests, Dias [28] applied empirical expressions for Kser and Fv,R from Eurocode 5 in the models and enabled a more comprehensive application of descriptive models for timber–concrete connections made with dowel-type fasteners. Based on the material and geometric properties of the connection, predictive models were mainly developed in the wooden connections field [32,33,34,35], while such models do not exist for timber–concrete connections. Researchers Heine and Dolan [36], Laszlo Erdodi [37], Dominguez [38], and Kobel [39] made outstanding contributions to the characterization of timber connections. They proposed numerical predictive models to determine the whole load–slip curve based on the European Yield Model (EYM), taking into account the non-linear behavior of the timber and the fastener (determined from experimental tests). These models have shown very good results and, with certain modifications, could also be suitable in the case of timber–concrete connections.
In order to adequately predict pronounced non-linear behavior and to evaluate the mechanical properties of the composite connection required for the calculation of the TCC structures, a novel approach is needed. Therefore, there is a need for a non-linear connection model that will cover most of the abovementioned problems. This paper presents a suitable numerical model for describing the complete load–slip curve. The model is based on the mechanical characteristics of timber, concrete, and fasteners and includes the non-linear behavior of composite connections. The proposed model describes the complete behavior of the composite connection. It was developed so that instead of the assumed idealized behavior of the applied materials, their non-linear behavior is taken into account, which was determined by simple experimental tests.

2. Proposed Numerical Model

Apart from the model proposed by Dias [28], which is descriptive, no model has been developed that predicts the load–slip relationship over the entire domain in timber–concrete connections based on the material and geometric characteristics of the connection. Therefore, those above-mentioned predictive numerical models [36,37,38,39] represent the basis for the formation of a new model.

2.1. Background Theory of the Numerical Model

The EYM is a commonly accepted procedure for determining the capacity of timber connections based on the yield theory assumptions that consider the rigid-plastic behavior of the fasteners and timber. It represents the equilibrium conditions on an isolated segment of the fastener for possible types of failure in the connection. It implies that the failure of the connection is ductile and caused by the yield of the fastener or timber failure. Requirements for minimum dimensions of timber elements as minimum edge and spacing distances of fasteners usually guarantee this condition. It should also be noted that the load-carrying capacity of the connection determined by the EYM represents the final value of the capacity of the connection at failure caused by the yield of the fastener and timber failure. This theory does not cover the gradual development of timber and fastener damage, so quantifying displacement is impossible.
Additionally, the EYM theory implies a specific state when a failure in the connection occurs, i.e., a condition in which the simultaneous formation of plastic hinges and yield of timber and crushing of concrete arise under the assumption shown in Figure 1a. For a particular type of connection failure (failure modes), the EYM determines the position of the maximum moments in the fastener, which also represents the location of the plastic hinges or the location of the fastener rotation. Based on the assumed development of the connection failure (Figure 1a), it is possible to determine the complete load–slip curve (Fv,R-δ) on the considered displacement domain (Figure 1b).
As the failure modes are correlated with the geometrical characteristics of the connection and the mechanical properties of the timber and the fastener, it is possible to predict the failure mode that will occur (Figure 2).
In the case when the connection is assembled with different materials (represented by the embedment strength of the materials fh,tfh,c), the transition from Mode II to Mode III must satisfy the requirements regarding the embedment depth of the screw in timber (lt,req) and concrete (lc,req) [40]:
l t l t , r e q = 1.15 · 2 · β 1 + β + 2 · M y f h , t · d
l c l c , r e q = 1.15 · 2 1 + β + 2 · M y f h , c · d
By using timber screws, the effect of the screw thread in timber provides a certain degree of anchoring, which creates the rope effect (Figure 3), i.e., due to the slip in the connection and prevented extraction of the screw from timber, a tensioning force occurs in the screw. Prevented extraction of the screw leads to the adhesion of timber and concrete, and an additional normal force, Ht, occurs at the contact, where its horizontal component (μ·Ht) contributes to an increase in the load-carrying capacity of the connection. Due to the tension force in the screw Fax(θ), its horizontal component also contributes to an increase in the load-carrying capacity of the connection by activating a larger timber crushing zone (outside the plastic hinges zone). In CEN/TS 19103, the contribution of the rope effect to an increase in the load capacity of fasteners (where anchoring is provided) is considered according to the recommendation given in Eurodoce 5 (8.2.2(2)), i.e., as an additional force Frope (min{Fax/4;100%·FJoh}) in relation to the Johansen part (FJoh), where Fax is the withdrawal capacity in timber, and 1/4 represents the friction coefficient μ = 0.25 (corresponding to friction between two timber elements) [41]. In the timber–concrete connection, according to experimental tests, the coefficient of friction between concrete and timber is μ = 0.45–0.66 [42,43]. Using the simplified EYM model and based on the experimental tests, Hao Du et al. [44] concluded that the angle of rotation of the screw has a significant part in contributing to the load-carrying capacity due to the rope effect. The coefficient of friction is in the range of 0.35–0.50 (μ = sin(90 − θ)) for an angle θ that varies from 20 to 30 degrees (for the slip in connection up to 15 mm). According to previous experimental research [3], for a slight slip in the connection (below the yield load point), the contribution of the load-carrying capacity due to the rope effect on the connection load-carrying capacity is negligible. On the contrary, when large slips occur, the rope effect significantly increases the final load-carrying capacity. In the expressions according to Eurocode 5, as the rope effect is determined regardless of the angle of the screw rotation, the force contribution due to the rope effect and friction cannot be perceived and determined in the final load-carrying capacity of the connection.
For the model formulation, the following assumptions and approximations were introduced:
  • Timber and concrete may be considered deformable bodies into which the fastener is pressed. The fastener rests on a deformable elastic foundation in timber and concrete and can be represented as a beam on which embedded stresses act (fh,t and fh,c);
  • The gradual development of the position of plastic hinges or the location of rotation of the fastener can be estimated by the EYM model for Mode III, i.e., considering the condition of equilibrium on an isolated screw segment, as shown in Figure 4c;
  • It is assumed that the embedment stresses (fh,t and fh,c) are uniformly distributed in the zone between the plastic hinges;
  • When pressed into timber and concrete, the fastener behaves like a rigid body, except for an infinitesimally small region where plasticization of the fastener (plastic hinge) occurs, around which rotation is achieved. This assumption results in a linear function of displacement in timber and concrete [36,39];
  • In the zones outside plastic hinges, the fastener remains in its initial position since, according to EYM, the fastener outside plastic hinges is balanced, while the deformations perpendicular to the axis of the fastener can be ignored;
  • The mechanical characteristics of timber (embedment strength in timber → fh − Δ) and the fastener (fastener yield moment → Myθ) are presented through non-linear behavior, i.e., determined experimentally. These characteristics are required to describe the behaviors of timber and the fastener as a function of slip in the shear plane;
  • Deformation in concrete is idealized through the ratio β = fh,c/fh,t, i.e., through the function of displacement in timber δc = δt/β, taking into account the deformations in concrete (concrete crushing as a local failure effect adopted as fh,c = 4·fc [45]);
  • The withdrawal capacity in timber was introduced through a sine function of the angle of rotation of the fastener θ (Equation (3), Figure 3). The withdrawal capacity Fax,max is determined via experimental testing or obtained using empirical expressions according to Eurocode 5. It only applies to failure modes involving fastener deformation (plastic hinges);
    F a x θ = F a x , m a x · sin θ
  • The fastener is perpendicular to the shear plane;
  • Although a plastic barrier (plastic sheet) was placed in order to minimize the effects of friction on the response of the connection and to protect the timber from the influence of moisture from concrete, the effect of friction in the connection (due to rope effect) was introduced into the calculation [43]. The coefficient of friction between timber and concrete adopted in this paper was assumed as μ = 0.45, for the case of applying plastic sheet, according to the suggestion from [43].

2.2. Analytical Procedure and Calculation Technique

The analytical procedure for determining the load-carrying capacity of a connection is divided into four phases:
Phase I—Determination of the failure mode of connection:
Applying Equations (1) and (2) determines the failure mode depending on the required embedment depth lreq. Expressions for the required embedment depth were derived based on the two considered failure modes (Mode II and Mode III), demonstrating the transition from one failure mode to another (Figure 2).
Phase II—Determination of the embedment depth in timber and concrete:
For the obtained failure mode determined in Phase I and using the EYM (taking into account the above simplifications and the equilibrium state on an isolated segment of the fastener), the positions of the maximum bending moments in the fastener are obtained, i.e., the crushing depth in timber bt (Equation (4)) and concrete bc = bt/β, Figure 4. Under the assumption of the rigid-to-plastic behavior of timber and the fastener, the first estimation of the position was determined for the maximum values of the embedment strength (fh,t,max) and yield moment of the fastener (My,max), which should be determined through experimental tests or calculated using the available empirical expressions given in Eurocode 5.
b t = 2 · M y + F v , a x θ · S f h , t · d e f f · 2 · β 1 + β
Phase III—Determination of the angle of rotation of the fastener:
For a given displacement S of the connection, the indentations in timber and concrete (δt and δc) are determined, as shown in Figure 4a. Based on the calculated indentation values and estimated crushing depths (bt and bc from Phase II) in timber and concrete, the angle θ is determined from the geometry of the considered deformed screw and according to the expression given for the angle θ shown in Figure 4a.
Phase IV—Determination of load-carrying capacity of a connection:
As the first estimation for bt and bc in Phase II was carried out based on fh,t,max and My,max, in order to estimate the angle θ, the final crushing depths btk and bck, which are used in Equation (6), must be iteratively determined using experimentally determined functions fh − Δ and Myθ. The final crushing depths btk and bck are determined using Equation (4) and applying the values determined in Phases II and III. If btk and bck values differ by more than 1% from the estimated bt and bc from Phase II, it is necessary to repeat the determination of the rotation angle, taking the recently calculated values (btk and bck) as input data for Phase III. Due to the impossibility of including the share of force Fh,ax(θ) in the model due to the unknown displacement in timber outside the plastic hinge, to satisfy the equilibrium condition (Figure 4c), Fh,ax(θ) is added to Fy,R afterward. As tension force occurs in the screw from the beginning, the force at the shear plane Fμ = μ·Fv,ax(θ) represents the share of friction in the connection. The friction force does not affect the moment equilibrium, which implies that Fμ can be added to Fy,R afterward. The load-carrying capacity of the connection is determined according to Equations (5) and (6).
F v , R = F y , R + μ · F v , a x θ + F h , a x θ
F y , R = b t k · f h , t · d e f f = b c k · β · f h , t · d e f f
The above equations represent the algorithm for calculating the load-carrying capacity of a connection for one screw. For obtaining the complete load–slip curve on the entire domain, an iterative procedure must be carried out for predefined displacements S in connection. The maximum displacement δmax on the considered domain in the proposed model in the amount of 15 mm was adopted according to Dias [28], which stated that for this slip, the maximum load is always reached, and the assumption is only valid for the TCC connection made with dowel-type fasteners. For a more straightforward interpretation of the proposed procedure, a summary of the calculation process can be seen on the flow chart in Figure 5.

3. Materials and Methods

3.1. Material Properties

The first step required to obtain the complete load–slip curve of the connection with the proposed numerical model is the determination of the mechanical properties of applied materials. Based on the proposed algorithm, two sets of material properties are needed for the analysis. The first set of material properties refers to the determination of mechanical properties of materials, i.e., the ultimate strength of screw fu, the withdrawal capacity in timber Fax, and the compressive strength of concrete fc. The second set involves determining the non-linear relationship of the embedment in timber (fh − Δ) and yield bending moment in screw (Myθ).
In this paper, shear test specimens were made of timber lag screws 10 mm × 150 mm (DIN 571 hexagon head timber screws, Grade 5.6, as specified by the manufacturer), lightweight aggregate concrete LWAC LC45/50 (density class D1.8) and glulam GL24 h (ρm = 477 kg/m3). The structural lightweight concrete was made with a combination of riverbed aggregate (size 0/4 and 8/16 mm) and Liapor HD 1/8 lightweight aggregate incorporating 40% (by volume) of fresh concrete. The hardened concrete had a dry density of 1740 kg/m3, compressive strength of fc = 55 MPa and modulus of elasticity E = 25.3 GPa. The timber test specimens were tested after conditioning at relative humidity of 65 ± 5% at a constant temperature of 20 ± 2 °C, leading to a timber moisture content of around 12%. Experimental tests were carried out for the first set of material properties, and in Table 1, the mechanical characteristics of the materials and the applied test standards are shown.
For the second set of material properties, in order to determine the non-linear relationship between the timber embedment strength fh − Δ and the yield bending moment of the screw Myθ, experimental tests were carried out according to EN 383 [49] and ISO 10984-1 [50], respectively. The test setups are presented in Figure 6a and Figure 7a, respectively.
To determine the embedment strength of timber parallel to grain, according to EN 383, the dimensions of the timber specimens were 30 × 60 × 140 mm, where the thickness t = 30 mm was adopted to minimize the effect of screw bending during the test. Five specimens were made for this test (Figure 6b). After processing the timber specimens, holes with a diameter din of 7 mm were drilled, and screws were inserted into the timber specimens using a wrench. This ensured a direct contact between the core diameter of the screw and timber, as will be achieved in a composite connection.
All timber specimens were monotonically loaded parallel to grain with a preloading cycle, while the load was adjusted following the loading procedure given in EN 383. The load increased or decreased at a constant rate of loading-head movement through a universal hydraulic press with a capacity of 2000 kN (P-250, Milaform-Service, Neftekamsk, Russia). The applied load was controlled using a 50 kN compression load cell with a measurement accuracy of 0.01 kN (CZL312A-5t, SAH Electronics, Belgrade, Serbia). Indentation of screws in timber specimens is measured using a Linear Variable Differential Transformer (LVDT) with a maximum deflection range of 10 mm and an accuracy of 0.01 mm (NOVOTECHNIK TR-0010). The measurement was made at half the timber specimen height, i.e., LVDTs were attached to the external steel arm while its tip rested on a small steel L profile screwed into the timber specimen at the level of the center line of the screw. Simultaneous data acquisition from LVDTs and load cells was conducted at a rate of 10 Hz using the “SPIDER 8” acquisition system. The test was stopped either when the maximum load was reached or when the deformation was 5 mm.
The embedment–slip curves (with labels FH1Z to FH5Z, the same as the name of the tested timber specimens) were obtained from the recorded measurements (load and indentation) and presented in Figure 6c. For the purpose of the proposed numerical model, for all tested embedment–slip curves, the initial slips and preload cycles (from 0.4·Fmax,est to 0.1·Fmax,est) were cut out from the curves presented in Figure 6c. In this paper, the effective diameter of the screw deff = 7.7 mm (deff = 1.1·din) was used to determine the embedment strength. The reason for taking the effective diameter of the screw arrived from the suggestion given in Eurodoce 5, i.e., in the case of screws, it is necessary to consider the thread’s positive effect using the effective diameter of the screw [51]. When the effective diameter of the screw is considered in the proposed expression for fh in Eurocode 5 (fh,EC5 = 36.1 MPa), a better agreement with the experimentally determined mean value of the embedment strength fh,mean = 35.3 MPa is obtained. Additionally, the same reference diameter of the screw must be used to derive the embedment strength of timber and load-carrying capacity given with Johansen’s expressions [51].
In ISO 10984-1, two suggested methods (Method A and B) are given to determine the yield moment and noted that the choice is left to the user to choose the method, whichever is most relevant. The experimental study conducted in [34] shows significant differences between the two test methods. These differences are reflected in terms of the test setup and determination of the load necessary to determine the yield moment. Within the EN rules of safety, which avoids reaching the actual collapse of connection, the determination of yield moment according to suggested methods leads to the yield moment that lies between the elastic and plastic moment according to the mechanical relationship [52]. This fact indicates that only the outer areas of the cross-section of the screw are under plastic deformation, and the plastic bending moment is only partially utilized in real connections, i.e., the plastic bending capacity of the screw in connection will never reach its maximum values. As mentioned above, for the proposed model, it is essential to determine the real behavior of the yielding moment of the screw in the connection, i.e., the relationship Myθ cannot be calculated using a simple empirical expression. Nevertheless, it must be determined in a bending test. For this purpose, both methods can be applied. Depending on the slenderness of the fasteners and the method of estimating the rotation angle, Method B was applied in this paper.
The procedure in accordance with ISO 10984-1 (Method B) was applied to determine the yield bending moment of the screw, via three-point bending test, with a distance of supports equal to L0 = 80 mm (Figure 7b). The test was carried out on 5 specimens, where the force was applied at a speed of 2 mm/min. Using recorded data of the screw displacement in the middle of the span u and the applied force P, the relationships Myθ (with labels My1 to My5, for all tested screws, Figure 7c) were determined according to Equation (7) for the yield moment, and Equation (8) for the rotation angle.
In this study, for the comparison with Eurocode 5, the yield bending moment of the screw was obtained using both methods. According to Method B, the yield bending moment My,0.05% = 41,369 Nmm was obtained using the yield force P0.05% = 2.04 kN, determined by the 5% offset method, which uses the core diameter of the screw (din = 7 mm) [51]. According to Method A, yield moment My,max = 46,613 Nmm was determined based on the maximum force Pmax = 2.33 kN. Using the empirical expression (8.14), given in Eurocode 5, the yield moment is My,EC5 = 42,071 Nmm (with deff = 7.7 mm and fu = 695 MPa). The deviation of −9.7% in relation to My,max is a consequence of the adopted assumption of ultimate strength in bending fub = 0.8·fu used for bolts and dowels in Eurocode 5 [26], while a deviation of −1.7% to My,0.05% shows good agreement with experimental one. Based on the yield moment My,max, the ultimate strength of the screw from the bending test can be estimated with Equation (9), based on the plastic moment according to the mechanical relationship via assuming full plastification of the screw cross-sections. In Equation (9), the effective screw diameter deff is applied [53]. The ultimate strength of the screw in bending was fub,max = 612 MPa (~0.9·fu). The ultimate strength in bending is in the range of the yield strength in tension fy and ultimate tensile strength fu. In order to allow a more realistic estimation of the maximum yield moment, if the recommendation according to [54] was adopted (fub,eff = 0.9·fu = 625 MPa, for fu > 450 MPa), the yield moment My,eff = 47,593 Nmm would have a deviation of +2.1%, which is in the interval of the experimentally obtained yield moment (My,min = 45,954 Nmm and My,max = 47,425 Nmm). The ultimate effective strength of the screw in bending fub,efffub,max agrees with experimental values.
M y = P · L 0 / 4
θ = 2 · tan 1 2 · u L 0
M p l = f u b · d e f f 3 / 6                 f u b = 6 · M y / d e f f 3

3.2. Push-Out Test Specimens, Loading Procedure, and Results

3.2.1. Description of Shear Test Specimens

For the proposed model validation, all push-out composite connection specimens were made according to the materials mentioned in Section 3.1 (Table 1). According to CEN/TS 19103, the suggestion of the specimen configuration, testing protocol, and determination of mechanical properties of the TCC connection is given in Annex C. The suggestion that test specimen dimensions can be defined based on a minimum admissible edge and spacing distances allowed for each specific fastener according to Eurocode 5 is unsuitable for TCC connections. As there are no specific criteria regarding the influence of test specimens on the experimental characterization of a TCC connection, the author of [55] proposed two configurations of the symmetrical type of connection (concrete–timber–concrete and timber–concrete–timber), as well as minimum dimensions of the specimens from the requirements of ductile behavior of the connection. Considering the mentioned problems, in this paper, the dimensions of the timber part of the specimens were derived from the criteria of minimum screw distances according to Eurocode 5, adopted as double the minimum value prescribed for timber connections. The screw distance of at least 150 mm (a1 ≥ 13·deff·n0.4 = 131 mm) was adopted based on the criterion of the effective number of fasteners in 1 row to provide that neff = n = 2 [56]. Doubling the value of minimum distances was adopted in order to ensure that the influence of the group of fasteners, the uneven redistribution of the shear force between the fasteners, and the possible splitting of timber perpendicular to the fibers is avoided, which is implicitly contained by the effective number of fasteners according to Eurocode 5. To ensure the ductile-type failure mode in the connection and the development of two plastic hinges, the embedment depth of the screw in timber must be at least lt = 9·din [11]. Considering the embedment depth of the screw in concrete, the minimum depth must be lc = 3·ds [9,57]. According to the recommended values, the screw (Φ10/150 mm) satisfies the minimum installation criteria since the depth in timber was lt = 100 mm and concrete lc = 50 mm.
In this study, double shear test specimens were adopted (concrete–timber–concrete), consisting of 1 central timber (160 × 240 × 400 mm) and 2 side concrete (160 × 100 × 400 mm) members connected with 2 timber screws 10 × 150 mm in each shear plane at a distance of 150 mm (Figure 8a). Concrete was reinforced with steel mesh reinforcement Φ4/50 mm in the bottom zone at 20 mm from the timber member (Figure 8b). Three shear test specimens were made for this study (Figure 8c). Before forming the test specimens and placing them in the formwork, the timber member was conditioned at 20 °C/65%, and during the shear test, the timber humidity was 10.9%. After pre-drilling the holes with a core diameter of 7 mm (din) in timber, the screws were installed manually with a wrench. A plastic sheet was placed between timber and concrete, which reduces the effects of friction and prevents water penetration into the wood from concrete.

3.2.2. Test Setup and Loading Procedure

Short-term push-out tests were conducted according to EN 26891 [58]. In Figure 9a, the shear test setup is shown. The force was applied to a timber member through a universal hydraulic press with a capacity of 2000 kN (P-250, Milaform-Service, Neftekamsk, Russia). The applied force was controlled using a 200 kN compression load cell with a measurement accuracy of 0.01 kN (CZL110D, SAH Electronics, Belgrade, Serbia). Displacements between timber and concrete were measured using a Linear Variable Differential Transformer (LVDT) with a maximum deflection range of 50 mm and an accuracy of 0.01 mm (SHTSEIKO KTR-50 mm). For each test specimen, two LVDTs were placed near the shear plane to register the displacements as precisely as possible and to observe any uneven movements. The measurement was made at half the model height, i.e., at half the distance between the screws. LVDTs were attached to the timber member of the test specimen, while its tip rested on a small steel L profile screwed into the concrete. Simultaneous data acquisition from LVDTs and load cells was conducted at a rate of 10 Hz using the “SPIDER 8” acquisition system.
The shear tests were conducted based on the loading protocol from EN 26891 (Figure 9b). At first, the test specimens were loaded up to 40% of the estimated ultimate load (Fest). The applied load was kept constant for 30 s. After that, the load was reduced to 10% Fest and kept constant for 30 s. The load was increased to 70% Fest with a rate of 20% Fest per minute. Finally, the displacement control mode was adopted to apply the further load until the occurrence of either failure or a 15 mm displacement.

3.2.3. Push-Out Test Results

Shear test specimens were tested according to a protocol given in EN 26891 that requires the estimated maximum load Fmax,est based on which the loading is carried out. In this study, the first test specimen was tested to accurately assess the maximum load (Fmax). The other two test specimens were tested according to the adopted loading procedure in Figure 9b. Although the first test specimen was used for a more accurate estimation of the maximum load, using Johansen’s expression for Mode III and experimental values of the mechanical characteristics of the material, the calculated estimated load (in the interval of 40.7 to 55.6 kN) was in agreement with the first tested specimen. As EN 26891 prescribes that the test should be carried out when the ultimate load is reached or when the slip is 15 mm, the first test specimen was tested up to a slip of ~30 mm, while the remaining 2 test specimens were tested up to the slip of 15 mm. This ensures that the maximum capacity is reached and shows the composite connection behavior beyond the proposed slip. The experimentally obtained load–slip curves for all tested specimens (with labels LWAC-1 to LWAC-3) are shown in Figure 10.
To show the failure mode of the screw and the damage to timber and concrete, parts of the timber and concrete were removed after the test. All tested specimens showed ductile behavior, i.e., during the screw indentation in the timber, a maximum embedment strength was reached, and two plastic hinges were formed in the screw (Figure 11a), as was expected at the given outset in the design of specimens to ensure a ductile-type failure of the connection. Although the first specimen was tested for slips beyond 15 mm (Figure 11b), it can be observed that there is no crushing in concrete around the screw, shear failure of the screw, or timber splitting. Figure 11c shows that the screw is pressed into the concrete and that the proposed approximation of neglecting the deformations in concrete is inappropriate. This fact implies that attention should be paid to this assumption, especially when LWAC with lower density classes than D1.6 are applied [12].

4. Validation of the Proposed Numerical Model

4.1. Numerical Simulations

In order to determine and evaluate the behavior of the composite connection using the proposed numerical model, a code was written in MATLAB according to the given algorithm. Based on the experimental curves for the embedment strength of timber and yield moment of the screw, numerical curves for simulations of the TCC connection were formed. Numerical curves for embedment strength fh − Δ (with labels fh1z to fh5z) were obtained based on the assumption that hardening after the yielding was negligible for the indentations above 5 mm, i.e., the horizontal branch of the curve was added to the experimentally determined curves (Figure 12a). As the experimentally obtained yield moment curves Myθ (My1 to My5) are uniform in behavior, the mean curve My-NUM (Figure 12b), with a horizontal branch after the maximum moment, was adopted for the numerical curve of the screw yield moment. Five numerical simulations (with names SIM(fh1z) to SIM(fh5z)) were carried out to validate the proposed model, in which the numerical embedment curves fh − Δ (fh1z to fh5z) were varied, while the mean curve of the yield moment of the screw Myθ (My-NUM) was the same for all simulations. All simulations were carried out up to the slip of 15 mm, and Figure 12c shows superimposed experimentally (LWAC-1 to LWAC-3) and numerically (SIM(fh1z) to SIM(fh5z)) obtained load–slip curves Fv,Rδ to compare the results.
In addition to the numerically determined curves shown in Figure 12c, the following diagrams for simulation 3 (SIM(fh3z)) are superimposed and shown in Figure 13a: bt (crushing depth in timber), bc (crushing depth in concrete), θ (screw rotation angle), My (yield moment of the screw), and fh,t (embedment strength in timber) as functions of the connection displacement δ in the shear plane. To perceive the contribution of forces to the increase in the load-carrying capacity of the connection (in addition to the final value Fv,R, and Johansen’s part FJoh without rope effect), the following contributions of forces are shown in Figure 13b:
  • The influence of the axial force in the screw Fax,f = Frope′ + Frope″, where Frope′ = Fy,RFJoh is the force that increases the indentation zone in timber when withdrawal capacity is activated (increasing bt), while Frope″ = Fh,ax(θ) is the force that represents a further expansion of the indentation zone beyond the plastic hinge in timber;
  • The effect of friction force Fμ = μ·Fv,ax(θ);
  • The combined effect of axial and friction forces Frope = Fax,f + Fμ, i.e., the “rope effect”.
Observing the local effects in the connection behavior (Figure 13a), it can be seen that after 4 mm of slip, the crushing depth in timber bt rapidly increases due to the rope effect. In contrast, the crushing depth in concrete is approximately constant bc = 5 mm (see Figure 12c). The crushing depth bt, derived from the equilibrium condition, shows the activation of a larger zone of timber that resists the movement/rotation of the screw due to the activation of the axial force and the change in the position of the plastic hinge in timber. In reality, changing the position of the plastic hinge is not possible, and this can be interpreted in a way that the plasticization and the moment increase in the screw spreads in a larger zone around the plastic hinge, which can be seen in Figure 11a,b. Additionally, Figure 11a,b show that even for displacements larger than 15 mm, the position of the plastic hinge remains unchanged and corresponds to a length of 25–35 mm obtained within the numerical simulation. The assumption that the fastener remains in its initial position in the zones outside the plastic hinge was experimentally confirmed for displacements up to 15 mm (Figure 11a). It can be noticed that for a displacement of δ = 15 mm, considering all simulations, the angle of rotation is θ ≈ 21°, and in agreement with the experimentally measured value of θ ≈ 22° (Figure 11a). It can also be seen that plastification of the screw was achieved with the formation of two plastic hinges, where the maximum value of yield moment My,max is reached. From Figure 13a, it can be seen that for a slip greater than 2 mm, due to the formation of the second plastic hinge in timber and to reaching the yield point of the screw My and timber fh,y, the connection shows a pronounced non-linear behavior, which is expected for these fasteners.
In Figure 13b, a significant increase in the load-carrying capacity of the connection due to the rope effect can be observed. Due to the slip increase in the connection and activation of the withdrawal capacity in timber, the gradual development of the rope force contribution Fax,f can be observed. It can be noted that this contribution gradually develops before the yield point of the connection is reached and that their further increase is approximately linear. As the friction force Fμ is the function of Fv,ax(θ) and friction coefficient μ (which has been adopted as constant), its contribution is developed from the very beginning and follows the trend of axial force in the screw Fax(θ). For the maximum rotation angle of screw θ = 22°, the contribution of the rope effect to the increase in the bearing capacity of the connection according to the numerical simulation is Frope = 6593 N (Figure 13b, purple dashed curve), which corresponds to the force that represents the limit in Johansen’s expression for the rope effect Frope = min{100%·FJoh = 6438 N; μ·Fax,max = 0.45 · 15,525 = 6986 N} = 6438 N. Based on the results, the rope effect significantly contributes to the connection’s load-carrying capacity. The proposed approximation of axial force in the screw as a function of the rotation angle θ (sine function) shows a good agreement with experimental results and assumptions from Eurocode 5. For a slip of δ ≈ 19 mm (Figure 10), the maximum load-carrying capacity of the connection occurs, followed by a decrease in load capacity due to the screw being pulled out of timber, as well as the formation of plastic hinges in the screw (Figure 11b). Due to the unknown axial tension force in the screw, as a relationship between the withdrawal force in timber and slip in the shear plane, as well as due to the impossibility of taking it into account through a kinematic relation [39], for displacements above 15 mm, the numerical model would continue to increase the force. This fact represents one of the limitations of the proposed numerical model. This issue is noted as one of the crucial problems in the design of timber connections with configurations that use fasteners with threads, i.e., when the contribution of the rope effect is essential to determining the connection resistance [59]. By comparing numerical and experimental curves, numerical simulations showed that if the material characteristics of timber and the fastener are known (represented in the form of fh − Δ and Myθ), it is possible to simulate the non-linear behavior of the timber–concrete connection on the whole domain (δ = 0–15 mm) by using the proposed numerical model.

4.2. Validation of the Proposed Model via an Example from the Literature

In order to additionally validate the proposed numerical model, an experimental study of TCC connections conducted by Stevanović in his Ph.D. thesis [60] was selected. All available input data required for the numerical model were taken from the [60]. The double-shear push-out test specimens were used. Test specimens were formed of 1 central timber (100 × 200 × 350 mm) and 2 side concrete (70 × 300 × 350 mm) members connected with 2 timber screws Φ10 × 150 mm in each shear plane at a distance of 100 mm (see Figure 3.10 in [60]). Three replicants of the shear test specimen were made and tested according to EN 26891. Based on the conducted experimental tests, the load-carrying capacity of the connection was Fv,R,mean = 10.65 kN (Fv,R,min = 10.47 kN and Fv,R,max = 10.78 kN), while the slip modulus was Ks,mean = 2580 N/mm (Ks,min = 2140 N/mm and Ks,max = 2860 N/mm). For details of experimental research conducted by Stevanović, see Figure 4.4 for Fv,Rδ curves (test specimens Z1–Z3) and Table 4.1 for values Fv,R and Ks in [60].
Since experimental tests for embedment strength in timber fh − Δ and yield moment of the screw Myθ were not conducted, they were reconstructed by using exponential functions [33], as shown in Figure 14a,b. The input data (fh,t and Kf) for the exponential function fh − Δ were determined according to the suggestions from [61], i.e., using the expression from Figures 6 and 7 provided in [61], Equations (10) and (11), which was based on timber density ρm. Input data (My and Kb) for the exponential function Myθ were determined according to the suggestions from this paper based on conducted tests on the screws. Yield moment My was determined as the effective value according to the expression in Table 2, and Kb was determined as the mean value of the initial slope of the Myθ curves from Figure 7c. The withdrawal capacity in timber Fax was estimated according to expression (8.3), given in Eurocode 5. The coefficient of friction was adopted as μ = 0.45. The input data are summarized in Table 2.
f h , t = 0.0936 · ρ m 8.454
K f = 0.2183 · ρ m 60.513
The numerical simulation results are shown in Figure 14 and Figure 15. In order to compare numerical load–slip curves with experimental ones, Figure 14c shows superimposed experimental (black lines for experimentally tested specimens ZP1 to ZP3 [60]) and numerical (red lines, for Fv,R and Fv,eff,R(neff)) load–slip curves Fv,Rδ. It should be noted that in Figure 14c, the red line (Fv,R) represents the Fv,Rδ curve of the connection without the influence of the group effects (effective number of fasteners in a line parallel to the grain on the load-carrying capacity of the connection), while the dashed line (Fv,eff,R(neff)) takes into account the spacing between the screws through the effective number of screws, neff = 0.93, reducing the load-carrying capacity of the connection. With this simple suggestion from Eurocode 5, as well as observations and conclusions reported in [53], in the proposed numerical model, an attempt was made to consider the reduction in the load-carrying capacity of the connection through the effective number of screws. This recommendation refers to the final value of the load-carrying capacity and that further investigation of this influence on the complete Fv,Rδ curve is necessary to improve the prediction by the proposed numerical model. Comparing the numerical results with the deformation of the screw (see Figure 4.8 in [60]), the distance of plastic hinges bt,NUM = 20–30 mm, and the maximum angle of rotation of the screw θ = 22.6° follows the experimental values (bt,EXP ≈ 30 mm and θ ≈ 24°, according to [60]). Based on results from the presented numerical simulation, it can be concluded that the numerical model can sufficiently accurately capture the load–slip response, although not all experimental data were available, i.e., embedment strength in timber fh − Δ, the yield bending moment of the screw Myθ, and the withdrawal capacity in timber Fax.

4.3. Discussion of Results—Load-Carrying Capacity and Slip Modulus of the TCC Connection

In order to compare the numerical, theoretical, and experimental values of the load-carrying capacity and slip modulus of the TCC connection, Table 3 shows the estimations of the proposed model compared to theoretical values (CEN/TS 19103) and the experimentally determined values considering the mean values. Based on the comparison presented in Table 3, the load-carrying capacity of the connection obtained by the numerical model gives a better estimation compared to the recommendations by CEN/TS 19103. This deviation is a consequence of the proposal for the rope effect in CEN/TS 19103 given in the form Fax/4. This effect should be revised or taken in the form of μ·Fax, where if there is a plastic sheet or if timber is in direct contact with concrete, the coefficient of friction between timber and concrete should be experimentally obtained or taken according to the recommendation from [43]. If the value μ = 0.45 was adopted, then the load-carrying capacity of the connection according to Eurocode 5 would have been estimated better Fv,R,corr compared to the experimental one. It should be noted that in the case of Stevanović, the load-carrying capacity of the connection should include the influence of the group effects in order to take into account the reduction in the load-carrying capacity of the connection through the effective number of screws.
Compared to the experimental values, the mean value of the slip modulus determined from the numerical curves, according to EN 26891, shows a good prediction compared to overestimated theoretical values according to CEN/TS 19103 (Table 3). It should be noted that the diameter of the screw is not precisely defined in CEN/TS 19103, which may result in an overestimation of the Kser value. According to suggestions from the literature, effective diameter deff is used to estimate the slip modulus in this paper. Analyzing the expression given in CEN/TS 19103, for TCC connections Kser,TCC, as well as the case for timber-timber connection (TTC) Kser,TTC, the slip modulus for TCC connection has large deviations to the experimental values in comparison to the slip modulus for the TTC connection. Although in the case of concrete with a higher strength class that was used in this paper, based on the experimental results, the recommendation by CEN/TS 19103 to double the value of the slip modulus is not justified. Using numerical curves, slip modulus should be determined for 40% of the ultimate load, as in codes and standards, slip modulus for serviceability limit state design Kser is derived as the ratio of the respective load (0.4·Fv,R) and the instantaneous deformation (v40%) [10]. Compared to the EN 26891 methodology for determining Ks, in this study, the slip modulus was determined as a secant value at 40% of the ultimate load (Ks,0.4Fmax). This determination method is adopted because there is no initial slip, i.e., all initial slips have already been removed from the embedment curves when the numerical curves were formed, while EN 26891 takes this into account when determining the slip modulus. It should be noted that the load–slip curves show a pronounced non-linear behavior, which significantly affects the determination of the slip modulus value. As discussed in the paper [62], it can be concluded that a limit of 0.4·Fv,R is a sensitive parameter to obtaining the slip modulus of the connections. Thus, the evaluation of the slip modulus should be performed with caution. It was shown that the proposed numerical model has advantages in determining and evaluating the slip modulus of the TCC connection, so it is recommended that for TCC connections with screws, the slip modulus should be determined using the proposed numerical model. The load-carrying capacity and stiffness values of the TCC connection, determined by the proposed numerical model, agree with the experimental values, and maximum deviations are within acceptable limits according to [24].

5. Conclusions

To improve the prediction and assessment of the stiffness and load-carrying capacity of the TCC connection, in this paper, a numerical model of the connection suitable for determining the complete load–slip curve is proposed. The model was developed to consider the real non-linear behavior of timber and fasteners, which is determined via simple experimental tests. Since there are no detailed criteria in the CEN/TS 19103 regarding the geometry of test specimens and their influence on the response of TCC connections, all recommendations and observations were considered based on a literature review from the field of composite connections. Suggestions of minimum dimensions of shear test specimens, the optimal spacing between screws, and its effect on the load-carrying capacity of the connection, as well as the occurrence of a possible failure mode of the screw in the connection, are discussed in the paper. According to these suggestions, shear test specimens in this study were made. For the numerical model validation, experimental tests were carried out at the level of the applied materials, as well as on the specimens of the composite connection. For additional validation of the proposed model, experimental research of the TCC connection taken from the available literature was used.
Through numerical simulations, analysis of the obtained results, and comparison with experimental and theoretical values, it was shown that it is possible to simulate the non-linear behavior of the TCC connection using the proposed numerical model. It can be concluded that the numerical simulation carried out on the TCC specimen adopted from [60] proved to be successful, although all the experimental data needed for the numerical model was not available. This proves that in the absence of experimental tests, it is possible to use the generic functions for fh − Δ and Myθ, i.e., using exponential functions in which the input parameters are determined empirically [61,63], or to evaluate the material behavior better through experimental testing. This aspect is of great importance in order to extend the application of the proposed numerical model, which requires simple experimental tests at the material level. The stiffness and load-carrying capacity values of the connection estimated with the numerical model agree with the experimental ones, and the deviations are minor compared to the recommendations and overestimated values given by the CEN/TS 19103. Investigating and updating the recommendations related to slip modulus is necessary because it overestimates the experimental values by more than +140%. These large deviations result from the adopted assumptions of negligible deformation in concrete, the unclear diameter of the fastener with threads, and the significant influence of timber properties. It can be noted that the connection response strongly depends on the adopted embedment fh − Δ curve and that the foundation modulus significantly influences the initial slope of the Fv,Rδ curve, i.e., slip modulus. For the load-carrying capacity of the connection, the proposed expression in Eurocode 5 is in reasonable agreement with the experimental values, with a deviation of about −20%. The difference in these values is due to the rope effect behavior based on a simple relationship between timber withdrawal capacity and friction coefficient, where μ is assumed to be 0.25 or 1/4 of the maximum withdrawal capacity in timber. In these situations, where the prediction and evaluation of the stiffness and load-carrying capacity of the TCC connection are essential, as well as there is a need for a complete curve in the case of non-linear analysis of TCC structures, it is suggested to use the proposed numerical model. Additionally, in addition to the simulation of the complete curve, the numerical model can determine local effects such as the crushing depth in timber and concrete, the rotation angle of the screw, and the participation of forces in the load-carrying of the connection for any displacements in the connection.
In addition to the mentioned possibilities, it should be noted that there are certain limitations in the proposed model. They are related to the adopted assumptions about the axial force Fax(θ) in the screw approximated by a sinusoidal function, constant friction coefficient, and the influence of the group effect of the screw (through the effective number neff) on the reduction in the load-carrying capacity of the connection. Additionally, this paper presents a numerical approach that solves the problem of the ductile failure of a composite connection with two plastic hinges when a fastener is placed perpendicular to the shear plane. The possibility is left to upgrade the model in order to cover other failure modes of the connection, as well as to take into account the behavior of LWAC with lower strengths that have a brittle behavior, i.e., the effects such as cone expulsion and pull-out. Future investigation should focus on more extensive experimental and theoretical research considering these limitations to improve the predictions of the proposed numerical model and to extend it in order to cover other dowel-type fasteners and types of concrete.

Author Contributions

Conceptualization, D.M., A.R., V.V. and A.S.-Ć.; methodology, D.M., A.R. and V.V.; software, D.M. and Đ.J.; validation, D.M. and A.R.; formal analysis, D.M. and A.Č.; investigation, L.K. and A.Č.; resources, D.M., A.R. and V.V.; data curation, D.M. and A.Č.; writing—original draft preparation, D.M., A.R. and V.V.; writing—review and editing, D.M., A.R. and V.V.; visualization, D.M. and L.K.; supervision, A.R. and V.V.; project administration, A.R. and V.V.; funding acquisition, D.M., A.R. and V.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research has been supported by the Ministry of Science, Technological Development and Innovation through project no. 451-03-47/2023-01/200156, “Innovative scientific and artistic research from the FTS domain”.

Data Availability Statement

The data presented in this study are available on request.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Lowercase symbols and abbreviations
a1minimum admissible spacing of screws within one row parallel to timber grain [mm]
bcestimated initial crushing depth in concrete [mm]
btestimated initial crushing depth in timber [mm]
bckfinal crushing depth in concrete [mm]
btkfinal crushing depth in timber [mm]
dfastener diameter [mm]
deffeffective diameter of screw [mm] (deff = 1.1·din)
dininner thread root diameter of screw [mm]
doutouter thread diameter of screw [mm]
dsshank diameter of screw [mm]
faxwithdrawal strength of timber perpendicular to the grain [N/mm2]
fccompressive strength of concrete (i.e., compressive strength of concrete based on the uniaxial compressive strength of cube with dimensions 150 × 150 × 150 mm) [MPa]
fckcharacteristic compressive strength of concrete determined in accordance with EN 1992-1-1 (based on the uniaxial compressive strength of cylinder) [MPa]
fhembedment strength in timber [MPa] according to EN 1995
fh,cembedding strength of confined concrete fh,c = 4·fc [45] [MPa](fh,c = 3·fck CEN/TS 19103 [9], fh,c = fc [3], fh,c < 4·fck CEB-FIB [64], fh,c = 4–5·fc [65,66])
fh,tembedment strength in timber [MPa]
fuultimate tensile strength of the screw [MPa]
fubbending yield strength of the screw [MPa]
fub,effultimate effective strength of the screw in bending [MPa]
lcpenetration depth of screw in concrete [mm]
lc,reqrequired embedment depth of screw in concrete [mm]
ltpenetration depth of screw in timber [mm]
lt,reqrequired embedment depth of screw in timber [mm]
nnumber of fasteners in one row parallel to timber grain
neffeffective number of fasteners in one row parallel to timber grain
urecorded data of the screw displacement in the middle of the span [mm]
Uppercase case symbols and abbreviations
Fv,R,corrproposal of corrected load-carrying capacity of the TCC connection [N]
Fmax,estestimated maximum load [N] according to EN 383 and EN 26891
Faxaxial withdrawal capacity of fastener in timber [N]
Fax(θ)the tension force in the screw depending on the screw rotation angle θ [N]
Fax,faxial force in the screw Fax,f = Frope’ + Frope’’ [N]
FJohload-carrying capacity of the TCC connection per shear plane per fastener according to Johansen’s part without rope effect [N]
Fropecontribution of the rope effect to an increase in the load-carrying capacity of fasteners according to EN 1995 or combined effect of axial and friction forces (i.e., the “rope effect”) in the proposed model [N]
Frope’force that increases the indentation zone in timber when withdrawal capacity is activated [N]
Frope’’force that represents a further expansion of the indentation zone beyond the plastic hinge in timber [N]
Fv,Rload-carrying capacity of the TCC connection per shear plane per fastener [N]
Fv,eff,Rload-carrying capacity of the TCC connection per shear plane per fastener with the influence of the group effects (i.e., effective number of fasteners neff) [N]
Fy,Rembedment force on timber crushing depth btk [N]
Fμforce at the shear plane that represents the share of friction in the connection [N]
Htadditional normal force occurs at the contact when prevented extraction of the screw leads to the adhesion of timber and concrete [N]
Kbinitial slope of yield moment function Myθ [Nmm/deg]
Kffoundation modulus in timber [N/mm3]
Ksslip modulus [N/mm] according to EN 26891
Ks,0.4Fmaxslip modulus determined as a secant value at 40% of the ultimate load [N/mm]
Kserslip modulus for serviceability limit states [N/mm] according to EN 1995
Kser,TTCslip modulus for TTC in serviceability limit states [N/mm] according to EN 1995
Kser,TCCslip modulus for TCC in serviceability limit states [N/mm] according to CEN/TS 19103
L0bearing point spacing [mm]
Myyield bending moment of the screw [Nmm]
My,effeffective yield bending moment of the screw [Nmm]
Sthe assigned displacement/slip in the connection in the numerical model [mm]
Greek symbols
βratio between the embedment strength of the members in the TCC connection (fh,c/fh,t)
δslip in shear plane between timber and concrete [mm]
δcindentation in concrete [mm]
δtindentation in timber [mm]
δmaxmaximum displacement on the considered domain in the proposed model (δmax = 15 mm)
Δindentation of the screw in timber specimen in experimental test [mm]
θthe angle between the deformed axis of the screw and the timber grain direction [°]
μfriction coefficient
ρmmean density of timber [kg/m3]

References

  1. Schänzlin, J.; Dias, A.M.P.G. Design of Timber-Concrete-Composite Structures. In Proceedings of the 4th International Conference on Timber Bridges, Biel/Bienne, Switzerland, 9–12 May 2022; Volume 174, pp. 342–344. [Google Scholar] [CrossRef]
  2. Ceccotti, A. Composite Concrete-Timber Structures. Prog. Struct. Eng. Mater. 2002, 4, 264–275. [Google Scholar] [CrossRef]
  3. Dias, A.M.P.G. Mechanical Behaviour of Timber-Concrete Joints. Ph.D. Thesis, University of Coimbra, Coimbra, Portugal, 2005. [Google Scholar]
  4. Girhammar, U.A.; Pan, D.H. Exact Static Analysis of Partially Composite Beams and Beam-Columns. Int. J. Mech. Sci. 2007, 49, 239–255. [Google Scholar] [CrossRef]
  5. Stevanović, B. Eksperimentalna i Teorijska Analiza Spregnutih Nosača Drvo-Beton Izvedenih Mehaničkim Spojnim Sredstvima. Materijali i Konstrukcije 2004, 47, 29–46. Available online: https://scindeks.ceon.rs/article.aspx?artid=0543-07980402029S (accessed on 20 April 2023.).
  6. EN 1995-1; Design of Timber Structures—Part 1-1: General—Common Rules and Rules for Buildings Eurocode. European Committee for Standardization (CEN): Brussels, Belgium, 2008.
  7. EN 1995-2; Design of Timber Structures—Part 2: Bridges. European Committee for Standardization (CEN): Brussels, Belgium, 2004.
  8. EN 1994-1; Design of Composite Steel and Concrete Structures—Part 1-1: General Rules and Rules for Buildings. European Committee for Standardization (CEN): Brussels, Belgium, 2004.
  9. CEN/TS 19103 Eurocode 5; Design of Timber Structures—Structural Design of Timber-Concrete Composite Structures—Common Rules and Rules for Buildings. European Committee for Standardization (CEN): Brussels, Belgium, 2021.
  10. Jockwer, R.; Caprio, D.; Jorissen, A. Evaluation of Parameters Influencing the Load-Deformation Behaviour of Connections with Laterally Loaded Dowel-Type Fasteners. Wood Mater. Sci. Eng. 2022, 17, 6–19. [Google Scholar] [CrossRef]
  11. Symons, D.; Persaud, R.; Stanislaus, H. Slip Modulus of Inclined Screws in Timber–Concrete Floors. Proc. Inst. Civ. Eng.-Struct. Build. 2010, 163, 245–255. [Google Scholar] [CrossRef]
  12. Kozarić, L.; Kukaras, D.; Prokić, A.; Bešević, M.; Kekanović, M. Slip Modulus of Screws in Timber and Lightweight Concrete Composite Structures. BioResources 2018, 13, 6021–6032. [Google Scholar] [CrossRef]
  13. Appavuravther, E.; Vandoren, B.; Henriques, J. Behaviour of Screw Connections in Timber-Concrete Composites Using Low Strength Lightweight Concrete. Constr. Build. Mater. 2021, 286, 122973. [Google Scholar] [CrossRef]
  14. Du, H.; Hu, X.; Xie, Z.; Wang, H. Study on Shear Behavior of Inclined Cross Lag Screws for Glulam-Concrete Composite Beams. Constr. Build. Mater. 2019, 224, 132–143. [Google Scholar] [CrossRef]
  15. EN 1992-1; Design of Concrete Structures—Part 1-1: General Rules and Rules for Buildings. European Committee for Standardization (CEN): Brussels, Belgium, 2004.
  16. Zhang, Y.; Raftery, G.M.; Quenneville, P. Experimental and Analytical Investigations of a Timber–Concrete Composite Beam Using a Hardwood Interface Layer. J. Struct. Eng. 2019, 145, 04019052. [Google Scholar] [CrossRef]
  17. Khorsandnia, N.; Valipour, H.; Bradford, M. Deconstructable Timber-Concrete Composite Beams with Panelised Slabs: Finite Element Analysis. Constr. Build. Mater. 2018, 163, 798–811. [Google Scholar] [CrossRef]
  18. Tao, H.; Yang, H.; Zhang, J.; Ju, G.; Xu, J.; Shi, B. Nonlinear Finite Element Analysis on Timber-Concrete Composite Beams. J. Build. Eng. 2022, 51, 104259. [Google Scholar] [CrossRef]
  19. Liu, C.; Ferche, A.C.; Vecchio, F.J. Modelling Short-Term Monotonic Response of Timber–Concrete Composite Structures. Can. J. Civ. Eng. 2022, 49, 201–211. [Google Scholar] [CrossRef]
  20. Dias, A.M.P.G.; Van de Kuilen, J.W.G.; Lopes, S.; Cruz, H. A Non-Linear 3D FEM Model to Simulate Timber–Concrete Joints. Adv. Eng. Softw. 2007, 38, 522–530. [Google Scholar] [CrossRef]
  21. Oudjene, M.; Meghlat, E.-M.; Ait-Aider, H.; Batoz, J.-L. Non-Linear Finite Element Modelling of the Structural Behaviour of Screwed Timber-to-Concrete Composite Connections. Compos. Struct. 2013, 102, 20–28. [Google Scholar] [CrossRef]
  22. Wusqo, U.; Awaludin, A.; Setiawan, A.F.; Irawati, I.S. Study of Laminated Veneer Lumber (LVL) Sengon to Concrete Joint Using Two-Dimensional Numerical Simulation. J. Civ. Eng. Forum 2019, 5, 275. [Google Scholar] [CrossRef]
  23. De Santis, Y.; Sciomenta, M.; Spera, L.; Rinaldi, V.; Fragiacomo, M.; Bedon, C. Effect of Interlayer and Inclined Screw Arrangements on the Load-Bearing Capacity of Timber-Concrete Composite Connections. Buildings 2022, 12, 2076. [Google Scholar] [CrossRef]
  24. Dias, A.M.P.G.; Martins, A.R.D.; Simões, L.M.C.; Providência, P.M.; Andrade, A.A.M. Statistical Analysis of Timber–Concrete Connections—Mechanical Properties. Comput. Struct. 2015, 155, 67–84. [Google Scholar] [CrossRef]
  25. Töpler, J.; Buchholz, L.; Lukas, J.; Kuhlmann, U. Guidelines for a Finite Element Based Design of Timber Structures and Their Exemplary Application on Modelling of Beech LVL. Buildings 2023, 13, 393. [Google Scholar] [CrossRef]
  26. Blass, H.J.; Bienhaus, A.; Krämer, V. Effective Bending Capacity of Dowel-Type Fasteners. In PRO 22: International RILEM Symposium on Joints in Timber Structures; RILEM Publications: Stuttgart, Germany, 2001; pp. 71–80. [Google Scholar]
  27. Khorsandnia, N.; Valipour, H.; Crews, K. Structural Response of Timber-Concrete Composite Beams Predicted by Finite Element Models and Manual Calculations. Adv. Struct. Eng. 2014, 17, 1601–1621. [Google Scholar] [CrossRef]
  28. Dias, A.M.P.G. Analysis of the Nonlinear Behavior of Timber-Concrete Connections. J. Struct. Eng. 2012, 138, 1128–1137. [Google Scholar] [CrossRef]
  29. Djoubissié Denouwé, D.; Messan, A.; Fournely, E.; Bouchair, A. Influence of Interlayer in Timber-Concrete Composite Structures with Threaded Rebar as Shear Connector-Experimental Study. Am. J. Civ. Eng. Archit. 2018, 6, 38–45. [Google Scholar] [CrossRef]
  30. Zhu, W.; Yang, H.; Liu, W.; Shi, B.; Ling, Z.; Tao, H. Experimental Investigation on Innovative Connections for Timber–Concrete Composite Systems. Constr. Build. Mater. 2019, 207, 345–356. [Google Scholar] [CrossRef]
  31. Wang, Z.; Wei, Y.; Li, N.; Zhao, K.; Ding, M. Flexural Behavior of Bamboo–Concrete Composite Beams with Perforated Steel Plate Connections. J. Wood Sci. 2020, 66, 4. [Google Scholar] [CrossRef]
  32. Foschi, R.O.; Bonac, T. Load-Slip Characteristics for Connections With Common Nails. Wood Sci. 1977, 9, 118–123. [Google Scholar]
  33. Kermani, A.; Goh, H.C.C. Load-Slip Characteristics of Multi-Nailed Timber Joints. Proc. Inst. Civ. Eng.-Struct. Build. 1999, 134, 31–43. [Google Scholar] [CrossRef]
  34. Coste, G. The Assessment and Applications of a New Connector Type for Use in Timber Structural Systems. Ph.D. Thesis, Edinburgh Napier University, Edinburgh, UK, 2010. [Google Scholar]
  35. Liu, Y.; Wang, Y.; Zhang, Y.; Chen, M.; Nie, X. Force–Displacement Relations of Bolted Timber Joints with Slotted-in Steel Plates Parallel to the Grain. J. Wood Sci. 2020, 66, 83. [Google Scholar] [CrossRef]
  36. Heine, C.P.; Daniel Dolan, J. A New Model to Predict the Load-Slip Relationship of Bolted Connections in Timber. Wood Fiber Sci. 2001, 33, 534–549. [Google Scholar]
  37. Erdodi, L.; Bódi, I. Numerical Determination of the Slip Modulus of Dowel-Type Timber Joints. Pollack Period. 2007, 2, 35–44. [Google Scholar] [CrossRef]
  38. Domínguez, M.; Fueyo, J.G.; Cabezas, J.A. Design and Calculation of Dowel-type Joints in Timber Structures. Behaviour of Load against Displacement. In Proceedings of the 3rd International Conference on Mechanical Models in Structural Engineering, Seville, Spain, 24–26 June 2015; pp. 199–212. [Google Scholar]
  39. Kobel, P. Dowel-Type Connections in Beech LVL. Ph.D. Thesis, ETH Zurich, Zürich, Switzerland, 2019. [Google Scholar]
  40. DIN EN 1995-1-1/NA:2013-08; National Annex—Nationally Determined Parameters Eurocode 5: Design of Timber Structures—Part 1-1: General—Common Rules and Rules for Buildings. European Committee for Standardization (CEN): Berlin, Germany, 2013.
  41. Bjetka, I.; Blass, H.J. Joints with Inclined Screws. In Proceedings of the Working Commission W18—Timber Structure, Kyoto, Japan, 16–19 September 2002. [Google Scholar]
  42. Dias, A.M.P.G.; Cruz, H.; Lopes, S.; Van de Kuilen, J.W.G. Experimental Shear-Friction Test on Dowel Type Fastener Timber-Concrete Joints. In Proceedings of the 8th World Conference on Timber Engineering (WCTE), Lahti, Finland, 14–17 June 2004. [Google Scholar]
  43. Mirdad, M.A.H.; Chui, Y.H. Stiffness Prediction of Mass Timber Panel-Concrete (MTPC) Composite Connection with Inclined Screws and a Gap. Eng. Struct. 2020, 207, 110215. [Google Scholar] [CrossRef]
  44. Du, H.; Hu, X.; Jiang, Y.; Wei, C.; Hong, W. Load-Carrying Capacity of Self-Tapping Lag Screws for Glulam-Lightweight Concrete Composite Beams. BioResources 2019, 14, 166–179. [Google Scholar] [CrossRef]
  45. Dias, A.M.P.G.; Lopes, S.; Van de Kuilen, J.W.G.; Cruz, H. Load-Carrying Capacity of Timber–Concrete Joints with Dowel-Type Fasteners. J. Struct. Eng. 2007, 133, 720–727. [Google Scholar] [CrossRef]
  46. EN 1382; Timber Structures—Test Methods—Withdrawal Capacity of Timber Fasteners. European Committee for Standardization (CEN): Brussels, Belgium, 2016.
  47. EN 12390-3; Testing Hardened Concrete—Part3: Compressive Strength of Test Specimens. European Committee for Standardization (CEN): Brussels, Belgium, 2009.
  48. EN ISO 6892-1; Metallic Materials—Tensile Testing—Part 1: Method of Test at Room Temperature. European Committee for Standardization (CEN): Brussels, Belgium, 2009.
  49. EN 383; Timber Structures—Test Methods. Determination of the Embedding Strength and Foundation Values for Dowel Type Fasteners. European Committee for Standardization (CEN): Brussels, Belgium, 1993.
  50. ISO 10984–1; Timber Structures—Dowel Type Fasteners—Part 1: Determination of Yield Moment. International Organization for Standardization: Geneva, Switzerland, 2009.
  51. Stamatopoulos, H.; Massaro, F.M.; Qazi, J. Mechanical Properties of Laterally Loaded Threaded Rods Embedded in Softwood. Eur. J. Wood Wood Prod. 2022, 80, 169–182. [Google Scholar] [CrossRef]
  52. Aquino, C.D.; Rodrigues, L.G.; Branco, J.M.; Gomes, W.J.S. Statistical Correlation Investigation of a Single-Doweled Timber-to-Timber Joint. Eng. Struct. 2022, 269, 114810. [Google Scholar] [CrossRef]
  53. Gečys, T.; Bader, T.K.; Olsson, A.; Kajėnas, S. Influence of the Rope Effect on the Slip Curve of Laterally Loaded, Nailed and Screwed Timber-to-Timber Connections. Constr. Build. Mater. 2019, 228, 116702. [Google Scholar] [CrossRef]
  54. Blass, H.J.; Colling, F. Load-Carrying Capacity of Dowelled Connections. In Proceedings of the International Network on Timber Engineering Research, Sibenik, Croatia, 24–27 August 2015; pp. 115–129. [Google Scholar]
  55. Carvalho, E.P. Avaliação Do Comportamento Mecânico Ao Cisalhamento de Conectores de Barras de Aço Com Ganchos Em Vigas Mistas de Madeira Laminada Colada e Concreto. Ph.D. Thesis, Escola de Engenharia da Universidade Federal de Minas Gerais, Belo Horizonte, Brazil, 2008. [Google Scholar]
  56. Moshiri, F. Structural Behaviour of Timber Concrete Composite Connections and Floors Utilising Screw Connectors. Ph.D. Thesis, University of Technology Sydney, Sydney, Australia, 2014. [Google Scholar]
  57. Gelfi, P.; Giuriani, E.; Marini, A. Stud Shear Connection Design for Composite Concrete Slab and Wood Beams. J. Struct. Eng. 2002, 128, 1544–1550. [Google Scholar] [CrossRef]
  58. EN 26891; Timber Structures—Joints Made with Mechanical Fasteners—General Principles for the Determination of Strength and Deformation Characteristics. European Committee for Standardization (CEN): Brussels, Belgium, 1991.
  59. Domínguez, M.; Fueyo, J.G.; Villarino, A.; Anton, N. Structural Timber Connections with Dowel-Type Fasteners and Nut-Washer Fixings: Mechanical Characterization and Contribution to the Rope Effect. Materials 2022, 15, 242. [Google Scholar] [CrossRef] [PubMed]
  60. Stevanović, B. Ponašanje spregnutih nosača tipa drvo-beton izvedenih mehaničkim spojnim sredstvima pri eksploatacionom i graničnom opterećenju. Ph.D. Thesis, Građevinski fakultet Univerziteta u Beogradu, Beograd, Serbia, 2003. [Google Scholar]
  61. Glišovíc, I.; Stevanovíc, B.; Kočetov-Mišulić, T. Embedment Test of Wood for Dowel-Type Fasteners. Wood Res. 2012, 57, 639–650. [Google Scholar]
  62. Dias, A.M.P.G.; Cruz, H.; Lopes, S.; Van de Kuilen, J.W.G. Stiffness of Dowel-Type Fasteners in Timber-Concrete Joints. Proc. Inst. Civ. Eng. Struct. Build. 2010, 163, 257–266. [Google Scholar] [CrossRef]
  63. Mirdad, M.A.H.; Khan, R.; Chui, Y.H. Analytical Procedure for Timber-Concrete Composite (TCC) System with Mechanical Connectors. Buildings 2022, 12, 885. [Google Scholar] [CrossRef]
  64. Comite Euro-International du Beton CEB. CEB-FIP 1990. CEB FIP Model Code 1990: Design Code, CEB-FIB-MC90, CEB-FIB; Lausanne: Thomas Telford, Switzerland, 1993. [Google Scholar]
  65. Gelfi, P.; Giuriani, E. Modelo Teorico Del Legame Costitutivo per Le Connessioni a Piolo. In Studi e Ricerche Vol.9. Corso di Perfezionamento per le Costruzioni in Cemento Armato Fratelli Pesenti; Politecnico di Milano: Milan, Italy, 1987; pp. 323–341. [Google Scholar]
  66. Vintzēleou, E.N.; Tassios, T.P. Mathematical Models for Dowel Action under Monotonic and Cyclic Conditions. Mag. Concr. Res. 1986, 38, 13–22. [Google Scholar] [CrossRef]
Figure 1. The behavior of timber–concrete-bonded connection: (a) gradual development of failure and (b) Fv-δ curve on the complete displacement domain.
Figure 1. The behavior of timber–concrete-bonded connection: (a) gradual development of failure and (b) Fv-δ curve on the complete displacement domain.
Buildings 13 01171 g001
Figure 2. Prediction of the failure modes dependent on the embedment depths of a fastener in the timber and concrete.
Figure 2. Prediction of the failure modes dependent on the embedment depths of a fastener in the timber and concrete.
Buildings 13 01171 g002
Figure 3. The behavior of the connection due to the rope effect: forces and stresses in the timber–concrete connection in the shear plane.
Figure 3. The behavior of the connection due to the rope effect: forces and stresses in the timber–concrete connection in the shear plane.
Buildings 13 01171 g003
Figure 4. The EYM model approaches for failure mode III: (a) mechanics of the yield mode III of the TCC connection, (b) failure mechanism, and (c) free-body diagram of the screw.
Figure 4. The EYM model approaches for failure mode III: (a) mechanics of the yield mode III of the TCC connection, (b) failure mechanism, and (c) free-body diagram of the screw.
Buildings 13 01171 g004
Figure 5. Summary flowchart for determining the load-carrying capacity of the TCC connection.
Figure 5. Summary flowchart for determining the load-carrying capacity of the TCC connection.
Buildings 13 01171 g005
Figure 6. Determination of the embedment properties parallel to the grain: (a) the setup of the embedment test, (b) timber specimens, and (c) fh − Δ curves obtained experimentally.
Figure 6. Determination of the embedment properties parallel to the grain: (a) the setup of the embedment test, (b) timber specimens, and (c) fh − Δ curves obtained experimentally.
Buildings 13 01171 g006
Figure 7. Determination of the yield moment of screw: (a) The setup of yield moment determination, (b) a scheme and dimensions of the three-point bending tests, and (c) My − θ curves obtained experimentally.
Figure 7. Determination of the yield moment of screw: (a) The setup of yield moment determination, (b) a scheme and dimensions of the three-point bending tests, and (c) My − θ curves obtained experimentally.
Buildings 13 01171 g007
Figure 8. Details of push-out specimens: (a) dimensions of TCC specimens, (b) reinforcement of concrete (steel mesh), and (c) TCC specimen after hardening of concrete.
Figure 8. Details of push-out specimens: (a) dimensions of TCC specimens, (b) reinforcement of concrete (steel mesh), and (c) TCC specimen after hardening of concrete.
Buildings 13 01171 g008
Figure 9. Push-out test: (a) test setup and (b) adopted loading procedure for the shear tests.
Figure 9. Push-out test: (a) test setup and (b) adopted loading procedure for the shear tests.
Buildings 13 01171 g009
Figure 10. The load–slip curves from the push-out tests.
Figure 10. The load–slip curves from the push-out tests.
Buildings 13 01171 g010
Figure 11. Failure mode of the push-out specimens: (a,b) crushing deformation of the timber and bending deformation of the lag screw (for δ = 15 and 30 mm) and (c) crushing zone of concrete.
Figure 11. Failure mode of the push-out specimens: (a,b) crushing deformation of the timber and bending deformation of the lag screw (for δ = 15 and 30 mm) and (c) crushing zone of concrete.
Buildings 13 01171 g011
Figure 12. Experimental and numerical curves: (a) embedment strength in timber, (b) yield moment of the screw, and (c) load–slip curves of the connection (with 4 screws).
Figure 12. Experimental and numerical curves: (a) embedment strength in timber, (b) yield moment of the screw, and (c) load–slip curves of the connection (with 4 screws).
Buildings 13 01171 g012
Figure 13. Numerical simulation 3 (SIM(fh3z)): (a) local effects and (b) participation of forces.
Figure 13. Numerical simulation 3 (SIM(fh3z)): (a) local effects and (b) participation of forces.
Buildings 13 01171 g013
Figure 14. Experimental [60] and numerical curves: (a) embedment strength in timber, (b) yield moment of the screw, and (c) load–slip curves of the connection.
Figure 14. Experimental [60] and numerical curves: (a) embedment strength in timber, (b) yield moment of the screw, and (c) load–slip curves of the connection.
Buildings 13 01171 g014
Figure 15. Numerical simulation: (a) local effects and (b) participation of forces.
Figure 15. Numerical simulation: (a) local effects and (b) participation of forces.
Buildings 13 01171 g015
Table 1. Mechanical properties of the first set of materials.
Table 1. Mechanical properties of the first set of materials.
MaterialsPropertiesMean ValueNo. of TestsApplied Standard
Glulam GL24hAxial withdrawal capacity Fax15,525 N3EN 1382 [46]
Concrete LC45/50Cubic compressive strength fc55 MPa3EN 12390-3 [47]
Screw Φ10/150 mmUltimate tensile strength fu695 MPa5EN ISO 6892-1 [48]
Table 2. Mechanical properties of materials of TCC connection according to [60].
Table 2. Mechanical properties of materials of TCC connection according to [60].
MaterialsPropertiesMean Value
Glulam
(approx. GL24h)
Axial withdrawal capacity Fax13,090 N
Density of timber ρm476 kg/m3
Embedment strength fh,t *36.1 MPa
Foundation modulus Kf *43.4 N/mm3
Concrete MB30 (approx. C25/30)Cubic compressive strength fc46 MPa
Screw Φ10/150 mmUltimate tensile strength fu590 MPa
Yield moment My **My,eff = 0.9·fu·deff3/6 = 40,403 Nmm
Initial slope of yield moment function Kb **13,867 Nmm/deg
* fh,t and Kf were calculated according to Equations (10) and (11). ** My and Kb were estimated according to the suggestions from this paper.
Table 3. Comparison of the load-carrying capacity and slip modulus of the TCC connection.
Table 3. Comparison of the load-carrying capacity and slip modulus of the TCC connection.
MethodMean Values 1Deviations [%]Mean Values 2Deviations [%]
Load-carrying capacity
Fv,R [N]
Experimental (δ = 15 mm)13,240-10,650-
Numerical (δ = 15 mm)13,4001.211,560/12,390 f8.5/16.4
TheoreticalCEN/TS 19103 a10,491−20.88773/9400 f−17.6/−11.7
Fv,R,corr b13,218−0.211,217/12,019 f5.3/12.9
Slip modulus
K [N/mm]
Experimental2889-2580-
NumericalKs (EN 26891)2670−7.62732/2980 f5.9/15.5
Ks,0.4Fmax c2825−2.23058/3279 f18.5/27.1
TheoreticalKser,TTC d348820.7347734.8
Kser,TCC = 2·Kser,TTC e6975141.46954169.5
1 Experiment from this study: deff = 7.7 mm, ρm = 477 kg/m3, fh,t,mean = 35.3 MPa, Fax = 15,525 N, fh,c = 220 MPa, My = 46,613 Nmm, μ = 0.45. 2 Experiment from [60]: deff = 7.7 mm, ρm = 476 kg/m3, fh,t = 36.1 MPa, Fax = 13,089 N, fh,c = 184 MPa, My = 40,403 Nmm, μ = 0.45. a Fv,R = FJoh + Fax/4; b Fv,R,corr = FJoh + Frope, where FJoh = 2 · β / 1 + β · 2 · M y · f h , t · d e f f and Frope = min{100%·FJoh; μ·Fax}. c Secant value determined from the load–slip curve at 40% of the ultimate load at slip δ = 15 mm. d For timber-timber connection (TTC): Kser,TTC = ρm1.5·d/23; e For TCC connections: Kser,TCC = 2·ρm1.5·d/23 (where d = deff). f With or without effective number of fasteners neff = 0.93, Fv,eff,R/Fv,R, where Fv,eff,R = neff · Fv,R.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Manojlović, D.; Rašeta, A.; Vukobratović, V.; Čeh, A.; Kozarić, L.; Jovanović, Đ.; Starčev-Ćurčin, A. Simulation of Load–Slip Capacity of Timber–Concrete Connections with Dowel-Type Fasteners. Buildings 2023, 13, 1171. https://doi.org/10.3390/buildings13051171

AMA Style

Manojlović D, Rašeta A, Vukobratović V, Čeh A, Kozarić L, Jovanović Đ, Starčev-Ćurčin A. Simulation of Load–Slip Capacity of Timber–Concrete Connections with Dowel-Type Fasteners. Buildings. 2023; 13(5):1171. https://doi.org/10.3390/buildings13051171

Chicago/Turabian Style

Manojlović, Dragan, Andrija Rašeta, Vladimir Vukobratović, Arpad Čeh, Ljiljana Kozarić, Đorđe Jovanović, and Anka Starčev-Ćurčin. 2023. "Simulation of Load–Slip Capacity of Timber–Concrete Connections with Dowel-Type Fasteners" Buildings 13, no. 5: 1171. https://doi.org/10.3390/buildings13051171

APA Style

Manojlović, D., Rašeta, A., Vukobratović, V., Čeh, A., Kozarić, L., Jovanović, Đ., & Starčev-Ćurčin, A. (2023). Simulation of Load–Slip Capacity of Timber–Concrete Connections with Dowel-Type Fasteners. Buildings, 13(5), 1171. https://doi.org/10.3390/buildings13051171

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop