Numerical Analysis of the Bending Properties of Cathay Poplar Glulam

This paper presents the formulae and finite element analysis models for predicting the Modulus of Elastic (MOE) and Modulus of Rupture (MOR) of Cathay poplar finger-jointed glulam. The formula of the MOE predicts the MOE of Cathay poplar glulam glued with one-component polyurethane precisely. Three formulae are used to predict the MOR, and Equation (12) predicts the MOR of Cathay poplar glulam precisely. The finite element analysis simulation results of both the MOE and MOR are similar to the experimental results. The predicted results of the finite element analysis are shown to be more accurate than those of the formulae, because the finite element analysis considers the glue layers, but the formulae do not. Three types of typical failure modes due to bending were summarized. The bending properties of Cathay poplar glulam were compared to those of Douglas fir glulam. The results show that Cathay poplar glulam has a lower stiffness, but a marginally higher strength. One-component polyurethane adhesive is shown to be more effective than resorcinol formaldehyde resin adhesive for Cathay poplar glulam. This study shows that Cathay poplar has the potential to be a glulam material in China.


Introduction
Modern wooden structures have been used in many areas around the world. Glued-laminated timber (glulam) is a type of engineered wood product that can preserve the natural beauty of the source wood during construction [1,2]. Structural glulam has been widely used in Japan, North America and Europe [3], but Chinese domestic wood species are rarely used in the construction of structural glulam manufactured in China. It is thus important to study Chinese domestic wood species that could be applicable for use in structural glulam.
Cathay poplar (Populus cathayana) is one of the most abundant, fast-growing species in northern China and is cold resistant, fast growing, straight grained, medium and evenly textured and easy to manufacture. Cathay poplar is thus an excellent material for construction components [4,5].
Finite element analysis (FEA) is an efficient research method used in the wooden materials field. By establishing a finite element analysis model and post-processing, visual predictions can be obtained without specimen failures. In 1980, Foschi and Barrett put forward the Foschi-Barrett simulation model, which used the Monte Carlo method, to predict the Modulus of Rupture (MOR) of glulam. The Ala Tabiei [6] model simulated the nonlinearity of wood, and the Moses D.M. [7] model predicted the nonlinear deformation and brittle fracture of wood composite material. The Serrano [8] model is also a three-dimensional nonlinear FEA model, which predicted the properties of the glue layer, the size effect and tensile fracture properties when the wood was connected with bolts. Gao [9] proposed the FEA model of the deformation of timber-framed plywood panel dome structures, and the simulated results were found to be in good agreement with the results of the experiments under similar loading conditions. Frese [10][11][12] simulated the bending and tensile strengths of spruce and beech glulam. Qiu [13] simulated the crack propagation behavior of Chinese larch, and the numerical results correlated well with the experimental results. Based on the stiffness profile of each board and its location within the glulam, Kandler [14] established a numerical finite element model, which is able to predict the effective glulam stiffness with high accuracy.
The formula model is also commonly used to make numerical simulation of glulam. Kohler [15], Fink [16,17] and Jockwer [18][19][20] have published contributions about the formula model. Kohler presented a probabilistic approach for modelling the tensile strength and tensile stiffness properties of timber boards and finger joint connections of knot clusters based on the results of the experimental investigation. Fink presented a probabilistic approach for modelling the load-bearing capacity, the bending stiffness and the failure modes of glulam. In Jockwer's research, the prediction of the load-carrying capacity of notched beams was developed based on experiments and theory, and the load at initial cracking of the reinforced notches tested can be predicted with the model. Besides, the impact of knots and grain deviations on the fracture perpendicular to the grain of timber is analyzed by means of numerical models.
The Modulus of Rupture (MOE) and MOR are both paramount properties of glulam. Four formulae were derived, and finite element analysis models were established to predict the MOE and MOR of Cathay poplar finger-jointed glulam. In addition, the bending properties of Cathay poplar glulam were examined and compared to those of Douglas fir glulam. A comparison of the test and simulation results was also performed.

Material Properties
Cathay poplar was obtained from the Beijing Longshun Wood Market, Beijing, China. Douglas fir was imported from North America (Can-Coast Development Corp, Vancouver, BC, Canada). The moisture content of the Cathay poplar and the Douglas fir ranged from 9% to 12%. The physical and mechanical properties of the Cathay poplar and Douglas fir lumber are shown in Table 1. Resorcinol formaldehyde resin adhesive (RF, Prefere 5837W, Dynea, Gaoyao, Guangdong, China) and one-component polyurethane adhesive (PU, HB S309, Purbond, Shanghai, China) were used to manufacture the glulam specimen. The resin content of RF and PU is 100% and 45%~75%, respectively; the curing time of RF and PU is 75 min and 60 min, respectively.

