# Effects of AC Coils Parameters on Transduction Efficiency of EMAT for Steel Plate Inspection

Authors: Xiaochun Song; Yawu Wang
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

In order to improve the transduction efficiency of electromagnetic acoustic transducer (EMAT) for steel plates inspection, the constitutive equation of magnetostrictive material was theoretically derived and simplified while the magnetostrictive force is parallel to the material surface. Based on the multiphysics field FEM, the effects of such excitation parameters as current, frequency, and pulse number in AC coils on magnetostrictive strain were mainly simulated, and the influence of the coil with different winding shapes on magnetostrictive strain was also analyzed. The simulation and experimental results indicate that magnetostrictive strain increases with a continuously increasing excitation current density, but it decreases with the increase of the frequency and pulse number of AC currents. Moreover, on condition that loop length and AC currents are held constant, spiral type coils have higher transduction efficiency than homocentric squares and figure-of-eight coils.


---

## Body

## 1. Introduction

Since the EMAT guided-wave technique can overcome such limitations as low-efficiency and high-cost of some conventional NDT techniques and realize the noncontact, long distance, and quick inspection, it shows some significant advantages for large-scale components inspection. And a typical EMAT, which is used to inspect the plate, mainly includes AC coils and bias magnets. When the high-frequency current is applied to the AC coil, the coil induces a time-varying magnetic field in the plate. And by coupling the alternating magnetic field in AC coils and the bias magnetic field provided by permanent magnets, an ultrasonic guided-wave is generated in the steel plate via the magnetostrictive effect. Therefore, the excitation parameters of AC coils are the key to influence the performance of EMAT.Due to the distinguishing characteristics of EMAT over other conventional approaches, simulation and optimization of the EMAT have been extensively studied by many researchers. Dutton et al. improved the efficiency of EMAT by adjusting the coil turns [1]. Jian et al. studied the lift-off effect of EMAT by establishing coil equivalent circuit [2] and also analyzed the influences of alternating magnetic field generated by AC current on the precision of EMAT [3, 4]. Yang et al. simulated generation mechanism of the electromagnetic acoustic wave in aluminum plates and analyzed its propagation characteristics [5]. Park et al. investigated the feasibility to inspect internal defects in thin aluminum plate [6, 7], and so on. Generally, the EMAT transduction mechanism is different with different materials; for example, the approach only relies on Lorentz force for nonferromagnetic materials, and magnetostrictive force plays leading role in generation of electromagnetic acoustic wave for ferromagnetic materials. Owing to the complexity of the magnetomechanical coupling, it is hard to describe the energy conversation mechanism of material magnetostriction accurately via current finite-element method relying on the analogy between piezoelectricity and piezomagnetism [8, 9]. With the high-efficiency numerical method being used widely, the FEM based on multiphysics field coupling has become a very useful technique to simulate the magnetostrictive effect. And by using COMSOL multiphysics software, the effects of such excitation parameters as current, frequency, and pulse number in AC coils on magnetostrictive strain were mainly simulated, and the influence of the coil with different shapes on magnetostrictive strain was also analyzed. Moreover, the simulation results are verified by experiments with the help of a self-developed EMAT guided-wave system.


## 2. Constitutive Equations of Ferromagnetic Materials

