Harmonic Distortion Measurements for Non-Destructive Evaluation of Steel Strips

TL;DR

This paper reviews recent advances in using harmonic distortion measurements for non-destructive evaluation of steel materials.

Contribution

The paper provides a state-of-the-art review focusing on recent simulation and inversion method advancements in the field.

Findings

Harmonic distortion measurement is a widespread technique for evaluating ferromagnetic materials.

Recent progress in simulation and data processing can enhance the accuracy of these measurements.

Sensor design and signal interpretation remain empirical challenges in the field.

Abstract

There is an increasing demand for better and more accurate non-destructive evaluation and testing techniques for ferromagnetic materials, in particular, steel structures, which span diverse industrial sectors. Among the available methods, the measurement of the harmonic content and the corresponding distortion due to the material non-linearity is one of the most widespread approaches, with significant implications in the steel industry. Although these approaches have been in use for several decades with considerable success, crucial aspects, including the sensor design, the fine-tuning of the method, and, most importantly, the processing and interpretation of the acquired signals, remain, to a certain degree, empirical. Meanwhile, significant progress in simulation and data processing techniques has been reported in recent years, whose potential in the further development of the existing measuring technology cannot be overstated. In view of this rapidly changing landscape, the present review article aims to provide an overview of the state-of-the-art of the harmonic distortion measurement technique from the perspective of the recent advances in simulation and inversion methods.

Figures

Full text

1. Introduction

Non-destructive material evaluation methods serve two basic purposes: i. retrieving characteristic material properties, and ii. assessing the degree of degradation with respect to a nominal state in a non-destructive way, i.e., without altering the structural integrity. They are introduced to replace, or at least minimise, expensive destructive tests, i.e., experiments involving local piece deformations (as in the case of the hardness test) or the removal of small specimens for further laboratory examination. It is well-understood that the latter family of methods is subjected to a number of practical and economical restrictions, which exclude the possibility of in situ examination.

When speaking of ferromagnetic materials, a key property that non-destructive measurements rely on is the strong correlation of the magnetic field with the material microstructure. Since magnetic field is an easily experimentally accessible quantity, one can understand that its measurement in the vicinity of the piece can deliver rich information about the state of the material. There are different available magnetic measurements, each one focussing one or more features of the magnetisation curve, depending on the material characteristics that one needs to reveal. Hristoforou, Ktena, Vourna, et al. used differential permeability measurements for the nondestructive assessment of residual stresses [1,2,3]. Residual stresses are also well correlated with the incremental permeability, as was demonstrated by Uchimotto, Duscharne et al. [4,5,6,7]. The magnetic Barkhausen noise (MBN) proves to be a very sensitive tool for hardening detection and characterisation, as is shown in a number of works [8,9,10,11,12,13,14,15]. Martinez-de-Guerrenu used a number of different identifiers derived from hysteresis curves, such as the coercive field, the remanence, and the hysteresis losses, to monitor the recovery process during annealing [16,17,18].

Many of the cited works concern laboratory measurements involving specially prepared specimens, which encounter a number of serious limitations when coming to in situ non-destructive measurements. An important step towards a systematic and robust exploitation of different magnetic identifiers suitable for outdoor/industrial environments was taken by Dobmann and coworkers during the 1980s. Their pioneering contributions resulted in the introduction of the 3MA (Micromagnetic-, Multiparameter-, Microstructure-, and Stress-Analysis) measuring concept [19,20,21], which was subsequently turned into a commercial device [22]. The 3MA system combines three different magnetic measurements, namely, harmonic distortion, incremental permeability, and MBN with a classical eddy-current measurement, which are then supplied as inputs to a dedicated regressor to provide predictions for different properties like yield and tensile strength, hardness, etc. Similar commercial solutions, albeit ones focused on specific measurements, are now offered by other companies such as the PropertyMon [23] or the IMPOC system [24].

The usefulness of such measurements has been demonstrated after years of operation in more than one steel production line and a number of associated research projects [25,26,27,28,29,30,31]. Despite the very satisfactory results so far, the rationale behind these technologies remains very pragmatic, heavily relying on the existence of calibration curves. There is a certain need for more quantitative results that will benefit from the recent advances in the domain of numerical simulation and signal processing.

In the present review article, we are interested in the particular case of the harmonic distortion measurement, which we try to examine from the above-mentioned perspective. Starting with the general mathematical statement of the characterisation problem and the definition of the inputs and outputs, as well as the involved mathematical operators, we proceed to a presentation of the underlying physical concepts and the typical experimental setups used for these kinds of measurements. This introductory section is complemented by a brief presentation of the three commercial implementations mentioned above: the 3MA, the PropertyMon, and its predecessor, the HACOM system.

In the second section, the focus shifts to the study of the material law. Different parametric models for the calculation of the magnetic hysteresis are discussed and classified according the the number of input parameters. This piece of information allows one to have an image of the problem dimensionality, which is an important aspect when we try to construct inverse models. Afterwards, the spectral decomposition of the operator is discussed, in an effort to give a qualitative understanding of how the hysteresis parameters are linked to the harmonic content. At the end of this section, the harmonic distortion variable is introduced, which is a quantity that provides an alternative figure of merit for magnetic materials.

The third section is dedicated to the processing of the measured signals. To facilitate understanding, some simple analytical and semi-analytical models for the calculation of the field response with the three considered setups are examined, followed by a detailed discussion of the model inversion.

The main ideas discussed in the article are numerically illustrated by means of a representative example. At the end of the paper, a short account of the correlation with the microstructure and the mechanical properties that are the main quantities of interest for real-world applications are given.

2. The Harmonic Distortion Measurement

2.1. Formal Statement of the Characterisation Problem

A ferromagnetic material is characterised by its magnetisation, which is a function of the externally applied magnetic field $[eqn]$ and the material microstructure, the latter being mathematically expressed by a set of appropriate parameters $[eqn]$ :

$$[eqn]$$

The dimension of $[eqn]$ and the physical meaning of the involved parameters depend upon the model used to describe $[eqn]$ . For this reason, $[eqn]$ is intentionally kept unspecified until the introduction of the hysteresis’ explicit mathematical expression. All we need to know at this level is that $[eqn]$ provides a complete characterisation of the material.

The conceptual image of the magnetic characterisation measurement is schematically depicted in Figure 1. Our primary task is to determine $[eqn]$ in terms of a set of suitable electromagnetic observables, herein referred to as $[eqn]$ . In other words, we seek to establish the relation $[eqn]$ . The measurement process itself can be understood as the concatenation of two operators: the material operator $[eqn]$ , which translates the material parameters $[eqn]$ to the magnetisation function (which is nothing more than the numerical evaluation of (1) under the given operational conditions of the setup),

$$[eqn]$$

and the measurement operator, which associates the material magnetisation state to the observables,

$$[eqn]$$

$[eqn]$ is a linear operator that only depends on the setup geometry, and its expression can be calculated by solving the underlying electromagnetic problem for the magnetic field around the piece and thereupon calculating its convolution with the transfer function of the sensor. This is the task of the electromagnetic solver, which is, in a general case, a numerical solver such as the finite element method (FEM) or the finite difference/integration method (FDM/FIT). In the case of simple geometries, this calculation can be also carried out analytically. Some standard approximate analytical solutions for three typical setups are given in a later section for completeness.

The determination of the material operators $[eqn]$ can be thus formally obtained by calculating the inverse operator

$$[eqn]$$

In the particular case of harmonic measurements, one usually does not proceed to the end of the chain by inverting $[eqn]$ but is usually satisfied with the information hidden in the spectral information of the measurement; otherwise stated, one works with the ratio of the higher harmonics to the fundamental one.

We assume that a sinusoidal excitation of radial frequency $[eqn]$ is exerted upon the specimen. There are two main excitation modes: either by controlling the magnetic field or the mean magnetic flux that passes through the specimen cross-section (in the case of uniaxial configurations). Whatever the excitation mode, all remaining (non-imposed) magnetic variables will contain the basic frequency with all the higher harmonics, a direct result of the field interaction with the non-linear material operator. Thereupon, the resulting expression for the magnetic field reads

$$[eqn]$$

with $[eqn]$ . The fact that only odd harmonics appear in (5) is a well-known experimental fact and reflects the point symmetry property of the hysteresis. We shall comme to this point in the following section. Note that, in the special case of an H-imposed excitation, $[eqn]$ .

The corresponding time signals for the magnetisation and the measurements will admit similar decompositions, namely,

$$[eqn]$$

and

$$[eqn]$$

Specific attention should be drawn here that the ratios between the harmonics of the measured quantities are not identical to those of the material operator and, hence, the latter should not be confused with the former. In fact, for a rigorous extraction of the material information, one needs to remove the spectral deformation introduced by the transfer function of the measurement disposition. This de-convolution operation is the task of the inverse electromagnetic solver:

$$[eqn]$$ $$[eqn]$$

Coming back to Figure 1, one can also discern two side branches. The lower-right branch relates the mechanical parameters of the material with the experimental observables $[eqn]$ . Lacking a rigorous mathematical framework, the usual way to perform this prediction is via regression, which is denoted in the figure via the action of the $[eqn]$ operator. The construction of the appropriate regression operator, however, is not a simple task since it necessitates a number of (destructive) calibration measurements, which are not always available.

The lower-left branch associates the input model parameters $[eqn]$ with actual microstructural features like the grain size and the dislocation density. Also, here, a complete theory is missing at the moment. Instead, a number of empirical formulas that establish certain links between the microstructure and basic hysteresis characteristics exist in the literature, which also require extensive calibration. Following the complete chain from the lower-left to the lower-right permits (at least theoretically) the establishment of an indirect assessment of the destructively determined mechanical parameters with the underlying microstructure, which is of strong interest in many industrial applications. More on this point is discussed in the last part of this article.

2.2. Measurement Principle, Basic Setups and Experimental Observables

Magnetic methods for the characterisation of ferromagnetic materials exploit the interaction between an inductor (or a system of these) and the tested specimen, which is brought in close interaction. The working principle is schematically explained in Figure 2.

According to Ampère’s equation, the inductor will induce a magnetic field in the surrounding area, including the piece itself:

$$[eqn]$$

where $[eqn]$ is the current density in the inductor windings and $[eqn]$ is the induced source field.

The piece will interact with $[eqn]$ , with the resulting total magnetic field after the interaction $[eqn]$ being given by the relation

$$[eqn]$$

where $[eqn]$ is the electrical conductivity of the medium and $[eqn]$ stands for the magnetic induction, which is locally related with the magnetic field via the magnetic constitutive relation:

$$[eqn]$$

$[eqn]$ being the magnetic permeability of the free space and $[eqn]$ the magnetisation of the medium. The latter contains all the physics of the material and, thereupon, its microstructural information. In non-magnetic regions like the air $[eqn]$ , and (12) reduces to a simple linear relation.

The solution to (11), together with the constitutive relations in (12) and the boundary conditions at the piece interfaces, will provide the sought $[eqn]$ and $[eqn]$ values inside the piece and the air. The solution of the Maxwell equations in the presence of a ferromagnetic medium is a difficult problem, which is the subject of the whole research area of computational electromagnetism. For this reason, we cannot proceed to a detailed discussion in the limited space of this article. A number of representative articles are given instead (recall once more that the literature is huge and the following list does not claim by any means exhaustiveness): [32,33,34,35,36,37,38,39,40,41,42,43,44]. In Section 4.2, a semi-analytical solution is sketched for the case of an axisymmetric problem, which is meant to serve as example solution for these kinds of problems.

From the above discussion, it is inferred that the total field in the vicinity of the material contains indirect information about the material microstructure through the employment of $[eqn]$ . This information can be retrieved by reading the field response at some point close to the piece. This can be achieved in two ways, either by directly measuring the magnetic field using a magnetic sensor like the Hall effect, giant magnetoresistance/magnetoimpedance (GMR/GMI), flux-gate, etc., or by measuring the time derivative of the magnetic flux density using the exciting or a secondary inductor. In the second case, our observable will be the electromotive force (EMF) at the inductor’s terminals $[eqn]$ , which is given by Faraday’s law:

$$[eqn]$$

where N is the number of turns and $[eqn]$ is the magnetic flux crossing the inductor cross-section, namely,

$$[eqn]$$

The typical setups for conducting this type of measurement are illustrated in Figure 3. In the first-from-the-left configuration, a homogeneous magnetic field is induced in the tested specimen by means of a long solenoid, and the flux variation is measured using a secondary coaxial to the solenoid coil. This is the simplest configuration, which suffers, however, from a major drawback: its application is restricted to narrow specimens (bars or ribbons) and, hence, this can be considered as a destructive technique. A very popular variation of the solenoid setup is the Epstein frame [45].

The second configuration involves a number of induction coils with a specific phase relation between them, which serves in guiding the magnetic field lines inside the tested specimen. It is suitable for testing planar specimens of arbitrary lateral dimensions, which is the major advantage of this method. In cases where only one side is accessible, the method is still applicable, yet less efficient, since the absence of flux concentrators (ferrite cores) and the less-favorable coil arrangement result in a significant number of stray fields. These types of arrangements are also referred to as open-circuit setups. An industrial realisation can be found under the name of the HACOM and its successor, the ProperyMon device [23].

The third configuration involves a closed magnetic circuit in the form of a ferromagnetic or ferrite yoke. For an even better magnetic coupling to the tested piece, a second yoke can be added to the opposite surface. This variation is known as a single-sheet (or single-strip, in the case of narrow specimens) tester [45]. It is the configuration that implements the 3MA device, another popular industrial device for the magnetic testing of ferromagnetic pieces [22].

2.3. Industrial Realisations: 3MA, HACOM, and PropertyMon Systems

In the previous paragraph, mention was made to three popular commercial devices that have been successfully used for the quality control of the steel strips in production lines via the in-line monitoring of the magnetic properties: the 3MA, the HACOM, and the PropertyMon system. All three devices perform harmonic distortion measurements (alongside others) and are capable of predicting properties like the yield $[eqn]$ and tensile strength $[eqn]$ , using integrated multi-parameter regression tools.