Grading of Sawn Lumber
Visual grading and the MOE test method were used in this study. A total of 52 pieces of Cathay poplar sawn lumbers with dimensions of 60 (width (W)) mmˆ800 (length (L)) mm was tested, where half of the pieces were 20 mm thick and the other half were 30 mm thick. The lumber was graded based on GB/T (Recommended National Standards of China) 26899-2011 [21]. The concentrated knot's diameter ratio (CKDR, the ratio of the diameters of all knots to the lumbers' width in a 150 mm-long range) of Grades I, II and III lumber is not more than 20%, 30% and 40%, respectively. Figure 1 shows that the MOE decreases as the ratio increases. When the ratio is less than 30%, the MOE of the 20 mm-thick lumber is higher than that of the 30 mm-thick (T) lumber. However, when the ratio is in the range of 30%~40%, the linear fit of the MOE of the 30 mm-thick lumber is higher than that of the 20 mm-thick lumber.

Grading of Sawn Lumber
Visual grading and the MOE test method were used in this study. A total of 52 pieces of Cathay poplar sawn lumbers with dimensions of 60 (width (W)) mm × 800 (length (L)) mm was tested, where half of the pieces were 20 mm thick and the other half were 30 mm thick. The lumber was graded based on GB/T (Recommended National Standards of China) 26899-2011 [21]. The concentrated knot's diameter ratio (CKDR, the ratio of the diameters of all knots to the lumbers' width in a 150 mm-long range) of Grades I, II and III lumber is not more than 20%, 30% and 40%, respectively. Figure 1 shows that the MOE decreases as the ratio increases. When the ratio is less than 30%, the MOE of the 20 mm-thick lumber is higher than that of the 30 mm-thick (T) lumber. However, when the ratio is in the range of 30%~40%, the linear fit of the MOE of the 30 mm-thick lumber is higher than that of the 20 mm-thick lumber.

Manufacturer of Specimen
Glulam specimens were manufactured based on GB/T 26899-2011. All of the glulam specimens were 3000 (L) mm × 60 (W) mm × 120 (T) mm. The length of finger joints was 21.5 mm, and the width was 8.0 mm. The distribution of finger joints in glulam satisfies the construction requirement of glulam. The positions of finger joints were at least 30 cm away from each adjacent layer.
The specimen parameters are shown in Table 2. To improve the bending strength of glulam [22], Grades I and II lumber were placed on the outermost layers of the 4-layer glulam and the outermost and outer layers of the 6-layer glulam. Grade III lumber was placed on the inner layers ( Figure 2).

Manufacturer of Specimen
Glulam specimens were manufactured based on GB/T 26899-2011. All of the glulam specimens were 3000 (L) mmˆ60 (W) mmˆ120 (T) mm. The length of finger joints was 21.5 mm, and the width was 8.0 mm. The distribution of finger joints in glulam satisfies the construction requirement of glulam. The positions of finger joints were at least 30 cm away from each adjacent layer.
The specimen parameters are shown in Table 2. To improve the bending strength of glulam [22], Grades I and II lumber were placed on the outermost layers of the 4-layer glulam and the outermost and outer layers of the 6-layer glulam. Grade III lumber was placed on the inner layers ( Figure 2).

Equipment and Experiment
A four-point load static bending test (i.e., bending test Method A) was performed based on GB/T 26899-2011 at a loading speed of 8 mm/min. An LVDT (Linear Variable Differential Transformer) sensor was set at the midspan of each specimen. (Figure 3). The MOE for bending and MOR are calculated by the following Equation: (1) where ∆ is the difference between the upper and lower loads at the proportional limit, l is the span of the specimen between supports, s is the span of the specimen between the loading points, Δy is the corresponding midspan deflection of ∆ , b is the width of the specimen, h is the depth of specimen and Pmax is the maximum load.

Failure Modes
Elastic deformation occurred in the specimens during the initial stage of loading. As the load increased, plastic deformation occurred in the specimens. At this stage, many small compression wrinkles appeared in the specimens. The bending rigidity of the specimen decreased marginally, and the deformations increased significantly. Finally, the bottom laminates reached their ultimate tensile stresses. The finger joints of the bottom laminate or defects in the timber failed first, which led to brittle failure of the specimen. The three typical failure modes observed are shown in Figure 4.
The failure modes were connected to the finger joint to some degree. Structural glulam requires a large number of finger joints, so studies of finger joint fatigue strength are very important [23]. There are quite a few research works aimed at the relationship between the bending strength of glulam and finger joint strength.