As the direction of bias magnetic field is parallel to the surface of the material and perpendicular to the direction of excitation current, the magnetostrictive force is parallel to the material surface. Thus, the magnetostrictive constitutive equations [10] can be written as
(1)
ε
x
=
1
E
σ
+
λ
s
(
1
-
σ
/
σ
s
)
M
w
s
2
M
x
2
-
θ
λ
s
[
M
x
4
-
M
w
s
4
(
1
-
σ
/
σ
s
)
4
]
M
w
s
4
,
(2)
H
x
=
1
k
f
-
1
(
M
x
M
s
)
-
λ
s
(
2
-
σ
/
σ
s
)
μ
0
M
s
2
M
x
W
+
4
θ
λ
s
σ
[
M
x
3
-
σ
s
M
w
s
3
(
1
-
σ
/
σ
s
)
3
]
μ
0
M
w
s
4
,
where (1) describes the magnetostrictive effect and (2) describes the inverse magnetostrictive effect. And the first item in (1) represents the inherent elasticity of materials and shows the linear relationship between materials strain and the applied stress, the second reflects the effects of the magnetic domain rotation on the elasticity due to stress anisotropy, and the last shows the influence of magnetomechanical coupling on the elasticity. And ε
x is response variables, λ
s is saturation magnetostriction, M
w
s is magnetic domain wall magnetization, M
x is magnetization intensity, μ
0 is the magnetic permeability of free space and is equal to 4
π
×
10
-
7 Hm−1, M
s is saturation magnetization, k and θ represent relaxation factor, k is related to saturation magnetic field, θ describes relations of M
x and the initial magnetization, f
(
x
)
=
coth
⁡
⁡
(
x
)
-
x
-
1, σ represents the applied stress, and σ
s represents the saturation stress.If the applied stress of the material is zero, the magnetostrictive constitutive relations can be simplified as follows:(3)
ε
x
=
M
x
2
λ
s
M
w
s
2
-
θ
λ
s
(
M
x
4
-
M
w
s
4
)
M
w
s
4
,
H
x
=
1
k
f
-
1
(
M
x
M
s
)
-
2
λ
s
μ
0
M
w
s
2
M
x
.


## 3. Effects of Excitation Parameters on Magnetostrictive Strain

### 3.1. Finite Element Model of EMAT for Steel Plates Inspection

A finite element model of EMAT which is used to inspect steel plates, shown as Figure1, is established by COMSOL multiple-physics software, which mainly consists of two permanent magnets and a homocentric-squares coil. The geometric dimensions and the properties parameters of the detected steel plate are shown as Table 1, and its B-H curve is shown as Figure 2. Material Nd-Fe-B N50 is chosen for the bias magnetic pole, and its magnetic properties parameters are shown as Table 2. In the mesh generation, such rectangular regions as permanent magnets and steel plates are divided by mapped-swept method firstly; then the rectangular grid is transformed into free tetrahedral elements using insert diagonal edges method, where the maximum element size is 5.28 mm, and the minimum element size is 0.384 mm. And other regions are meshed by adaptive free tetrahedral method, where the maximum element size is 14.4 mm, the minimum element size is 2.69 mm, the maximum element growth rate is 1.6, the curvature factor is 0.7, and the resolution of narrow regions is 0.4. The expression of excitation current in the homocentric-squares coil is shown as
(4)
i
=
0.5
×
I
×
1
-
cos
⁡
⁡
(
2
×
π
×
f
×
t
)
6
×
cos
⁡
⁡
(
a
×
π
×
f
×
t
)
,
where I represents current, f is excitation frequency, and a represents coefficient of pulse number.Table 1
Geometric dimensions and properties parameters of the detected steel plate.
Parameters
Length
Width
Thickness
Mass damping parameter
Stiffness damping parameter
Relative permittivity
Electrical conductivity
λ
s
M
ws
θ
k
Values
96 [mm]
70 [mm]
6 [mm]
300 [1/s]
3.2
e
-
5 [S]
1
4.032
e
6 [S/m]
1.5
e
6 [A/m]
1
e
6 [A/m]
1
4.3
e
-
4Table 2
Geometric parameters and magnetic properties parameters of permanent magnets.
Parameters
Bias magnet length
Bias magnet width
Bias magnet thickness
Electrical conductivity
Relative permeability
Residual flux density
Relative permittivity
Values
50 [mm]
20 [mm]
15 [mm]
1 [S/m]
1.26
0.38 [T]
1Figure 1
FEM of EMAT for plates inspection.Figure 2
B-H curve of steel plates.


### 3.2. Effects of AC Current Density on Magnetostrictive Vibration

According to (4), the initial values of loaded current in AC coils are specified as I
=
10 A, f
=
120 kHz, and a
=
2, respectively. Using of the COMSOL finite element analysis software, the magnetostritive displacements in steel plates are calculated with different excitation parameters in order to analyze the influences of excitation parameters on the magnetostrictive strain. The results can be used to set excitation parameters in accordance with geometric dimensions of the detected steel plate and obtain the best transduction efficiency in practical inspection.Figure3(a) shows the influences of the excitation current on magnetostritive displacements, where the current density is from 1 to 5 A, respectively. It can be seen that magnetostrictive vibration increases continuously with the increasing of AC current intensity. Figure 3(b) shows the relation between the magnetostritive displacements and the current intensity for different steel plate thickness (
D
), and it is nearly linear if the excitation frequency and pulse number are held constant. Moreover, it is found that the magnetostritive displacement decreases gradually with a continuously increasing plate thickness.Effects of excitation current on magnetostrictive strain.
(a)
AC current-magnetostrictive strain
(b)
AC current-magnetostriction with different plate thickness


### 3.3. Effects of Excitation Frequency on Magnetostrictive Vibration

Figure4 shows the influences of the excitation frequency on magnetostritive displacements if the other parameters are held constant, where the AC current frequency is from 100 to 500 kHz, respectively. It is found that both the magnetostrictive vibration amplitude and the magnetostrictive period decrease with the increase of the current frequency loaded in AC coils. A possible cause of this phenomenon is that the vibration period could shorten simultaneously and cause the times of particle vibration passing balance point per unit time to increase as the excitation frequency increases. Moreover, since the material damping increases with the decrease of particle vibration period, the magnetostrictive amplitude decreases. From Figure 4(b), it can be seen that the magnetostrictive amplitude decreases with the excitation frequency. And it decreases sharply at first, but when the excitation frequency exceeds a certain value and the plate undergoes gradual magnetostrictive saturation, and the increment of magnetostrictive displacement becomes slow. Additionally, on condition that excitation parameters are held constant, the magnetostritive displacement decreases with the plate thickness.Effects of excitation frequency on magnetostrictive strain.
(a)
Current frequency-magnetostrictive strain
(b)
Frequency-magnetostriction with different plate thickness


### 3.4. Effects of AC Current Pulse Number on Magnetostrictive Vibration

Figure5(a) shows the relations of the AC current pulse number on magnetostritive displacements, where the pulse number is from 2 to 6, respectively. It is found that the increase of pulse number could decrease the transduction efficiency and cause the amplitude of magnetostrictive vibration to decrease. A possible reason of this phenomenon is that the times of particle vibration passing balance point per unit time could increase simultaneously and cause the velocity of magnetostrictive vibration to increase as the pulse number increases. And because the material damping increases with the decrease of particle vibration velocity, the magnetostrictive amplitude decreases. From Figure 5(b), it can be seen that the magnetostrictive amplitude decreases sharply at first with the pulse number, and when the pulse number exceeds 4 and the plate undergoes gradual magnetostrictive saturation, the increment of magnetostrictive displacement becomes slow. Similarly, on condition that excitation parameters are held constant, the magnetostritive displacement also decreases with the plate thickness.Effects of pulse number on magnetostrictive strain.
(a)
Pulse number-magnetostrictive strain
(b)
Pulse number-magnetostriction with different plate thickness


## 3.1. Finite Element Model of EMAT for Steel Plates Inspection

A finite element model of EMAT which is used to inspect steel plates, shown as Figure1, is established by COMSOL multiple-physics software, which mainly consists of two permanent magnets and a homocentric-squares coil. The geometric dimensions and the properties parameters of the detected steel plate are shown as Table 1, and its B-H curve is shown as Figure 2. Material Nd-Fe-B N50 is chosen for the bias magnetic pole, and its magnetic properties parameters are shown as Table 2. In the mesh generation, such rectangular regions as permanent magnets and steel plates are divided by mapped-swept method firstly; then the rectangular grid is transformed into free tetrahedral elements using insert diagonal edges method, where the maximum element size is 5.28 mm, and the minimum element size is 0.384 mm. And other regions are meshed by adaptive free tetrahedral method, where the maximum element size is 14.4 mm, the minimum element size is 2.69 mm, the maximum element growth rate is 1.6, the curvature factor is 0.7, and the resolution of narrow regions is 0.4. The expression of excitation current in the homocentric-squares coil is shown as
(4)
i
=
0.5
×
I
×
1
-
cos
⁡
⁡
(
2
×
π
×
f
×
t
)
6
×
cos
⁡
⁡
(
a
×
π
×
f
×
t
)
,
where I represents current, f is excitation frequency, and a represents coefficient of pulse number.Table 1
Geometric dimensions and properties parameters of the detected steel plate.
Parameters
Length
Width
Thickness
Mass damping parameter
Stiffness damping parameter
Relative permittivity
Electrical conductivity
λ
s
M
ws
θ
k
Values
96 [mm]
70 [mm]
6 [mm]
300 [1/s]
3.2
e
-
5 [S]
1
4.032
e
6 [S/m]
1.5
e
6 [A/m]
1
e
6 [A/m]
1
4.3
e
-
4Table 2
Geometric parameters and magnetic properties parameters of permanent magnets.
Parameters
Bias magnet length
Bias magnet width
Bias magnet thickness
Electrical conductivity
Relative permeability
Residual flux density
Relative permittivity
Values
50 [mm]
20 [mm]
15 [mm]
1 [S/m]
1.26
0.38 [T]
1Figure 1
FEM of EMAT for plates inspection.Figure 2
B-H curve of steel plates.


