Next Article in Journal
Analysis and Multi-Objective Optimization for Reducing Energy Consumption and Improving Surface Quality during Dry Machining of 304 Stainless Steel
Next Article in Special Issue
Magnesia (MgO) Production and Characterization, and Its Influence on the Performance of Cementitious Materials: A Review
Previous Article in Journal
Data Merging of AE Sensors with Different Frequency Resolution for the Detection and Identification of Damage in Oxide-Based Ceramic Matrix Composites
Previous Article in Special Issue
Concrete for Precast Blocks: Binary and Ternary Combination of Sewage Sludge Ash with Diverse Mineral Residue
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of the Drying Method of Pine and Beech Wood on Fracture Toughness and Shear Yield Stress

1
Department of Manufacturing and Production Engineering, Faculty of Mechanical Engineering, Gdańsk University of Technology, Gabriela Narutowicza Street 11/12, 80233 Gdańsk, Poland
2
InnoRenew CoE, Livade 6, 6310 Izola, Slovenia
3
Andrej Marušič Institute, University of Primorska, Muzejski trg 2, 6000 Koper, Slovenia
4
HS Hydromech, Wybickiego 21, 83050 Lublewo Gdańskie, Poland
5
The Institute of Fluid-Flow Machinery, Polish Academy of Sciences, Fiszera 14, 80231 Gdansk, Poland
*
Author to whom correspondence should be addressed.
Submission received: 23 September 2020 / Revised: 16 October 2020 / Accepted: 19 October 2020 / Published: 21 October 2020

Abstract

:
The modern wood converting processes consists of several stages and material drying belongs to the most influencing future performances of products. The procedure of drying wood is usually realized between subsequent sawing operations, affecting significantly cutting conditions and general properties of material. An alternative methodology for determination of mechanical properties (fracture toughness and shear yield stress) based on cutting process analysis is presented here. Two wood species (pine and beech) representing soft and hard woods were investigated with respect to four diverse drying methods used in industry. Fracture toughness and shear yield stress were determined directly from the cutting power signal that was recorded while frame sawing. An original procedure for compensation of the wood density variation is proposed to generalize mechanical properties of wood and allow direct comparison between species and drying methods. Noticeable differences of fracture toughness and shear yield stress values were found among all drying techniques and for both species, but only for beech wood the differences were statistically significant. These observations provide a new highlight on the understanding of the effect of thermo-hydro modification of wood on mechanical performance of structures. It can be also highly useful to optimize woodworking machines by properly adjusting cutting power requirements.

1. Introduction

Trees contain a high amount of water during their lifetime that is necessary for their survival. Consequently, logs delivered to the sawmill for the downstream conversion are wet, frequently reaching absolute moisture content above 100%, depending on the wood species, season, storing conditions or time after harvesting. Though, a majority of wood products are used in a dry state (10–15%) because of the wood tendency of reaching moisture equilibrium with the surrounding air. Drying of wood is, therefore, an important step in the manufacturing process of any wood product [1]. Several methods of drying wood are commonly used in the modern wood industry with air drying, kiln drying and vacuum drying recognized as the most frequent [1,2,3,4]. Alternative methods using super-critical drying, higher temperature schedules, low pressures, assistance of microwaves, among others, were developed to shorten time of this process [5,6]. Nevertheless, a consequence of each drying event is a change of the native material properties that can occur to a different extent depending on the wood species and process settings—particularly exposure time and treatment temperature.
Herrera-Díaz et al. [7] analyzed an effect of the air temperature during selected drying processes on the mechanical properties of pine wood, including changes to the modulus of elasticity (MOE) and modulus of rupture (MOR). Two drying methods (air and kiln) as well as thermal modification of radiata pine wood were investigated with a conclusion that values of MOR decreased with increase of the process temperatures. Conversely, the MOE was unchanged independent of the process type. A similar analysis for spruce wood dried in lower temperatures was performed by Oltean et al. [8]. The result of this study revealed that a significant effect on the mechanical properties (MOE, MOR, bending strength) was noticed for wood drying temperatures above 80 °C. Roszyk et al. [9] reported that high treatment temperatures of wood affect the MOE and Rc (relative/compressive strength) in different directions related to the fibers but do not influence the intrinsic anisotropy of mechanical parameters. Furthermore, Borrega and Kärenlampi [10] noted that the kinetics of the drying process affect an alteration of the mechanical properties of dried wood to the same extent as drying temperature. An effect of the drying method combined with the process temperature on physical properties of dried wood was studied by several researchers [3,11,12,13,14,15]. It was also reported that drying of wood may result in an increase of the material porosity, especially after delignification process, as investigated by Vitas et al. [16]. An effect of the drying method on the parallel compressive strength of bamboo was reported by Arantes et al. [17]. Three methods of drying, including air drying, fire-drying and kiln drying, were compared. A slight reduction of the compressive strength parallel to the fiber direction was detected in the case of the bamboo kiln drying. Generally, the mechanical properties of the wood fall sharply as a result of temperature evenly in all directions, as demonstrated by Fonseca and Barreira [18] and Fonseca et al. [19] on the example of several representative wood species.
Application of a vacuum as an addition to the classical wood drying changes the mechanisms of moisture transport and, consequently, kinetics of the process and resulting material modifications [20,21]. Blanchet et al. [22] reported a small negative effect of the vacuum drying process on the mechanical properties of wood, when compared to conventional oven drying. However, noticeable chemical changes in wood polymers as well as its hygroscopic properties because of high-temperature heat treatment were reported by several researchers [23,24,25,26]. These chemical changes are expected, therefore, to affect other material characteristics including strength and mechanical properties. Indeed, a gradual reduction of poplar’s mechanical strength was observed by Sandak et al. [27] when increasing the vacuum thermal treatment time, temperature and, consequently, the process intensity. An elevated temperature of wood drying, and its thermal modification, resulted in the higher fragmentation and granularity of sawdust obtained during sawing, milling or sanding processes [28,29,30,31]. This indicates an extensive change in the fracture properties of such treated wood.
It can be concluded, therefore, that diverse drying methods affect (to varying degrees) the mechanical, physical and chemical properties of dried wood. It is a common understanding that rapid reduction of dynamic mechanical properties is the first and most noticeable consequence of drying or thermal treatment. It is related to increased brittleness of the resulting materials. Hence, it can be hypothesized that drying affects the related mechanical properties, particularly fracture toughness R and shear yield stress τγ. Both properties are relatively difficult to be estimated in routine material characterization tests but are highly relevant for the engineered structure safety assurance as defined in building codes. R and τγ are also fundamental for the proper estimation of cutting forces and corresponding cutting power. Optimization of cutting processes is highly relevant when considering the overall economy of production and manufacturing process in woodworking and other sectors. An accurate forecasting of the power demand for the cutting of wood allows optimal choice of the process configuration, including the number of saws and their spacing in frame saws, multi-circulars saws or tandem band saws. This ensures an effective use of machines as well as prevents overloading that can result in the damage of tools or other losses associated with long downtime in the manufacturing process [32,33,34].
A consistent model of cutting power requirements based on mechanical properties of wood (both fracture toughness R and shear yield stress τγ) has been proposed in the proceeding works [35,36,37,38]. This model is comprehensive for diverse cutting configurations and can be used to predict the cutting forces acting on different process kinematics or tools of various geometries [38,39,40,41].
The goal of this study was to determine variations within fracture toughness and shear yield stresses of hard and soft woods due to differences in the applied drying procedures. The mentioned material properties could be useful for engineers designing wooden structures as well as for the proper estimation of the cutting power requirements.

