Experimental and numerical evaluation of electromagnetic acoustic ...

10 downloads 152 Views 787KB Size Report
Sep 1, 2011 ... Electromagnetic Acoustic Transducers (EMATs) are an attractive ... electromagnetic counterparts, i.e. ultrasonic waves in the testpiece.
NDT&E International 45 (2012) 32–38

Contents lists available at SciVerse ScienceDirect

NDT&E International journal homepage: www.elsevier.com/locate/ndteint

Experimental and numerical evaluation of electromagnetic acoustic transducer performance on steel materials R. Ribichini a,n, F. Cegla a, P.B. Nagy a,b, P. Cawley a a b

UK Research Centre in NDE, Department of Mechanical Engineering, Imperial College, London, SW7 2AZ, UK School of Aerospace Systems, University of Cincinnati, Cincinnati, OH 45221, USA

a r t i c l e i n f o

a b s t r a c t

Article history: Received 12 December 2010 Received in revised form 14 August 2011 Accepted 20 August 2011 Available online 1 September 2011

Electromagnetic Acoustic Transducers (EMATs) are an attractive alternative to standard piezoelectric probes in a number of applications thanks to their contactless nature. EMATs do not require any couplant liquid and are able to generate a wide range of wave-modes; however these positive features are partly counterbalanced by a relatively low signal-to-noise ratio and by the dependence of EMAT performance on the material properties of the test object. A wide variety of steel materials is employed in many industrial applications, so it is important to assess the material-dependent behaviour of EMATs when used in the inspection of different types of steel. Experimental data showing the performance of bulk shear wave EMATs on a wide range of steels is presented, showing the typical range of physical properties encountered in practice. A previously validated Finite Element model, including the main transduction mechanisms, the Lorentz force and magnetostriction, is used to evaluate the experimental data. The main conclusion is that the Lorentz force is the dominant transduction effect, regardless of the magnitude and direction of the bias magnetic field. Differently from magnetostriction, the Lorentz force is not significantly sensitive to the typical range of physical properties of steels, as a consequence the same EMAT sensor can be used on different grades of ferritic steel. & 2011 Elsevier Ltd. All rights reserved.

Keywords: Electromagnetic Acoustic Transducers Magnetostriction Lorentz force Steel

1. Introduction Electromagnetic Acoustic Transducers (EMATs) are able to generate and detect ultrasonic waves thanks to contactless electromagnetic coupling with the test object, rather than with mechanical coupling, as in standard piezoelectric probes [1–4]. This feature makes EMATs an attractive alternative to piezoelectric transducers in all those applications where contactless inspections are required, for example when high temperature or moving objects are to be tested. Moreover, EMATs can excite a wide range of wave-modes and can be employed as a standard for ultrasonic calibration. However, EMATs have some disadvantages: the signal-to-noise ratio is relatively low compared to standard transducers and their performance depends significantly on the material properties of the inspected sample. A wide range of different kinds of steel materials, with different physical properties is employed in modern engineering. The variation of EMATs performance with material properties represents a major concern for practical applications, since it raises the question whether the same EMAT probe can be successfully used to inspect different kinds of steel, or if transducers optimized for each steel grade have to be developed. For instance, the results of pulse-echo

n

Corresponding author. E-mail address: [email protected] (R. Ribichini).

0963-8695/$ - see front matter & 2011 Elsevier Ltd. All rights reserved. doi:10.1016/j.ndteint.2011.08.007

tests performed with the same EMAT transducer on different steel grades (with no oxidation) is presented in Fig. 1. When employed on ferromagnetic materials such as steel, EMATs exploit mainly two different types of transduction mechanisms: the Lorentz force and magnetostriction. When an eddy current density, Je, induced in a conductive sample by a driving electric current, interacts with a static bias magnetic flux density, B produced by a magnet, the Lorentz force arises: f ¼ Je  B,

ð1Þ

where f is a force per unit volume. On the other hand, magnetostriction is due to the fact that ferromagnetic domains tend to align along the direction of the total magnetic field, causing a net mechanical strain [5,6]. When the magnetic field has a timevarying component this strain can be exploited to launch ultrasonic waves. Both transduction mechanisms have mechanical to electromagnetic counterparts, i.e. ultrasonic waves in the testpiece induce an electrical signal in the transducer so it can be used as a receiver. While the Lorentz force mechanism is linear and relatively insensitive to material properties such as electric conductivity s and relative magnetic permeability mr, magnetostriction is highly non-linear, depends significantly on the physical properties of the sample and is a function of the applied magnetic field, stress state, magneto-mechanical loading history and surface conditions [7].

