ARTIFICIAL NEURAL NETWORKS TO ESTIMATE THE PHYSICAL-MECHANICAL PROPERTIES OF AMAZON SECOND CUTTING CYCLE WOOD

Timber from the second cutting cycle may make up the majority of future crop volumetric. However, there are few studies of the physical and mechanical properties of this timber, which are important to support the consolidation of new species. This study aimed to use Artificial Neural Networks to estimate the physical and mechanical properties of wood from the Amazon, based on basic density. The properties were: shrinkage (tangential, radial and volumetric), static bending, parallel and perpendicular to the fiber compression, parallel and transverse to the fibers, Janka hardness, traction, splitting and shear. The estimate followed the tendency of the data observed for the tangential, radial and volumetric shrinkage. The network estimated the mechanical properties with significant accuracy. Distribution of errors, static bending, parallel compression and perpendicular to the fiber compression also showed significant accuracy. Artificial Neural Networks can be used to estimate the physical and mechanical properties of wood from Amazon species.


INTRODUCTION
The period between one selective logging harvest and another is called the cutting cycle.During this period, forests have time to regrow harvested stock.These recovery periods are determined by Brazilian legislation (Brasil 2006).
Many studies demonstrated the feasibility of the second cutting cycle in relation to carbon stock recovery.They also showed that the first harvest helps the recovery of the forest volume and basal area, boosting the growth of the remaining forest, as reported in studies by Reis et al. (2010), Reis et al. (2014) and Reis et al. (2015).
The Tapajós National Forest completed its second monitored cutting cycle in 2014.It is the first forest in the world to achieve this outcome, demonstrating that selective wood harvesting in native forests is sustainable.According to Reis et al. (2010), over 28 years the forest managed to recover its total volumetric stock.This showed that the second cutting cycle in this experimental location is possible, mainly due to the dominant stock of species that were not used in the first harvest.Second cycle wood comes from species that were favored by forest harvesting and increased their density in the area.Even with little commercial attractiveness, it may have a significant market, increasing manager production and subsequently, income.
According to Adeodato et al. (2011), studies have indicated that the commercialization of second line wood (with only potential or no consolidated markets) may increase revenue per hectare by 40% to 50% in forestry management.There is a great effort from research organizations to analyze the physical characteristics of wood whose species are unknown in the market or which have low commercial value, to replace those with higher commercial value (Almeida et al. 2010).
In order to promote new species, we require a better understanding of the physical-mechanical properties of wood, to back its introduction into the market.The determination of physical and mechanical wood properties is important to decide on its correct use.This will increase the number of available species and reduce pressure on the most common commercial species.It also helps in forestry management, since the growth in the number of potential species makes it possible to increase production per area unit and to adopt more harvesting options.
To estimate these properties, Artificial Neural Networks (ANN) have been used.These are computational models inspired by the nervous system of living beings.They can be defined as a massive set of parallel processing units, characterized by artificial neurons linked by a great number of interconnections.ANN's have a natural ability to store experimental information making it available for subsequent use (Silva et al. 2010).Various studies have been conducted using ANN's in forest areas, providing successful estimates, in both even and uneven-aged populations (Ashraf et al. 2015, Diamantopoulou et al. 2015, Reis et al. 2016, Richards et al. 2008) as well as in wood technology areas (Avramidis et al. 2006, Avramidis and Iliadis 2005, Esteban et al. 2009, Tiryaki and Aydın 2014) ANN's show advantages in terms of regression analyses, due to their inherent properties and capacities, such as: non-linearity, learning capacity, generalization abilities, adaptability, fault tolerance and data organization (Silva et al. 2010).Since they present data such as the physical and mechanical wood properties, which have many intra and interspecies variations, ANN's may improve the adjustments to estimate these properties.
The present study seeks to use Artificial Neural Networks to estimate the physical and mechanical properties of second cutting cycle timber from the Amazon starting from basic density.
In 1975, a forest census was realized and in 1979, 64 wood species were harvested, removing 72.5 m 3 ha -1 .In 1981, 36 permanent plots (50x50m) were randomly positioned in the study area.In these plots, all dbh ≥ 5 cm trees were botanically identified in loco and samples were collected for subsequent identification at the IAN Herbarium of the Embrapa Eastern Amazon (Reis et al. 2010).
In 2014, the harvesting area was divided into two treatments: the first cycle area where harvesting had not taken place had 39 ha (T1), while the second cycle area where harvesting had occurred in 1979 had 31,5 ha (T2).The two areas totaled 70,5 ha.Thirty-five species were harvested from the two areas, with an average intensity of 3,76 trees ha -1 and 15,3324 m 3 ha -1 .
The physical and mechanical property data of the timber was obtained from studies about species Artificial neural networks to..: Reis et al. from the Amazon region.Sources were publications by the Forest Product Laboratory (FPL 2016).
In order to estimate wood properties, basic density (kg.m -3 ) was considered an independent variable.For data regarding wood shrinkage, literature meeting the Pan-American Commission for Technical Standards (Copant) was considered.To obtain mechanical properties, all references consulted used test samples with a 12% moisture content, also in agreement with COPANT standards.
Artificial Neural Networks (ANN) were used to estimate the physical and mechanical variables of timber with no added commercial value, belonging to the species list from the survey conducted in 2014 in Flona do Tapajós.
Data was divided into two groups, one for ANN training (85%) and the other for the generalization-test (15%); this was to evaluate the generalization capacity of the networks.The training group was not part of the ANN test group.
Three thousand nine hundred ANN's were trained using basic density as an input.3600 ANN's were trained with simple outputs and 300 more were trained with multiple outputs (Table 1).ANN training consists of the application of an ordained step set that aims to adjust the weights and thresholds of its neurons.Such adjustment processes, also known as learning algorithms, seek to fine tune the network so that its responses are closer to the desired values (Silva et al. 2010).
Training was based on the generalized delta rule, also known as the backpropagation algorithm, applied in multi-layer feed forward networks -Multilayer Perceptron (MLP).
The backpropagation algorithm is performed through successive applications of two specific phases.According to Silva et al. (2010) and Reis et al. (2018) the first phase to be applied is called "forward propagation", in which signals{x 1 , x 2 ,..., x n } of a sample of the training set are inserted into the network inputs and are propagated layer by layer until the respective outputs are produced.The application of this phase only aims to obtain the network responses, taking into account only current values of synaptic weights and thresholds of its neurons, which will remain unchanged during each execution.Subsequently, the responses produced by the outputs of the network are compared with the respective desired responses that are available.The respective deviations (errors) between the desired responses and those produced by the output neurons are then calculated.This will then be used to adjust the weights and thresholds of all their neurons.Because of these error values, the second phase of the backpropagation method, called "reverse propagation" (backward), is then applied.The changes (adjustments) of the synaptic weights and thresholds of all the neurons of the network are executed during this phase (Silva et al. 2010).
In summary, the successive applications of the forward and backward phases make the synaptic weights and thresholds of the neurons automatically adjust at each iteration, implying a gradual decrease in the sum of the errors produced by the responses of the network compared with the desired ones (Silva et al. 2010, Reis et al. 2018).
For training, the Intelligent Problem Solver (IPS) tool of the Statistica 13 software (StaSoft Inc 2016) was used to analyze the activation functions of the intermediate and output layers (Identity, Logistics, Hyperbolic Tangent and Exponential).The initial weights of the networks were randomly generated, and the stopping criterion of the algorithm occurred when the mean square error or the cross-entropy error began to increase.When this happened, training was interrupted.Network training continued up to 10000 cycles as long as the error was decreasing (Reis et al. 2018).
During training, only one hidden layer was used, and the interval of the neuron number in this layer was defined by the Fletcher-Gloss method, where the interval of the neuron number is defined according to the number of input and output variables.This was defined by the Equation (1): where n is the network input number, n 1 is the neuron quantity in the hidden layer and n 2 is the neuron quantity in the output layer.
The maximum neuron number was defined by the Fletcher-Gloss method and the minimum was a single neuron.This was done to analyze the generalization capacity of a simpler topology.The maximum defined by the method was to avoid excessive input data memorization (overfitting) or insufficient information extraction during training (underfitting).
In order to choose and compare ANN's, the correlation statistics between estimated values and bibliographic values and the root-mean-square error (RMSE) were used.The graphic analysis of percentage error dispersion (Error %) in relation to the observed values was also used (Equation 2): (2) where : estimated values by ANN; Y: bibliographic values.

The correlation indicates the intensity of the relationship between estimated and observed values.
The closer to 1, the higher the correlation between the variables.The correlation was calculated using the Equation (3): (3) where : estimated values by ANN; Y: observed values; S 2 : variance and Cov: covariance.As for the root-mean-square error (RMSE), the lower the estimates, the more reliable the training or generalization.In order to calculate RMSE for the comparison between ANN's, the Equation ( 4) was used:
The network estimate for tangential and volumetric shrinkage was not as precise as the tested mechanical properties.Observing the correlation matrix (Table 3), the correlation between basic density and radial shrinkage was not significant.
( ) Shrinkage results were not as precise as the tested mechanical properties, because shrinkage is strongly influenced by the anatomical structure of the wood.According to studies by Zhang and Zhong (1992) analyzing the relation between anatomical structures and physical-mechanical properties of wood, fiber diameter is the anatomical characteristic that most influences tangential and radial shrinkage, followed by vessel and fiber proportion, which influence the radial and tangential shrinkage, respectively.To improve the accuracy of the model in relation to the shrinkage of the wood, laboratory tests and inclusion of wood structure variables, such as vessel and fiber length, should be carried out.
Wood shrinkage is also affected by other factors, such as the microfibril angle of the cell wall central secondary layer (Chauhan andWalker 2006, Poubel et al. 2011) and the structural characteristics of wood that are resistant to shrinkage, such as ray quantity in the radial direction (Poubel et al. 2011).Despite this, the estimate could follow the tendency of the observed data for tangential, radial and volumetric shrinkage.
Regarding error distribution, we noticed that the tangential, radial and volumetric shrinkage tended to estimate lower errors for the highest values (Figure 1).The network precisely estimated wood mechanical properties.Estimated data followed the tendency of the observed data.It was also possible to observe the high correlation between basic density and mechanical properties, which presented significant values (Table 3).
For all mechanical properties analyzed, increased density values were mirrored by increased property values (Figure 2).This is caused by the high correlation (Table 3) between density and other mechanical properties.Armstrong et al. (1984) demonstrated that density as an independent variable explains mechanical properties (FEM, FRM and FParC) with a correlation coefficient above 0,70 for a 95% reliability level, in green and dry wood conditions.Araújo (2007) also obtained significant results between wood density and mechanical properties.Other studies such as those by Abruzzi et al. (2012) and Protásio et al. (2012) also observed a relationship between density and mechanical properties.Density is one of the most important properties for recommended wood usage (Leite et al. 2016) and(Peres et al. 2012).It is one of the properties that provides the most information regarding wood characteristics and is related to its resistance.
With the results obtained, it is possible to observe that all mechanical properties studied have a correlation with basic wood density.Mechanical wood properties depend on basic density (Lobão et al. 2004).
The variability of most mechanical wood properties may be estimated based on density variations.The properties of resistance to static flexion and its elasticity module are correlated with density, which in turn is correlated to cell dimension.The relationship between density and mechanical properties may be changed by the presence of extractives that are added to wood mass, mainly increasing wood resistance to axial shrinkage (Lobão et al. 2004, Panshin andZeeuw 1980).
Regarding error distribution, the properties: static flexion, fiber parallel compression and fiber perpendicular compression had lower errors and, consequently, higher precision (Figure 2).Beltrame et al. (2010) showed that mechanical properties define wood behavior when submitted to mechanical stress, allowing comparison to other woods with known properties.

Studies conducted by
Regarding static flexion, fiber parallel compression, fiber perpendicular compression and fiber parallel shearing, the correlation result was high and the error result low (Table 2), indicating that the mentioned properties may be estimated with artificial neural networks using basic density.According to Lisboa et al. (1993), these properties are extremely important to obtain information about possible final uses.
Artificial neural networks to..:Reis et al.

Figure 1 .
Figure 1.Basic wood density in relation to the estimated tangential, radial and volumetric shrinkage, and error% in relation to the estimated variables.

Figure 2 .
Figure 2. Modulus of rupture (MR) and modulus of elasticity (ME), estimated parallel and perpendicular to the fiber compression in relation to the basic wood density and error% in relation to the estimated variables.

Table 1 .
Input and output variables used for the training of the ANN's and numbers of trainings performed.

Table 2 .
Characteristics, input and output variables of trained and validated ANNs.
a Bd: