# Development of Testing Platform and Comparison Studies for Wood Nondestructive Testing

Authors: Jian Li; Yiming Fang; Jiyong Tang; Hailin Feng; Xiongwei Lou
Journal: Journal of Electrical and Computer Engineering
 (2016)

---

## Abstract

Stress wave based techniques have been developed for evaluating the quality of the wooden materials nondestructively. However the existing techniques have some shortcomings due to the significant variation of the wood properties and are now in need of updating. There are also stress wave based instruments which have been widely used for nondestructive testing of wood. But most of them are inflexible and unsuitable for the tentative studies. This paper proposed and implemented a wood nondestructive testing platform based on NI virtual instrument. Three wood nondestructive testing methods, including peak time interval measurement, cross-correlation, and spectrum analysis, were also tested on this platform with serious decay sample, early decay sample, and defect-free sample. The results show that new methods can be verified easily and the researches of wood nondestructive testing will be accelerated with the designed platform.


---

## Body

## 1. Introduction

Wood nondestructive testing technology can predict the wood internal condition without destruction and has been utilized to detect the mechanical properties of standing trees [1] or the decay and holes of wood [2]. By far, many wood nondestructive testing techniques have been studied, such as X-ray detection [3], ultrasonic [4], mechanical stress [5], and acoustic emission [6].Among them, the stress wave detection method attracts more attention due to its insensitivity to the shape and the size of the wood, harmlessness to human body, and easiness to use [7]. Ross et al. showed that the flight time of the stress wave can be used to predict whether the wood was defective [8]. Li et al. investigated the stress wave velocity pattern in the cross sections of black cherry trees and then employed this pattern to detect the internal defect [9]. All of the above work have proven the effectiveness of stress wave technique for the detection of wood defect.During the past decades, many commercial instruments sprang up, such as Fakopp TreeSonic acoustic device [10] and Model 239A Stress Wave Timer [11]. Other instruments have also been reported [12, 13]. However, the functions of those instruments are always predefined by manufacturers and cannot be altered to satisfy the user’s different requirements. In fact, the propagation of the stress wave is quite complex due to the anisotropy and inhomogeneity of the wood material, and most instruments like the above-mentioned cannot satisfy the practical requirement due to the unobtrusive precisions. Therefore, efforts should be made to investigate more advanced signal processing methods after the stress wave signals are collected. So a flexible testing platform is needed to meet this demand in the pilot phase of the tentative study.The NI virtual instrument has a major advantage that the hardware and the software can be redesigned or improved easily and quickly to meet the variable demand [14]. In this study, we design a wood nondestructive testing platform using the NI virtual instrument. Subsequently, we verify the testing platform by comparing three methods which have potentialities to be applied to wood nondestructive testing in the future.


## 2. Platform for Nondestructive Testing of Wood

### 2.1. Hardware Setup

The hardware system of the testing platform was composed of piezoelectric acceleration sensors, charge amplifier, data acquisition card, and computer. The structure can be illustrated in Figure1. The piezoelectric acceleration sensors were employed to sense the stress wave and fixed on the wood by a nail. The output signal of the sensor was amplified by the charge amplifier, and then the signal was input into NI data acquisition card. Finally, the analogy signal was converted into digital form and sent to computer through universal serial bus (USB).Figure 1
The structure of hardware system.It is very important to select reasonable sensors in the implementation of nondestructive testing wood. A good example is that the frequency band should be wider than that of the stress wave. In this work, the piezoelectric acceleration sensors were selected from Beidaihe Rand Technology Co. Ltd., which has good performance with its wide broadband, low impedance output, antijamming ability, and low noise. And the matched charge amplifier was also manufactured by this company.The Data Acquisition Equipment was selected in terms of the performance in collection efficiency, reliability of the data collection, and convenience of the field use. The equipment NI USB-6259 is a multifunction DAQ module with a high-speed USB interface, which can provide high accuracy at high-speed sampling. The Plug and Play function also makes it very easy to install the hardware system. The MX driver and measurement services software of NI-DAQ provide a simple configuration and programming interface. Together with the DAQ Assistant, NI USB-6259 can help users to save the development time.


### 2.2. Implementation of Signal Acquisition and Data Storage

