Next Article in Journal
YOLO-LWNet: A Lightweight Road Damage Object Detection Network for Mobile Terminal Devices
Previous Article in Journal
Validity and Reliability of a Wearable Goniometer Sensor Controlled by a Mobile Application for Measuring Knee Flexion/Extension Angle during the Gait Cycle
Previous Article in Special Issue
Real-Time Laser Interference Detection of Mechanical Targets Using a 4R Manipulator
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of the Properties of Eddy Current Sensors Based on Their Equivalent Parameters

Faculty of Automatic Control, Electronics and Computer Science, Silesian University of Technology, Akademicka 16, 44-100 Gliwice, Poland
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(6), 3267; https://doi.org/10.3390/s23063267
Submission received: 29 December 2022 / Revised: 2 February 2023 / Accepted: 25 February 2023 / Published: 20 March 2023

Abstract

:
This paper presents a practical way of using the method of evaluating the metrological properties of eddy current sensors. The idea of the proposed approach consists of employing a mathematical model of an ideal filamentary coil to determine equivalent parameters of the sensor and sensitivity coefficients of tested physical quantities. These parameters were determined on the basis of the measured value of the real sensor’s impedance. The measurements were carried out with an air-core sensor and an I-core sensor while they were positioned at different distances from the surface of tested copper and bronze plates. An analysis of the influence of the coil’s position in relation to the I core on the equivalent parameters was also carried out, and the interpretation of the results obtained for various sensor configurations was presented in a graphical form. When equivalent parameters and sensitivity coefficients of examined physical quantities are known, it is possible to compare even very different sensors with the employment of one measure. The proposed approach makes it possible to make a significant simplification of the mechanisms of calibration of conductometers and defectoscopes, computer simulation of eddy current tests, creating the scale of a measuring device, and designing sensors.

1. Introduction

Eddy current testing is commonly used in many industries for the nondestructive evaluation of metal products. The tests are carried out using a single sensor in the form of a coil [1,2,3,4,5,6] or a system consisting of several coils constituting one sensor [7,8,9,10,11,12]. The basic idea of such an inspection is to utilise the phenomenon of electromagnetic induction. The alternating current flowing in the coil creates a magnetic field. If there is a conductive element within this field, eddy currents are induced in the conductor, creating a secondary magnetic field. According to Lentz’s rule, the secondary magnetic field is directed opposite to the primary field and thus causes a change in the impedance of the coil. This property makes it possible to detect defects in conductive objects [13,14,15,16,17,18,19] and determine their geometric dimensions [20,21,22,23,24,25], magnetic permeability [26,27,28], or electrical conductivity [29,30,31]. The changes in the impedance components of the sensor, and therefore the sensitivity of the measurement, strongly depends on the frequency of the current that excites the eddy currents and on the geometry of the sensor. In order to select the eddy current sensor construction that is optimal for a specific application, it is necessary to carry out many calculations and tests at the stage of its design. It is also important to properly select the frequency of the current exciting the coil since it affects not only the eddy currents penetration depth but also the metrological properties of the measuring instrument, including its sensitivity. Another important solution that needs to be implemented in the measuring instrument is the mechanism eliminating the influence of the lift-off effect [32,33] and the inaccuracy of the sensor position in relation to the tested element.
In all these applications, it is most convenient to employ the simplest possible coil model, which will significantly reduce the computational complexity and make it easier for the designer to select the optimal geometry of the sensor and test parameters. The mathematical model for a single-turn coil positioned over a conductive half-space was developed in [34,35]. The authors found in their research that, in many respects, it is much more convenient to use the model of an ideal filamentary coil (Figure 1). Such a coil contains n infinitely thin turns concentrated in a circle of radius r0, positioned at a distance h0 from the tested surface. Any coil used in eddy current testing can be associated with the corresponding equivalent parameters r0, h0 of a filamentary coil while the number of turns is maintained the same [36]. For this purpose, it is enough to determine such values of equivalent radius r0 and equivalent height h0 so that the impedance of the filamentary coil and the measured impedance of any real coil are equal. A mathematical model containing analytical expressions for the impedance of a filamentary coil placed over a two-layer conductive half-space [37] was derived using the truncated region eigenfunction expansion (TREE) method [38]. The algorithm in Matlab enabling the effective determination of the parameters r0, h0 was presented in [39] and subsequently used to detect defects with a filamentary coil [40]. The mathematical model for the system consisting of two coils was developed in [41,42,43,44,45,46]. In [47], it was proved that the change in the frequency and the change of the conductivity of the plate do not affect the change of the equivalent parameters. This convenient property was used by the authors while designing a sensor for a station for detecting cracks in clutches made with the method of powder metallurgy. The complicated shape of the tested elements necessitated the use of I-core sensors [48,49,50,51,52,53] with a small diameter. This core acted as a magnetic circuit and made it possible to move the coil away from the tested element. Such a thin core may touch the surface located in deep slots of the tested object.
In this work, the method of evaluating the properties of eddy current sensors with the employment of equivalent parameters r0, h0, described in [47], was used and extended. Section 2 presents the already-known analytical final formulas for the change in the impedance of the filamentary coil. Then, for the first time, a sensitivity analysis using a filamentary coil was performed for I-core sensors, which made it possible to determine the following:
the sensitivity of the resistance component to lift-off,
the sensitivity of the inductance component to lift-off,
the sensitivity of the resistance component to the electrical conductivity of the material under test,
the sensitivity of the inductance component to the electrical conductivity of the material under test,
the sensitivity of the resistance component to the magnetic permeability of the material under test,
the sensitivity of the inductance component to the magnetic permeability of the material under test.
The characteristics of the changes in sensitivity coefficients have been shown in the graphs. An algorithm for determining the sensitivity coefficients of the examined parameter in relation to changes in the resistance and reactance of the filamentary coil has also been presented. The Results section describes the measurements and calculations in detail. Equivalent parameters were determined for the air-core sensor and the I-core sensor. The analysis of what the changing of the distance of the sensor from the tested surface (lift-off) presented in [47] brings about was significantly extended. The novelty here is the detailed study of the influence of the coil position in relation to the core on the values of equivalent parameters, which were performed for various configurations. Subsequently, the obtained results were presented for 16 configurations in a graphical form, making it possible to observe the relationship between the geometry and position of the sensor and the value of equivalent parameters. In the Discussion section, the six most important conclusions for sensor constructors were formulated and discussed, thus enabling the practical application of the developed approach. A comparative analysis that describes how to correctly interpret the obtained values of equivalent parameters was also carried out.