Equipment and Experiment
A four-point load static bending test (i.e., bending test Method A) was performed based on GB/T 26899-2011 at a loading speed of 8 mm/min. An LVDT (Linear Variable Differential Transformer) sensor was set at the midspan of each specimen. (Figure 3). The MOE for bending and MOR are calculated by the following Equation: MOE " ∆P pl´sq`2l 2`2 ls´s 28 ∆ybh 3 MOR " 3P max pl´sq 2bh 2 (2) where ∆P is the difference between the upper and lower loads at the proportional limit, l is the span of the specimen between supports, s is the span of the specimen between the loading points, ∆y is the corresponding midspan deflection of ∆P, b is the width of the specimen, h is the depth of specimen and P max is the maximum load.

Equipment and Experiment
A four-point load static bending test (i.e., bending test Method A) was performed based on GB/T 26899-2011 at a loading speed of 8 mm/min. An LVDT (Linear Variable Differential Transformer) sensor was set at the midspan of each specimen. (Figure 3). The MOE for bending and MOR are calculated by the following Equation: (1) where ∆ is the difference between the upper and lower loads at the proportional limit, l is the span of the specimen between supports, s is the span of the specimen between the loading points, Δy is the corresponding midspan deflection of ∆ , b is the width of the specimen, h is the depth of specimen and Pmax is the maximum load.

Failure Modes
Elastic deformation occurred in the specimens during the initial stage of loading. As the load increased, plastic deformation occurred in the specimens. At this stage, many small compression wrinkles appeared in the specimens. The bending rigidity of the specimen decreased marginally, and the deformations increased significantly. Finally, the bottom laminates reached their ultimate tensile stresses. The finger joints of the bottom laminate or defects in the timber failed first, which led to brittle failure of the specimen. The three typical failure modes observed are shown in Figure 4.
The failure modes were connected to the finger joint to some degree. Structural glulam requires a large number of finger joints, so studies of finger joint fatigue strength are very important [23]. There are quite a few research works aimed at the relationship between the bending strength of glulam and finger joint strength.

