2.1. Shear Specimens Preparation and Fabrication
Seven flatsawn eastern white pine (Pinus strobus) boards obtained from a local timber supplier near Syracuse, NY were used for this project. The boards were NeLMA [9] finish grade with actual dimensions of 22 mm thickness, 197 mm width, and 3.65 m length. These boards had been kiln-dried and planed. They were stored in ambient conditions in the testing laboratory for 30 days, and the moisture content of all the boards was measured with a Delmhorst (BD 2100) pin-type moisture meter before fabrication. The average measured moisture content of the boards was approximately 11 ± 1%. These boards were planed to the desired thickness of 19 mm to assure uniform thickness and a fresh surface for correct adhesive performance. The boards were then cut to the desired length and orientation for the angled layers.
First, six two-layer panels with 19mm-thick layers, a width of 197mm, and a length of 3.65m with angle-ply orientations of 0°, ± 10°, ± 20°, ± 40°, ± 70°, and 90° were fabricated. The layers were cut to the desired size from the straight grain and defect-free parts of the boards to eliminate the possible effects of knots on the results. However, because of the NeLMA finish grade’s allowance of small knots, and the process of cutting the boards in different grain orientations, a minimal number of knots were included in some of the test specimens. Next, three replicates of each grain orientation panel were cut from the full-size panels to 38mm thickness, 152mm width, and 457mm length (Fig. 2). The width and length of the specimens in this study followed ASTM D2718 [8]. Loctite HB X452 Purbond adhesive (Henkel Corporation, Bridgewater, NJ) with a recommended spread rate of 180 g/m2 was used to glue the top and bottom of the lamina (face-glued). The grain orientations of the boards were alternated in the adjacent layer to equalize the effect of shear failure on the overall performance. A vacuum press set to maintain a negative pressure of 0.82 MPa (24 in Hg) was used to apply sufficient pressure for the gluing process. This pressure was found to be sufficient based on previously tested specimens showing no glue failures. The boards were under pressure for 42 hours meeting the minimum requirement of 24 hours and were then released and cut to the final dimension shear blocks.
As Fig. 3 illustrates, to perform the two-plate shear tests, the shear specimens were glued to two steel plates with dimensions of 13- mm thickness, 152- mm width, and 559- mm length with one side knife-edged at 20° on the top and bottom layers.
The steel plates were first heated in an oven at 150°C for one hour to provide dry metal surfaces and to improve the gluing process. Then a Loctite PL Premium Fast Grab Polyurethane adhesive with an approximate spread rate of 1500 g/m2 and an assembly time of 20 minutes was used to glue the specimens to the steel plates. As shown in Fig. 4, the specimens glued to the steel plates were placed under a pressure of 1.03 MPa using a Tinius Olsen universal testing machine for 48 hours at 25°C to ensure sufficient bonding.
After each test, the steel plates were prepared by removing wood residue from the previous test. Then they were scrubbed with a clean bristle brush with acetone solvent and then washed under cold running tap water and then placed in the oven again as described above.
2.2. Two-plate Planar Shear Test
The planar shear test, as defined in the European CLT standard EN 16351 − 2021 [11] and as a modified method of the shear compression test in ASTM D2718 [8], was conducted using a compression test fixture configuration with a Young Universal Mechanical Testing Machine. The test was performed to measure the angled and rolling shear modulus (Gα relative to the grain orientation of the angle-ply layers). As shown in Fig. 5, the experimental setup provides uniform shear stress in the layers and guarantees angled and rolling shear failure. The 10° angle (ϴ) of inclination was determined based on the length and thickness of the specimens [12]. Each shear specimen was loaded in compression until failure and the load-displacement curve was recorded to calculate the shear modulus (Gα). The shear modulus (Gα) was calculated using the slope between 40% and 50% of the peak load using Eq. 1,
\({\text{G}}_{{\alpha }}=\frac{t}{\text{l}\times \text{b}} \times \frac{\text{P}}{\varDelta } \times \text{cos}\text{ϴ}\) Eq. 1
where l is the length of inclined height and b is the width of the specimen, t is the thickness of the cross-layer, P/Δ is the slope of the load-displacement linear curve between 40% and 50% range from the first loading step (zero point) to max load, and ϴ is the angle between the shear plane and force direction which here is 10°.
2.3. Finite Element Method
The research version of the software ANSYS® Workbench 2022 [13] was used to provide numerical finite element analysis (FEA) for insight into the shear stress distribution in each specimen during the two-plate (planar) shear test. The boundary conditions of one of the knife-edge steel plates were assumed fixed and the other knife-edge of the steel plate was subjected to the compressive load in the form of displacement.
The specimens were modeled using linear elastic orthotropic wood materials. The nine independent elastic constants are required to define the mechanical response of an orthotropic material in the simulation. These constants are three elastic moduli (E), three Poisson’s ratios (ν), and three shear moduli (G). The actual test configuration was modeled for each shear specimen with two steel plates with a modulus of elasticity and Poisson’s ratio respectively set to 200 GPa and 0.30. For the wood material property values input in the simulations, the longitudinal elastic modulus value (MoE, EL) for eastern white pine was obtained from non-destructive simply supported bending tests on three eastern white pine boards (span to depth ratio > 15 [8]). The values for two constants of GLR (G0) and GTR (G90) were also acquired from the two-plate shear tests in the study. The USDA Wood Handbook[14] provides elastic ratios with respect to the longitudinal elastic modulus (EL) for some softwood and hardwood species. However, there are no available published values or elastic ratios for the constants of ER, GTL, ET, and the three required Poisson’s ratios for eastern white pine. To predict these elastic constants using the measured longitudinal (parallel to the grain) modulus of elasticity, linear regressions (confidence level of 95%) were performed on each constant of the species in the USDA Wood Handbook with available elastic ratios. The regressions were performed with Microsoft Excel[14] using the longitudinal MoE as the predictor. The regressions showed R2 ranging from 80–99% P < 0.05), indicating a statistically significant correlation between MoE and the response elastic constant (Appendix). Final elastic ratios and values for the material in FEM simulations have been reported in Table 1 and Table 2.
Table 1. Investigated Ratios of Eastern White Pine Mechanical Properties for the ANSYS Simulations Input
Species
|
Poisson's ratio
|
Shear Modulus
|
Young’s Modulus
|
νLT
|
νLR
|
νTR
|
GLT/ EL
|
ET/ EL
|
ER/ EL
|
Eastern White Pine
|
0.233
|
0.223
|
0.220
|
0.0693
|
0.0547
|
0.0896
|
Table 2. Final Material Properties of Eastern White Pine for ANSYS Simulations
Species
|
Poisson's ratio
|
Shear Modulus (MPa)
|
Young’s Modulus (MPa)
|
Specific Gravity
|
νLT
|
νLR
|
νTR
|
GLR
|
GLT
|
GTR
|
ET
|
ER
|
EL
|
SG
|
Eastern White Pine
|
0.233
|
0.223
|
0.220
|
465.0
|
520.8
|
37.5
|
411.1
|
673.4
|
7,515.3
|
0.35
|
FEM ANSYS Workbench [13] simulations were performed under the “Static Structural” analysis system on a 3D model of each shear specimen with grain orientations of 0°, ± 10°, ± 20°, ± 40°, ± 70°, and 90°. The shear blocks were modeled with the properties and configurations such as material, loading, and boundary conditions identical to those tested. Mesh development was performed under mechanical physics preference (program-controlled element order) by subdividing the mesh size to the rectangular elements with the size of 3.5mm for each solid geometry. The mesh elements used for the sample specimens were higher order elements: SOLID186, a 3-D 20-node element was used for solids, TARGE170, a 3-D target surface for the associated contact element of CONTA174, a 3-D 8-node surface-to-surface element. The mesh size was refined to ensure mesh convergence. Table 3 shows the details of the mesh sizing and the refinement quality of the models.
Table 3
Mesh Sizing and Details for ANSYS Modeling
Group
|
Bounding Box Diagonal Size (mm)
|
Average Surface Area (cm2)
|
Minimum Edge Length Size (mm)
|
Numbers of Nodes
|
Number of Elements
|
Specimen 0°
|
607
|
292
|
12.7
|
572,624
|
121,432
|
Specimen ± 10°
|
574
|
151
|
12.7
|
541,574
|
113,460
|
Specimen ± 20°
|
566
|
150
|
12.7
|
536,634
|
112,416
|
Specimen ± 40°
|
566
|
150
|
12.7
|
566,792
|
119,460
|
Specimen ± 70°
|
566
|
121
|
12.7
|
560,352
|
117,672
|
Specimen 90°
|
607
|
135
|
12.7
|
576,989
|
121,432
|
For each simulation, four coordinate systems were defined: one global coordinate system aligned with the testing machine load heads) two local coordinate systems aligned to the two diagonal layers with ± α° rotation, and one additional coordinate system parallel to the 10° angle of inclination (theta symbol here) of the test set-up. The absolute shear strain was measured with respect to this last coordinate system. The simulation’s details and the details of different coordinate systems designed in the modeling are shown in Fig. 6.
2.4. Evaluating the Shear Stiffness of Specimens Using Analytical Models
Finally, theoretical investigations were performed to evaluate the accuracy of Hankinson’s criterion and the RMS method by experimental and FEM results. The first tested model (Eq. 2) was initially developed by Hankinson in 1921 [16] and is still commonly used to describe uniaxial strength and/or stiffness with respect to the grain angle. Hankinson’s formula has been one of the most useful theories to predict the mechanical properties of solid wood when its grain orientation is at a diagonal angle. Also, a modified model was assessed using the Root-Mean-Square (RMS) method as shown in Eq. 3, which was applied by Gupta et al [6] to find the relationship between the shear strength and the grain angle.
\({\text{G}}_{{\alpha }}= \frac{{\text{G}}_{0 }\times {\text{G}}_{90}}{{\text{G}}_{0}\times {\text{sin}}^{2}{\alpha } +{\text{G}}_{90}\times {\text{cos}}^{2}{\alpha }}\) Eq. 2
\({\text{G}}_{{\alpha }}=\sqrt{{G}_{0}^{2}\times {\text{cos}}^{2}{\alpha }+{G}_{90}^{2}\times {\text{sin}}^{2}{\alpha } }\) Eq. 3
where, α is the angle of the grain orientation off the main axis, Gα is the shear modulus of the diagonal layer, G0, and G90 are shear moduli parallel and perpendicular to the grain.
This study has evaluated the accuracy of analytical models for predicting the shear stiffness values of angle-ply layers compared to the experimental tests and FEM results. Because these theoretical models use different methods to interpolate intermediate values, it is unclear which theory can provide an accurate relationship over the entire grain angle range [17].