2. Materials and Methods

The filamentary coil with n turns concentrated in a circle of radius r0 at a distance h0 from the conductive surface is shown in Figure 2. It was assumed that the thickness of the plate is greater than the depth of the eddy current penetration, which makes it possible to replace the plate in the mathematical model by means of half-space. The change in the coil impedance resulting from bringing it close to a material with conductivity ϭ and magnetic permeability μr was formulated as follows [36]:
Δ Z = j ω π   r 0 μ 0 μ r n 2 β 0 c ( λ ) e α β λ J 1 2 ( β λ ) d λ ,
where:
α = 2 h 0 r 0 ,
β = r 0 ω μ 0 μ r σ ,
c ( λ ) = λ λ 2 + j λ + λ 2 + j ,
ω is the angular frequency, μ0 is the permeability of free space, and J1(x) is the Bessel function of the first kind of order 1.
When the impedance of the coil positioned at some distance from the conductor is denoted as Z0, and the impedance of the coil close to the conductor as Z, the change in impedance ΔZ = ZZ0 can be represented as:
Δ Z = R - R 0 + j ω ( L L 0 ) = R - R 0 - j ω ( L 0 L ) .
The changes in the impedance components were written in such a way so that both the change in resistance ΔR = r = RR0 and the change in inductance ΔL = l = L0L were positive for non-ferromagnetic material. By determining the real part and the imaginary part of the impedance, the change in resistance r and the change in inductance l were subsequently calculated [36].
r = R R 0 = n 2 ω   π μ 0 μ r r 0   φ ( α , β ) ,
l = L 0 L = n 2 π μ 0   μ r r 0   χ ( α , β ) ,
where:
φ ( α , β ) = Re { j β 0 c ( λ ) e α β λ J 1 2 ( β λ ) d λ } ,
χ ( α , β ) = Im { j β 0 c ( λ ) e α β λ J 1 2 ( β λ ) d λ } .
Using the exact differential equations, the Functions (8) and (9) can be written using the parameters α and β, defined in Equations (2) and (3).
Δ φ ( α , β ) = φ α Δ α + φ β Δ β ,
Δ χ ( α , β ) = χ α Δ α + χ β Δ β ,
where:
Δ α = 2 r 0 Δ h 0 ,
Δ β = β 2 σ Δ σ .
According to Equations (6) and (7), the changes in the components of the coil depend on the changes in the functions ϕ (α,β) and χ (α,β). The changes in resistance and inductance of the coil can be written in the following form:
Δ r = n 2 π β 2 r 0 σ Δ φ ,
Δ   l = n 2 π μ 0 μ r r 0 Δ χ .
Expressions (10)–(15) can be used in the sensitivity analysis to determine the influence of a given parameter on the change in the resistance and the inductance of a filamentary coil. The differentiation of Equations (6) and (7) brought about the determination of the sensitivity coefficients that define:
the change in the coil resistance resulting from changing the distance of its location from the tested surface (16),
the change in the coil resistance caused by changing the conductivity of the tested element (17),
the change in the coil inductance resulting from changing the distance of its location from the tested surface (18),
the change in the coil inductance caused by changing the conductivity of the tested element (19).
Δ r h 0 = Δ r Δ h 0 = 2 n 2 π β 2 r 0 2 σ φ α ,
Δ r σ = Δ r Δ σ = n 2 π β 3 2 r 0 σ 2 φ β ,
Δ l h 0 = Δ l Δ h 0 = 2 n 2 π μ 0 χ α ,
Δ l σ = Δ l Δ σ = n 2 π μ 0 r 0 β 2 σ χ β .
To denote the sensitivity, three-letter symbols have been adopted. The symbol Δr means that the sensitivity is defined as the ratio of the change in resistance r to the change in a parameter of the tested element. Analogously, in the case of Δl, it is the ratio of the change in inductance l to the change in a given parameter. The third letter of the sensitivity designation is the symbol of a parameter of the tested element.
Determining the sensitivity coefficients for changes in magnetic permeability is more complicated because μr occurs not only in Equation (3) but also in Equations (6) and (7). At first, μr was obtained from Equation (3) and substituted into Equations (6) and (7). In the next step, the differentiation of Equations (6) and (7) brought about the determination of the expressions for the change in resistance and inductance of the coil.
Δ r = [ 2 π n 2 ω μ 0 μ r σ φ ( α , β ) + π n 2 r 0   ω μ 0 μ r φ ( α , β ) β ]   Δ β ,
Δ l = [ 2   π n 2 μ 0 μ r ω   σ χ ( α , β ) + π n 2 r 0   μ 0 μ r χ ( α , β ) β ]   Δ β .
Then, the dependence Δβ = f(μr) was determined from Equation (3).
Δ β = r 0 2 ω   σ μ 0 μ r Δ μ r .
In the final step, Equation (22) was substituted into Equations (20) and (21) and the sensitivity coefficients for magnetic permeability were obtained.
Δ   r μ r = Δ   r Δ μ r = [ 2 π n 2 ω μ 0 μ r σ   φ ( α , β ) + π n 2 r 0   ω μ 0 μ r φ ( α , β ) β ]   r 0 2   ω σ μ 0 μ r ,
Δ l μ r = Δ   l Δ μ r = [ 2 π n 2 μ 0 μ r ω σ χ ( α , β ) + π n 2 r 0 μ 0 μ r χ ( α , β ) β ]   r 0 2 ω σ μ 0 μ r .
Expressions (23) and (24) can be written in a compact form:
Δ   r μ r = π n 2 r 0 ω μ 0 φ ( α , β ) + π n 2 r 0 2 ω 3 μ 0 3 μ r σ 2 φ ( α , β ) β ,
Δ l μ r = π n 2 r 0 μ 0 χ ( α , β ) + π n 2 r 0 2 ω μ 0 3 μ r σ 2 χ ( α , β ) β .
In order to present the characteristic of the changes in sensitivity coefficients, they were presented as a function of parameter β for several values of α (Figure 3). This choice results from the fact that the conductivity of the tested material in the expression for the change in the coil impedance (1) occurs only in the form of parameter β, and the distance h0 only in the form of parameter α.
Parameter βm is the maximum of function ϕ (α,β), and it determines the frequency for which the resistance of the sensor does not change along with the change in the conductivity. As distance h0 increases, all sensitivity coefficients (10)–(13) decrease, and the sensitivity of the instrument also decreases. The sensitivity coefficients Δ and Δ depend greatly and non-monotonically on the conductivity, frequency, and radius r0. This means that only the optimal selection of these parameters will ensure the best sensitivity of the sensor. Therefore, at the stage of designing a conductometer, flaw detector, and other eddy current devices, it is crucial to carefully select the geometry of the sensor and the frequency of the current. According to the authors, it is convenient to employ the method developed for this purpose, which consists in carrying out what follows:
(1)
In the first step, the components of impedance Z0 of the sensor positioned away from the conductive material should be measured. During the measurement, it is best to use the same frequency value as in the planned tests.
(2)
The second measurement is made for the same sensor after bringing it close to a conductive material of known conductivity σ and known magnetic permeability µr. The impedance Z should be measured at the same frequency as in step 1.
(3)
On the basis of the measured impedance values, the changes in resistance ΔR = RR0 and the changes in inductance ΔL = L0L caused by bringing the sensor close to the conductor are determined. In the case of the computer simulation of the tests or at the stage of designing the sensor, the changes in impedance components may be calculated with, for example, the finite element method or analytical expressions.
(4)
After obtaining the changes in the impedance components of the sensor, with the employment of the mathematical model of the filamentary coil, substitute parameters r0, h0 are calculated. At this stage, it is also possible to determine coefficients of the sensitivity of the analysed parameter to the change in the resistance and the change in the filamentary coil’s reactance.
When the values of r0, h0 and appropriate sensitivity coefficients are known, it is possible to precisely determine the properties of the sensor. The use of the same method of evaluation based on equivalent parameters for different sensors with different geometries makes it possible to compare them according to one criterion, which greatly simplifies the analysis process.