Failure Modes
Elastic deformation occurred in the specimens during the initial stage of loading. As the load increased, plastic deformation occurred in the specimens. At this stage, many small compression wrinkles appeared in the specimens. The bending rigidity of the specimen decreased marginally, and the deformations increased significantly. Finally, the bottom laminates reached their ultimate tensile stresses. The finger joints of the bottom laminate or defects in the timber failed first, which led to brittle failure of the specimen. The three typical failure modes observed are shown in Figure 4.
The failure modes were connected to the finger joint to some degree. Structural glulam requires a large number of finger joints, so studies of finger joint fatigue strength are very important [23]. There are quite a few research works aimed at the relationship between the bending strength of glulam and finger joint strength.  In terms of experimental research, Ohuchi [24,25] came to a conclusion based on bending tests of hinoki (Chamaecyparis obtusa) and sugi finger-jointed laminae. In these studies, to examine the optimum adhesive condition in the finger-joint part that influences the strength properties of the large-scale finger-jointed laminae, the block shear tests with block specimens that assumed the finger-joint part and glued under various adhesive conditions was performed. In conclusion, for glulam with good strength properties, the evaluation of finger-joint properties is important. In addition, Frese [26] tested the characteristic bending strength of beech glulam and characteristic finger joint bending strength. The conclusion was proposed that the bending strength of glulam depends on both board strength and finger joint strength.
In terms of the glulam strength model, Frese [11] put forward a glulam model, where its bending strength depends both on the board tensile strength and finger joint tensile strength. This model could determine the influence of the board and finger joint strength on glulam bending strength, which proved the relationship between the bending strength of glulam and finger joint strength. As a result, the poor characteristic bending strength values of the 38 test beams were mainly caused by too low requirements for boards and finger joints, which verified the model. In addition, Frese [10,27] proposed a predictive model for characteristic glulam bending strength based on tension tests of the boards and the finger joints. The experimental values of them were used as input data for the strength models to calculate the predictive values of the glulam. This model predicted the characteristic glulam bending strength precisely. Comparisons between the experimental data and the analytical results from the computer model show a good agreement. The predictive values for the characteristic glulam bending strength and the experimental ones differed not more than 6%. This indicated that the bending strength of glulam was determined by finger joint strength to a large extent.
Therefore, based on the above research, the following conclusions were proposed. Shear failure of the finger joint led to the failures of PR4 and PR6. The bending strengths of PR4 and PR6 were indicative of the bonding strength of a finger joint with RF, but did not fully describe the mechanical properties of Cathay poplar.
Shear failure of the finger joint and splitting of timber caused the failure of PO6. The bending strength of PO6 was indicative of both the bonding strength of a finger joint with PU and the In terms of experimental research, Ohuchi [24,25] came to a conclusion based on bending tests of hinoki (Chamaecyparis obtusa) and sugi finger-jointed laminae. In these studies, to examine the optimum adhesive condition in the finger-joint part that influences the strength properties of the large-scale finger-jointed laminae, the block shear tests with block specimens that assumed the finger-joint part and glued under various adhesive conditions was performed. In conclusion, for glulam with good strength properties, the evaluation of finger-joint properties is important. In addition, Frese [26] tested the characteristic bending strength of beech glulam and characteristic finger joint bending strength. The conclusion was proposed that the bending strength of glulam depends on both board strength and finger joint strength.
In terms of the glulam strength model, Frese [11] put forward a glulam model, where its bending strength depends both on the board tensile strength and finger joint tensile strength. This model could determine the influence of the board and finger joint strength on glulam bending strength, which proved the relationship between the bending strength of glulam and finger joint strength. As a result, the poor characteristic bending strength values of the 38 test beams were mainly caused by too low requirements for boards and finger joints, which verified the model. In addition, Frese [10,27] proposed a predictive model for characteristic glulam bending strength based on tension tests of the boards and the finger joints. The experimental values of them were used as input data for the strength models to calculate the predictive values of the glulam. This model predicted the characteristic glulam bending strength precisely. Comparisons between the experimental data and the analytical results from the computer model show a good agreement. The predictive values for the characteristic glulam bending strength and the experimental ones differed not more than 6%. This indicated that the bending strength of glulam was determined by finger joint strength to a large extent.
Therefore, based on the above research, the following conclusions were proposed. Shear failure of the finger joint led to the failures of PR4 and PR6. The bending strengths of PR4 and PR6 were indicative of the bonding strength of a finger joint with RF, but did not fully describe the mechanical properties of Cathay poplar.
Shear failure of the finger joint and splitting of timber caused the failure of PO6. The bending strength of PO6 was indicative of both the bonding strength of a finger joint with PU and the mechanical properties of Cathay poplar. Compared to the RF, the Cathay poplar glued with PU performs better.
Splitting of timber resulted in the failures of DO4 and DO6. The bending strengths of DO4 and DO6 were indicative of the mechanical properties of Douglas fir, but did not fully describe the bonding strength of a finger joint with PU.

Experimental Results
The load displacement curves of the mean value of the specimens from each series are shown in Figure 5 as solid lines, and the single results are shown as dotted lines. Compared to the Douglas fir glulam, the Cathay poplar glulam has a lower stiffness, but a slightly higher strength. The PU is shown to be more effective than the RF when used to glue the Cathay poplar glulam specimens. With the same thickness, the strength of the specimen increased when the number of layers increased. Splitting of timber resulted in the failures of DO4 and DO6. The bending strengths of DO4 and DO6 were indicative of the mechanical properties of Douglas fir, but did not fully describe the bonding strength of a finger joint with PU.

Experimental Results
The load displacement curves of the mean value of the specimens from each series are shown in Figure 5 as solid lines, and the single results are shown as dotted lines. Compared to the Douglas fir glulam, the Cathay poplar glulam has a lower stiffness, but a slightly higher strength. The PU is shown to be more effective than the RF when used to glue the Cathay poplar glulam specimens. With the same thickness, the strength of the specimen increased when the number of layers increased.   Table 3. Based on these results, the following conclusions are proposed.    Table 3. Based on these results, the following conclusions are proposed. Splitting of timber resulted in the failures of DO4 and DO6. The bending strengths of DO4 and DO6 were indicative of the mechanical properties of Douglas fir, but did not fully describe the bonding strength of a finger joint with PU.

