Next Article in Journal
Effect of Heat Curing Method on the Mechanical Strength of Alkali-Activated Slag Mortar after High-Temperature Exposure
Next Article in Special Issue
A Study on Fatigue State Evaluation of Rail by the Use of Ultrasonic Nonlinearity
Previous Article in Journal
Effect of the Nano-Ca(OH)2 Addition on the Portland Clinker Cooking Efficiency
Previous Article in Special Issue
Enhanced Sensing and Data Processing System for Continuous Profiling of Pavement Deflection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical and Experimental Research on Identifying a Delamination in Ballastless Slab Track

1
Shanghai Institute for Advanced Communication and Data Science, School of Communication and Information Engineering, Shanghai University, Shanghai 200444, China
2
School of Urban Railway Transportation, Shanghai University of Engineering Science, Shanghai 201620, China
*
Author to whom correspondence should be addressed.
Materials 2019, 12(11), 1788; https://doi.org/10.3390/ma12111788
Submission received: 3 April 2019 / Revised: 12 May 2019 / Accepted: 31 May 2019 / Published: 2 June 2019
(This article belongs to the Special Issue Non-destructive Inspection)

Abstract

:
This paper aims to adopt the total focusing method (TFM) and wavenumber method for characterizing a delamination in ballastless slab track. Twelve dry point contact (DPC) transducers located at the upper surface of the slab track compose a linear array. These transducers are employed to actuate shear waves, which are suitable for identifying the delamination. The technique of removing the surface wave has been implemented for only retaining the scattered wave caused by the delamination and the reflected wave from the bottom of bed plate. Numerical and experimental results demonstrate that the delamination and bottom of the bed plate can be identified by the proposed methods. Furthermore, the near-surface pseudomorphism is distinctly restrained after removing the surface wave. Compared to TFM, the wavenumber method has the great advantages of improving computational performance and lateral resolution. However, they have no significant difference in the longitudinal resolution. Furthermore, it has been confirmed that the lateral resolution can be affected by the amount of transducers. This paper can provide valuable suggestions on improving the computational performance and the imaging accuracy when we identify a delamination in ballastless slab track.

1. Introduction

The ballastless track has been widely applied in high-speed railways [1]. Compared to the conventional ballast track, the ballastless track has the advantages of high structural stability, good stiffness uniformity, long service life, and less maintenance work [2,3]. Nevertheless, the complex natural environment, rain erosion, train load, and thermal load may lead to the ballastless track bearing some structural diseases, such as cracks, air voids, and delaminations [4,5,6]. These diseases may not be easily recognized by visual investigation.
Among these diseases, the delamination located between multi-layers is frequently observed in ballastless track. The dynamic behavior of slab track might be affected by the smaller delamination. However, the larger delamination will dramatically reduce the safety performance, and even result in the potential failure of slab track, such as fracture and transfiguration. It is significant to obtain the location and depth of the delamination. Many researchers have proposed various nondestructive testing (NDT) techniques, such as ground penetrating radar (GPR) [7,8], active infrared thermography (IRT) technique [9,10], and ultrasonic testing (UT) method [11].
Hong et al. utilized high frequency GPR system for periodically monitoring reinforcement corrosion in concrete, and visualizing the accumulative corrosion progress [12]. However, the identification of void diseases in a ballastless track using GPR is difficult. This is mainly due to the fact that the echoes caused by the void will be affected by the reflected wave occurring on the rebars surface. Void diseases are especially located at the lower part of rebars.
Tran et al. proposed the square pulse thermography (SPT) to identify the subsurface delamination in concrete [13]. They investigated the effect that steel bars exerted on the detectability of delamination. The results showed that the absolute contrast of a delamination was reduced slightly under the effect of steel bars. In addition, the test result was relevant to the width-to-depth ratio of the delamination. Similarly, Huh et al. applied long pulsed thermography (LPT) to characterize the delamination in a concrete slab [14]. They employed the non-dimensional prefactor k to predict the depth of a delamination and evaluated the effects of steel bars. The active infrared thermography was an effective tool for detecting the delamination with the lower water/cement ratio. It was verified that the delamination with a large size or being located near the surface was more easily detected by using IRT [15]. However, IRT is unsuited for detecting the deep delamination. Moreover, IRT is strongly influenced by environmental conditions.
The UT method is a commonly used NDT technique and has shown great potential for structural health monitoring (SHM). Damage imaging methods based on elastic waves are multifarious, such as synthetic aperture focusing technique (SAFT) [16,17], total focusing method (TFM) [18], ultrasonic phased array [19], reverse time migration (RTM) [20], and ultrasonic tomography [21].
It is well known that the high frequency-thickness product leads to multimodal guided waves coexisting [22]. Under this condition, the mode identification and separation will become difficult. As a result, it is not easy to identify the damage existing in multi-layers concrete structure using guided waves. Li et al. adopted guided waves to evaluate the integrity of CFRP (Carbon Fiber Reinforced Polymer)-reinforced concrete structures [23]. Their work indicated that debonding indices were only applicable for quantifying the severity of debonding. Song et al. demonstrated that interface bonding conditions of steel–clad concrete structures could be quantitatively evaluated using attenuation characteristics of guided waves [24].
Considering the above reasons, many researchers tend to apply bulk waves for detecting damage in concrete structure. Schabowicz found that the ultrasonic tomography technique can not only be used to determine the location and size of a defect, but also to estimate the thickness of a concrete member [25]. One minor flaw was that the components of the received signals had not yet been analyzed. Beniwal et al. investigated the ultrasonic imaging of concrete using scattered elastic wave modes [26]. Compared with the traditional TFM, the weighted sum technique and statistical technique were conducive to improve the detection performance. However, the near-surface pseudomorphism could not be ignored as the surface wave with high amplitude was contained in the received signal.
For the purpose of maintaining the safety and stability of the high-speed railway, the TFM and wavenumber method based on full matrix capture (FMC) were used for identifying a delamination in ballastless slab track. In this paper, 12 DPC transducers located at the upper surface of the slab track were employed to actuate shear waves [27]. The main work was focused on suppressing the near-surface pseudomorphism, improving the computational performance, and retaining a relatively favorable accuracy.
This paper is organized with seven sections, including this Introduction. Section 2 illustrates the theory of TFM and the wavenumber method. Section 3 presents the numerical results. Experimental results are explicated in Section 4. Section 5 concludes the discussion. Section 6 summarizes the main work. Section 7 provides the recommendation for future research.

