Research Article  Open Access
Jing Zhang, ZhiFang Liu, Yong Xu, MaiLi Zhang, LiuCheng Mo, "Cyclic Behavior and Modeling of Bolted Glulam Joint with Cracks Loaded Parallel to Grain", Advances in Civil Engineering, vol. 2021, Article ID 6612886, 16 pages, 2021. https://doi.org/10.1155/2021/6612886
Cyclic Behavior and Modeling of Bolted Glulam Joint with Cracks Loaded Parallel to Grain
Abstract
Under varying humidity and temperature conditions, with the constraint of metal fasteners to wood shrinkage, cracks along the bolt lines are generally observed in bolted glulam joints. A threedimensional (3D) numerical model was established in software package ANSYS to investigate the cyclic behavior of bolted glulam joints with local cracks. A reversed cyclic loading was applied in the paralleltograin direction. The accuracy of numerical simulation was proved by comparison with fullscale experimental results. Typical failure modes were reproduced in the numerical analysis with the application of wood foundation zone material model and cohesive zone material model. The effect of crack number and length on the hysteretic behavior of bolted glulam joints was quantified by a parametric study. It was found that initial cracks impair the peak capacity and elastic stiffness of bolted glulam joints significantly. More decrease in capacity was observed in joints with more cracks, and longer cracks affect elastic stiffness more dramatically. Moreover, with the existence of initial cracks, the energy dissipated and equivalent viscous damping ratio of bolted joints are reduced by 24% and 13.3%, respectively.
1. Introduction
Connections are always recognized as the weakest part of timber structures. It was reported that nearly 80% of structural collapses arise from connections [1]. The mechanical behavior of joints significantly affects resistance, durability, and energy dissipation of timber structures. Doweltype joints, an effective fastening technique extensively applied in actual engineering, have attracted tremendous research interest. Extensive experimental and numerical studies on the static [2–5] and cyclic behavior [6–8] of doweltype joints have been conducted and available in the literature. The loadcarrying capacity and failure modes of bolted joints are related to material behavior, geometry, and loading type. Relatively large energy dissipation and ductile behavior could be achieved with small diameter bolts, and slender bolts are preferred in the bolted glulam joints in order to avoid brittle failure modes [9].
In terms of modeling the mechanical behavior of bolted glulam joints, several methods have been proposed by previous researchers. An analytical method named BeamonFoundation model was first applied [10] and modified [11] to predict the monotonic strength and hysteretic behavior of nailed and bolted wood connections. To represent the dynamic response of timber connections, hysteresis models such as BoucWenBaberNoori model [12], Florence model [13], and Pinching4 [14] were employed in structural analysis. A threedimensional (3D) finiteelementbased model is a recent and efficient approach to estimate the behavior of bolted glulam joint. Complete strain and stress distributions can be obtained and observed from numerical analysis. In the numerical modeling, a reasonable definition of material behavior is particularly important to achieve relatively accurate numerical results. To take into account the localized crushing behavior of wood underneath the fasteners, a foundation zone around the fasteners with weakened material properties was defined in a twodimensional numerical model [15], and it was extended to 3D FE models by Hong [16]. To consider the progressive failure of wood in bolted glulam joints, several failure criteria were proposed by previous researchers. In the 3D numerical model of a singlebolt connection developed by Moses and Prion [17], anisotropic plasticity material model was defined for wood in compression and Weibull Weakest Link Theory was applied to consider brittle failure. It was indicated that the forcedisplacement relationship was predicted well, while the ultimate loads were underestimated. Hill yield criterion and maximum stress failure criterion were selected by Kharouf [18] when establishing a numerical model of timber joints, and it was found that failure modes and deformation predicted by numerical analysis fit well with experimental results. Hoffman failure criterion was applied by Xu et al. [19] to take into account the damage evolution of wood. TsaiWu failure criterion is able to consider the difference between tension and compression strengths of wood [20], while stress interaction coefficients defined in the criterion are difficult to determine [21]. The application of cohesive zone material model into bolted glulam joints has been attempted by several researchers [22–24]. With the definition of mixed mode damage onset and propagation criteria, it was proved to be a powerful tool to simulate splitting and row shear commonly observed in bolted joints.
Extensive researches have been conducted on mechanical behavior and strength enhancement of bolted glulam joints as above. However, little attention was paid to the mechanical behavior of bolted glulam joints with local cracks. Shrinkage cracks are generally observed in bolted glulam joints due to their sensitivity to varying relative humidity of the environment. Shrinkage and swelling strains of wood are restrained with the existence of steel fasteners in the joint area, and cracking of wood is easily caused by perpendiculartograin restrained stress [25]. Previous researches have proved that the capacity, stiffness, and ductility of bolted glulam joints are decreased by initial cracks under monotonic loading, and failure modes of bolted glulam joints switched from ductile failure with bolt yielding to brittle failure of wood [26, 27]. It is essential to further investigate the influence of initial cracks on the cyclic behavior of bolted glulam joints. In respect of numerical modeling, the complexity of contact problems and material behavior provides challenges of convergence. Even though some 3D numerical analysis has been conducted on the monotonic behavior of bolted glulam joints, few attempts have been conducted to model the cyclic behavior of timber joints by 3D numerical analysis.
To reveal the influence law of cracks on the cyclic behavior of bolted glulam joints intuitively, a 3D numerical model was established in this paper. Hill yield criterion and wood foundation model were applied to simulate the behavior of wood in compression and wood around fasteners, respectively. Brittle failure of wood and propagation of initial cracks were modeled with the application of cohesive zone material model. Several contact pairs were defined to model the interaction between different parts in the joint, and reversed cyclic loading was applied in the paralleltograin direction. The experimentally observed failure modes and hysteretic curves were used to verify the numerical results, and the influence of different crack patterns was explored by a parametric study.
2. Development of FiniteElementBased Model
2.1. Material Modeling
To model the mechanical behavior of wood in compression, the transversely isotropic plastic material model was defined with the application of the Hill yield criterion, which had been incorporated in business software ANSYS. To take into account the embedment behavior of wood underneath bolts, a foundation zone was set around bolts with a radius of 1.8time bolt diameter. The generalized Hill plasticity model with anisotropic hardening rule was applied to represent the mechanical behavior of wood foundation. 1% of initial modulus was taken as tangent modulus of hardening segment. With embedment tests loaded in perpendicular and paralleltograin directions, a bilinear relationship between load per unit length and deformation was obtained, as shown in Figure 1, and nominal yield stress and strain of foundation zone were calculated by
(a)
(b)
where and are the nominal yield strain and stress (MPa) of wood foundation, respectively; and are the deformation (mm) and yield load (N/mm) in Figure 1. Further calibration was conducted to derive effective foundation parameters. Shear modulus and yield strain were calculated according to a theoretical formula given in [28], as presented in
where is the initial shear modulus; is the effective modulus of wood foundation; is Poisson’s ratio; is the shear yield strain; and is the intercept of the second linear portion in the bilinear stressstrain relationship. For bolts and steel plate existing in bolted glulam joints, elastic modulus and Poisson’s ratio were taken as 210,000 MPa and 0.3, respectively. The isotropic plastic material model with a bilinear constitutive relationship was defined, as shown in Figure 2.
To simulate the propagation of initial cracks and brittle failure (i.e., splitting and row shear) of wood in bolted glulam joints, surfacetosurface contact pairs with elements CONTA174 and TARGE170 were defined on the predicted cracking path. To characterize the constitutive relationship of the interface, the cohesive zone material model with damage initiation and growth criteria was implemented in the numerical model. As shown in Figure 3, a bilinear traction separation law was defined, and mode I and II fractures caused by perpendiculartograin tensile stress and longitudinal shear stress were taken into account. In Figure 3, and represent normal and tangential contact stiffness, respectively. A linear softening segment was observed after peak point corresponding to material strengths and , and complete debonding occurs when reaching critical values of relative displacement or . Mixed mode fracture was considered in the analysis with a combined energy criterion, expressed as follows:
(a)
(b)
where and indicate fracture energies in modes I and II and and represent respective critical values, which are obtained from Double Cantilever Beam (DCB) test and End Notched Flexure (ENF) test, respectively, with the application of compliance combination method. In this numerical model, a small fictitious viscosity is introduced to avoid convergence difficulties.
2.2. FiniteElement Model
A 3D finiteelement model of bolted glulam joints was established in FEsoftware ANSYS. Eightnode solid element SOLID185 was applied to embody the joints. To model the interaction between bolts and steel plate, timber and bolts, timber and steel plate, and timber and steel gasket, surfacetosurface contact pairs were defined with CONTA174 and TARGE170 elements. For bolttosteel plate and woodtosteel contacts, coefficients of friction were taken as 0.001 and 0.3, respectively [29]. An example of the distribution zone of different material models was presented in Figure 4. The transversely isotropic plastic material model with the Hill yield criterion was induced to represent the compressive behavior of wood in the joint. Wood foundation zone material model was applied to simulate the embedment behavior of wooden parts surrounding the bolts, and the radius of the cylinder was set to be , where indicates the bolt diameter. The mechanical behavior of bolts and steel plate was modeled with a bilinear isotropic plastic constitutive relationship. To simulate the splitting and plug shear failure of wood in the joint area, crack growth paths along the bolt lines were predicted based on experimental observations [23], as shown in Figure 5. The width of wooden parts plugged out by bolt load was assumed to be , and was determined based on the friction between bolt and wood as given in [30]. For bolted glulam joints with initial cracks, hard contact conditions were defined to simulate the opening and closure of initial cracks. All the material constants applied in the analysis were listed in Table 1.
(a)
(b)
(a)
(b)