Experimental Results
The load displacement curves of the mean value of the specimens from each series are shown in Figure 5 as solid lines, and the single results are shown as dotted lines. Compared to the Douglas fir glulam, the Cathay poplar glulam has a lower stiffness, but a slightly higher strength. The PU is shown to be more effective than the RF when used to glue the Cathay poplar glulam specimens. With the same thickness, the strength of the specimen increased when the number of layers increased.   Table 3. Based on these results, the following conclusions are proposed.   Firstly, a comparison of the different wood species was performed. The MOR of PR6 was 2.4% higher than that of DO6, and the MOE of DO6 was 70.6% higher than that of PO6. In this study, the Cathay poplar glulam had a lower stiffness, but a marginally higher strength compared to those of the Douglas fir glulam. The ultimate tensile strain of PO6 was 0.0053, and the maximum compressive strain was 0.0063; these values were both higher than those of Douglas fir glulam.
Secondly, a comparison of the different adhesives was performed. Based on the experimental results, the MOR of PO6 was 2.6% higher than that of PR6, and the MOE of PR6 was 32.9% higher than that of PO6. For the Cathay poplar glulam, the PU is shown to be more effective than the RF.
Lastly, a comparison of the different numbers of layers was performed. The MOR of PR6 was 1.2% higher than that of PR4, and the MOR of DO6 was 6.4% higher than that of DO4. The MOE of PR6 was 14.3% higher than that of PR4, and the MOE of DO6 was 4.3% higher than that of DO4. With the same thickness, the bending strength increased as the number of layers increased.

Checking Calculation
The bending strength, the shear strength and the deflection of the glulam specimens were checked.
The load of a room was assumed to be carried by sample PO6. Based on the characteristic load combination in GB 50009-2012 [28], the design value of the load was 4000 N. According to the four-point load static bending test method, the design value of the maximum moment M and maximum shear strength V were 1.68ˆ10 6 N¨mm and 2000 N based on the principle of the balance of force and moment, respectively. The shear diagram and bending moment diagram are shown as Figure 7.  Firstly, a comparison of the different wood species was performed. The MOR of PR6 was 2.4% higher than that of DO6, and the MOE of DO6 was 70.6% higher than that of PO6. In this study, the Cathay poplar glulam had a lower stiffness, but a marginally higher strength compared to those of the Douglas fir glulam. The ultimate tensile strain of PO6 was 0.0053, and the maximum compressive strain was 0.0063; these values were both higher than those of Douglas fir glulam.
Secondly, a comparison of the different adhesives was performed. Based on the experimental results, the MOR of PO6 was 2.6% higher than that of PR6, and the MOE of PR6 was 32.9% higher than that of PO6. For the Cathay poplar glulam, the PU is shown to be more effective than the RF.
Lastly, a comparison of the different numbers of layers was performed. The MOR of PR6 was 1.2% higher than that of PR4, and the MOR of DO6 was 6.4% higher than that of DO4. The MOE of PR6 was 14.3% higher than that of PR4, and the MOE of DO6 was 4.3% higher than that of DO4. With the same thickness, the bending strength increased as the number of layers increased.

Checking Calculation
The bending strength, the shear strength and the deflection of the glulam specimens were checked. The load of a room was assumed to be carried by sample PO6. Based on the characteristic load combination in GB 50009-2012 [28], the design value of the load was 4000 N. According to the four-point load static bending test method, the design value of the maximum moment M and maximum shear strength V were 1.68 × 10 6 N•mm and 2000 N based on the principle of the balance of force and moment, respectively. The shear diagram and bending moment diagram are shown as Figure 7.  The bending strength was checked by Equation (3). S is the section modulus, which was 1.44ˆ10 5 mm 3 , and M was 1.68ˆ10 6 N¨mm.
The value of checking calculation was 11.67 MPa, which was lower than the design value of bending strength f m (17 MPa). In conclusion, PO6 achieves the required bending strength in GB 50009-2012.

Checking Calculation of the Horizontal Shearing Strength
The shearing strength was checked by Equation (4). Q is the first moment of area, which was 1.08ˆ10 5 mm 3 ; I is the second moment of area (i.e., moment of inertia), which was 8.64ˆ10 6 mm 4 ; b is the width of the specimen, which was 60 mm; and V was 2000 N.
The value of checking calculation was 0.42 MPa, which was lower than design value of horizontal shearing strength f m (2.4 MPa). In conclusion, PO6 achieves the required horizontal shearing strength in GB 50009-2012 [28].

Checking Calculation of the Deflection
The deflection was checked by Equation (5). δ 1 is the design value of the deflection, which was 8.64 mm, and δ is the calculated deflection, which was calculated by Equation (6).
The MOE of PO6 was 11,119 MPa, so the value of checking calculation was 8.14 mm, which was lower than the design value of deflection. In conclusion, PO6 satisfies the deflection requirements in GB 50009-2012 [28]. Therefore, the bending strength, shear strength and deflection of PO6 pass the checking calculations.