The 3MA device was originally based on the contributions of Dobmann et al. during the 1980s and has undergone a number of further developments/improvements thereafter; it is currently commercialised by the IZFP institute in Saarbrücken. It offers four types of magnetic measurements, namely, harmonic analysis, magnetic Barkhausen noise (MBN), incremental permeability, and multi-frequency eddy-current measurements. Its setup involves a closed magnetic circuit approach, combined with secondary coils for performing the inductive measurements and a magnetic sensor for the acquisition of the tangential magnetic field component. It is delivered with an integrated software module, which permits data analysis and regression tools for the prediction of mechanical properties like $[eqn]$ , hardness, etc. An example of 3MA installation for the monitoring of steel quality in a production line is given in Figure 4.

The HACOM instrument is an industrial realisation of the second configuration shown in Figure 3, involving an intelligent arrangement of air-cored coils. Similarly, with the 3MA, it is combined with a dedicated software for performing acquisition, data analytics, and regression for the estimation of the key mechanical/technological properties of the tested specimen.

PropertyMon can be considered as the successor of HACOM (which is no longer commercially available), based on the same principle with its predecessor, yet offering essential improvements in terms of hardware and software tools. PropertyMon is commercialised by Primetals Technologies Austria GmbH, in Linz, Austria. A typical plant integration of the product is given in Figure 4.

3. Magnetic Properties

3.1. Parametric Models for the Description of Magnetic Hysteresis

The magnetisation function $[eqn]$ has been introduced in (1) without giving details of its form or its properties. In this section, we are interested in specifying the properties of the magnetisation function and providing approximate expressions for its calculation.

In its general form, $[eqn]$ is a non-linear, hysteretic, and anisotropic function. Since, however, working with a fully vectorised version of the magnetisation is complicated (and, for many practical situations, unnecessary), it is much more convenient to work with the special case of scalar hysteresis, i.e., the magnetic law only applies to the amplitudes of the magnetic variables:

$$[eqn]$$

with $[eqn]$ and $[eqn]$ . This simplification is a good approximation in situations where the magnetic field is predominantly in a specific direction. This is certainly the case in single-strip configurations and applies equally well in open-circuit setups such as the PropertyMon, where the magnetic field is intentionally guided parallelly to the specimen surface.

There is a plethora of magnetic hysteresis models in the literature, such as the Preisach model [46], the Jiles–Atherton model [47,48], the Stoner–Wohlfarth model [46], the Hauser (or energetic) model [49], the Mel’gui model [50], and the dry-friction model [51], just to mention some of the most well-known among them. One may also mention the Rayleigh model [52], but, since this model is only applicable in the weak-field region, it has been deliberately omitted from the above list. It can be, however, used as a constitute of hybrid formulations for the improvement of the calculations at the low-field domain, where most of the above models behave rather poorly [53].

An important question that has to be answered when dealing with characterisation is the minimum number of parameters that one needs to capture the basic hysteresis features. The larger the number of parameters the model relies on, the better the fit with the experimental data, but the more difficult its identification becomes (recall that the parameter estimations are achieved via the solution of an inverse problem).

The absolute minimum of parameters that one needs, when dealing with a non-linear model, is two. These can be the slope at zero field $[eqn]$ (initial susceptibility) and the high-field asymptote (magnetisation at saturation) $[eqn]$ , or the coercive field $[eqn]$ and the corresponding susceptibility value $[eqn]$ . Both approximations provide a partial description of the full hysteresis curve, which can be useful when we are only interested in a particular subdomain. In particular, a parametrisation based on $[eqn]$ describes the high-field regime well, which is dominated by reversible domain rotations, whereas the second parametrisation based on the local information around the coercive field $[eqn]$ is useful when working with step hysteresis loops or when we are primarily interested in the behaviour around the coercive point. This approximation is also known in mechanics as the backlash model, and this term is also used herein. The two approximations are illustrated in Figure 5(left) for the case of dual-phase (DP) steel. The measured maximum curve of the material is also given for comparison. For this particular illustration, the Langevin model has been used for the high-field approximation, but other models can perform equally well. The explicit expression of the Langevin expression reads

$$[eqn]$$

where the a parameter is related to the initial susceptibility via the relation $[eqn]$ .

In order to be able to reconstruct the closed loop, one needs at least one additional parameter that provides the information of the loop width, expressed by the coercive field. A simple expression for the hysteresis loop was proposed by Ponomarev in the 1980s [54]:

$$[eqn]$$

where

$$[eqn]$$

$[eqn]$ is the saturation magnetisation, $[eqn]$ the coercivity, $[eqn]$ is a parameter describing the squareness of the hysteresis loop, and $[eqn]$ stands for the magnetic field amplitude at the reversal point. The negative/positive sign applies to the ascending/descending branch.

The $[eqn]$ parameter can be expressed by the more familiar remanent magnetisation $[eqn]$ by setting $[eqn]$ in (17) and solving for $[eqn]$ , which yields

$$[eqn]$$

In addition to its striking simplicity, the Ponomarev model succeeds in delivering a fairly good description of the experimental data for a large family of soft industrial steels, where the irreversible magnetisation mechanisms due to wall movement are predominant. The comparison of the model prediction with our example DP steel is shown in Figure 5(center). It is worth noting that the given model is also interesting because, although originally introduced as a heuristic expression, it can be derived from first principles under certain assumptions, as shown in [55]. In the same paper, it is also shown that the skewness parameter $[eqn]$ is linked to the material disorder. The derivation of (17) is given in Appendix A for completeness.

A generalisation of Ponomarev’s expression was proposed by Mel’gui in [50], which introduces a supplementary empirical term in order to improve the model accuracy in the high-field region, where domain rotational mechanisms become important. This improvement comes, however, at the cost of two additional parameters. The Mel’gui expression reads

$$[eqn]$$

where $[eqn]$ is the Ponomarev expression given in (17). The two additional parameters introduced in Melgui’s experssion are the initial susceptibility $[eqn]$ and the weighting coefficient k, defined as

$$[eqn]$$

$[eqn]$ stands for the coercive field, $[eqn]$ is the remanence, and $[eqn]$ is defined as the virgin curve value that corresponds to the coercive field, namely, $[eqn]$ . The model performance is illustrated in Figure 5(right). We can easily observe the better fit for the higher field values.

A summary of the discussed models and their tuning parameters is given Table 1.

3.2. Spectral Analysis of the Hysteresis Operator

After having introduced the parametric model for the description of the hysteresis operator, the next question that we need to address is how the model parameters are associated to the spectrum of the magnetisation function (or its derivative, the differential susceptibility), which is the property of our concern.

In order to proceed with this discussion, we need to select a working hysteresis model. The model of choice has to be as simple as possible but, an the same time, representative of the main material features. Ponomarev’s model makes a good candidate given its simple form, the restricted number of parameters, and the very good approximation of real hysteretic curves that it provides.

Although, the magnetisation spectrum can be easily computed numerically, it is instructive to try to derive approximate analytic relations, involving explicit expressions of the material parameters $[eqn]$ , and $[eqn]$ . This can be achieved by Taylor series development, which yields the following expression for the lowest harmonics, assuming a sinusoidal excitation of the form $[eqn]$ :

$$[eqn]$$

The derivation of (22) and the explicit expressions of the series coefficients are given in Appendix B.

The results of the previous analysis can be complemented by numerical studies with the Ponomarev model in Expression (17) for variable $[eqn]$ and $[eqn]$ and by calculating the harmonic coefficients of the resulting magnetisation signal via Fourier transform. As study material, we consider the same IF steel grade used for the studies in [55]. The material parameters are obtained by the identification of the experimental curves, following the standard procedure, and they are summarised in Table 2. We then let $[eqn]$ and $[eqn]$ vary in the interval 10–500% with respect to their nominal value. The resulting bunch of curves obtained with the considered values are demonstrated in Figure 6.

The relative variation of the first three harmonics with respect to their value obtained for the original curve is illustrated in Figure 7. From a comparison of the curves, it turns out that $[eqn]$ variations result in a stronger variation of the higher harmonics than the ones when varying $[eqn]$ . In other words, the sensitivity of the harmonic measurements with respect to the coercive field is stronger that that for the form factor $[eqn]$ . In this particular example, the curves corresponding to $[eqn]$ variation have been multiplied by a factor of 100.

3.3. The Harmonic Distortion Identifier K

Instead of directly using the harmonic amplitudes $[eqn]$ , many practitioners prefer working with two derived indicators: the coercive field $[eqn]$ and the harmonic distortion factor K, defined as the ratio of the mean of the higher harmonic amplitudes $[eqn]$ to that of the basic harmonic $[eqn]$ :

$$[eqn]$$

The advantage of using the harmonic distortion identifier is that it directly integrates the information from several harmonics, thus offering better robustness to noise. Assuming that we have estimated $[eqn]$ independently, the knowledge of K together with $[eqn]$ is theoretically sufficient to identify the normalised Ponomarev’s model (and, indirectly, the material itself).

In the previous paragraph, the dependency of the first three harmonics from the $[eqn]$ and $[eqn]$ parameters has been demonstrated both analytically and numerically. From the dimensionality of the parameters space, it turns out that one needs to measure at least two harmonics for the identification of the $[eqn]$ parameters. We can, therefore, infer that $[eqn]$ , K deliver sufficient information for the reconstruction of the normalised hysteresis curve.

On the contrary, since $[eqn]$ is canceled out as a common factor in all terms in (23), its determination requires an additional independent measurement. The simplest choice is with proceed to an eddy-current measurement, which provides us the complementary information of the initial magnetic susceptibility $[eqn]$ (under the assumption that there is zero remanence at the measurement point). Using the Langevin expression at zero-field in (16), we easily derive the relation initial susceptibility $[eqn]$ . The a parameter can be estimated by fitting the normalised Langevin curve to the high-field part of the (also normalised) Ponomarev model, which leaves us with only one unknown parameter: $[eqn]$ .

4. Measurement Processing and Inversion

4.1. Semi-Analytical Approximations for the Inverse Operator

Given the complexity of the field problem, the calculation of the transfer function $[eqn]$ is carried out, in most cases, numerically. Nevertheless, reasonable analytical (or semi-analytical) approximations can be obtained for simplified versions of the setups in Figure 3. The interest in such approximations is not only academic. Fast approximate solutions can be used to give a coarse qualitative image of the setup response, and they can be combined with minimisation algorithms for the solution of the inverse problem, which tend to be numerically cumbersome when using numerical solutions.

Stating with the solenoid, the calculation of both $[eqn]$ and $[eqn]$ is very straightforward. The variable of interest here is the voltage at the leads of the receiving coil:

$$[eqn]$$

where $[eqn]$ are the linear density of the solenoid windings (in turns/m), $[eqn]$ is the magnetic permeability of air, $[eqn]$ are the cross-section and the number of turns in the receiver coil, and I the excitation current. Equation (24) can be trivially inverted, giving

$$[eqn]$$

In the case of the closed magnetic circuit in Figure 3, a simple coarse approximation can be obtained by using the well-known magnetic circuit model. The very low magnetic reluctivity of the yoke allows us to assume that the magnetic flux is preferentially contained inside the yoke volume, with only a weak proportion of the magnetic streamlines straying away. Furthermore, since the yoke core has a very low conductivity, the eddy-current effect can be ignored, which implies that the magnetic field distribution throughout the yoke cross-section is nearly homogeneous.

At a first order, we may also assume that the same approximations also hold inside the tested specimen, though less accurately than in the yoke. We may additionally assume that the field is approximately constant along the yoke and along the specimen, which is equivalent with stating the the demagnetisation field due to the geometry boundaries is negligible. These kinds of models, where the spatial variation of the involved quantities is discarded, are classified in the literature as 0D or lumped models. The above-described assumptions of the working model are schematically depicted in Figure 8.

The application of Ampère’s law along the dashed integration path of Figure 8 yields

$$[eqn]$$

where $[eqn]$ are the magnetic field in the yoke and the specimen, respectively, $[eqn]$ are the corresponding path lengths, $[eqn]$ is the number of turns in the excitation coil, and I the feed current.

It is also reasonable to assume that the $[eqn]$ curve of the yoke is very narrow, with low coercivity $[eqn]$ and high saturation field, which allows us to approximate it with a linear magnetic material of high relative permeability $[eqn]$ . This assumption is very close to reality for many commercial yokes. One can thus write

$$[eqn]$$

Imposing the continuity of the magnetic flux density across the yoke poles, we end up with the expression

$$[eqn]$$

The solution is the common point of (28) and $[eqn]$ , which is obtained by iteration.

Our observable is the tangential magnetic field $[eqn]$ , which is measured by means of a magnetic field sensor (Hall element, GMR, etc.). Assuming a linear response function for the sensor, we can write

$$[eqn]$$

Both simplified models examined above share the common assumption that the entire magnetic flux passes through the piece. The effect of the piece conductivity has been also tacitly ignored, resulting in a homogeneous field in the piece cross-section. Should the fringing effect at the piece interfaces be non-negligible or should the conductivity play an important role, the field distribution become more complicated, which requires the use of more sophisticated approaches. Even in such situations, however, semi-analytical solutions, where the field is developed in a series, remain tangible. This is the case, for example, in the open-circuit setup shown in Figure 3 if a single pair of coils is involved.

The development of the semi-analytical solution for this case is described in [43]. Here, we simply sketch the main ideas of the approach in order to give the coarse image.

The magnetic field in this case is 2D, which makes it more convenient to work with a uniaxial magnetic potential instead of the magnetic field itself. The common choice is the magnetic vector potential, which, in the case of an axisymmetrical problem, reduces to a single-component potential along the azimuth direction $[eqn]$ , namely, $[eqn]$ , with $[eqn]$ standing for the corresponding unit vector. The magnetic potential is related to the magnetic induction via the relation

$$[eqn]$$