To study the cyclic behavior of bolted glulam joints, a gap of 30 mm was set between the top surface of steel plates and the slot, as presented in the side view of Figure 4. The bottom surface of the slottedin steel plate was fixed, and a reversed cyclic loading was applied on the top end of glulam specimens. The loading protocol recommended in EN 12512 standard [31] was applied, as presented in Figure 6. The parameter required for cyclic loading protocol was determined as yield displacement under monotonic loading. The amplitudes of loading cycles ranged from to , the single cycle was considered for the amplitudes and , and three cycles were repeated for other amplitude levels.
3. Model Verification
To validate the feasibility of the numerical model, a bolted glulam joint with five replicates was manufactured and tested. The configuration of tested joints is presented in Figure 7. The width and thickness of glulam members are 200 mm and 100 mm, respectively. The slottedin steel plates are 8 mm in thickness with a nominal yield stress of 235 MPa. Two bolts with a strength grade of 8.8 and a diameter of 12 mm are included in the joint, which has a nominal yield strength of 640 MPa according to Chinese code [32]. As can be seen from Figure 8, under the cyclic loading, obvious bending deformation of bolts was observed in the bolted glulam joints.
The numerical model of the tested joint was established and the results obtained from numerical analysis were compared to experimental findings for verification purposes. The failure modes of bolted joints predicted by numerical analysis are shown in Figure 8. Bolt bending deformation and embedment deformation of wood observed in the joint are quite similar to the experimental phenomenon. The comparison between the test obtained and the numerical predicted hysteretic curves was conducted and presented in Figure 9. As can be seen from the figure, the shape of the forcedisplacement hysteretic loops predicted by numerical analysis fit well with the test results. The deviation in displacement levels is caused by the difference between experimental obtained and numerical predicted yield displacement under monotonic loading, which was related to the variability in wood properties and bolt clearance existing in test specimens. The mechanical parameters such as elastic stiffness and peak loading capacity are calculated and listed in Table 2. It can be seen that elastic stiffness obtained from the test is lower than the numerical results due to wood properties variation in the tested specimens. The experimentally obtained peak loadbearing capacity is slightly larger than the numerically predicted values, which is believed to relate to the discrepancy between the nominal and actual yield strength of bolts. The enclosed area of hysteretic loops was calculated to obtain the accumulative energy dissipation of the joints under reversed cyclic loading. As listed in Table 2, the numerical predicted energy dissipation is 1766.4 J, which is only 7% larger than the test value. It validates the reasonability of applying the numerical model to simulate the cyclic behavior of bolted glulam joints.