In this work, the flexible testing platform provides a variety of analyzing methods. What is more, the system can be easily redesigned with functions altered. Therefore, we employed the PC-DAQ which, based on Graphical Programming, can facilitate the hierarchical and modular design.As mentioned above, both the maximum amplitude and the maximum frequency of the output signal were within the ranges of the sensors. Therefore, we connected the amplifier to the ports of the DAQ card directly. Ai0 and Ai1 were set at the sampling channels. The amplitude range was −10 V to 10 V. The differential mode was selected. The sampling frequency was 500 KHz.The signal can be saved in many file formats, such as excel file, text file, and binary file. Among them, the binary file has the fastest accessing speed and can meet the requirement of high-speed data acquisition. Besides, it holds high performance on compatibility, which enables the data to be analyzed in other software. Therefore, all the signals recorded were saved as binary files.The Labview code for implementing signal acquisition and data storage is shown in Figure2.Figure 2
Labview code of signal acquisition and data storage.


## 2.1. Hardware Setup

The hardware system of the testing platform was composed of piezoelectric acceleration sensors, charge amplifier, data acquisition card, and computer. The structure can be illustrated in Figure1. The piezoelectric acceleration sensors were employed to sense the stress wave and fixed on the wood by a nail. The output signal of the sensor was amplified by the charge amplifier, and then the signal was input into NI data acquisition card. Finally, the analogy signal was converted into digital form and sent to computer through universal serial bus (USB).Figure 1
The structure of hardware system.It is very important to select reasonable sensors in the implementation of nondestructive testing wood. A good example is that the frequency band should be wider than that of the stress wave. In this work, the piezoelectric acceleration sensors were selected from Beidaihe Rand Technology Co. Ltd., which has good performance with its wide broadband, low impedance output, antijamming ability, and low noise. And the matched charge amplifier was also manufactured by this company.The Data Acquisition Equipment was selected in terms of the performance in collection efficiency, reliability of the data collection, and convenience of the field use. The equipment NI USB-6259 is a multifunction DAQ module with a high-speed USB interface, which can provide high accuracy at high-speed sampling. The Plug and Play function also makes it very easy to install the hardware system. The MX driver and measurement services software of NI-DAQ provide a simple configuration and programming interface. Together with the DAQ Assistant, NI USB-6259 can help users to save the development time.


## 2.2. Implementation of Signal Acquisition and Data Storage

In this work, the flexible testing platform provides a variety of analyzing methods. What is more, the system can be easily redesigned with functions altered. Therefore, we employed the PC-DAQ which, based on Graphical Programming, can facilitate the hierarchical and modular design.As mentioned above, both the maximum amplitude and the maximum frequency of the output signal were within the ranges of the sensors. Therefore, we connected the amplifier to the ports of the DAQ card directly. Ai0 and Ai1 were set at the sampling channels. The amplitude range was −10 V to 10 V. The differential mode was selected. The sampling frequency was 500 KHz.The signal can be saved in many file formats, such as excel file, text file, and binary file. Among them, the binary file has the fastest accessing speed and can meet the requirement of high-speed data acquisition. Besides, it holds high performance on compatibility, which enables the data to be analyzed in other software. Therefore, all the signals recorded were saved as binary files.The Labview code for implementing signal acquisition and data storage is shown in Figure2.Figure 2
Labview code of signal acquisition and data storage.


## 3. Tentative Methods for Nondestructive Testing of Wood

In the earlier studies, the flight time (or the propagation velocity) of the stress wave has usually been used to predict the internal conditions of the wood, such as Young’s modulus [15] and internal decay [16]. Therefore, it plays a key role in wood nondestructive testing to measure the propagation time of the stress wave accurately. In this study, we employed peak time interval measurement method and cross-correlation method to accurately obtain propagation time of stress wave in wood.The spectral analysis has also been proven to be one of the most effective techniques for detecting internal defect of wood nondestructively. Yu et al. evaluated the wood by analyzing the frequency spectrum of stress wave captured by transmitting and receiving sensors on the wood. The result showed that higher frequency components of the stress wave will be attenuated more quickly in wood, and the stress wave will be attenuated more quickly in decay wood than normal wood [17].


### 3.1. Peak Time Interval Measurement Method

In this work, the time interval between the first peak of transmitted signal and received signal was used as the propagating time of stress wave in the wood. However, there are many problems which hinder the implementation of this method. For example, the peak is prone to misjudgment due to the impulse noise and the weak amplitude of the signal. Therefore, the Median Filtering VI, provided by the Labview software, was firstly employed to reduce the impulse noise. Then, a Peak Detector VI was developed to locate the peaks. The threshold was set as 0.2 V. Finally, the appearing times of the peaks were sent to the Subtraction VI and the flight time of the stress wave can be computed. Figure3 illustrated the Labview code of flight time measurement.Figure 3
Labview code of peak time interval measurement.


### 3.2. Cross-Correlation Method

Cross-correlation method was often used to estimate the signal propagating time [18]. An accurate evaluation can be obtained especially when the two signals are very similar.Assuming that the transmitted signal and the received signal arex
1
(
k
) and x
2
(
k
), respectively, where k is the sampling index. The flight time τ
^ is given by the following function:(1)
τ
^
=
arg
⁡
m
a
x
m
R
x
y
m
,where R
x
y
(
m
) is the cross-correlation function of x
1
(
k
) and x
2
(
k
) and can be expressed as follows:(2)
R
x
y
m
=
1
N
-
m
∑
k
=
0
N
-
m
-
1
x
1
k
x
2
k
-
m
,where N defines the length of both x
1
(
k
) and x
2
(
k
).Because the noise will diminish the similarity of two stress wave signals, we firstly employed the IIR filter VI to eliminate noise. The pass-band used here is 100 Hz to 5 KHz, which is the same as in the literature [17]. Subsequently, the cross-correlation VI was used to compute the cross-correlation value of the two input signals. Finally, the Statistics Express VI was used to locate the maximal correlation value. The flight time was then computed using the location of the maximal correlation. Figure 4 depicts the implementation of this method.Figure 4
Labview code of correlation method.


### 3.3. Analysis of Stress Wave Signal in Frequency Domain

At first, the Create Waveforms VI was developed to cut the stress wave signal data into continuous segments, with the duration being 1 ms. Then the Power Spectrum VI was employed to transform the stress wave signal into the frequency domain. The Labview implementation is shown in Figure5.Figure 5
Labview code of spectrum analysis.


## 3.1. Peak Time Interval Measurement Method

In this work, the time interval between the first peak of transmitted signal and received signal was used as the propagating time of stress wave in the wood. However, there are many problems which hinder the implementation of this method. For example, the peak is prone to misjudgment due to the impulse noise and the weak amplitude of the signal. Therefore, the Median Filtering VI, provided by the Labview software, was firstly employed to reduce the impulse noise. Then, a Peak Detector VI was developed to locate the peaks. The threshold was set as 0.2 V. Finally, the appearing times of the peaks were sent to the Subtraction VI and the flight time of the stress wave can be computed. Figure3 illustrated the Labview code of flight time measurement.Figure 3
Labview code of peak time interval measurement.


## 3.2. Cross-Correlation Method

Cross-correlation method was often used to estimate the signal propagating time [18]. An accurate evaluation can be obtained especially when the two signals are very similar.Assuming that the transmitted signal and the received signal arex
1
(
k
) and x
2
(
k
), respectively, where k is the sampling index. The flight time τ
^ is given by the following function:(1)
τ
^
=
arg
⁡
m
a
x
m
R
x
y
m
,where R
x
y
(
m
) is the cross-correlation function of x
1
(
k
) and x
2
(
k
) and can be expressed as follows:(2)
R
x
y
m
=
1
N
-
m
∑
k
=
0
N
-
m
-
1
x
1
k
x
2
k
-
m
,where N defines the length of both x
1
(
k
) and x
2
(
k
).Because the noise will diminish the similarity of two stress wave signals, we firstly employed the IIR filter VI to eliminate noise. The pass-band used here is 100 Hz to 5 KHz, which is the same as in the literature [17]. Subsequently, the cross-correlation VI was used to compute the cross-correlation value of the two input signals. Finally, the Statistics Express VI was used to locate the maximal correlation value. The flight time was then computed using the location of the maximal correlation. Figure 4 depicts the implementation of this method.Figure 4
Labview code of correlation method.


## 3.3. Analysis of Stress Wave Signal in Frequency Domain

At first, the Create Waveforms VI was developed to cut the stress wave signal data into continuous segments, with the duration being 1 ms. Then the Power Spectrum VI was employed to transform the stress wave signal into the frequency domain. The Labview implementation is shown in Figure5.Figure 5
Labview code of spectrum analysis.


## 4. Experiments and Results

To test the performance of the developed platform, experiments were carried out on threeCamphor samples, as shown in Figure 6. The average diameter of the wood samples was 30 cm and the average height was 40 cm. TheCamphor samples were characterized by different sizes of defects. The first had serious decay area. The second one had early decay area and the last one was a sound wood.Figure 6
Photos of the wood samples used in this work. (a–c): (a) serious defect, (b) early defect, and (c) sound sample.
(a)
(b)
(c)Typical stress waves are shown in Figure7. The blue dashed line plots the transmitted signal and the red solid line corresponds to the received signal.Figure 7
Transmitted and received signals.


### 4.1. Results of Time Domain Analysis