MOE of the Cathay Poplar Glulam
The MOE for bending of the Cathay poplar glulam was calculated based on that of sawn lumber using a formula that will now be described. The accuracy of the formula is also verified below.
Wang and Chang [30] showed that the strength of laminate was based on the bending stiffness of laminates (e j I j = the bending stiffness of the j laminate). The bending stiffness of a glulam specimen was the sum of the bending stiffness of each laminate layer under the non-glued layers, as shown in Equation (7) [31]: where E b is the MOE of the glulam specimen; I is the section moment of inertia of the glulam specimen; and e j I j is the MOE multiplied by the section moment of inertia of the j laminate (j = 1, 2, 3 . . . , 6), respectively. Based on the parallel axis theorem from material mechanics [32], the section moment of inertia of laminates can be transformed into the section moment of inertia of the central plane of the glulam specimen as Equation (8): where I j0 is the section moment of inertia of the j-th laminate to the central plane of the glulam specimen; A j is the sectional area of the j-th laminate; and d j is the distance between the j-th laminate and the central plane of the glulam specimen. Equation (8) was then substituted into Equation (7): The mean values of the theoretical results of Equation (9) are shown in Table 4. The error between the test and theoretical results of PO6 was 2.89%; thus, the formula predicts the MOE of PO6 accurately. However, the error of PR4 and PR6 is somewhat larger; this is probably because the shear failure of finger joint led to the failures of PR4 and PR6, and thus, the bending properties of PR4 and PR6 did not fully describe the mechanical properties of the Cathay poplar. Therefore, the calculated MOE of PO6 is more precise than that of PR4 and PR6 based on existing data. Moreover, Equation (9) would be validated with more data by further experiments to test its accuracy.

MOR of the Cathay Poplar Glulam
The MOR of glulam is typically significantly better compared to that of sawn lumber due to the distribution of defects in the lumber [33]. The MOR of glulam is closely related to its material properties. To explore the relationship between the MOR of glulam and the tensile strength or the MOR of the lumber used, three formulae were derived.