3. Results

In the measurements and calculations, a sensor in the form of a cylinder-shaped coil containing 700 turns was applied (Figure 4). The outer diameter of the coil was 25 mm, whereas the height was 14.8 mm. In the centre of the coil, there was a hole with a diameter of 8.1 mm, wherein an I core made of material F1001 with magnetic permeability µ = 3500 was placed. The employed core had a diameter of 8 mm and a height of 29.3 mm. The sensor impedance was measured using the Agilent E4980A precision LCR meter with an accuracy of ±0.1% for frequencies f = 10 kHz and f = 100 kHz. Each measurement was performed three times, and the arithmetic mean of the obtained values was calculated. In this way, the effect of random measurement error was reduced. Such an error is also influenced by the precision of the core installation and the positioning of the sensor in relation to the tested element. The influence of these parameters on the measurement result may be noticeable primarily for small values of the lift-off distance. In order to minimise this influence, washers were utilised, which ensured: the stabilisation of individual elements of the sensor and a position parallel in relation to the tested surface. The geometric dimensions of the coil and the core were measured using a micrometre screw gauge with an accuracy of ±0.05 mm.
In the first stage of the experiment, the impedance components of the air-core sensor positioned at different distances from the conductive surface, which were copper (σ = 58.42 MS/m) and bronze (σ = 10.47 MS/m) plates, were measured. The electrical conductivity of the plates was measured using a Foerster Sigmatest 2.069 instrument with an absolute accuracy of ±0.5% of the measured value. The thickness of the plates was greater than the depth of the eddy currents penetration. At the same time, the surface of each plate was large enough so that no edge effect was observed. In the next stage of the experiment, a core was placed in the coil. By changing the position of the coil in relation to the core and changing the distance of the sensor from the conductive surface, the impedance components of the sensor were measured. The measurements were carried out for 16 test configurations utilising the air-core sensor (Table 1) and for 32 configurations using the I-core sensor (Table 2). The lift-off distance of the sensor from the tested surface is marked as h. In the case of the sensor without a core, it is the distance from the lower edge of the coil to the surface of the tested element. If the sensor has a core, it is the distance from the core to the surface. In such a case, the distance from the bottom edge of the coil to the bottom edge of the core is hc. For each measurement, equivalent parameters r0, h0 of the filamentary coil were determined using the algorithm developed in [39]. The value of the filamentary coil impedance determined with the employment of the equivalent parameters was different from the measured value of the sensor impedance by less than 0.01%. The integration in Formula (1) was performed using Gauss–Kronrod quadrature implemented in Matlab. The absolute error tolerance of the integration in the numerical procedure was 1 × 10−8. The changes in frequency f and conductivity σ caused very small changes in the values of r0, h0—much smaller than the error of the method for measuring the impedance components. For this reason, Table 1 and Table 2 also show the average value of the equivalent parameters, which turned out to be very useful in the final assessment of the properties of the sensors.
The interpretation of the obtained results is shown in Figure 5, Figure 6, Figure 7 and Figure 8. Next to the dimensioned sketch of the sensor and the tested element, equivalent parameters r0, h0 are presented. Due to this, it is easy to compare the influence of positioning the coil in relation to the core and the influence of moving the sensor away from the tested surface on the value of the equivalent parameters and, thus, on the metrological properties of the sensors.
Table 1. Equivalent parameters of the air-core sensor.
Table 1. Equivalent parameters of the air-core sensor.
Average Value
h
[mm]
hc
[mm]
f
[kHz]
σ
[MS/m]
ΔR
[Ω]
ΔX
[Ω]
ΔL
[mH]
r0
[mm]
h0
[mm]
r0
[mm]
h0
[mm]
Figure
Number
0-1058.425.8159.880.958.306.138.306.13Figure 5a and Figure 8a
-10058.4220.07643.61.028.296.12
-1010.4711.4451.250.828.316.14
-10010.4744.87615.30.988.306.12
4-1058.421.6421.920.358.3710.128.3710.13Figure 5b and Figure 7a
-10058.425.55231.50.378.3610.12
-1010.473.3619.510.328.3710.14
-10010.4712.56223.70.368.3710.12
8-1058.420.589.760.168.4214.248.4214.24Figure 5c
-10058.421.94101.90.168.4214.23
-1010.471.238.920.148.4314.25
-10010.474.4499.20.168.4214.23
16-1058.420.122.910.058.5022.458.5022.45Figure 5d
-10058.420.4030.030.058.5022.45
-1010.470.272.740.048.5022.46
-10010.470.9329.470.058.5022.45
Table 2. Equivalent parameters of the I-core sensor.
Table 2. Equivalent parameters of the I-core sensor.
Average Value
h
[mm]
hc
[mm]
f
[kHz]
σ
[MS/m]
ΔR
[Ω]
ΔX
[Ω]
ΔL
[mH]
r0
[mm]
h0
[mm]
r0
[mm]
h0
[mm]
Figure
Number
001058.4225.74362.65.7715.35.0815.305.12Figure 6a and Figure 8b
10058.4284.6938126.0715.55.27
1010.4755.04325.15.1815.04.99
10010.47195.336955.8815.45.14
041058.4217.452734.3415.66.7115.686.76Figure 7b and Figure 8c
10058.4257.1528544.5515.86.91
1010.4737.68247.63.9415.66.64
10010.47123.227754.4215.76.78
081058.4211.571963.1215.68.515.618.46Figure 8d
10058.4237.820433.2515.88.71
1010.4725.16179.32.8615.68.23
10010.4787.5619913.1715.58.39
401058.4210.58155.92.4813.37.5513.357.56Figure 6b
10058.4235.4616382.6113.47.56
1010.4722.1140.52.2413.47.56
10010.4780.7615882.5313.47.56
441058.427.79127.92.0414.09.2814.059.29Figure 7c
10058.4225.9713362.1314.09.28
1010.4716.47116.51.8614.19.31
10010.4759.3612992.0714.09.28
801058.424.1173.791.1713.611.8913.6111.91Figure 6c
10058.4213.68768.11.2213.611.88
1010.478.7467.841.0813.611.98
10010.4731.30749.11.1913.611.87
881058.422.4652.150.8314.615.5514.6315.53Figure 7d
10058.428.13539.40.8614.615.51
1010.475.3348.610.7714.715.54
10010.4718.7528.10.8414.615.53
1601058.420.9423.650.3814.421.1114.3721.11Figure 6d
10058.423.09243.30.3914.421.08
1010.472.0722.30.3614.421.15
10010.477.15238.90.3814.421.10

