Original Article

Determination of Factors Influencing the Coefficient of Internal Friction of Natural and Stained Oak

Vladimir SHAMAEV1,https://orcid.org/0009-0002-1696-6152, Alexandr RUSSU1, Ilya MEDVEDEV1, Nikolay TRUBNIKOV1, Varvara DRUZYANOVA2, Natalia ZADRAUSKAITE3, Sergey REVYAKO4
Author Information & Copyright
1Department of Wood Science, Voronezh State University of Forestry and Technologies named after G.F. Morozov, Voronezh 394087, Russia
2Department of Operation of Road Transport and Auto Repair, Northeastern Federal University named after M.K. Ammosov, Yakutsk 677000, Russia
3Department of Timber Industry and Materials Processing, Northern (Arctic) Federal University named after M.V. Lomonosov, Arkhangelsk 163000, Russia
4Department of Environmental Engineering Machines, Novocherkassk Engineering and Reclamation Institute, Don State Agrarian University, Novocherkassk 346428, Russia
Corresponding author: Vladimir SHAMAEV (e-mail: shamaev122@outlook.com)

Copyright 2025 The Korean Society of Wood Science & Technology. This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Jun 13, 2024; Revised: Jul 17, 2024; Accepted: Aug 24, 2024

Published Online: Jan 25, 2025

ABSTRACT

This study aims to determine the influence of density, temperature, humidity, and processing time in ultrasonic and pulsed magnetic fields on the wood of natural and stained oak on the magnitude of their internal friction. The logarithmic decrement of damping, obtained by free flexural vibrations, was used to assess the magnitude of internal friction. As a result of the performed research, it was established that the most significant influence on the increase in the logarithmic decrement magnitude of natural oak is exerted by ultrasonic treatment (34%–41%), followed by an increase in humidity (31%–37%), treatment with pulsed magnetic field (21%–29%), with temperature increase having the least effect (6%–7%). As for stained oak, the most significant influence on the increase in the logarithmic decrement magnitude is exerted by treatment with a pulsed magnetic field (17%–42%), followed by an increase in humidity (10%–41%), ultrasonic treatment (11%–37%), with temperature increase having the least effect (6%–9%). Mathematical dependencies have been developed that allow for the prediction of the magnitude of the logarithmic decrement of damping of natural and stained oak wood, depending on density, temperature, and humidity, as well as processing time with ultrasonic and pulsed magnetic fields.

Keywords: internal friction; natural oak; stained oak; logarithmic decrement of damping; ultrasonic; pulsed magnetic field

1. INTRODUCTION

Internal friction in wood, also known as damping, is the general term for mechanisms converting elastic energy into heat, playing a role as a damping factor. Alongside the speed of sound, internal friction is considered one of the vibrational properties of wood and depends on its orientation (longitudinal, radial, and tangential; Afzalinia and Roberge, 2007; Radhamani and Balakrishnan, 2023; Weller and Wert, 1983).

Scientific literature indicates a substantial number of publications dedicated to investigating internal friction in wood (Ha et al., 2023; Kalasee et al., 2023). Studies have been conducted to explore the relationships between the characteristics of natural wood and internal stresses, viscoelastic and vibrational properties, hysteresis damping, the dependency of vibrational properties on anatomical structure, and geometric, and mechanical properties of wooden structures (Alkadri et al., 2018; Brémaud et al., 2012; Clair, 2012; Grigorev et al., 2014; Kunickaya et al., 2019a, 2019b; Labonnote et al., 2013; Lukina et al., 2022). Additionally, attempts have been made to study the frequency characteristics and modes of natural vibrations and correlations between wood moisture content and damping behavior (Göken, 2018; Ouis, 2022). Efforts are underway to adapt the findings in the field of internal friction studies to develop non-destructive density determination techniques, enhance the mechanical efficiency of wood materials, create metamaterials, design wooden structures with high mechanical energy dissipation, increase resource conservation in wood panel production, and enhance key consumer values in the production and operation of frictional products and tribological setups (Fitzgerald et al., 2021; Gao et al., 2017; Gonzalez et al., 2022; Yin et al., 2022).

The findings indicated that the internal friction tanδ value exhibited a threefold higher vulnerability to electromagnetic compatibility (EMC) alterations in comparison to the E′/γ value (Brémaud and Gril, 2021). The radial (R) properties, namely tanδR, showed much greater sensitivity than the longitudinal (L) properties, resulting in a noticeable increase in anisotropy as the EMC rose. The impact of EMC on the elastic modulus (E′/γ) in the L direction was quite minor when compared to the inherent fluctuations in wood. Nevertheless, the impact of EMC on the damping factor, tanδ, was as substantial as the inherent fluctuations seen in high-quality spruce. The vibrational properties displayed hysteresis concerning relative humidity (RH) but showed little hysteresis concerning equilibrium moisture content. The relationship between tanδ and EMC is significantly affected by the precise duration required for stabilization to happen after equilibrium moisture content is reached.

The empirical data obtained from the research conducted by academics provides evidence supporting the existence of a certain moisture threshold that limits the occurrence of internal friction (Senalik and Farber, 2021). The internal friction experiences a rise on both sides of the minimum value as the moisture content approaches zero or reaches the fiber saturation limit. The temperature is a clear determinant of the moisture content at which the lowest internal friction is seen. Likewise, there exist certain temperatures at which the occurrence of internal friction is reduced, and the determination of the lowest level of internal friction is contingent upon the moisture content. An increase in the temperatures at which little internal friction occurs is seen when the quantity of moisture drops. In the case of wood with high moisture content, there exists a consistent inclination for the internal friction to diminish with rising temperatures.