Assuming a harmonic excitation of the form $[eqn]$ , the resulting state equation for the nth harmonic of the magnetic potential is the diffusion equation, which in the cylindrical coordinates system ( $[eqn]$ ) reads:

$$[eqn]$$

where $[eqn]$ is the air permeability, $[eqn]$ is the Kronecker delta, and J is the current density in the coil. Note that J is zero everywhere except inside the coil. $[eqn]$ stands for the material magnetisation, and it is only nonzero inside the plate. Although it is not a priori known, it stands in the right-hand side of (31) since it acts as a magnetic source. Its value is estimated via the successive solution of (31) by requiring the fulfillment of the material law (12), according to the standard fixed-point technique. (The reader may notice here that $[eqn]$ is not a uniaxial field in the 2D case considered here, which raises questions about the adequateness of the scalar hysteresis approach adopted throughout this article. In fact, it can be shown that the normal to the plate $[eqn]$ component is much weaker that the tangential one $[eqn]$ , which makes the field solution inside the plate quasi-uniaxial. This is a peculiarity of the special coil arrangement in this setup).

Due to the odd symmetry of the field with respect to the z-axis, we set the $[eqn]$ point to be in the middle of the plate and we truncate the domain by considering only the $[eqn]$ axis. To be consistent with the problem symmetry, we also need to assume a perfect magnetic conductor condition at $[eqn]$ . In order to simplify the spectral expressions, we assume an artificial boundary at a certain radial distance from the coil $[eqn]$ . This truncation is justified by the fact that the field is rapidly decreasing as we move away from the source. The type of termination condition is irrelevant. For simplicity, we assume a perfect magnetic wall at $[eqn]$ .

The solution for the magnetic potential in the air region above the plate reads

$$[eqn]$$

where $[eqn]$ and $[eqn]$ are the series coefficients for the source term and reflection from the plate and $[eqn]$ are the first-order Bessel functions. The eigenvalues $[eqn]$ are determined by the boundary condition at the truncation surface:

$$[eqn]$$

The solution is more complicated inside the plate since, now, the right-hand side of (31) is nonzero. The series development for the homogeneous Helmholtz equations needs to be complemented with a second term that takes care of the magnetisation term. We can thus write

$$[eqn]$$

with

$$[eqn]$$

and the $[eqn]$ eigenvalues being given by

$$[eqn]$$

The $[eqn]$ coefficients stand for the special solution of (31), and they are calculated by projecting the right-hand side to the development basis:

$$[eqn]$$

where $[eqn]$ is the normalisation coefficient of the Bessel function basis. $[eqn]$ are obtained by observing (37) at $[eqn]$ , namely,

$$[eqn]$$

The $[eqn]$ , $[eqn]$ coefficients stand for the homogeneous solution, and they are determined by applying the continuity relations at the piece boundary. The final relations read

$$[eqn]$$

and

$$[eqn]$$

where R and T are Fresnel’s reflection and the transmission coefficient, respectively:

$$[eqn]$$ $$[eqn]$$

The source coefficients $[eqn]$ are obtained by solving the field problem with the coil in the area, and, therefore, they depend only upon the coil geometry and the excitation current. Their expression for the region below the coil and the plate (we are usually interested only in the field in the immediate neighbourhood of the plate interface; similar expressions that hold for the domains beside and over the coil can be readily obtained, yet their explicit expressions are not provided here for brevity—the interested reader can refer to [56]) is as follows:

$$[eqn]$$

where $[eqn]$ stand for the coil inner and outer radius, respectively, l is the coil thickness, and $[eqn]$ is the current density across its cross-section. The $[eqn]$ function stems from the integration over the coil cross-section and it has a closed-form expression in terms of Struve functions $[eqn]$ : $[eqn]$ .

The measurement in this case is either the voltage variation at the terminals of the coils or the tangential magnetic field at a specific observation position $[eqn]$ . In the second case, we can write, for the nth harmonic of the response,

$$[eqn]$$

4.2. Direct Extraction of the Coercive Field from the Magnetic Field

The direct measurement of the magnetic flux density $[eqn]$ is not possible in setups involving open or closed magnetic circuits. This inhibits the direct measurement of the coercive field $[eqn]$ . Alternative approaches for an approximate estimation of $[eqn]$ were proposed by Dobmann et al. [57,58] and also by Fillon et al. [59], where they demonstrated, either experimentally or via simulation, that $[eqn]$ is closely related to the zeros of the second time derivative of the tangential magnetic field. This point also proves to be correlated with the cross-over point of the higher magnetic field component (i.e., the signal we obtain by subtracting the basic harmonic).

A good theoretical justification of the first observation is given in [57]. Let us consider the definition relation of the differential susceptibility, which, using the chain rule, can be written as the quotient of the time derivatives of M and H:

$$[eqn]$$

The time derivative of the susceptibility reads

$$[eqn]$$

Given that the susceptibility becomes a maximum at the coercive field, and given that the magnetic field is stationary only at the turning points of the hysteresis, i.e., at the maximum of the excitations, we can safely assume that the denominator is nonzero in the region of interest. We can hence conclude that

$$[eqn]$$

We also know that, at the coercive point, the Bloch wall mobility also reaches a maximum. In the simple of a two-domain system separated by a single wall, it can be shown should be present here that the magnetisation rate is related with the velocity of the wall via the relation [60]:

$$[eqn]$$

where $[eqn]$ is the mean wall velocity and L is its length. From (48), and by assuming that the rate of the mean wall velocity is stationary at the coercive point $[eqn]$ , we can state that the coercive field is an inflection point for the magnetisation. Let $[eqn]$ be the time instance when reaching the coercive field, i.e., $[eqn]$ . Then,

$$[eqn]$$

Combined with (47), the last condition yields

$$[eqn]$$

In other words, the coercive point is an inflection point for the magnetic field as well.

Conditions (49) and (50) are naturally satisfied in case of an induction-driven system: the magnetic flux density B is the experimentally controlled variable (since the external magnetic field is the physical intensive variable that is accessible to the experimentalist, an external closed-loop system is necessary for the realisation of a B-driven system). Indeed, from the constitutive equation, we have

$$[eqn]$$

$[eqn]$ at $[eqn]$ by definition. For a sinusoidal magnetic induction, it follows that the second derivative is also zero at $[eqn]$ , and, taking (47) into account, we obtain

$$[eqn]$$

The first term is nonzero, which naturally yields (50). It is clear that Condition (50) cannot be fulfilled when the magnetic field is externally imposed (e.g., for a solenoid excitation). In that case, we must conclude that (49) cannot hold and that the mean wall velocity is not at its exact maximum there. In other cases between the two extremes of pure sinusoidal B or H, both fields dispose higher harmonics, and we can thus conclude that (49) and (50) are approximately fulfilled [57,59]. These are the cases where there is an M-H coupling due to the presence of a demagnetising field caused by the piece boundaries.

The second approximation for the coercive field is based on the location of the zeros of the higher harmonic components. It is known from theory that the discontinuity of the normal magnetisation component at the piece interfaces acts as an effective magnetic surface-charge density:

$$[eqn]$$