4. Discussion

The essence of the proposed approach lies in correctly interpreting the values of parameters r0, h0, which unambiguously characterise the metrological properties of eddy current sensors. Based on long-running studies and the results presented in this paper, the authors formulated and then discussed a number of conclusions that enable the effective and efficient application of the developed approach.
(1)
The lower the value of h0, the greater the sensitivity of the sensor.
The intensity of the eddy currents induced in the tested material has a decisive influence on the sensitivity of the sensor. The smaller the distance h0, the closer the filamentary coil is to the conductor, so the intensity of the eddy currents increases (Figure 3). It is well known that the core concentrates the magnetic flux around itself, thus increasing the sensitivity of the sensor. How large this increase is can be assessed based on the value of parameter h0.
(2)
The higher the value of parameter r0, the greater the changes in the impedance components caused by the change in the measured quantities.
Both in the process of designing the sensor and while selecting the optimal test parameters, it is important to obtain the highest value of r0 and, thus, the highest possible sensitivity of the sensor.
(3)
The employment of the core brings “the coil magnetically closer to the tested material”.
The use of the core results in a decrease in the value of h0, which corresponds to bringing the coil in contact with the conductor. Parameter h0 may be treated as a measure of the efficiency of the core performance.
(4)
The higher the value of parameter r0, the larger the area of the tested element where eddy currents that change the components of the coil impedance are induced.
Parameter r0 depends primarily on the inner radius and the outer radius of the real coil. Along with an increase in the value of these radii, the coverage of eddy currents induction in the conductor also increases. This property should be utilised appropriately according to the application of the sensor. In the case of testing large surfaces, such as plates or metal sheets, designers aim to obtain the most extensive inspection coverage. However, there are applications, such as rod testing, coin sorting or crack detection in small-diameter cylinders, where a short inspection coverage is preferable. This requirement results from the fact that various factors disturbing the measurement occur in close proximity or that there is the likelihood of an edge effect.
(5)
The higher the value of hc, the lower the sensitivity of the sensor.
When height hc increases due to moving the coil “up” along the core, the value of equivalent parameters increases. The sensitivity of the sensor decreases because a small increase in r0 cannot compensate for a very large increase in h0.
(6)
The number of turns of the coil does not affect the value of equivalent parameters.
Changing the number of coil turns does not change the value of r0, h0. The parameter specifying the number of turns appears as n2 both in the formula for the impedance of the real coil and in the formula for the filamentary coil (1). This property constitutes an important advantage over the analysis with a single-turn coil. Only in the case of a filamentary coil is it possible that the number of turns does not impact the results of the comparison and evaluation of the properties of eddy current sensors.
The above conclusions were used to perform an exemplary comparative analysis. The equivalent parameters determined for the air-core sensor placed on the tested surface are r0 = 8.3 mm and h0 = 6.13 mm. After having inserted the core, the value of parameter h0 decreases (h0 = 5.12 mm), thus increasing the sensitivity of the sensor. At the same time, parameter r0 increases from 8.3 mm to 15.3 mm, which causes an additional increase in sensitivity due to the increase in the magnetic field range. Similarly, it is possible to analyse the effect of moving the sensor away from the tested surface. The air-core sensor that is brought into contact with the tested surface has an equivalent height of h0 = 6.13 mm. After moving the sensor 4 mm away from the tested surface, h0 = 10.13 mm, i.e., Δh = 10.13 − 6.13 = 4 mm. If the sensor is 8 mm away from the tested surface, h0 = 14.2 mm, i.e., Δh = 14.2 − 6.13 = 8.07 mm. Moving the air-core sensor away from the surface of the tested element by a particular value increases h0 by the same value. In the case of the I-core sensor, the changes in parameter h0 are slightly different. Such a sensor, when brought into contact with the tested surface, takes the value of parameter h0 = 5.12 mm. When the sensor is moved 4 mm away from the tested surface, h0 = 7.56 mm, i.e., Δh = 7.56 − 5.12 = 2.44 mm. Moving the sensor away from the tested surface by 8 mm brings about an increase of the value of parameter h0 = 11.91 mm, i.e., Δh = 11.91 − 5.12 = 6.79 mm. The beneficial effect of the core is apparent. For small distances, moving the sensor away increases the equivalent height h0, however, by a value that is smaller than the actual distance. Thus, the sensitivity of the I-core sensor decreases less significantly than that of the air-core sensor.