The mechanical characteristics of sawdust and wood chips were measured at five different moisture levels (Stasiak et al., 2015), which is valuable information for process design and management. The relationship between the moisture content of the material and the elastic modulus has been firmly established. The elastic modulus of dry wood chips was measured to be 0.33 MPa, but at 50% humidity, it was found to be 0.25 MPa. The mean angle of internal friction for sawdust was 27°, whereas for wood chips it was roughly 33°. The angle of internal friction ranged from 34° for sawdust to 42° for wood chips. An inverse correlation was detected between the internal friction angle and the moisture content of sawdust. The authors showed that the angles of internal friction from wood were higher than from straw (Swietochowski et al., 2018).

Table 1 demonstrates that the wooden chips have the highest angles in comparison with wooden pellets (Wu et al., 2011).

Table 1. Internal friction angles
Material of wood Size (mm) Angle (°)
Pellets < 6 37
< 8 41
< 12 35
Chips < 20 45
< 40 47
< 100 47
Download Excel Table

The research was undertaken by the authors to investigate the internal friction angle and friction coefficient of various surface types, with rapeseed serving as the model material (Xu et al., 2019). A direct shear apparatus was used for the purpose of their inquiry. The researchers furthermore developed prediction models via the process of fitting the empirical data. The rapeseed exhibited an internal friction angle ranging from 24 to 35 degrees. During normal stress conditions, it exhibited a decline and an increase in response to moisture levels.

The research measured the dynamic elastic modulus and internal friction of beech wood at various equilibrium moisture levels within a temperature range of −130°C to 25°C (Sellevold et al., 1975). The response curves show a notable displacement at around –90°C for wood samples with moisture levels equivalent to or above the fiber saturation threshold. The transition is analogous to that seen in other microporous materials that can absorb water, such as solidified cement paste and permeable glass. Lowering the moisture content in the wood reduces the magnitude of the transition and leads to its occurrence at higher temperatures. We hypothesize that the shift is caused by the gradual solidification, sometimes referred to as the “glass transition,” of the water that has been adsorbed.

In the study Russu et al. (2023), prepared samples of natural and modified wood, cut from the stem part of birch, and were subjected to ultrasonic and pulsed magnetic field treatments in radial and tangential directions. The highest level of approximation was achieved in testing samples of pendulous birch natural wood (R2 = 0.98) for the radial direction.

Analysis of existing studies indicates that research in the direction of studying internal friction in wood is being actively pursued. At the same time, it is notable that there is a lack of comprehensive and precise understanding of the physical nature of dissipation and relaxation processes in wood, hindering the development of physical and mathematical descriptions of regularities. This is relevant both theoretically and practically, as it allows for establishing correlations between the properties of internal material structures, obtaining new insights into processes within the substance structure, and their influence on the main macroscopic material characteristics.

The objective of this study is to obtain experimental data and mathematical dependencies of the magnitude of internal friction on factors such as density, temperature, humidity, ultrasonic treatment, and pulsed magnetic field, as well as to develop mathematical models that adequately describe them, to find optimal solutions.

For this research, natural and stained oak wood is chosen as the experimental material. The selection of natural and stained oak wood is motivated by limited data available on these materials in the specific segment of the scientific and technical literature, as well as their undeniable value for the end consumer.

2. MATERIALS and METHODS

2.1. Materials

For the experiment, specimens of healthy natural wood (Quercus robur) were used, obtained from the Voronezh region, Russian Federation. Specimens of stained oak were obtained from wood harvested in the Medveditsa River (a tributary of the Oka River) in the Tver region (40 km from the city of Tver), Russian Federation.

The physical and mechanical properties of the wood species used are presented in Table 2.

Table 2. Mechanical properties of wood
Type Air-dried density ρ (kg/m3) Impact strength (kJ/m2) Tensile strength static bending (MPa)
Moisture Moisture
12% 30% 12% 30%
Natural oak (Quercus robur) 625 76 65 103 66
Stained oak wood 775 38 32 58 35
Download Excel Table

The appearance of the specimens was assessed visually. Quality control was conducted during specimen preparation, ensuring that specimens with cracks, knots, decay, and other defects were not included. For the experiment, specimens of the following geometric dimensions were cut: length 500 mm, width 42 mm, thickness 2.5 mm.

The surface roughness of wood specimens can significantly affect their mechanical properties, including damping. To minimize the impact of surface roughness on experimental results, it is essential to process the specimens to achieve a consistent and acceptable level of roughness. We used fine sandpaper to smooth the surface of the wood specimens [start with a coarser grit (120) and progress to finer grits (220) to achieve a smooth finish]. After processing, we inspected the surface of the specimens using surface roughness testers to ensure the desired roughness level (Ra below 10 micrometers) according to ISO 21920-2:2021 Geometrical product specifications (GPS) – Surface texture: Profile.

Prior to testing, the specimens were weighed, and their density was determined. A minimum of three specimens were required for each condition. The visual appearance of the specimens is illustrated in Fig. 1.

wood-53-1-33-g1
Fig. 1. Appearance of samples of the tested wood. 1, 2: oak natural, 3, 4: oak stained. The top-bottom or bottom-top direction corresponds to the 1: oak natural radial direction, 2: oak natural tangential direction, 3: oak stained radial direction, 4: oak stained tangential direction. Left-right is the lengthwise direction of the lumber, and the wood fibers are in the longitudinal direction.
Download Original Figure
2.2. Methods

When determining internal friction using the method of free decaying flexural vibrations, the following relationship is employed:

Q 1 = δ / π = 1 π ln ( A n A n + 1 )
(1)

Where Q-1 represents internal friction, δ denotes the logarithmic decrement of damping, and An and An+1 are the amplitudes of neighboring vibration periods.

Thus, the magnitude of internal friction and the logarithmic decrement of damping can be equivalently utilized to describe experimental dependencies obtained through the method of free flexural vibrations.

