2024
26
28032022
29112023
Fatih Kurul fatihkurul@istanbul.edu.tr
Istanbul University - Cerrahpaşa, Turkey
Ömer Asım Şişman
Gebze Technical University, Turkey
Türker Dündar
Istanbul University - Cerrahpaşa, Turkey
Abstract: For the mechanical characterization of Turkish fir and black pine, 400 board specimens with 22 mm× 50 mm × 420 mm were visually graded according to TS 1265 standard. Nondestructive tests were here upon performed using the stress wave method. After specimens were intentionally tested under flatwise bending to research the applicability as an alternative to tension and edgewise bending tests in European strength grading system. According to analyses of variance, the mean values of MOR and MOE differed in four groups at a p<0,05 significance level for visually graded boards. High correlations were found between MOR-MOE (R2=0,837) for fir and MOR-MOE (R2=0,776) for black pine. In addition, correlations of MOR-Knot rate for fir and black pine were respectively R2=0,669 and R2=0,660 showing the effectiveness of flatwise bending tests with the visual grading standard. For nondestructive tests, the mean values of the dynamic modulus of elasticity were very close in between fir and black pine grades while the usage of defect-free density performed better than the density of the whole specimen. Higher strength classes were found for black pine boards (Class 1= C40, Class 2= C27 and Class 3= C22) compared to fir boards (Class 1= C24, Class 2= C22 and Class 3= C18), respectively. Moreover, a simplified nonlinear material model was proposed for numerical modelling, and the results were found in good agreement in terms of the bending stiffness, strength, and deformation capacity of boards especially for class 1 and class 2 in both softwood species.
Keywords: Board, flatwise bending, finite element method, mechanical properties, stress wave, visual grading.
ARTÍCULO
Mechanical characterization of visually graded boards from turkish fir and black pine by nondestructive and destructive tests
Fatih Kurul fatihkurul@istanbul.edu.tr
Ömer Asım Şişman
Türker Dündar
Received: 28 March 2022
Accepted: 29 November 2023
Timber, as a natural and sustainable construction material, has variations in mechanical and physical properties due to anisotropy, defects and biological characteristics (Köhler et al. 2007). National standards such as DIN 4074-1 (2012), BS 4978: 2007+A2 (2017), UNI 11035-1 (2010) and UNE 56544 (2011) have been developed for visual grading of timber products based on knots, the slope of grain, strength, stiffness, density, and other visible parameters as described in EN 1912 (2012b). Timber products with a CE marking, representing the guaranteed resistance, are called strength graded timber that can be used to produce structural elements (Burawska-Kupniewska et al. 2020).
In addition to visual grading, nondestructive tests have been recently developed to evaluate the dynamic modulus of elasticity such as vibration, ultrasonic and stress wave methods and used in practice for an efficient grading of structural timber (Divós and Kiss 2010, Sanabria et al. 2011, Barriola et al. 2021).
Vibration method was prominent, however, the stress wave method was found to be more practical among nondestructive methods (Esteban et al. 2009, Barriola et al. 2020). Stress wave technique for wood quality assessment is based on the measurement of the velocity of a stress wave propagation generated by a shock. This technique is also applied for the production of engineered wood products (Bucur 2006, Ross 2015). According to Arriaga et al. (2022), velocities measured with the ultrasonic and stress wave methods had similar average value.
Moreover, the strength classes of structural timber are evaluated by machine grading that measure one or more mechanical and physical properties to establish its indicating properties (Mvolo et al. 2021, Vega et al. 2019). Kandler et al. (2018) experimentally tested glued laminated timber (GLT) beams to evaluate the impact of knots on bending strength and stiffness utilizing finite element model (FEM) with the image-processing.
For cross-laminated timber (CLT) panels produced from Australian pine, Navaratnam et al. (2020) predicted bending and shear performance using FEM analysis which was validated by experimental test data. Pang et al. (2021) investigated the effects of knots on the flexural performance of CLT panels to enhance the potential of using low-grade boards.
Furthermore, Crovella et al. (2019) verified shear analogy approach to evaluate bending stiffness of CLT panels made from softwood and hardwood. Rais et al. (2021) investigated European beach boards using 3D fiber orientation via laser technique for tensile strength prediction. The shift from structural timber to engineered wood products such as GLT post & beam and CLT wall or floor panels increased the speed and height of timber structures.
In recent years, the use of visual grading has gained momentum in Turkey for sustainable development and constructing timber buildings using green-certified products. Görgün and Dündar (2018) examined the compatibility between visual grades and the mechanical properties of timber joists from Turkish black pine determined by using nondestructive methods and edgewise bending tests. Moreover, Guntekin and Bulbul (2014) and Guntekin et al. (2013) used the stress wave method for the determination of the mechanical properties of timber joists from Turkish black and red pine, respectively.
Black pine is the second most common pine species following Turkish red pine, constituting 31 % of the softwood forest corresponding to 19 % of the total forestland in Turkey (OGM 2021). On the other hand, fir forest constitutes 3 % of the forestland in Turkey. Both species are also regarded as raw materials with a high potential for engineered wood products for structural usage in the timber industry.
Before the lamination process, boards shall be strength graded and declared as tabulated strength class in accordance with EN 14080 (2013) for GLT and EN 16351 (2015) for CLT. The entire system of strength grading and loading protocols in Europe are focused in edgewise bending or tension tests although boards used in engineered wood products are stressed in flatwise bending (Stapel and van de Kuilen 2014). This implies the need for research to verify the applicability of flatwise bending testing for the European standardization with C classes.
The motivation of this study was to determine the mechanical characterization of visually graded boards using flatwise bending tests to check whether the strength classes in EN 338 (2016b) obtained by edgewise bending tests are similar for Turkish fir and black pine timber. For the nondestructive testing, the stress wave method is investigated to evaluate the potential applicability and effectiveness for grading. Furthermore, a simplified nonlinear material model was proposed for the numerical modelling of boards in structural applications like flooring or a layer definition of CLT panels and GLT beams.
In this study, flatwise bending tests are arranged for 400 specimens with 22 mm ×50 mm ×420 mm (depth·with ·length of boards) including 200 fir (Abies nordmanniana subsp. Bornmuelleriana) from Bolu region and 200 black pine (Pinus nigra Arnold. subsp. pallasiana) from Balıkesir region from Turkey as described in Table 1. The length of specimens is selected as 19 times the depth of boards, which is shorter than edgewise bending tests. The moisture contents of boards were measured using the oven dry method by following the procedure defined in EN 13183-1 (2002). The average moisture contents were found as 11,4 % both fir and black pine.
In the outline of the research, specimens are first visually sorted as shown in Figure 1. Then, boards are tested by the stress wave method due to its easy applicability and the lack of experimental data in the literature. Next, boards are destructed using a flat wise bending test setup. Afterwards, a numerical model is proposed. Finally, the destructive and nondestructive results are analyzed statically for the comparison and discussed in detail.
All specimens used in this study were classified according to TS 1265 (2012). For structural usage, softwood timber products were sorted into three grades: class 1 with high load-bearing capacity, class 2 with moderate load-bearing capacity, and class 3 with low load-bearing capacity. On the other hand, the boards that did not match any of the grading requirements were referred to as rejected. Although the knot area ratio (KAR) is used for visual grading before the lamination process in the standard of ASTM D3737 (2018), TS 1265 (2012) requires calculating the knot diameter ratio (KDR) for each board specimen. KDR is defined as the knot diameters divided by twice the width of the board. Therefore, knot dimensions were visually checked and measured through the length of specimens. As shown in Figure 2a, knot diameter is measured by the surface which is parallel to the board. Moreover, the knot clusters are calculated by the sum of the KDRs of the knots existing in any 150 mm length of a piece of boards as shown in Figure 2b. Visual grading criteria for structural timber boards were defined in Table 2 according to TS 1265 (2012).
Moreover, the slope of grain, annual ring width and other defects were measured as specified in EN 1309-3 (2018b). Then, all specimens were cut according to the location of critical knots to be placed in between the loading points. Apart from visual grading, the angle between the flat surface of the samples and annual rings, αwas also measured for each specimen. Then, specimens were grouped as 0°-30° (tangential),30°-60° (rift) and 60°-90° (radial) as shown in Figure 3 according to EN 844 (2019).
Following the visual grading, the dynamic modulus of elasticity (MOEd) was determined for each specimen by using stress wave method with a portable microsecond stress wave timer. Fakopp Microsecond Timer (23 kHz) has two piezoelectric-type transducers with 60 mm long spikes as shown in Figure 4. For the application, the source is selected as a simple hammer impact and then the time of flight is measured. Then, stress wave velocity (m/s) was calculated using the distance between the transducers (L, m) and the time of flight taken from the device (∆t, μs) by the following equation (Equation 1). After that, dynamic modulus of elasticity was calculated with two different methods. MOEd1 (MPa) included the whole specimen’s density (ρ1, kg/m3) as a standard method but the density of defect-free piece was proposed by MOEd2 (MPa) which is cut from the same specimen (ρ2, kg/m3) by the following equations (Equation 2 and Equation 3). Meanwhile, whole specimens’ density (ρ1) was calculated by volume and mass however defect-free density (ρ2) was calculated after the destructive tests.
(1)
(2)
(3)After the nondestructive tests, 400 boards with 22 mm ×50 mm ×420 mm were tested flatwise in a four-point bending test setup as stipulated in EN 408 (2012a) with a 396 mm span length. Tests were carried out with the universal testing machine, which was equipped with a load cell of 100 kN. All specimens were tested under flatwise bending, and the deformations (w) were measured at the center of the span. Thereby, global modulus of elasticity (MOE) and modulus of rupture (MOR) were determined by following Equation 4 and Equation 5.
(4)
(5)Where: MOE is the global modulus of elasticity (MPa), MOR is the modulus of rupture (MPa), l is the span distance (mm), a is the distance between a loading position and the nearest support in a bending test (mm), F2-F1 is the load difference at 10 % and 40 % of maximum load (N), w2-w1 is the deflection difference at 10 % and 40 % of maximum load (mm), G is the shear modulus assumed as 650 MPa, b is the with (mm) and h is the height of samples (mm). The defect-free density (d2) of each specimen at the moisture content during the test was determined from a piece of the board cut as close as the location of the fracture after the experiments. Then, the pieces were immediately weighed with an accuracy of 0,01 g and their dimensions were measured with an accuracy of 0,01 mm. Furthermore, the moisture contents were found with the same samples by following the procedure defined in EN 13183-1 (2002).
According to EN 384 (2018a), there are several adjustments required for obtaining characteristic values:
The experimental values for the modulus of elasticity and the density of specimens that were not at the reference moisture content, were adjusted by using the following formulas (Equation 6 and Equation 7):
(6)
(7)Where: u is the moisture content at testing (8 % ≤ u ≤ 18 %), and uref is the reference moisture content (normally uref = 12 %).
MOR shall be adjusted to 150 mm depth by dividing with the factor kh as described in the formula. Therefore, the bending strength values obtained from the samples with a depth of nominal 22 mm were adjusted to 150 mm depth by dividing by the factor kh (Equation 8).
(8)MOE shall be adjusted to the modulus of elasticity E0 by using the following formula (Equation 9):
(9)After completing the required adjustments, the 5-percentile strength values f05,i, defect-free density d05,i and the mean stiffness values were determined for each grade of fir and black pine species as stimulated in EN 14358 (2016a) where bending strength and density were assumed as logarithmically distributed and modulus of elasticity was assumed as normally distributed. Thus, the parametric method was used to calculate the 5-percentile values of bending strength and density. The coefficient of subsample was neglected because specimens were delivered from one region only.
Numerical studies focused previously on the mechanical characterization of timber such as bending strength and stiffness (Berg et al. 2019, Olsson et al. 2013), tensile strength and stiffness (Mitsuhashi et al. 2008, Kohler et al. 2013). In this study, a simplified nonlinear material model was adopted to capture both elastic and plastic zones for modelling the out-of-plane bending behavior of boards in structural applications. Hence, the material models of fir and black pine species were generated for each grade using the stress distribution as shown in Figure 5. The distinct characters of tension parallel-to-grain and compression parallel-to-grain of timber under bending were expressed separately in the material model. For each visual grade, ultimate tension stress (ft,0,u) of boards was taken as equal to the mean values of ultimate bending stress fm according to Nwokoye (1972). Bending behavior of boards in tension stress was defined as linear brittle material where the mean values of the elastic modulus Emean were obtained from the experimental tests. However, bending behavior in compression stress was defined as ideal elasto-plastic material using Emean and the ultimate stress of compression parallel-to-grain (fc,0,u) was assumed to be reached at 0,003 mm/mm based on (Frese et al. 2012). In the numerical models, the support conditions of a simple beam and vertical loads according to four-point bending were firstly defined in SAP 2000 software (2000) as described in the test setup as shown in Figure 6. Then, boards were created using 2D layered/nonlinear shell elements with the simplified nonlinear material models for each grade and species. Finally, the vertical loads were increased step by step to find out peak displacements and load-carrying capacities of boards under bending. The tensile and compressive stress distributions can be also seen in Figure 6.
The destructive and nondestructive results were used for the analysis of variance (ANOVA) to detect whether any significant differences exist statistically within density values, MOR, MOE, velocity, MOEd1 and MOEd2 considering the visual grades of the specimens based on TS 1265 (2012). Average values for moisture content, ρ1, ρ2, MOR, MOE, velocity and MOEd1 and MOEd2 were given separately for each visual grade in Table 3 and Table 4, for fir and black pine species respectively.
Furthermore, the standard deviation values were given in parentheses. Afterwards, Duncan test was applied to the variables that had significant differences and grouped as shown with lower-case letters over the mean values.
In the visual grades, ANOVA was performed for each variable to determine whether there was any difference as equal to a significance level of p<0,05. According to Table 3 and Table 4, a significant difference was not found in ρ2 as expected for the defect-free pieces compared to ρ1for fir and black pine species, respectively. However, significant differences were found in MOR and MOE. For both species, the mean values of MOR and MOE decreased in lower visual classes as it was anticipated for visual grading.
Using stress wave method, three groups were found between velocity values as class 1, class 2-3 and rejected for both species. Due to the location and large size of edge knots as shown in Figure 7, the prediction of velocity varied significantly in board specimens other than class 1 that may lead to unconservative results. According to mean values of MOE, class 1 and 2 were higher than mean values of dynamic modulus of elasticity but class 3 was almost equal for both species.
In the studies of Van Duong and Matsumura (2018) and Van Duong and Ridley-Ellis (2021), however, the modulus of elasticity values obtained by three-point bending tests were about 15 % lower than the dynamic modulus of elasticity values measured by the stress wave method for small clear wood samples with 20 mm ×20 mm ×320 mm length. The difference between MOE and MOEd could be attributed to tested specimens with defects and the four-point bending test setup used in this paper.
For each visual grade, the difference between MOEd1 and MOEd2 was resulted from the density of the defect-free pieces. The mean values of MOEd1 were higher than the mean values of MOEd2 which was caused by the increased weight with higher density of knots especially for class 3 and rejected samples. Hence, MOEd2 groups went parallel with velocity remarkably for fir and black pine and more effective to detect rejected samples. For fir, the ratios between MOE and MOEd2were 0,79 for class 1; 0,85 for class 2; 0,92 for class 3 and 1,24 rejected samples, respectively.
On the other hand, the ratios between MOEd1and MOE were 0,81 for class 1; 0,88 for class 2;0,99 for class 3 and 1,39 for rejected specimens in fir, respectively. In black pine, the ratios between MOEd2 and MOE were 0,82 for class 1; 0,89 for class 2; 0,96 for class 3 and 1,14 for rejected samples, respectively where the ratios between MOE and MOEd1were 0,84 for class 1; 0,93 for class 2 and 1,01 for class 3 and for 1,26 rejected samples, respectively. Please note that MOEd1 and MOEd2 overestimated the modulus of elasticity of rejected samples while defect-free density performed better compared to whole specimens’ density.
In addition, no significant difference was found for the angle, α in destructive and nondestructive results according to Duncan tests according to Table 5. Therefore, the mechanical properties of boards were not affected remarkably by the angle, α under bending.
The frequency histograms of specimens and regression matrix of mechanical properties were plotted using MOR-MOE, MOR-knot rate, MOE-knot rate, MOR-velocity, MOE-velocity and knot rate-velocity. Linear regression matrices for MOR, MOE, knot rate and velocity were shown in Figure 8 and Figure 9 for fir and black pine respectively. In the matrix, histograms were located at the points where the variables coincided with each other. For fir, there were very strong correlations in MOR-MOE (R2=0,837), MOR-Knot rate (R2=0,669) and MOE-Knot rate (R2=0,616).
Furthermore, strong correlations were found in MOR-Velocity (R2=0,430), MOE-Velocity (R2=0,554) and a moderate correlation was found in Velocity-Knot rate (R2=0,307). For black pine, there were very strong correlations in MOR-MOE (R2=0,776), MOR-Knot rate (R2=0,660) and MOE-Velocity (R2=0,700). Similarly, strong correlations were found in MOR - Velocity (R2=0,489) and MOE-Knot rate (R2=0,474). In addition, there was a moderate correlation in Velocity-Knot rate (R2=0,253). For clarity and better understanding, the relations of MOEd1 and MOEd2 were presented separately using correlation matrix.
For sub-groups of visual grading, correlation matrices for MOR, MOE, knot rate and velocity as well as MOEd1 and MOEd2 were given in Table 6 and Table 7 for fir and black pine, respectively. For fir specimens, Table 6 shows that by using proposed defect-free density instead of the density of whole specimen, the overall correlations of MOR-MOEd2 and MOE-MOEd2 increased by 90 % and 56 % compared to MOR-MOEd1 and MOE-MOEd1, respectively. Likewise, the overall correlations of MOR-MOEd2 and MOE-MOEd2 increased by 46 % and 19 % compared to MOEd1 for black pine, respectively. More over, the overall correlation of velocity-MOEd2 increased by 25 % for fir and 9 % for black pine compared to MOEd1, respectively.
The strength class CXX defines the bending strength of edgewise bended samples in the strength grading of Europe whereas specimens were intentionally tested flat wise in this paper to research the potential for the applicability. After the 5 % percentile values of bending strength, modulus of elasticity and density were calculated, strength classes were declared in accordance with CEN EN 338 (2016b) for fir in Table 8 and black pine in Table 9.
Although the characteristic bending strength of fir and black pine were respectively determined as 39 MPa and 43,2 MPa in class 1, lower strength classes are attained due to the low density and modulus of elasticity for both species: C24 for fir and C40 for black pine. Similarly, boards with class 2 were attained as C22 and C27 whereas class 3 boards were attained as C18 and C22 for fir and black pine, respectively. On the other hand, rejected grades could not be included in any strength class due to the limited capacity of bending strength and modulus of elasticity for both species.
To compare the strength classes obtained by flatwise bending with edgewise bending, Anatolian black pine and Caucasian and taurus fir species of Turkey presented in CEN EN 1912 (2012b) were used. Higher resistances were found for black pine boards (Class 1= C40, Class 2= C27 and Class 3= C22) in flatwise bending tests compared to strength classes using edgewise bending tests where Class 1= C35, Class 2= C24 and Class 3= C18. Similarly, same resistances were found for fir boards in Class 1= C24 and Class 3= C18 however a small decrease is found for Class 2 as equal to C22 in flatwise bending compared to C24 in edgewise bending tests.
The failure mode of boards can be summarized generally in three stages: (I) the stresses of beam section were in the linear elastic zone with low deflections on the middle span, and (II) linear elastic point was exceeded in the load-deflection curve and plastic deformations increased with the compression strains started to accumulate in the upper side of the board close to the loading head, then (III) the ultimate deflection was reached with brittle tension rupture of timber grains at the bottom of the board due to high tensile demands under bending as shown in Figure 10. However, the tension rupture of timber grains occurred earlier in grades with high KDR values. Therefore, the second stage was developed rapidly with very limited plasticity in class 3 and rejected specimens. For instance, the mean ultimate deformation capacity of boards in rejected class was found lower than other strength classes for both species.
In addition, the transition point that defines the start of plastic region and softening was found for each specimen by monitoring the change in the slope of the load-deflection curves with a correlation coefficient of 0,99 after exceeding linear load limit beyond 40 % of the peak load. The mean values of force, Ft and deflection, Dt were given in Table 10 at the transition point for fir and black pine grades. While class number gets better, the mean values of force and deformation were increased in both species.
The experimental load-deflection curves of specimens were displayed for fir and black pine species respectively in Figure 11 and Figure 12 with the numerical models generated for each visual grade. From the comparison of load-deformation curves and numerical models, the simplified nonlinear material model was able to capture sufficiently not only the linear behavior but also the plastic behavior with a sudden drop in strength after the tensile rupture of timber parallel-to-grain under bending.
After the experimental tests, the mean values of the ultimate deflection, Dmax,m at the mid-span when 80 % of the ultimate load was reached at the descending part of the load-deformation curve and the mean values of the ultimate load, Fmax,m obtained in tests were also calculated and given in Table 11 for fir and black pine visual grades to compare the accuracy of the numerical modelling. Table 11 illustrates that the numerical model estimated the bending strength accurately with a difference of 0,24 % and 2,77 % respectively for class 1 and class 2 but higher differences were found for class 3 and rejected for fir grades.
In black pine grades, the numerical model caught the bending strength with a maximum difference of 3,97 % in visual classes except for the rejected specimens. Although the deformation capacities were estimated with high accuracy for class 1 and class 2, the simplified material model overestimated deformation capacities of the class 3 for both species. As a result of the high KDR values, the failure mode tended to be more brittle compared to specimens with fewer defects leading unconservative deflection results in class 3. Finally, the proposed material model should not be utilized for the rejected specimens which are forbidden to be used in structural elements.
In conclusion, 400 board specimens from fir and black pine species were visually graded according to TS 1265 (2012) and stress wave method was applied for nondestructive tests. Then, boards were tested under flatwise bending and load-deflection relationships were given detailly for each grade. After the statistical analysis, the mean values of MOR and MOE differed to three visual classes and rejected at a p<0,05 significance level.
Hence, high correlations between modulus of rupture and knot ratio indicated that limiting parameters of visual grading standard were suitable for boards.
To optimize visual grading, KDR calculation method in the existing standard could be improved especially for edge knots. Further more, the mechanical properties of boards under bending were not affected statistically by the angle between the flat surface of the samples and annual rings.
After comparing the visual grading and the stress wave method, it was discovered that visual grading was more reliant to sort board specimens into classes since the difference in the mean values of MOE between visual grades were more remarkable than the mean values of MOEd. For the determination of dynamic modulus of elasticity, however, the usage of the defect-free density performed better than whole specimen’s density specially to detect rejected samples. Overall, flexural stiffness obtained by destructive tests were 14,6 % for fir 10,9 % for black pine higher than nondestructive test results except for the rejected samples. The difference could be attributed to the small length of tested specimens and the transverse shear effects in the calculation of the global modulus of elasticity as stimulated by the latest version of CEN EN 408 (2012a)
The strength classes of boards were determined using destructive test results according to CEN EN 338 (2016b) with the required adjustments. The strength class assignments using flatwise bending were presented only for demonstration which should be considered carefully.
Since the specimens were taken from one region Bolu for fir and Balıkesir for black pine, no adjustment factors for the sampling were applied. Higher resistances were found for black pine boards in flatwise bending tests compared to edgewise bending classes defined in CEN EN 1912 (2012b). However, same resistances are found for fir boards in class 1 and class 3 except for class 2 with a small decrease. Therefore, more tests using flatwise bending are required to apply for the standardized system as boards used in engineered wood products are stressed in flatwise.
Lastly, the failure mode of boards under bending was observed by the rupture of timber fibers parallel-to-grain due to high tensile stresses for both species. With the proposal of a simplified nonlinear material model, experimental and numerical model results were found in good agreement in terms of bending strength and deformation capacity especially for class 1 and class 2 to model boards in structural applications. For class 3, the usage of linear material models should be appealed to be safer and more conservative in terms of the deformation capacity. For future studies, an experimental campaign is planned to produce and test cross-laminated timber panels by using the structural timber boards from Turkish fir and black pine.
♠Corresponding author: fatihkurul@istanbul.edu.tr
