Research Article  Open Access
Huazhang Zhou, Xiaoqiang Zhou, "FE Modeling for Bolted Wood Connection Using a Porous Constitutive Model", Advances in Civil Engineering, vol. 2020, Article ID 6621333, 8 pages, 2020. https://doi.org/10.1155/2020/6621333
FE Modeling for Bolted Wood Connection Using a Porous Constitutive Model
Abstract
According to the facts of localized crushing failure of bolt groove in wood connection with enough end distance and the threephase composites of wood with solid (wood substance), water, and gas, a confined compression test for the wood cylinder was conducted for achieving constitutive relation under the complex stress state in wood groove. A porous constitutive model was developed according to the confined compression experiments. Then, the constitutive model was implemented in a finite element modeling of mental doweltype fasters in woodtowood connections to analyse the loadcarrying capacity parallel to the grain. Through changing the thicknesses of centre members and side members of wood connections made of a similar wood species, Pinus Sylvestris var. Mongolica, the effects of thickness combinations of centre members and side members on the failure modes and loadcarrying capacity of bolted wood connection including numerical simulations and experiments were compared. The failure modes, including the yielding of centre member, the yielding of side member, and the yielding of the bolt, as well as the rigid rotation of the bolt, all reappeared by the finite element modeling with the porous constitutive model. The predicted deformation shapes and loaddisplacement relations of bolted wood connections were compared with experimental ones, and good correlations were observed. This paper presents a new approach to simulate the local embedment crushing of bolt groove in wood connections.
1. Introduction
Finite element (FE) numerical simulation of the localized problem in wood structure is always a difficult issue due to the inevitable differences between FE results and experiment results. The reasons causing the differences involve the high variable mechanical properties of wood, differences between clear wood and structural wood, size effect, defects of wood, effect of moisture content, duration of load, and even the test method for material parameters for the localized compression wood. Hence, to reduce the differences between FE simulation and experiment, some hypotheses were adopted in FE modeling to calibrate the mechanical parameters. For example, Moses reduced the initial modulus to cover the workhardening behaviour of wood, and a multiplier 0.01 was used in the modeling [1]. Guan adopted a modification factor of the empirical modulus to simulate the stiffness of the flat nail embedment test, where a multiplier 0.25 was used [2]. Barrett developed a FE modeling using a foundation material model and limit foundation zone around bolt hole of 4.5 times of diameter of the bolt to simulate the localized problem of nailed connections in wood [3]. Sandhaas developed a continuum damage mechanics model including eight types of failure modes to simulate the performance of timber joints with the slottedin steel plate, where the size effect for fracture energy was taken into consideration by means of the characteristic element length [4]. By utilizing the modification factors, the simulation results were consistent with those of corresponding experiments. However, the modification factors are always empirical coefficients which are dependent on the experiences of sophisticated engineers.
Difference induced by the test method for mechanical parameter is also one source for the differences between simulation results and experiment results. The mechanical parameters for the constitutive model in FE modeling were usually determined through uniaxial compression or tension tests. However, in bolted wood connection, the stress state in wood groove is not uniaxial one. Localized wood groove is restrained by lateral wood, which usually deforms along the direction parallel to the grain including shear deformation and crushing deformation components. However, a stressstrain relation from the uniaxial test cannot descript accurately the stress state in wood groove. On the other hand, the dowelbearing strength test specified in ASTM D576497a only specifies the average embedment strength of groove [5]. The embedment strength is also different with the strength from the uniaxial test due to different stress states. Hence, a confined compression test of the wood cylinder was conducted to obtain the real deformationcompression relation in wood groove. A new approach to simulate the localized issue of wood groove in bolted connection by applying the porous constitutive model for wood groove was established to mitigate the dependency on the experiences.
A threedimensional (3D) FE modeling developed in this paper is based on the phenomena of failure process and failure model of wood groove in bolted connection. During the contacting between bolt and wood groove, localized wood groove is compressed and wood void becomes less and less. If without enough end distance (7d) for the connection, wood groove will crack along the direction of load and shear failure will happen. However, for the connection with enough end distance, only localized wood will be crushed and then plastic hinge may appear in the bolt. In the process of groove crushing, the embedment stiffness of connection changes with the changing of the contact state. At the beginning of contact, embedment stiffness is small. Then, it becomes bigger and nearly keeps a constant. When groove crush is happening, embedment stiffness increases steeply which means that groove becomes compact and the wood void becomes small. In the numerical modeling, the effect of embedment stiffness and strength will be reflected in the porous constitutive by the corresponding relationship between void ratio of wood and pressure on the porous medium.
2. Porous Constitutive Model for Wood
Wood can be seemed as a threephase medium which is composed of wood solid, liquid, and gas. The wood substance components of wood are the product of wood parenchyma. The liquid component is usually water, and the gas component is usually air. The gaps between the wood substance are called voids. The voids may contain air, water, or both. The void ratio e is the proportion of the volume of voids with respect to the volume of wood substance [6]. Volume of wood changes due to the water expelled from the voids. Mechanical behaviour of wood is also largely dependent on the proportion of wood solids versus voids. When wood is localized compressed, crushing failure is one of typical failure modes and its void ratio becomes less and less. In the elastic stage, the relationship between void ratio and compression stress can be described by a differential equation, the change of the void ratio and the change of logarithm of the equivalent pressure stress [7], shown in the following equation:where is the void ratio of wood; is the equivalent pressure stress; , , and are stress components in the directions of radius, tangent, and longitude, respectively; is the value of the logarithmic bulk modulus, a dimensionless material parameter, which is the slope of the curve between void ratio and the logarithm of the equivalent pressure stress when wood is confined compressed, and it describes the detailed shape of void ratio vs. stress; and is the tensile strength of wood.
After using the equivalent pressure stress, the total stress and deviator stress have the following relationships [7]:where is the total stress; is the deviator stress; is a unit matrix; is the shear modulus; and is the deviator strain [7], which can be expressed aswhere is the total strain and is the logarithmic volumetric strain; and , where V is the volume of wood and V_{0} is the initial volume of wood. Then, the stress and strain can be solved by the above equations.
3. Confined Compression Test of Wood
A setup for the confined compression test was devised to conduct the experiments for obtaining the relationship between void ratio of wood and equivalent pressure in the direction of parallel to grain, shown as in Figure 1(a). A wood cylinder made of the species of Pinus Sylvestris var. mongolica was inserted into the steel block through a round hole bored in the block. Then, the wood cylinder was compressed by a steel cylinder which was fit to the hole. The diameters of the wood cylinder and the hole are both 20 mm. After sanded by the sand paper, the wood cylinder was inserted in the hole. During the test, the compression velocity of the cylinder was controlled at a rate of 0.5 mm/min [5, 8] and the compression load and displacement were recorded. The information for wood cylinders is shown in Table 1.
(a)
(b)