2. Principles of the TFM and the Wavenumber Method

2.1. Sound Field Transmission and Reception

The transmission and reception of the ultrasonic phased array is introduced based on the two-dimensional (2-D) medium. The schematic diagram used for describing the process of sound field transmission and reception is shown in Figure 1. Linear array elements are located at the upper surface of the material, while a defect is embedded into the material. The X axis is parallel to the material’s surface, and the Z axis is perpendicular to the material’s surface. The elements are spaced at a fixed pitch Δ u , and each element has a width of D. The coordinate origin is at the center of the array.
Generally, the signals for each transmitter–receiver pair within the array are recorded in space and time. The full matrix capture data [28] is shown in Figure 2. S i j denotes the signal which is obtained under transmitter element i and receiver element j. The three-dimensional (3-D) matrix corresponds to the transmitter and receiver positions and the temporal axis. ( u , 0 ) and ( v , 0 ) are the Cartesian coordinates of the transmitter and receiver elements. ( x , z ) is the Cartesian coordinate of the defect. r 1 is the distance between the transmit element and the defect, and r 2 is the distance between the receive element and the defect. Naturally, r 1 and r 2 can be expressed as:
r 1 = ( x u ) 2 + z 2
r 2 = ( x v ) 2 + z 2

2.2. Principle of the TFM

TFM is a type of ultrasonic imaging and it has been widely used for Non-destructive Evaluation (NDE). The first step of implementing TFM algorithm is discretizing the target region into a grid. Then, the travelled time is calculated by the travelled distance and the group velocity [29]. This operation can be denoted as:
T i j ( x , z ) = ( x x i ) 2 + ( z z i ) 2 + ( x x j ) 2 + ( z z j ) 2 C g
where, ( x i , z i ) is the Cartesian coordinates of the transmitting sensor. ( x j , z j ) is the Cartesian coordinates of the receiving sensor ( x , z ) is the Cartesian coordinates of the gird. c g is the group velocity. T i j ( x , z ) denotes the travelled time.
Finally, the amplitude of the signal’s envelope (relative to the transducer pair) associated to this travelled time is extracted and added to the pixel intensity [30]. This operation is repeated for each transmitter-receiver pair.
I ( x , z ) = | i = 1 N j = 1 N S i j ( T i j ( x , z ) ) |
here, I ( x , z ) is the intensity of the pixel located at the Cartesian coordinates ( x , z ) . N denotes the amount of sensors. S i j is the amplitude of the signal’s envelope.

2.3. Principle of the Wavenumber Method

The wavenumber method based on FMC is adapted from the theories described by Hunter et al. and Moghimirad et al. [31,32]. The algorithm is expatiated as follows. First of all, the received echo signal e ( t , u , v ) is given by E ( ω , u , v ) in frequency domain.
E ( ω , u , υ ) = P ( ω ) f ( x , z ) g ( ω , x u , z ) g ( ω , x v , z ) d x d z
where P ( ω ) is the Fourier transform of the transmitted signal P ( t ) . f ( x , z ) denotes a continuous distribution of the scatterer. The 2-D free-space Green’s function can be expressed as [33]:
g ( ω , x , z ) = j 4 π exp ( j | z | k 2 k x 2 + j k x x ) k 2 k x 2 d k x
Naturally, the temporal spectrum of the received echo signal can be rewritten as:
E ( ω , u , v ) = P ( ω ) 1 ( 4 π ) 2 exp ( j k u u + j k v v ) k 2 k u 2 k 2 k v 2 [ f ( x , z ) × exp ( j ( k u + k v ) x j ( k 2 k u 2 + k 2 k v 2 ) z ) × d x d z ] d k u d k v
The next step is to obtain the frequency-wavenumber spectrum E ( ω , k u , k v ) by applying the 2-D fast Fourier transform [34] of E ( ω , u , v ) .
E ( ω , k u , k v ) = P ( ω ) F ( k u + k v , k 2 k u 2 + k 2 k v 2 ) ( 4 π ) 2 k 2 k u 2 k 2 k v 2
here, F is the spatial Fourier transform of the scattering function f ( x , z ) . The process of transforming the E ( ω , u , v ) to E ( ω , k u , k v ) can also be accomplished by using the 3-D fast Fourier transform [35]. The system model of the incident wavenumber k u is given by:
F ( k x , k z | k u ) = ( 4 π ) 2 S 1 { k 2 k u 2 k 2 k v 2 E ( ω , k v | k u ) }
where the operator S 1 { } denotes the Stolt mapping, which can be obtained from the Equations (10)–(12). E ( ω , k v | k u ) denotes the 2-D slices for each value of k u . In this case, the system model in wavenumber domain is obtained by summing together F ( k x , k z | k u ) . The operation given by Equation (13) is for the purpose of reducing the side lobes.
k x = k u + k v
k z = k 2 k u 2 + k 2 k v 2
k = ± k z 4 + 2 ( k u 2 + ( k x k u ) 2 ) k z 2 + k u 4 + ( k x k u ) 4 2 k u 2 ( k x k u ) 2 2 k z
F ^ ( k x , k z ) = F ( k x , k z | k u ) d k u
Finally, the image of the scatterer distribution can be reconstructed by applying the 2-D inverse Fourier transform (IFFT2) and this operation is given by Equation (14). It is important to note that IFFT2 is based on the uniform baseband grid. However, the raw baseband grid is inhomogeneous as the Stolt mapping is a kind of non-linear coordinate transforms [36,37]. Therefore, f ^ ( k x , k z ) is obtained by implementing the complex-valued linear interpolation on F ^ ( k x , k z ) . The wavenumber method based on FMC is summarized and its flow chart is shown in Figure 3.
f ^ ( x , z ) = 1 ( 2 π ) 2 f ^ ( k x , k z ) exp ( j k x x + j k z z ) d k x d k z

3. Numerical Simulation

3.1. Finite Element Model

A commercial software package called PZFlex was used to create the 2-D finite element model. PZFlex is the time domain finite element method (FEM). The complex acoustic wave propagation can be simulated using this method [38].
The ballastless slab track is a typical concrete structure with three layers. It consists of the track slab, the self-compacting concrete and the bed plate from top to bottom. It should not be ignored that the subgrade is under the bed plate. The cross-section of the ballastless slab track in longitudinal direction is shown in Figure 4. Parameters of the finite element model are listed in Table 1. The upper and lower interfaces of the finite element model are with the free boundary conditions, while the left and right interfaces are with the absorbing boundary conditions. The sizes of the slab track, the self-compacting concrete, the bed plate and the subgrade are 400 mm × 200 mm, 400 mm × 60 mm, 400 mm × 170 mm, and 400 mm × 200 mm, respectively. The delamination with the size of 300 mm × 30 mm is located at the boundary between the self-compacting concrete and the bed plate. Twelve transducers are located at the upper surface of the slab track, and the central distance between adjacent transducers is 30 mm. The vibration direction of the transducer is along the horizontal direction, which is suitable to actuate the shear wave. The parameters of the ultrasonic transducer array are listed in Table 2.
Considering that the damping existed in slab track can not be ignored under the higher frequency, the central frequency of the excitation source is set at 50 kHz. Figure 5 shows the excitation signal, which is a three-cycle sinusoidal tone burst modulated by the Hanning window. In our test, the shear wave was about 2,466 m/s. Therefore, the wavelength was 49.320 mm ( λ = 2466/50 = 49.320 mm). It is conducive to the detection of delamination when the value of wavelength remains in a reasonable range. The sampling interval was only 1 mm ( λ / 50 ). Apparently, the sampling interval satisfies the sampling theorem and can provide a good spatial resolution.
The process of the excitation and reception should be in accordance with the FMC. In other words, each transducer has a chance to act as the excitation, and the acoustic signals are received by all the transducers. Naturally, a total of 144 signals will be received as the linear array consists of twelve transducers.

3.2. FEM Results

3.2.1. Snapshots of the Wavefield

The interactions between wavefield and the delamination have been analyzed by many researchers. Lin et al. described a layered structure with an embedded delamination and illustrated the screenshots of the wave propagation [39]. The reflection and refraction seemed complex at the overlay-specimen interface as the wave propagation velocities were obviously different in multi-layers structure. Figure 6 shows the snapshots of the ballastless slab track in longitudinal direction. These snapshots are generated by the FEM program when the sixth transducer acts as the excitation source. It should be noted that the wavefield in the subgrade has not been shown as we focus on the ballastless slab track. The incident, reflected, and transmitted waves can be observed. Figure 6d indicates that the scattered wave from the delamination is clear while the reflected wave from the bottom is very weak.

3.2.2. FEM Results

Figure 7 shows the simulated signal in the time- and frequency-domain. It can be clearly found that the amplitude of surface wave is much greater than that of the scattered wave caused by the delamination and the reflected wave from bottom of the bed plate. In fact, it is only necessary to retain the scattered and reflected waves as they are beneficial to the detection. A normal method is to subtract the signal acquired in the pristine structure from that acquired in the damaged structure. The process of removing surface wave for multiple signals is shown in Figure 8.
The simulation results based on TFM are shown in Figure 9. The rectangle with white dotted line represents the location of delamination. The delamination and the bottom of bed plate both were obscure in Figure 9a as the pseudomorphism caused by the surface wave was conspicuous. It could be found that the near-surface pseudomorphism was distinctly suppressed after removing the surface wave. Thus, the delamination and the bottom of bed plate were clearly seen in Figure 9b. As a whole, the location of identified delamination matches well with its actual location. However, the identified delamination is shorter than its actual size, revealing that the lateral resolution of the TFM must be improved.
The simulation results based on wavenumber method are shown in Figure 10. Benefiting from removing the surface wave, the delamination and the bottom of bed plate were well identified. Furthermore, the identified delamination area matched well with its actual size. Compared to TFM, the wavenumber method took less time in obtaining the intensity imagery while maintaining a relatively accurate location and better lateral resolution.

4. Experimental Verification

4.1. Experimental Setup

The experimental system is composed of DPC transducers, an ultrasonic device (JPR-600C, Japan Probe, Yokohama, Japan), a pre-amplifier (PR-60A-06, Japan Probe, Yokohama, Japan) and the signal acquisition interface. The DPC transducer is a piezoelectric ceramic-based transducer with the central frequency of 50 kHz, and its bandwidth is about ±60% of the central frequency. The pre-amplifier was used to improve the signal-to-noise ratio (SNR) and its gain was set at 60 dB. The sampling frequency was set at 1 MHz. It is a suitable time to record the signal when the waveform tends to be stable. Each measurement was repeated 15 times in order to reduce the measurement error. The signal acquisition interface can display the detailed information of received signals, such as the sample interval, trigger point, amplitude, etc. The experimental conditions were in agreement with the simulation conditions. In fact, the ultrasonic device only has two channels. Nevertheless, the FMC still can be successfully completed. More specifically, one channel was used for transmitting. The other one was used for picking up direct and reflected waves. First of all, the excitation source was located at the fixed position. Then, the receiver was removed with a fixed interval after the signal was recorded. Thus, multi-channel received signals could be successfully obtained by using two transducers. The experimental setup is shown in Figure 11.