at the coercive point $[eqn]$ , which implies that $[eqn]$ . Thereupon, we can conclude that the demagnetisation field, which is produced by this surface charge, is approximately zero. The total magnetic field inside the specimen is the sum of two terms: the externally imposed magnetic field $[eqn]$ and the demagnetisation field $[eqn]$ produced by the piece discontinuities,

$$[eqn]$$

Assuming an ideal current source excitation, the external field $[eqn]$ is proportional to the current and, hence, only consists of the basic harmonic. $[eqn]$ , on the contrary, contains higher harmonic components, which are approximately zero at the coercive point. Therefore, we expect that the higher-harmonic part of H, comprising solely the $[eqn]$ contribution, will be nearly zero at the coercive point. In the case of voltage excitation, the coil current will also contain a small portion of higher harmonics, and $[eqn]$ will also comprise higher-harmonic components, which makes the location of $[eqn]$ less precise.

An illustration of the $[eqn]$ , as well as the higher-harmonic part as a function of time, can be found in [59] and is reproduced in Figure 9. In the same plot, we also draw the variation of the higher-harmonic component of the magnetic field. The numerical results have been obtained by using the 0D approximation of the magnetic circuit presented in Section 4.2 for the closed magnetic circuit problem. One can readily observe the closeness of the two points. The numerical rresults were produced using the Ponomarev model for the sample with $[eqn]$ kA/m, $[eqn]$ kA/m, and $[eqn]$ kA/m. For more details of the considered setup, the reader is referred to the original paper [59].

Figure 10 demonstrates the correlation of the coercive field approximations obtained by the two methods with the real value of the coercive field. Both estimators, i.e., the second-derivative zeros $[eqn]$ and the cross-over points of the higher harmonics $[eqn]$ , demonstrate a good correlation and, hence, can be used for the prediction of $[eqn]$ from tangential field measurements after calibration.

4.3. Numerical and Statistical Analysis of Forward and Inverse MF

Operator

The study of the correlation between physical observable quantities, model parameters, and material properties may provide a deeper physical insight into the model behaviour over a broad set of configurations. Nevertheless, data from experiments are often scarce, which limits the ability to draw meaningful conclusions on problems involving variability in the model factors to be evaluated (e.g., optimisation of a design, assessment of performances). The employment of numerical solvers has been demonstrated to enhance the analysis stage, based on experimental data, in terms of robustness, efficiency, and cost saving. Consequently, the integration of physics-based models with statistical methods has emerged as a prevailing methodology for addressing certain material characterisation problems.

In this context, the numerical modelling of the $[eqn]$ operator is pivotal in conducting advanced statistical studies on configuration that are pertinent to a specific problem. This includes highlighting potential unobserved or underestimated correlations between material parameters ( $[eqn]$ ), setup parameters, and observables ( $[eqn]$ ). As illustrated in Figure 1, the operators in (2) and (3) are essential in this regard if they can be numerically computed, enabling working with the constitution of datasets of simulations. Toward this end, the numerical counterpart of the total transfer $[eqn]$ operator can be modelled as

$$[eqn]$$

where, from now on, we refer to as the numerical operator. $[eqn]$ can then be utilised as an online operator to address optimisation tasks or to generate a synthetic dataset composed of finite number (e.g., N) of samples $[eqn]$ to accelerate advanced numerical studies. In the latter case, based on a given dataset $[eqn]$ , the numerical model itself is replaced by a suitable machine learning (ML) model (also called metamodel or surrogate model), based on a given dataset, with the aim of avoiding further calls to the physics-based model. Notwithstanding the reliance of $[eqn]$ rely on a physics-based or ML model, its potential for substantial exploitation can be employed to address optimisation problems, with the objective being the estimation of an unknown set of material parameters $[eqn]$ based on some observed experimental measures $[eqn]$ via the minimisation of a given cost function, for example,

$$[eqn]$$

Furthermore, the application of the $[eqn]$ operator can be exploited for the computation of statistical indices, aiming to retrieve the impact of the hysteresis model parameter on the observables. That is, the numerical counterpart of the $[eqn]$ operator replaces the needed for the measures to come from the total transfer for computing the Global Sensitivity Analysis (GSA) indices such as the first-order Sobol’ ( $[eqn]$ ), the $[eqn]$ -sensitivity measures ( $[eqn]$ ), etc. [61].

As stated in (56), the inversion of $[eqn]$ values based on the candidate observable $[eqn]$ can be cast as an optimisation problem to be minimised iteratively. The inversion problem can be alternatively formulated as an inverse problem to be solved via ML methods. That is, one can define the numerical counterpart of the equivalent of the inverse total transfer $[eqn]$ operator as

$$[eqn]$$

where the notation $[eqn]$ indicates the explicit mapping performed by the operator, without distinguishing between the separated inversion of the $[eqn]$ or $[eqn]$ , as outlined in (4). The selection of the most suitable ML schema is contingent upon the volume of data and the available computational resources, in addition to the anticipated outcomes. In principle, kernel machine-based algorithms, such as support vector machines, Gaussian process regressors, and kernel ridge regression, are suitable for addressing these tasks, in conjunction with deep neural network approaches.

5. Numerical Example

5.1. Harmonic Distortion Measurements in an Open-Circuit Configuration

The theoretical discussion of the previous sections can be better understood by a fully worked-out numerical example. The considered configuration is a basic open-circuit setup, based on the same working principle as that of the HACOM and the PropertyMon systems (both HACOM and PropertyMon systems are realisations of the open-circuit setup; the main differences from this example are the number and geometry of the inductors).

The considered configuration is shown in Figure 11. The details of the geometry are given in [62] and are not repeated here for the sake of brevity. The studied material is a typical structural steel that can be practically considered as isotropic. The hysteresis model used in the study is the five-parameter Mel’gui model presented above. The hysteresis model is identified using the measured major hysteresis curve in a specialised characterisation experiment according to the standard procedure. A detailed description of the latter can be found in the same article [62].

The two inductors are fed with a sinusoidal voltage with $[eqn]$ phase difference in order to enhance the tangential magnetic field component that passes through the plate and thereupon amplifies the material effect. Two excitation frequencies have been examined: 10 Hz and 20 Hz. Note that the frequency must be kept low in order to minimise the eddy-current effects, which hinder the magnetic field penetration in the interior of the specimen.

The field problem is solved numerically using the finite integration technique (FIT). A major difficulty when treating material with hysteresis consists of the multi-branch behaviour of the material curve at the reversal points of the feed current, which may destabilise the numerical solver. A further challenge stems from the low-field inner loops, for which most parametric models (including the Mel’gui model considered herein) behave poorly. Both issues can be remedied by considering the hysteresis law in an augmented parametric space, as explained in [53,62].

Figure 12 shows a comparison between the measured and simulated field response for a given position of the field sensor. In the left plot, there are the measured signals of the tangential magnetic field for four different excitation voltages. It is easily observed that the signal deviates from the sinusoidal form with increasing field intensity. The corresponding spectra obtained for the maximum level of excitation at both frequencies are shown in the right-hand plot.