The BING device was used (Traore et al., 2010). The test rests on two rigid supports (nylon threads 17 cm long and 0.35 mm in diameter) or elastic supports (Fig. 2). Whenever possible, the supports are placed on the vibration nodes of the fundamental frequency. The specimen is impacted at one end with a 13-mm-diameter steel ball weighing 9 g (bending or compression vibration depending on the orientation of the impact). A microphone is placed at a distance of 2 cm from the other end of the specimen. The microphone is unidirectional, type EMU4535 (Lem Industries, Bucine, Italy). The analog-to-digital converter is a Picoscope 3224 (Pico Technology, Eaton Socon, UK). It possesses a 12-bit resolution. The antialiasing filter was set to the Nyquist-Shannon frequency. Recordings were not performed in an anechoic room.

wood-53-1-33-g2
Fig. 2. Experimental setup of bending (transversal) and compression (longitudinal) free-vibration device.
Download Original Figure

The experimental setup for generating free decaying oscillations is depicted in Fig. 3. The temperature and humidity of the working environment for conducting experiments were maintained at 22°C and 68%, respectively.

wood-53-1-33-g3
Fig. 3. Laboratory equipment. 1: computer, 2: stabilizer, 3: carrier, 4: magnet, 5: dynamometer, 6: sample, 7: three-position tool, 8: electro transformer, 9: electro corrector, 10: ampermeter, 11: protective tool, 12: digital desktop.
Download Original Figure

The development of the laboratory equipment was carried out based on the findings of the study (Brancheriau et al., 2010). In this research, the authors were able to demonstrate that the type of support significantly influences the measurement of temporary damping; measurements of temporary damping obtained from flexural free vibration are linearly correlated with measurements obtained from compressional vibration.

The sample is securely fastened cantilever-style using a three-position clamping fixture on a sturdy base. Bending occurs in the horizontal plane. To induce and establish initial bending stress, an electromagnet was utilized (Fig. 3, position 4). The response of the wooden sample to the magnetic field is facilitated by a thin strip (0.02 mm thick), weighing no more than 0.6 g, made of permalloy (alloy 79NM) with magnetic-soft properties and high coercive force. The strip is attached to a small portion of the free end of the sample (Fig. 3, position 6). The magnitude of the initial force required to induce free decaying bending oscillations, ranging from 10 to 50 N, is controlled by the magnetic field of the electromagnet (Fig. 3, position 4). The distance from the fastening point to the measurement point is 35 mm and can vary from 0 to 100 mm. The parameters measured during the experiment are listed in Table 3.

Table 3. Measurement parameters during the experiment
Parameter Step Interval
Bending force 10 N 0–50 N
Temperature 20°C 20°C–100°C
Moisture 10% 0%–40%
Duration of ultrasonic treatment 300 s 0–1,200 s
Duration of treatment with pulsed magnetic field 30 s 0–120 s
Download Excel Table

The technical specifications of the ultrasonic bath are provided in Table 4.

Table 4. The technical specifications of ultrasonic bath
Name Unit Value
Operating frequency kHz 24.0 ± 7.5%
Ultrasonic power W 80
Download Excel Table

The technical specifications of the pulsed magnetic field generator are presented in Table 5.

Table 5. Generator specifications pulsed magnetic field
Name Unit meas. Meaning
Maximum magnetic field induction T 0.3
Pulse current kA 2
Number of pulses during processing 3,000 / 6,000
Pulse frequency Hz 100
Pulse duration 10−6 s 10

T: tesla, kA: kiloAmpere, Hz: hertz, s: second.

Download Excel Table

Thus, the testing methods, measuring instruments, and devices used in the experiment comply with the necessary regulatory documents, and the maximum relative measurement error did not exceed 1%.

The experiment is based on a full factorial design. We will seek a mathematical description of the experiment in the form of polynomials, taking into account the nonlinearity of factor interactions. In general, when the number of varied factors is equal to k, the model of the object has the form:

y = b 0 + i = 1 k b i x i + i , j = 1 ; i < j k b i j x i x j + i , j , l = 1 ; i < j < l k b i j l x i x j x l + + x 1 x 2 · · x k
(2)

Where bm…n are the theoretical regression coefficients.

In this case, 5 factors will be used for variation: x1: density (ρ), x2: temperature (T), x3: humidity (w), x4: ultrasonic treatment (u), x5: pulsed magnetic field (m).

For 5 factors, the number of regression coefficients is 21. Table 6 presents the natural and formalized factors.

Table 6. Natural and normalized factors
Natural values Normalized values
F-r i Xi min Xi max Xi 0 xi min xi max xi 0
ρ 1 625 kg/m3 775 kg/m3 700 kg/m3 −1 1 0
T 2 293K 373K 333K −1 1 0
w 3 0% 30% 15% −1 1 0
u(t) 4 0 s 300 s 600 s −1 1 0
m(t) 5 0 s 30 s 60 s −1 1 0

F-r: factor, ρ: density, T: temperature, w: humidity, u(t): ultrasonic treatment, m(t): pulsed magnetic field.

Download Excel Table

Let’s determine the sample size (the number of experiment repetitions) n:

n = t 2 s 2 Δ 2
(3)

Where t is the Student's coefficient, s2 is the sample variance, Δ is the magnitude of the maximum error.

Taking into account the degrees of freedom f = 9 and the confidence probability p = 0.95, we obtain the significance level q = 0.05 and the value of Student’s coefficient t = 2.262. Let’s assume s2> = 5·10−5 and set Δ = 5·10−2. Then, the calculation yields n = 10, 23..., which, when rounded to the nearest integer, gives the minimum number of repetitions for each experiment (the minimum sample size) as 10.

To conduct a regression analysis, the following requirements must be met: homogeneity of variances (Cochran’s criterion), i.e., the absence of any variances that significantly exceed all others; absence of gross errors (Grubbs’ criterion); homogeneity of variances ensures the orthogonality and rotatability of the design.

In this study, a one-sided Grubbs’ criterion (Smirnov-Grubbs criterion) was utilized, indicating the relative deviation of the tested quantity y from its mean value in units of the SD (Göken, 2018):

G r = | y y ¯ | s
(4)

Where Gr is the Grubbs’ criterion, y is the tested value, ȳ is the sample mean, s is the sample SD.

According to the verification results, the condition should be met:

G r < G r max
(5)

Where Grmax is the theoretical (tabulated) value of the Grubbs’ criterion.

At a significance level of 0.05 and a sample size of 10, the value of Grmax = 2.176.

The hypothesis testing for the homogeneity of variances (dispersion characteristics) can be conducted using Cochran’s criterion since the number of experiment repetitions in each case is the same:

G = s max 2 i = 1 N s i 2
(6)

To accept the hypothesis, it is necessary to meet the requirement:

G < G 0 , 05 ; 9 ; 32
(7)

At a significance level of 0.05, with 32 experiments and a sample size of each experiment being 10, the value of G0.05;9;32 = 0.0902.

After verifying the homogeneity of variances, we proceed to compute the variance of experiment reproducibility S2y:

s y 2 = 1 N i = 1 N s i 2
(8)

Subsequently, utilizing the obtained experimental data, the model coefficients bi…n are computed using known formulas (Göken, 2018).

Following the calculation of the model coefficients, their significance is assessed using Student’s criterion:

| b i | > Δ b i = t 0.05 ; f s y 2 N
(9)

After calculating the model coefficients and verifying their significance, the adequacy variance is determined using the Fisher F criterion Fад:

s a d 2 = 1 N l i = 1 N ( y ¯ i y i * ) 2
(10)
F a d = s a d 2 s y 2 < F 0.05 ; N l ; N l
(11)

Where l is the number of significant coefficients in the regression equation, ȳi is the mean response value in the i-th experiment, y*i is the calculated response value according to the regression equation in the i-th experiment.

3. RESULTS and DISCUSSION

3.1. Experimental investigation

As the investigated response function model involves five factors, we first evaluate the role of each factor separately using functions in a two-dimensional plane.

The density factor is represented by two values of natural (625 kg/m3) and steamed (775 kg/m3) oak wood. Clearly, the graph represents a straight line. Visualization of the experimental data for the other four factors is presented in Figs. 48.

wood-53-1-33-g4
Fig. 4. Influence of bending force on δ. 1: bog oak tangential direction, 2: bog oak radial direction, 3: natural oak radial direction, 4: natural oak tangential direction.
Download Original Figure

In the initial stage, it was necessary to determine at what external force magnitude, based on the sample sizes and the strength limit of the corresponding species of natural or modified wood, further experimental work should be conducted. Therefore, measurements of the logarithmic decrement of damping δ were performed for all samples in the range of external force from 10 N to 50 N.

Fig. 4 clearly shows that at the magnitude of 30 N, there is a change in the curvature slope of the graph. Attention is drawn to the fact that between the changes in the δ value in the radial and tangential directions, as seen in Fig. 4, a certain similarity is established. This is associated with the increase in stress in the sample, resulting in approaching the limit of mechanical strength.

For natural oak, the increase in δ in the radial direction was approximately 20%, while in the tangential direction, it was around 30%. For steamed oak, the pattern in the radial direction was similar to that of natural oak, with an increase in δ of about 23% (Fig. 4), whereas in the tangential direction, it was significantly lower, around 11%. The curve trend for stained oak in the tangential direction indicated that besides mechanical factors, the change in the microstructure also influenced the magnitude of the logarithmic decrement of damping. A comparison of the graphs of natural oak and steamed oak (Fig. 4) revealed that for the former, the graph in the radial direction was higher than in the tangential direction, while for the latter, it was the opposite. The reasons for this phenomenon lie in the change in the microstructure and molecular composition of steamed oak, particularly the increased presence of iron.

A characteristic feature observed in all wood samples used in the experiment is the attainment of a stable plateau region in δ at an external bending force magnitude of 40 N or higher (Fig. 4). The most distinctive inflection is observed at an external bending force of 30 N. Since no sharp transitions in the magnitude of δ were observed within the range of external bending forces for all wood samples, it can be inferred that the variation of this quantity is continuous.

Based on the analysis of the curves in Fig. 4, an external bending force magnitude of 40 N was selected for further investigation.

When assessing the influence of the temperature field (Fig. 5), for natural oak, the increase in δ in the radial direction was approximately 7%, while in the tangential direction, it was around 6%. For steamed oak samples, the increase in δ in the radial direction was approximately 6%, while in the tangential direction, it was around 9%. Notably, the tangential direction of steamed oak exhibited the highest magnitude of δ compared to other samples.

wood-53-1-33-g5
Fig. 5. Effect of temperature on δ. 1: bog oak tangential direction, 2: bog oak radial direction, 3: natural oak radial direction, 4: natural oak tangential direction.
Download Original Figure

In the case of increasing positive temperatures, the thermal effect on wood leads to an expansion in volume, thereby serving as a factor in the growth of internal stresses. It was experimentally determined that there is also a slight increase in the magnitude of δ. This was the observed pattern for the specimens investigated by us.

Moisture content significantly influences the technical properties of wood. For quantitative assessment, we utilize the RH of wood expressed in percentages. In this study, humidity varied within the range of 0% to 40%. The experimental results are presented in Fig. 6. For all specimens, a significant increase in the magnitude of δ and a transition to a “plateau” was observed at a humidity level of 30% or higher. It is evident that at 30%, this represents the point of hygroscopic moisture, hence further increases in humidity do not affect the mechanical properties of wood.

wood-53-1-33-g6
Fig. 6. Influence of moisture on δ. 1: bog oak radial direction, 2: bog oak tangential direction, 3: natural oak radial direction, 4: natural oak tangential direction.
Download Original Figure

When assessing the influence of humidity (Fig. 6), for natural oak, the growth of δ in the radial direction was approximately 37%, and in the tangential direction, it was about 31%. For samples of steamed oak, the growth of δ in the radial direction was approximately 41%, and in the tangential direction, it was around 10%, which is nearly three times less than that for natural oak wood.

Thus, the influence of humidity on the magnitude of the logarithmic decrement of damping is more significant than the influence of the temperature field. Moreover, in some cases, the relative increase reaches a threefold growth. Only for the tangential direction of steamed oak wood does the influence of temperature and humidity contribute approximately equally to the overall effect.

Senalik and Farber (2021) describe the moisture level at which internal friction is minimal. Internal friction increases on both sides of this minimum as moisture content decreases to zero or reaches the fiber saturation point. Temperature determines moisture content and minimum internal friction. The minimum internal friction at a moisture content is approximately 6% at 23°C; at –20°C, it is 14%; and at 70°C, it is 4%. The minimum temperature at 90°C is not defined and is observed in the absence of moisture. There are also temperatures of minimum internal friction and temperatures of maximum internal friction, which vary depending on moisture content. As moisture content decreases, the temperature of minimum internal friction increases. Internal friction significantly increases at temperatures above 0°C and moisture content above 10%, accompanied by a strong positive interaction with moisture content. Internal friction decreases with increasing temperature for very dry wood. The logarithmic decrement of hot moist wood is approximately 0.1, while for hot dry wood, it is less than 0.02. Cold wood, regardless of moisture content, will be acceptable.

Brémaud and Gril (2021) conducted experiments on two different wood kinds, namely softwood (spruce) and hardwood (maple), in both the longitudinal (L) and radial (R) orientations. The R properties, namely tanδR, showed much greater sensitivity (2–3 times) than the L properties, resulting in a noticeable increase in anisotropy as the EMC rose.

The propagation of ultrasonic waves in a medium ensures the periodicity of compression and rarefaction processes of the medium’s components. The influence of ultrasound was found to be quite significant for all tested samples (Fig. 7). Considering the sample processing technology, the moisture content of the samples was maintained at 20 ± 1%. Additional measurements of processing duration were conducted for 20 minutes to determine the behavior of the logarithmic decrement curve, taking into account the maximum duration of ultrasonic treatment, which was 15 minutes. The experiment revealed that when the duration of ultrasonic treatment exceeded 15 minutes, the curves reached a plateau.

wood-53-1-33-g7
Fig. 7. Influence of the ultrasonic field on δ. 1: bog oak radial direction, 2: bog oak tangential direction, 3: natural oak radial direction, 4: natural oak tangential direction.
Download Original Figure

When assessing the influence of ultrasound (Fig. 7), for natural oak, the growth of δ in the radial direction was approximately 41%, and in the tangential direction, it was approximately 34%. For samples of bog oak, the corresponding growth of δ in the radial direction was approximately 37%, and in the tangential direction, it was approximately 11%, which is nearly three times less than for natural oak wood. This relationship is analogous to the influence of moisture discussed above. It should be noted that the ultrasonic effect had a slightly greater impact on the magnitude of δ in the samples compared to moisture. This is due to the physics of the wave process, which significantly affects the macroscopic and microscopic characteristics of the medium, forming ordered structures in a certain way.

The graphs illustrating the influence of weak pulsed magnetic fields (PMF) on the magnitude of δ are presented in Fig. 8.

wood-53-1-33-g8
Fig. 8. Influence of a pulsed magnetic field on δ. 1: bog oak radial direction, 2: bog oak tangential direction, 3: natural oak tangential direction, 4: natural oak radial direction.
Download Original Figure

The magnetic field was applied in a series of 1,500 to 6,000 symmetric triangular pulses with a duration of τ = 10 μs, an amplitude of B0 = 0.3 T (tesla), and a repetition frequency of 10 ms (milliseconds). During the exposure, the samples were positioned such that the wood fibers were parallel to the direction of the field.

A characteristic feature for all samples was the plateauing of the δ curves within 60 seconds after the onset of magnetic field treatment.

When assessing the influence of the pulsed magnetic field (Fig. 8), for natural oak, the growth of δ in the radial direction was approximately 21%, and in the tangential direction, it was around 29%. For samples of stained oak, the growth of δ in the radial direction was approximately 42%, nearly twice as much as for natural oak, and in the tangential direction, it was about 17%.

The pulsed magnetic field, by affecting the internal microstructure of the wood, significantly influences the increase in the magnitude of the logarithmic decrement of damping within the first 60 seconds of treatment. The greatest enhancing effect on δ (~42%) was observed in samples of the radial direction of stained oak. Thus, overall, the influence of the pulsed magnetic field on the radial direction is noticeably greater than on the tangential direction.

Table 7 presents the percentage increase values of the logarithmic decrement of damping for natural and stained oak wood in the radial and tangential directions when varying temperature, humidity, as well as processing time under ultrasonic and pulsed magnetic fields.

Table 7. Increase in the magnitude (% relative) of the logarithmic decrement of damping with variations in different factors
Indicator Natural oak Stained oak
Rad. Tang. Rad. Tang.
Temperature (20°C–100°C) 7 6 6 9
Humidity (0%–40%) 37 31 41 10
Ultrasonic treatment (0–1,200 s) 41 34 37 11
Pulsed magnetic field treatment (0–120 s) 21 29 42 17

Rad.: radial, Tang.: tangential.

Download Excel Table

The data presented in Table 7 indicate that the greatest influence on the increase in the value of the logarithmic decrement of natural oak is exerted by ultrasonic treatment (34%–41%), followed by an increase in humidity (31%–37%), treatment with a pulsed magnetic field (21%–29%), and lastly, temperature elevation (6%–7%). As for the stained oak, the greatest influence on the increase in the value of the logarithmic decrement is attributed to treatment with a pulsed magnetic field (17%–42%), followed by an increase in humidity (10%–41%), ultrasonic treatment (11%–37%), and lastly, temperature elevation (6%–9%).

The empirical dependencies obtained by us are consistent with the results of the study by Russu et al. (2023), where it was established that humidity, as well as treatment in ultrasonic and pulsed magnetic fields, have the most significant impact on increasing the value of the logarithmic decrement of natural and modified birch.

Ultrasonic and pulsed magnetic field treatments significantly impact the internal structure and channels of wood, thereby affecting its internal friction. Ultrasonic treatment increases the logarithmic decrement of damping by 34%–41% for natural oak through microstructural changes that enhance viscoelastic properties. Pulsed magnetic field treatment similarly increases damping by 17%–42% for stained oak, potentially altering its magnetic and structural properties. These treatments modify the wood's internal channels and structure, converting elastic energy into heat, and are influenced by environmental factors like humidity and temperature, as demonstrated through mathematical models predicting internal friction based on these variables.

3.2. Development of regression models

Based on the proposed methodology, assessments of the obtained experimental data were conducted using the Grubbs' criterion. The experimental value of the Grubbs’ criterion, Gr = 2.116 < Grmax = 2.176. This indicates that there are no significant outliers in the experimental data.

The hypothesis of homogeneity of variances was tested using Cochran’s criterion: the experimental value G = 0.0468 < G0,05;32;10 = 0.0865. Therefore, the hypothesis of homogeneity of variances is accepted.

The reproducibility variance of the experiment is S2y = 8.999·10−3. The significance level is Δbi = 28.8·10−3.

After performing calculations and considering the significance level, we obtain the model - the final normalized regression equation for the logarithmic decrement of damping δr, (where r denotes the radial direction and, t the tangential direction), linearly associated with the internal friction magnitude:

δ r = π Q r 1 = 2.912 0.237 x 1 + 0.057 x 2 + 0.379 x 3 + 0.435 x 4 + 0.253 x 5 0.042 x 1 x 4 0.031 x 3 x 4 0.033 x 1 x 2 x 5 + 0.040 x 1 x 3 x 4 0.034 x 1 x 2 x 3 x 4
(12)

Since the number of significant coefficients in the regression equation is l = 11, we find F0.05;21;10 = 2,764. Adequacy check Fad=4.2521038.999103=0.472<F0.05;21;10 has demonstrated that with a probability of p = 0,95, it can be asserted that the regression equation adequately approximates the response function.

Similarly, we obtain the second model for the tangential direction. All the data are compiled into a comprehensive Table 8.

Table 8. Regression models
Radial direction Tangential direction
p 0.95 0.95
Gr 2.116 < 2.176 2.033 < 2.176
No major misses No major misses
G 0.0468 < 0.0865 0.0444 < 0.0865
The hypothesis of homogeneity of variances Accepted Accepted
S 2 y 8.999·10−3 9.535·10−3
Δbi 28.8·10−3 31.6·10−3
δ δ r = π Q r 1 = = 2.912 0.237 x 1 + + 0.057 x 2 + 0.379 x 3 + + 0.435 x 4 + 0.253 x 5 0.042 x 1 x 4 0.031 x 3 x 4 0.033 x 1 x 2 x 5 + 0.040 x 1 x 3 x 4 0.034 x 1 x 2 x 3 x 4 δ τ = π Q τ 1 = 2.559 + 0.051 x 1 + + 0.074 x 2 + 0.304 x 3 + + 0.356 x 4 + 0.131 x 5 + + 0.078 x 1 x 2 + 0.073 x 1 x 5 0.063 x 1 x 2 x 5 + 0.063 x 1 x 3 x 5 0.054 x 1 x 2 x 3 x 4
l 11 11
Fad F a d = 4.252 10 3 8.999 10 3 = 0.472 < 2.764 F a d = 7.704 10 3 9.535 10 3 = 0.808 < 2.764
The model is adequate The model is adequate

p: probability, Gr: Grubbs’ criterion, G: Cochran’s criterion, S2y: the reproducibility variance of the experiment, Δbi: the significance level, δ: the logarithmic decrement of damping, l: the number of significant coefficients in the regression equation, Fad: Fisher F criterion.

Download Excel Table

Thus, experimental dependencies in the form of graphs and regression models have been obtained to describe the physico-mechanical interaction between external factors and wood material.

A comparative analysis of the regression coefficients of models for the radial direction of wood and the tangential direction revealed that density exerts an inverse influence in the radial direction compared to the tangential direction, with a difference of 4.6 times in absolute terms. Density as a factor is inconvenient for controlling the wood diagnostics process since it is an intrinsic material parameter. Nevertheless, its consideration is necessary, as comparing its magnitude for the radial direction reveals that it is only 6% less than the coefficient of influence of the impulse magnetic field factor. The temperature factor for the tangential direction is significantly higher by 30% than for the radial direction. The moisture factor is more significant for the radial direction than for the tangential direction by 25%. The ultrasonic field factor is more significant for the radial direction than for the tangential direction by 22%. The impulse magnetic field factor is more significant for the radial direction than for the tangential direction by 93%, nearly twice as much.

At the same time, comparing the coefficients reflecting the joint influence of factors, we notice that on average, this group effect is 86% greater in the tangential direction than in the radial direction. Investigating the reasons for the relationship between the joint action of two or more factors is a promising task that goes beyond the scope of our study at this stage and will be considered in the future.

Since the experiment was conducted under identical conditions, it is possible to rank the factors according to their degree of significance and technological feasibility. The greatest influence on the test samples was exerted by humidity and the ultrasonic field, followed by a slightly lesser influence from the impulse magnetic field. The thermal field had the weakest impact. Therefore, ultrasonic and impulse magnetic fields can be proposed as the most technologically feasible control parameters, while humidity, temperature, and density can be suggested as additional parameters.

Thus, it can be inferred that employing the internal friction magnitude as a criterion for assessing the structural quality of wood represents a promising idea.

Our findings align with the results of other researchers’ studies dedicated to exploring the relationship between the dissipative, acoustic, and mechanical characteristics of wood for the production of musical instruments (Hossen et al., 2018; Obataya et al., 2020; Shirmohammadi et al., 2021; Sproßmann et al., 2017). Extensive documentation exists about the impact of the microfibril angle inside the layer of the tracheid, specifically in relation to the L direction, on the internal friction value and Young’s modulus of wood (Ono and Norimoto, 1983). In addition, it should be noted that our research is consistent with the results of recent research in this field of deep manufacturing of wood (Gong et al., 2021; Hwang and Oh, 2021, 2024).

Internal friction in wood, the conversion of elastic energy into heat during deformation, is crucial for its processing and utilization. High internal friction enhances damping properties, essential for applications like musical instruments, where it reduces unwanted resonance, and in construction materials, where it dampens vibrations, improving structural stability. In furniture and sports equipment, higher internal friction means better shock absorption and longevity. Understanding and controlling internal friction through treatments like ultrasonic and pulsed magnetic fields, as explored in this article, optimizes wood’s mechanical properties and processing efficiency, reducing tool wear and energy consumption. This research provides valuable methods to enhance wood’s performance for specific applications, making it more reliable and effective across various industries, thereby highlighting its practical significance.

4. CONCLUSIONS

The study aimed to investigate the influence of density, temperature, humidity, and processing time under ultrasonic and pulsed magnetic fields on the internal friction of natural and stained oak wood. The internal friction magnitude was determined using the logarithmic decrement obtained by the free flexural vibration method. A full factorial design was employed for experimental planning. The study revealed that the logarithmic decrement for natural oak was δt = 1.40 for the tangential direction and δt = 1.34 for the radial direction, while for stained oak, it was δt = 1.59 and δt = 1.78, respectively. Ultrasonic treatment (34%–41%), humidity increase (31%–37%), pulsed magnetic field processing (21%–29%), and temperature elevation (6%–7%) had the most significant impact on increasing the logarithmic decrement of natural oak. Pulsed magnetic field processing (17%–42%), humidity increase (10%–41%), ultrasonic treatment (11%–37%), and temperature elevation (6%–9%) had the greatest influence on increasing the logarithmic decrement of stained oak.

Mathematical dependencies have been developed that allow predicting the magnitude of the logarithmic decrement of attenuation in natural and stained oak wood based on density, temperature, humidity, and processing time under ultrasonic and pulsed magnetic fields.

Further research should focus on increasing the number of examined samples of various species of both natural and modified (not only stained) wood, as well as on investigating wood with different defects, primarily cracks.

CONFLICT of INTEREST

No potential conflict of interest relevant to this article was reported.

ACKNOWLEDGMENT

Not applicable.

REFERENCES

1.

Afzalinia, S., Roberge, M. 2007. Physical and mechanical properties of selected forage materials. Canadian Biosystems Engineering 49: 2.23-2.27.

2.

Alkadri, A., Carlier, C., Wahyudi, I., Gril, J., Langbour, P., Brémaud, I. 2018. Relationships between anatomical and vibrational properties of wavy sycamore maple. IAWA Journal 39(1): 63-86.

3.

Brancheriau, L., Kouchade, C., Brémaud, I. 2010. Internal friction measurement of tropical species by various acoustic methods. Journal of Wood Science 56(5): 371-379.

4.

Brémaud, I., El Kaïm, Y., Guibal, D., Minato, K., Thibaut, B., Gril, J. 2012. Characterisation and categorisation of the diversity in viscoelastic vibrational properties between 98 wood types. Annals of Forest Science 69: 373-386.

5.

Brémaud, I., Gril, J. 2021. Moisture content dependence of anisotropic vibrational properties of wood at quasi equilibrium: Analytical review and multi-trajectories experiments. Holzforschung 75(4): 313-327.

6.

Clair, B. 2012. Evidence that release of internal stress contributes to drying strains of wood. Holzforschung 66(3): 349-353.

7.

Fitzgerald, D., Sinha, A., Miller, T.H., Nairn, J.A. 2021. Axial slip-friction connections for cross-laminated timber. Engineering Structures 228: 111478.

8.

Gao, S., Wang, X., Wiemann, M.C., Brashaw, B.K., Ross, R.J., Wang, L. 2017. A critical analysis of methods for rapid and nondestructive determination of wood density in standing trees. Annals of Forest Science 74: 1-13.

9.

Göken, J., Fayed, S., Schäfer, H., Enzenauer, J. 2018. A study on the correlation between wood moisture and the damping behaviour of the tonewood spruce. Acta Physica Polonica A 133(5): 1241-1260.

10.

Gong, D.M., Shin, M.G., Lee, S.H., Byeon, H.S., Park, H.M. 2021. Dynamic property of cross-laminated woods made with temperate seven species. Journal of the Korean Wood Science and Technology 49(5): 504-513.

11.

Gonzalez, S., Chacra, E., Carreño, C., Espinoza, C. 2022. Wooden mechanical metamaterials: Towards tunable wood plates. Materials & Design 221: 110952.

12.

Grigorev, I.V., Grigorev, G.V., Nikiforova, A.I., Kunitckaia, O.A., Dmitrieva, I.N., Khitrov, E.G., Pásztory, Z. 2014. Experimental study of impregnation birch and aspen samples. BioResources 9(4): 7018-7026.