4.2. Experimental Results

4.2.1. Removing Surface Wave

Filtering and removing surface wave are required for experimental signals. First of all, the received signals are filtered by a Chebyshev band-pass filter [40] with a bandwidth of 20–80 kHz. Consequently, the noise is successfully removed and the filtered signal is shown in Figure 12.
Generally, it is very hard to obtain the same surface wave in the pristine structure and the damaged structure as the real ballastless slab track is nonhomogeneous. Therefore, the method of removing surface wave is different from the case of FEM. In our test, the surface wave, the scattered wave from the delamination, and the reflected wave from the bottom were segregative in the time-domain. Therefore, the surface wave can be separated in the time-domain. The surface wave velocity of the track slab denoted by C r can be obtained by calculating the following equation [41]:
C r = 0.87 + 1.12 u 1 + u E ρ ( 1 2 ( 1 + u ) )
where, u is the Poisson’s ratio. ρ is the density. E is elasticity modulus. The values of these parameters have been listed in Table 1. The surface wave velocity of the slab track is about 2249 m/s. The travelled time of surface wave T 1 can be written as:
T 1 = ( x j x i ) 2 + ( z j z i ) 2 C r
T 2 = T 1 + T 0
where, T 0 is the duration of the original excitation signal envelope. T 2 is the deadline of the surface wave. The received signal is multiplied by a very small ratio coefficient before T 2 and filtered in the frequency-domain. In this paper, the process of removing surface wave is illustrated under the condition that No.1 transducer is selected as the excitation source. Figure 13 and Figure 14 reveal the processing procedure for a single signal and multiple signals, respectively. It must be pointed out that the surface wave can be removed in frequency-wavenumber domain if the wave packets are overlapped in the time-domain.

4.2.2. Experimental Results

Figure 15 and Figure 16 show the experimental results based on TFM and the wavenumber method, respectively. The experimental results were in good agreement with the simulation results. It can be found that a relatively high accuracy was achieved if the surface wave was removed successfully. Meanwhile, it is reconfirmed that the wavenumber method can achieve a better lateral resolution.

5. Discussion

The TFM and wavenumber method have been used for characterizing the delamination in a ballastless slab track. Two methods are compared in the computational performance, lateral resolution, and longitudinal resolution. In addition, the difference between simulation and experimental results has been pointed out. Moreover, the delamination and the bottom of bed plate are compared in the imaging results. Finally, the influence caused by the amount of transducers has been analyzed.
(1) The TFM and the wavenumber method are compared in the computational performance.
The wavenumber method can provide with the superior computational performance. We consider a linear array with J transducers. Each signal is recorded by K time samples. The output image is assumed with M × N pixels. The computational performances for the TFM and the wavenumber method can be evaluated by [31,32]:
C T F M = O ( M N J 2 )
C ω k = O ( J 2 K log 2 ( J 2 K ) + M N J + M N log 2 ( M N ) )
where, the Landau notation O ( ) is used to denote the order of computational complexity. C T F M and C ω k denote the computational complexities of TFM and wavenumber method, respectively. In this paper, the relevant parameters were set as: M = 490, N = 330, J = 12, and K = 2048. Therefore, the computational complexity ratio of TFM to wavenumber method was C T F M / C ω k = 2.306 theoretically. In practice, the time cost of using TFM is given by T T F M = 4.192 s. While the time cost of using the wave number method is given by T ω k = 2.011 s. The times were measured on a personal computer with Core i7 3.40 GHz CPU and 24 GB RAM. It is easy to obtain the ratio T T F M / T ω k 2 . 084 . The theoretical and experimental results confirm that the wavenumber method can provide with the superior computational performance. Especially, the detection is implemented in a large area or the array consists of many transducers.
(2) The TFM and the wavenumber method are compared in the lateral resolution and the longitudinal resolution.
In this paper, the actual length and thickness of the delamination were 300 mm and 30 mm, respectively. Figure 17 shows the profile curves of experimental results. The size of delamination is estimated by setting a threshold at the level of −6 dB and results are shown in Table 3. The measured delamination length is only 201 mm by using the TFM. However, this value reached up to 234 mm by using the wavenumber method. The measured delamination thickness is 45 mm by using the TFM, and this value was 49 mm by using the wavenumber method.
It is easy to draw a conclusion that the wavenumber method is better than the TFM in the lateral resolution, but there is no significant difference in the longitudinal resolution. One reasonable explanation is that a mathematically rigorous solution to the inverse problem is obtained in wavenumber method, while the TFM only adopts the delay-and-sum beam-forming. Therefore, the wavenumber method can achieve better accuracy. It must be pointed out that the wavenumber method requires a linear array, while the TFM is applicable with a flexible transducer array.
(3) The difference between simulation results and experimental results is analyzed.
Experimental results show that the delamination and the bottom of bed plate can be easily identified without removing the surface wave, while they are obscure in simulation results. This discovery is obtained by comparing Figure 9a with Figure 15a, or comparing Figure 10a with Figure 16a. To explain this phenomenon, the simulated signal and the experimental signal were compared in time-domain. Figure 18a shows that the amplitude of surface wave is much higher than that of the delamination scattered wave and the bottom reflected wave. However, Figure 18c shows that their difference in amplitude is not so serious. Figure 18b,d is the time-frequency spectra for the simulated signal and experimental signal, respectively. The time-frequency spectra have further confirmed that the pseudomorphism caused by the surface wave is serious under the simulation condition.
(4) The location of the delamination is precise while the bottom of bed plate can not be identified very well.
On the one hand, it is not appropriate to assume the propagation speed as a constant during the process of implementing algorithms. In fact, the ballastless slab track is a typical multi-layers structure. We should pay more attention to the fact that the wave propagation velocities in multi-layers are different. More specifically, the shear wave velocity is about 2466 m/s in the slab track. This value is about 2377 m/s in the self-compacting concrete layer and the bed plate. Although the difference in velocity is not obvious, however, it still can affect the test quality at the bottom of bed plate.
On the other hand, the imagery quality is relevant to the attenuation coefficient of materials. In other words, the amplitude of the scattered wave caused by the delamination is higher than that of the reflected wave from the bottom because of the attenuation. Consequently, the delamination is conspicuous while the bottom of bed plate is obscure.
(5) The amount of transducers has a great influence on the test results.
In order to make the result be more intuitive, all the analysis was based on the experimental signals, which have been removed surface waves. The test results by applying different amount of transducers are shown in Figure 19 and Figure 20. It could be found that the identified delamination length was gradually shortening with the decrease of transducers. Meanwhile, the expansion of the residual pseudomorphism was inevitable. Generally, the more the transducers are, the better test results can be achieved.