R. Ribichini et al. / NDT&E International 45 (2012) 32–38

33

Am m plit u de [a rrb.]

Table 1 List of the steel samples under investigation.

0

60

Designation

C [%]

Other elements [%]

Notes

EN8 EN16 EN24 EN36 EN3 EN32B BO1 AISI 304 L80a L80b L80SS TN80cr3 J55 CS70

0.32–0.40 0.30–0.40 0.30–0.40 0.10 0.16–0.24 0.13–0.18 0.90–1.00 0.08 0.25–0.30 0.25–0.30 0.25–0.30 0.25–0.30 0.40–0.50 0.65–0.75

0.80 Mn 1.50 Mn, 0.25 Mo 0.60 Mn, 0.25 Mo, 1.50 Ni, 1.20 Cr 0.50 Mn, 3.50 Ni, 0.90 Cr 0.70 Mn 0.80 Mn 1.20 Mn, 0.50 Cr, 0.50 W, 0.22V 9.00 Ni, 19.00 Cr 1.40 Mn. 0.12 Cu, Mo, Cr, Ti 1.40 Mn. 0.12 Cu, Mo, Cr, Ti Mo, Cr, Ti Mo, Cr, Ti 1.00 Mn, 0.17 Cr, 0.09 Cu, Mo, Ni 0.70 Mn

Mild steel Hardenable Hardenable Hardenable Mild steel Mild steel Tool steel Austenitic Pipe steel Pipe steel Pipe steel Pipe steel Pipe steel Pipe steel

120

Time [µs] Fig. 1. Pulse-echo signals from a spiral coil EMAT on a range of different steel grade samples with equal thickness: from top down EN24, AISI 304, CS70, TN80cr3, L80. Oxide layers can significantly increase the signal amplitude due to magnetostriction.

For this reason, it is fundamental to determine, which transduction mechanism dominates for a given EMAT configuration as it affects the behaviour of the transducer when used on materials with different properties. Previous research has established that magnetostriction is the leading phenomenon in those EMAT configurations where the bias field is parallel to the surface of the sample [3]. However, when the static field is normal to the sample, some authors state that the Lorentz force dominates [3,8,9], while others [4,10] claim that magnetostriction is the major effect for most practical cases. This paper presents an experimental study of bulk shear wave EMAT performance on a wide range of steel materials commonly used in engineering. The experimental results are compared with simulations from a previously validated Finite Element model in order to obtain a physical interpretation of the data. An analysis of the relative importance of the transduction mechanisms is performed and practical conclusions are drawn.

2. The experimental study The steel grades under investigation are among the most commonly used in modern engineering, ranging from mild steel to tool and alloy steel, and including pipe steel and an austenitic steel (AISI 304). Two samples of the same grade (L80) were also included for reference. The materials tested are listed in Table 1. All the samples have the same dimensions: 70  30  4 mm. 2.1. Conductivity and permeability measurement The electrical conductivity s and the relative magnetic permeability mr, of each sample were measured with the alternating current potential drop (ACPD) technique. A pair of electrodes injects an alternating current in the testpiece and a second pair of electrodes measures the resulting potential drop; the resistance can then be computed as the real part of the ratio between the potential difference and the current. Resistance varies with frequency due to the electromagnetic skin depth effect; once the geometric configuration of the probe and the thickness of the sample are known, analytical solutions [11] can be employed to compute the couple {s, mr} that minimizes the root mean square error between theoretical and experimental data. The resistance of each sample was measured in the frequency range between 2 and

Table 2 Measured electromagnetic properties. Relative magnetic permeability was estimated via ACPD technique and a Feritscope instrument. The electric conductivity was measured with ACPD. The static magnetic flux density B when the spiral coil EMAT is applied to each sample, measured at the centre of the transducer and on the surface of the sample, is also given. Data accuracy 75 mT. Designation

r [MS/m]

lr, ACPD

lr, Ferritscope

