Next Article in Journal
Relative Importance of Landscape and Climate Factors to the Species Diversity of Plant Growth Forms along an East Asian Archipelago
Next Article in Special Issue
The Influence of High-Temperature and -Pressure Treatment on Physical Properties of Albizia falcataria Board
Previous Article in Journal
Modeling Potential Changes in Rare Species Habitat from Planned Timber Harvest in Minnesota, USA
Previous Article in Special Issue
Prediction of Mechanical Properties of Thermally Modified Wood Based on TSSA-BP Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Temperature and Process Duration as a Method for Predicting the Mechanical Properties of Thermally Modified Timber

1
Key Lab of State Forest and Grassland Administration on “Wood Quality Improvement & High Efficient Utilization”, School of Forestry & Landscape Architecture, Anhui Agricultural University, Hefei 230036, China
2
Department of Wood Science and Technology, Faculty of Technical Engineering, University of Bihać, 77000 Bihać, Bosnia and Herzegovina
3
Department of Wood Science, Biotechnical Faculty, University of Ljubljana, 1000 Ljubljana, Slovenia
4
Department of Natural Resources and Earth Science, Shahrekord University, Shahrekord 64165478, Iran
*
Author to whom correspondence should be addressed.
Forests 2022, 13(2), 217; https://doi.org/10.3390/f13020217
Submission received: 6 January 2022 / Revised: 24 January 2022 / Accepted: 26 January 2022 / Published: 31 January 2022
(This article belongs to the Special Issue Thermal Modification of Wood: Process and Properties)

Abstract

:
This study aims to investigate the influence of thermal modification (TM) on the physical and mechanical properties of wood. For this purpose, the experimental part focused on selected influential parameters, namely temperature, residence time, and density, while the four-point bending strength is obtained as the output parameter. The obtained experimental data are stochastically modeled and compared with the model created by genetic programming (GP). The classical mathematical analysis obtained treatment parameters in relation to the maximum bending strength (T = 187 °C, t = 125 min ρ = 0.780 g/cm3) and compared with the results obtained by genetic algorithm (GA) (T = 208 °C, t = 122 min, and ρ = 0.728 g/cm3). It is possible to obtain models that describe experimental results well with stochastic modeling and evolutionary algorithms.

1. Introduction

Wood is one of the oldest construction materials, and new wood-based materials continue to develop and are being successfully introduced into the engineering and construction marketplace [1]. Wood modification can be defined as improving wood properties by creating new materials. In the last twenty years, there has been a notable increase in the use of modified wood for a variety of applications [2]. Among various modification options, thermal modification (TM) is reported the most important in terms of traded volumes in Europe. Up to 500,000 m3 of TM wood is placed on the market annually [3].
Therefore, modified wood is frequently advertised as a new species of wood [4,5,6,7,8]. Thermally modified timber can improve specific properties such as resistance to biological pests, fire, radiation, water, aggressive chemicals, dimensional changes or mechanical loads [9], and higher brittleness determining its use in the interior (e.g., furniture, cladding, flooring, saunas) and exterior applications (e.g., garden furniture, children’s playgrounds, noise barriers) [4,10,11].
Despite its properties, wood demonstrates some limitations in the exterior environment use. The main disadvantages of wood, such as poor dimensional stability and biological degradation or deterioration, are mainly due to the nature of the main polymers of the wood cell wall, primarily due to the abundance of free hydroxyl groups (OH) [12,13,14].
The study [15,16] states that thermal modification improved dimensional stability and fungal attack resistance, while the TM process reduced the mechanical properties of wood. Thermal modification can be presented as partial pyrolysis in a reaction chamber with almost anoxic conditions. This process results in a changed chemical composition of the modified wood. Among the three key wood components, the hemicelluloses are being the most affected [17]. The first signs of hemicellulose modification appear at rather low temperatures. Decomposition starts at 140 °C and α-cellulose at 150 °C [18]. Lignin is more resilient to elevated temperatures [19]. One of the key results of thermal modification is the reduction of readily available hydroxyl groups, which affect the sorption properties of wood [20]. The equilibrium moisture content of TM wood is therefore considerably lower than that of non-modified wood when determined under the same climatic conditions [21,22]. It should be considered that the water exclusion efficacy is one of the key factors that determine the wood performance in above-ground conditions [23]. The effect of the thermal modification depends on the modification parameters. Shuang-Yan et al. stated that the chemical components are the main parameters affecting the mechanical properties of wood fibers [24]. The thermal modification of wood at different temperatures and durations decreased mechanical properties and increased dimensional stability and biological durability [25,26,27,28,29,30,31,32].
Physical and mechanical properties, modulus of rupture (MOR), modulus of elasticity (MOE), impact bending, compression strength, hardness in tangential and radial direction were determined in research [33,34].
The thermal modification process is influenced by various parameters, is including temperature, wood density, and residence time. Wood processing and testing technologies, which have been applied for many years, can modernize the application of knowledge through modeling, simulation, optimization, process theory, and computer technology [35,36,37]. The experimental results were obtained according to relevant standards for determining some physical and mechanical properties of timber EN 408 [38] and EN 789 [39], and using generally accepted scientific methods based on stochastic modeling and optimization [40,41,42,43,44]. Regarding genetic programming versus stochastic modeling, GP with high selection pressure finds a solution earlier, i.e., it converges faster without compromising the quality of the resulting solution, since faster convergence causes a higher probability of getting trapped in the local optimum [45]. On the other hand, if the selection pressure of the model is too low, computer time is spent on useless iterations because the convergence is too slow. Thus, this state has to be avoided with proper programming and specification of limit conditions [45,46,47].
This study attempts to predict optimal parameters of wood thermal treatment (modification temperature and duration) of wood products based on bending strength using stochastic modeling and genetic programming.

