Measuring density and Young’s modulus of a log through the vibration test without measuring its weight

Abstract

This study examines a simple estimation method for the measurement of the mass of an on-site log through a vibration test. In this method, rather than the log itself, a cut end portion is weighed. For this purpose, the vibration method with additional mass (VAM) was applied to Sitka spruce (Picea sitchensis Carr.) circular truncated cones (model log) and Japanese cedar (Cryptomeria japonica D. Don) logs. Longitudinal vibration tests were performed on the circular truncated cones with/without an additional mass. Furthermore, the cut end portions of the circular truncated cones and logs were used as the virtual additional mass in the VAM. From the results of the vibration test using specimens with/without the concentrated mass, it is possible to estimate the mass of a circular truncated cone with 10% error by the VAM. The cut end portion of a circular truncated cone could be used as the virtual mass in the VAM. From the experimental and theoretical results, to maintain high estimation accuracy, the specimen length must not be too short as shown in our previous study for a specimen with constant cross-sectional shape. The cut end portion of the logs could be used for the virtual mass of the VAM.

Introduction

When logs are bought in log markets, stable qualities of the logs, in other words, appropriate classifying based on log shapes such as diameter, length and warp is required. Since timber quality is affected by its Young’s modulus, classifying logs based on the Young’s modulus enables to enrich lineups of log markets. Evaluating the Young’s modulus of a log before sawmilling allows it to be sold at a high value. However, the Young’s modulus of logs cannot be estimated from the evaluation based on the shape.

In classifying lines at log markets and large-scale sawmilling factories, the shape of logs on a conveyor is measured and the logs are consequently classified. Although there is a classifying process having a scale and a tapping system to measure the mass and Young’s modulus of the logs, the process reduces the production efficiency. Because a log is large and heavy, handling it is difficult, and determining the output of a scale is time consuming. Hence, a method to measure the mass of a log in short time is desirable.

Because logs have an extra length of about 100 mm with respect to the nominal length, use of the extra length portion was examined in this study. It is possible that the density of the cut extra length portion is regarded as the density of a whole log. However, the measurement of the accurate volume of the cut end portion (the extra length portion) will be difficult because the shape of the butt end of a log is irregular due to swelled butt, notches and backcuts in felling. On the other hand, an inexpensive scale with high accuracy that is easy and safe to handle can be used for the cut end portion that is much smaller than the main body of a log. Weighing the cut end portion does not cause a serious reduction of the production efficiency.

A method for measuring the mass, density, and Young’s modulus of a specimen without needing to weigh the specimen was developed based on a frequency equation that incorporates the effect of an additional mass attached to a wooden bar [1,2,3,4,5,6]. This method is referred to as the vibration method with additional mass (VAM) in this study. This method utilizes decreases in the resonance frequency by attaching additional mass. The ratio of the resonance frequency with the additional mass to that without it is used in the frequency equation incorporating the effect of the concentrated mass attached to a specimen and its position. The mass ratio (concentrated mass/specimen) is subsequently calculated using the frequency equation.

To analyze the potential application of the VAM for practical purposes, a series of parameters were investigated, including the suitable mass ratio (additional mass/specimen) [7], the connection between the additional mass and specimen [7], the crossers’ position for the piled lumber [8], the specimen moisture content [9], and the bending vibration generation method [10]. The VAM could perhaps be used to assess the deterioration of the cross beams of timber guardrails [11] and the mass of the piled lumber [12]. The estimation accuracy of the VAM decreased with an increase in the effect of shear and rotatory inertia, and it could be corrected [13].

Our previous study [14] showed that a cut end portion could be used as the virtual mass in the VAM. In this study, we investigate whether the use of a cut end portion in the VAM can be applied to a log. First, the VAM was applied to the circular truncated cone used as the model of a log. Next, a cut end portion was used for the virtual additional mass in the VAM using a circular truncated cone and log.

Theory

Longitudinal vibration for a circular truncated cone with a concentrated mass

A circular truncated cone of radii of r₁ and r₂ (r₁ $\ge$ r₂, k is the ratio of r₂ to r₁, k = r₂/r₁) and length l with concentrated mass M placed at (x₁, x₂) = (al, 0) shown in Fig. 1 is considered. Herein, longitudinal vibration under the free-free condition is considered.

Fig. 1

Circular truncated cone with additional mass

A circular truncated cone is prepared by rotating segments AB and BC about the x axis as shown in Fig. 1. The equations for the segments AB and BC are:

Segment AB:

$${y}_{1}={p}_{1}{x}_{1}+{q}_{1}\left(0\le {x}_{1}\le al\right). $$

(1)

$${p}_{1}=\frac{\left(k-1\right){r}_{1}}{l}. $$

(2)

$${q}_{1}={r}_{1}. $$

(3)

Segment BC:

$${y}_{2}={p}_{2}{x}_{2}+{q}_{2}\left(0\le {x}_{2}\le bl\right). $$

(4)

$${p}_{2}={p}_{1}. $$

(5)

$${q}_{2}={r}_{1}\left\{\left(k-1\right)a+1\right\}. $$

(6)

The longitudinal displacements u₁(x₁, t) and u₂(x₂, t) are represented as

$${u}_{1}\left({x}_{1},t\right)={v}_{1}\left({x}_{1}\right){\mathrm{sin}}{\omega }_{n}t. $$

(7)

$${v}_{1}\left({x}_{1}\right)=\frac{{C}_{11}}{{p}_{1}{x}_{1}+{q}_{1}}{\mathrm{sin}}\left(\sqrt{\frac{\rho }{E}}{\omega }_{n}{x}_{1}+{C}_{21}\right), $$

(8)

$${u}_{2}\left({x}_{2},t\right)={v}_{2}\left({x}_{2}\right){\mathrm{sin}}{\omega }_{n}t, $$

(9)

$${v}_{2}\left({x}_{2}\right)=\frac{{C}_{12}}{{p}_{2}{x}_{2}+{q}_{2}}{\mathrm{sin}}\left(\sqrt{\frac{\rho }{E}}{\omega }_{n}{x}_{2}+{C}_{22}\right), $$

(10)

where C₁₁, C₂₁, C₁₂, and C₂₂ are constants, t, ρ, E, and ω (equal to 2πf, where f is the resonance frequency) are time, density, Young’s modulus, and angular frequency, and the subscript n is the mode number.

Because the axial force does not exist at each extremity under a free-free condition:

$$\left\{\begin{array}{c}{x}_{1}=0: \displaystyle \frac{\partial {u}_{1}}{\partial {x}_{1}}=0\\ \\ \displaystyle {x}_{2}=bl: \frac{\partial {u}_{2}}{\partial {x}_{2}}=0\end{array}\right.$$

(11)

Because both parts of the bar are connected and the difference in the axial force in each bar is equal to the inertia force exerted by the concentrated mass at (x₁, x₂) = (al, 0)

$$\left\{\begin{array}{c}{u}_{1}={u}_{2}\\ ES({x}_{2}) \displaystyle \frac{\partial {u}_{2}}{\partial {x}_{2}}-ES({x}_{1})\frac{\partial {u}_{1}}{\partial {x}_{1}}=M\frac{{\partial }^{2}{u}_{1}}{\partial {t}^{2}}=M\frac{{\partial }^{2}{u}_{2}}{\partial {t}^{2}}\end{array}\right.$$

(12)

where S is cross-sectional area.

From Eqs. (7)–(12)

$$\left\{{\left(k-1\right)}^{2}+k{m}_{n}^{2}\right\}{\mathrm{tan}}{m}_{n}-{m}_{n}{\left(k-1\right)}^{2}+\frac{1}{3}\mu {m}_{n}\frac{{k}^{2}+k+1}{{\left\{\left(k-1\right)a+1\right\}}^{2}}\frac{\left\{{m}_{n}+\left(k-1\right){\mathrm{tan}}a{m}_{n}\right\}\left\{{km}_{n}+\left(k-1\right){\mathrm{tan}}b{m}_{n}\right\}}{1-{\mathrm{tan}}a{m}_{n}{\mathrm{tan}}b{m}_{n}}=0$$

(13)

where µ is the ratio of the concentrated mass to the mass of the circular truncated cone and is defined as

$$\mu =\frac{M}{\frac{1}{3}\uppi \left({r}_{1}^{2}+{r}_{1}{r}_{2}+{r}_{2}^{2}\right)l\rho }=\frac{3M}{\uppi \left({k}^{2}+k+1\right){r}_{1}^{2}l\rho }$$

(14)

Equation (13) is the frequency equation for a circular truncated cone with a concentrated mass.

• i. If µ  = 0, then Eq. (13) is as follows [15]:

  $$\left\{{\left(k-1\right)}^{2}+k{m}_{n}^{2}\right\}{\mathrm{tan}}{m}_{n}-{m}_{n}{\left(k-1\right)}^{2}=0$$

  (15)

• ii. If k = 1, then Eq. (13) is as follows [3]:

  $${\mathrm{sin}}{m}_{n}+\mu {m}_{n}{\mathrm{cos}}{am}_{n}{\mathrm{cos}}{bm}_{n}=0.$$

  (16)

Longitudinal vibration testing method without weighing the specimen

The resonance frequency, represented by f_n0 (0: value without the additional mass), is expressed as follows:

$${f}_{n0}=\frac{{m}_{n0}}{2\uppi l}\sqrt{\frac{E}{\rho }}$$

(17)

where m_n0 is a constant depending on k for a circular truncated cone from Eq. (15).

The resonance frequency is decreased experimentally by attaching the additional mass, whereas, the dimensions, density, and Young’s modulus are not altered. Hence, it can be said that m_n0 changes to m_n. The resonance frequency after attaching the additional mass is expressed as follows:

$${f}_{n}=\frac{{m}_{n}}{2\uppi l}\sqrt{\frac{E}{\rho }}$$

(18)

From Eqs. (17) and (18)

$${m}_{n}=\frac{f_{n}}{f_{n0}}m_{n0}$$

(19)

The measured resonance frequencies f_n0 and f_n are substituted in Eq. (19) to calculate m_n, and the calculated m_n is substituted in Eq. (13) to calculate µ. A specimen mass and density can be obtained by substituting the calculated µ, the concentrated mass, and the dimensions of the bar in Eq. (14). The Young’s modulus can be calculated by substituting the estimated density, resonance frequency without the concentrated mass, and dimensions of the bar in Eq. (17) [1,2,3,4,5,6].

The above steps represent the calculation procedure for the VAM. The mass of the specimen is not required for the calculations.

Materials and methods

Specimens

Air-dried, 1000-mm length, circular truncated cones (larger diameter, smaller diameter) = (100 mm, 90 mm), (100 mm, 80 mm), (100 mm, 70 mm), (100 mm, 60 mm), (100 mm, 50 mm), (100 mm, 40 mm), and (100 mm, 30 mm) were produced from Sitka spruce (Picea sitchensis Carr.) with a cross section of 110 mm square and a length of 1500 mm without a pith. The circular truncated cone is a model of a log and differs from an actual log in that the circular truncated cone has the slope grain. The larger and smaller diameter ends are referred as to butt end and top end, respectively, in this study. Two specimens were made for each dimension. The specimens were conditioned at 20 °C and 65% relative humidity. The vibration tests for the Sitka spruce circular truncated cones were conducted under the same conditions.

In addition, six Japanese cedar (Cryptomeria japonica D. Don) logs with bark having a length of 3000 mm and the larger diameter of about 190—260 mm and the smaller diameter of about 160—210 mm planted in the Forestry and Forest Products Research Institute (Tsukuba, Japan) were used as specimens. Their tree ages were 16, 32, and 41 years. Japanese cedar logs were used within two days after felling. The vibration tests for the Japanese cedar logs were conducted in a sawmill building without controlling temperature and relative humidity. The cut end portions of the logs were weighed before the vibration tests and oven-dried at 105 °C.

Vibration tests

Free-free longitudinal vibration tests were conducted for the following three purposes.

First, to experimentally prove the theory of the VAM for the circular truncated cones mentioned above, the free-free longitudinal vibration tests were conducted on the circular truncated cones with the length of 1000 mm with/without the concentrated mass. The concentrated mass was attached to the butt end (x = 0) or the top end (x = l). Iron plates with dimensions 25 mm × 15 mm × 2 mm, mass 5.4 g (M1) and that with the dimensions of 25 mm × 25 mm × 2 mm, mass 10.0 g (M2) were used for the concentrated masses. The concentrated masses M3, M4, M5, and M6 were made by laminating 2, 3, 4, and 5 M2 iron plates using the double-sided tapes in the thickness direction, respectively.

Second, to experimentally examine the use of cut end portions in the VAM, the free-free longitudinal vibration tests were conducted on the circular truncated cones. One specimen was cut from the butt end and the other was cut from the top end. All the specimens were cut at the distances of 10 mm until the length was 700 mm. The longitudinal vibration tests were performed, while the specimens were shortened stepwise.

Third, to experimentally examine the application of the VAM to the logs, the free-free longitudinal vibration tests were conducted. Three logs were cut from the butt end and the other three were cut from the top end of a specimen. All the logs were cut at the distances of 50 mm until the length was 2850 mm. The longitudinal vibration tests were performed, while the specimen was shortened stepwise.

The free-free longitudinal vibration tests were performed using the following procedure: the test bar was placed on a support near its center point. Next, longitudinal vibration was generated by tapping the RT-plane of one end in the L-direction using a hammer with an iron-head, and the first-mode vibration of the specimen was detected using a microphone (PRECISION SOUND LEVEL METER 2003, NODE Co., Ltd, Tokyo, Japan) placed at the other end. The signal was processed using a Fast Fourier Transform (FFT) digital signal analyzer (Multi-Purpose FFT Analyzer CF-5220, Ono-Sokki, Co., Ltd., Yokohama, Japan) to observe the high-resolution resonance frequencies. The diagram of the experimental setup and the photograph of the vibration test for the logs are shown in Figs. 2 and 3, respectively.

Fig. 2

Diagram of the experimental setup

Fig. 3

Vibration test for the logs

Results and discussion

The mean (standard deviation) density and Young’s modulus of the initial specimens before cutting were 480 (23) kg/m³, 15.45 (0.98) GPa (Sitka spruce), and 878 (121) kg/m³, 5.13 (1.61) GPa (Japanese cedar), respectively. The mean (standard deviation), maximum, minimum moisture content of the cut end portion of Japanese cedar log were 172 (34), 238, and 133%, respectively.

The estimation accuracy of the VAM was expressed using the mass ratio (estimated specimen mass through the VAM/measured specimen mass) M_VAM/M₀ and the estimation accuracy of VAM was considered to be sufficiently high when 0.9 ≤ M_VAM/M₀ ≤ 1.1 in this study. Since the density and Young’s modulus are proportional to mass and density, respectively, the ratio of the estimated density through the VAM to the measured density, and the ratio of the estimated Young’s modulus through the VAM to the Young’s modulus measured through the normal vibration test without the concentrated mass are equal to M_VAM/M₀.

Table 1 shows the accuracy of the VAM for the circular truncated cones using the decrease in the resonance frequency by attaching the concentrated mass. The range of µ was from 0.000156 to 0.0347, and M_VAM/M₀ was from 0.9 to 1.1 in several µ values. Hence, it was possible to estimate the masses of the circular truncated cones with 10% error by the VAM. The estimation accuracy of VAM for the smaller concentrated mass was higher than that for the larger concentrated mass regardless of k. This tendency is explained in the same way as the uniform cross-section specimen using the relationship between m_n and µ calculated from the frequency equation incorporating the effect of the concentrated mass attached to a specimen and its position [3, 5]. m_n decreases monotonically with respect to µ, and the relationship between m_n and µ is convex downward as shown in Fig. 4. Hence, the high estimation accuracy of VAM is likely to be obtained for the larger absolute value of the differential coefficient of the m_n− µ relationship i.e. for the smaller µ because the effect of any errors in m_n caused by measuring the resonance frequency on µ can be ignored. There was not a clear tendency between the estimation accuracy of VAM and k.

Table 1 Estimation acccuracy of the vibration method with additional mass (VAM) M_VAM/M₀ for the circular truncated cones when the concentrated mass was placed at the butt end or the top end

Fig. 4

Relationship m₁ – µ in Eq. (13)

Figures 5 and 6 demonstrate the estimation accuracy of the VAM when the butt end portion or the top end portion of the circular truncated cones was cut. The initial properties (dimensions, mass and resonance frequency) for the length of 1000 mm and properties after cutting correspond to the values with/without the concentrated mass in the section “Longitudinal vibration testing method without weighing the specimen”, respectively. For accuracy, µ was calculated by

Fig. 5

Experimental estimation acccuracy of the vibration method with additional mass (VAM) M_VAM/M₀ (M_VAM: Estimated mass by the VAM, M₀: Mesasured mass) for the circular truncated cones when the cut butt end portion was used for the virtual additional mass

Fig. 6

Experimental estimation acccuracy of the vibration method with additional mass (VAM) M_VAM/M₀ (M_VAM: Estimated mass by the VAM, M₀: Mesasured mass) for the circular truncated cones when the cut top end portion was used for the virtual additional mass

$${\mu }_{\mathrm{M}}=\frac{{M}_{\mathrm{i}}-{M}_{\mathrm{cut}}}{{M}_{\mathrm{cut}}}$$

(20)

where M_i and M_cut are the initial specimen mass, and specimen mass after end cutting, respectively.

m_n0 was obtained by substituting k in Eq. (15). Equation (15) was solved using a software [16].

M_VAM/M₀ decreased when the specimens were shortened. Several results deviated from the M_VAM/M₀—length curve at the beginning of the shortening the specimen because the change in the resonance frequency was small. These tendencies were similar to those of the specimens with constant cross-sectional shapes [14].

The estimation accuracy reduced with the decreasing specimen length (Figs. 5 and 6). The effect of the decrease in the specimen length on the accuracy of estimation was explored theoretically.

The specimens can be divided into two parts: the end portion and the main body. The ratio of the mass (end portion/main body) is obtained. Dimensions of the main body and end portion are shown in Table 2.

Table 2 Dimensions of the main body and the end portion of the circular truncated cones

When the range of $0\le x\le al$ in Fig. 1 is the end portion, i.e. when the butt end is cut, the radius ratio (top end/butt end) of the main body is

$${k}_{\mathrm{EC}}=\frac{{r}_{2}}{{r}_{al}}=\frac{k}{\left(1-k\right)b+k}$$

(21)

where the subscript EC means “end cut”.

The volume ratio (end portion/main body) is

$${\mu }_{\mathrm{v}}=\frac{\frac{\uppi }{3}al\left({r}_{1}^{2}+{r}_{1}{r}_{al}+{r}_{al}^{2}\right)}{\frac{\uppi }{3}bl\left({r}_{al}^{2}+{r}_{al}{r}_{2}+{r}_{2}^{2}\right)}=\frac{1-b}{b}\frac{{\left\{\left(1-k\right)b+k\right\}}^{2}+\left(1-k\right)b+k+1}{{\left\{\left(1-k\right)b+k\right\}}^{2}+\left\{\left(1-k\right)b+k\right\}k+{k}^{2}}.$$

(22)

When the range of $al\le x\le l$ in Fig. 1 is the end portion, i.e. when the top end is cut, the radius ratio (top end/but end) of the main body is

$${k}_{\mathrm{EC}}=\frac{{r}_{al}}{{r}_{1}}=\left(k-1\right)a+1.$$

(23)

The volume ratio (end portion/main body) is

$${\mu }_{\mathrm{v}}=\frac{\frac{\uppi }{3}bl\left({r}_{al}^{2}+{r}_{al}{r}_{2}+{r}_{2}^{2}\right)}{\frac{\uppi }{3}al\left({r}_{1}^{2}+{r}_{1}{r}_{al}+{r}_{al}^{2}\right)}=\frac{1-a}{a}\frac{{\left\{\left(k-1\right)a+1\right\}}^{2}+\left\{\left(k-1\right)a+1\right\}k+{k}^{2}}{{\left\{\left(k-1\right)a+1\right\}}^{2}+\left(k-1\right)a+2}.$$

(24)

Substituting the initial values before cutting the end (resonance frequency: f_ni, constant: m_ni, and length: l_i) for Eq. (18):

$${f}_{n\mathrm{i}}=\frac{{m}_{n\mathrm{i}}}{2\uppi {l}_{\mathrm{i}}}\sqrt{\frac{E}{\rho }}$$

(25)

Substituting the values after cutting the end (resonance frequency: f_ncut, constant: m_ncut, and length: l_cut) for Eq. (17),

$${f}_{n\mathrm{cut}}=\frac{{m}_{n\mathrm{cut}}}{2\uppi {l}_{\mathrm{cut}}}\sqrt{\frac{E}{\rho }}$$

(26)

Substituting f_ni and f_ncut for Eq. (19), m_n to be substituted for Eq. (13) is

$${m}_{n}=\frac{f_{n\mathrm{i}}}{f_{n\mathrm{cut}}}m_{n\mathrm{cut}}=\frac{l_{\mathrm{cut}}}{l_{\mathrm{i}}}m_{n\mathrm{i}}$$

(27)

If the µ calculated by substituting m_n from Eq. (27) for Eq. (13) is written as µ_cut, the accuracy of the VAM is

$$\frac{{M}_{\mathrm{VAM}}}{{M}_{0}}=\frac{{\mu }_{\mathrm{V}}}{{\mu }_{\mathrm{cut}}}.$$

(28)

The ratio M_VAM/M₀ is plotted against the specimen length ratio (after/before end cutting) (Figs. 7 and 8). The experimental estimations had high overlap with the theoretical estimations. To maintain high estimation accuracy, a specimen length must not be too short as shown in our previous study for a specimen with the constant cross-sectional shape [14]. From the theoretical analysis, the estimation accuracy was more than 0.90 when the specimen length ratio was more than 0.855 for the butt end cut and 0.826 for the top end cut.

Fig. 7

Theoretical estimation acccuracy of the vibration method with additional mass (VAM) M_VAM/M₀ (M_VAM: Estimated mass by the VAM, M₀: Mesasured mass) for the circular truncated cone when the cut butt end portion was used for the virtual additional mass

Fig. 8

Theoretical estimation acccuracy of the vibration method with additional mass (VAM) M_VAM/M₀ (M_VAM: Estimated mass by the VAM, M₀: Mesasured mass) for the circular truncated cone when the cut top end portion was used for the virtual additional mass

Table 3 shows the mass of the Japanese cedar logs estimated by the VAM. The calculation procedure is the same as that for the circular truncated cones. It can be said that M_VAM/M₀ was 0.9 or more. Hence, it is possible that the mass of a log can be simply estimated using a cut end portion with the sufficient accuracy. Since the estimation accuracy was sufficiently high when the cutting length was from 50 to 150 mm, the high estimation accuracy will be obtained by cutting the aforementioned extra length.

Table 3 Estimation acccuracy of the vibration method with additional mass (VAM) for logs when the butt end or top end were cut

We consider the application of the VAM to on-site logs. Since m_n0 depends on k, it is necessary to contain a process to calculate an m_n0 or an m_n0—k table in the VAM. Although the logs were placed on a support near their center point, i.e., the nodal position of the first mode of the longitudinal vibration, the logs will be supported by at least two points on-site. Hence, the number of supporting points (the resonance mode of the longitudinal vibration to be used) and materials to be used for the supporting points need to be examined in the future.

Conclusions

Longitudinal vibrations were undertaken and the VAM was applied to the circular truncated cones and logs. The following results were obtained.

• 1. From the results of the vibration test using specimens with/without the concentrated mass, it was possible to estimate the mass of a circular truncated cone with 10% error by the VAM.

• 2. The cut butt end and the cut top end portions of a circular truncated cone could be used as the virtual mass in the VAM.

• 3. From the experimental results and theoretical results, in order to maintain high estimation accuracy, a specimen length must not be too short as shown in our previous study for a specimen with the constant cross-sectional shape.

• 4. The cut butt end and the cut top end portions of a log could be used as the virtual mass in the VAM. Hence, it is possible to simply measure the mass of log without weighing the log itself on site through a vibration test.