6. Conclusions

The paper presents the TFM and wavenumber method for the delamination localization in ballastless slab track. Both simulation results and experimental results demonstrate that the delamination and the bottom of the bed plate can be identified by the proposed methods. The conclusions can be drawn as follows:
(1)
The test result has been distinctly improved as the near-surface pseudomorphism is restrained by applying the technique of removing the surface wave.
(2)
Compared to TFM, the wavenumber method can provide with the better computational performance and the lateral resolution. However, they have no significant difference in the longitudinal resolution.
(3)
It can be observed that the delamination is conspicuous while the bottom of bed plate is obscure in the image. The explanation could be the inhomogeneous propagation speed and the attenuation existed in concrete.
(4)
The test result is relevant to the amount of transducers. Generally, the more the transducers are, the better test result can be achieved.
(5)
The low frequency air-coupled ultrasonic inspection combined with Lamb wave is viewed as a preliminary step for obtaining the location of the delamination in the transversal direction.
(6)
Compared with the classical delay-and sum (DAS) method, the wavenumber method can greatly reduce the computational time while retaining a relatively favorable accuracy.
The proposed methods can provide useful information for the early warning of the structural failure. It must be noted that the accurate shear-wave velocity and removing surface wave are very crucial for the delamination detection.

7. Future Work

Detection with a large area will be taken into consideration for future studies. In practice, the delamination location is usually unknown in a ballastless slab track. For this reason, a large amount of data must be processed if the proposed methods are employed for the detection with a large area. Naturally, consuming much time is inevitable. In addition, it is not flexible, as the DPC transducers still have to contact the surface of the detected object. Fortunately, air-coupled ultrasonic inspection combined with Lamb wave is a promising non-contact method for the rapid and convenient inspection [42,43,44]. More specifically, the low frequency air-coupled ultrasonic inspection combined with Lamb wave is viewed as a preliminary step for obtaining the location of the delamination in the transversal direction. On this basis, the proposed methods are implemented on the small scanning region where delamination is likely to exist. Consequently, delamination can be localized both in the transversal and longitudinal directions. The detailed solutions will be introduced in future work.

Author Contributions

G.F. wrote this article and performed the FEM simulation; H.Z. (Haiyan Zhang) was the principle researcher on this project and revised the manuscript; W.Z. and H.Z. (Hui Zhang) designed the experiment; X.C. analyzed the experimental data and provided many valuable suggestions.

Funding

This research was funded by the National Natural Science Foundation of China (Grant Nos. 11874255 and 11674214), and the Key Technology R&D Project of Shanghai Committee of Science and Technology (Grant No. 16030501400).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yu, Z.; Xie, Y.; Shan, Z.; Li, X. Fatigue performance of CRTS III slab ballastless track structure under high-speed train load based on concrete fatigue damage constitutive law. J. Adv. Concr. Technol. 2018, 16, 233–249. [Google Scholar] [CrossRef]
  2. Chen, L.; Chen, J.; Wang, J. Calculation of Reasonable Tension Value for Longitudinal Connecting Reinforcement of CRTS II Slab Ballastless Track. Appl. Sci. 2018, 8, 2139. [Google Scholar] [CrossRef]
  3. Shan, Y.; Zheng, S.; Zhang, X.; Luo, W.; Mao, J.; Kong, D. Fatigue Performance of the CA Mortar Used in CRTS I Ballastless Slab Track under Simulated Servicing Condition. Materials 2018, 11, 2259. [Google Scholar] [CrossRef] [PubMed]
  4. Cao, S.; Zhai, S.; Dai, F.; Deng, S.; Wang, W. Hydrodynamic Pressure and Velocity Distributions in the Interlayer Crack of Ballastless Track under High-Speed Train Load: A Theoretical Analysis. Math. Probl. Eng. 2018, 3187979. [Google Scholar] [CrossRef]
  5. Yang, R.; Li, J.; Kang, W.; Liu, X.; Cao, S. Temperature Characteristics Analysis of the Ballastless Track under Continuous Hot Weather. J. Transp. Eng. Pt A-Syst. 2017, 143, 04017048. [Google Scholar] [CrossRef]
  6. Liu, J.; Wen, X.; Zhang, Z.; Li, S.; Zeng, Z. Influence of the Stabilizer on Interfacial Bonding Behavior of Cement Asphalt Mortar in Slab Ballastless Track. J. Mater. Civ. Eng. 2018, 30, 04018245. [Google Scholar] [CrossRef]
  7. Tasker, L.; Karrech, A.; Shragge, J.; Josh, M. Time-lapse monitoring of internal alteration of a concrete structure using ground penetrating radar. Constr. Build. Mater. 2018, 191, 300–310. [Google Scholar] [CrossRef]
  8. Cassidy, N.J.; Eddies, R.; Dods, S. Void detection beneath reinforced concrete sections: The practical application of ground-penetrating radar and ultrasonic techniques. J. Appl. Geophys. 2011, 74, 263–276. [Google Scholar] [CrossRef]
  9. Omar, T.; Nehdi, M.L.; Zayed, T. Infrared thermography model for automated detection of delamination in RC bridge decks. Constr. Build. Mater. 2018, 168, 313–327. [Google Scholar] [CrossRef]
  10. Hiasa, S.; Birgul, R.; Catbas, F.N. A data processing methodology for infrared thermography images of concrete bridges. Comput. Struct. 2017, 190, 205–218. [Google Scholar] [CrossRef]
  11. Ongpeng, J.; Oreta, A.W.C.; Hirose, S. Characterization of Damage Using Ultrasonic Testing on Different Types of Concrete. Mater. Eval. 2018, 76, 1532–1541. [Google Scholar]
  12. Hong, S.; Wiggenhauser, H.; Helmerich, R.; Dong, B.; Dong, P.; Xing, F. Long-term monitoring of reinforcement corrosion in concrete using ground penetrating radar. Corrosion Sci. 2017, 114, 123–132. [Google Scholar] [CrossRef]
  13. Tran, Q.H.; Huh, J.; Mac, V.H.; Kang, C.; Han, D. Effects of rebars on the detectability of subsurface defects in concrete bridges using square pulse thermography. NDT E Int. 2018, 100, 92–100. [Google Scholar] [CrossRef]
  14. Huh, J.; Mac, V.H.; Tran, Q.H.; Lee, K.Y.; Lee, J.I.; Kang, C. Detectability of Delamination in Concrete Structure Using Active Infrared Thermography in Terms of Signal-to-Noise Ratio. Appl. Sci. 2018, 8, 1986. [Google Scholar] [CrossRef]
  15. Rocha, J.H.A.; Póvoas, Y.V.; Santos, C.F. Detection of Delaminations in Sunlight-Unexposed Concrete Elements of Bridges Using Infrared Thermography. J. Nondestruct. Eval. 2019, 38, 8. [Google Scholar] [CrossRef]
  16. Tong, J.H.; Chiu, C.L. Elastic-wave-based synthetic aperture focusing technique imaging system for detecting voids inside concrete structures. Struct. Control. Health Monit. 2017, 24, e1956. [Google Scholar] [CrossRef]
  17. Ganguli, A.; Rappaport, C.M.; Abramo, D.; Fascetti, S.W. Synthetic aperture imaging for flaw detection in a concrete medium. NDT E Int. 2012, 45, 79–90. [Google Scholar] [CrossRef]
  18. Tseng, P.Y.; Chang, Y.F.; Lin, C.M.; Nien, W.J.; Chang, C.H.; Huang, C.C. A Study of Total Focusing Method for Ultrasonic Nondestructive Testing. J. Test. Eval. 2013, 41, 557–563. [Google Scholar] [CrossRef]
  19. Huggett, D.J.; Dewan, M.W.; Wahab, M.A.; Okeil, A.; Liao, T.W. Phased Array Ultrasonic Testing for Post-Weld and OnLine Detection of Friction Stir Welding Defects. Res. Nondestruct. Eval. 2017, 28, 187–210. [Google Scholar] [CrossRef]
  20. Grohmann, M.; Mueller, S.; Niederleithinger, E.; Sieber, S. Reverse time migration: Introducing a new imaging technique for ultrasonic measurements in civil engineering. Near Surf. Geophys. 2017, 15, 242–258. [Google Scholar]
  21. Choi, H.; Popovics, J.S. NDE Application of Ultrasonic Tomography to a Full-Scale Concrete Structure. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2015, 62, 1076–1085. [Google Scholar] [CrossRef] [PubMed]
  22. Khalili, P.; Cawley, P. Excitation of Single-Mode Lamb Waves at High-Frequency-Thickness Products. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2016, 63, 303–312. [Google Scholar] [CrossRef] [PubMed]
  23. Li, J.; Lu, Y.; Guan, R.; Qu, W. Guided waves for debonding identification in CFRP-reinforced concrete beams. Constr. Build. Mater. 2017, 131, 388–399. [Google Scholar] [CrossRef]
  24. Song, H.; Popovics, J.S. Characterization of steel-concrete interface bonding conditions using attenuation characteristics of guided waves. Cem. Concr. Compos. 2017, 83, 111–124. [Google Scholar] [CrossRef]
  25. Schabowicz, k. Ultrasonic tomography-The latest nondestructive technique for testing concrete members Description, test methodology, application example. Arch. Civ. Mech. Eng. 2014, 14, 295–303. [Google Scholar] [CrossRef]
  26. Beniwal, S.; Ghosh, D.; Ganguli, A. Ultrasonic imaging of concrete using scattered elastic wave modes. NDT E Int. 2016, 82, 26–35. [Google Scholar] [CrossRef]
  27. Hoegh, k.; Khazanovich, L. Extended synthetic aperture focusing technique for ultrasonic imaging of concrete. NDT E Int. 2015, 74, 33–42. [Google Scholar] [CrossRef] [Green Version]
  28. Caulder, A. Full Matrix Capture and Total Focusing Method: The Next Evolution in Ultrasonic Testing. Mater. Eval. 2018, 76, 591–597. [Google Scholar]
  29. Bazulin, E.G. Reconstruction of reflector images using the C-SAFT method with account for the anisotropy of the material of the test object. Russ. J. Nondestr. Test. 2015, 51, 217–226. [Google Scholar] [CrossRef]
  30. Lin, L.; Cao, H.; Luo, Z. Total focusing method imaging of multidirectional CFRP laminate with model-based time delay correction. NDT E Int. 2018, 97, 51–58. [Google Scholar] [CrossRef]
  31. Hunter, A.J.; Drinkwater, B.W.; Wilcox, P.D. The Wavenumber Algorithm for Full-Matrix Imaging Using an Ultrasonic Array. IEEE Trans. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2008, 55, 2450–2462. [Google Scholar] [CrossRef] [PubMed]
  32. Moghimirad, E.; Villagómez Hoyos, C.A.; Mahloojifar, A.; Asl, B.M.; Jensen, J.A. Synthetic Aperture Ultrasound Fourier Beam formation Using Virtual Sources. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2016, 63, 2018–2030. [Google Scholar] [CrossRef] [PubMed]
  33. Chew, W.C. Waves and Fields in Inhomogeneous Media; IEEE Press: Piscataway, NJ, USA, 1994. [Google Scholar]
  34. Fateri, S.; Boulgouris, N.V.; Wilkinson, A. A two-dimensional fast Fourier transform using incremented frequency measurement for guided wave analysis. Insight 2014, 56, 499–504. [Google Scholar] [CrossRef]
  35. Lee, C.; Park, S. Flaw Imaging Technique for Plate-Like Structures Using Scanning Laser Source Actuation. Shock Vib. 2014, 9, 1–14. [Google Scholar] [CrossRef]
  36. Li, Z.; Wang, J.; Liu, Q.H. Interpolation-Free Stolt Mapping for SAR Imaging. IEEE Geosci. Remote Sens. Lett. 2014, 5, 926–929. [Google Scholar] [CrossRef]
  37. Moghimirad, E.; Mahloojifar, A.; Mohammadzadeh, A.B. Computational Complexity Reduction of Synthetic Aperture Focus in Ultrasound Imaging Using Frequency-Domain Reconstruction. Ultrasonic Imaging 2016, 38, 175–193. [Google Scholar] [CrossRef] [PubMed]
  38. Fan, G.; Zhang, H.; Zhang, H.; Zhu, W.; Chai, X. Lamb Wave Local Wavenumber Approach for Characterizing Flat Bottom Defects in an Isotropic Thin Plate. Appl. Sci. 2018, 8, 1600. [Google Scholar] [CrossRef]
  39. Lin, S.; Shams, S.; Choi, H.; Azari, H. Ultrasonic imaging of multi-layer concrete structures. NDT E Int. 2018, 98, 101–109. [Google Scholar] [CrossRef]
  40. Losada, R.A.; Pellisier, V. Designing IIR Filters with a Given 3-dB Point. IEEE Signal Process. Mag. 2005, 22, 95–98. [Google Scholar] [CrossRef]
  41. Petro, J.T.; Kim, J. Detection of delamination in concrete using ultrasonic pulse velocity test. Constr. Build. Mater. 2012, 26, 574–582. [Google Scholar] [CrossRef]
  42. Sunarsa, T.Y.; Aryan, P.; Jeon, I.; Park, B.; Liu, P.; Sohn, H. A Reference-Free and Non-Contact Method for Detecting and Imaging Damage in Adhesive-Bonded Structures Using Air-Coupled Ultrasonic Transducers. Materials 2017, 10, 1402. [Google Scholar] [CrossRef] [PubMed]
  43. Chen, J.; Yang, C.; Guo, Q. Evaluation of surface cracks of bending concrete using a fully non-contact air-coupled nonlinear ultrasonic Technique. Mater. Struct. 2018, 51, 104. [Google Scholar] [CrossRef]
  44. Asadollahi, A.; Khazanovich, L. Analytical reverse time migration: An innovation in imaging of infrastructures using ultrasonic shear waves. Ultrasonics 2018, 88, 185–192. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Transmission and reception of the ultrasonic phased array.
Figure 1. Transmission and reception of the ultrasonic phased array.
Materials 12 01788 g001
Figure 2. Schematic diagram of full matrix capture.
Figure 2. Schematic diagram of full matrix capture.
Materials 12 01788 g002
Figure 3. Flow chart of the wavenumber method based on the full matrix capture.
Figure 3. Flow chart of the wavenumber method based on the full matrix capture.
Materials 12 01788 g003
Figure 4. Finite element model.
Figure 4. Finite element model.
Materials 12 01788 g004
Figure 5. The excitation signal in the (a) time-domain and (b) frequency-domain.
Figure 5. The excitation signal in the (a) time-domain and (b) frequency-domain.
Materials 12 01788 g005
Figure 6. Snapshots of wavefield. (a) 66 μs; (b) 107 μs; (c) 127 μs; and (d) 152 μs. The rectangular area delimited by the white dotted lines denotes the delamination. The rectangular area delimited by the black dotted lines represents the self-compacting concrete.
Figure 6. Snapshots of wavefield. (a) 66 μs; (b) 107 μs; (c) 127 μs; and (d) 152 μs. The rectangular area delimited by the white dotted lines denotes the delamination. The rectangular area delimited by the black dotted lines represents the self-compacting concrete.
Materials 12 01788 g006
Figure 7. The simulated signal in the (a) time-domain and (b) frequency-domain.
Figure 7. The simulated signal in the (a) time-domain and (b) frequency-domain.
Materials 12 01788 g007
Figure 8. Removing surface wave for multiple signals. (a) Simulated signals with surface waves and (b) simulated signals without surface waves.
Figure 8. Removing surface wave for multiple signals. (a) Simulated signals with surface waves and (b) simulated signals without surface waves.
Materials 12 01788 g008
Figure 9. Simulation results based on the TFM (total focusing method). (a) Result with surface wave and (b) result without surface wave. The rectangle with white dotted line represents the location of a delamination.
Figure 9. Simulation results based on the TFM (total focusing method). (a) Result with surface wave and (b) result without surface wave. The rectangle with white dotted line represents the location of a delamination.
Materials 12 01788 g009
Figure 10. Simulation results based on the wavenumber method. (a) Result with surface wave and (b) result without surface wave. The rectangle with white dotted line represents the location of a delamination.
Figure 10. Simulation results based on the wavenumber method. (a) Result with surface wave and (b) result without surface wave. The rectangle with white dotted line represents the location of a delamination.
Materials 12 01788 g010
Figure 11. Experimental setup. (a) Cross-section of the delamination in transversal direction and (b) cross-section of the delamination in longitudinal direction.
Figure 11. Experimental setup. (a) Cross-section of the delamination in transversal direction and (b) cross-section of the delamination in longitudinal direction.
Materials 12 01788 g011
Figure 12. Band-pass filtering of an experimental signal. (a) Original signal; (b) spectrum of the original signal; (c) filtered signal; and (d) spectrum of the filtered signal.
Figure 12. Band-pass filtering of an experimental signal. (a) Original signal; (b) spectrum of the original signal; (c) filtered signal; and (d) spectrum of the filtered signal.
Materials 12 01788 g012
Figure 13. Removing surface wave for a single signal. (a) Experimental signal with surface wave and (b) experimental signal without surface wave.
Figure 13. Removing surface wave for a single signal. (a) Experimental signal with surface wave and (b) experimental signal without surface wave.
Materials 12 01788 g013
Figure 14. Removing surface wave for multiple signals. (a) Experimental signals with surface waves and (b) experimental signals without surface waves.
Figure 14. Removing surface wave for multiple signals. (a) Experimental signals with surface waves and (b) experimental signals without surface waves.
Materials 12 01788 g014
Figure 15. Experimental results based on the TFM. (a) Result with surface wave and (b) result without surface wave. The rectangle with white dotted line represents the location of a delamination.
Figure 15. Experimental results based on the TFM. (a) Result with surface wave and (b) result without surface wave. The rectangle with white dotted line represents the location of a delamination.
Materials 12 01788 g015
Figure 16. Experimental results based on the wavenumber method. (a) Result with surface wave and (b) result without surface wave. The rectangle with white dotted line represents the location of a delamination.
Figure 16. Experimental results based on the wavenumber method. (a) Result with surface wave and (b) result without surface wave. The rectangle with white dotted line represents the location of a delamination.
Materials 12 01788 g016
Figure 17. Profile curves of experimental results. (a) Lateral profile curve and (b) longitudinal profile curve.
Figure 17. Profile curves of experimental results. (a) Lateral profile curve and (b) longitudinal profile curve.
Materials 12 01788 g017
Figure 18. Comparison between the simulated signal and the experimental signal. (a) Simulated signal in time domain; (b) simulated signal in time-frequency domain; (c) experimental signal in time domain; and (d) experimental signal in time-frequency domain.
Figure 18. Comparison between the simulated signal and the experimental signal. (a) Simulated signal in time domain; (b) simulated signal in time-frequency domain; (c) experimental signal in time domain; and (d) experimental signal in time-frequency domain.
Materials 12 01788 g018
Figure 19. Experimental results based on the TFM. (a) Twelve transducers; (b) ten transducers; (c) eight transducers; and (d) six transducers. The rectangle with white dotted line represents the location of a delamination.
Figure 19. Experimental results based on the TFM. (a) Twelve transducers; (b) ten transducers; (c) eight transducers; and (d) six transducers. The rectangle with white dotted line represents the location of a delamination.
Materials 12 01788 g019
Figure 20. Experimental results based on the wavenumber method. (a) Twelve transducers; (b) ten transducers; (c) eight transducers; and (d) six transducers. The rectangle with white dotted line represents the location of a delamination.
Figure 20. Experimental results based on the wavenumber method. (a) Twelve transducers; (b) ten transducers; (c) eight transducers; and (d) six transducers. The rectangle with white dotted line represents the location of a delamination.
Materials 12 01788 g020
Table 1. Parameters of the finite element model.
Table 1. Parameters of the finite element model.
PropertySlab TrackSelf-Compacting Concrete LayerBed PlateSubgrade
Concrete strengthC60C40C40/
Density (kg/m3)2500250025001850
Elasticity modulus (Gpa)36.532.532.525
Poisson’s ratio0.20.20.20.35
Thickness (mm)20060170200
Velocity of shear wave2466237723772237
Table 2. Parameters of the ultrasonic transducer array.
Table 2. Parameters of the ultrasonic transducer array.
PropertyValue
Amount of transducers12
Central distance between adjacent transducers (mm)30
Central frequency of the excitation signal (kHz)50
Bandwidth by the level −6 dB (kHz)20–80
Table 3. The size of delamination evaluated by TFM and wavenumber method.
Table 3. The size of delamination evaluated by TFM and wavenumber method.
MethodDelamination Length (mm) Delamination Thickness (mm)
TFM20145
Wavenumber method23449

Share and Cite

MDPI and ACS Style

Fan, G.; Zhang, H.; Zhu, W.; Zhang, H.; Chai, X. Numerical and Experimental Research on Identifying a Delamination in Ballastless Slab Track. Materials 2019, 12, 1788. https://doi.org/10.3390/ma12111788

AMA Style

Fan G, Zhang H, Zhu W, Zhang H, Chai X. Numerical and Experimental Research on Identifying a Delamination in Ballastless Slab Track. Materials. 2019; 12(11):1788. https://doi.org/10.3390/ma12111788

Chicago/Turabian Style

Fan, Guopeng, Haiyan Zhang, Wenfa Zhu, Hui Zhang, and Xiaodong Chai. 2019. "Numerical and Experimental Research on Identifying a Delamination in Ballastless Slab Track" Materials 12, no. 11: 1788. https://doi.org/10.3390/ma12111788

APA Style

Fan, G., Zhang, H., Zhu, W., Zhang, H., & Chai, X. (2019). Numerical and Experimental Research on Identifying a Delamination in Ballastless Slab Track. Materials, 12(11), 1788. https://doi.org/10.3390/ma12111788

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

Article Metrics

Back to TopTop