4. Parametric Study
To investigate the effects of crack patterns on the cyclic behavior of bolted glulam joints, a parametric study was conducted with the consideration of different bolt configurations. A singlebolt line with two bolts was included in Joint 1, as shown in Figure 10(a). The crosssectional size of the glulam member is 120 x 90 mm, and the bolts are 12 mm in diameter with a nominal yield strength of 480 MPa. The geometric size of the glulam member in Joint 2 is identical to the tested joint, and the bolts are 16 mm in diameter with a nominal yield strength of 640 MPa, as presented in Figure 10(b). Different crack numbers and locations were considered in the analysis. The detailed information of crack patterns is shown in Figure 11. For crack patterns SB and ST in Joint 1, different crack lengths are considered and initial crack is cut through the bottom and top bolt, respectively. Initial cracks set on both sides of the steel plate are included in Joint 1 with crack patterns DB and DT. For crack patterns LS and LD in Joint 2, initial cracks along the left bolt line are set on single and both sides of the slottedin steel plate, respectively. Initial cracks along both right and left bolt lines are included in the joints with crack patterns RLS and RLD.
(a)
(b)
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
4.1. Failure Modes and Hysteretic Curves
For initially perfect Joint 1, splitting failure of wood is observed under positive loading (i.e., tension load), as shown in Figure 12(a). Bolt yielding occurred under negative loading, as can be seen from the von Mises stress distribution of the middle cross section in Figure 12(b). In terms of the hysteretic curve, an obvious pinch phenomenon is observed from Figure 13(a), which is caused by the brittle failure of wood and plastic deformation of bolts. Under positive loading, a sudden drop of loadcarrying capacity occurred at a loading amplitude level of 2.0 mm due to the splitting failure of wood. Compared to the mechanical behavior of the joint under positive loading, a more ductile performance is observed under negative loading. The joint with crack pattern DT exhibited a more pitched hysteretic behavior when compared with the initially perfect joint, as shown in Figure 13(b). The peak loading capacity and elastic stiffness are reduced by initial cracks. Under positive loading, a small decrease in capacity occurred when the displacement reaches 2 mm, which is caused by the propagation of initial cracks.
(a)
(b)
(a)
(b)
(c)
(d)
To investigate the effect of different numbers of cracks on the hysteretic behavior of bolted glulam joints, a comparison is conducted between Joint 1 with crack patterns SB and DB as presented in Figure 13(c). As can be seen from the figure, for Joint 1SB with a single crack on one side of the steel plate, the brittle failure of wooden parts on the other side leads to the decrease of the loadcarrying capacity from 30.7 kN to 19.2 kN. Compared to Joint 1SB, a relatively lower peak capacity on the positive side is observed in Joint 1DB with two initial cracks. Different crack lengths are considered in Figure 13(d), and the hysteretic curve of Joint 1DB is compared to that of 1DT. As can be seen from the figure, crack lengths only affect the elastic stiffness of bolted joints under positive loading. The initial crack in Joint 1DB propagates from bottom bolt to top bolt at the displacement of 1 mm, and a slight decrease of capacity is observed in the figure. After this crack propagation, the same cracking conditions are obtained in these two joints, and thus the hysteretic curve of Joint 1DB coincides with Joint 1DT under further application of displacement load.
The main failure modes of initially perfect Joint 2 include the embedment failure of wooden parts around bolt holes and brittle failure of wood in the joint area, as shown in Figure 14. Obvious embedment deformation is observed before the occurrence of splitting failure of wood, as presented in Figure 14(a), while the bolts remain stiff and straight, and no significant bolt bending deformation is observed. The von Mises stress distribution of the middle section is given in Figure 14(c). It is indicated that only a slight portion of the bolts reaches the yield strength and the rest performs an elastic behavior.
(a)
(b)
(c)
Compared to Joint 1, a more significant pinch phenomenon is observed in the hysteretic curve of Joint 2, as shown in Figure 15(a), which is caused by the unrecoverable embedment deformation and splitting failure of wood. When the relative displacement of the initially perfect joint reaches 3.1 mm, splitting failure occurred along the bolt lines, and the capacity of the joint under positive loading decreases from 71 kN to 46 kN. In the loading cycle with an amplitude level of 6.0 mm, the peak capacity of the joint under positive loading is close to 60 kN. The capacity drops quickly in the unloading process, and the residual plastic deformation is around 4.3 mm. Nearly zero capacity is observed when the displacement declines from 4.3 mm to −2.5 mm under reversed loading due to the significant embedment deformation. By comparison with the initially perfect joint, the peak loading capacity and elastic stiffness of the joint with crack pattern RLD are impaired by initial cracks, as can be observed from the positive side of Figure 15(a). The propagation of cracks occurred at the displacement of 2.9 mm. In Figure 15(b), a comparison is conducted between Joint 2 with crack patterns LD and RLD to explore the influence of different numbers of cracks. It can be seen that when compared with Joint 2RLD including four cracks, a larger peak capacity is observed in Joint 2LD with two cracks along the left bolt line. A drop of capacity in Joint 2LD is observed at the displacement of 4.2 mm due to the splitting failure of wood along the right bolt line.
(a)
(b)
4.2. Skeleton Curves
The skeleton curves of Joint 1 are presented in Figure 16, which are obtained by connecting peak points of each primary loading cycle in the forcedisplacement hysteretic loops. Different numbers and lengths of initial cracks are considered in the comparison. As can be seen from Figure 16, the positive and negative sides of skeleton curves can be divided into three segments, respectively. A linear loaddisplacement relationship is observed in the first segment. Under positive loading, the stiffness of the linear segment is reduced by initial cracks when compared with the initially perfect joint, which is caused by the influence of crack opening on the interaction between bolts and wooden parts. The second segment is nonlinear, and stiffness attenuation is observed due to the occurrence of bolt yielding. After the peak point, a significant decrease in loadbearing capacity is observed from skeleton curves, which is caused by a brittle failure of wood in the joint area.
Mechanical parameters of Joint 1 under positive and negative loading were calculated and listed in Table 3, including elastic stiffness, yield force, yield displacement, peak capacity, ultimate displacement, and ductility ratio. The yield point was estimated based on ASTM D5764 [33], as shown in Figure 17. A straight line is drawn parallel to the linear segment and the offset is set as 5% of bolt diameter. The intersection point between the line and the skeleton curve is determined as the yield point. For envelope curves with the declining branch, the ultimate displacement is defined as the deformation corresponding to 80% of peak strength and the ductility ratio is calculated as ultimate displacement divided by yield displacement.
 