Equation Referenced EN 1194 [34]
Based on many experimental results, EN (European Norm) 1194 (Europe Standard Committee) has shown that the relationship between the MOR of the bending of glulam and the tensile strength of the lumber used is linear. A laminate coefficient λ is thus introduced and is expressed in the following two Equation: f b,gl,k " 7.35`1.12 f t,lam,k λ " 1.12`7.35{ f t,lam,k where f b,gl,k is the characteristic MOR of glulam and f t,lam,k is the characteristic value of the tensile strength of the lumber. The laminate coefficient λ is inversely proportional to the characteristic value of the tensile strength of the lumber f t,lam,k . For lumber that has higher tensile properties, its laminate coefficient is smaller than lumber that has lower tensile properties. Therefore, the bending properties of glulam made of worse lumber have greatly improved compared to glulam made of better lumber.
The correlation coefficient r = 0.945 and the two formulae are only suitable for glulam with a thickness of 600 mm. If the dimensions are different, the MOR should be multiplied by the volume factor k [35,36], where F b is the MOR of glulam with other dimensions, F b0 is the MOR of glulam with standard dimensions, V 0 is the standard volume (130 mmˆ300 mmˆ6400 mm), V is the measured volume and k is the volume factor = 0.076 [35].

Regression Analysis
Based on the bending and tensile results, regression analysis is used to describe the relationship between the MOR of glulam and the tensile strength of the lumber used. The tensile strengths of the lumber mean the mean experimental values of 20 mm-thick lumber and the mean experimental values of 30 mm-thick lumber. In this formula, the tensile strength of the lumber is an independent variable, and the MOR of glulam is a dependent variable. Based on the scatter diagram drawn by Origin (Figure 8), the relationship appears to be linear. The correlation coefficient r = 0.945 and the two formulae are only suitable for glulam with a thickness of 600 mm. If the dimensions are different, the MOR should be multiplied by the volume factor k [35,36], where Fb is the MOR of glulam with other dimensions, Fb0 is the MOR of glulam with standard dimensions, V0 is the standard volume (130 mm × 300 mm × 6400 mm), V is the measured volume and k is the volume factor = 0.076 [35].

Regression Analysis
Based on the bending and tensile results, regression analysis is used to describe the relationship between the MOR of glulam and the tensile strength of the lumber used. The tensile strengths of the lumber mean the mean experimental values of 20 mm-thick lumber and the mean experimental values of 30 mm-thick lumber. In this formula, the tensile strength of the lumber is an independent variable, and the MOR of glulam is a dependent variable. Based on the scatter diagram drawn by Origin (Figure 8), the relationship appears to be linear.

Equation Based on the Parallel Axis Theorem
Equation (14) was derived to describe the relationship between the MOR of glulam and the MOR of the lumber used. For symmetrical glulam, the MOR is expressed by the following Equation: σ (9) where σ is the MOR of glulam, E is the MOE for the bending of glulam and ε is the elongation of a unit length. Equation (8) in Section 4.1.1 was substituted into Equation (13) and Equation (14) to derive:

Equation Based on the Parallel Axis Theorem
Equation (14) was derived to describe the relationship between the MOR of glulam and the MOR of the lumber used. For symmetrical glulam, the MOR is expressed by the following Equation: σ " Eε (9) where σ is the MOR of glulam, E is the MOE for the bending of glulam and ε is the elongation of a unit length. Equation (8) in Section 4.1 was substituted into Equation (13) and Equation (14) to derive:

Comparison between Test and Theoretical Results
The mean values of the theoretical results of Equations (10), (12) and (14) are shown in Table 5. In terms of existing data, Equation (12) predicts the MOR well, but the errors of Equations (10) and (14) are somewhat large. Moreover, these three formula would be validated with more data by further experiment to test their accuracy. The theoretical results are shown to be 28.03% higher than the corrected experimental results on average. This is probably because Equation (10) was derived based on the experimental results of European glulam, and their mechanical properties have many differences from that of Cathay poplar glulam.
Equation (12) was derived based on the experimental results, and the theoretical results are similar to the experimental results. Equation (12) can predict the MOR of Cathay poplar glulam.
Equation (14) shows that the error between the test and calculated values of PR4 is somewhat larger. This is because of the softening of a finger-joint bond. This deviation leads to 23.04% lower theoretical results compared to the experimental results.

Finite Element Analysis
The MOE and MOR of Cathay poplar glulam are calculated using the finite element program ANSYS, Version 14.5. Figure 9 shows the meshed finite element model with the loads and constraints. The "8 nodes solid 185" element was chosen, and the material properties were defined based on the experimental results of Cathay poplar glulam. Cathay poplar is orthotropic, and the glue layers (the MOE is 3000 MPa, and Poisson's ratio µ is 0.37) are isotropic. Knots were not taken into consideration when the model was established. The compressive zone is assumed to be ideally elastoplastic, and the tensile zone is assumed to be ideally elastic. The calculation does not stop until the tensile stress in the center of an element lies in between a range˘0.5% of the tensile strength of the board section [10].

. Comparison between Test and Theoretical Results
The mean values of the theoretical results of Equations (10), (12) and (14) are shown in Table 5. In terms of existing data, Equation (12) predicts the MOR well, but the errors of Equations (10) and (14) are somewhat large. Moreover, these three formula would be validated with more data by further experiment to test their accuracy. The theoretical results are shown to be 28.03% higher than the corrected experimental results on average. This is probably because Equation (10) was derived based on the experimental results of European glulam, and their mechanical properties have many differences from that of Cathay poplar glulam.
Equation (12) was derived based on the experimental results, and the theoretical results are similar to the experimental results. Equation (12) can predict the MOR of Cathay poplar glulam.
Equation (14) shows that the error between the test and calculated values of PR4 is somewhat larger. This is because of the softening of a finger-joint bond. This deviation leads to 23.04% lower theoretical results compared to the experimental results.

Finite Element Analysis
The MOE and MOR of Cathay poplar glulam are calculated using the finite element program ANSYS, Version 14.5. Figure 9 shows the meshed finite element model with the loads and constraints. The "8 nodes solid 185" element was chosen, and the material properties were defined based on the experimental results of Cathay poplar glulam. Cathay poplar is orthotropic, and the glue layers (the MOE is 3000 MPa, and Poisson's ratio μ is 0.37) are isotropic. Knots were not taken into consideration when the model was established. The compressive zone is assumed to be ideally elastoplastic, and the tensile zone is assumed to be ideally elastic. The calculation does not stop until the tensile stress in the center of an element lies in between a range ±0.5% of the tensile strength of the board section [10].  In addition, four elastic cushion blocks are set at two of the supports and two loading points to mitigate stress concentrations and to create a more accurate simulation. Based on He [37] and Tuohuti [38], the elastic cushion blocks are isotropic. The MOE is 1000 MPa, and Poisson's ratio µ is 0.2.
The mean values of the simulation results of the MOE for bending and MOR of Cathay poplar glulam are shown in Table 6. For the MOE, the simulation results are all higher than the experimental results. The mean error of PR4 is 8.98%, while that of PR6 is 7.67%; the error of PO6 is lower (2.79%). Thus, the simulation of PO6 is more accurate than those of PR4 and PR6.
For the MOR, the simulation results are all lower than the experimental results. The mean error of PR4 is´5.33%, that of PR6 is´5.39% and that of PO6 is´3.89%. Thus, the simulation of PO6 is more accurate than those of PR4 and PR6.
Therefore, the simulations of the MOE and MOR of PO6 are better than those of PR4 and PR6. Figure 10a compares the experimental results of the MOE for bending and MOR with the calculation and simulation results.
In addition, four elastic cushion blocks are set at two of the supports and two loading points to mitigate stress concentrations and to create a more accurate simulation. Based on He [37] and Tuohuti [38], the elastic cushion blocks are isotropic. The MOE is 1000 MPa, and Poisson's ratio μ is 0.2.
The mean values of the simulation results of the MOE for bending and MOR of Cathay poplar glulam are shown in Table 6. For the MOE, the simulation results are all higher than the experimental results. The mean error of PR4 is 8.98%, while that of PR6 is 7.67%; the error of PO6 is lower (2.79%). Thus, the simulation of PO6 is more accurate than those of PR4 and PR6.
For the MOR, the simulation results are all lower than the experimental results. The mean error of PR4 is −5.33%, that of PR6 is −5.39% and that of PO6 is −3.89%. Thus, the simulation of PO6 is more accurate than those of PR4 and PR6. Therefore, the simulations of the MOE and MOR of PO6 are better than those of PR4 and PR6. Figure 10a compares the experimental results of the MOE for bending and MOR with the calculation and simulation results. For the MOE, the calculation and simulation results of PO6 approach the experimental results, but the errors between the results of PR4 and PR6 are somewhat larger. Therefore, the formulae and simulation model are both suitable to predict the MOE of Cathay poplar glulam glued with PU.
For the MOR, the error between the test and simulation results is small, but the error between the test and theoretical results is somewhat larger. Therefore, the results of the finite element analysis are more accurate than those of the formulae. This is probably because the glue layers were considered in the simulation models, but the formulae were derived without considering the glue layers.

Conclusions
Based on the experiments, calculations and finite element analysis, the following conclusions can be drawn.
Cathay poplar glulam has a lower stiffness, but slightly higher strength compared to Douglas fir glulam in this experiment. The ultimate tensile strain of Cathay poplar glulam glued with For the MOE, the calculation and simulation results of PO6 approach the experimental results, but the errors between the results of PR4 and PR6 are somewhat larger. Therefore, the formulae and simulation model are both suitable to predict the MOE of Cathay poplar glulam glued with PU.
For the MOR, the error between the test and simulation results is small, but the error between the test and theoretical results is somewhat larger. Therefore, the results of the finite element analysis are more accurate than those of the formulae. This is probably because the glue layers were considered in the simulation models, but the formulae were derived without considering the glue layers.

Conclusions
Based on the experiments, calculations and finite element analysis, the following conclusions can be drawn.
Cathay poplar glulam has a lower stiffness, but slightly higher strength compared to Douglas fir glulam in this experiment. The ultimate tensile strain of Cathay poplar glulam glued with polyurethane is 0.0053, and its maximum compressive strain is 0.0063, both of which are higher than those of Douglas fir glulam. This study indicates that Cathay poplar has the potential to be glulam material in China. In terms of the adhesive, the MOR of six-layer Cathay poplar glulam glued with one-component polyurethane adhesive is 2.6% higher than that of six-layer Cathay poplar glulam glued with resorcinol formaldehyde resin adhesive, and the MOE of the latter is 32.9% higher than that of the former.
A formula is derived to predict the MOE of Cathay poplar glulam based on the MOE of its laminates. The results indicate that the formula could predict the MOE of six-layer Cathay poplar glulam glued with one-component polyurethane adhesive. Three formulae are used to predict the MOR of Cathay poplar glulam based on the MOR and the tensile strength of its laminates. Equation (12) is suitable to predict the MOR of Cathay poplar glulam. The simulation results of both the MOE and MOR of Cathay poplar glulam are similar to the experimental results. The simulation of Cathay poplar glulam glued with one-component polyurethane is shown to be more accurate than those of Cathay poplar glulam glued with resorcinol formaldehyde resin adhesive. The error of the MOE is shown to be 2.8%, and that of the MOR is shown to be´3.9%.