13.

Ha, Y.S., Lee, H.J., Lee, S.J., Shin, J.A., Song, D.B. 2023. A study on floor impact sound insulation performance of cross-laminated timber (CLT): Focused on joint types, species and thickness. Journal of the Korean Wood Science and Technology 51(5): 419-430.

14.

Hossen, M.F., Hamdan, S., Rahman, M.R. 2018. Investigation of the acoustic properties of chemically impregnated Kayu malam wood used for musical instrument. Advances in Materials Science and Engineering 2018: 7829613.

15.

Hwang, J.W., Oh, S.W. 2021. Bending strength of board manufactured from sawdust, rice husk and charcoal. Journal of the Korean Wood Science and Technology 49(4): 315-327.

16.

Hwang, J.W., Oh, S.W. 2024. Density profile and sound absorption capability of ceramics manufactured from sawdust, chaff and charcoal: Effect of carbonization temperature and mixing ratio. Journal of the Korean Wood Science and Technology 52(3): 234-242.

17.

Kalasee, W., Lakachaiworakun, P., Eakvanich, V., Dangwilailux, P. 2023. Sound absorption of natural fiber composite from sugarcane bagasse and coffee silver skin. Journal of the Korean Wood Science and Technology 51(6): 470-480.

18.

Kunickaya, O.A., Shadrin, A.A., Burmistrova, O.N., Markov, O.B., Gasparyan, G.D., Davtyan, A.B., Lapshina, M.L., Sleptsova, N.A., Ustinova, V.V., Kruzhilin, S.N. 2019a. Wood treatment with hydro impact: A theoretical and experimental study. Bulgarian Journal of Agricultural Science 25(Suppl. 2): 158-166.

19.

Kunickaya, O.A., Shadrin, A.A., Kremleva, L.V., Mueller, O.D., Ivanov, V.A., Bederdinova, O.I., Kruchinin, I.N., Burgonutdinov, A.M., Zakharova, O.I., Struchkov, N.A. 2019b. Modeling of the processes of the modification of the current volume warming by drainage and pressing. Bulgarian Journal of Agricultural Science 25(Suppl. 2): 167-177.

20.

Labonnote, N., Rønnquist, A., Malo, K.A. 2013. Modified hysteretic damping model applied to Timoshenko timber beams. Computers & Structures 121: 22-31.

21.

Lukina, A., Lisyatnikov, M., Martinov, V., Kunitskya, O., Chernykh, A., Roschina, S. 2022. Mechanical and microstructural changes in post-fire raw wood. Architecture and Engineering 7(3): 44-52.

22.

Obataya, E., Zeniya, N., Endo-Ujiie, K. 2020. Effects of seasoning on the vibrational properties of wood for the soundboards of string instruments. The Journal of the Acoustical Society of America 147(2): 998-1005.

23.

Ono, T., Norimoto, M. 1983. Study on Young’s modulus and internal friction of wood in relation to the evaluation of wood for musical instruments. Japanese Journal of Applied Physics 22(4R): 611-614.

24.

Ouis, D. 2022. Experimental modal analysis of a palm tree log under radial vibrational excitation. Wood and Fiber Science 54(1): 1-13.

25.

Radhamani, R., Balakrishnan, M. 2023. NiTi shape memory alloy: Unraveling the role of internal friction in passive damping: A review. Materials Today Communications 37: 107276.

26.

Russu, A.V., Shamaev, V.A., Razinkov, E.M., Zimelis, A. 2023. Internal friction investigation of the natural and compressed birch (Betula pendula Roth) wood. Forestry Engineering Journal 13(1): 236-256.

27.

Sellevold, E.J., Radjy, F., Hoffmeyer, P., Bach, L. 1975. Low temperature internal friction and dynamic modulus for beach wood. Wood and Fiber Science 7(3): 162-169.

28.

Senalik, C.A., Farber, B. 2021. Mechanical Properties of Wood. In: Wood Handbook: Wood as an Engineering Material, Ed. by Wiemann, W.C. U.S. Department of Agriculture, Forest Service, Forest Products Laboratory, Madison, WI, USA.

29.

Shirmohammadi, M., Faircloth, A., Redman, A. 2021. Assessment of sound quality: Australian native hardwood species for guitar fretboard production. European Journal of Wood and Wood Products 79(2): 487-497.

30.

Sproßmann, R., Zauer, M., Wagenführ, A. 2017. Characterization of acoustic and mechanical properties of common tropical woods used in classical guitars. Results in Physics 7: 1737-1742.

31.

Stasiak, M., Molenda, M., Bańda, M., Gondek, E. 2015. Mechanical properties of sawdust and woodchips. Fuel 159: 900-908.

32.

Swietochowski, A., Dabrowska, M., Lisowski, A. 2018. Friction properties of pellets made of wood and straw. Engineering for Rural Development 23: 1815-1820.

33.

Traore, B., Brancheriau, L., Perre, P., Stevanovic, T., Diouf, P. 2010. Acoustic quality of vene wood (Pterocarpus erinaceus Poir.) for xylophone instrument manufacture in Mali. Annals of Forest Science 67: 815.

34.

Weller, M., Wert, C.A. 1983. Internal friction of coal and of other natural macromolecular solids. Journal de Physique Colloques 44(C9): C9-191-C9-196.

35.

Wu, M.R., Schott, D.L., Lodewijks, G. 2011. Physical properties of solid biomass. Biomass and Bioenergy 35(5): 2093-2105.

36.

Xu, Q., Cheng, X., Chen, X. 2019. Models for predicting frictional properties of rapeseed. International Agrophysics 33(1): 61-66.

37.

Yin, W., Lu, H., Zheng, Y., Tian, Y. 2022. Tribological properties of the rotary friction welding of wood. Tribology International 167: 107396.