5. Conclusions

This paper presents a practical way of using the method of evaluating the metrological properties of eddy current sensors with the use of equivalent parameters of a filamentary coil. The presented sensitivity analysis makes it possible to obtain expressions that define the influence of a selected parameter of the sensor-tested element system on changes in the coil impedance components. The developed geometric interpretation of changes in sensitivity coefficients allows a better understanding of the characteristic of changes in the analysed sensitivity coefficients. It is the first time that the question of how the position of the coil in relation to the core affects the equivalent parameters has been examined. In the measurements, a sensor with a removable I-core was used. Due to this, it was possible to observe the influence of the core on the metrological parameters of the sensor. The measurements were carried out for different frequency values using thick copper and bronze plates. The sensor was positioned at different distances from the tested surface, while the position of the coil with respect to the core was also changed. The measured values of the impedance changes were employed to determine the equivalent parameters. The results of the calculations confirmed that parameters r0, h0 may be successfully used to evaluate the properties of eddy current sensors.
The presented method for describing the sensors using two numbers that are equivalent parameters of r0, h0 is easy to perform and provides a number of benefits. The usefulness of the proposed concept consists in enabling the comparison of very different sensors with the employment of the same measure. It is enough to determine the sensitivity coefficients of the examined physical quantities with the application of the universal filamentary coil model. Then, equivalent parameters, which are the dimensions of the ideal filamentary coil, are determined for the real sensors. On the basis of the once-determined values of r0, h0 and sensitivity coefficients, the metrological properties of the sensors can be assessed. The proposed method can be utilised for many applications, such as the designing of sensors, selection of test parameters, computer simulations, interpretation of the obtained results, etc.
In this article, the properties of I-core sensors were analysed due to the fact that they were used in the designed device. In further work, sensors with different geometry will also be employed to detect various types of defects. In addition, tests of conductors of other shapes, including roughness and porosity of the tested surface, are being planned.

Author Contributions

Conceptualization, L.D.; methodology, L.D.; software, G.T.; validation, L.D. and G.T.; formal analysis, L.D. and G.T.; investigation, L.D. and G.T.; resources, L.D. and G.T.; data curation, L.D. and G.T.; writing—original draft preparation, L.D. and G.T.; writing—review and editing, L.D. and G.T.; visualization, L.D.; project administration, G.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tytko, G. An eddy current model of pot-cored coil for testing multilayer conductors with a hole. Bull. Pol. Ac. Tech. 2020, 68, 1311–1317. [Google Scholar]
  2. Faria, A.; Marques, L.; Ferreira, C.; Alves, F.; Cabral, J. A Fast and Precise Tool for Multi-Layer Planar Coil Self-Inductance Calculation. Sensors 2021, 21, 4864. [Google Scholar] [CrossRef] [PubMed]
  3. Luo, Y.; Theodoulidis, T.; Zhou, X.; Kyrgiazoglou, A. Calculation of AC resistance of single-layer coils using boundary-element method. IET Electr. Power Appl. 2021, 15, 1–12. [Google Scholar] [CrossRef]
  4. Wang, S.-C.; Xie, B.-R.; Huang, S.-M. Design and Analysis of Small Size Eddy Current Displacement Sensor. Sensors 2022, 22, 7444. [Google Scholar] [CrossRef]
  5. Tytko, G. Measurement of multilayered conductive discs using eddy current method. Measurement 2022, 204, 112053. [Google Scholar] [CrossRef]
  6. Liu, Z.; Li, Y.; Ren, S.; Ren, Y.; Abidin, I.M.Z.; Chen, Z. Pulse-Modulation Eddy Current Evaluation of Interlaminar Corrosion in Stratified Conductors: Semi-Analytical Modeling and Experiments. Sensors 2022, 22, 3458. [Google Scholar] [CrossRef]
  7. Dehui, W.; Qisheng, S.; Xiaohong, W.; Fan, Y. Analytical model of mutual coupling between rectangular spiral coils with lateral misalignment for wireless power applications. IET Power Electron. 2018, 11, 781–786. [Google Scholar] [CrossRef]
  8. Betta, G.; Ferrigno, L.; Laracca, M.; Rasile, A.; Sangiovanni, S. A novel TMR based triaxial eddy current test probe for any orientation crack detection. Measurement 2021, 181, 109617. [Google Scholar] [CrossRef]
  9. Lu, M.; Meng, X.; Huang, R.; Peyton, A.; Yin, W. Analysis of Tilt Effect on Notch Depth Profiling Using Thin-Skin Regime of Driver-Pickup Eddy-Current Sensor. Sensors 2021, 21, 5536. [Google Scholar] [CrossRef]
  10. Lu, M.; Chen, L.; Meng, X.; Huang, R.; Peyton, A.; Yin, W. Thickness Measurement of Metallic Film Based on a High-Frequency Feature of Triple-Coil Electromagnetic Eddy Current Sensor. IEEE Trans. Instrum. Meas. 2021, 70, 1–9. [Google Scholar] [CrossRef]
  11. Mirzaei, M.; Ripka, P.; Grim, V.; Jelinek, L.; Svatos, J. Thickness evaluation of hollow nonmagnetic cylinders utilizing a motional eddy current. Measurement 2022, 189, 110463. [Google Scholar] [CrossRef]
  12. Wang, D.; Zhang, Z.; Yin, W.; Chen, H.; Ma, H.; Zhou, G.; Zhang, Y. Linear Characteristics of the Differences in Phase Tangents of Triple-Coil Electromagnetic Sensors and Their Application in Nonmagnetic Metal Classification. Sensors 2022, 22, 7511. [Google Scholar] [CrossRef] [PubMed]
  13. Ali, K.B.; Abdalla, A.N.; Rifai, D.; Faraj, M.A. Review on system development in eddy current testing and technique for defect classification and characterization. IET Circuits Devices Syst. 2017, 11, 330–343. [Google Scholar] [CrossRef] [Green Version]
  14. Yu, Y.; Gao, K.; Liu, B.; Li, L. Semi-analytical method for characterization slit defects in conducting metal by eddy current nondestructive technique. Sens. Act. A Phys. 2020, 301, 111739. [Google Scholar] [CrossRef]
  15. Xu, P.; Zhu, C.L.; Zeng, H.M.; Wang, P. Rail crack detection and evaluation at high speed based on differential ECT system. Measurement 2020, 166, 108152. [Google Scholar] [CrossRef]
  16. Ha, N.; Lee, H.-S.; Lee, S. Development of a Wireless Corrosion Detection System for Steel-Framed Structures Using Pulsed Eddy Currents. Sensors 2021, 21, 8199. [Google Scholar] [CrossRef]
  17. Xie, F.; Wu, X.; Zhang, H.; Zhang, X. Research on pulsed eddy current automatic detection technology for weld defects of pressure vessel cylinder. Measurement 2021, 176, 109081. [Google Scholar] [CrossRef]
  18. Fan, X.; Chen, T.; He, Y.; Du, J. Influence of spatial winding distribution of flexible eddy current sensor on quantitative monitoring of subsurface cracks. Measurement 2021, 178, 109382. [Google Scholar] [CrossRef]
  19. Farag, H.E.; Toyserkani, E.; Khamesee, M.B. Non-Destructive Testing Using Eddy Current Sensors for Defect Detection in Additively Manufactured Titanium and Stainless-Steel Parts. Sensors 2022, 22, 5440. [Google Scholar] [CrossRef]
  20. Wang, H.; Li, W.; Feng, Z. Noncontact thickness measurement of metal films using eddy-current sensors immune to distance variation. IEEE Trans. Instrum. Meas. 2015, 64, 2557–2564. [Google Scholar] [CrossRef]
  21. Li, W.; Ye, Y.; Zhang, K.; Feng, Z. A thickness measurement system for metal films based on eddy-current method with phase detection. IEEE Trans. Ind. Electron. 2017, 64, 3940–3949. [Google Scholar] [CrossRef]
  22. Huang, R.; Lu, M.; Zhang, Z.; Zhao, Q.; Xie, Y.; Tao, Y.; Meng, T.; Peyton, A.; Theodoulidis, T.; Yin, W. Measurement of the radius of metallic plates based on a novel finite region eigenfunction expansion (FREE) method. IEEE Sens. J. 2020, 20, 15099–15106. [Google Scholar] [CrossRef]
  23. Hu, G.; Huang, R.; Lu, M.; Zhou, L.; Yin, W. Measurement of radius of a metallic ball using eddy current testing based on peak frequency difference feature. Measurement 2021, 184, 109876. [Google Scholar] [CrossRef]
  24. Meng, X.; Lu, M.; Yin, W.; Bennecer, A.; Kirk, K.J. Evaluation of Coating Thickness Using Lift-Off Insensitivity of Eddy Current Sensor. Sensors 2021, 21, 419. [Google Scholar] [CrossRef] [PubMed]
  25. Wang, H.; Huang, J.; Liu, L.; Qin, S.; Fu, Z. A Novel Pulsed Eddy Current Criterion for Non-Ferromagnetic Metal Thickness Quantifications under Large Liftoff. Sensors 2022, 22, 614. [Google Scholar] [CrossRef] [PubMed]
  26. Desjardins, D.; Krause, T.W.; Clapham, L. Transient eddy current method for the characterization of magnetic permeability and conductivity. NDT E Int. 2016, 80, 65–70. [Google Scholar] [CrossRef]
  27. Lu, M.; Huang, R.; Yin, W.; Zhao, Q.; Peyton, A. Measurement of permeability for ferrous metallic plates using a novel lift-off compensation technique on phase signature. IEEE Sens. J. 2019, 19, 7440–7446. [Google Scholar] [CrossRef] [Green Version]
  28. Wen, D.; Fan, M.; Cao, B.; Xue, Z.; Wang, P. A PEC Thrice Subtraction Method for Obtaining Permeability Invariance Feature in Conductivity Measurement of Ferromagnetic Samples. Appl. Sci. 2019, 9, 2745. [Google Scholar] [CrossRef] [Green Version]
  29. Yu, Y.; Zhang, D.; Lai, C.; Tian, G. Quantitative approach for thickness and conductivity measurement of monolayer coating by dual-frequency eddy current technique. IEEE Trans. Instrum. Meas. 2017, 66, 1874–1882. [Google Scholar] [CrossRef]
  30. Wang, Z.; Yu, Y. Thickness and Conductivity Measurement of Multilayered Electricity-Conducting Coating by Pulsed Eddy Current Technique: Experimental Investigation. IEEE Trans. Instrum. Meas. 2019, 68, 3166–3172. [Google Scholar] [CrossRef]
  31. Wang, D.; Zhang, Z.; Yin, W.; Chen, H.; Ma, H.; Zhou, G.; Zhang, Y. Conductivity Classification of Multi-Shape Nonmagnetic Metal Considering Spatial Position Drift Effect with a Triple-Coil Electromagnetic Sensor. Sensors 2022, 22, 5694. [Google Scholar] [CrossRef]
  32. Tashiro, K.; Wakiwaka, H.; Minakata, K.; Kimura, T.; Nakamura, Y. A novel eddy current method for magnetic plate identification with elimination of lift-off effect. J. Jpn. Soc. Appl. Electromagn. Mech. 2019, 27, 165–168. [Google Scholar] [CrossRef] [Green Version]
  33. Lu, M.; Meng, X.; Huang, R.; Chen, L.; Peyton, A.; Yin, W. Liftoff tolerant pancake eddy-current sensor for the thickness and spacing measurement of nonmagnetic plates. IEEE Trans. Instrum. Meas. 2021, 70, 1–9. [Google Scholar] [CrossRef]
  34. Cheng, D.H.S. The reflected impedance of a circular coil in the proximity of a semi-infinite medium. IEEE Trans. Instrum. Meas. 1965, IM-14, 107–116. [Google Scholar] [CrossRef] [Green Version]
  35. Dodd, C.V.; Deeds, W.E. Analytical solutions to eddy-current probe-coil problems. J. Appl. Phys. 1968, 39, 2829–2838. [Google Scholar] [CrossRef] [Green Version]
  36. Simankova. L. Mathematical presentation of impedance variation of a coil cause by the measured object. TESLA Electron. 1971, 4, 112–118. [Google Scholar]
  37. Tytko, G.; Dziczkowski, L. Fast calculation of the filamentary coil impedance using the truncated region eigenfunction expansion method. Appl. Comput. Electromagn. Soc. J. (ACES) 2018, 33, 1461–1466. [Google Scholar]
  38. Theodoulidis, T.P.; Kriezis, E.E. Eddy Current Canonical Problems (with Applications to Nondestructive Evaluation); Tech Science Press: Forsyth, GA, USA, 2006. [Google Scholar]
  39. Tytko, G. Fast Method of Calculating the Air-Cored Coil Impedance Using the Filamentary Coil Model. PIER M 2020, 91, 101–109. [Google Scholar] [CrossRef]
  40. Tytko, G. Locating defects in conductive materials using the eddy current model of the filamentary coil. J. Nondestruct. Eval. 2021, 40, 66. [Google Scholar] [CrossRef]
  41. Babic, S.; Akyel, C. Calculating mutual inductance between circular coils with inclined axes in air. IEEE Trans. Magn. 2008, 44, 1743–1750. [Google Scholar] [CrossRef]
  42. Poletkin, K.V.; Korvink, J.G. Efficient calculation of the mutual inductance of arbitrarily oriented circular filaments via a generalisation of the kalantarov-zeitlin method. J. Magn. Magn. Mater. 2019, 483, 10–20. [Google Scholar] [CrossRef] [Green Version]
  43. Ishida, K.; Itaya, T.; Tanaka, A.; Takehira, N. Mutual inductance of arbitrary-shaped coils using shape functions. IET Sci. Meas. Technol. 2019, 13, 1085–1091. [Google Scholar] [CrossRef]
  44. Kumar, S.S.; Lee, C.H.; Mamleyev, E.; Poletkin, K. Calculation of Mutual Inductance Between Circular and Arbitrary Shaped Filaments: Segmentation Method. In Proceedings of the ACTUATOR 2022. International Conference and Exhibition on New Actuator Systems and Applications, Mannheim, Germany, 29–30 June 2022. [Google Scholar]
  45. Pirinççi, N.; Altun, H. A New Analytical Study on Mutual Inductance Calculations for Wireless Power Transfer Using Magnetic Vector Potential. IEEE Trans. Magn. 2022, 58, 8002714. [Google Scholar] [CrossRef]
  46. Poletkin, K.; Babic, S. Magnetic stiffness calculation for the corresponding force between two current-carrying circular filaments arbitrarily oriented in the space. J. Magn. Magn. Mater. 2022, 568, 170067. [Google Scholar] [CrossRef]
  47. Dziczkowski. L.; Tytko, G. A method for comparing the metrological properties of eddy current probes. IEEE Trans. Instrum. Meas. 2021, 70, 6004906. [Google Scholar]
  48. Zhu, Y.; Chen, B.; Luo, Y.; Zhu, R. Inductance calculations for coils with an iron core of arbitrary axial position. Electromagnetics 2019, 39, 99–119. [Google Scholar] [CrossRef]
  49. Kucheryavskiy, S.; Egorov, A.; Polyakov, V. Coupling Analytical Models and Machine Learning Methods for Fast and Reliable Resolution of Effects in Multifrequency Eddy-Current Sensors. Sensors 2021, 21, 618. [Google Scholar] [CrossRef]
  50. Zhang, S. Analytical Model of an I-core Coil for Nondestructive Evaluation of a Conducting Cylinder below an Infinite Plane Conductor. Meas. Sci. Rev. 2021, 21, 99–105. [Google Scholar] [CrossRef]
  51. Tytko, G. Eddy current testing of small radius conductive cylinders with the employment of an I-core sensor. Measurement 2021, 186, 110219. [Google Scholar] [CrossRef]
  52. Vasic, D.; Rep, I.; Špikic, D.; Kekelj, M. Model of Magnetically Shielded Ferrite-Cored Eddy Current Sensor. Sensors 2022, 22, 326. [Google Scholar] [CrossRef]
  53. Zhang, S.; Ye, C. Model of Ferrite-cored Driver-pickup Coil Probe Application of TREE Method for Eddy Current Nondestructive Evaluation. ACES J. 2022, 37, 632–638. [Google Scholar] [CrossRef]
Figure 1. Filamentary coil placed over a conductive cylinder.
Figure 1. Filamentary coil placed over a conductive cylinder.
Sensors 23 03267 g001
Figure 2. Rectangular cross–section of the filamentary coil located above the conductive half–space.
Figure 2. Rectangular cross–section of the filamentary coil located above the conductive half–space.
Sensors 23 03267 g002
Figure 3. Sensitivity of resistance and inductance of the filamentary coil to conductance variations of the examined workpiece material as well as to variations of the distance between the measuring coil and the examined surface. The graphs are plotted as functions of the generalised β parameter and for several values of α.
Figure 3. Sensitivity of resistance and inductance of the filamentary coil to conductance variations of the examined workpiece material as well as to variations of the distance between the measuring coil and the examined surface. The graphs are plotted as functions of the generalised β parameter and for several values of α.
Sensors 23 03267 g003
Figure 4. I-core sensor and plates made of bronze (B) and copper (C) used in measurements.
Figure 4. I-core sensor and plates made of bronze (B) and copper (C) used in measurements.
Sensors 23 03267 g004
Figure 5. Equivalent parameters of the air-core sensor for height h: (a) h = 0 mm, (b) h = 4 mm, (c) h = 8 mm, (d) h = 16 mm.
Figure 5. Equivalent parameters of the air-core sensor for height h: (a) h = 0 mm, (b) h = 4 mm, (c) h = 8 mm, (d) h = 16 mm.
Sensors 23 03267 g005
Figure 6. Equivalent parameters of the I-core sensor for height hc = 0 mm: (a) h = 0 mm, (b) h = 4 mm, (c) h = 8 mm, (d) h = 16 mm.
Figure 6. Equivalent parameters of the I-core sensor for height hc = 0 mm: (a) h = 0 mm, (b) h = 4 mm, (c) h = 8 mm, (d) h = 16 mm.
Sensors 23 03267 g006
Figure 7. Equivalent parameters of the air-core sensor and the I-core sensor: (a) h = 4 mm, (b) h = 0 mm and hc = 4 mm, (c) h = 4 mm and hc = 4 mm, (d) h = 8 mm and hc = 8 mm.
Figure 7. Equivalent parameters of the air-core sensor and the I-core sensor: (a) h = 4 mm, (b) h = 0 mm and hc = 4 mm, (c) h = 4 mm and hc = 4 mm, (d) h = 8 mm and hc = 8 mm.
Sensors 23 03267 g007
Figure 8. Equivalent parameters of the air-core sensor and the I-core sensor, for height h = 0 mm: (a) without the core, (b) hc = 0 mm, (c) hc = 4 mm, (d) hc = 8 mm.
Figure 8. Equivalent parameters of the air-core sensor and the I-core sensor, for height h = 0 mm: (a) without the core, (b) hc = 0 mm, (c) hc = 4 mm, (d) hc = 8 mm.
Sensors 23 03267 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Dziczkowski, L.; Tytko, G. Evaluation of the Properties of Eddy Current Sensors Based on Their Equivalent Parameters. Sensors 2023, 23, 3267. https://doi.org/10.3390/s23063267

AMA Style

Dziczkowski L, Tytko G. Evaluation of the Properties of Eddy Current Sensors Based on Their Equivalent Parameters. Sensors. 2023; 23(6):3267. https://doi.org/10.3390/s23063267

Chicago/Turabian Style

Dziczkowski, Leszek, and Grzegorz Tytko. 2023. "Evaluation of the Properties of Eddy Current Sensors Based on Their Equivalent Parameters" Sensors 23, no. 6: 3267. https://doi.org/10.3390/s23063267

APA Style

Dziczkowski, L., & Tytko, G. (2023). Evaluation of the Properties of Eddy Current Sensors Based on Their Equivalent Parameters. Sensors, 23(6), 3267. https://doi.org/10.3390/s23063267

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