2. Materials and Methods

2.1. Materials

Scots pine (Pinus sylvestris L.) and beech (Fagus silvatica L.) species were used for preparation of experimental samples. The green wood, originated from Pomeranian District of Poland, was not exposed to any intended drying processes nor thermal treatment before preparation of experimental samples. Boards from the sawmill were exposed to four different drying processes following standard industrial protocols as well as laboratory prototypes. Ten blocks with dimensions W = 60 mm × Hp = 60 mm × Lp = 700 mm (width×height×length, respectively) were cut-out from randomly selected boards. The wood density, considered here as a ratio of air-dry wood mass to its volume, was measured separately on each block and results are summarized in Table 1.
The moisture content MC of wood was determined with resistance-based moisture content meter WRD 100 (TANEL SJ, Gliwice, Poland) with declared accuracy of ±2%. The drying of wood was performed in four alternative scenarios, with measured initial moisture content MC = 40% for pine and MC = 70% for beech.
A schematic representation of four drying process configurations is summarized in Figure 1 and most relevant technical details are described below.

2.1.1. Modified Air-Drying Process

The first method for investigation of the changes to the wood induced by drying was a classic air-drying process [1]. In that case, the batch of green wood samples was stored for two years, outdoors under a roof, at the Gdansk University of Technology campus, assuring proper ventilation and shadowing. After initial phase of natural drying, the wood moisture content reached approximately MC = 16% for both species. In order to reduce the MC to the usual level of indoor use, all the samples were additionally conditioned in the laboratory conditions assuring constant air temperature of 20 °C and relative humidity of 55%. The second phase of air drying took three months, and the resulting MC was reduced to 10%. Such level of MC was assumed as similar to that obtained when drying wood according to state-of-the-art industrial solutions.

2.1.2. Conventional Kiln Drying Process

The second batch of experimental samples was dried in the industrial kiln O.S Panto 120/F (PANTO, Szczytno, Poland) installed at the wood processing company PHU Drew-Met from Kiełpino, Poland. Experimental samples were mixed with a similar batch of wood that was used for the production needs of the factory. The drying process implemented followed the routine drying schedules of the company in varying temperatures of 30 °C to 55 °C for beech and 35 °C to 75 °C for pine. The drying process control assured adaptation of air parameters following MC changes of processed wood. It took 14 and 37 days to dry pine and beech samples, respectively. It is a standard procedure implemented in the PHU Drew-Met company to cool down the kiln and condition the wood before the conclusion of the drying process. The final moisture content of both investigated woods was MC = 10%.

2.1.3. Vacuum Kiln Drying Process

An alternative to classic kiln drying, widely implemented in the wood industry, is vacuum kiln drying that was also tested in this research. The third batch of experimental samples was processed at the same company as kiln drying (PHU Drew-Met, Kiełpino, Poland) using the industrial vacuum kiln SP-5 (LAC S.R.O., Židlochovice, Czech Republic). The wood was heated by direct contact with hot plates and conduction [3]. The drying schedule followed routine settings of the company. The vacuum inside of the kiln was 100 mbar (10 kPa) with temperature of 65 °C for 88 h when drying pine wood. The process became two-steps in the case of beech with the initial pre-drying phase of 55 °C at 10 kPa (84 h), followed by the final drying phase of 65 °C at 6 kPa for 65 h.

2.1.4. Warm Air–Steam Mixture Experimental Drying Process

The fourth batch of samples was dried in the experimental kiln developed at the Gdansk University of Technology (GUT). The technical details regarding that innovative process are presented in the work of Baranski [42]. The method uses heated steam as a drying medium and resulting materials are considered as a hybrid thermally treated wood. The specific settings of the drying schedule were adopted taking into consideration available recommendations [43,44]. The process consisted of two-stages where temperature of the medium in the kiln increased from 65 °C (stage one) to 80 °C at the second stage. The first stage lasted until the wood reached fiber saturation point FSP. In the case of pine, it was approximately 32 h, and for beech 125 h. Total time of drying process for pine was approximately 90 h and in the case of beech 300 h. The relative humidity RH of the heated steam was constant along the process (RH = 80%), while the drying medium flow velocity was 2.5 m∙s−1. The whole batch of wooden boards was cooled down before opening the kiln to reduce the induced stresses.

2.2. Machinability Tests

Experimental cuttings were performed on the PRW15M sash gang saw with a hybrid dynamically balanced driving system and elliptical teeth trajectory movement. The concept of the machine was developed at Gdańsk University of Technology [45] and prototype manufactured by REMA S.A. (Reszel, Poland). The use of electric power (active and passive) during idling and working cycles was continuously monitored with the power converter PP54 (LUMEL S.A., Zielona Góra, Poland). The data were recorded with a time stamp and further processed to determine energetic effects of cutting. A detailed list of sawing machine settings and used tool characteristics is summarized in Table 2.
The mean value of feed per tooth fz for a sash gang saw was calculated as in Equations (1) and (2) [31,41]:
f z = 1000 v f t p n F H F
v f = L p t c
where: vf—feed speed (m·min−1), tp—tooth pitch (mm), Lp—length of the sample (m), HF—saw frame stroke (mm), nF—number of strokes of saw frame per min (spm) and tc—cutting time (min) necessary to process sample of the length Lp.
The average cutting power Pc was calculated as the difference of the mean total power PT and the average idle power Pi [41,46], as expressed in Equation (3):
P c = P T P i ,
The average idle power Pi of the frame saw PRW15-M was determined each time before initiation of the proper cutting cycle. It allowed minimization of an effect of the varying temperature of the machine components (such as hydraulic oil, gear boxes, etc.) on the energetic effects corresponding directly to the cutting process. The average cutting power in a working stroke Pcw was calculated as in Equation (4), following the works [41,46]:
P c w = 2 P c ,