## 3.2. Effects of AC Current Density on Magnetostrictive Vibration

According to (4), the initial values of loaded current in AC coils are specified as I
=
10 A, f
=
120 kHz, and a
=
2, respectively. Using of the COMSOL finite element analysis software, the magnetostritive displacements in steel plates are calculated with different excitation parameters in order to analyze the influences of excitation parameters on the magnetostrictive strain. The results can be used to set excitation parameters in accordance with geometric dimensions of the detected steel plate and obtain the best transduction efficiency in practical inspection.Figure3(a) shows the influences of the excitation current on magnetostritive displacements, where the current density is from 1 to 5 A, respectively. It can be seen that magnetostrictive vibration increases continuously with the increasing of AC current intensity. Figure 3(b) shows the relation between the magnetostritive displacements and the current intensity for different steel plate thickness (
D
), and it is nearly linear if the excitation frequency and pulse number are held constant. Moreover, it is found that the magnetostritive displacement decreases gradually with a continuously increasing plate thickness.Effects of excitation current on magnetostrictive strain.
(a)
AC current-magnetostrictive strain
(b)
AC current-magnetostriction with different plate thickness


## 3.3. Effects of Excitation Frequency on Magnetostrictive Vibration

Figure4 shows the influences of the excitation frequency on magnetostritive displacements if the other parameters are held constant, where the AC current frequency is from 100 to 500 kHz, respectively. It is found that both the magnetostrictive vibration amplitude and the magnetostrictive period decrease with the increase of the current frequency loaded in AC coils. A possible cause of this phenomenon is that the vibration period could shorten simultaneously and cause the times of particle vibration passing balance point per unit time to increase as the excitation frequency increases. Moreover, since the material damping increases with the decrease of particle vibration period, the magnetostrictive amplitude decreases. From Figure 4(b), it can be seen that the magnetostrictive amplitude decreases with the excitation frequency. And it decreases sharply at first, but when the excitation frequency exceeds a certain value and the plate undergoes gradual magnetostrictive saturation, and the increment of magnetostrictive displacement becomes slow. Additionally, on condition that excitation parameters are held constant, the magnetostritive displacement decreases with the plate thickness.Effects of excitation frequency on magnetostrictive strain.
(a)
Current frequency-magnetostrictive strain
(b)
Frequency-magnetostriction with different plate thickness


## 3.4. Effects of AC Current Pulse Number on Magnetostrictive Vibration

Figure5(a) shows the relations of the AC current pulse number on magnetostritive displacements, where the pulse number is from 2 to 6, respectively. It is found that the increase of pulse number could decrease the transduction efficiency and cause the amplitude of magnetostrictive vibration to decrease. A possible reason of this phenomenon is that the times of particle vibration passing balance point per unit time could increase simultaneously and cause the velocity of magnetostrictive vibration to increase as the pulse number increases. And because the material damping increases with the decrease of particle vibration velocity, the magnetostrictive amplitude decreases. From Figure 5(b), it can be seen that the magnetostrictive amplitude decreases sharply at first with the pulse number, and when the pulse number exceeds 4 and the plate undergoes gradual magnetostrictive saturation, the increment of magnetostrictive displacement becomes slow. Similarly, on condition that excitation parameters are held constant, the magnetostritive displacement also decreases with the plate thickness.Effects of pulse number on magnetostrictive strain.
(a)
Pulse number-magnetostrictive strain
(b)
Pulse number-magnetostriction with different plate thickness