B [mT]

EN8 EN16 EN24 EN36 EN3 EN32B BO1 AISI 304 L80a L80b L80SS TN80cr3 J55 CS70

4.12 3.71 3.80 3.03 4.47 4.46 4.01 1.39 4.54 4.54 4.19 2.61 4.06 3.77

92 52 65 99 128 108 90 1 70 61 67 86 137 59

170 110 132 142 166 150 157 1 143 139 126 167 164 100

753 769 774 762 772 772 777 433 780 760 768 756 766 747

400 Hz and a two-variable fit with the analytical formula was performed to deduce the electromagnetic properties (Table 2). The measured conductivity of the ferritic steels falls within the range s A [2.5; 4.5] MS/m while the conductivity of austenitic AISI 304 is 1.39 MS/m. The measured permeabilities for the ferritic steel samples were between 50 and 140, while AISI 304, being nonferromagnetic, has approximately unit relative permeability. The magnetic permeability was also measured with a Fischer Technology Feritscope MP30E-S. This instrument measures an engineering parameter, the equivalent ferrite content, from which permeability can be estimated using an approximated relationship found in the literature [12]. While ACPD employs low-intensity currents, in the order of a few milliamperes, the Feritscope induces much larger currents in the sample. The EMATs used in the experimental study were driven by an approximately 10 A peak to peak pulse, and the currents used by the Feritscope are closer to the actual experimental conditions than the ACPD ones, however, this instrument gives much less accurate values, reported in Table 2. 2.2. Magnetostriction measurement The magnetostrictive curves, i.e. magnetostrictive strain against magnetic field strength, of four steel grades (EN3, EN24, EN32B, BO1) were measured in order to determine the magnetostrictive parameters to be fed in the numerical model. In each measurement, a small sample (30  20  1 mm) was placed in the

R. Ribichini et al. / NDT&E International 45 (2012) 32–38

2.3. EMAT wave amplitude measurement Two commercial transducers (Sonemat Ltd.) have been used: a spiral coil EMAT and a linear racetrack coil EMAT. Both transducers generate shear waves, with radial and linear polarization, respectively, propagating in the bulk of the material. The static magnetic field is normal to the surface of the sample and is due to a permanent magnet (NdFeB), while the coil generates eddy current and dynamic magnetic fields parallel to the surface of the sample. The transducers are driven by a broadband pulse, whose centre frequency is around 2.5 MHz. The result of a typical pulse-echo test is shown in Fig. 3 (a): the ultrasonic pulse travels across the thickness of the sample and the reflections from its back-wall are received by the transducer. For each type of transducer, five

Magnetostrictiv Strain, [ppm]

10

0

-10

EN32B EN32

-20

BO1 EN24 nickel

-30

-40 0

20

40

60

80

Static Magnetic Field, H [kA/m] Fig. 2. Magnetostriction curves of four steel grades and industrially pure nickel.

A Amplitu ude [V V]

0.3

-0.3

0

20

Time [µs] 0.3

Vp-pp [V]

air gap of a magnetic circuit. Two electromagnets driven by an adjustable DC current generated the bias field, the resulting magnetic field being proportional to the driving current. The magnetic flux density generated at the surface of the sample, in a direction parallel to the surface (Bair) was measured by using a Hall gaussmeter (GM04, Hirst Magnetic Instruments). The magnetic field strength inside the material can then be estimated by acknowledging that Hair ¼Bair/m0 and that the boundary conditions for H prescribe the continuity of the tangential component at the boundary between two media, so Hsteel ffi Hair. Since the magnetostrictive strain to be measured is relatively small (less than 8 ppm) four strain gages (Kyowa) in a full bridge configuration were employed. Two gages on the opposite arms of the Wheatstone bridge were parallel to the static bias field, while the other two gages were perpendicular to it. Two gages orthogonal to each other were on each side of the sample; this configuration maximizes the sensitivity to the strain in the bias field direction while cancelling out any bending strain or thermal expansion strain. The resulting magnetostriction curves, shown in Fig. 2, are consistent with data available in the literature [5,7,10,13]. For comparison, the magnetostriction curve of industrially pure (99.0%) nickel is also shown [14]. In the steel samples the application of a magnetic field initially causes a positive strain (i.e. an expansion) along the direction of the field. The deformation reaches a maximum for Ho20 kA/m and turns into a compressional strain for higher bias fields. Even though the shapes of the four curves are similar, the position and magnitude of the maxima differ significantly for each grade because of the presence of alloy elements and due to thermal treatments. Conversely, nickel shows a monotonic contraction whose amplitude is significantly larger than the strain observed in any steel.