Note. Mechanical parameters are calculated from the positive and negative sides of skeleton curves, respectively. 
As can be seen from Figure 16 and Table 3, the peak loading capacity and initial stiffness of Joint 1 under positive loading are impaired by initial cracks, while negligible influence is observed on the mechanical performance under negative loading. Under positive loading, the peak capacity of Joint 1SB and Joint 1DB is 31.4 kN and 29.4 kN, respectively. Compared to the initially perfect joint, the decreasing ratio of capacity can be up to 21.8%, which is calculated by dividing the decrease by the corresponding value of the initially perfect joint. Further, it is indicated that more cracks lead to a more significant reduction in peak loading capacity. In terms of elastic stiffness, when compared to Joint 1SB, lower stiffness is observed in Joint 1DB with two initial cracks.
Different crack lengths are considered in Joint 1DB and Joint 1DT. As can be seen from Table 3, a similar peak capacity is found for these two joints, which indicated that crack lengths have little influence on the peak capacity of bolted joints. The elastic stiffness of Joint 1DB and Joint 1DT is 24.5 kN/mm and 17.2 kN/mm, respectively. Under positive loading, a longer crack reduces the elastic stiffness of bolted joint more significantly. When compared with an initially perfect joint, the decreasing ratio of elastic stiffness can be up to 47.4% with the existence of two cracks through the top bolt. In terms of ductility ratio, it can be found from Table 3 that the ductility ratio of the initially perfect joint is lower than the joints with initial cracks. Brittle failure of wood under external load leads to a sudden drop of capacity in the initially perfect joint, while a smoothing skeleton curve is realized in Joint 1DB and Joint 1DT with the occurrence of crack propagation instead of brittle failure of wood. The ductility ratio of Joint 1DB and Joint 1DT is 2.85 and 2.32, respectively.
The skeleton curves and mechanical parameters of Joint 2 are presented in Figure 18 and Table 4, respectively. Similar to Joint 1, the reduction in the peak loading capacity and elastic stiffness is observed with the existence of initial cracks. As can be seen from the table, the peak capacity of Joint 1LS and Joint 1RLD under positive loading is 68.4 kN and 60.1 kN, respectively. Compared to an initially perfect joint with a capacity of 71.4 kN, the decreasing ratio of capacity can be up to 15.8%. With initial cracks on both bolt lines, the elastic stiffness of bolted joints declines from 49.6 kN/mm to 40.2 kN/mm with a decreasing ratio of 19.0%. Compared to Joint 1, the influence of initial cracks is less significant on the peak capacity and elastic stiffness of Joint 2, while initial cracks impair the ductility behavior of Joint 2, contrary to the observed trend in Joint 1. Similar ultimate displacement is observed in the joints with or without initial cracks, as shown in Table 4.
 
