1. Introduction
Eddy current testing (ECT) has been widely applied in the detection of surface and subsurface defects in conductive materials. It utilizes a magnetic probe to generate an alternating magnetic field, which induces an eddy current in the specimen and measures the secondary magnetic field produced by the eddy current. Physical anomalies and discontinuities in material properties existing in the specimen can then be identified by analyzing the probe output signals [
1]. The ECT probe plays roles as the actuating source and also the receiving sensor for damage detection. Throughout the literature, different ECT probe types have been reported and their operation modes are generally classified into four categories—absolute, differential, reflection and hybrid [
1,
2,
3]. Reflection probes have separate excitation and detection sensors, which can be individually optimized for their intended applications [
3,
4]. In the past few decades, modern magnetic sensors such as Hall devices; anisotropic, giant and tunnel magneto resistive (AMR, GMR and TMR) sensors and superconducting quantum interference devices (SQUIDs) are becoming alternatives of detection coils due to their intrinsic advantages over induction coils, including their higher sensitivity at low frequency and smaller sensing footprint for small defects [
5,
6,
7]. However, for the excitation purpose, a wire-wound coil, no matter air-cored, ferrite-cored or printed, is still the primary and most common choice.
As the ECT probe’s capability and accuracy for defect detection correlate with the noise level of eddy current signals, one main objective in the ECT probe design is to achieve an output with high signal-to-noise ratio (SNR) [
8,
9,
10]. This can be achieved by improving the probe coil configurations and excitation modes, examples of which have been extensively reported in previous studies [
2,
10,
11,
12,
13].
Following the notion of the figure-of-eight coil used in deep transcranial magnetic stimulation [
14,
15], a so-called “figure-8-shaped” focusing probe was introduced by Zhang et al. [
16,
17]. It consists of two coils tilted oppositely like a butterfly over the specimen. Compared with the circular pancake coil, it can induce unidirectional eddy currents with greater current density and focus them in a smaller region in the conductor. Their study also indicates that the eddy current density reaches its maximum beneath the probe center at which the eddy current produced by circular pancake probes tends to be zero. This type of probe shows great sensitivity to the defect perpendicular to the eddy current but would lose its detectability when defects orient in parallel along the eddy current direction [
18].
To cope with the practical situation that the defect orientation is not known in advance or is unpredictable, two methods have been developed—rotationally scanning on the suspicious zone [
19,
20,
21] and developing a rotating current excitation [
22,
23,
24,
25,
26]. They have the same underlying physics for the sensing mechanism but the latter is superior to the former because it avoids the inevitable noise caused by mechanically rotating and is suitable for a C-scan or even a line scan. In general, this method employs two orthogonally placed coils driven by two current excitations with 90 degrees phase difference to result in a uniform and rotating eddy current (EC) field in the specimen. By doing this, the probe shows the same sensitivity to arbitrarily oriented defects. Yang et al. [
22] proposed a printed-circuit-board (PCB)-based rotating field EC-GMR probe for the detection of embedded cracks in riveted multilayer structures. Each side of the PCB was fabricated with a planar coil, resulting in a structural lift-off deviation for the two excitation coils. Thus, a compensation method was required to form a circularly rotating field. This probe structure was then further developed by Ye et al. [
23,
24], including GMR arrays to eliminate the impact of background field and applying TMR to replace GMR for increasing the spatial resolution. In a similar way, Repelianto et al. [
25] presented a rotating uniform eddy current probe with two pairs of orthogonally installed rectangular excitation coils and a small detection coil. Since the stacked excitation coils have quite different lift-off distances, they are well adjusted with different numbers of turns and excitation currents to make a rotating field. Recently, Wang et al. [
26] presented a probe with a printed right-triangle excitation coil and two TMR arrays placed on the triangle’s two legs. In order to make an equivalent rotating field, acquired TMR data were fused by a series of complex algorithms including shifting, rotating and multifrequency mixing. In short, to generate a rotating field, additional processing like excitation adjustment or data fusion after acquiring the EC signals is required for state-of-the-art EC probe design.
This work proposed a novel rotating field EC probe by integrating the figure-8-shaped coil focusing scheme. Two pairs of focusing coils are orthogonally arranged at the same lift-off, which thus avoids additional excitation adjustment for generating a rotating EC field and also makes the field concentrated with high density. A prototype probe with a GMR sensor as the receiver was developed and its capacity for arbitrary orientation defects was investigated using finite element analysis and an experiment test. The GMR sensor was mounted at the probe center without increasing the probe’s design lift-off. The experimental signal had superior SNR and its amplitude can be used to evaluate the depth of arbitrarily oriented defect, demonstrating the effectiveness of the proposed probe.
3. Model-based Study of Defect Detection
In this section, 3D finite element analysis (FEA) is performed by using ANSYS Workbench software to study the effectiveness of the presented probe for detecting arbitrary orientation defects. The unidirectional focusing probe is also studied for giving comparative results. The Solid236 element with 20 nodes is selected to model all the entities including the probe, the specimen and the air region. It is capable of modeling electromagnetic fields based on the edge-based magnetic vector potential method. In this method, the magnetic vector potential
and electric scalar potential
are used to represent the electromagnetic field. The basic equations involved are expressed as Maxwell–Faraday’s law, Maxwell–Ampere’s law, constitutive relations and potential functions:
where the various quantities are defined as
: Electric field intensity (V/m)
: Magnetic flux density (T)
: Magnetic field intensity (A/m)
: Total conduction current density (A/m2)
: Applied source current density (A/m2)
: Induced eddy current density (A/m2)
: Velocity current density (A/m2)
: Time (s)
: Permittivity (F/m)
: Magnetic permeability (H/m)
: Conductivity (S/m)
Substituting Equation (8) into Equation (4) results in the following equation:
Invoking Equation (9),
can be written as
Finally, by substituting Equations (11), (6) and (7) into Equation (5), the governing equations are obtained:
Figure 3a,b depict the FEA model. For clarity, the elements of air properties are not plotted except those inside the defect. The geometry parameters are defined here—2
β, the flare angle between two coils of the sub-probe;
lc, the distance from the coil center to the specimen surface;
dc, the center distance between two coils;
θ, the angle from the symmetric line of one sub-probe to the defect lengthwise direction.
The specimen under test is a carbon steel plate with parameters listed in
Table 1. A groove with the size of 10 mm × 2 mm × 2 mm (length, width, depth) is set on the central top surface of the specimen to simulate the defect.
Table 2 lists the geometric parameters of the probe coil.
Figure 4 shows the distribution of the induced eddy current on the surface of the specimen with the rotating field probe. The eddy current under the probe center is oriented towards varying directions as the time changes. It is dynamically formed in a counterclockwise direction in a whole excitation cycle
T and indeed exhibits a spatial rotating feature. This observation is consistent with the prediction described in
Figure 2b and therefore proves the effectiveness of the rotating field probe.
Figure 5 compares the distribution of normal magnetic flux density on the specimen surface with the groove lain in eddy currents induced by the focusing probe and the rotating field probe.
Figure 5a and b respectively show the cases with eddy current flow parallel and perpendicular to the groove length direction. It can be seen that the normal magnetic field only arises at the two edges along the current flow. In contrast, in the rotating field as shown in
Figure 5c, the normal magnetic field is produced alongside all four edges and with uniform density distribution, which provides omni-directional information about the groove and therefore makes the groove able to be sensed from any orientation. Furthermore, by checking the color variation against the legend, the magnetic flux density is found to decrease sharply from the groove edge to its outer vicinity, which is attributed to the focusing effect and can contribute towards enhancing the defect-related signal and suppress the influence of the background field.
The simulated normal magnetic flux density
is a complex quantity and can be expressed in terms of the magnitude
and the phase
as follows:
where Re and Im are the real and imaginary part of
, respectively.
3.1. Detection of Arbitrary Orientation Defects
To simulate different defect orientations in a FEA model, the probe is rotated step-by-step to change the value of angle θ. In this way, the meshes of the defect region can be kept exactly the same in different simulation models, thus ensuring the comparability of the simulation results. Considering the structure symmetry, a set of θ from 0 to 90 degrees by steps of 15 degrees is used for the focusing probe model. For rotating the focused field ECT probe, a set of θ from 0 to 165 degrees (θ = 180° is the same as the case of 0°) by steps of 15 degrees is used.
Simulation data is processed by three steps to obtain the signal. First, a 30 mm long data path along the defect length direction is aligned, with 2.5 mm distance off the plate and middle point above the defect center. Data sets representing the normal magnetic flux density mapped to this path are extracted. Afterwards, all the data are subtracted by a data set acquired from a defect-free model, yielding a differential signal with the notation diff. . Finally, the magnitude and phase angle of diff. are calculated by using formulas (14) and (15).
The excitation frequency is a key parameter in eddy current testing. For a given conductive material, the frequency is not only concerned with the skin depth but also with the probe signal magnitude. An optimal frequency exists at which the maximum signal is retrieved for a fixed-size defect [
27].
Figure 6 gives the simulated results by using the rotating probe excited by harmonic currents with the same density of 1.67 × 10
6 A/m
2 but different frequencies. As seen in
Figure 6a, the magnitude signal curve exhibits two symmetric peaks whose positions correspond to the two lengthwise ends of the defect. This phenomenon agrees with the fact that the defect edge causes a severe distortion to eddy current flow and thus triggers a significant normal magnetic field.
Figure 6b shows that the probe signal magnitude reaches the maximum at the frequency of 1 kHz as the frequency varies from 1 Hz to 100 kHz. Considering the need for detecting a deeper defect, a compromise is made and the 100 Hz is selected as the excitation frequency.
The simulated signal magnitudes of the focusing and the rotating probes for detecting defects with various orientations are shown in
Figure 7a,b, respectively. As the defect orientation varies (i.e.,
θ increases), the peak value of the focusing probe signal decreases significantly, while that of the rotating focused field probe remains almost unchanged. From this comparison, it can be concluded that the strategy of integrating the focusing coil into the rotating field probe works well and the hybrid probe shows uniform sensitivity to arbitrary orientation defects in the specimen.
Further, the signal data with maximum peaks, arising at the defect edge as analyzed above, is extracted and its phase angles are calculated.
Figure 8 shows the variation of phase angle with the defect orientation for using the rotating probe. The fitted line indicates a descending linear relationship between them, which suggests a potential method for determining the defect orientation.
3.2. Effect of Flare Angle
As an important factor, the effect of the flare angle 2
β on the performance of the rotating focused field probe is studied.
Figure 9 shows the simulation result for 2
β varying from 0 to 360 degrees by steps of 15 degrees. Note the magnitude represents the curve’s peak magnitude such as that plotted in
Figure 7. It can be seen that the signal magnitude slowly increases to reach its maximum and afterwards decreases until the end. For the given probe, the peak value appears at 2
β = 150°. This peak point is not fixed but depends on the coil geometry. As the coil height
h decreases from 18 mm to 5 mm, the flare angle corresponding to the peak point (circled on the curve) shows a gradual decrease as well. The curve of
h = 18 mm does not have values at 2
β = 0° and 360° because in either case the adjacent coils of the probe will interfere with each other. When the coil height increases, the curves will be more incomplete.
It has to be pointed out that the change rule is presented based on the given probe geometry with fixed coil center distance
dc and fixed coil center lift-off
lc as shown in
Figure 3a. When the flare angle changes, the distance from the bottom of the coil to the specimen also changes. For a single tilted coil, the eddy current is stronger when the coil is closer to the specimen [
28]. For the presented probe with four tilted coils, the concerned rotating field is the superposition of the field generated by four coils. To accommodate adjacent coils, an adequate space between each pair of oppositely tilted coils is required. In this case, the eddy current density is related not only to the distance from the coil to the specimen but also to the space between opposite coils. Since eddy current produced by a single tilted coil tends to focus underneath the near side of the coil [
28], two focusing coils with a smaller space between their inner sides will result in a stronger net eddy current field. However, when the flare angle changes, the distance from the coil on the near side of the specimen (denoted by
l0) and the space between opposite coils (denoted by
d0) will change in opposite directions, as described by three typical cases in
Figure 10 when the flare angle is 0°, 120° and 180°, respectively. When the flare angle is 0°,
l0 is the minimum value and
d0 is the maximum. The opposite situation occurs when the flare angle changes to 180°. Then, the intermediate case as shown in
Figure 10b has the potential to generate stronger eddy current underneath the probe center. When the flare angle exceeds 180°, the inner side of the coils comes up while the outside goes down, as is the case with a flare angle of 240°, shown in
Figure 10d. Taking the case with a 180° flare angle as the reference, due to the overall higher values of
l0, eddy current generated by the probes with flare angles greater than 180° will be weaker than eddy current generated by the probes with rival flare angles, for instance, the case with a 240° flare angle versus the case with a 120° flare angle. This explains the curves in
Figure 9 with the peak values emerging before 180° and with steeper left rising parts than right descending parts.
4. Experimental Verification
Experimental study was also conducted to verify the probe performance.
Figure 11 shows the composition, prototype probe and the specimen of the experiment system. The probe has a flare angle of 150 degrees and its coil parameters are the same as those given in
Table 2. The specimen is a Q235 steel plate with four prefabricated grooves, with the same length and width but different depths, as shown in
Figure 11c. Due to the machining error, the actual depths of the grooves are 2.24 mm, 4.80 mm, 5.94 mm and 8.64 mm, measured by a Vernier caliper with an accuracy of 0.02 mm, which deviate respectively from their design values of 2 mm, 4 mm, 6 mm and 8 mm. A sine wave of 100 Hz frequency and 0° initial phase is generated by a function generator and output through two channels. One is directly amplified by the power amplifier A, while the other is sent to the phase shifter module to form a 90-degree phase shift before being amplified by the power amplifier B. After that, the two amplified sinusoidal currents are fed into the sub-probe A and B respectively. The two power amplifiers are of the same type and with identical performances.
Considering the low excitation frequency, a GMR AA002-02 sensor from NVE Corporation [
29] was used as the detection sensor. It has a linear range from 1.5 to 10.5 Oe (1 Oe = 0.1 mT in air) and presents the highest sensitivity of [3.0, 4.2] mV/V-Oe in NVE AA00x series [
29]. The GMR sensor was mounted in the probe center surrounded by the excitation coils with its sensitivity axis along the normal direction of the specimen surface. This compact structure helps decrease the probe design lift-off. A small circuit board for mounting the GMR was used. Since the GMR sensor operates in a unipolar mode, a small permanent magnet is placed near the sensor to bias the sensor in the middle of the linear range [
30]. Taking all factors into account, the total lift-off of the GMR sensor is around 2.5 mm. The GMR output is conditioned by a custom-made preamplifier and then digitalized by a NI data acquisition system linked with a personal computer for subsequent signal processing. The magnitude and phase angle are extracted by using a digital lock-in amplifier program coded in LabVIEW.
The excitation current after the power amplifier has an amplitude of 1 A, which matches the corresponding current density used in the simulation. The phase shifter module is homemade and can generate 0 to 360 degrees shift at the frequency range of a few hertz up to 1 kHz. It employs a two-stage shifter structure and provides two multi-turn trimpots with high resolution for manual adjustment to obtain a high accuracy of phase shift. The custom-made preamplifier functions as a signal filter and an amplifier whose cut-off frequency and gain factor can be configured as needed. During the experiment, the cut-off frequency and gain factor were set as 1 kHz and 10, respectively.
Figure 12 presents the schematic circuit diagrams of the phase shifter and the preamplifier. Note that R18 and R19 in
Figure 12a represent the trimpots, while S1 to S4 in
Figure 12b denote the DIP switches for setting cut-off frequency and gain factor.
First, the probe was placed with the GMR sensor right above the width edge of the defect. In order to simulate the variation of defect orientation, the probe was in situ rotated step by step on the specimen surface.
Figure 13a and b show the experiment signal magnitude and phase angle versus the defect orientation for four different depth defects. For comparison, simulation results are also shown. The conversion of simulation data from mT to V are based on the following relationship—the response curve of GMR AA002-02 at room temperature is tested for 100 Hz applied magnetic field and its sensitivity in the linear range is determined as 3.2 mV/V-Oe. The gain factor of the preamplifier is 10 and thus, the equivalent voltage is 0.32 V per 1 mT. It can be seen that the experiment results are basically in agreement with the simulation results—as the defect orientation varies, the magnitude holds stable while the phase angle exhibits a linear variation. The observed slight deviations between experiment and simulation results are tolerable considering various kinds of influence factors including the linear and homogeneous material assumption made in the simulation, the GMR sensor calibration error and the phase shifting caused by the analog electronic devices involved in experiments. Furthermore, according to
Figure 13a, as the defect depth increases, the signal magnitude increases as well, while the signal phase shows no significant change. All the above findings indicate that the signal magnitude can be used to evaluate the defect depth, regardless of the defect orientation.
A line scan using the presented probe is also carried out on the specimen. The scan path is along the lengthwise center line of four defects. To avoid the edge effect, the beginning and ending points are 50 mm away from the edges, so the total scanning distance is 300 mm. All scan points are operated manually with a step of 2 mm.
Figure 13c shows the scanning results. From left to right, the four amplitude-rising signals correspond well with the four defects with depth of 2.24, 4.80, 5.94 and 8.64 mm, respectively. Good SNR was also observed. Small signal fluctuations are probably due to the inaccurate manual operation and interferences caused by the intrinsic non-uniform magnetic permeability of the carbon steel.
5. Conclusions
A novel rotating focused field eddy-current probe for the detection of arbitrary orientation defects is presented. The magnetic field focusing method of the figure-8-shaped coil is integrated into the generation of the rotating excitation field. The hybrid excitation strategy gives the designed probe advantages in both structure and function. Instead of a stacked structure, the probe employs a symmetric and equal lift-off configuration for all excitation coils and the receiving GMR sensor, thus making a small design lift-off and avoiding the additional excitation adjustment encountered in existing rotating field probes. By optimizing the probe’s flare angle, the probe shows the capability of focusing eddy current at the defect area while weakening it at the background area, which contributes to a defect signal with a high signal-to-noise ratio. Feasibility of the probe to detect arbitrarily oriented defects was studied using simulation and validated by experiment tests on a carbon steel plate. A line scan over defects with different depth was also performed to demonstrate the probe’s focusing effect.
It should be noted that, in this paper, only the effect of the flare angle was preliminarily analyzed for the fixed-coil-geometry probe. Other parameters including the coil shape, coil dimensions, the distance between each pair’s coils are factors related to the probe’s performance. Multi-parameter optimization is pending in future work. In addition, more work will be done to further validate and refine the probe performances, including developing a C-scan presentation to image the defect and performing more experiments to study the detection accuracy for fine cracks and deep buried defects.