Seven wood cylinders with different initial density and different initial void ratio were compressed parallel to grain. Figure 1(b) shows the typical crushing failure mode of the cylinder. The wood cylinders were shortened, and diagonal cracks appeared. Figure 2 shows the curves of stress versus elongation parallel to grain from the confined compression tests of wood cylinders.
Based on the relationship among three components including wood substance, water, and void [9], the initial void ratio of wood can be expressed bywhere , , and are initial void ratio, moisture content, and density of wood, respectively; is the density of the wood solids, for different wood species, g/cm^{3}, and here an average value [6] was applied; and S, l_{0}, and m_{0} are the crosssection area, initial length, and mass of wood cylinder specimens, respectively. The initial void for wood cylinder specimens is shown in Table 1.
For the confined compression test, the longitudinal strain is not zero, while the tangent strain and radial strain are both zero, which can be derived by the following equation:
Then, the radial stress and tangent stress can be expressed by the longitudinal stress, shown in the following equations:
For the tested wood species, the elastic ratio and Poisson’s ratio [10] are shown in Table 2. After applying the elastic ratio and Poisson’s ratio, the equivalent stress of the confined compression wood cylinder can be solved, , , and , respectively.

The change of void ratio of the wood cylinder during compressing can be solved by the compression deformation, and its initial void ratio and height can be expressed aswhere is the change of void ratio of wood, is the initial length of the wood cylinder, and is the compressive deformation of the wood cylinder. After calculating the equivalent stress and , curves of the equivalent stress versus the change of void ratio can be obtained, shown as in Figure 3.
(a)
(b)
Comparing the curves from the experiments of seven wood cylinders in Figure 3(a), the curves of different cylinders with different initial void ratios are close to each other except the specimen S5. However, an average curve of equivalent stress versus change of void ratio can reflect the compression behaviour of wood with the boundary condition of confined restrain. Based on equation (1), the logarithmic bulk modulus of wood was calculated by the average curve. From the average curve shown in Figure 3(b), it can be found that the logarithmic bulk modulus is not a constant, which reflects three stages of the contact between wood groove and bolt. At the beginning of compression, the logarithmic bulk modulus decreases rapidly until bolt and wood groove contact closely. Then, it nearly keeps a constant. Finally, it climbs rapidly at the same time plasticity is yielded in wood groove. Here, a regressed relationship between the logarithmic bulk modulus and the change of void ratio was expressed by equation (10) using a polynomial based on the average curve:
4. FE Modeling and Validation of Embedment Crushing of Wood Groove
4.1. 3D FE Modeling for Bolted Wood Connection
A 3D FE modeling for bolted connection in wood was created using C3D8 solid elements in ABAQUS [7]. Surface contact between wood groove and steel bolt was taken into consideration with the coefficient of friction of 0.15 [11]. The material property of perfect porous elastoplasticity was assigned to the wood groove zone around the three times of bolt diameter according to crushing deformation in experiment results and FE calculation. A Drucker–Prager plasticity was assigned to wood in the groove zone. The material angle of friction of wood was set as 11°, and dilation angle was zero [4]. An orthotropic elasticity property was assigned to other zones of wood. A perfect elastoplasticity property was assigned to the steel bolt. The requisite material properties used in FE modeling are listed in Table 3. Here, the properties of clear wood were assigned to the pin groove due to the fact that the wood without defects was usually chosen to make joints and only localized crushing occurred at the joint zone with enough end distance (7d). The dependent relationship of the logarithmic bulk modulus on change of void ratio was implemented in a user defined subroutine UFIELD in which a user defined field variable was defined as the void ratio of each integrate point in the FE software ABAQUS [7].
 