## 4. Effects of AC Coil Shapes on Transduction Efficiency of EMAT

In order to study the influence of AC coil types on transduction efficiency of EMAT, a 235 mm long enameled wire is wound in such three shapes as spiral type coils, figure-of-eight, and homocentric squares. According to the above analysis, an AC current with the same excitation parameters is loaded in the AC coils, and some calculated results are shown in Figure6. It is found that the maximum magnetostrictive displacement of spiral type coils, figure-of-eight, and homocentric squares are 2.63703
×
10
-
9 mm, 9.47149
×
10
-
10 mm, and 1.89573
×
10
-
9 mm, respectively. Obviously, on condition that loop length and excitation current are constant, the spiral type coil has higher transduction efficiency than figure-of-eight coils and homocentric squares.Magnetostrictive displacement with different coil types.
(a)
Magnetostriction of spiral type coils
(b)
Magnetostriction of figure-of-eight coilsFigures7 and 8 show the influence of excitation parameters on transduction efficiency of EMAT with spiral type coils and figure-of-eight coils, respectively, where the loaded AC current density is from 1 to 5 A, the excitation frequency is from 100 to 500 kHz and the pulse number is from 2 to 6. It is found that the generated magnetostriction of AC coils with different shapes show the same influence with variation of excitation parameters. That is to say, the transduction efficiency of different AC coils increases with an increasing of AC current density, but it decreases with an increasing of the excitation frequency and pulse numbers.Magnetostriction of spiral type coils.
(a)
Magnetostriction-current density
(b)
Magnetostriction-excitation frequency
(c)
Magnetostriction-pulse numberMagnetostriction of figure-of-eight coils.
(a)
Magnetostriction-current density
(b)
Magnetostriction-excitation frequency
(c)
Magnetostriction-pulse number


## 5. Experimental Analysis

According to the above analysis, in the experimental condition, a steel plate with 6 mm thickness, 80 mm width, and 400 mm length is selected as the inspected sample, and a through hole with 3 mm diameter is drilled on the surface of the plate. An EMAT, whose AC coils are spiral shapes and bias pole pieces are two 20 mm wide and 15-mm-height permanent magnets, is applied to inspect the plate and to verify the simulation results. Adjusting the button of current, frequency, and pulse number of the exciting signal source as shown in Figure9, some testing signals can be obtained and shown in Figures 10, 11, and 12, respectively. Exacting the maximum amplitude of defect signals, the effects of AC excitation parameters on the echo signals are shown as Figure 13. It can be seen that the magnetostrictive vibration increases continuously with the applied current density, but the magnetostrictive vibration decreases with the applied excitation frequency and pulse number. Compared with the FEM calculated results shown in Figures 3(b), 4(b), and 5(b), it is found that the experimental results are basically consistent with simulation data, which indicates that the finite element model is correct and simulation results are reliable.Figure 9
Exciting signal source.Testing signals with different current density.
(a)
I
=
3 A
(b)
I
=
4 A
(c)
I
=
5 ATesting signals with different excitation frequency.
(a)
f
=
100 kHz
(b)
f
=
200 kHz
(c)
f
=
300 kHzTesting signals with different pulse numbers.
(a)
a
=
2
(b)
a
=
3
(c)
a
=
4Testing signals voltage with different AC coils parameters.
(a)
AC current-testing signals voltage
(b)
Current frequency-testing signals voltage
(c)
Pulse number-testing signals voltage


## 6. Conclusions

Based on the constitutive equation of magnetostrictive material, a finite element model of EMAT for steel plates inspection is established by COMSOL multiple-physics software, the effects of AC excitation parameters on magnetostrictive displacements have been mainly simulated, and the influence of AC coil shapes on EMAT transduction efficiency are also analyzed. Moreover, successful comparisons have been obtained by a validation experimental case for a steel plate with a through hole in the center. Works in progress are focusing on optimizing transmitters or receivers of EMAT guided-wave, which can provide theoretical and methodological basis for choosing the excitation parameters in practical EMAT applications.



---

*Source: 712418-2014-09-28.xml*