1. Introduction
Wood is a viable material that is suitable for use in different applications. In civil construction, it has been highly recommended worldwide because of its sustainability since it is a renewable natural resource with low energy consumption and a great contribution to reducing greenhouse gases [
1,
2].
Currently, there are structural calculation programs that allow for precise dimensioning of structures, even the most complex ones. However, the accuracy of these tools depends on the insertion of data in a complete and coherent way with the type of material (isotropic, orthotropic, and transversely isotropic), making the acquisition of these data essential. In addition to civil construction, wood has great importance in the manufacture of musical instruments, which is another example of an area for which knowledge of elastic properties is essential [
3,
4,
5,
6,
7,
8,
9]. Nevertheless, as wood is a complex material due to its anisotropic and heterogeneous characteristics, obtaining knowledge of its mechanical properties is still a challenge [
10,
11].
By reducing the dimensions of specimens, it is possible to obtain more rectilinear growth rings and, thus, get closer to the condition of orthotropy with three main axes of symmetry (longitudinal—L, radial—R, and tangential—T) that respond differently to the action of loads. The L axis is parallel to the direction of the grain (fiber or tracheid); the R axis is perpendicular to the grain in the direction of the radius; and the T axis is perpendicular to the grain but tangents the growth rings. Thus, for the complete characterization of this material, one must determine twelve elastic constants (nine of them independent): three longitudinal moduli of elasticity (EL, ER, ET), three shear moduli of elasticity (GLR, GLT, GRT), and six Poisson ratios (μLR, μRL, μRT, μTR, μLT, μTL).
Ultrasound technology allows one to determine, in a non-destructive way, all the elastic constants of orthotropic materials, such as wood, and its theoretical framework is discussed by Bucur [
12], Musgrave [
13], and Hearmon [
14]. Over the years, with the development of equipment, several studies have been conducted corroborating the effectiveness and accuracy of ultrasound technology in the characterization and inference of wood properties [
15,
16,
17,
18,
19,
20,
21,
22,
23]. For the complete characterization of wood with ultrasound, it is possible to use several specimen formats, such as spherical [
10,
24], prismatic [
3,
20,
25,
26,
27], polyhedral [
12,
20,
22], and multifaceted discs [
12,
20]. Polyhedral specimens with 26 faces were initially proposed by Francóis [
28] and polyhedral specimens with 18 faces were proposed by Arts [
28,
29].
In the process of the complete characterization of wood by ultrasound, some authors report difficulties in determining the terms outside the main diagonal of the stiffness matrix [
12,
20] since a distortion of the wavefront, caused by the inclination of the wood growth rings, can be observed. According to Bucur and Archer [
27], reducing the size of the specimen favors the approximation of orthotropy since the curvature of the growth rings can be neglected, while in larger sections, the curvatures of the growth rings are significant [
30]. On the other hand, the reduction in the specimen must be compatible with the maintenance of plane wave propagation in infinite media, and for this, one must make the reduction in dimensions compatible with the frequency of the transducer used [
12]. Considering the identification of interferences in the inclination of the growth rings in the results of ultrasound characterization and the greater feasibility of obtaining specimens with more rectilinear growth rings, reducing the size of the specimen, this research aimed to verify whether the reduction in the polyhedral specimen compared to the traditionally used one (26 faces starting from a 70 mm edge cube) changes the results of the elastic parameters obtained in the characterization of wood by ultrasound. To achieve the objectives, the experimental design involved the determination of the elastic parameters of three wood species using polyhedral specimens of 26 faces and dimensions traditionally used in the literature (starting from 70 mm edge cubes), which were considered as references, as well as polyhedra of 26 faces of reduced dimension and polyhedra of 18 faces, both starting from edge cubes with approximately half the traditional polyhedron dimension (≅35 mm). To minimize interference caused by other parameters, dimension reductions were made to continue guaranteeing the theoretical framework of wave propagation in infinite media.
Thus, the main novelty of this study is the analysis of the behavior of the results of the elastic characterization of wood by ultrasound using specimens in a polyhedral format with reduced dimensions and fewer faces than those traditionally used. This analysis offers a safer basis for the modification of a traditionally used and known specimen and, as a consequence, reduces the amount of material used in the manufacture of specimens and simplifies the process of producing polyhedra. The simplification of the process, in turn, allows for an increase in the accuracy of the specimen’s production with better-directed growth rings.
2. Materials and Methods
The experimental scheme of the research, presented in detail in this topic, is shown in
Figure 1.
2.1. Material
Three wood species were selected for the tests: Apuleia leiocarpa, Dinizia excelsa, and Cedrelinga catenaeformis. The choice of species was made to cover a range of densities that would allow for greater generalization of the results.
For each species, three pieces of 0.10 × 0.20 × 1.0 m³ were obtained, which were purchased with moisture already around 12% (
Table 1). The pieces had their weight and volume measured for the calculation of densities (
Table 1) and monitoring of moisture by contact measurement equipment (Merlin, PM1-E, Tumeltsham, Austria).
For wood, a material with great variability within the same species and even within the same tree, comparisons with samples from other authors that aim at the validation of tests are made only based on the order of magnitude of the results. According to the Institute of Technological Research [
31],
Apuleia leiocarpa and
Dinizia excelsa have average apparent densities (with 12% moisture) of 830 kg·m
−3 and 1090 kg·m
−3, respectively, which are values with mean differences of around 2.5% from those obtained in this study. In the case of
Cedrelinga catenaeformis, the difference was 22% compared to that obtained by Dias and Lahr [
32], of 566 kg·m
−3. This species was also the one with the highest density variability between the pieces (
Table 1). Despite the greater difference obtained for
Cedrelinga catenaeformis, the order of magnitude of the densities obtained in this research is compatible with the one obtained by other authors.
2.2. Acquisition of Polyhedral Specimens
The polyhedra that composed the research sample and that were used in the ultrasound tests were obtained from the 0.10 × 0.20 × 1.0 m
3 pieces (
Table 1). For each of the species, 6 polyhedra of each shape and dimension adopted in the research were made (
Figure 2), totaling 54 polyhedra, including 18 polyhedra with 26 faces (P26N) with dimensions traditionally used in the literature, 18 polyhedra with 26 faces of reduced dimension (P26A), and 18 polyhedra with 18 faces (P18).
The polyhedra were taken from cubes oriented in the three main directions of the wood (L or 1, R or 2, and T or 3), with the growth rings as rectilinear and parallel as possible in relation to the edge of the cube. The polyhedron with 26 faces of traditionally used dimensions [
17,
18,
19,
20,
22] was obtained from a 70 mm edge cube (
Figure 3c). The polyhedron with 26 faces and reduced dimensions was made from a 35 mm edge cube (
Figure 3a), half the size of the cube edge (70 mm) that originated the traditionally used polyhedron (
Figure 3c). The dimension of the cube edge that originated the polyhedron with 18 faces (
Figure 3b) was adopted to maintain, as the minimum dimension of the faces, the diameter of the 1.0 MHz frequency transducer (12 mm), ensuring that the transducer could be circumscribed to the face of the polyhedron, avoiding the wall effect during wave propagation [
12]. The 1.0 MHz frequency transducer was adopted as a reference for the proposal of the dimensions of the specimens because this frequency level is considered the maximum that allows for the maintenance of a wavelength compatible with the dimensions of the anatomical elements of the wood, minimizing signal attenuation problems [
12].
2.3. Measurement of the Growth Ring Inclination on the Face of the Polyhedron
Considering that the reduction in the dimension of the polyhedra is proposed to test the hypothesis that it allows for a reduction in the inclination of the growth rings compared to the radial or tangential face and that this reduction minimizes wave dispersion effects, improving the purest propagation in the axes or planes, the inclination of the rings was measured. To measure the inclination of the growth rings, it was first necessary to sand the cross section (RT) of each polyhedron until it was possible to clearly identify the rings, using sandpaper weights of 120 and 400. Initially, the samples were sanded in an orbital sander (DeWalt, DWE6411, Uberaba-MG, Brazil), and then, to allow a finer finish, a polisher (Panambra, DP-10, São Paulo-SP, Brazil) was used. The cross section (RT) of the polyhedra was then photographed using a camera (NIKON D7000 with AF-S DX VR Zoom-Nikkor 18–105 mm ƒ 3.5–5.6 lens) fixed on a tripod. The samples were positioned so that the images had the radial and tangential directions parallel to the Y and X axes (
Figure 4). The images were then opened in an app (Paint) where the direction of the growth ring (blue line in
Figure 4) and the radial (R) and tangential (T) directions (red lines in
Figure 4) were demarcated, the red line crossing the cross section being indicative of the direction of propagation of the ultrasound wave (
Figure 4). The lines (blue and red) were used to determine the angle α (
Figure 4) between the direction of the growth rings and the direction of propagation of the ultrasound wave using the ImajeJ software (1.53t).
2.4. Ultrasound Tests on Polyhedra
Tests on polyhedral specimens were performed with ultrasound equipment (EP1000, Olympus, Waltham, MA, USA) with 1 MHz frequency compression and shear wave transducers (
Figure 5). The coupling between the transducer and the material to be tested is a factor of great importance for the accuracy of the results, and the use of viscous material is indicated to eliminate the air between the transducers and the surface of the evaluated specimen. The coupler also has the function of reducing the impedance difference in the wave passage to the material under inspection, minimizing the attenuation of the signal as a function of wave deviation (Snell’s law); this research adopted starch glucose, which, according to Gonçalves et al. [
17], presents good results for the shear wave in wood.
2.5. Calculation of Wave Propagation Velocities
For the calculation of longitudinal velocities (V
11, V
22, and V
33), Equation (1) was used, in which the path length corresponds to the distance between the faces of the transducers and the wave propagation times in directions 1 (longitudinal), 2 (radial), and 3 (tangential), respectively, are obtained using the longitudinal transducer.
where V = velocity of wave propagation in a given direction (m·s
−1); L = wave path length (m); t = wave propagation time in a given direction (µs).
Considering the same directions but with the use of the shear transducer, the transversal velocities (V12, V13, V21, V31, V32, and V23) were calculated by obtaining the wave path time (t12, t13, t21, t31, t32, and t23) with propagation in a certain direction (first numerical index) and its polarization in the perpendicular direction (second numerical index).
For the determination of the velocities corresponding to the wave propagation outside the axes of symmetry, the wave propagation time was obtained using shear transducers, but on the faces representing axes inclined 45° with respect to each plane.
2.6. Calculation of Stiffness [C] and Flexibility [S] Matrices
With the velocity data, the stiffness matrix [C] (Equation (2)) was determined using Christoffel’s equations (Equations (3)–(6)). The density values in the equations were those obtained for the pieces used to make the polyhedra (
Table 1). The stiffness coefficients of the main diagonal (C
11, C
22, C
33, C
44, C
55, and C
66) were obtained by the Christoffel equation (Equation (3)) with the longitudinal velocities (V
11, V
22, and V
33) propagating and polarizing in the direction of the main axes of symmetry (longitudinal, radial, and tangential:
Figure 6). To calculate coefficients C
44, C
55, and C
66, the velocities whose wave propagation times on the axes were obtained with the shear transducer were used, with wave propagation on one axis and its polarization on the perpendicular axis, with the numbering related to the axes (propagation/polarization) having the following nomenclature: 44 = propagation on axis 2 and polarization on axis 3 or propagation on axis 3 and polarization on axis 2 (plane 2–3); 55 = propagation on axis 1 and polarization on axis 3 or propagation on axis 3 and polarization on axis 1 (plane 1–3); and 66 = propagation on axis 1 and polarization on axis 2 or propagation on axis 2 and polarization on axis 1 (plane 1–2). In the same way (by solving the Christoffel equation), the 3 terms outside the diagonal (C
12, C
13, and C
23) can be obtained. For this, the propagation must occur inclined to the axes of symmetry (45° in this research). The general formulas, deduced from the Christoffel tensor for determining the constants in planes 1–2 (LR), 2–3 (RT), and 1–3 (LT), are given by Equations (4)–(6).
where i = 1, 2, 3, 4, 5, and 6; ρ = density of the material (kg·m
−3); V = velocity of wave propagation in the considered direction (m·s
−1).
where α = angle (45°); n
1 = cos α, n
2 = sin α, and n
3 = 0 (α is taken with respect to axis 1) (plane 1–2) for Equation (4); n
1 = cos α, n
3 = sin α, and n
2 = 0 (α is taken with respect to axis 1) (plane 1–3) for Equation (5); n
2 = cos α, n
3 = sin α, and n
1 = 0 (α is taken with respect to axis 2) (plane 2–3) for Equation (6).
By inverting the matrix [C]
−1, it was possible to obtain the flexibility matrix [S] (Equation (7)), which is associated with the elastic parameters of the material (longitudinal moduli of elasticity—E, shear moduli of elasticity—G, and Poisson’s ratio—υ).
2.7. Verification of Basic Aspects of Wave Propagation in Infinite Media
As we used polyhedral specimens of reduced dimensions (18 faces and 26 faces reduced), it is important to analyze whether the theoretical aspects regarding wave propagation in infinite media are being met. For this propagation to occur in infinite media, the distance between transducers (dimension of the piece in the direction of wave propagation) must be a few times greater than the wavelength (λ). For wood, some researchers [
12,
33,
34] indicated values from 2 λ to 5 λ for this condition to be met. The minimum dimension of the piece is linked to the frequency of the transducer (f) and to the propagation velocity of the wave in the material (V) since λ is given by the ratio between these two parameters (λ = V/f). When this theoretical condition is violated, the results of the elastic parameters obtained by ultrasound can be affected [
12,
26,
34].
Because of the attenuation of the wave in the wood, the effective frequency (FEF) is different from the nominal frequency of the transducer (f), which, in the case of this research, is 1 MHz. To obtain the effective frequency (FEF), one must obtain the time difference (Δt) between two successive peaks of the wave (FEF = 1/Δt) and therefore the signal must be evaluated (example in
Figure 7). Considering the frequency of the transducer and the dimensions of the specimens adopted in this research, it was found that, even in the most critical cases (shorter path length and higher velocity), at least two complete waves passed through the material (
Table 2), thus meeting the theoretical conditions of infinite media.
2.8. Analysis of Results
The results of the inclination of the growth rings, wave propagation velocities (longitudinal and shear) in different directions and planes, and elastic parameters were initially evaluated for normality to ensure the feasibility of using parametric statistics. The multiple range test was used for the comparisons of the results for the traditional polyhedra with 26 faces, the reduced polyhedra with 26 faces, and the polyhedra with 18 faces.