2.3. Methodology for Determination of Material Properties from the Cutting Test

The average total cutting power in the working stroke PcT for a single saw blade in the sash gang can be determined by means of the cutting forces model proposed by Atkins [47,48]. This methodology was adopted for the case of cutting wood on the frame sawing machine by authors Orlowski et al. [35], Chuchala et al. [39] and Sinn et al. [41] and is summarized in Equation (5):
P c T = P c w + P a c + P d u l l = m H p τ γ γ S t Q t p h v c + m H p R S t Q t p v c + P a c + P d u l l ,
Equation (5) consists of four components. The first term ❶ describes the internal work of plasticity along the shear plane, where τγ is shear yield stress and γ is the shear strain along the shear plane. The value of γ can be calculated according to Equation (6), assuming that Φc corresponds to the shear angle:
γ = cos γ f cos ( Φ c γ f ) sin Φ c ,
where: γf—tool side rake angle.
The second component ❷ defines internal work required for the formation of new surface, where R corresponds to the fracture toughness or specific work of material separation. The coefficient of friction correction Q appears in both ❶ and ❷, representing an effect of friction between tool rake face and separated material. Q is computed according to Atkins [47,48] and Orlowski et al. [35] and can be represented as in Equation (7):
Q = 1 sin β μ sin Φ c cos ( β μ γ f ) cos ( Φ c γ f ) ,
where: βμ = tan−1μ is a friction angle (rad) directly related to the coefficient of friction μ.
The other parameters included in the first two terms of Equation (5) are as follow: m—number of saws in the gang, Hp—height of cutting material, St—overall set (kerf width), h—uncut chip thickness, vc—cutting speed.
The power needed for acceleration of chips Pac is also included in Equation (5) as component ❸. It can be described as variation function of the mass flow and tool velocity [35]. It is important to notice that the value of chip acceleration power Pac is calculated globally for the sawing process and is not directly related to the number of working teeth. As the contribution of Pac on the overall cutting power PcT is negligible [41,46], it is not considered in analyses performed for the needs of this research.
The last component ❹ of Equation (5) corresponds to the excessive energy use in case of improper chip formation related to the dullness of the cutting edge. It is an important component of the energetic balance of the real-world processes, and it explains an increase of the cutting forces observed along the tool life and increase of dullness. However, assuring an appropriate sharpness of the tool, component ❹ can be ignored, especially at the initial phase of the tool use. Summarizing, the following assumptions were made for the needs of this research:
  • only freshly sharpened blades were used in cutting tests: component ❹ of Equation (5) = 0;
  • chip acceleration power Pac is omitted: component ❸ of Equation (5) = 0;
  • values of the shear angle Φc were calculated following the approach of Merchant [46], considering specific cutting zone geometry as well as coefficient of friction. This procedure is valid when deviations in the shear angle assigned to inherent material properties can be neglected for larger values of uncut chip thicknesses, as used in this experiment;
  • the value of friction coefficient µ= 0.6 was adopted following the work of Glass and Zielinka [49];
  • the effect of lateral forces on the power consumption can be omitted when cutting in straight direction and with minimal saw deviation [50].
As a consequence, it is possible (by implementing the above listed assumptions) to express Equation (5) as a linear regression function (Equation (8)):
Pcw(h) = c1h + c0,
In that case, c1 and c0 correspond to the slope and intercept, respectively. An independent variable of the regression is the uncut chip thickness h. It has become possible, therefore, to determine values of fracture toughness R and shear yield stress τγ by matching the regression Equation (8) with the experimental data from the cutting tests. A similar approach was reported for diverse materials and cutting kinematics [40,46,51]. However, in the case of the frame sawing process investigated here, values of both fracture toughness R and shear yield stresses τγ are determined for cutting perpendicular to wood fibers direction ΦG–vc (case 90−90 according to Kivimaa [52]). The detailed cutting configuration is presented in Figure 2.
The mathematical procedure for computation of above material properties is expressed in Equations (9) and (10) for τγ and R, respectively:
τ γ = c 1 Q z a γ S t v c ,
R = c 0 Q z a S t v c ,
where za = Hp/tp – number of teeth in contact with the kerf (average).

3. Results and Discussion