Note. Mechanical parameters are calculated from the positive and negative sides of skeleton curves, respectively. 
4.3. Stiffness Degradation
Under reversed cyclic loading, the stiffness of bolted glulam joints declined gradually with increased loading amplitudes with wood crushing or bolting yielding. The secant stiffness of primary loading cycles was calculated as follows:
where and represent the maximum positive load and corresponding displacement of loading cycle ; and indicate the maximum negative force and corresponding displacement of cycle . The change of secant stiffness with the increase of displacement is presented in Figure 19. It can be seen that the secant stiffness of the initially perfect joint is the largest in the same loading cycle compared to joints with initial cracks. Also, more numbers of cracks and longer cracks lead to more decrease in the secant stiffness of joints.
(a)
(b)
At the early loading stage, the gradual growth of secant stiffness is observed with full contact interaction. For initially perfect Joint 1, the largest secant stiffness of 36.5 kN/mm is obtained at the displacement of 0.5 mm, as shown in Figure 19(a). It decreases to 32 kN/mm when the amplitude of the loading cycle rises to 1.0 mm, which is caused by the development of bolt yielding. A deep decrease of secant stiffness is observed with the amplitude level growing to 2 mm due to the occurrence of brittle failure of wood. For Joint 1 with crack pattern DT, the peak secant stiffness is around 26 kN/mm, which is 28.8% lower than the initially perfect joint. A relatively flat segment was observed after the peak point, and a little decrease of secant stiffness occurs when the amplitude increases from 0.5 mm to 1.0 mm. It is indicated that the damage accumulation such as bolt yielding and wood locally compression is not significant at the loading stage.
Similar to Joint 1, a first growing tendency of secant stiffness is observed in Joint 2, as shown in Figure 19(b). The largest secant stiffness of the initially perfect joint is 50.9 kN/mm obtained at the displacement of around 1 mm. The peak value of Joint 2RLD is 45.9 kN/mm, which is decreased by 9.8% with the existence of initial cracks. A rapid drop of secant stiffness is observed after the peak point due to the occurrence of embedment failure of wood. Compared to Figure 19(a), it can be found that the decrease of peak secant stiffness in Joint 2 is much lower than that observed in Joint 1, which indicates that initial cracks have a less dramatic effect on the secant stiffness of bolted joints with embedment failure mode.
4.4. Energy Dissipation
To estimate the accumulative energy dissipation of the joint under reversed cyclic loading, the enclosed area of the hysteretic loop was calculated for each loading cycle based on numerical integration, and the accumulative energy dissipation was presented in Figure 20. As can be seen from Figure 20(a), the energy dissipated in Joint 1 is negligible during the loading cycles with a magnitude lower than 0.5 mm, which indicated that the materials behave elastically at this deformation level. With the displacement increasing from 1.0 mm to 4.0 mm, bolt yielding and obvious bolt bending deformation are observed in the joint, and the energy dissipated grows rapidly. The total energy dissipated in initially perfect Joint 1 is around 677 J. Compared to joints with crack pattern DB and DT, the joints with single crack dissipated more energy at the displacement lower than 4 mm, while similar total energy dissipated of around 514 J is obtained at the end of the loading process. The energy dissipated of bolted glulam joints with initial cracks is only 76% of the initially perfect joint.
(a)
(b)
For initially perfect Joint 2, a similar increasing tendency is observed and the total energy dissipated is 932 J. As can be found from Figure 20(b), lower energy dissipated is observed in the joints with more cracks. The total energy dissipated of Joint 2RLD is 764 J, and the decreasing ratio is 18% compared to the initially perfect joint.
Equivalent viscous damping ratios have been calculated at different amplitude levels as follows:where represents the enclosed area of the hysteretic loop for each primary loading cycle and equals the area of two triangles (), as shown in Figure 21. The relationships between equivalent viscous damping ratios and the displacement are presented in Figure 22.
(a)
(b)
It can be seen from Figure 22(a) that the equivalent viscous damping ratios of Joint 1 are quite low at displacement of 0.5 mm. With the displacement rising from 0.5 mm to 2.0 mm, the damping ratio of the initially perfect joint increases to 0.18 linearly with damage accumulation, and a damping ratio of 0.196 is achieved at the displacement of 3.0 mm. For Joint 1 with initial cracks, the damping ratios increase as the displacement grows from 0.5 mm to 2.0 mm, and a damping ratio of around 0.17 is obtained. The damping ratio is reduced by 13.3% with the existence of initial cracks. After the peak point, a drop of equivalent viscous damping ratio is observed with further application of displacement loading.
As can be seen from Figure 22(b), a similar changing tendency of equivalent viscous damping ratio is observed in Joint 2 before the displacement reaches 2.0 mm. The equivalent viscous damping ratio of initially perfect Joint 2 is lower than that of Joint 1. A peak damping ratio of 0.114 is obtained at the displacement of 2.0 mm and the damping ratio decreases with further increase in loading amplitude. The damping ratio of Joint 2LS and Joint 2RLD is 0.112 and 0.109, respectively. Less significant influence of initial cracks on equivalent viscous damping ratio is observed in Joint 2 when compared to Joint 1.
5. Conclusion
In the paper, a 3D numerical model was developed to investigate the influence of initial cracks on the cyclic behavior of bolted glulam joints under paralleltograin loading. Cohesive zone material law was applied to simulate the propagation of initial cracks and brittle failure of wood. With the application of the Hill yield criterion and wood foundation zone model, the local crushing behavior of wood was reproduced by numerical results. The feasibility of the numerical model was verified by comparison with fullscale experimental results, and different crack patterns and bolt configurations were further considered in a parametric study.
It was found that peak capacity and elastic stiffness of bolted glulam joints were reduced with the existence of initial cracks. More decrease in capacity was observed in joints with more cracks, and longer cracks affect elastic stiffness more dramatically. For Joint 1 with splitting failure occurring under positive loading, the decreasing ratios can be up to 21.8% and 47.4%, respectively. The ductility behavior of Joint 2 with controlled failure mode of wood embedment is impaired by initial cracks. Moreover, under reversed cyclic loading, less energy dissipated is observed in bolted glulam joints with initial cracks, which is only 76% of energy dissipated in the initially perfect joint. Further, the equivalent viscous damping ratio of bolted glulam joint is reduced by 13.3% with the existence of initial cracks.
Data Availability
The data used to support the findings of this study are included within the article and available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors gratefully acknowledge the National Natural Science Foundation of China (Grant no. 51908147) for supporting this research.
References
 R. Y. Itani and K. F. Faherty, Structural Wood Research: StateofTheArt and Research Needs, ASCE, Reston, VA, USA, 1984.
 T. V. P. Caldeira, N. Dourado, A. M. P. De Jesus, M. F. S. F. De Moura, and J. J. L. Morais, “Quasistatic behavior of momentcarrying steelwood doweled joints,” Construction and Building Materials, vol. 53, pp. 439–447, 2014. View at: Publisher Site  Google Scholar
 E.M. Meghlat, M. Oudjene, H. AitAider, and J.L. Batoz, “A new approach to model nailed and screwed timber joints using the finite element method,” Construction and Building Materials, vol. 41, pp. 263–269, 2013. View at: Publisher Site  Google Scholar
 C. L. Santos, A. M. P. De Jesus, J. J. L. Morais, and B. F. C. Fontoura, “An experimental comparison of strengthening solutions for doweltype wood connections,” Construction and Building Materials, vol. 46, pp. 114–127, 2013. View at: Publisher Site  Google Scholar
 P. Zarnani and P. Quenneville, “Strength of timber connections under potential failure modes: an improved design procedure,” Construction and Building Materials, vol. 60, no. 16, pp. 81–90, 2014. View at: Publisher Site  Google Scholar
 I. Gavric, M. Fragiacomo, and A. Ceccotti, “Cyclic behaviour of typical metal connectors for crosslaminated (CLT) structures,” Materials and Structures, vol. 48, no. 6, pp. 1841–1857, 2015. View at: Publisher Site  Google Scholar
 A. Jorissen and M. Fragiacomo, “General notes on ductility in timber structures,” Engineering Structures, vol. 33, no. 11, pp. 2987–2997, 2011. View at: Publisher Site  Google Scholar
 A. Lokaj and K. Klajmonová, “Dowel type joints of round timber exposed to static and cyclic tension forces,” Procedia Engineering, vol. 114, pp. 240–247, 2015. View at: Publisher Site  Google Scholar
 M. Popovski, “Seismic performance of braced timber frames,” University of British Columbia, Vancouver, Canada, 2000, Doctoral dissertation. View at: Google Scholar
 T. L. Wilkinson, “Analysis of nailed joints with dissimilar members,” Journal of the Structural Division, vol. 98, no. 9, pp. 2005–2013, 1972. View at: Publisher Site  Google Scholar
 R. O. Foschi and H. Prion, “Behaviour of heavy timber connections under cyclic loads,” Report, University of British Columbia, Vancouver, B.C., Canada, 1996, Thesis. View at: Google Scholar
 J. Xu and J. D. Dolan, “Development of nailed wood joint element in ABAQUS,” Journal of Structural Engineering, vol. 135, no. 8, pp. 968–976, 2009. View at: Publisher Site  Google Scholar
 Y.L. Shen, J. Schneider, S. Tesfamariam, S. F. Stiemer, and Z.G. Mu, “Hysteresis behavior of bracket connection in crosslaminatedtimber shear walls,” Construction and Building Materials, vol. 48, pp. 980–991, 2013. View at: Publisher Site  Google Scholar
 Z. Shu, Z. Li, M. He, X. Zheng, and T. Wu, “Seismic design and performance evaluation of selfcentering timber moment resisting frames,” Soil Dynamics and Earthquake Engineering, vol. 119, pp. 346–357, 2019. View at: Publisher Site  Google Scholar
 R. O. Foschi, “Loadslip characteristic of nails,” Wood Science, vol. 7, pp. 69–76, 1974. View at: Google Scholar
 J. P. Hong, “Threedimensional nonlinear finite element model for singe and multiple doweltype wood connections,” University of British Columbia, Vancouver, B.C., Canada, 2007, Doctoral Dissertation. View at: Google Scholar
 D. M. Moses and H. G. L. Prion, “Stress and failure analysis of wood composites: a new model,” Composites Part B: Engineering, vol. 35, no. 3, pp. 251–261, 2004. View at: Publisher Site  Google Scholar
 N. Kharoufa, G. McClure, and I. Smith, “Elastoplastic modelling of wood bolted connections,” Computers and Structures, vol. 81, no. 811, pp. 747–754, 2003. View at: Publisher Site  Google Scholar
 B. H. Xu, A. Bouchaïr, M. Taazount, and E. J. Vega, “Numerical and experimental analyses of multipledowel steeltotimber joints in tension perpendicular to grain,” Engineering Structures, vol. 31, no. 10, pp. 2357–2367, 2009. View at: Publisher Site  Google Scholar
 M. PattonMallory, P. Pellicane, and F. Smith, “Qualitative assessment of failure in bolted connections: tsaiWu criterion,” Journal of Testing and Evaluation, vol. 26, no. 5, pp. 497–505, 1998. View at: Google Scholar
 M. PattonMallory, P. J. Pellicane, and F. W. Smith, “Modeling bolted connections in wood: review,” Journal of Structural Engineering, vol. 123, no. 8, pp. 1054–1062, 1997. View at: Publisher Site  Google Scholar
 N. Dourado, F. G. A. Silva, and M. F. S. F. De Moura, “Fracture behavior of woodsteel dowel joints under quasistatic loading,” Construction and Building Materials, vol. 176, pp. 14–23, 2018. View at: Publisher Site  Google Scholar
 J. Zhang, M.J. He, and Z. Li, “Numerical analysis on tensile performance of bolted glulam joints with initial local cracks,” Journal of Wood Science, vol. 64, no. 4, pp. 364–376, 2018. View at: Publisher Site  Google Scholar
 C. L. Santos, A. M. P. De Jesus, J. J. L. Morais, and J. L. P. C. Lousada, “Quasistatic mechanical behaviour of a doubleshear single dowel wood connection,” Construction and Building Materials, vol. 23, no. 1, pp. 171–182, 2009. View at: Publisher Site  Google Scholar
 J. Sjödin and C.J. Johansson, “Influence of initial moisture induced stresses in multiple steeltotimber dowel joints,” Holz als Roh und Werkstoff, vol. 65, no. 1, pp. 71–77, 2006. View at: Publisher Site  Google Scholar
 M. He, J. Zhang, and Z. Li, “Influence of cracks on the mechanical performance of dowel type glulam bolted joints,” Construction and Building Materials, vol. 153, pp. 445–458, 2017. View at: Publisher Site  Google Scholar
 J. Zhang, M. He, and Z. Li, “Mechanical performance assessment of bolted glulam joints with local cracks,” Journal of Materials in Civil Engineering, vol. 30, no. 6, 2018. View at: Publisher Site  Google Scholar
 E. P. Saliklis, T. J. Urbanik, and B. Tokyay, “Bilinear modelling of cellulosic orthotropic nonlinear materials,” Journal of Pulp and Paper Science, vol. 29, no. 12, pp. 407–411, 2003. View at: Google Scholar
 B. H. Xu, M. Taazount, A. Bouchaïr, and P. Racher, “Numerical 3D finite element modelling and experimental tests for doweltype timber joints,” Construction and Building Materials, vol. 23, no. 9, pp. 3043–3052, 2009. View at: Publisher Site  Google Scholar
 J. L. Jensen and P. Quenneville, “Fracture mechanics analysis of row shear failure in dowelled timber connections,” Wood Science and Technology, vol. 44, no. 4, pp. 639–653, 2009. View at: Publisher Site  Google Scholar
 BSI, Timber StructuresTest MethodsCyclic Testing of Joints Made with Mechanical Fasteners, BSI, London, UK, 2001.
 BSI, Mechanical Properties of Fasteners Bolts, Screws and Studs, BSI, London, UK, 2000.
 ASTM, Standard Test Method for Evaluating DowelBearing Strength of Wood and Woodbase Products, ASTM, West Conshohocken, PA, USA, 2013.
Copyright
Copyright © 2021 Jing Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.