2. Materials and Methods

2.1. Description of the Resources Needed for the Experiment

Experimental research related to this work is based on the thermal modification of samples in the chamber for thermal modification, performed at Biotechnical Faculty in Ljubljana. Unmodified wood samples of European beech (Fagus sylvatica), Linden (Tilia sp.), and Silver fir (Abies alba) for bending process, with the average density of 0.675 g/cm3, 0.472 g/cm3, and 0.366 g/cm3, respectively. For modified samples, the average density for Beech was 0.655 g/cm3, for Linden 0.455 g/cm3, and 0.353 g/cm3 for Fir. The test samples had a minimum length of 19 times the depth of the section. Ten samples without visible defects, full cross-section, free from knots and resin pockets, with the dimensions of 380 mm × 50 mm × 20 mm, were prepared for bending strength test. Before the TM process, all samples were dried to a constant weight, over 24 h in a drying chamber from the “Kambič” manufacturer, at a temperature of 103 °C. Drying was carried out to determine the oven-dry mass for further mass-loss calculations. After drying, the samples were cooled in a desiccator and weighed in an unmodified state. They were then transferred to a vacuum reactor, where the thermal modification process was performed according to a commercial process (Silvapro®, Silvaprodukt, Ljubljana, Slovenia) [48]. After the specimens were placed in a vacuum-pressure reactor, about 95% vacuum was achieved in the chamber, while the absolute pressure was about 5 kPa. The parameters of the heat treatment process for testing the mechanical and physical properties of European beech (Fagus sylvatica), Linden (Tilia sp.), and silver fir (Abies alba) are following the parameters used in commercial processes [48]. The heating of the samples goes through five phases in which the samples were heated to maximum temperatures (170 °C, 180 °C, 195 °C, 210 °C, 220 °C) and treated with different maximum durations of the process (78, 120, 180, 240, 276 min). In the next 12 h, the samples were cooled and weighed.
Measurement of four-point breaking force of unmodified and modified specimens was performed on a testing machine SHIMADZU type SIL-50KNAG at the Faculty of Technical Engineering in Bihać. The specimens were made and tested as recommended in the standard BAS EN 408 + A1 (Figure 1) [49].
Bending force (Fb) testing was performed at four points of ten specimens for each experimental point. The specimens were always positioned so that they were subjected to maximum load. The test specimens, having a minimum length of approximately 19 times the depth of the section, were simply supported and symmetrically loaded in bending at two points throughout approximately 18 times the depth. Figure 1 shows the maximum force measured within the loading points.

2.2. Experimental Procedure and Setup

The experiment was conducted by using a central composition design with five levels of the three main independent parameters: process temperature (X1 = T), process duration (X2 = t), and density (X3 = ρ ), as shown in Table 1 [35]. The overall number of experiments conducted of bending strength:
N = 2 k p + 2 k + n 0 = n k + n α + n 0 = 20  
where: n k number   of   variables ,   n α   number   of   symmetrically   set   points ,     n 0   number   of   repetition   points   in   the   plan   center . In the specific case k = 3   (changeable variables:   T ,   t ,   ρ ) , n 0 = 6 , and n α = 6 ( n α = 2 k ).

2.3. Polynomial Form of Bending Strength Function

Based on the experiment, can present the functional relationship between the response of bending strength and the investigated independent parameters, by the following polynomial form of mathematical model:
Y F = b 0 + b 1 X 1 + b 2 X 2 + b 3 X 3 + b 11 X 1 2 + b 22 X 2 2 + b 33 X 3 2 + b 12 X 1 X 2 + b 13 X 1 X 3 + b 23 X 2 X 3 + b 123 X 1 X 2 X 3  
where states that:
Y F —dependent variable corresponding to the physical or experimental value of bending strength
X 1 ,   X 2 ,   X 3 —independent variables corresponding to the physical values of the processing regime T ,   t ,   ρ
b 0 ,   b 1 ,   b 2   ,   b 3 ,   b 11 ,   b 22 ,   b 33 —regression coefficient of the mathematical model.
The formalized description of Y F = f   T ,   t ,   ρ   changes to mathematical Y F = f   X 1 ,   X 2 ,   X 3 , which describes the influence of their interactions X 12 ,   X 11 ,   X 22 on the output value Y in addition to the influence on all three sizes X 1 ,   X 2 ,   X 3 individually. Table 2 lists the experimental design matrix or coded form of input parameters.
Regression coefficients are calculated as:
b 0 = a 1 j = 1 N y j + a 2 i = 1 k j = 1 N X i j 2 y j
b i = a 3 j = 1 N x i j y j     for   i = 1 ,   2 ,   3 , k
b i m = a 4 j = 1 N X j x m j y j ,   for   1 i   m k
b i i = a 5 j = 1 N X i j 2 y j + a 6 j = 1 k j = 1 N X i j 2 y j + a 7 j = 1 N y j ,   i = 1 ,   2 ,   3 ,
The mathematical analysis provided the model with ten significant coefficients in coded form   b 0 , b 1 ,   b 2 ,   b 3 ,   b 12 , b 13 ,   b 23 ,   b 11 ,   b 22 ,   b 33 . The dispersion is homogenous; thus, the experiment can continue. The mathematical analysis provided the model with ten significant coefficients in coded form:
Y = 8.402 0.1430 X 1 0.1247 X 2 + 0.2015 X 3 0.1515   X 1 2 0.0936 X 2 2 0.1151 X 3 2 0.1172 X 1 X 2 0.0835 X 1 X 3 0.0750 X 2 X 3
After the transformation of Equation (7), the bending strength model as a function of the process temperature (T), process time (t), and density (ρ), has the following physical form:
Y = 27.9482 + 2.2505   t 0.0936   t 2 + 0.2982 T 0.00782   tT 0.00067 T 2 + 16.0142   ρ 0.5   t ρ 0.0371   T ρ 5.1156   ρ 2
A multiple regression coefficient R = 0.97 was obtained, indicating a good interdependence of the mathematical model parameters. The mathematical model describes the four-point bending of thermally modified wood accurately and reliably within the space covered by the experiment.