The data are from experiments. 
4.2. Test of Single Bolted Wood Connection
To validate the FE modeling with the perfect elastoplastic porous constitutive model, 12 lateral resistance tests of single bolt connection of wood made of the species of Pinus sylvestris var. mongolica were conducted in the direction of parallel to grain according to ASTM D5652 [8]. The dimensions of centre member and side member of connections are all 150 mm in the direction of parallel to grain and 80 mm in the direction of perpendicular to grain, shown as in Figure 4. The diameters of the bolt and groove are 10.5 mm and 13 mm, respectively. The side distance is 40 mm and end distance is 100 mm for every connections. Four groups of thickness combinations for connection specimens, three samples per combination, were prepared, whose detail thicknesses are shown in Table 4.

4.3. Validation of FE Modeling
The deformations of all tested connections from different thickness combinations are shown in Figure 5. For each combination, one typical deformed connection and all three deformed bolts are exhibited in the same figure, including these from the tested connection and numerical modeling. The deformation of FE modeling shows good agreement with the experiment observation. Localized crushing failure of pin groove and plastic hinges in the steel bolts appears in different connections. It is obvious that crushing failure did appear at the compressive pin groove, and stress contours for combination 1 are shown in Figure 6. The stress concentration phenomenon and local crushing deformation occur under the contact surface between bolt and groove. The number of plastic hinges in bolts depends on the thickness ratio between the centre member and side member. For the combination 2, the thickness of side member is the thinnest among the four combinations, crushing failure appeared in side member, and the rigid body deformation happened in the bolt. Hence, there was no plastic hinge in the bolt in combination 2. Meanwhile, in combination 4, there were three plastic hinges in the bolt and two hinges in each shear plane. The FE modeling with perfect elastoplastic porous constitutive exhibits the practicability to simulate the deformation of wood connection including localized deformation of groove and bolt’s deformation.
(a)
(b)
(c)
(d)
(a)
(b)
The loaddisplacement curves for all tested connections and FE simulations are shown in Figure 7. All curves showed the changes of stiffness of bolted connection. The stiffness becomes less and less with the increase in deformation when the localized wood and steel bolt come into plasticity. The ultimate load from FE results was very close to the maximum loads among three experiment results. While the diameter of the bolt is large and wood members are thin, the yield of groove also can be simulated. The six yield modes for single shear and four yield modes for double shear can all be modeled by the FE modeling. Comparing the deformations of experiments and FE simulations, although the initial stiffness by FE modeling is a little bigger than those of experiments, the overall simulated curves are general accordant with those of experiment curves. The possible reason is that the contact gap between bolt and wood groove brings additional deformation. Hence, the initial stiffness in simulation is more than that of experiment.
(a)
(b)
(c)
(d)
5. Conclusions
A method to simulate the localized crushing failure of wood groove in bolted wood connection applying a perfect elastoplasticity porous constitutive model was developed. According to the threephase composites of wood with solid, water, and gas, wood void ratio expression was derived. A confined compression of the wood cylinder was conducted for void ratio versus stress relationship. The tests can simulate the real stress state in wood groove and provide the FE parameter of wood species for numerical simulation. A relationship between wood void ratio and compression stress was established. The empirical constitutive model for localized crushing wood was validated by the good agreement between experiments and numerical simulations for bolted wood connection with enough end distance (7d). The FE modeling also can predict the load bearing capacity of the all yield modes for single shear and double shear. Its application may be expanded to simulate complex bolted connection in wood structures, such as multiple fastener connections.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors thank the financial support of National Key Research and Development Program of China (no. 2019YFD1101001).
References
 D. M. Moses, “Constitutive and analytical models for structural composite lumber with applications to bolted connections,” University of British Columbia, Vancouver, Canada, 2000, Ph.D. Dissertation. View at: Google Scholar
 T. Zhou and Z. W. Guan, “A new approach to obtain flat nail embedding strength of doublesided nail plate joints,” Construction and Building Materials, vol. 25, no. 2, pp. 598–607, 2011. View at: Publisher Site  Google Scholar
 H. JungPyo and D. Barrett, “Threedimensional finiteelement modelling of nailed connections in wood,” Journal of Structural Engineering, vol. 136, no. 6, pp. 715–722, 2010. View at: Publisher Site  Google Scholar
 S. Carmen, “Mechanical behaviour of timber joints with slottedin steel plates,” Technology University Delft, Delfts, Netherlands, 2012, Ph.D. Dissertation. View at: Google Scholar
 ASTM D 576497a, Standard Test Method for Evaluating DowelBearing Strength of Wood and WoodBased Products, American Society for Testing and Materials, West Conshohocken, PA, USA, 2002.
 F. F. P. Kollmann, A. Wilfred, and J. R. Cote, Principles of Wood Science and Technology, 1 Solid Wood, SpringerVerlag New Yourk Inc, New York, NY, USA, 1968.
 D. S. Simulia, ABAQUS 6.14. ABAQUS Theory Guide, Dassault Systèmes Simulia Corp., Johnston, RI, USA, 2014.
 ASTM D 565295, Standard Test Methods for Bolted Connections in Wood and Wood–Based Products, American Society for Testing and Materials, West Conshohocken, PA, USA, 2002.
 M. D. Wood, Soil Behaviour and Critical State Soil Mechanics, Cambridge University Press, Cambridge, UK, 1990.
 D. W. Green, E. Winandy Jerrold, and E. Kretschmann David, Mechanical Properties of Wood, Wood as an Engineering Material, Forest Products Laboratory, Madison, WI, USA, 1999.
 W. Simpson and T. W. Anton, Physical Properties and Moisture Relations of Wood: Wood Handbook, Wood as an Engineering Material, Forest Products Laboratory, Madison, WI, USA, 1999.
Copyright
Copyright © 2020 Huazhang Zhou and Xiaoqiang Zhou. 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.