The peak time interval measurement method and the cross-correlation method mentioned above were carried on with our platform. Table1 lists the flight time obtained from the three samples using different methods. In order to reduce the influence of some random factors, the flight time listed in Table 1 is the average of 10 repetitions.Table 1
Flight time obtained by two methods (μs).
Serious defect
Early defect
Defect-free
Peak method
504
160
48
Correlation method
428
144
42It can be seen from Table1 that the flight time obtained by two methods are very close. They are both sensitive to the degree of the defect. It took about 45 μs for stress wave to travel in the defect-free sample and approximately 3 times longer to travel in the sample with early defect. The time of flight corresponding to the seriously sample is 9 times of those in response to the sound sample. Therefore, both of the two methods have the potentiality of application in nondestructive testing of wood and further studies are needed to demonstrate their performance.It also can be observed that the value obtained employing correlation method is a little less than that employing peak method. Literature [19] reported that current time of flight techniques provide stiffness measurements that are higher than those obtained employing other techniques. Therefore, it can be expected that the correlation method has the potential to improve the accuracy of stiffness testing employing time of flight.


### 4.2. Results of Frequency Domain Analysis

With the Labview code of the spectrum analysis, we transformed the signal recorded on the received sensor into frequency domain. Figure8 shows the obtained spectra. From (a, b, c), the spectra corresponded to sample with serious defect, sample with early defect, and sound sample, respectively.Figure 8
Spectra of signals recorded on different samples. (a–c): (a) sample with serious defect, (b) sample with early defect, and (c) sound sample.
(a)
(b)
(c)The results of stress wave signal analysis show that if the size of defect is larger the low frequency part of received stress wave is more, because the high-frequency stress wave will attenuate due to the timber defect. And the results also indicate spectra analysis method can be used to distinguish timbers with large defect, small defect, and defect-free by the frequency distribution of stress wave.


## 4.1. Results of Time Domain Analysis

The peak time interval measurement method and the cross-correlation method mentioned above were carried on with our platform. Table1 lists the flight time obtained from the three samples using different methods. In order to reduce the influence of some random factors, the flight time listed in Table 1 is the average of 10 repetitions.Table 1
Flight time obtained by two methods (μs).
Serious defect
Early defect
Defect-free
Peak method
504
160
48
Correlation method
428
144
42It can be seen from Table1 that the flight time obtained by two methods are very close. They are both sensitive to the degree of the defect. It took about 45 μs for stress wave to travel in the defect-free sample and approximately 3 times longer to travel in the sample with early defect. The time of flight corresponding to the seriously sample is 9 times of those in response to the sound sample. Therefore, both of the two methods have the potentiality of application in nondestructive testing of wood and further studies are needed to demonstrate their performance.It also can be observed that the value obtained employing correlation method is a little less than that employing peak method. Literature [19] reported that current time of flight techniques provide stiffness measurements that are higher than those obtained employing other techniques. Therefore, it can be expected that the correlation method has the potential to improve the accuracy of stiffness testing employing time of flight.


## 4.2. Results of Frequency Domain Analysis

With the Labview code of the spectrum analysis, we transformed the signal recorded on the received sensor into frequency domain. Figure8 shows the obtained spectra. From (a, b, c), the spectra corresponded to sample with serious defect, sample with early defect, and sound sample, respectively.Figure 8
Spectra of signals recorded on different samples. (a–c): (a) sample with serious defect, (b) sample with early defect, and (c) sound sample.
(a)
(b)
(c)The results of stress wave signal analysis show that if the size of defect is larger the low frequency part of received stress wave is more, because the high-frequency stress wave will attenuate due to the timber defect. And the results also indicate spectra analysis method can be used to distinguish timbers with large defect, small defect, and defect-free by the frequency distribution of stress wave.


## 5. Conclusion

In this paper, we designed a wood nondestructive testing platform based on NI virtual instrument which can collect, convert, and store the stress wave data. Based on the testing platform, the performance of three methods, including peak time interval measurement method, cross-correlation method, and spectrum analysis method, was compared with each other for the nondestructive testing of wood. The hardware of the platform and the Labview codes for implementing the three methods were described in detail. The results indicated that the designed platform has the flexibility which enable the tentative ideas to be verified easily. All of the three methods have potential to be used for the nondestructive testing of wood. Timbers with large defect and small defect and that are defect-free can be distinguished.This study proves that this platform is feasible and has a great potential in the development of new methods for detecting wood defects in the radial direction. Further research into defect detection in the longitudinal direction is recommended. Moreover, efforts will be made to investigate new methods to evaluate the mechanical behavior of the wood nondestructively with this platform.



---

*Source: 6279869-2016-08-11.xml*