2.4. Modelling by Genetic Programming

A genetic model of the four-point bending force is developed using GPdotNET software [45]. The output variable is the bending force of thermally modified wood (Table 3). The population size is expected to correlate positively with the quality and rate of convergence of the solution. Mutation probability is a key parameter that allows GP new genetic material to be introduced into the population and obtain better solutions. Therefore, a significant complex effect of mutation on the solution’s quality is expected. The number of generations should be positively correlated with the quality of the solution and negatively correlated with the program execution rate. Figure 2 shows a model of the bending force of thermally modified wood in the form of an expression tree. Table 3 lists the experimental results and the results calculated with the obtained GP model of bending force.
GP model of bending force obtained for a specific GP configuration depends on the following parameters: population size, mutation probability, number of generations, and window size. These parameters can take on a wide range of values, which can significantly affect GP solutions’ performance and quality. The influence of window size and other parameters on the operation of the genetic algorithm can be found in the literature [45]. The following GP parameters were used during modeling:
-
Set of functions F = {+, −, *, /}
-
Input variable vectors X1, X2, X3, Set T= {T, t, ρ}—set of terminals, temperature (T), time (s), and density (ρ).
-
R-makes a set of randomly generated constants that can be found in the expression (r1 = 0.517009973526; r2 = 0.930670022964478)
-
Size of population G = 500,
-
Initial depth of binary wood 5,
-
Depth of wood at mutation and crossing 8,
-
Probability of crossing 90%,
-
Probability of mutation 5%,
-
Probability of reproduction 20%,
-
Selection method, Elite selection,
-
Method of initialization of mixed population ‘’ramped half and half’’
-
Number of iterations (evolutions).
-
The criterion function of chromosome goodness-of-fit testing (computer programs) is defined by multiple regressions, as follows:
R = 1 ( y i y i ˜ ) 2 ( y i y i ¯ ) 2
where states that:
y i —values of experimental data,
y i ¯ —mean value of experimental values
y i ˜ —model values.
In addition to a graphical representation of the model (Figure 2), the GPdotNET program provides the possibility to export the model to other programs in the form of an analytical expression. is shown in the following analytical expression represents the analytical expression of the GP model of the bending strength.
F s = X 1 2 · X 3 6 + r 1 2 ( X 3 6 + X 1 2 · 2 + X 3 2 + X 1 · X 3 5 · 2 + X 3 2 r 2 + X 1 2 + r 1 2 · X 1 4 + r 2 · X 1 4 + X 1 5 · X 3 6
Based on the obtained results, the experiment and the GP model are presented in Figure 3 and Table 3. The presented diagram shows the quality distribution of the solution of the set of parameters of thermally modified wood for bending. The diagrams suggest that the algorithm consistently finds solutions close to the optimum for these sets of parameters. In other words, by changing the input parameters, it is known up to what final breaking force the thermally modified wood can withstand the load. The correlation coefficient between the model and the validation set is 0.99. The data acquired in this way show that the model obtained by the genetic programming method is reliable and can be used in further research.

2.5. Optimization of Bending Strength of Thermally Modified Wood by Classical Mathematical Analysis Method

The obtained mathematical model (8) should be optimized so that the model parameters such as T , t , and ρ can assume optimal values, when the target function, expressed over the maximum four-point bending process, gets the maximum value F c = F m a x .
The objective function F c = F s = F s X 1 ,   X 2 ,   X 3 for the area 1.682 < X i < 1.682 gets maximum   F c = F m a x for coded values X 1 = X 10 ,   X 12 = X 20 , X 3 = X 30 or physical values T = T 0 ,   t = t 0 ,   ρ = ρ 0 . Determining the maximum bending strength extreme values is about differentiating the mathematical model (8) according to the expression (11).
The following condition is used for the determination of the optimum function F c = F m a x :
  F x j = 0   ;   j = 1 , 2 , 3 , , n
By solving the system of equations, we can obtain the coordinates of a stationary point:
X 10 = 0.4975 ,   X 20 = 0.8943   i   X 30 = 1.3471
These are nonlinear equations and solutions ( X 10 , X 20 , X 30 ). Iterative methods obtain stationary points. If there is a minimum, maximum, bending, or saddle point in the stationary point, it is necessary to examine the sign of the determinants ∆1, ∆2, and ∆3, as follows:
1 = F x 1 x 1 = 2 F x x 1 x 1 = 2 F x x 1 2  
Δ 2 = d e t F x 1 x 1 F x 1 x 2 F x 2 x 1 F x 2 x 2 = d e t 2 F x x 1 x 1 2 F x x 1 x 2 2 F x x 2 x 1 2 F x x 2 x 2  
Δ 3 = d e t F x 1 x 1 F x 1 x 2 F x 1 x 3 F x 2 x 1 F x 2 x 2 F x 2 x 3 F x 3 x 1 F x 3 x 2 F x 3 x 3 = d e t 2 F x x 1 x 1 2 F x x 1 x 2 2 F x x 1 x 3 2 F x x 2 x 1 2 F x x 2 x 2 2 F x x 2 x 3 2 F x x 3 x 1 2 F x x 3 x 2 2 F x x 3 x 3
Maximum bending strength function F i = F i m a x is determined with the expressions (12)–(15), as follows:
1 = 0.3030 < 0 ;   2 = 0.0430 > 0 ;   3 = 0.00836 < 0
Based on the calculated values of the determinants ∆i and criteria (12)–(15), the function of the maximum bending strength for X 10 = 0.4975 ,   X 20 = 0.8943   and   X 30 = 1.3471 has a maximum value of FB = FB max. The maximum bending strength was obtained for the values of thermal modification parameters T = 187 °C, t = 126 min and ρ = 0.780 g/cm3. Figure 3 shows the maximum bending force that depends on the input parameters. In contrast, the whole diagram represents the behavior of the maximum bending force and the relationship between the two input variables. The extreme point Fs max (T, t, ρ ) shows the intensity of the change of the bending force depending on the individual change of each of the analyzed bending process parameters. Mathematical models of the maximum force are obtained and determined with two variables through the optimal point, with an orthogonal cross-section:
Y t , ρ = 4.266 + 0.784 t 0.093 t 2 + 9.050 ρ 0.500 t ρ 5.115 ρ 2
  Y T , ρ = 23.624 + 0.281 T 0.000674 T 2 + 14.961 ρ 0.037 T ρ 5.11
  Y T , t = 18.552 + 1.859 t 0.0936 t 2 + 0.269 T 0.0078 t T 0.00067 T 2
Orthogonal cross-section of the objective function Fs = Y (8) through the optimal point, enable mathematical models of the maximum bending force, determined by one physical variable:
Y T = 15.052 + 0.252 T 0.00067 T 2 ,
Y t = 8.216 + 0.393 t 0.0936
Y p = 5.503 + 7.998 p 5.115 p 2
Figure 4 shows the curves of the thermal modification parameters for the maximum four-point bending concerning the different input variables and optimal points.
The application of stochastic modeling, based on the statistical method of defining a mathematical model, confirms the determined dependence of the maximum four-point bending force and the optimal parameters of the thermal modification process. The obtained mathematical model is adequate for the optimization of the thermal modification parameters if the maximum bending force is taken as the criterion of optimality.

2.6. Optimization of Bending Strength of Thermally Modified Wood by Genetic Algorithm

Genetic algorithms use stochastic search methods, meaning that multiple population goodness in a sequence can remain approximately the same over generations. After a while, a generation with superior goodness will undoubtedly emerge. This way of working with genetic algorithms can cause a problem in defining the objective function. A common practice in applying genetic algorithms is to complete the genetic algorithm after a certain number of generations (iterations). The best chromosomes in the population are then studied. If no satisfactory solution is obtained, the genetic algorithm is restarted. The mathematical model (10) of thermally modified four-point bending wood obtained by genetic programming was derived according to a genetic algorithm. A gene representation over real numbers was used for the genetic algorithm. The GA parameters used for optimization are as follows:
-
population size 500,
-
number of iterations 272,
-
probability of mutation 5%,
-
crossing probability 90%,
-
probability of reproduction 20%,
-
rank selection method.
The maximum value of the mathematical model was calculated for T = 208 °C, t = 121 min, and ρ = 0.728 g/cm3, where the maximum bending force of the thermally modified wood was F = 5950 N. The comparison of this result with the maximum value obtained by the regression analysis, where the maximum force is F = 6438 N, indicates that the forces are approximately equal.

3. Results and Discussion

3.1. Comparative Results of Bending Strength Experiment and Modelling

Genetic programming is a very efficient method of model making, and unlike classical regression analysis, there are no restrictions on the degree of polynomials. Based on the concept of computer programs, this method is an intuitive and simple way to build a model over experimental data that can then be used in future research in various wood tests as a starting point. The results showed that at maximum flexural fracture force, genetic programming can yield results that are closer to experimental data than classical regression analysis, that is, better results than stochastic modeling (classical mathematical analysis). Comparative results obtained by stochastic and genetic modeling of the maximum bending strength of thermally modified wood are presented, i.e., the reliability of using mathematical models without conducting expensive experimental research. Table 3 lists the comparative results of the bending experiment and the results obtained according to stochastic mathematical modeling and genetic modeling.
Table 3 summarizes that the results of the maximum four-point bending force, obtained by the stochastic mathematical model and the genetic model, are almost indistinguishable, which further confirms the adequacy of the obtained models (Figure 5).
Figure 5 represents the results obtained by methods of mathematical modeling and GP are approximate values concerning the experimental results of the maximum four-point bending force.
Not many models are developed that enable the optimization of wood properties based on density-independent of wood species. The majority of the current approaches are based on considering two parameters, modification temperature, and duration of the modification process [50]. These models are based on the mass loss of the timber after thermal modification and are predominately focused on durability. Up to our knowledge, this is the first approach to address the mechanical properties of wood with models based on the density of wood, modification temperature, and modification duration.

3.2. The Comparison of the Optimal Results

The obtained models (8) and (10) present a very good base for finding the optimal parameters to increase process efficiency, which is verified by the confirmation test. The presented optimization approaches provide the same optimal values for the TM parameters and give accurate results (as indicated by the confirmation test) with a slight deviation between each other (Table 4).
The modeling provides a helpful tool for optimizing thermally modified wood properties. A considerable part of the chambers for thermal modification is somewhat of small capacity. This enables the producers to adopt the process based on the consumers’ quality (density) and the expectation (needs). Thus, the respective solutions enable the higher added value of the thermally modified wood and better competitiveness on the market.

4. Conclusions

By mathematical modeling of the four-point bending process, have obtained a model in a coded form that has ten significant coefficients, which means that the relationship between (T), (t), and the wood ( ρ ) plays a significant role in the adequacy of the mathematical model of bending strength.
The study shows that stochastic modeling and evolutionary algorithms can provide reliable results that are consistent with experimental data. By optimizing the obtained models of maximum bending force, using the method of classical mathematical analysis and evolutionary algorithm, optimal values for T , t , and ρ were obtained. The optimal parameters of the thermal modification process by classical mathematical analysis and genetic algorithms have approximately the same value and represent an ideal treatment that can be applied to the wood.
There is the possibility of further research in finding a model for predicting the bending force of fracture which will be based on additional parameters such as the influence of humidity, or the behavior of the timbers’ different conditions of exploitation, using GA in the input parameters optimization process.

Author Contributions

Conceptualization, M.H. and A.H.; methodology, M.B., M.P. and D.K.; software, R.H.; validation, R.H. and D.H.; formal analysis, M.H., R.H. and M.P.; investigation, R.H., M.H., M.P., D.H.; resources, R.H., M.H. and M.P.; data curation, R.H., M.B., D.K. M.H. and A.H.; writing—original draft preparation, R.H., M.P., M.H., D.K. and A.H.; writing—review and editing, R.H., M.B. and D.C.; visualization, R.H. and D.H.; supervision, M.H., M.P.; project administration, R.H. and A.H.; funding acquisition, R.H. and D.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by “Natural Science Foundation of Anhui Province (No. 2008085QC130)”, “Bilateral project between the Republic of Slovenia and Bosnia and Herzegovina’’, and the research program “P4-0015 Wood and lignocellulosic composites” financed by the Slovenian Research Agency (ARRS).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data included in this study are available upon request by contact with the corresponding author.

Acknowledgments

The authors gratefully acknowledge the support of the “Natural Science Foundation of Anhui Province (No. 2008085QC130)”, “Bilateral project between the Republic of Slovenia and Bosnia and Herzegovina’’, and the research program “P4-0015 Wood and lignocellulosic composites” financed by the Slovenian Research Agency (ARRS).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kržišnik, D.; Grbec, S.; Lesar, B.; Plavčak, D.; Šega, B.; Šernek, M.; Straže, A.; Humar, M. Durability and Mechanical Performance of Differently Treated Glulam Beams during Two Years of Outdoor Exposure. Drv. Ind. 2020, 71, 243–252. [Google Scholar] [CrossRef]
  2. Goli, G.; Marcon, B.; Fioravanti, M. Poplar wood heat treatment: Effect of air ventilation rate and initial moisture content on reaction kinetics, physical and mechanical properties. Wood. Sci. Technol. 2014, 48, 1303–1316. [Google Scholar] [CrossRef]
  3. Jones, D.; Sandberg, D.; Goli, G.; Todaro, L. Wood Modification in Europe: Processes, Products; University of Florence: Florence, Italy, 2018; Volume 41. [Google Scholar]
  4. Hill, C.A. Wood Modification: Chemical, Thermal and Other Processes, 3rd ed.; John Wiley & Sons, USA: Hoboken, NJ, USA, 2007; pp. 99–100. [Google Scholar]
  5. Esteves, B.; Carmo, J.; Nunes, L. Commercialisation and production of modified wood. In Proceedings of the European Conference on Wood Modification, Portugal, Lisbon, 10–12 March 2014. [Google Scholar]
  6. Kržišnik, D.; Lesar, B.; Thaler, N.; Humar, M. Micro and material climate monitoring in wooden buildings in sub-Alpine environments. Constr. Build. Mater. 2018, 166, 188–195. [Google Scholar] [CrossRef]
  7. Popescu, C.M.; Jones, D.; Kržišnik, D.; Humar, M. Determination of the effectiveness of a combined thermal/chemical wood modification by the use of FT-IR spectroscopy and chemometric methods. J. Mol. Struct. 2020, 1200, 127–133. [Google Scholar] [CrossRef]
  8. Sandberg, D.; Kutnar, A.; Karlsson, O.; Jones, D. Wood Modification Technologies Principles, Sustainability, and the Need for Innovation, 1st ed.; Taylor & Francis Group: Oxfordshire, UK, 2021. [Google Scholar]
  9. Reinprecht, L. Wood Deterioration, Protection and Maintenance, 1st ed.; John Wiley & Sons: Hoboken, NJ, USA, 2016. [Google Scholar]
  10. Boonstra, M.J.; Blomberg, J. Semi-isostatic densification of heat-treated radiata pine. Wood Sci. Technol. 2007, 41, 607. [Google Scholar] [CrossRef] [Green Version]
  11. Hasanagić, R. Thermal modified wood: A critical review. In Proceedings of the RIM 2019—Development and Modernization of Production, Sarajevo, Bosnia and Herzegovina, 10–14 September 2019. [Google Scholar]
  12. Schmidt, O. Wood and Tree Fungi, 1st ed.; Springer: Berlin/Heidelberg, Germany, 2006; p. 100. [Google Scholar]
  13. Van den Bulcke, J.; De Windt, I.; Defoirdt, N.; De Smet, J.; Van Acker, J. Moisture dynamics and fungal susceptibility of plywood. Int. Biodeterior. Biodegrad. 2011, 65, 708–716. [Google Scholar] [CrossRef] [Green Version]
  14. Eaton, R.A.; Hale, M.D.C. Wood: Decay, Pests, and Protection; Chapman & Hall: London, UK; New York, NY, USA,, 1993. [Google Scholar]
  15. Tjeerdsma, B.F.; Stevens, M.; Militz, H. Durability Aspects of Hydrothermal Treated Wood; International Research Group on Wood Preservation: Stockholm, Sweden, 2000. [Google Scholar]
  16. Boonstra, M.J.; Tjeerdsma, B. Chemical analysis of heat treated softwoods. Holz Roh- Werkst. 2006, 64, 204–211. [Google Scholar] [CrossRef]
  17. Hill, C.A.S. Wood modification: An update. BioResources 2011, 2, 918–919. [Google Scholar]
  18. Zauer, M.; Kowalewski, A.; Sproßmann, R.; Stonjek, H.; Wagenführ, A. Thermal modification of European beech at relatively mild temperatures for the use in electric bass guitars. Eur. J. Wood Wood Prod. 2016, 74, 43–48. [Google Scholar] [CrossRef]
  19. Korošec, R.C.; Lavrič, B.; Rep, G.; Pohleven, F.; Bukovec, P. Thermogravimetry as a possible tool for determining modification degree of thermally treated Norway spruce wood. J. Therm. Anal. Calorim. 2019, 98, 189–195. [Google Scholar] [CrossRef]
  20. Tjeerdsma, B.F.; Boonstra, M.; Pizzi, A.; Tekely, P.; Militz, H. Characterisation of thermally modified wood: Molecular reasons for wood performance improvement. Holz Roh- Werkst. 1998, 56, 149–153. [Google Scholar] [CrossRef]
  21. Brito, A.F.; Calonego, F.W.; Bond, B.H.; Severo, E.T.D. Color changes, EMC, and biological resistance of thermally modified yellow poplar. Wood Fiber Sci. 2018, 50, 439–446. [Google Scholar] [CrossRef]
  22. Kozakiewicz, P.; Drozdzek, M.; Laskowska, A.; Grześkiewicz, M.; Bytner, O.; Radomski, A.; Zawadzki, J. Effects of thermal modification on selected physical properties of sapwood and heartwood of black poplar (Populus nigra L.). BioResources 2019, 14, 8391–8404. [Google Scholar]
  23. Meyer-Veltrup, L.; Brischke, C.; Alfredsen, G.; Humar, M.; Flæte, P.O.; Isaksson, T.; Brelid, P.L.; Westin, M.; Jermer, J. The combined effect of wetting ability and durability on outdoor performance of wood: Development and verification of a new prediction approach. Wood Sci. Technol. 2017, 51, 615–637. [Google Scholar] [CrossRef]
  24. Shuang-Yan, Z.; Fei, B.-H.; Yu, Y.; Cheng, H.-T.; Wang, C.-G. Effect of the amount of lignin on tensile properties of single wood fibers. For. Sci. Pract. 2013, 15, 56–60. [Google Scholar]
  25. Militz, H. Thermal Treatment of Wood: European Processes and Their Background; International Research Group on Wood Preservation: Stockholm, Sweden, 2000. [Google Scholar]
  26. Rahimi, S.; Singh, K.; DeVallance, D. Effect of different hydrothermal treatments (steam and hot compressed water) on physical properties and drying behavior of yellow-poplar (Liriodendron tulipifera). For. Prod. J. 2019, 69, 42–52. [Google Scholar] [CrossRef]
  27. Le, X.P.; Satoshi, S.; Yukie, S. Effects of heat treatment on brittleness of Styrax tonkinensis wood. J. Wood Sci. 2007, 53, 181–186. [Google Scholar]
  28. Esteves, B.M.; Pereira, M.H. Wood modification by heat treatment: A review. BioResources 2008, 4, 370–404. [Google Scholar] [CrossRef]
  29. Kaboorani, A.; Englund, K.R. Water sorption and mechanical performance of preheated wood/thermoplastic composites. J. Compos. Mater. 2010, 45, 1423–1433. [Google Scholar] [CrossRef]
  30. Gokhan, G.; Suleyman, K.; Deniz, A.; Ilter, B. The density, compression strength and surface hardness of heat treated hornbeam (Carpinus betulus L.) wood. Maderas Cienc. Tecnol. 2009, 1, 61–70. [Google Scholar]
  31. Otávio, J.; Brito, J.; Nivaldo, G.; Geraldo, B.; Júnior, A.; Pessoa, M.; Paulo, M.H. The density and shrinkage behavior of Eucalyptus grandis wood suubmitted to difrent temperatures of termoretification. Cerne 2015, 12, 182–188. [Google Scholar]
  32. Hasanagić, R.; Sauradipta, G.; Bajramović, E.; Hasanagić, A. Mechanical properties changes in fir wood (abies sp.), linden wood (tilia sp.), and beech wood (fagus sp.) subjected to various thermal modification process conditions. IOP Conf. Ser. Mater. Sci. Eng. 2021, 1208, 012025. [Google Scholar] [CrossRef]
  33. Bilgin, I.; Abdullah, B. Effects of industrial heat treatment on some physical and mechanical properties of iroko wood. Drvna industrija 2017, 68, 229–2369. [Google Scholar]
  34. Vasiliki, K.; Barboutis, I.; Vasileios, V. Influence of thermal treatment on mechanical strength of Scots pine (Pinus sylvestris L.) wood. Wood Res. 2014, 59, 373–378. [Google Scholar]
  35. Jurković, M. Mathematical Modelling and Optimization of Machining Processes, 2nd ed.; Faculty of Engineering, University of Rijeka: Rijeka, Croatia, 2006. [Google Scholar]
  36. Jurković, M. Matematičko Modeliranje Nžinjerskih Procesa i Sistema, 1st ed.; Univerzitet u Bihaću-Tehnički fakultet: Bihać, Bosnia & Herzegovina, 1999. [Google Scholar]
  37. Dean, A.; Morris, M.; Stufken, J.; Bingham, D. Handbook of Design and Analysis of Experiments, 1st ed.; Chapman and Hall/CRC: New York, NY, USA, 2015; p. 591. [Google Scholar]
  38. Montgomery, D.C. Design and Analysis of Experiments, 3rd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2017. [Google Scholar]
  39. EN 408. In Timber and Glued Laminated Timber-Determination of Some Physical and Mechanical Properties.; Swedish Institute for Standards: Stockholm, Sweden, 1995.
  40. EN 789; Timber structures—Test methods—Determination of Mechanical Properties of Wood-Based Panels; European Committee for Standardization: Brussels, Belgium, 2004.
  41. Jenkel, C.; Leichsenring, F.; Graf, W.; Kaliske, M. Engineering Structures, Stochastic modelling of uncertainty in timber engineering. Elsevie 2015, 99, 296–310. [Google Scholar]
  42. Hodžić, A.; Hasanagić, R. Matematičko modeliranje sile loma na istezanje masivnog drveta produženog sa zupčastim vezom. In Proceedings of the Zbornik Međunarodnog Kongresa o Procesnoj Industriji—Procesing, Beograd, Srbija; 2017. [Google Scholar]
  43. Hasanagić, R.; Hodžić, A.; Jurković, M. Modelling and optimization of tensile break force of solid wood elements lengthened by finger joint. J. Adhes. Sci. Technol. 2020, 34, 1013–1027. [Google Scholar] [CrossRef]
  44. Hasanagić, R.; Hodžić, A.; Bjelić, A. Mathematical modelling of the surface roughness in orthogonal cutting of the massive wood. Wood Mater. Sci. Eng. 2019, 16, 74–80. [Google Scholar] [CrossRef]
  45. Hrnjica, B.; Danandeh, M.A. Optimized Genetic Programming Applications: Emerging Research and Opportunities, 1st ed.; IGI Global: Hershey, PA, USA, 2018; p. 190. [Google Scholar]
  46. Hasanagić, R. Modeling and prediction of fracture force to tighten solid wood elements by genetic programming. Tehnika 2018, 73, 653–657. [Google Scholar] [CrossRef]
  47. Gondia, A.; Ezzeldin, M.; El-Dakhakhni, W. Mechanics-guided genetic programming expression for shear-strength prediction of squat reinforced concrete walls with boundary elements. J. Struct. Eng. 2020, 11, 04020223. [Google Scholar] [CrossRef]
  48. Rep, G.; Pohleven, F.; Kosmerl, S. Development of the Industrial kiln for Thermal Wood Modification by a Procedure with an Initial Vacuum and Commercialisation of Modified Silvapro Wood. In Proceedings of the 6th European Conference on Wood Modification, Ljubljana, Slovenia, 17–18 September 2012; Jones, M.H., Militz, M.P.D., Petrič, H., Pohleven, M.F., Eds.; Slovenia University of Ljubljana: Ljubljana, Slovenia, 2012. [Google Scholar]
  49. BAS EN 408+A1. In Timber Structures-Structural Timber and Glued Laminated Timber-Determination of Some Physical and Mechanical Properties; Institute for Standardization BIH: Sarajevo, Bosnia and Herzegovina, 2013.
  50. Altgen, M.; Welzbacher, C.; Humar, M.; Willems, W.; Militz, H. Bestimmung der Behandlungsqualitat von Thermoholz mithilfe von Schnellverfahren. Holztechnologie 2012, 6, 44–49. [Google Scholar]
Figure 1. Bending test: (a) 4-point bending procedure, (b) rectangular cross-sectional sample dimensions [49].
Figure 1. Bending test: (a) 4-point bending procedure, (b) rectangular cross-sectional sample dimensions [49].
Forests 13 00217 g001
Figure 2. The model of bending strength in the form of an expression tree.
Figure 2. The model of bending strength in the form of an expression tree.
Forests 13 00217 g002
Figure 3. Graphic representation of the function: (a) (YT, ρ maximum bending force depending on the duration of the process (t) and wood density ( ρ ), (b) (YT, ρ ) maximum bending force depending on the process temperature (T) and wood density ( ρ ), (c) (YT, t) maximum bending force depending on the process temperature (T) and the duration of the process (t).
Figure 3. Graphic representation of the function: (a) (YT, ρ maximum bending force depending on the duration of the process (t) and wood density ( ρ ), (b) (YT, ρ ) maximum bending force depending on the process temperature (T) and wood density ( ρ ), (c) (YT, t) maximum bending force depending on the process temperature (T) and the duration of the process (t).
Forests 13 00217 g003
Figure 4. Graphic representation of the function of: (a) YT thermal modification parameters on bending strength depending on the temperature of the thermal modification process (T), (b) Yt thermal modification parameters on bending strength depending on the duration of the thermal modification process (t), (c) Y ρ thermal modification parameters on bending strength depending on the wood density ( ρ ).
Figure 4. Graphic representation of the function of: (a) YT thermal modification parameters on bending strength depending on the temperature of the thermal modification process (T), (b) Yt thermal modification parameters on bending strength depending on the duration of the thermal modification process (t), (c) Y ρ thermal modification parameters on bending strength depending on the wood density ( ρ ).
Forests 13 00217 g004
Figure 5. Comparative results of experiments (blue curve), stochastic model (black curve), and GP model (red curve) in bending.
Figure 5. Comparative results of experiments (blue curve), stochastic model (black curve), and GP model (red curve) in bending.
Forests 13 00217 g005
Table 1. Levels of bending strength parameters.
Table 1. Levels of bending strength parameters.
Parameters/LevelsLowestLowCentreHighHighest
Coding–classical experimental design−1.682−1011.682
Temperature (°C) X1 = T170180195210220
Process duration (min) X2 = t78120180240276
Density (g/cm3)
X3 = ρ
0.3300.4300.5800.7300.830
Coding-orthogonal array (Xi)−1.682−1011.682
Table 2. Modeling matrix.
Table 2. Modeling matrix.
N Species *T °Ct min ρ
g/cm3
X0X1X2X3X1X2X1X3X2X3X1X2X3X12X22X32
1 Fir1801200.431−1−1−1111−1111
2 Fir2101200.4311−1−1−1−111111
3 Fir1802400.431−11−1−11−11111
4 Fir2102400.43111−11−1−1−1111
5 Beech1801200.731−1−111−1−11111
6 Beech2101200.7311−11−11−1−1111
7 Beech1802400.731−111−1−11−1111
8 Beech2102400.7311111111111
9–14 Linden1951800.5810000000000
15 Linden1701800.581−α000000(−α)200
16 Linden2201800.581α000000α200
17 Linden195780.5810−α000000(−α)20
18 Linden1952760.5810α000000α20
19 Fir1951800.33100−α000000(−α)2
20 Beech1951800.83100α000000α2
* Ten samples were measured for each experimental point.
Table 3. Results of bending strength breaking force on the bending of the TMD experiment and model.
Table 3. Results of bending strength breaking force on the bending of the TMD experiment and model.
N *Experimental
Results Y [N]
Standard
Deviations
Results Per Models
Stochastic
Model(YR) [N]
Genetic Model(GP) [N]
126652.0525222569
230814.2128303146
328733.4028872928
421644.0223591889
561506.1051836438
652994.0041655053
751107.0643944880
826536.4722102977
944035.3244594283
1042014.1142374461
1148565.1543964621
1244384.0645594161
1346063.2743794282
1443945.8449594299
1540153.1136954054
1622005.0422832602
1741233.7542204628
1829764.1027742905
1927615.0622932897
2039336.2445194502
* Ten samples were tested for each experimental point.
Table 4. The comparison of the optimal results obtained.
Table 4. The comparison of the optimal results obtained.
MethodOptimal TM Parameters
T (°C)t (min)ρ (g/cm3)
Classic mathematical analysis1931260.780
Genetic algorithm (GA)1971210.728
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chu, D.; Hasanagić, R.; Hodžić, A.; Kržišnik, D.; Hodžić, D.; Bahmani, M.; Petrič, M.; Humar, M. Application of Temperature and Process Duration as a Method for Predicting the Mechanical Properties of Thermally Modified Timber. Forests 2022, 13, 217. https://doi.org/10.3390/f13020217

AMA Style

Chu D, Hasanagić R, Hodžić A, Kržišnik D, Hodžić D, Bahmani M, Petrič M, Humar M. Application of Temperature and Process Duration as a Method for Predicting the Mechanical Properties of Thermally Modified Timber. Forests. 2022; 13(2):217. https://doi.org/10.3390/f13020217

Chicago/Turabian Style

Chu, Demiao, Redžo Hasanagić, Atif Hodžić, Davor Kržišnik, Damir Hodžić, Mohsen Bahmani, Marko Petrič, and Miha Humar. 2022. "Application of Temperature and Process Duration as a Method for Predicting the Mechanical Properties of Thermally Modified Timber" Forests 13, no. 2: 217. https://doi.org/10.3390/f13020217

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop