Next Article in Journal
Industrial Soft Sensor Optimized by Improved PSO: A Deep Representation-Learning Approach
Previous Article in Journal
Using Multi-Dimensional Dynamic Time Warping to Identify Time-Varying Lead-Lag Relationships
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Crack Detection of Threaded Steel Rods Based on Ultrasonic Guided Waves

1
College of Civil Engineering and Architecture, Zhejiang University, Hangzhou 310058, China
2
Center for Balance Architecture, Zhejiang University, Hangzhou 310058, China
3
College of Computer Science and Technology, Zhejiang University, Hangzhou 310027, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(18), 6885; https://doi.org/10.3390/s22186885
Submission received: 14 August 2022 / Revised: 4 September 2022 / Accepted: 7 September 2022 / Published: 12 September 2022
(This article belongs to the Special Issue Sensing Technologies for Damage Detection)

Abstract

:
Fatigue cracks are typical damage of threaded steel rods under dynamic loads. This paper presents a study on ultrasonic guided waves-based, fatigue-crack detection of threaded rods. A threaded rod with given sizes is theoretically simplified as a cylindrical rod. The propagation characteristics of ultrasonic guided waves in the cylindrical rod are investigated by semi-analytical finite element method and the longitudinal L(0, 1) modal ultrasonic guided waves in low frequency band is proposed for damage detection of the rod. Numerical simulation on the propagation of the proposed ultrasonic guided waves in the threaded rod without damage shows that the thread causes echoes of the ultrasonic guided waves. A numerical study on the propagation of the proposed ultrasonic guided waves in the threaded rod with a crack on the intersection of the smooth segment and the threaded segment shows that both linear indexes ( R f and ARS) and nonlinear indexes ( β r e and β ) are able to detect the crack. A constant-amplitude tensile fatigue experiment was conducted on a specimen of the threaded rod to generate fatigue cracks in the specimen. After every 20,000 loading cycles, the specimen was tested by the proposed ultrasonic guided waves and evaluated by the linear indexes and nonlinear indexes. Experimental results show that both the linear and nonlinear indexes of the ultrasonic guided waves are able to identify the crack before it enters the rapid growth stage and the nonlinear indexes detect the crack easier than the linear indexes.

1. Introduction

Threaded steel rods are widely used in engineering structures as load-bearing members or connections such as bridge suspenders, tendons and bolts. These members are usually under high stress and vulnerable to defects and fatigue cracks under periodic dynamic loads [1]. Failure of them will result in safety accidents and economic losses. Nondestructive damage detection of in-service threaded steel rods is a useful way to avoid sudden failure of them.
There are various nondestructive damage detection approaches, such as magnetic particle testing [2,3], eddy current sensing [4,5], acoustic emission [6,7], vibration-based inspection [8,9,10], and ultrasonic testing [11,12,13,14] for steel members. The magnetic-particle method relies on the interaction between magnetic powder and defects. Although its efficiency is improved with modern intelligent magnetic-particle machines [2], the detection process is still cumbersome due to unavoidable procedures such as surface cleaning, magnetization and demagnetization [3]. In addition, it has different sensitivity to the defection in a different magnetization direction, which makes it is difficult to find circumferential cracks paralleling the threads [14]. The eddy current method, which is based on the change of voltage or current in sensing coil, is efficient and sensitive to circumferential cracks. However, the defect at the root of the threads usually causes special signals that are difficult to process [4] and the signal lags severely as the detection depth increases [5]. The acoustic emission method, which judges and locates damage by capturing the high-frequency elastic waves generated by the release of energy when damage occurs, has high spatial resolution and sensitivity [6]. It can be used to detect metal cracks even if the acoustic emission is relatively weak [7]. However, there are some weaknesses in this method that need our attention such as complex cabling and high sensitivity to environmental noise [6]. The vibration-based inspection method, which detects damage through changes of dynamic properties, has been used to detect and localize multiple cracks in uniform rods [8] and non-uniform rods [9]. However, it is vulnerable to the noise, and external factors such as inconsistent boundary conditions and environment changes will affect the stability and accuracy of the detection [10].The ultrasonic method is well-developed and has high-efficiency [13]. It mainly uses straight probes and small-angle oblique, identifying damage from the waveform changes [11,12]. However, the ultrasonic signal will be interfered with by boundary and variable cross-section contour boundary, resulting in disturbance wave, deformation wave, etc., which lowers the effectiveness of damage detection [14]. Moreover, Single probe angle, detection space limitation and poor sonic accessibility also bother the users [14]. Ultrasonic guided waves (UGWs) are a kind of ultrasonic waves that propagate over the entire cross section of the components. Compared with the traditional high-frequency ultrasonic testing, it has lower frequency and is more sensitive to frequency. Its low attenuation and high propagation distance make it possible for long-range and hidden-area detection [15,16].
The UGW method has already been applied for damage inspection of pipes [17], composite plates [18], cables [19], and other types of structural components [20]. In particular, applying the UGW method to damage detection of rod members also has been intensively studied. Shoji et al. experimentally studied the UGW-based detection of anchor rods embedded in soil, and successfully identified the artificial V-shaped necking defect [21,22]. Stepinski et al. proved that time-frequency analysis of UGWs could be used to assess the integrity of rock bolts [23]. Rong et al. analyzed the pulse echo of UGWs in rock bolts, and found that the wave mode L(0, 1) is more sensitive to the area of damage than other wave modes [24]. Zhao et al. found that the guided wave echoes of the corrosion pit in the low-frequency longitudinal wave mode are well correlated with the position and cross-sectional area of the pits [25]. Amjad et al. studied the time-of-flight (TOF) of low-frequency longitudinal mode UGWs in steel rods with various degrees of corrosion, and found that the corrosion level can be quantified by the change in TOF [26].
The influence of threads on the propagation characteristics of UGWs is theoretically insignificant, since the size of threads is much smaller than the wavelength of low-frequency UGWs [27,28]. Therefore, in the numerical simulation of UGWs propagation in threaded rods, especially for externally threaded rods (such as anchor rods [23,27] and rebars [26]), they are mostly simplified as ideal rods. However, for the real threaded rods, especially for the internally threaded rods, the threads cause a change in the cross-section and thus result in significant changes in echoes [29]. To the authors’ knowledge, few studies on UGW-based detection of internally threaded rods have been reported. Meanwhile, most of the previous studies on UGW-based damage detection of threaded rods focused on inspection of artificial damages by analyzing the linear characteristic variations of guided waves [21,22,23,24,25,26,27]. However, the first thread of the internally threaded rod is prone to fatigue due to stress concentration [30], while the linear characteristics of guided waves (e.g., wave speed, echo amplitude, etc.) are not sensitive to micro-fatigue cracks [15,31]. The detection of the fatigue cracks in the threaded rod using UGW-based method is still a technological challenge and worthy of being investigated. This paper carries out a study towards this challenge.
The layout of the paper is as follows: Section 2 presents the propagation characteristics and nonlinear characteristics of guided waves in rods. In Section 3, numerical simulations are carried out to explore the feasibility of using guided waves to detect cracks in threaded rods. Section 4 gives an experimental validation of the UGW-based detection of fatigue cracks in threaded rods. The conclusions of the study are drawn in Section 5.

2. UGWs in Cylindrical Rod

2.1. Propagation Characteristics

The threaded rod can be simplified as a cylindrical rod considering that the size of thread is much smaller than the wavelength of the UGWs, and thus the UGWs in the threaded rod can be theoretically idealized as UGWs in the cylindrical rod [14,27,28]. Figure 1 presents an infinitely long cylinder with a radius of r, which is a homogeneous isotropic elastic medium, and its central axis coincides with the z-axis. There are three groups of modes for the UGWs propagating in the cylinder, i.e., longitudinal mode, torsional mode and flexural mode, which are labeled as L(0, m), T(0, m) and F(n, m), respectively, where n is the harmonic order and m is a sequence number. They are distinguished from each other by the vibration patterns of the particles.
The frequency equation of the UGWs in the cylindrical rod can be obtained by substituting different stress-displacement boundary conditions into the differential equation of motion and the wave equation, and then the modes of UGWs in the cylindrical rod can be determined by solving the frequency equation [32]. Substituting the sectional size and material properties, including diameter d, density ρ , elastic modulus E , Poisson’s ratio v , etc., into the frequency equation, the wavenumber k at different frequencies f can be obtained by semi-analytical finite element (SAFE) method [33]. Then, the dispersion curves of the group velocity ( c g ) of UGWs in the cylindrical rod can be determined by the relationship that c g = d ( 2 π f ) d k . Taking the steel cylinder with a diameter of 27 mm as an example, the dispersion curves of it are obtained by substituting diameter d = 27 mm, density ρ = 7850 kg/m3, elastic modulus E = 206 GPa, and Poisson’s ratio v = 0.3 into the above procedures, as shown in Figure 2. It is clearly shown that the group velocity of each mode varies with frequency, i.e., there is frequency dispersion.
The L(0, 1) modal UGWs at low frequencies (0–150 kHz) are frequently used in the previous studies and the linear characteristic variation of UGWs, including the amplitude and the wave velocity, are used to inspect the damage of rod members [21,22,23,24,25,26,27]. It also can be observed in Figure 2 that, at a given frequency of 60 kHz, the L(0, 1) modal UGWs has a larger group velocity than the T(0, 1) modal UGWs and the F(1, 1) modal UGWs. This observation means that the L(0, 1) modal signal reaches the damaged position of the rod significantly earlier than the signals of other two modes, which is conducive to mode identification and long-distance detection. The cut-off frequency of L(0, 2) modal UGWs is 134 kHz.
The displacement components of L(0, 1) modal UGWs at 60 kHz are plotted in Figure 3, which shows that there is no circumferential displacement uθ and the absolute value of the axial displacement uz is larger than the radial displacement ur. The displacement distribution of L(0, 1) modal UGWs in the cylindrical rod is simple, and its energy is concentrated in the axial direction. These features indicate that L(0, 1) modal UGWs are sensitive to circumferential defects and large radial (depth) defects, but not to axial defects. Besides, in cylindrical rods, the axisymmetric L modal UGWs is easier to be excited than the T modal UGWs and the asymmetric F modal UGWs [15,16]. Therefore, the L(0, 1) modal UGWs at low frequency (≤134 kHz) were chosen for damage detection of threaded rods.
The wavelength of the L(0, 1) modal UGWs at low frequencies is too large to detect tiny cracks effectively [34]. The linear characteristics of the low-frequency L(0, 1) modal UGWs are not sensitive enough for early warning of fatigue cracks [15,31]. Therefore, the nonlinear characteristics of the low-frequency L(0, 1) modal UGWs are also investigated and used for fatigue crack detection of threaded rods in this study.

2.2. Nonlinear Characteristics

2.2.1. Nonlinear Parameters

When UGWs travel in a medium with fatigue cracks, higher order harmonic waves are generated due to Contact Acoustic Nonlinearity (CAN) [35]. For micro-cracks with rough surface, energy dissipation due to micro-contact and friction may also cause higher order harmonic waves [36].
The threaded rod is stressed by the contact forces between the threads and the connected components. The first thread, which is the transition between the smooth segment and the threaded segment, is prone to fatigue due to stress concentration. Local nonlinear damage occurs with the initiation and propagation of micro-cracks. In this case, the “CAN” can be simplified as the elastic modulus change at micro-crack z 1 [37], where the stress σ 1 is given by
σ 1 ( z 1 ) = E 1 ( z 1 ) ε 1 = ( 1 β m a x ( u ( z 1 ) z ) ) E 0 u ( z 1 ) z + β E 0 ( u ( z 1 ) z ) 2
where E 0 is the initial elastic modulus; E 1 is the average value of the elastic modulus after the crack z 1 occurs; u 1 is the displacement at the micro-crack z 1 ; ε 1 is the strain at the micro-crack z 1 ; and β is the nonlinear parameter. The second term on the right side of Equation (1) is the second harmonic waves (SHW). In order to ensure the generation of SHW, the UGWs should oscillate and generate strain at the micro-crack. In other words, the inequation u ( z 1 ) z 0 must be satisfied, which is also known as main mode vibration at the contact surfaces of micro-cracks.
Substituting Equation (1), u 1 = A 1 L sin ( k z 1 ω t ) , max ( u ( z 1 ) z ) = A 1 L k cos ( k z 1 ω t ) and the solution of SHW u 1 H = f ( t ) sin ( 2 k z 1 2 ω t ) + g ( t ) cos ( 2 k z 1 2 ω t ) into the wave equation and ignoring the third and higher order components, SHW can be further simplified as
u 1 H = β c k 2 ( A 1 L ) 2 t 8 cos ( 2 k z 1 2 ω t )
where ω = 2 π f a , where f a is the central frequency of excitation signal; t is the time; and β = A 1 H ( A 1 L ) 2 8 k 2 z is usually referred as the acoustic nonlinear parameter, which can also be derived from the nonlinear Hooke’s law combined with the 1-D wave equation, representing the material nonlinearity [38].
In practice, the relative acoustic nonlinear parameter β and the relative amplitude of the second harmonic β r e given in Equations (3) and (4) are usually used to scale the acoustic nonlinearity.
β = A 1 H ( A 1 L ) 2
β r e = A 1 H A 1 L
where   A 1 L and A 1 H are the amplitude of the fundamental guided waves and SHW in the frequency domain, respectively. The relative acoustic nonlinearity parameter β is adopted to quantify the nonlinearity induced by local micro-cracks. The relative amplitude of the second harmonic β r e is derived based on a bi-linear stiffness model [39], which considers the change of elastic modulus under different pressure and tension. Compared to β r e , β amplifies the amplitude of variation at the fundamental frequency.

2.2.2. Central Frequency of Excitation Signal

Considering that there may be only one end face of the threaded rod suitable for the coupling of the probe, the pulse-echo method is adopted to excite and receive low-frequency UGW signals. In this case, the L(0, 1) mode is the main mode of UGWs propagating in the threaded rod and the displacement of UGWs is mainly distributed in the axial direction. The fatigue crack at the root of the thread mainly grows in both radial and circumferential directions. In other words, the contact surface of micro-crack is perpendicular to the axial direction, which means the L-modal UGWs will cause opening and closing vibration of micro-crack.
Noting that the cut-off frequency of the L(0, 2) mode is 134 kHz, when the second harmonic frequency ( 2 f a ) is greater than 134 kHz, the L(0, 2) mode may be generated. In order to reduce the influence of other modes, the maximum excitation frequency is set as 65 kHz. Further, considering that too small excitation frequency will reduce the sensitivity of the signal to micro-cracks and will cause a long detection blind zone near the excitation end, the range for the central frequency of excitation is finally set at 30–65 kHz.

3. Numerical Study

3.1. UGWs Propagation in Undamaged Threaded Rod

3.1.1. Finite Element Model

A 700 mm-long threaded steel rod with a diameter of 27 mm is considered in this section. There is a 200 mm threaded segment between point B and C, as shown in Figure 4. The pitch of the threads is 3 mm and the depth of the threads is 1.62 mm. The material properties of the steel rod are the same as those given in Section 2.
The commercial finite element software Abaqus [40] was used to carry out the simulation. The mesh size Δ x is suggested to be not larger than λ 20 , where λ = 2 π k is the wavelength of the guided wave, to capture the fluctuation effect between the elements and then to ensure the accuracy and the convergence of the numerical computations [41,42]. Since larger frequency corresponds to smaller wavelength and central frequencies not larger than 70 kHz will be used, the wavelength λ = 70.1 mm corresponding to the central frequency f a = 70 kHz is used to determine the upper bound λ/20 = 3.5 mm. So, the global element size of the numerical model was set as Δ x = 2   mm and the model was meshed with eight-node brick elements (C3D8R) through adaptive meshing.
To simulate the pulse-echo excitation, an axial displacement excitation was applied on a circle area with a diameter of 18 mm at end A of the threaded rod (Figure 4), and a node 2 mm away from the end was selected as the receiver. In this way, there would only be axial and radial displacements in the rod, but no circumferential displacements. In other words, only the L(0, 1) modal guided waves would be excited.
The sinusoidal waves with Hanning window at the central frequencies f a of 30 kHz, 40 kHz, 50 kHz, 60 kHz and 70 kHz were input as the excitation, respectively. Here the maximum excitation frequency is set as 70 kHz to check the influence of L(0, 2) mode to the nonlinear effect of guided waves. The excitation signal is expressed as follows:
f ( t ) = { A ( ( 1 cos ( 2 π f a t n ) ) sin ( 2 π f a t ) 0 0 < t < τ t > τ
where τ = n f a is the excitation time; n is the number of modulation cycle; and A is the amplitude of the signal. Here n = 5 and A = 0.5 were used. With the increasing of n, the energy leakage of the spectrum decrease, which means that more energy is concentrated to central frequency f a . Meanwhile, the time-history signal is more complex due to the increase in the excitation time τ. Figure 5 shows a typical normalized excitation signal with the central frequency f a = 6 0   kHz , whose energy is concentrated to the central frequency f a , which is beneficial to reduce the dispersion phenomenon of the guided wave. Each signal propagated in the rod for 0.0032 s and this process was tracked by the Explicit Solver embedded in Abaqus.
The guided wave signals used in this paper are all normalized as follows:
a ¯ ( t ) = a ( t ) μ σ μ = 1 N i = 1 N a ( t i ) σ 2 = 1 N 1 i = 1 N ( a ( t i ) μ ) 2
where a ( t ) is the guided wave time-history signal; μ is the mean of a ( t ) ; σ is the mean square error; and   N is the number of sampling points. After the normalization, the mean amplitude of the signals becomes zero and the mean square error becomes one.

3.1.2. Signal Analysis

As shown in Figure 6, thread-induced, cross-sectional changes cause some signals to be reflected and these reflected signals are called thread echo. With the increasing of the propagation time, the thread echo is enhanced by multiple reflections at threaded cross-sections. As a result, the waveform becomes more complex and the end echo gradually decreases due to energy loss.
The group velocity c g , the thread reflection coefficient R f , and the location of reflection cross-section L o c f as defined in Equations (7)–(9) can be determined by picking up the excitation signal, the thread echo, and the peak amplitude of the end echo and their propagation times.
c g = 2 L t 1 t 0
R f = H f H 0
L oc f = c g 2 ( t f t 0 )
where H 0 and H f are the peak amplitudes of the excitation signal and the thread echo, respectively; and t 0 , t f and t 1 are the propagation times corresponding to the peak of the excitation signal, the thread echo and the end echo, respectively. Here the time history signal of the first round-trip is selected to calculate the three parameters defined in Equations (7)–(9). The obtained results corresponding to different excitation frequencies are plotted in Figure 7, Figure 8 and Figure 9.
It can be seen from Figure 7 that, in the given frequencies, the group velocity obtained by the numerical simulation of the threaded rod is close to the theoretical result of the simplified smooth rod. To some degree, it justifies the simplification that a threaded rod can be treated as a smooth rod with the same diameter in the analytically solving of the dispersion curve.
It is found that the location of the reflection cross-section is close to the location of the first thread (Figure 8). It indicates that the interface between the smooth segment and the threaded segment is the main cause of the thread echo.
Figure 9 shows that the thread reflection coefficient reaches maximum at 60 kHz. It is also found that the location of reflection cross-section is closest to the location of the first thread at 60 kHz (Figure 8). These observations indicate that the UGWs at 60 kHz is most sensitive to the cross-sectional change at the first thread.
To further investigate the effect of the threads to the propagation of the UGWs, the frequency versus wavenumber relationship obtained by applying 2D Fourier transform on the axial displacement signals of a set of nodes along the length of the rod [43] is plotted in Figure 10. It shows that the wavenumber of the signal propagating in the threaded rod coincides with the wavenumber of L(0, 1) modal UGWs, which indicates that the L(0, 1) mode is the main mode of the signal traveling in the rod. It also shows that there are slight non-axisymmetric bending modes F(1, 1) and F(1, 2) because of the non-axisymmetric thread segment.

3.2. UGWs Propagation in Threaded Rod with Crack

3.2.1. Modeling of “Breathing” Crack

To mimic the fatigue crack, the “surface-to-surface contact” interaction in Abaqus was applied to each surface of the crack, and the normal behavior and the tangential friction coefficient of the interaction were set as “hard” contact and 0.1, respectively.
There are 6 cases of “CAN” area considered to mimic fatigue cracks with different radial depths. They are named as C1-C5 and CB6, and the corresponding “CAN” areas are given in Table 1. Among them, the cracks in cases C1-C5 are fully “breathing” and the whole cracked area has “CAN”. The crack in case CB6 is deemed as open wider than others and there is no “breathing” effect at the partial area close to the surface. Note, that the cases C5 and CB6 have the same cracked area. Illustrations for the numerical modeling of fully “breathing” crack and partially “breathing” crack are shown in Figure 11.

3.2.2. Signal Analysis

The time domain envelopes of full signals are obtained by normalizing and Hilbert transforming the guided wave signals in different cases successively. The Hilbert transform was performed to the difference of normalized guided wave signals under undamaged state and damaged state to obtain time domain envelopes of residual signals.
As shown in Figure 12, the amplitude of the crack echo increases and the amplitude of the end echo decreases with the increase in the crack area. It is notable that there are two obvious wave packets before the first end echo. The first wave packet is the L(0, 1) modal crack echo, and the second one is the F(1, 1) modal echo, which is generated by the guided wave conversion. Although with the same crack area as case C5, the amplitude of the crack echo of CB6 is larger than that of C5 and the amplitude of the end echo of CB6 is smaller than that of C5. This finding indicates that the energy reflected by the full “breathing” crack is smaller than that reflected by the partial “breathing” crack and more energy transmits through the full “breathing” crack than the partial “breathing” crack.
Again, the first round-trip echo in time domain is selected to calculate the linear parameters including the group velocity difference Δ c g   of the L(0, 1) mode between the damaged state and undamaged state, the reflection coefficient of the crack R f , and the location of the crack L o c f .
As shown in Figure 13, the group velocities of L(0, 1) modal UGWs in the damaged rod is lower than that in the undamaged rod at central exciting frequency of 30–70 kHz except for 50 KHz. As the crack area increases, the absolute differences between the group velocities of the damaged rod and the undamaged rod gradually increases. The reflection coefficient of the crack at different central exciting frequencies generally increases as the crack area increases (Figure 14). It is also worth noting that the reflection coefficient at 30 kHz and 40 kHz increases very slowly as the crack area increases (Figure 14), which indicates that UGWs at these frequencies are not sensitive to cracks. When the central exciting frequency is over 50 kHz, the reflection coefficient increases significantly after case C4 (Figure 14). Figure 15 shows that the UGWs at 50 kHz are more suitable for damage localization than the others.
The sum of residual square of the main frequency signal in time domain ARS, as expressed by Equation (10), is an important index for the change of signal in time domain due to damage. It is the accumulated square of the difference between the amplitudes of time-history signals in undamaged state and damaged state.
A R S = a ˜ 2 = t = 0 t = T I M E ( a B a s e , t a D a m , t ) 2
where a B a s e , t is the amplitude of the signal in the undamaged component at the central frequency at time t; and a D a m , t is the amplitude of the signal in the damaged component at central frequency at time t.
As shown in Figure 16, the ARS increases with the increasing of the crack area and the central frequency. In particular, the ARS increases very slightly with the increasing of the crack area at 30 kHz, while it increases much significantly with the increasing of the crack area at frequencies above 30 kHz, especially from the case C4 to CB6. Besides, the CB6 has a higher ARS than the C5 under the same conditions.
To explore the influence of the change of the “CAN” area on the nonlinear characteristics of the UGW, the short-time Fourier transform was performed to the signals and the result was shown in Figure 17.
It can be found that with the increasing of the crack area, the amplitude of the main frequency spectrum slightly decreases. On the contrary, the amplitude of the second order frequency increases and the SHW is significantly enhanced by the increasement of the crack area. It is noteworthy that the second-order frequency curve of CB6 is close to that of the undamaged case. This indicates that with the appearance of separated area, the nonlinearity of the crack decreases and thus the SHW reduces.
The peak amplitude in the central frequency f a band ( f a 5 kHz, f a + 5 kHz) and in the second-order frequency band (2 f a 5 kHz, 2 f a + 5 kHz) are denoted by   A 1 L and   A 1 H , respectively. The relative acoustic nonlinear parameter β and the relative amplitude of the second harmonic β r e are obtained using Equations (3) and (4), as shown in Figure 18 and Figure 19.
Figure 18 and Figure 19 show that, within the given range of central exciting frequencies, the nonlinear parameters β and β r e obviously increase from the undamaged case to C5 with the increasing of the nonlinear damage area of “CAN”. This finding indicates that the nonlinear parameters can be used to detect the nonlinear damage. However, for the case CB6, which has a same area of “CAN” to C5, the nonlinear indexes drop to a level close to the case C1. The observation reveals that the nonlinear indexes characterized the magnitude of the area of “CAN”, but not the whole area of crack. When f a = 70 kHz, the amplitudes of β r e and β are smaller than those at 60 kHz except for β in cases C1-C3, which shows that in general the existence of L(0, 2) mode reduce the nonlinear effect of guided waves traveling in rods. It is also worth noting that the nonlinear indexes β r e and β are more sensitive to the change of the “CAN” crack area at 60 kHz than at other frequencies.
According to the above numerical study, it can be found that when the central frequency of the excitation is 60 kHz, both the linear characteristics of UGWs, such as the reflection coefficient R f of the main frequency signals and the accumulated value ARS of the square residual of the main frequency signals, and the nonlinear characteristics of UGWs, such as the relative amplitude of the second harmonic β r e and the relative acoustic nonlinear parameter β , obviously vary with the change of the crack area. Therefore, the central frequency of excitation signal is set as 60 kHz in the following experimental study.

4. Experiment

4.1. Specimen and Setup

A specimen, the same as the threaded rod considered in the numerical study, was used in the experiment. The detailed geometry of it is shown in Figure 20. A groove-shaped notch with a length of about 13 mm, a width of about 1 mm, and a depth of about 2 mm was machined at the interface between the smooth segment and the threaded segment (Figure 20). A round hole with a diameter of about 1 mm and a depth of about 0.5 mm was machined in the middle of the groove notch. The groove notch and round hole were used as initial defects to initiate and accelerate the formation and growth of cracks.
The experimental setup is shown in Figure 21. The Instron8805 test machine was used to carry out the constant-amplitude tensile fatigue testing. A length of 90 cm of each end of the specimen was clamped by the test machine with a clamping force of 4.5 MPa. A piezoelectric probe, which works as both actuator and receiver in the plus-echo mode, was bonded at the end A of the threaded rod by using Loctite EA E-05CL adhesive. The piezoelectric probe was connected to the MSGW30 ultrasonic guided wave detector [44] through a data cable, and the detector was connected to the PC via WIFI.
A periodic sinusoidal tensile force whose maximum and minimum are 131.0 kN and 39.3 kN, respectively, was applied to the specimen by the test machine at a frequency of 8 Hz. This cyclic load caused a stress periodically varying from 68.6 MPa to 228.6 MPa in the specimen. After every 20 k (20,000) cycles of loading, the specimen was unloaded. Then sinusoidal waves with 3-cycle and 10-cycle Hanning windows at a central exciting frequency of 60 kHz were input as excitation signals in sequence. For each excitation signal, 20 echo signals were collected by the transducer at a sampling frequency of 2 MHz. Each sampling lasted for 0.004 s. During the experiment, the environment temperature was kept constant at 13 °C. Besides, a total of 80 initial signals were collected as baseline signals before cyclic loading.

4.2. Fatigue Crack Development

The vertical displacement of the upper clamp of the test machine throughout the experiment is shown in Figure 22. It shows that the amplitude of the displacement keeps constant before 330 k cycles of loading and then obviously rises till to 337.5 k cycles when the rod abruptly breaks at the cross-section with initial defect due to rapid growth of a fatigue crack.
The development of the fatigue crack during the cyclic load is shown in Figure 23. At the initial state, there was no visible crack around the initial defect (Figure 23a). After 320 k cycles of loading, a tiny crack through the artificial hole was observed (Figure 23b) and rapidly expanded (Figure 23c–e). The rod was broken after 337.5 k cycles of loading. The broken section was perpendicular to the longitudinal axis of the rod and crack expanding area and sudden breaking area can be recognized on the broken section (Figure 23f).

4.3. Experimental Results

In order to further remove the random noise, the guided wave signals obtained in the experiment are filtered by the continuous wavelet transform [45] and then normalized according to Equation (6). Figure 24 illustrates the filtered echo signal collected by the piezoelectric probe before the loading test, when sinusoidal waves with a 3-cycle Hanning window at 60 kHz are used as the excitation signal.
The reflection coefficient R f of the baseline signal is 0.078. The group velocity of the UGWs determined based on the baseline signal is 4548 m/s which is 2.1% smaller than the theoretically prediction and 1.4% smaller than the numerically prediction of the group velocity of the L(0, 1) modal UGWs. The geometric and material deviations of the specimen and the random vibrations of the bonding layer might cause the difference [46,47]. This general agreement in the group velocity verifies that the L(0, 1) mode is the main mode of the baseline signal.
Figure 25 shows the envelopes of the filtered fundamental echo signal and residual echo signal after various loading cycles. It is found that the amplitudes of the crack echo and its residual increase rapidly after 325 k cycles, while those of the end echo decreases rapidly. It is also found that the second echo changes more significantly with the increasing of load cycles than the first echo (Figure 25).
The reflection coefficient R f and the accumulative residual squares ARS of the echo signals are determined and plotted in Figure 26 and Figure 27 where the threshold for crack identification is set as the mean of the corresponding parameter determined from the baseline signals plus three times of the standard deviation.
It shows that both the reflection coefficient R f and the accumulative residual squares ARS exceed the threshold and monotonically increase after 320 k loading cycles. This verifies that both the reflection coefficient R f and the accumulative residual squares ARS are effective indexes for damage detection of the threaded rod. Moreover, the ARS slightly exceeds the threshold after 260 k loading cycles when there is no crack visually observed in the experiment. This indicates that the ARS may be used as an early warning index for the fatigue crack of threaded rods.
The frequency spectrums obtained by performing short-time Fourier transform on the echo signals are shown in Figure 28. In the early loading stage, there is no significant difference in the frequency spectrum between the test and baseline signal. After 335 k loading cycles, the amplitudes in the fundamental frequency band begin to shift and decrease obviously. The peaks in the second harmonic frequency band increase as the loading cycles increase. However, after 335 k loading cycles, the peaks of the second harmonic frequency spectrum decreases compared to that after 330 k loading cycles. The reason for this phenomenon is that the wider opening of the crack after 335 k loading cycles reduces the “CAN” area, which is similar to the case CB6 in the numerical study.
The relative amplitude of the second harmonic β r e and the relative acoustic nonlinear parameter β of the second harmonic are shown in Figure 29 and Figure 30. The thresholds for these two indexes are also defined as the mean of the corresponding parameter determined from the baseline signals plus three times of the standard deviation. It can be found that β r e and β stably exceed their thresholds after 160 k loading cycles and 240 k loading cycles, respectively. Then, they are fluctuating upward and reach their peaks at 320 k loading cycles. After that, they continuously decrease as the number of loading cycles increases, and fall below the threshold at 335 k loading cycles. It can be seen that the nonlinear indexes β r e and β are able to warn the damage much earlier than the linear indexes R f and ARS.

4.4. Influence of Number of Modulation Cycles

To investigate the influence of the number of modulation cycles on the echo signals, sinusoidal waves with 10-cycle Hanning window at 60 kHz was also used as excitation signal. With the increasing modulation cycles of the excitation signals, the spectral energy is more concentrated and the bandwidth is smaller, which is conducive to the capture of nonlinear characteristics based on the spectrum. However, it also means longer duration of the excitation signal, causing the echo signals to be overlapped by the excitation signals, which makes it is difficult to obtain linear characteristics. Therefore, only the nonlinear characteristics are investigated in this chapter.
The frequency spectrums corresponding to the 10-cycle excitation signal are shown in Figure 31. Comparing Figure 31a with Figure 28a, it can be noted that, with the increment of the number of modulation cycles, the energy of UGWs is more concentrated to the central frequency, the amplitudes of the fundamental spectrums increase greatly, and the frequency shift becomes more obvious. Comparing Figure 31b with Figure 28b, it can be found that the amplitudes of the second harmonic frequency spectrums also increase, but not as much as the fundamental spectrums.
The nonlinear indexes β r e and β corresponding to the second harmonic corresponding to the 10-cycle excitation signal are shown in Figure 32 and Figure 33. Comparing Figure 32 and Figure 33 with Figure 29 and Figure 30, respectively, it is notable that the variation trends of β r e and β under the two excitation signals with different number of modulation cycles are similar. It can also be found that β r e and β corresponding to the 10-cycle excitation signal are more stable than those corresponding to the 3-cycle excitation signal, especially after their first pass of the thresholds. Meanwhile, both β r e and β corresponding to the 10-cycle excitation signal stably exceed the thresholds after 200 k cycles, which is earlier than those corresponding to the 3-cycle excitation signal. It is also worth noting that both β r e and β corresponding to the 10-cycle excitation signal do not fall below the thresholds after 335 k cycles. This finding indicates that the nonlinear indexes become more sensible to wider cracks, as the number of modulation cycles increase.

5. Conclusions

This paper investigates the feasibility of UGW-based fatigue crack detection of threaded rods through analytical, numerical and experimental studies. The main conclusions that can be drawn on the basis of the studies are given as follows:
(1)
The longitudinal L(0, 1) modal guided wave has a simple displacement distribution at low frequency when traveling in the rod. It is easy to be excited and the energy of it is concentrated in the axial direction of the rod. Hence, it is suitable for damage detection of rods.
(2)
Via numerical simulations, it is found that the thread-induced cross-sectional changes will cause thread echoes. The linear and nonlinear characteristics of the guided wave obviously change with the growth of cracks, which shows a great potential to use these indexes for damage detection. It is also found that these indexes generally have the highest sensitivity to cracks at central exciting frequency of 60 kHz in the experimental frequency region.
(3)
The experiment on UGW-based fatigue crack detection of a threaded rod under cyclic tensile load shows that the reflection coefficient R f is able to warn the crack when it is visible and the accumulative residual squares ARS is able to warn the crack before it is visible. The spectrum-based nonlinear damage indexes, i.e., the relative amplitude of second harmonic β r e and the relative acoustic nonlinear parameter β , are generally able to give an earlier warning on the fatigue crack than the linear indexes R f and ARS.
(4)
Increasing the number of modulation cycles of excitation signals from 3 to 10 improves the stability of the nonlinear indexes β r e and β , and their sensibility to wider cracks.
It should be pointed out that it is mainly a proof-of-concept study based on an ideal case in laboratory. If an in-service threaded rod is considered, there may be impossible to obtain baseline signal from undamaged state and there may be more noise in the echo signals due to the imperfections of the rod and the environment. Improving the proposed approach to be applicable to in-service threaded rods is the ultimate goal of the authors and will be the topic of the authors’ future work.

Author Contributions

Conceptualization, X.X.; Formal analysis, K.P.; Funding acquisition, X.X.; Investigation, K.P. and Y.Z.; Methodology, K.P.; Project administration, X.X. and Y.L.; Resources, X.X. and Y.L.; Software, K.P.; Supervision, X.X., J.H. and Y.L.; Validation, K.P.; Visualization, K.P.; Writing—original draft, K.P. and Y.Z.; Writing—review & editing, X.X., J.H. and Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China under grant 2021YFF0501002 and by the National Natural Science Foundation of China under grant U21A20462 and 61872285.

Acknowledgments

This work was supported by Center for Balance Architecture of Zhejiang University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Casanova, F.; Mantilla, C. Fatigue Failure of the Bolts Connecting a Francis Turbine with the Shaft. Eng. Fail. Anal. 2018, 90, 1–13. [Google Scholar] [CrossRef]
  2. Bakunov, A.S.; Kurozaev, V.P.; Kudryavtsev, D.A.; Bronnikov, V.K.; Kravchenko, V.G. Increasing the Reliability of Magnetic-Particle Testing by Means of a yM∧K-10 Automated Unit for Magnetic Fluorescent-Penetrant Inspection of Pipe End Faces. Russ. J. Nondestruct. Test. 2004, 40, 311–316. [Google Scholar] [CrossRef]
  3. The Chinese Society of NDT. Magnetic Particle Testing, 2nd ed.; China Machine Press: Beijing, China, 2004. [Google Scholar]
  4. Chen, G.; Zhang, W.; Zhang, Z.; Jin, X.; Pang, W. A New Rosette-like Eddy Current Array Sensor with High Sensitivity for Fatigue Defect around Bolt Hole in SHM. NDT E Int. 2018, 94, 70–78. [Google Scholar] [CrossRef]
  5. Im, K.-H.; Lee, S.-G.; Kim, H.-J.; Song, S.-J.; Woo, Y.-D.; Ra, S.-W.; Lee, H.-H. Experimental Approach and Simulation-Based Design of Eddy Current Sensors for Inspection of Vehicular Bolts. J. Korean Soc. Manuf. Technol. Eng. 2015, 24, 294–301. [Google Scholar] [CrossRef]
  6. Yao, Y.; Tung, S.-T.; Glisic, B. Crack Detection and Characterization Techniques-An Overview. Struct. Control Health Monit. 2014, 21, 1387–1413. [Google Scholar] [CrossRef]
  7. Klepka, A.; Staszewski, W.J.; Jenal, R.B.; Szwedo, M.; Iwaniec, J.; Uhl, T. Nonlinear Acoustics for Fatigue Crack Detection Experimental Investigations of Vibro-acoustic Wave Modulations. Struct. Health Monit. 2012, 11, 197–211. [Google Scholar] [CrossRef]
  8. Shifrin, E.I. Identification of a Finite Number of Small Cracks in a Rod Using Natural Frequencies. Mech. Syst. Signal Process. 2016, 70, 613–624. [Google Scholar] [CrossRef]
  9. Shifrin, E.I. Inverse Spectral Problem for a Non-uniform Rod with Multiple Cracks. Mech. Syst. Signal Process. 2017, 96, 348–365. [Google Scholar] [CrossRef]
  10. Tarunpreet, S.; Shankar, S. Damage Identification Using Vibration Monitoring Techniques. Mater. Today Proc. 2022. [Google Scholar] [CrossRef]
  11. Suh, D.M.; Kim, W.W. A New Ultrasonic Technique for Detection and Sizing of Small Cracks in Studs and Bolts. J. Nondestruct. Eval. 1995, 14, 201–206. [Google Scholar] [CrossRef]
  12. Wagle, S.; Kato, H. Ultrasonic Evaluation of Fatigue Crack Appearing at Bolt Joint of Aluminum Alloy Plate. In Proceedings of the JSME Annual Meeting, Tokyo, Janpan, 2 August 2008; pp. 443–444. [Google Scholar] [CrossRef]
  13. DL/T 694-2012; The Technical Guide of Ultrasonic Inspection for High-temperature Tight Bolts. China Electric Power Press: Beijing, China, 2012.
  14. Wei, S. Study on Ultrasonic Guided Wave Testing Technology for Wind Power Bolts. Master Thesis, Hubei University of Technology, Wuhan, China, 2021. (In Chinese). [Google Scholar]
  15. Mitra, M.; Gopalakrishnan, S. Guided Wave Based Structural Health Monitoring: A Review. Smart Mater. Struct. 2016, 25, 053001. [Google Scholar] [CrossRef]
  16. He, C.; Zheng, M.; Lv, Y.; Deng, P.; Zhao, H.; Liu, X.; Song, G.; Liu, Z.; Jiao, J.; Wu, B. Development, Applications and Challenges in Ultrasonic Guided Waves Testing Technology. Chin. J. Sci. Instrum. 2016, 37, 1713–1735. [Google Scholar] [CrossRef]
  17. Niu, X.; Marques, H.R.; Chen, H.P. Sensitivity Analysis of Circumferential Transducer Array with T(0,1) Mode of Pipes. Smart Struct. Syst. 2018, 21, 761–776. [Google Scholar] [CrossRef]
  18. Garg, M.; Sharma, S.; Sharma, S.; Mehta, R. Non-contact Damage Monitoring Technique for FRP Laminates Using Guided Waves. Smart Struct. Syst. 2016, 17, 795–817. [Google Scholar] [CrossRef]
  19. Schaal, C.; Bischoff, S.; Gaul, L. Damage Detection in Multi-wire Cables Using Guided Ultrasonic Waves. Struct. Health Monit. 2016, 15, 279–288. [Google Scholar] [CrossRef]
  20. Wang, K.; Cao, W.; Su, Z.; Wang, P.; Zhang, X.; Chen, L.; Guan, R.; Lu, Y. Structural Health Monitoring of High-speed Railway Tracks Using Diffuse Ultrasonic Wave-based Condition Contrast: Theory and Validation. Smart Struct. Syst. 2020, 26, 227–239. [Google Scholar] [CrossRef]
  21. Shoji, M.; Higashi, Y. Fundamental Study on Guided Wave Testing of Cylindrical Bars Embedded in Soil. In Proceedings of the 2014 IEEE International Ultrasonics Symposium, Chicago, IL, USA, 3–6 September 2014; pp. 1404–1407. [Google Scholar] [CrossRef]
  22. Shoji, M.; Hirata, A. Ultrasonic Guided Wave Testing of Anchor Rods Embedded in Soil. In Proceedings of the 2016 IEEE International Ultrasonics Symposium, Tours, Frace, 18–21 September 2016. [Google Scholar]
  23. Stepinski, T. Time-Frequency Analysis of Guided Ultrasonic Waves Used for Assessing Integrity of Rock Bolts. In Smart Materials and Nondestructive Evaluation for Energy Systems; SPIE: Portland, OR, USA, 2017; Volume 10171, p. 101710J1-10. [Google Scholar] [CrossRef]
  24. Rong, X.; Lin, P.; Liu, J.; Yang, T. A New Approach of Waveform Interpretation Applied in Nondestructive Testing of Defects in Rock Bolts Based on Mode Identification. Math. Probl. Eng. 2017, 2017, 7920649. [Google Scholar] [CrossRef] [Green Version]
  25. Zhao, J.; Durham, N.; Abdel-Hadi, K.; McKenzie, C.A.; Thomson, D.J. Acoustic Guided Wave Techniques for Detecting Corrosion Damage of Electrical Grounding Rods. Measurement 2019, 147, 106858. [Google Scholar] [CrossRef]
  26. Amjad, U.; Yadav, S.K.; Kundu, T. Detection and Quantification of Diameter Reduction Due to Corrosion in Reinforcing Steel Bars. Struct. Health Monit. 2015, 14, 532–543. [Google Scholar] [CrossRef]
  27. He, C.; Sun, Y.; Wu, B.; Wang, X.; Liu, Z. Application of Ultrasonic Guided Waves Technology to Inspection of Bolt Embedded in Soils. Chin. J. Geotech. Eng. 2006, 28, 1144–1147. [Google Scholar]
  28. Beard, M.D.; Lowe, M.J.S.; Cawley, P. Development of a Guided Wave Inspection Technique for Rock Bolts. AIP Conf. Proc. 2002, 615, 1318–1325. [Google Scholar] [CrossRef]
  29. Tsai, Y.; Liu, S.; Light, G. Crack-Depth Effects in the Cylindrically Guided Wave Technique for Bolt and Pump-Shaft Inspections. Nucl. Eng. Des. 1992, 133, 77–80. [Google Scholar] [CrossRef]
  30. Du, Y.; Ouyang, Q.; Zhou, F. Study on Stress Concentration Factor of Bolt Screw. Eng. Mech. 2014, 31, 174–180. [Google Scholar]
  31. Su, Z.; Zhou, C.; Hong, M.; Cheng, L.; Wang, Q.; Qing, X. Acousto-Ultrasonics-Based Fatigue Damage Characterization: Linear versus Nonlinear Signal Features. Mech. Syst. Signal Process. 2014, 45, 225–239. [Google Scholar] [CrossRef]
  32. Rose, J.L. Ultrasonic Guided Waves in Solid Media; Cambridge University Press: New York, NY, USA, 2014. [Google Scholar]
  33. Hayashi, T.; Song, W.-J.; Rose, J.L. Guided Wave Dispersion Curves for a Bar with an Arbitrary Cross-Section, a Rod and Rail Example. Ultrasonics 2003, 41, 175–183. [Google Scholar] [CrossRef]
  34. Wu, B.; Qi, W.; He, C.; Zhou, J. Recognition of Defects on Steel Rod Using Ultrasonic Guided Waves Based on Neural Network. Eng. Mech. 2013, 30, 470–476. [Google Scholar]
  35. Jhang, K.-Y. Nonlinear Ultrasonic Techniques for Nondestructive Assessment of Micro Damage in Material: A Review. Int. J. Precis. Eng. Manuf. 2009, 10, 123–135. [Google Scholar] [CrossRef]
  36. Broda, D.; Staszewski, W.J.; Martowicz, A.; Uhl, T.; Silberschmidt, V.V. Modelling of Nonlinear Crack-Wave Interactions for Damage Detection Based on Ultrasound—A Review. J. Sound Vib. 2014, 333, 1097–1118. [Google Scholar] [CrossRef]
  37. Solodov, I.Y.; Krohn, N.; Busse, G. CAN: An Example of Nonclassical Acoustic Nonlinearity in Solids. Ultrasonics 2002, 40, 621–625. [Google Scholar] [CrossRef]
  38. Lee, Y.F.; Lu, Y.; Guan, R. Nonlinear Guided Waves for Fatigue Crack Evaluation in Steel Joints with Digital Image Correlation Validation. Smart Mater. Struct. 2020, 29, 035031. [Google Scholar] [CrossRef]
  39. Wang, K.; Li, Y.; Su, Z.; Guan, R.; Lu, Y.; Yuan, S. Nonlinear Aspects of “Breathing” Crack-Disturbed Plate Waves: 3-D Analytical Modeling with Experimental Validation. Int. J. Mech. Sci. 2019, 159, 140–150. [Google Scholar] [CrossRef]
  40. Smith, M. ABAQUS/Standard User’s Manual, Version 6.9; Dassault Systèmes Simulia Corp: Providence, RI, USA, 2009. [Google Scholar]
  41. Moser, F.; Jacobs, L.J.; Qu, J. Modeling Elastic Wave Propagation in Waveguides with the Finite Element Method. NDT E Int. 1999, 32, 225–234. [Google Scholar] [CrossRef]
  42. Datta, D.; Kishore, N.N. Features of Ultrasonic Wave Propagation to Identify Defects in Composite Materials Modelled by Finite Element Method. NDT E Int. 1996, 29, 213–223. [Google Scholar] [CrossRef]
  43. Alleyne, D.; Cawley, P. A Two-Dimensional Fourier Transform Method for the Measurement of Propagating Multimode Signals. J. Acoust. Soc. Am. 1991, 89, 1159–1168. [Google Scholar] [CrossRef]
  44. Zheda Jingyi. Available online: http://www.jingyitech.com/ (accessed on 13 November 2021).
  45. Hernandez-Salazar, C.D.; Herrejon, A.B.; Sanchez, J.I.A. Damage Detection in Multi-Wire Cables Using Continuous Wavelet Transform Analysis of Ultrasonic Guided Waves. In Proceedings of the 2009 Electronics, Robotics and Automotive Mechanics Conference (CERMA), Cernavaca, Mexico, 22–25 September 2009; pp. 250–255. [Google Scholar] [CrossRef]
  46. Tu, J.Q.; Tang, Z.F.; Yun, C.B.; Wu, J.J.; Xu, X. Guided Wave-based Damage Assessment on Welded Steel I-beam under Ambient Temperature Variations. Struct. Control Health Monit. 2021, 28, e2696. [Google Scholar] [CrossRef]
  47. Zhang, P.F.; Tang, Z.F.; Duan, Y.F.; Yun, C.B.; Lv, F.Z. Ultrasonic Guided Wave Approach Incorporating SAFE for Detecting Wire Breakage in Bridge Cable. Smart Struct. Syst. 2018, 22, 481–493. [Google Scholar] [CrossRef]
Figure 1. Infinitely long cylindrical rod with diameter of r.
Figure 1. Infinitely long cylindrical rod with diameter of r.
Sensors 22 06885 g001
Figure 2. Dispersion curves of a typical infinitely long cylindrical rod with diameter of 27 mm.
Figure 2. Dispersion curves of a typical infinitely long cylindrical rod with diameter of 27 mm.
Sensors 22 06885 g002
Figure 3. Displacement distribution of L(0, 1) modal UGWs in cylindrical rod at 60 kHz.
Figure 3. Displacement distribution of L(0, 1) modal UGWs in cylindrical rod at 60 kHz.
Sensors 22 06885 g003
Figure 4. Numerical model of threaded rod.
Figure 4. Numerical model of threaded rod.
Sensors 22 06885 g004
Figure 5. Excitation signal with central frequency of 60 kHz. (a) Excitation signal in time domain. (b) Excitation signal in frequency domain.
Figure 5. Excitation signal with central frequency of 60 kHz. (a) Excitation signal in time domain. (b) Excitation signal in frequency domain.
Sensors 22 06885 g005
Figure 6. Time history and envelope of UGWs with central frequency of 60 kHz propagating in threaded rod.
Figure 6. Time history and envelope of UGWs with central frequency of 60 kHz propagating in threaded rod.
Sensors 22 06885 g006
Figure 7. Group velocity under different excitation frequencies.
Figure 7. Group velocity under different excitation frequencies.
Sensors 22 06885 g007
Figure 8. Location of refection cross-section under different excitation frequencies.
Figure 8. Location of refection cross-section under different excitation frequencies.
Sensors 22 06885 g008
Figure 9. Thread reflection coefficients at different excitation frequencies.
Figure 9. Thread reflection coefficients at different excitation frequencies.
Sensors 22 06885 g009
Figure 10. Wavenumber diagram based on 2D Fourier transform at 60 kHz.
Figure 10. Wavenumber diagram based on 2D Fourier transform at 60 kHz.
Sensors 22 06885 g010
Figure 11. Modeling of “breathing” crack. (a) “CAN” area of fully “breathing” crack. (b) “CAN” area of partially “breathing” crack.
Figure 11. Modeling of “breathing” crack. (a) “CAN” area of fully “breathing” crack. (b) “CAN” area of partially “breathing” crack.
Sensors 22 06885 g011
Figure 12. Time domain envelopes of signals at 60 kHz. (a) Time domain envelopes of full signals. (b) Time domain envelopes of residual signals.
Figure 12. Time domain envelopes of signals at 60 kHz. (a) Time domain envelopes of full signals. (b) Time domain envelopes of residual signals.
Sensors 22 06885 g012
Figure 13. L(0, 1) mode group velocity difference.
Figure 13. L(0, 1) mode group velocity difference.
Sensors 22 06885 g013
Figure 14. Crack reflection coefficient diagram.
Figure 14. Crack reflection coefficient diagram.
Sensors 22 06885 g014
Figure 15. Crack location.
Figure 15. Crack location.
Sensors 22 06885 g015
Figure 16. Sum of residual square of main frequency signal in time domain.
Figure 16. Sum of residual square of main frequency signal in time domain.
Sensors 22 06885 g016
Figure 17. Spectrum at central exciting frequency of 60 kHz. (a) Main frequency spectrum. (b) Second-order frequency spectrum.
Figure 17. Spectrum at central exciting frequency of 60 kHz. (a) Main frequency spectrum. (b) Second-order frequency spectrum.
Sensors 22 06885 g017
Figure 18. Relative acoustic nonlinear parameter β at different cases.
Figure 18. Relative acoustic nonlinear parameter β at different cases.
Sensors 22 06885 g018
Figure 19. Relative amplitude of second harmonic β r e at different cases.
Figure 19. Relative amplitude of second harmonic β r e at different cases.
Sensors 22 06885 g019
Figure 20. Geometry of specimen.
Figure 20. Geometry of specimen.
Sensors 22 06885 g020
Figure 21. Photograph of experimental setup.
Figure 21. Photograph of experimental setup.
Sensors 22 06885 g021
Figure 22. Displacement of upper clamp of test machine.
Figure 22. Displacement of upper clamp of test machine.
Sensors 22 06885 g022
Figure 23. Crack growth and broken at defected cross-section. (a1) left side of hole; (a2) hole bottom; (a3) right side of hole; (b1) left side of hole; (b2) hole bottom; (b3) right side of hole; (c1) left side of hole; (c2) hole bottom; (c3) right side of hole; (d1) left side of hole; (d2) hole bottom; (d3) right side of hole; (e1) left side of hole; (e2) hole bottom; (e3) right side of hole; (f1) Upper broken surface (close to end D); (f2) Lower broken surface (close to end A).
Figure 23. Crack growth and broken at defected cross-section. (a1) left side of hole; (a2) hole bottom; (a3) right side of hole; (b1) left side of hole; (b2) hole bottom; (b3) right side of hole; (c1) left side of hole; (c2) hole bottom; (c3) right side of hole; (d1) left side of hole; (d2) hole bottom; (d3) right side of hole; (e1) left side of hole; (e2) hole bottom; (e3) right side of hole; (f1) Upper broken surface (close to end D); (f2) Lower broken surface (close to end A).
Sensors 22 06885 g023aSensors 22 06885 g023b
Figure 24. Filtered echo signal at initial baseline state.
Figure 24. Filtered echo signal at initial baseline state.
Sensors 22 06885 g024
Figure 25. Envelopes of fundamental and residual echo signals after various of loading cycles. (a) Envelope of echo signal; (b) Envelope of residual echo signal.
Figure 25. Envelopes of fundamental and residual echo signals after various of loading cycles. (a) Envelope of echo signal; (b) Envelope of residual echo signal.
Sensors 22 06885 g025
Figure 26. Reflection coefficient R f of filtered echo signal.
Figure 26. Reflection coefficient R f of filtered echo signal.
Sensors 22 06885 g026
Figure 27. Accumulative residual squares ARS of filtered echo signal.
Figure 27. Accumulative residual squares ARS of filtered echo signal.
Sensors 22 06885 g027
Figure 28. Frequency spectrums after various loading cycles. (a) Fundamental frequency spectrums around central frequency; (b) Second harmonic frequency spectrums around double of central frequency.
Figure 28. Frequency spectrums after various loading cycles. (a) Fundamental frequency spectrums around central frequency; (b) Second harmonic frequency spectrums around double of central frequency.
Sensors 22 06885 g028
Figure 29. Relative amplitude of second harmonic β r e versus loading cycles.
Figure 29. Relative amplitude of second harmonic β r e versus loading cycles.
Sensors 22 06885 g029
Figure 30. Relative acoustic nonlinear parameter β versus loading cycles.
Figure 30. Relative acoustic nonlinear parameter β versus loading cycles.
Sensors 22 06885 g030
Figure 31. Frequency spectrums after various loading cycles corresponding to 10-cycle excitation. (a) Fundamental frequency spectrums around central frequency; (b) Second-order frequency spectrums around double of central frequency.
Figure 31. Frequency spectrums after various loading cycles corresponding to 10-cycle excitation. (a) Fundamental frequency spectrums around central frequency; (b) Second-order frequency spectrums around double of central frequency.
Sensors 22 06885 g031
Figure 32. Relative amplitude of second harmonic β r e versus loading cycles corresponding to 10-cycle excitation signal.
Figure 32. Relative amplitude of second harmonic β r e versus loading cycles corresponding to 10-cycle excitation signal.
Sensors 22 06885 g032
Figure 33. Relative acoustic nonlinear parameter β versus loading cycles corresponding to 10-cycle excitation signa.
Figure 33. Relative acoustic nonlinear parameter β versus loading cycles corresponding to 10-cycle excitation signa.
Sensors 22 06885 g033
Table 1. Micro-crack settings of different areas.
Table 1. Micro-crack settings of different areas.
Working ConditionC1C2C3C4C5CB6
“CAN” area (mm2)11.126.136.549.266.911.1
Separated area (mm2)0.00.00.00.00.055.9
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Peng, K.; Zhang, Y.; Xu, X.; Han, J.; Luo, Y. Crack Detection of Threaded Steel Rods Based on Ultrasonic Guided Waves. Sensors 2022, 22, 6885. https://doi.org/10.3390/s22186885

AMA Style

Peng K, Zhang Y, Xu X, Han J, Luo Y. Crack Detection of Threaded Steel Rods Based on Ultrasonic Guided Waves. Sensors. 2022; 22(18):6885. https://doi.org/10.3390/s22186885

Chicago/Turabian Style

Peng, Kunhong, Yi Zhang, Xian Xu, Jinsong Han, and Yaozhi Luo. 2022. "Crack Detection of Threaded Steel Rods Based on Ultrasonic Guided Waves" Sensors 22, no. 18: 6885. https://doi.org/10.3390/s22186885

APA Style

Peng, K., Zhang, Y., Xu, X., Han, J., & Luo, Y. (2022). Crack Detection of Threaded Steel Rods Based on Ultrasonic Guided Waves. Sensors, 22(18), 6885. https://doi.org/10.3390/s22186885

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