Extrapolated Value 0 20

0

Time [µs] Fig. 3. (a) Signal received by an EMAT transducer in a pulse-echo test. The peak to peak amplitudes of the back-wall reflections have been interpolated via an exponential fit (b). It is then possible to estimate a theoretical attenuation-free amplitude for zero time of flight.

0.6

A Adjuste ed Signnal Am mplitudee [ V]

34

J55 TN80

L80SS

L80

EN36 EN16

0.4

EN3 EN32B

EN24

BO1 EN8 304

0.2

0.0 1

2

3

4

5

Fig. 4. Experimental EMAT amplitudes on different steels plotted against their electric conductivity. The amplitudes are attenuation compensated and squarerooted to account for the wave generation process only.

acquisitions per steel sample were taken, each resulting from the average of 1000 time traces. The peak to peak amplitudes of the first seven reflections were measured and were fitted with an exponential function, in order to extrapolate the theoretical amplitude for zero time of fight (Fig. 3 (b)). This is necessary in order to compensate both for diffraction effects and for the ultrasonic attenuation, which is different for each kind of steel. Since the tests used the EMAT in pulse echo mode, the square root of the values obtained was taken in order to account for the generation mechanism only, on the assumption that reciprocity holds. Being magnetostriction highly non-linear, it is in general non-reciprocal. However, a linearization can be employed when a small dynamic ~ is superimposed on a large static bias field H, such that field H ~ Hb H. This assumption is normally satisfied in EMATs, as the bias field due to the magnet is usually much larger than the timevarying field caused by the driving current. The experimental results are shown in Figs. 4 and 5 for the linear coil transducer. The adjusted signal amplitudes are plotted against electric conductivity s (Fig. 4) and against magnetic permeability

R. Ribichini et al. / NDT&E International 45 (2012) 32–38

A Adjuste ed Signnal Am mplitudee [ V]

0.6

0.4

0.2

0.0 40

70

100 130 Magnetic permeability,  r

160

Fig. 5. Experimental EMAT amplitudes on different steels plotted against their magnetic permeability as measured with ACPD technique (AISI 304 not shown in this graph as mr ¼ 1). The amplitudes are attenuation compensated and squarerooted to account for the wave generation process only.

(measured with ACPD technique), mr (Fig. 5). Error bars show the experimental standard deviations of the quantities under investigation for each steel grade. Analogous graphs were obtained from the spiral coil EMAT. The data show that the signal amplitudes do not have a large scatter and are not obviously correlated with the electric conductivity and magnetic permeability. Even using the permeabilities values measured with the Feritscope there is no better correlation between EMAT amplitudes and permeabilities. The only exception is the case of austenitic steel whose lower amplitude is due to the fact that since this material is not ferromagnetic the magnetic flux density is significantly smaller than in the case of ferromagnetic steels. Indeed, measurements indicated that B ¼ 410 mT for AISI 304, against an average of B ffi 770 mT for all the other samples (Table 2); this reduces the resulting amplitude by a factor of about 2, as the Lorentz force is linear in B (Eq. (1)). If we compensate the amplitude of AISI 304 for this effect, all the experimental points have similar amplitudes. This strongly suggests that the transduction is mainly due to the Lorentz force, whose magnitude does not depend significantly on conductivity or permeability; if magnetostriction were dominant, a much larger scatter would be expected because of the observed differences in the magnetostriction curves of the various grades. In order to test this hypothesis and shed light on the experimental results numerical simulations were carried out.