Experimental results from the series of cuttings performed on wood samples exposed to different drying procedures are summarized in Figure 3 and Figure 4 for Scots pine and beech, respectively. Each chart presents two test point groups that correspond to the mean value and standard deviations of measured cutting powers at two levels of feed speed vf. Values of feed speed corresponds to the basic geometrical parameter of the cutting process, i.e., uncut chip thickness h. The experimental results are clustered around values h = 0.11 mm and h = 0.22 mm. However, exact values of the feed per tooth (that correspond to h in the case of such gang saws) were determined individually for each processed sample based on recorded experimental data. The data fitting curve (linear regression), as well as regression equation with coefficient c1 and intercept c0, are provided in each chart. A relatively wide range of standard deviations noticed for similar values of feed speed can be observed in both figures. It is related to the high variance of chemical-physical properties native to biomaterials. In the case of experimental wood samples, it was associated to the within batch differences of wood density as well as change of mechanical properties induced by the drying process. Furthermore, the density variation within a single sample, which is associated with the early and late wood differences, presence of wood defects or other common irregularities of the wood tissue, are recorded during cutting tests, increasing even more the scatter of results. Even if all the care was given to assure a homogenous and defect-less set of samples used for different drying experiments, the average density would vary between samples (Table 1). Standard deviation ranges of density within batches in individual test groups corresponding to studied drying methods are summarized in Table 1.
The evident effect of wood density ρ [54,55] can be compensated by normalizing values of resulting mechanical parameters [41]. Diverse algorithms can be implemented for the data unification and further non-biased interpretation. The approach adopted for the need of this study included correction of the measured cutting power by the variation related to the average wood assessed individually for each processed board. Consequently, Equation (5) was revised as follows (Equation (11)), assuming that Pac and Pdull can be neglected for the cutting process on the frame saw:
P c T = P c w = m H p S t Q t p ( τ γ * γ h + R * ) v c ρ ,
Two novel meta variables were introduced here that correspond to normalized fracture toughness R* [J m kg−1] and normalized shear yield stress in the shear zone τ*γ [MPa m3 kg−1]. Equation (11) can, therefore, be expressed as a lineal equation by following the same logic as used for deriving Equation (8). In that case, the cutting power related to the wood density as a function of the chip thickness h is presented in Equation (12):
P c w * ( h ) = P c w ρ ( h ) = c 1 h + c 0 ,
The experimental values of the cutting power per one saw blade Pcw* normalized by density are shown in Figure 5 and Figure 6 for pine and beech, respectively. It is evident that the standard deviation ranges are noticeably reduced compared to the not-normalized results. It is confirmed by higher values of determination coefficients r2 in the majority of investigated cases. This phenomenon confirms the hypothesis of the direct correlation of the density and fracture properties of wood. However, it is also clear that considering density only does not lead to complete elimination of the experimental data scatter. It was previously observed by Chuchala et al. [56] that the effect of density upon the cutting power (cutting forces) is better explained if also considering the synergetic effect of varying cutting conditions, such as feed speed. It is clearly visible in Figure 3, Figure 4, Figure 5 and Figure 6 that the repeatability of results, observed as a minimized standard deviation ranges, is higher for low chip thicknesses.
Average values of R* and τ*γ for Scots pine are summarized in Table 3 together with standard deviations of the observed variation. The fracture toughness increases when wood is dried in other than air drying conditions. The air drying [1] is considered as a reference method that affects the native physical and chemical properties of wood to a small extent [2]. That leads to the conclusion that industrially dried wood becomes more brittle than wood not exposed to elevated temperatures. Conversely, the shear yield stress slightly decreases, making the Scots pine wood easier to deform. A heated-steam experimental drying appears to be the mildest as the fracture parameters were less altered. The relatively high range of the observed values variation is affected by rather low number of cutting tests implemented in this pilot research. Statistical analyses of the experimental results are summarized in Table 4. No significant differences were noticed in the case of pine wood for both the normalized fracture toughness R* and normalized shear yield stress τ*γ.
A very diverse effect of investigated drying techniques on the fracture properties of hardwood (beech) can be derived from data presented in Table 3. While the average shear yield stress (normalized by the wood density) was not large affected by different drying techniques, the fracture toughness evidently dropped when drying in industrial kilns. It is especially noticeable for wood dried according to experimental heat stem approach where R* reached negligible values. The wood of beech consists of a significant share of hemicelluloses that are very sensitive to thermo-hydro modification. It results in extensive hydrolysis when exposed to elevated temperatures in the presence of moisture. Even if the mass loss due to mild thermal modification during drying is negligible, the chemical reconfiguration of hemicelluloses, evidenced as a decrease of hygroscopicity, was observed [15,25,27,57]. It was also reported that thermolysis of hemicelluloses is highly inhibited when drying/thermal modification occurs in vacuum [25,27]. This was confirmed in this research where fracture toughness of vacuum dried wood decreased only slightly. The ANOVA test revealed significant differences between the determined normalized fracture toughness and normalized shear yield stress values for beech wood dried by four different methods (Table 4). However, no statistical differences were detected between air-dried and vacuum-dried beech wood.

4. Conclusions

The research presented allows us to derive the following conclusions:
  • The wood machining process on the frame saws can be adopted for straightforward determination of the fracture properties of materials by monitoring the power needed for cutting.
  • Wood density is an important factor affecting the mechanism of the chip formation and should, therefore, be considered when interpreting cutting test results. A straightforward way for the results normalization was presented in this report.
  • Technique of wood drying affects the physical-chemical properties that are revealed as a change of the fracture toughness and shear yield stress derived from the cutting tests.
  • The observed trend of fracture properties changes was different for examined softwood (Scots pine) and hardwood (beech). Statistically significant differences (α = 0.05) between the determined mechanical properties of wood dried by four different methods were observed only for beech. It is related to differences in microstructure of both species and diverse deterioration mechanisms when exposing wood to elevated temperatures.
  • The methodology presented here can be highly useful for accurate forecasting of the cutting power that is indispensable to optimize construction of woodworking machines and tools.
The results obtained revealed that drying wood with the use of high temperature air-steam mixture significantly reduces the fracture toughness of beech wood. However, it was not noticed for pine wood. Further studies and experimentation are therefore indispensable to properly understand mechanisms of material changes due to drying.

Author Contributions

Conceptualization, K.A.O., D.C. and T.M.; methodology, K.A.O., D.C. and T.M.; investigation, K.A.O., D.C. and T.M.; writing—original draft preparation, D.C., J.S. and K.A.O.; writing—review and editing, D.C., K.A.O., J.S., M.L. and T.O.; formal analysis, D.C., T.M., M.L. and T.O.; supervision, K.A.O., J.S. and M.L.; project administration, K.A.O. and D.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Science Centre, Poland, within the project MINIATURA 2, grant number 2018/02/X/ST8/03399, and Ministry of Science and Higher Education, Poland, for the maintenance of scientific and research equipment—PRW-15M frame saw, grant number 21/E-359/SPUB/SP/2019.

Acknowledgments