It is interesting to note the fast decrease in the amplitude with increasing harmonic index. This means that the essential (and exploitable) information must be sought in the first harmonics since higher-harmonics are expected to be under the noise threshold. Obviously, the higher the excitation signal, the better the signal-to-noise ratio that should be attained; however, the aforementioned steep decrease should be kept in mind.

5.2. Correlation with Hysteresis Parameters

Once the numerical model of the problem is established and validated, we can turn our attention to the relation of the material parameters with the measurement observables. In this sense, we exploit simulations based on a design of experiment (DoE), accounting for the measurement setup parametrisation in terms of excitation frequency, voltage excitations, and the Mel’gui model parameters. The variations in the Mel’gui parameters used for the establishment of the DoE are described in Table 3.

The main purpose of the analysis is to verify how Mel’gui parameters correlate with the harmonic amplitudes for the different measurement frequencies and the different excitation levels. Toward this end, the DoE given in Table 3 is scanned to generate a dataset composed of 766 samples containing a time signal of 100 timesteps, associated with the tangential and normal field components. By applying the Fourier transform to the output signals, the measurement spectrum is obtained, from which only the first three harmonics having the highest amplitudes (and, thus, being less impacted by noise) are considered.

For the formal analysis, the Pearson’s and Spearman’s correlation coefficients have been employed. The former one postulates an ‘almost’ linear dependence between analysed variables, whereas the latter exploits the monotonicity in the analysed factors [63]. The Pearson’s correlation coefficient is defined as

$$[eqn]$$

with the subscript i indicating the i-th sample and $[eqn]$ and $[eqn]$ referring to the mean values of x and y, respectively. The Spearman’s correlation coefficient is defined as

$$[eqn]$$

where $[eqn]$ and $[eqn]$ represent the rank variables. In this work, the samples $[eqn]$ , $[eqn]$ and the stochastic variables X and Y are defined based on the Mel’gui model factors and the observable, which are considered such that $[eqn]$ $[eqn]$ , V is the excitation voltage.

Toward this end, we discussed the results obtained from the point of view of correlations and qualitative plots in order to highlight the most meaningful outcomes observed. In Figure 13, we analysed the Spearman’s and Pearson’s correlations between Mel’gui model factors and the amplitude of the the spectrum at the first three harmonics. Overall, we notice that Pearson’s and Spearman’s correlation indexes agree with each other, regardless of the components considered, as shown in the first columns of each plot. Furthermore, regardless of the harmonic considered, $[eqn]$ and $[eqn]$ are the two factors that correlate more with the amplitude of the spectral components considered, whereas the weakest correlation across the harmonics is always between $[eqn]$ and $[eqn]$ .

Furthermore, it has been observed that $[eqn]$ does not correlate with the amplitude of the spectrum for the investigated frequencies and, thus, could be likely discharged from the model parameters, with minor changes in the measured quantities. It has been also observed that one can take advantage of the correlation analysis of Mel’gui parameters to choose the most-suitable harmonic contents to be analysed in order to possibly retrieve the Mel’gui model parameters when the inversion of the hysteresis parameters is targeted [64].

5.3. Sensitivity Analysis of Hysteresis Parameters Based on Harmonic Distortion Measurements

Sensitivity analysis is a statistical methodology that aims to assess the ‘sensitivity’ of one or multiple observables against the variation in a set of factors that are supposed to have an impact on these observables. The sensitivity analysis estimation enables one to identify and/or rank how a specific parameter leads to changes in the observables. In statistics, multiple methods have been developed in order to perform sensitivity analysis studies; the most advanced one requires one to establish a meaningful statistical set in order to draw meaningful and robust conclusions [61]. For the sake of brevity, it is worth mentioning that, based on the state-of-the-art literature, we consider in the following the use of the Sobol’ and $[eqn]$ -sensitivity measurement GSA methods, which are well-established in the applied statistics community to assess the uncertainties in complex numerical models [61]. One can show that the Sobol’ indices’ calculation, in the case of the estimation of the separated effects of the i-th set of a factor onto the observable—the so-called first-order Sobol’ indices—representing the effects of the i-th factor on the output, is defined as

$$[eqn]$$

where $[eqn]$ stands for the maximum value of the ML-based operator, evaluated on the samples of a given set of factors $[eqn]$ , $[eqn]$ , with $[eqn]$ representing the expectation and $[eqn]$ representing the variance, calculated without the contribution of the ith factor [65]. Sobol’ indices are suitable measurements of the sensitivity based on the variance decompositions. In certain classes of problems, as a function of the observable considered (e.g., the amplitude of specific harmonic components, the amplitude of the measured signal), variance decomposition may not the most suitable approach; it is then interesting to perform GSA based on moment-independent indices like the $[eqn]$ -sensitivity measures for the i-th factor, which are defined as

$$[eqn]$$

with $[eqn]$ referred to as the inner statistic or separation, which represents the surface of the intersection between the conditional and unconditional probability densities of the model output $[eqn]$ obtained for the ith factor value [66].

GSA indices rely on the calculation of scalar quantities, which are supposed to capture the overall variability of the signal, making them suitable for quantifying the impact of the hysteresis parameter model on the harmonic distortion measurements.

Looking at Figure 14, and comparing these results with the correlation plots shown in the previous section, one can notice that the ranking obtained via GSA is aligned with the results. Furthermore, since GSA calculation are based on the use of a metamodel, the demanding statistical calculation can be computed in a efficient way, leading to a factor ranking calculation that is less sensitive to the dataset used. On top of that, GSA enables the estimations of the uncertainties associated with the ranking too. In view of the inversion procedure, the analysis of GSA can qualitatively provide some insight on the ill-posedness of the inverse problem by analysing the ranking results.

5.4. Prediction of the Material Parameters via Direct Inversion of

MF Operator

The reliable and fast inversion of the $[eqn]$ operator is an open and ongoing research field that may further enhance the analysis capability based on measured data. Regardless of the inversion procedure used, i.e., iterative (56) or non-iterative (57) ones, one postulates the access to a set of possible solutions in order to retrieve the actual parameters from the observable. Considering non-iterative inversion, which mainly relies on the use of ML algorithms trained in a supervised framework, one can define a regressor operator $[eqn]$ , as defined in (57) based on a specific ML model. The choice of the algorithm typically depends on the number of samples available for training the ML model as well as the cardinality of the measures and the factors to be estimated. In the context of material parameter estimation, based on the Mel’gui model, the use of kernel machines, like the Gaussian process regressor, support vector machine, etc., or effective neural network architecture can be considered.

Having examined the correlations between the hysteresis parameters (inputs) and field measurements (outputs), we now shift the focus to the direct output inversion by the numerical calculation of the $[eqn]$ operator. Given the ill-posedness of this inversion operator, we restrict the study to a DoE, which consists of a Latin hypercube, with each side spanning a $[eqn]$ interval of the nominal parameter values (with the term nominal meaning the Mel’gui parameters obtained via identification). Based on this dataset, we consider a set of samples, i.e., the so-called training set, from which the ML is based, such that, in the case of a Gaussian process regressor, is given as $[eqn]$ [64]. Thus, to compute, by inversion, the model parameters $[eqn]$ (posterior mean) associated with an experimental measure $[eqn]$ are given as $[eqn]$ , with $[eqn]$ representing a suitable covariance matrix. The uncertainty associated with the predictions can be straightforwardly computed as $[eqn]$ . The inversion results for the hysteresis parameters are graphically depicted in Figure 15. Each point in these plots stands for the inverted input (hysteresis parameter) of a synthetic field signal, randomly chosen from the DoE. The predictions of the input parameters are compared against the true input parameters used for the calculation of the synthetic signal. Note that, the closer to the diagonal that the points lie, the better the predictions. The shadowed area around the parameters highlights the uncertainties (±a standard deviation) associated with the predictions. A cross-analysis of Figure 14 and Figure 15 readily highlights the fact that the less-sensitive parameters have a harder inversion task. Such results also highlight the fact that the analysis of GSA indecision via a forward metamodel based on single harmonic contents in observable signals can be considered in order to qualitatively foresee the ‘degree’ of the ill-posedness of the inversion problem.

In Figure 16, the hysteresis curve obtained upon the inversion of the experimental data in [62] is compared against the identified hysteresis curve of the material. Despite the preliminary nature of these results, they demonstrate the potential of moving a step further than the usual calibration curves used by the majority of state-of-the-art equipment and performing direct predictions of the hysteresis parameters based on open-circuit measurements. Similar conclusions are presented in [64], where the numerical inversion of impulse response measurements have been used to retrieve the material parameters.

6. Correlation with the Microstructure and the Mechanical Properties

Upto this point, the ferromagnetic materials have been characterised via the parameters of the hysteresis models (i.e., its magnetic properties), which provide only an indirect link with the material microstructure. For a nondestructive assessment of the latter, however, a more direct link is necessary.

Still, the dependence of the magnetic properties on the microstructure is an extremely complex problem, which has been a subject of intensive study since the 1940s, with the seminal contributions of Kersten and others, with a general quantitative relation still missing. Nevertheless, notable progress has been made in the understanding of the main mechanisms that influence the wall movements, which has resulted in several useful empirical relations. Among the most-important features that have an impact on the magnetic and mechanical properties, one can mention the solid-solution proportion, the dislocation density, the grain size, the precipitates and/or second phase domains, and the texture.

Dislocations act as impediments to Bloch wall movements. The majority of studies converge to the conclusion that the coercive field $[eqn]$ is proportional to the square root of the dislocation density $[eqn]$ [67]:

$$[eqn]$$

Grain boundaries represent two-dimensional lattice defects, which also act as domain wall impediments. There is again general agreement on the mathematical expression of this dependence, which admits the following inverse exponential form [67,68,69,70,71,72]:

$$[eqn]$$

where $[eqn]$ are constants, d stands for the mean grain size, and n is an exponent whose value lies in the interval $[eqn]$ . One must note that there is no general consensus on the exact value of the n exponent, with the majority of researchers tending to $[eqn]$ . Regardless of the n value, (59) expresses an inverse dependence on the grain size, that is, the smaller the grain, the magnetically harder the material, which is consistent with the physics.

Sablik and Dupré used the combined form of the above two relations, namely,

$$[eqn]$$

with $[eqn]$ , $[eqn]$ , and n being determined upon fitting with experimental data [71,72].

The effect of the precipitates is more complicated since the total effect depends on several factors like the type, quantity, size, distribution, and hardness of the particles. It is, therefore, harder to obtain a simple quantitative general expression. Existing attempts are based on the foreign-particle theory, which roots back to Kersten but is not considered here [73,74]. Second-phase aggregates are also expected to have a similar effect.

Finally, the existence of texture is responsible for anisotropy and, hence, can be indirectly taken into account by assuming the directional dependence of the magnetic characteristics.

The yield-strength dependence on the microstructure shows similar trends with the coercive field, which is plausible since the same lattice irregularities responsible for the wall-pinning during magnetisation also act as impediments to the lattice layer creep, which is the main mechanism of plasticity. In particular, the empirical expression that relates the yield strength with the microstructure characteristics admits the general form

$$[eqn]$$

where $[eqn]$ is for the Peierls–Nabarro stress (the tension to move dislocations at $[eqn]$ K), $[eqn]$ expresses the influence of the solid-solution atoms (i.e., substitution atoms like Si, Mn, or P or interstitial atoms like C or N), $[eqn]$ stands for the dislocation effect, $[eqn]$ is the influence of the grain side, and $[eqn]$ describes the contribution of the precipitates. $[eqn]$ and $[eqn]$ follow similar laws with the corresponding ones for the coercive field, namely,

$$[eqn]$$

and

$$[eqn]$$

known as Hall–Petch and Taylor expressions, respectively. The precipitate effect is shown to be analogous to the mean distance between precipitates.

The experimental demonstration of the dependence between material microstructure and harmonic measurements is not a simple task since one needs to perform measurement campaigns using sets of different, well-characterised microstructures. Such systematic studies can be found in a number of reports from past Research Fund of Coal and Steel (RFCS) projects [25,26,30]. An example of the variation of basic harmonic as a function of grain elongation and annealing time for different cold-rolled steel samples of DX56D-Z grade is shown in Figure 17.

The similarity between (58), (59) and (62), (59) implies a direct relation between $[eqn]$ and $[eqn]$ . This is indeed confirmed by the correlation of harmonic distortion measurements with destructive measurements of the yield strength. An example of such a correlation is given in Figure 18, where the HACOM measurements of the first and fifth harmonics for a number of samples of conventional IF are compared against destructive measurements of $[eqn]$ .

7. Conclusions

The harmonic distortion measurement of the magnetic field signal in ferromagnetic materials is a flexible non-destructive technique, which offers rich information about the material’s magnetic properties. The main features of these properties can be retrieved by a basic analysis of the spectral information, notably, the harmonic ratio and the harmonic distortion identifier, which is the main analysis approach used by available commercial systems (PropertMon, 3MA). In a second step, this information, combined with an appropriately calibrated regressor (using well-characterised reference specimens), can be applied for the estimation of mechanical and other technological material properties. A regression example is shown in Figure 18, where a clear relation between the harmonic content and the yield strength is demonstrated.

The recent advances in the simulation techniques and the rapidly transforming landscape of signal processing and inversion algorithms, including machine learning, aided by the tremendous progress in hardware architecture, pave new directions for the processing of measurement data. The results in Section 5.4 demonstrate the possibility to retrieve the hysteresis curves via the numerical inversion of harmonic distortion measurements. The quality of these predictions is expected to be further enhanced when enriched with other indicators (like the incremental permeability) in comprehensive data-fusion models, which can also benefit from the new possibilities offered by tools like artificial neural networks. For example, 3MA, already offers such models, and a lot of research is currently conducted in this direction.

New exciting perspectives involving multiscale simulations, which are capable of taking features of the grain structure directly into account, give hope for a new generation of powerful tools that will provide estimations of material features that are presently accessible only via microscopic methods [30,75].