3. Finite element simulations An EMAT numerical model has been developed using a Finite Element (FE) commercial software, COMSOL Multiphysics [14]. The program solves simultaneously the electrodynamic problem, accounting for eddy-current induction, and the elastic problem, accounting for wave generation. The magnetomechanical coupling is achieved by adding the Lorentz force and magnetostriction. The Lorentz force is implemented using its definition: the input force of the elastic problem, a mechanical effect, is caused by electrical quantities i.e. the vector product of eddy current density and static magnetic flux density (Eq. (1)). Modelling magnetostriction requires a modification of the constitutive equations in a way analogous to piezoelectricity: ( e~ ¼ SH s~ þdH~ , ð2Þ ~ B~ ¼ dT s~ þ ms H

35

~ are where e~ and s~ are the strain and stress tensors and B~ and H the magnetic flux density and the magnetic field strength, respectively. Dynamic components of physical quantities are here indicated with a tilde as opposed to the static components denoted by a bar. Together with the usual elastic compliance matrix SH (measured for constant H) and permeability matrix ms (measured at constant stress), coupling terms, proportional to the magnetostriction matrix d, (and its transpose dT) are present. The first equation accounts for direct magnetostriction, i.e. strain caused by the application of a magnetic field, whereas the second equation describes inverse magnetostriction, used in the detection process, when a stress produces magnetic flux density changes that can be picked up by the transducer. Under the large ~ it is possible to derive all the bias field approximation, Hb H, components of the magnetostrictive coupling matrix d from a single experimental curve of the magnetostrictive strain versus the applied bias field H, which is characteristic of the ferromagnetic material under investigation. All the non-zero terms of the matrix are either proportional to the derivative of the magnetostriction curve at the operation point, i.e. the bias field H, or to the ratio between the total magnetostrictive strain, e and the corresponding static bias field. The latter is usually the most relevant for magnetostrictive wave generation; Ogi and Hirao [4] have shown that it can be computed as: d15 ¼

3e H

ð3Þ

The magnetostriction model (and its underlying assumptions) has been quantitatively validated for SH0 wave generation in a nickel plate [14]. The model predicted the wave amplitude from first principles, without any adjustable parameters, with a 720% accuracy over a 200 kHz range and over a wide range of bias field strength. The model was used to help to understand the results of the experiments on the steel samples discussed above. An axisymmetric two-dimensional model in a cylindrical reference system {r,z,j} of an EMAT has been developed. The driving current in the coil is modelled as a zero cross-section current sheet, flowing in the circumferential direction above the metal, that induces eddy currents Jj. These interact with the vertical component of the static flux density Bz producing a Lorentz body force fr ¼ Jj  Bz in the radial direction that generates shear waves. Magnetostriction also contributes to the wave generation since it can be shown from Eq. (2) that shear strains e~ rz are produced by the dynamic magnetic field:

e~ rz pd15 H~ r

ð4Þ

where d15 is the magnetostrictive coupling matrix component ~ r is the radial component involved in shear wave generation and H of the dynamic magnetic field strength. The other magnetostriction contributions are proportional to the normal component of ~ z , and are considerably lower as the magnetic the dynamic field, H ~ below the coil is mostly parallel to the surface of the field H sample. The outer and inner diameters of the coil are 34 mm and 6 mm, respectively; the distance between the coil and the sample (lift-off) is 0.6 mm. The coil is driven by a 1 A current oscillating at a frequency f¼2 MHz. The mesh consists of approximatively 150,000 triangular elements. The elastic properties used were the same for all the grades of steel: Young’s modulus 200 GPa, Poisson’s ratio 0.33, mass density 7850 kg/m3. Just below the coil, full magnetostrictive constitutive equations are employed to simulate the transduction process. For a depth larger than a few skin depths d, i.e. 9z9 44d, the dynamic magnetic field becomes negligible and no transduction occurs. For this reason, purely elastic constitutive equation can be used to describe wave propagation saving significant computational time. In order to

R. Ribichini et al. / NDT&E International 45 (2012) 32–38

simulate the operation on a half-space, an absorbing region with finite damping constant surrounds the elastic domain, to avoid back-reflections from the boundaries of the model. The result of a typical FE simulation is shown in Fig. 6. The displacement amplitudes produced separately by the Lorentz force and by magnetostriction were computed for four of the steel samples (EN3, EN24, EN32B, BO1). In turns, the sole Lorentz force was applied, without any magnetostriction, and then the simulation was repeated with purely magnetostrictive effects and no Lorentz force, in order to evaluate the contribution of each mechanism. The magnetostrictive and magnetic properties were obtained from the experiments discussed above. In order to test the hypothesis that the Lorentz force is the dominant effect, the most favourable conditions for magnetostriction were taken into account to assess its maximum contribution. The magnetic permeabilities used in the simulations were those measured via ACPD, which are lower than those estimated with the Feritscope. Lower permeabilities imply a larger skin depth as d ¼ ðpf smÞ1=2 , ~ that is, there is a larger region where a significant dynamic field H is present. In other words, this means that the area over which Eq. (4) has to be integrated is wider hence the effect of magnetostriction is stronger. Moreover, the magnetic bias field in the material H, which determines the operation point cannot be estimated without a degree of uncertainty. This is a consequence of the fact that at the boundary between two media the perpendicular component of B (in our case Bz ) is continuous, while the perpendicular component of H is discontinuous. In other words, we know accurately the value of Bz from experimental data, but we can only estimate Hz using FE models. For the case under study it was found that Hz A [6, 15] kA/m. The maximum values of the magnetostrictive constant d15 falling in this range were considered to assess the largest possible impact of magnetostriction on wave generation. For the Lorentz force computations the static bias field Bz was assumed to be the same for all the samples and was set to the experimental value: Bz ¼ 770mT. Remembering that the values for magnetostriction are to be considered an ideal upper limit, the simulations indicate that for the investigated steels, the Lorentz force is the main transduction mechanism and that the contribution of magnetostriction is never larger than  30% of the

Lorentz force for three samples, and reaches  70% for EN24 (Fig. 7). For comparison, the simulations were also performed on nickel. The material properties used were: Young’s modulus 200 GPa, Poisson’s ratio 0.29, mass density 8900 kg/m3, electric conductivity 14.3 MS/m, and assuming the most favourable operation point for magnetostriction on nickel, i.e. 20 kA/m, relative permeability 24, and magnetostrictive constant d15 ¼4.09 nm/A. It has also to be noted that the bias magnetic flux density in nickel is B ¼ 600 mT due to magnetic saturation. Nickel is significantly more magnetostrictive than steel, thus in this case magnetostriction is the larger effect, the resulting displacement being 1.7 times the one due to the Lorentz force mechanisms. These results are summarized in Table 3. The predictions made for magnetostriction are essentially an upper limit; not only have we considered the maximum magnetostrictive constant for a given steel and the lowest measured permeability, also an implicit assumption has been made: that magnetostrictive constants are frequency independent. The magnetostriction curve of each material was measured in dc conditions, applying a static bias field and the resulting magnetostrictive constants were used for ac simulations. This assumption was made simply because assessing the frequency dependency of magnetostriction is a very complex experimental task, and in the literature there is a lack of dynamic magnetostriction properties. However, it is likely that when a dynamic magnetic field oscillating at frequencies in the order of hundreds of kilohertz is applied to a ferromagnetic material, not all the magnetic domains are able to follow the driving input, resulting in a reduction of the magnetostrictive coefficients. This hypothesis is strongly supported by the fact that magnetic permeability significantly decreases with frequency [15,16]; since permeability and magnetostriction are macroscopic effects caused by the same microscopic structures,

2.0

Displacemeent [arb.]

36

Lorentz Magnetostriction

1.5

1.0

0.5

0.0 EN32B

EN3

BO1

EN24

Ni

Fig. 7. Simulated displacements caused by the Lorentz force and magnetostriction in four steel grades and nickel. The amplitudes are not necessarily in phase. The same unit driving current oscillating at 2 MHz was used for all the simulations.

Table 3 Maximum magnetostrictive constants d15 of four steel samples in the range Hz A [6, 15] kA/m. The corresponding EMAT signal amplitudes (experimental), for the wave generation process only, are also shown. The last column displays the percentage ratio of the displacement caused by magnetostriction against the one due to the Lorentz force as predicted by the FE model for f¼ 2 MHz. Data on nickel are also shown for reference.

Fig. 6. Axisymmetric FE model of spiral coil EMAT. The displacement in the r direction generated by the transducer is represented by the colour plot. The dynamic magnetic field produced by the coil is represented by the contour lines.

Material

d15 [nm/A]

Exp. Amp. [OV]

MS/LOR

EN32B EN3 BO1 EN24 Nickel

1.30 1.44 1.23 1.71 4.09

0.411 0.417 0.427 0.413 –

27.5% 25.7% 35.6% 70.4% 173.5%

R. Ribichini et al. / NDT&E International 45 (2012) 32–38

some scatter in the experimental data. This is mainly due to the contribution of magnetostriction, together with experimental uncertainties in the measurement of magnetic flux density B and of the driving current I, which were quantified to 73–6% uncertainty of the signal amplitudes. From a practical point of view, since the measured amplitudes on different kinds of steel are similar, it is possible to use the same EMAT probe on a wide range of grades. Large amplitude variations have been observed in the field while inspecting steel components. Such variations are probably due to the presence of highly magnetostrictive oxide layers. In those cases, at the frequencies and permeabilities considered in this study, the transduction is mostly confined in the oxide layer, and magnetostriction is the dominant mechanism, as in the case of nickel, significantly increasing the overall signal level. It can be concluded that normal bias field EMATs do not show large variations in the performance when operating on steel with a range of different material properties, except when a highly magnetostrictive oxide layer is present.

i.e. magnetic domains, it is likely that the value of d15 used in our computation is overestimated. This is experimentally hinted at by the fact that there is no correlation between the magnetostrictive constants measured and the EMAT wave amplitudes.

4. Discussion

No rmalize d Amp l itude [a rb.]

The numerical and experimental results lead to the conclusion that the Lorentz force mechanism is the dominant one in steel, while magnetostriction plays a less significant role. This conclusion can be interpreted via the physics of the two transduction mechanisms. As long as the eddy current penetration depth is much smaller than the acoustic wavelength, it is found by integrating Eq. (1) that the total Lorentz force is proportional to R the total induced current: F L pB JdA, whereas the total magnetostrictive force is proportional to the integral of the dynamic R ~ magnetic field: F MS pd15 HdA. By using an electromagnetic FE model, or analytical solutions [17], we can compute the dependencies of these quantities on electrical conductivity and magnetic permeability. The results are shown in Fig. 8, normalized on the y-axis in order to show the relative variations of the integrals with the electromagnetic properties. The Lorentz force is not very sensitive to changes in s and mr because highly conductive materials show a shielding effect: the eddy currents tend to equal and mirror the driving current, regardless of their spatial distribution, which is governed by conductivity and permeability [3,5]. For this reason, the total eddy current, and thus the Lorentz force, is relatively insensitive to conductivity and permeability changes in highly conductive materials. On the other hand, magnetostriction is highly affected by s and mr because not only ~ along the does the distribution of the dynamic magnetic field H depth of the material change, but also its amplitude. This means that the integral of the magnetic field, and thus magnetostriction, is strongly affected by the electromagnetic properties of the material. The overall conclusion is that if the Lorentz force mechanism is dominant, a small variation of signal amplitudes with conductivity and permeability is to be expected, while if magnetostriction is the main transduction mechanism, large variations in the amplitudes should be observed. The relatively small variation of signal amplitudes in the experimental data supports the argument that Lorentz force is the dominant transduction mechanism for this EMAT configuration [3,8,9], in agreement with FE predictions. A purely Lorentz force mechanism would give virtually no variation with s and mr; however, there is

5. Conclusions Electromagnetic Acoustic Transducers operate on ferromagnetic materials via two physical phenomena: the Lorentz force and magnetostriction. Previous research on bulk shear wave EMATs has established that when the magnetic bias field is parallel to the surface of the sample magnetostriction is the dominant effect, while when it is normal to the surface diverging conclusions have been drawn. Some authors stated that the Lorentz force is the main effect [3,8,9], while others [4,10] claim that magnetostriction is up to two orders of magnitude larger than the Lorentz force. Experimental tests and numerical simulations undertaken in this study indicate that the Lorentz force is the largest transduction mechanism on steel materials, regardless of the level of magnetic bias field employed, while the Lorentz force and magnetostriction are of the same order in nickel. This finding is in contradiction with relatively recent claims [4,10], but agrees with previous studies [3,8,9]. This conclusion is of practical importance because, unlike magnetostriction, the Lorentz force is relatively insensitive to the range of material properties of steels. This implies that using the same EMAT probe on various grades is possible and yields similar performance. However, signals will increase when a highly magnetostrictive oxide is present so magnetostriction becomes significant, while the performance on austenitic steels is poorer than ferritic steels because of the reduced bias magnetic field.

1.0

1.0

0.8

0.8

0.6

0.6

00.4 4

0.4 04

0.2

0.2

0.0 1

2

3

4

Conductivity, σ [MS/m]

37

5

0.0 40

70

100

130

160

Magnetic permeability, r

R R ~ continuous line) plotted against (a) electric conductivity and (b) magnetic Fig. 8. Total induced current ( JdA, dashed line) and total dynamic magnetic field ( HdA, permeability. Since the Lorentz force is proportional to the total induced current and magnetostriction is proportional to the total dynamic magnetic field, these plots show the dependency of the two transduction mechanisms on material properties. The values on the y-axis are normalized to show the relative variations with s and mr.

38

R. Ribichini et al. / NDT&E International 45 (2012) 32–38

References [1] Dobbs ER. Electromagnetic generation of ultrasonic waves. In: Mason WP, Thurston RN, editors. Physical Acoustics, vol. X. New York: Academic Press; 1973. p. 127–89. [2] Maxfield BW, Fortunko CM. The design and use of electromagnetic acoustic wave transducers (emats). Mater Eval 1983;41:1399–408. [3] Thompson RB. Physical principles of measurements with emat transducers. In: Mason WP, Thurston RN, editors. Physical Acoustics, vol. XIX. New York: Academic Press; 1990. p. 157–200. [4] Hirao M, Ogi H. EMATs for Science and Industry: noncontacting Ultrasonic Measurements. Boston: Kluwer Academic Publishers; 2003. [5] Jiles D. Introduction to Magnetism and Magnetic Materials. London: Chapman and Hall; 1998. [6] Lee EW. Magnetostriction and magnetomechanical effects. Rep Prog Phys 1955;18(1):184–229. [7] Bozorth RM. Ferromagnetism. Princeton: Van Nostrand; 1951. [8] Wilbrand A. Emus-probes for bulk waves and Rayleigh waves. Model for sound field and efficiency calculations. In: Holler P, editor. New Procedures in Nondestructive Testing. Berlin: Springer; 1983. p. 71–82. [9] Wilbrand A. Quantitative modeling and experimental analysis of the physical properties of electromagnetic-ultrasonic transducers. In: Chimenti D, Thompson D, editors. Review of Progress in Quantitative Nondestructive Evaluation, vol. 7. New York: Plenum; 1987. p. 671–8. [10] Ogi H. Field dependence of coupling efficiency between electromagnetic field and ultrasonic bulk waves. J Appl Phys 1997;82:3940–9.

[11] Bowler JR, Bowler N. Theory of four-point alternating current potential drop measurements on conductive plates, Proceedings of the Royal Society A: mathematical. Physical and Engineering Sciences 2007;463(2079):817–36. [12] Yin W, Peyton AJ, Strangwood M, Davis CL. Exploring the relationship between ferrite fraction and morphology and the electromagnetic properties of steel. J Mater Sci 2007;42(16):6854–61. [13] Thompson RB. Mechanisms of electromagnetic generation and detection of ultrasonic lamb waves in iron–nickel alloy polycrystals. J Appl Phys 1977; 48(12):4942–50. [14] Ribichini R, Cegla F, Nagy PB, Cawley P. Quantitative modeling of the transduction of electromagnetic acoustic transducers operating on ferromagnetic media, IEEE Transactions on Ultrasonics. Ferroelectr and Freq Control 2010;57(12):2808–17. [15] Rose JH, Uzal E, Moulder JC. Magnetic permeability and eddy- current measurements. In: Chimenti D, Thompson D, editors. Review of Progress in Quantitative Nondestructive Evaluation, vol. 14. New York: Plenum; 1995. p. 315–22. [16] Bowler N. Frequency-dependence of relative permeability in steel. In: Chimenti D, Thompson D, editors. Review of Progress in Quantitative Nondestructive Evaluation, vol. 25. Brunswick, Maine, USA: AIP; 2006. p. 1269–76. [17] Chari M, Reece P. Magnetic field distribution in solid metallic structures in the vicinity of current carrying conductors, and associated eddy-current losses, power apparatus and systems. IEEE Trans Power App Syst 1974; 93(1):45–56.