The authors would like to acknowledge PhD. Eng. Jacek Baranski from Gdansk University of Technology, Poland, for sharing the drying kiln for research. The authors gratefully acknowledge the European Commission for funding the InnoRenew project [Grant Agreement # 739574] under the Horizon2020 Widespread-Teaming program, the Republic of Slovenia (investment funding from the Republic of Slovenia and the European Union European Regional Development Fund) and infrastructural ARRS program IO-0035. “Delivering fingertip knowledge to enable service life performance specification of wood”—CLICK DESIGN (N° 773324) supported under the umbrella of ERA-NET Cofund Forest Value by Ministry of Education, Science and Sport of the Republic of Slovenia. Forest Value has received funding from the European Union’s Horizon 2020 research and innovation program.

Conflicts of Interest

The authors declare no conflict of interest.

List of Symbols

BA
BC
BV
BS
sample codes for beech woodStoverall set (kerf width)
FSPfiber saturation pointTair temperature
HFsaw frame strokeWsample width
Hpsample heightc0c1coefficients of linear regression function
Lpsample lengthfzfeed per tooth
MCmoisture contenthuncut chip thickness
MOEmodulus of elasticitymsaws number
MORmodulus of rupturenFnumber of strokes of saw frame per min
Pacpower needed for acceleration of chipspair pressure
Pcaverage cutting powertccutting time
PcToverall cutting powertptooth pitch
Pcwaverage cutting power in a working strokevccutting speed
P*cwaverage cutting power in a working stroke normalized by densityvffeed speed
Pdullpower needed to chip formation while dullness of the cutting edgezanumber of teeth in contact with the kerf (average)
Piaverage idle powerβμfriction angle
PTmean total powerγshear strain along the shear plane
PA
PC
PV
PS
sample codes for pine woodγfrake angle
Qcoefficient of friction correctionμfriction coefficient
Rfracture toughnessρdensity
Rfracture toughness in perpendicular directionτγshear yield stress
R*fracture toughness in perpendicular direction normalized by densityτγshear yield stress in perpendicular direction
Rcrelative/compressive strengthτ*γshear yield stress in perpendicular direction normalized by density
RHrelative humidityΦcshear angle

References

  1. Williamson, C. The energy sector: A hidden goliath. In Primary Wood Processing, 2nd ed.; Walker, J.C.F., Ed.; Springer Science and Business Media B.V.: Dordrecht, The Netherlands, 2006; p. 596. [Google Scholar]
  2. Bond, B.H.; Espinoza, O. A Decade of Improved Lumber Drying Technology. Curr. For. Rep. 2016, 2, 106–118. [Google Scholar] [CrossRef]
  3. Espinoza, O.; Bond, B.H. Vacuum Drying of Wood—State of the Art. Curr. For. Rep. 2016, 2, 223–235. [Google Scholar] [CrossRef] [Green Version]
  4. Allegretti, O.; Rémond, R.; Perré, P. Nonsymmetrical drying tests—Experimental and numerical results for free and constrained spruce samples. Dry. Technol. 2018, 36, 1554–1562. [Google Scholar] [CrossRef]
  5. He, Z.; Yang, F.; Peng, Y.; Yi, S. Ultrasound—Assisted vacuum drying of wood: Effects on drying time and product quality. BioResources 2013, 8, 855–863. [Google Scholar] [CrossRef] [Green Version]
  6. Zhang, J.-W.; Liu, H.-H.; Yang, H.; Yang, L. Drying Characteristics of Eucalyptus urophylla × E. grandis with Supercritical CO2. Materials 2020, 13, 3989. [Google Scholar] [CrossRef]
  7. Herrera-Díaz, R.; Sepúlveda-Villarroel, V.; Pérez-Peña, N.; Salvo-Sepúlveda, L.; Salinas-Lira, C.; Llano-Ponte, R.; Ananias, R.A. Effect of wood drying and heat modification on some physical and mechanical properties of radiata pine. Dry. Technol. 2017, 36, 537–544. [Google Scholar] [CrossRef]
  8. Oltean, L.; Teischinger, A.; Hansmann, C. Influence of low and moderate temperature kiln drying schedules on specific mechanical properties of Norway spruce wood. Holz Als Roh Werkst. 2010, 69, 451–457. [Google Scholar] [CrossRef] [Green Version]
  9. Roszyk, E.; Stachowska, E.; Majka, J.; Mania, P.; Broda, M. Moisture-Dependent Strength Properties of Thermally-Modified Fraxinus excelsior Wood in Compression. Materials 2020, 13, 1647. [Google Scholar] [CrossRef] [Green Version]
  10. Borrega, M.; Kärenlampi, P.P. Three mechanisms affecting the mechanical properties of spruce wood dried at high temperatures. J. Wood Sci. 2009, 56, 87–94. [Google Scholar] [CrossRef]
  11. Ayadi, R.; Koubaa, A.; Braghiroli, F.L.; Migneault, S.; Wang, H.; Bradai, C. Effect of the Pyro-Gasification Temperature of Wood on the Physical and Mechanical Properties of Biochar-Polymer Biocomposites. Materials 2020, 13, 1327. [Google Scholar] [CrossRef] [Green Version]
  12. Bekhta, P.; Sedliacik, J.; Bekhta, N. Effect of Veneer-Drying Temperature on Selected Properties and Formaldehyde Emission of Birch Plywood. Polymers 2020, 12, 593. [Google Scholar] [CrossRef] [Green Version]
  13. Bliem, P.; Konnerth, J.; Frybort, S.; Gartner, C.; Mauritz, R.; Van Herwijnen, H.W. Influence of drying and curing parameters on phenol-formaldehyde impregnated wood veneers. J. Adhes. 2019, 96, 253–271. [Google Scholar] [CrossRef]
  14. Klement, I.; Vilkovský, P.; Vilkovská, T. The Effect of Contact-Drying on Physical Properties of European Beech (Fagus sylvatica L.). Forests 2020, 11, 890. [Google Scholar] [CrossRef]
  15. Nakagawa, M.; Umemura, K.; Kawai, S.; Kanayama, K. Influence of drying temperature on NO2 sorption ability of cedar timber. J. Wood Sci. 2020, 66, 1–9. [Google Scholar] [CrossRef]
  16. Vitas, S.; Segmehl, J.S.; Burgert, I.; Cabane, E. Porosity and Pore Size Distribution of Native and Delignified Beech Wood Determined by Mercury Intrusion Porosimetry. Materials 2019, 12, 416. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Arantes, A.K.; De Conti, C.; Sturaro, S.M.E.; De Conti, A.C. Effect of Drying Methods in the Resistance Strength Compression Parallel Test for Bamboo. Adv. Mater. Res. 2015, 1088, 664–668. [Google Scholar] [CrossRef]
  18. Fonseca, E.; Barreira, L. Experimental and Numerical Method for Determining Wood Char-Layer at High Temperatures Due to An Anaerobic Heating. Int. J. Saf. Secur. Eng. 2011, 1, 65–76. [Google Scholar] [CrossRef]
  19. Fonseca, E.; Coelho, D.; Barreira, L. Structural safety in wooden beams under thermal and mechanical loading conditions. Int. J. Saf. Secur. Eng. 2012, 2, 242–255. [Google Scholar] [CrossRef] [Green Version]
  20. Allegretti, O.; Brunetti, M.; Cuccui, I.; Ferrari, S.; Nocetti, M.; Terziev, N. Thermo—Vacuum modification of spruce (Picea abies Karst.) and fir (Abies alba Mill.) wood. BioResources 2012, 7, 3656–3669. [Google Scholar]
  21. Castro, G.; Rosso, L.; Allegretti, O.; Cuccui, I.; Cremonini, C.; Negro, F.; Zanuttini, R. Influence of thermo-vacuum treatment on bending properties of poplar rotary-cut veneer. iForest Biogeosci. For. 2017, 10, 161–163. [Google Scholar] [CrossRef] [Green Version]
  22. Blanchet, P.; Kaboorani, A.; Bustos, C. Understanding Effects of Drying Methods on Wood Mechanical Properties at Ultra and Cellular Levels. Wood Fiber Sci. 2016, 48, 1–12. [Google Scholar]
  23. Báder, M.; Németh, R.; Sandak, J.; Sandak, A. FTIR analysis of chemical changes in wood induced by steaming and longitudinal compression. Cellulose 2020, 27, 1–19. [Google Scholar] [CrossRef]
  24. Nge, T.T.; Tobimatsu, Y.; Yamamura, M.; Takahashi, S.; Takata, E.; Umezawa, T.; Yamada, T. Effect of Heat Treatment on the Chemical Structure and Thermal Properties of Softwood-Derived Glycol Lignin. Molecules 2020, 25, 1167. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Sandak, J.; Sandak, A.; Allegretti, O. Chemical Changes to Woody Polymers Due to High-Temperature Thermal Treatment Assessed with near Infrared Spectroscopy. J. Near Infrared Spectrosc. 2016, 24, 555–562. [Google Scholar] [CrossRef] [Green Version]
  26. Vidholdová, Z.; Sandak, A.; Sandak, J. Assessment of the chemical change in heat treated pine wood by near infrared spectroscopy. Acta Fac. Xylologiae 2019, 61, 31–42. [Google Scholar] [CrossRef]
  27. Sandak, A.; Allegretti, O.; Cuccui, I.; Rosso, L.; Castro, G.; Negro, F.; Cremonini, C.; Zanuttini, R. Thermo-Vacuum Modification of Poplar Veneers and its Quality Control. Bioresources 2016, 11, 10122–10139. [Google Scholar] [CrossRef] [Green Version]
  28. Očkajová, A.; Barcík, Š.; Kučerka, M.; Koleda, P.; Korčok, M.; Vyhnáliková, Z. Wood dust granular analysis in the sanding process of thermally modified wood versus its density. BioResources 2019, 14, 8559–8572. [Google Scholar] [CrossRef]
  29. Očkajová, A.; Kučerka, M.; Kminiak, R.; Krišťák, L.; Igaz, R.; Réh, R. Occupational Exposure to Dust Produced when Milling Thermally Modified Wood. Int. J. Environ. Res. Public Heal. 2020, 17, 1478. [Google Scholar] [CrossRef] [Green Version]
  30. Očkajová, A.; Kučerka, M.; Kminiak, R.; Rogoziński, T. Granulometric composition of chips and dust produced from the process of working thermally modified wood. Acta Fac. Xylologiae 2020, 62, 103–111. [Google Scholar] [CrossRef]
  31. Orłowski, K.A.; Chuchała, D.; Muziński, T.; Barański, J.; Banski, A.; Rogoziński, T. The effect of wood drying method on the granularity of sawdust obtained during the sawing process using the frame sawing machine. Acta Fac. Xylologiae 2019, 61, 83–92. [Google Scholar] [CrossRef]
  32. Nasir, V.; Cool, J. A review on wood machining: Characterization, optimization, and monitoring of the sawing process. Wood Mater. Sci. Eng. 2018, 15, 1–16. [Google Scholar] [CrossRef]
  33. Nasir, V.; Cool, J. Characterization, optimization, and acoustic emission monitoring of airborne dust emission during wood sawing. Int. J. Adv. Manuf. Technol. 2020, 109, 2365–2375. [Google Scholar] [CrossRef]
  34. Nadolny, K.; Kapłonek, W.; Sutowska, M.; Sutowski, P.; Myslinski, P.; Gilewicz, A. Experimental Studies on Durability of PVD-Based CrCN/CrN-Coated Cutting Blade of Planer Knives Used in the Pine Wood Planing Process. Materials 2020, 13, 2398. [Google Scholar] [CrossRef] [PubMed]
  35. Orlowski, K.A.; Ochrymiuk, T.; Atkins, A.; Chuchala, D. Application of fracture mechanics for energetic effects predictions while wood sawing. Wood Sci. Technol. 2013, 47, 949–963. [Google Scholar] [CrossRef] [Green Version]
  36. Orlowski, K.; Ochrymiuk, T.; Atkins, A. An Innovative Approach to the Forecasting of Energetic Effects While Wood Sawing. Drv. Ind. 2014, 65, 273–281. [Google Scholar] [CrossRef] [Green Version]
  37. Orlowski, K.A.; Ochrymiuk, T.; Sandak, J.; Sandak, J. Estimation of fracture toughness and shear yield stress of orthotropic materials in cutting with rotating tools. Eng. Fract. Mech. 2017, 178, 433–444. [Google Scholar] [CrossRef]
  38. Orlowski, K.A.; Ochrymiuk, T.; Hlásková, L.; Chuchala, D.; Kopecky, Z. Revisiting the estimation of cutting power with different energetic methods while sawing soft and hard woods on the circular sawing machine: A Central European case. Wood Sci. Technol. 2020, 54, 457–477. [Google Scholar] [CrossRef] [Green Version]
  39. Chuchala, D.; Ochrymiuk, T.; Orlowski, K.A.; Lackowski, M.; Taube, P. Predicting cutting power for band sawing process of pine and beech wood dried with the use of four different methods. Bioresources 2020, 15, 1844–1860. [Google Scholar] [CrossRef]
  40. Hlásková, L.; Orlowski, K.A.; Kopecký, Z.; Sviták, M.; Ochrymiuk, T. Fracture toughness and shear yield strength determination for two selected species of central European Provenance. Bioresources 2018, 13, 6171–6186. [Google Scholar] [CrossRef]
  41. Sinn, G.; Chuchała, D.; Orlowski, K.A.; Taube, P. Cutting model parameters from frame sawing of natural and impregnated Scots pine (Pinus sylvestris L.). Holz Als Roh Werkst. 2020, 78, 1–8. [Google Scholar] [CrossRef]
  42. Baranski, J. Moisture content during and after high- and normal-temperature drying processes of wood. Dry. Technol. 2017, 36, 751–761. [Google Scholar] [CrossRef]
  43. Barański, J.; Klement, I.; Vilkovská, T.; Konopka, A. High Temperature Drying Process of Beech Wood (Fagus sylvatica L.) with Different Zones of Sapwood and Red False Heartwood. Bioresources 2017, 12, 1861–1870. [Google Scholar] [CrossRef] [Green Version]
  44. Klement, I.; Vilkovská, T.; Baranski, J.; Konopka, A. The impact of drying and steaming processes on surface color changes of tension and normal beech wood. Dry. Technol. 2018, 37, 1490–1497. [Google Scholar] [CrossRef]
  45. Wasielewski, R.; Orlowski, K.A. Hybrid dynamically balanced saw frame drive. Holz Als Roh Werkst. 2002, 60, 202–206. [Google Scholar] [CrossRef]
  46. Orlowski, K.A.; Pałubicki, B. Recent progress in research on the cutting processes of wood. A review COST Action E35 2004–2008: Wood machining—Micromechanics and fracture. Holzforschung 2009, 63, 181–185. [Google Scholar] [CrossRef]
  47. Atkins, A. Modelling metal cutting using modern ductile fracture mechanics: Quantitative explanations for some longstanding problems. Int. J. Mech. Sci. 2003, 45, 373–396. [Google Scholar] [CrossRef]
  48. Atkins, A. Toughness and cutting: A new way of simultaneously determining ductile fracture toughness and strength. Eng. Fract. Mech. 2005, 72, 849–860. [Google Scholar] [CrossRef]
  49. Glass, S.V.; Zelinka, S.L. Moisture Relations and Physical Properties of Wood. In Wood Handbook—Wood as an Engineering Material; Centennial, Ed.; General Technical Report FPL-GTR-190; Department of Agriculture, Forest Service, Forest Products Laboratory: Madison, WI, USA, 2010; p. 508. [Google Scholar]
  50. Mohammadpanah, A.; Hutton, S.G. Dynamics Behavior of a Guided Spline Spinning Disk, Subjected to Conservative In-Plane Edge Loads, Analytical and Experimental Investigation. J. Vib. Acoust. 2016, 138, 041005. [Google Scholar] [CrossRef]
  51. Hlásková, L.; Kopecký, Z.; Solař, A.; Patočka, Z.; Mendel University. Cutting test as source of fracture toughness and shear yield strength for axial-perpendicular model of wood cutting. Wood Fiber Sci. 2019, 51, 58–68. [Google Scholar] [CrossRef]
  52. Kivimaa, E. Die Schnittkraft in der Holzbearbeitung. Holz Als Roh Werkst. 1952, 10, 94–108. [Google Scholar] [CrossRef]
  53. Orlicz, T. Obróbka Drewna Narzędziami Tnącymi, (in Polish: Wood Processing With Cutting Tools), 6th ed.; Skrypty SGGW-AR w Warszawie, Wydawnictwo SGGW-AR: Warszawa, Poland, 1988; p. 503. [Google Scholar]
  54. Kminiak, R.; Gaff, M. Roughness of Surface Created by Transversal Sawing of Spruce, Beech, and Oak Wood. Bioresources 2015, 10, 2873. [Google Scholar] [CrossRef]
  55. Stasiak, M.; Molenda, M.; Bańda, M.; Horabik, J.; Wiącek, J.; Parafiniuk, P.; Wajs, J.; Gancarz, M.; Gondek, E.; Lisowski, A.; et al. Friction and Shear Properties of Pine Biomass and Pellets. Materials 2020, 13, 3567. [Google Scholar] [CrossRef] [PubMed]
  56. Chuchala, D.; Orlowski, K.A.; Sandak, A.; Sandak, J.; Pauliny, D.; Barański, J. The Effect of Wood Provenance and Density on Cutting Forces While Sawing Scots Pine (Pinus sylvestris L.). Bioresources 2014, 9, 5349–5361. [Google Scholar] [CrossRef] [Green Version]
  57. Hermawan, A.; Sakagami, H.; Sobri, S.A.; Amini, M.H.M.; Ramle, S.F.M.; Rasid, S. The effects of drying temperatures on preservative retention and penetration of some Malaysian fast-growing species timbers. Dry. Technol. 2020, 1–10. [Google Scholar] [CrossRef]
Figure 1. Schematic of schedules for dying experimental beech samples with four studied processes. Note: air p—air pressure, air RH—air relative humidity, air T—air temperature, FSP—fiber saturation point.
Figure 1. Schematic of schedules for dying experimental beech samples with four studied processes. Note: air p—air pressure, air RH—air relative humidity, air T—air temperature, FSP—fiber saturation point.
Materials 13 04692 g001
Figure 2. Cutting speed directions when splitting orthotropic materials; axial || cutting along fibers (cutting direction 90°–0°), perpedicular ⊥ cutting across fibers (cutting direction 90°–90°), and cross-cutting ≠ direction (cutting direction 0°–90°) [38,52,53].
Figure 2. Cutting speed directions when splitting orthotropic materials; axial || cutting along fibers (cutting direction 90°–0°), perpedicular ⊥ cutting across fibers (cutting direction 90°–90°), and cross-cutting ≠ direction (cutting direction 0°–90°) [38,52,53].
Materials 13 04692 g002
Figure 3. Relation between cutting power per one saw blade and uncut chip thickness when sawing pine wood dried with different methods: PA—air drying, PC—conventional kiln drying, PS—warm heated-steam experimental drying, PV—vacuum kiln drying.
Figure 3. Relation between cutting power per one saw blade and uncut chip thickness when sawing pine wood dried with different methods: PA—air drying, PC—conventional kiln drying, PS—warm heated-steam experimental drying, PV—vacuum kiln drying.
Materials 13 04692 g003
Figure 4. Relation between cutting power per one saw blade and uncut chip thickness when sawing beech wood dried with different methods: BA—air drying, BS—warm heated-steam experimental drying, BC—conventional kiln drying, BV—vacuum kiln drying.
Figure 4. Relation between cutting power per one saw blade and uncut chip thickness when sawing beech wood dried with different methods: BA—air drying, BS—warm heated-steam experimental drying, BC—conventional kiln drying, BV—vacuum kiln drying.
Materials 13 04692 g004
Figure 5. Relation between normalized (by density) cutting power per one saw blade and uncut chip thickness when sawing pine wood dried with different methods: PA—air drying, PC—conventional kiln drying, PS—warm heated-steam experimental drying, PV—vacuum kiln drying.
Figure 5. Relation between normalized (by density) cutting power per one saw blade and uncut chip thickness when sawing pine wood dried with different methods: PA—air drying, PC—conventional kiln drying, PS—warm heated-steam experimental drying, PV—vacuum kiln drying.
Materials 13 04692 g005
Figure 6. Relation between normalized (by density) cutting power per one saw blade and uncut chip thickness when sawing beech wood dried with different methods: BA—air drying, BC—conventional kiln drying, BS—warm heated-steam experimental drying, BV—vacuum kiln drying.
Figure 6. Relation between normalized (by density) cutting power per one saw blade and uncut chip thickness when sawing beech wood dried with different methods: BA—air drying, BC—conventional kiln drying, BS—warm heated-steam experimental drying, BV—vacuum kiln drying.
Materials 13 04692 g006
Table 1. Density of tested pine and beech wood samples.
Table 1. Density of tested pine and beech wood samples.
Sample CodeDrying MethodWood Density ρ
(Mean Value and Standard Deviation)
[kg·m−3]
Pine Wood
PAmodified air drying515 ± 57.1
PCconventional kiln drying467 ± 26.4
PSheated-steam experimental drying490 ± 33.5
PVvacuum kiln drying474 ± 41.5
Beech Wood
BAmodified air drying723 ± 52.0
BCconventional kiln drying722 ± 36.9
BSheated-steam experimental drying706 ± 25.4
BVvacuum kiln drying757 ± 38.5
Table 2. Settings of the sash-gang saw, and its saw blade as used in the experimental cuttings.
Table 2. Settings of the sash-gang saw, and its saw blade as used in the experimental cuttings.
ParameterSymbolValueUnit
Machine Setting
number of strokes of saw frame per minnF685spm
saw frame strokeHF162mm
number of saws in the gangm5
average cutting speedvc3.69m·s−1
feed speedslowvf10.9m·min−1
fastvf21.9m·min−1
feed per toothslowfz10.11mm
fastfz20.22mm
Tool Setting
the sharp saw blades with stellate tipped teeth
overall set (kerf width)St2mm
saw blade thicknesss0.9mm
free length of the saw bladeL0318mm
blade widthb30mm
tooth pitchtp13mm
tool side rake angleγf9°
tool side clearanceαf14°
tension stresses of saws in the gangσN300MPa
Table 3. Fracture toughness R* and shear yield stress τ*γ normalized by density (average value and standard deviations) for of Scots pine and beech dried with different methods.
Table 3. Fracture toughness R* and shear yield stress τ*γ normalized by density (average value and standard deviations) for of Scots pine and beech dried with different methods.
Sample CodeDrying MethodNormalized Fracture Toughness R*
[J·m·kg−1]
Normalized Shear Yield Stress τ*γ
[MPa·m3·kg−1]
Pine Wood
PAair drying2.72 (±1.66)0.040 (±0.009)
PCconventional kiln drying3.69 (±1.54)0.036 (±0.007)
PSheated-steam experimental drying3.34 (±1.20)0.032 (±0.006)
PVvacuum kiln drying3.59 (±1.08)0.034 (±0.006)
Beech Wood
BAair drying2.42 (±0.65)0.062 (±0.003)
BCconventional kiln drying1.08 (±0.40)0.066 (±0.003)
BSheated-steam experimental drying0.39 (±0.21)0.068 (±0.002)
BVvacuum kiln drying1.70 (±1.03)0.064 (±0.003)
Table 4. Significance of differences between fracture toughness and shear yield stress of pine and beech samples dried with different methods (ANOVA) (α = 0.05).
Table 4. Significance of differences between fracture toughness and shear yield stress of pine and beech samples dried with different methods (ANOVA) (α = 0.05).
Pine Wood
Sample CodeSourceDFAdj SSAdj MSF-Valuep-ValueF-Critical
Fracture Toughness
PAbetween groups34.18581.39530.73470.54042.9604
PCwithin groups2751.27381.8990
PStotal3055.4596
PV
Shear Yield Stress
PAbetween groups30.0002859.51 × 10−51.94360.14642.9604
PCwithin groups270.0013214.89 × 10−5
PStotal300.001606
PV
Beech Wood
Fracture Toughness
BAbetween groups314.28894.76309.62250.000263.0280
BCwithin groups2311.38450.4950
BStotal2625.6734
BV
Shear Yield Stress
BAbetween groups30.0001344.46 × 10−53.17430.043333.0280
BCwithin groups230.0003231.4 × 10−5
BStotal260.000457
BV
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chuchala, D.; Sandak, J.; Orlowski, K.A.; Muzinski, T.; Lackowski, M.; Ochrymiuk, T. Effect of the Drying Method of Pine and Beech Wood on Fracture Toughness and Shear Yield Stress. Materials 2020, 13, 4692. https://0-doi-org.brum.beds.ac.uk/10.3390/ma13204692

AMA Style

Chuchala D, Sandak J, Orlowski KA, Muzinski T, Lackowski M, Ochrymiuk T. Effect of the Drying Method of Pine and Beech Wood on Fracture Toughness and Shear Yield Stress. Materials. 2020; 13(20):4692. https://0-doi-org.brum.beds.ac.uk/10.3390/ma13204692

Chicago/Turabian Style

Chuchala, Daniel, Jakub Sandak, Kazimierz A. Orlowski, Tomasz Muzinski, Marcin Lackowski, and Tomasz Ochrymiuk. 2020. "Effect of the Drying Method of Pine and Beech Wood on Fracture Toughness and Shear Yield Stress" Materials 13, no. 20: 4692. https://0-doi-org.brum.beds.ac.uk/10.3390/ma13204692

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