# A Simple Magnetostatic Sensing Method for Assessing the Local Hysteresis Properties in Ferromagnetic Sheet Materials

Authors: Ivan J. Garshelis; Guillaume Crevecoeur
Journal: Journal of Sensors
 (2012)

---

## Abstract

Global hysteretic properties of electrical steels can be measured using ring or strip samples, while the assessment of the local hysteretic properties is a much more difficult task since the measurement method needs to be very sensitive. This paper presents a new method wherein the intensity and spatial distribution of the magnetic field, arising from large gradients in the local magnetization, are measured. These large gradients are induced by the passage of a test sample through the steep gradient field of a small, proximate permanent magnet. Magnetic field measurements during both directions of motion provide information indicative of the hysteresis properties. We theoretically analyze these measurements and show experimentally that the measurements correlate well with all the significant aspects of conventional hysteresis loops. The results given in this paper are qualitative, and the method is both by its simplicity and its sensitivity to important hysteresis features a powerful means of magnetic nondestructive evaluation.


---

## Body

## 1. Introduction

Hysteresis loss is not only a critical factor in the selection of steels for use in electrical machines but, by its nature, can also provide significant information on the structural condition and/or magnetic anisotropy of ferromagnetic materials generally. With electrical steels, the energy loss itself is the usual parameter of prime interest, whereas one or both of its (typically) key components, namely, coercivity and the remanent induction, or at some other “standard” field intensity, provide the sought for information. There is a need for a technique that can give useful information about relationships between microstructure or internal stress with magnetic performance.Conventional measurements of hysteretic properties usually employ ring or strip samples, with the latter being necessary when anisotropy is being explored [1]. Indeed, the directional nonuniformity of the magnetic properties is measured by cutting samples at various angles from a sheet of material. Nondestructive determinations of local properties on the other hand, for example, to explore the relative structural damage associated with different processes used to cut electrical steel sheets, typically require either specially prepared samples [2], sophisticated techniques such as the needle probe method [3], or specialized apparatus such as with the drag force method [4]. All methods except drag force employ magnetic fields derived from electric currents.A new magnetostatic method for obtaining comparative measurements of hysteresis loss and its components in ferromagnetic sheet materials is described in this paper. The measurements will be shown to provide qualitative information, which correlates directly with all of the significant aspects of conventional hysteresis loops. Essential features of the measurement apparatus are shown in Figure1. The distance between the Motion Limits is called the “Stroke.” The field sensor is located at the center of the Stroke, generally as close to the SUT surface as its physical package allows (c≈0). The size of the Stroke is typically from 10 G to 20 G. Remark that the presented method can measure the directional nonuniformity of the magnetic properties in a nondestructive way.(a) Side view of basic arrangement. (b) Magnified view of the active elements showing the dimensional factors that contribute to the “Gap”,G·x is the longitudinal distance between the magnet and the field sensor.
(a)(b)Themodus operandi of the measurement method is as follows: the magnet, polarized normal to the sample under test (SUT) surface, initially located at one end of the Stroke is moved, at a speed slow enough (~3 mm/s) such that making it slower does not significantly alter the results, across the face of the SUT at a small fixed separation distance (0.05–2.5 mm) until it reaches the opposite Motion Limit at which time the motion is stopped. It is then moved, at the same speed and separation distance, back to its starting position. This “initializing” cycle of motion is repeated several times to obtain stable magnetization excursions in the SUT. During the next forward traversal of the Stroke, both the intensity of the field (in the longitudinal direction), HF(x), and the coincidental longitudinal distance between the magnet and the field sensor, x, are measured and recorded. After a momentary stop at the end of the Stroke, motion is resumed in the reverse direction, during which, the field, HR(x), and corresponding position are again measured and recorded. For each recorded data pair (H(x),x), the difference, HF(x)-HR(x)=D(x), is calculated. As will be seen, D(x) will have some salient features (maxima, minima, and others) that correlate with the SUT’s hysteresis loss and its components.


## 2. Theory

### 2.1. Assumptions

The inarguable simplicity of the method and apparatus notwithstanding, a detailed understanding of its operation is quite the opposite. Nevertheless, the operational basis and aqualitative relationship between features of D(x) and hysteresis can be made adequately clear by an analytical model built from the following assumptions.(1)
Theexistence of the signature features in D(x) reflects magnetic hysteresis in the sample; in the absence of hysteresis D(x)≡0.(2)
The field from the magnet and its spatial distribution are approximated well enough by that of an equal moment dipole located within the magnet body.(3)
A solution of the 2D problem is sufficient.(4)
The sample is assumed to be thin enough such that the field from the magnet, while varying with longitudinal position, is uniform throughout its thickness. This assumption ignores thus also radial components of the field from the magnet.(5)
The instantaneous local magnetization at points within the sample,M(x), is determined entirely by the instantaneous local field from the magnet, HL(x), and the history of changes in that field. It is recognized that HF(x) and HR(x) are each comprised of components from two sources, the magnet and ∇·M(x) within the sample. Nevertheless, effects of these latter fields on M(x) are ignored on the assumption that they modify only details of the D(x) signature features, not their presence.(6)
TheM-H characteristics of the sample can be defined by any function which results in closed, symmetrical, sigmoidal loops without concern for underlying physical sources, for example, domain wall pinning or anisotropy. Thus suitable loops can be constructed from purposefully modified Langevin functions.The authors are fully aware that this model is an approximation, but it enables to directly correlate simple hysteresis properties with theD(x) measurements. This is the major aim of the theory shown here. More advanced numerical techniques could be used for understanding in depth the measurements.


### 2.2. Analysis

Guided by the previous assumptions, the analysis proceeds as follows.(a)
Determine the variation withx of the longitudinal component of the field (HL) from the equivalent dipole source.(b)
Determine the sequence of field variation at underlying points in the sample during forward and reverse motion of the magnet, that is,HL(x)F and HL(x)R.(c)
Create families of hysteresis loops with variable loss densities and components, that is,M(HL) functions with different values of coercivity and remanence.(d)
DetermineM(x) for both directions of motion by transposing the sequence of field variation, that is, HL(x)F and HL(x)R onto the M(HL) loops created in (c).(e)
Determine∇·M(x) (= φ, the free-pole density) for each direction of motion. (This becomes a second source (in addition to m) of the sensed field.)(f)
Determine∇φ(x)=-H for each direction of motion, that is, HF(x) and HR(x).(g)
DetermineD(x)=HF(x)-HR(x), identify, and characterize salient features.(h)
Correlate feature characterizations from (g) with hysteresis loop features from (c).For the arrangement diagrammed in Figure2, the relationship between HL at point P (the location of the field sensor in Figure 1), and P’s distance from a dipole of moment m, is found (following from Cullity’s derivation [5]) as(1)HL=3mxG(x2+G2)5/2.Figure 2
Diagram showing the fieldH, present at a point P, distance r, and angle q (x and G) from dipole m.The longitudinal componentHL is directed to the right at points to the right of O (x>0) and to the left at points to the left of O (x<0), thus HL(x) = -HL(-x). Equation (1) with m=1 and G=1 is plotted in Figure 3. The numbers indicate noteworthy magnet (dipole) positions relative to a field sensor at x=0 during its motion in the forward, F, and reverse, R, directions of motion. When, for example, the magnet is at position 1, HL1 = −0.00435, at position 2, HL2 = −0.8587, at position 3, HL3 = 0, and so forth. Since the field distribution is effectively “attached” to the magnet, it moves together with the magnet. Determination of M(x) requires knowledge of the corresponding M-H relationships for the sample material. Assumption 6 permits suitable, analytically defined M-H functions to be created from modified Langevin functions.Figure 3
Variation in (normalized)HL with distance x. Note that HL is symmetrical around the origin, has peak values of ±0.8587 at x=±0.5, and becomes only ±0.00435 at x=±5.Ascending and descending limbs of hypothetical hysteresis loops are, respectively, generated from(2)Ma=coth⁡(k(H-Hc))-1k(H-Hc),(3)Md=coth⁡(k(H-Hc))-1k(H+Hc),
wherein the “constants” Hc and k provide means to vary the “coercive field” and “remanence ratio,” respectively, thereby to simulate different sample materials. Indeed, we use Hc to shift the function along the H axis, to the left (+Hc) to obtain the descending limb of the loop and to the right (-Hc) to obtain the ascending limb. Similarly, we use k to adjust the tilt of these loop limbs, which modifies the remanence ratio.A closed loop is formed by shiftingMa upward and Md downward by 1/2 the difference, Δ=Md-Ma, when H is at its peak value, HP, in (2) and (3). The loop will then consist of an ascending limb: (4)MA=Ma+Δ2,
and a descending limb:(5)MD=Md-Δ2.Hypothetical positive and negativeminor limbs from remanence (H=0) to H=±HP can be produced from(6)Mm+=MA+MD-MA1+f*H,Mm-=MD+MA-MD1-f*H,
where f(≪1) is selected to insure that  |Mm| increases continuously with increasing |H| while staying always between MA and MD. Of the three shape modifiers (Hc, k, and f), only Hc has physical dimensions (the same as those of H).From the above equations, it is possible to determine∇·M(x)=φ for both the forward and reverse directions: φF=(∂Mm+(x))/∂x and φR=(∂Mm-(x))/∂x, from which HF(x)=∇φF(x) and HR(x), and thusD(x), can be calculated.To avoid the need to deal with equations having ever growing numbers of ever more complex terms, and to provide means for graphically following the evolving analysis, we assign arbitrary (but as will be seen, arguably reasonable) values to the material-dependent parametersHc, k, and f, and set HP = 15. In this way, MA, MD, Mm+ (for 0≤H≤HP), and M-m (for 0≥H≥-HP) become numerical functions of H. For similar reasons and to make HLPeak = HP = 15, we assign the value 5/0.8587 = 17.47 to m in (1).As the magnet moves forward fromx=-5 to x=5, HL(0) and M(0) will follow the sequence shown in Figure 4(a), and during reverse motion (from x=5 to x=-5), the sequence shown in Figure 4(b). The corresponding variations of M(0)F and M(0)R with magnet position are shown in Figure 5(a), with the central regions of ∇(∇·M(0)F) and ∇(∇·M(0)R) shown in Figure 5(b). The starting points in Figure 5(a) are not exactly the same as the numbered points in Figure 4. This is because the field at the end of the motion in either direction does not reach zero at any finite values of x. Figure 4(a) presumes that the value of x at the motion extrema is large enough for the field starting and finishing points to be close enough to zero. In contrast, the start (and finish) of the plots (but not the motion) in Figure 5 clearly occur at finite (±8.8 G) distances from the sensor. Thus the end of the plot in one direction does not have the same magnetization as the start of the plot in the opposite direction.(a) Variation of field and resulting magnetization at the point on the SUT, which is directly over the field sensor during forward motion of the magnet. (b) Same for reverse motion. (Hc=2, k = 1, f=0.135).
(a)(b)(a) Variation ofM at the sensor position during forward (F) and reverse (R) motion. (b) Corresponding sensed field = ∇(∇·M(0)F) and ∇(∇·M(0)R).
(a)(b)Major hysteresis loops created from (4) and (5) for the values of Hc and k indicated are shown in Figures 6(a) and 7(a). Resultant values of D(x) are shown in Figures 6(b) and 7(b). The sizes of the central positive peak and the symmetrical negative peaks are seen to vary directly with the loop squareness, peak M, and Hc, as does hysteresis loss.(a) Hysteresis Loops for constantHc, variable k values. (b) Resulting D(x).
(a)(b)(a) Hysteresis Loops for constantk, variable Hc values. (b) ResultingD(x).
(a)(b)


## 2.1. Assumptions

The inarguable simplicity of the method and apparatus notwithstanding, a detailed understanding of its operation is quite the opposite. Nevertheless, the operational basis and aqualitative relationship between features of D(x) and hysteresis can be made adequately clear by an analytical model built from the following assumptions.(1)
Theexistence of the signature features in D(x) reflects magnetic hysteresis in the sample; in the absence of hysteresis D(x)≡0.(2)
The field from the magnet and its spatial distribution are approximated well enough by that of an equal moment dipole located within the magnet body.(3)
A solution of the 2D problem is sufficient.(4)
The sample is assumed to be thin enough such that the field from the magnet, while varying with longitudinal position, is uniform throughout its thickness. This assumption ignores thus also radial components of the field from the magnet.(5)
The instantaneous local magnetization at points within the sample,M(x), is determined entirely by the instantaneous local field from the magnet, HL(x), and the history of changes in that field. It is recognized that HF(x) and HR(x) are each comprised of components from two sources, the magnet and ∇·M(x) within the sample. Nevertheless, effects of these latter fields on M(x) are ignored on the assumption that they modify only details of the D(x) signature features, not their presence.(6)
TheM-H characteristics of the sample can be defined by any function which results in closed, symmetrical, sigmoidal loops without concern for underlying physical sources, for example, domain wall pinning or anisotropy. Thus suitable loops can be constructed from purposefully modified Langevin functions.The authors are fully aware that this model is an approximation, but it enables to directly correlate simple hysteresis properties with theD(x) measurements. This is the major aim of the theory shown here. More advanced numerical techniques could be used for understanding in depth the measurements.


## 2.2. Analysis

Guided by the previous assumptions, the analysis proceeds as follows.(a)
Determine the variation withx of the longitudinal component of the field (HL) from the equivalent dipole source.(b)
Determine the sequence of field variation at underlying points in the sample during forward and reverse motion of the magnet, that is,HL(x)F and HL(x)R.(c)
Create families of hysteresis loops with variable loss densities and components, that is,M(HL) functions with different values of coercivity and remanence.(d)
DetermineM(x) for both directions of motion by transposing the sequence of field variation, that is, HL(x)F and HL(x)R onto the M(HL) loops created in (c).(e)
Determine∇·M(x) (= φ, the free-pole density) for each direction of motion. (This becomes a second source (in addition to m) of the sensed field.)(f)
Determine∇φ(x)=-H for each direction of motion, that is, HF(x) and HR(x).(g)
DetermineD(x)=HF(x)-HR(x), identify, and characterize salient features.(h)
Correlate feature characterizations from (g) with hysteresis loop features from (c).For the arrangement diagrammed in Figure2, the relationship between HL at point P (the location of the field sensor in Figure 1), and P’s distance from a dipole of moment m, is found (following from Cullity’s derivation [5]) as(1)HL=3mxG(x2+G2)5/2.Figure 2
Diagram showing the fieldH, present at a point P, distance r, and angle q (x and G) from dipole m.The longitudinal componentHL is directed to the right at points to the right of O (x>0) and to the left at points to the left of O (x<0), thus HL(x) = -HL(-x). Equation (1) with m=1 and G=1 is plotted in Figure 3. The numbers indicate noteworthy magnet (dipole) positions relative to a field sensor at x=0 during its motion in the forward, F, and reverse, R, directions of motion. When, for example, the magnet is at position 1, HL1 = −0.00435, at position 2, HL2 = −0.8587, at position 3, HL3 = 0, and so forth. Since the field distribution is effectively “attached” to the magnet, it moves together with the magnet. Determination of M(x) requires knowledge of the corresponding M-H relationships for the sample material. Assumption 6 permits suitable, analytically defined M-H functions to be created from modified Langevin functions.Figure 3
Variation in (normalized)HL with distance x. Note that HL is symmetrical around the origin, has peak values of ±0.8587 at x=±0.5, and becomes only ±0.00435 at x=±5.Ascending and descending limbs of hypothetical hysteresis loops are, respectively, generated from(2)Ma=coth⁡(k(H-Hc))-1k(H-Hc),(3)Md=coth⁡(k(H-Hc))-1k(H+Hc),
wherein the “constants” Hc and k provide means to vary the “coercive field” and “remanence ratio,” respectively, thereby to simulate different sample materials. Indeed, we use Hc to shift the function along the H axis, to the left (+Hc) to obtain the descending limb of the loop and to the right (-Hc) to obtain the ascending limb. Similarly, we use k to adjust the tilt of these loop limbs, which modifies the remanence ratio.A closed loop is formed by shiftingMa upward and Md downward by 1/2 the difference, Δ=Md-Ma, when H is at its peak value, HP, in (2) and (3). The loop will then consist of an ascending limb: (4)MA=Ma+Δ2,
and a descending limb:(5)MD=Md-Δ2.Hypothetical positive and negativeminor limbs from remanence (H=0) to H=±HP can be produced from(6)Mm+=MA+MD-MA1+f*H,Mm-=MD+MA-MD1-f*H,
where f(≪1) is selected to insure that  |Mm| increases continuously with increasing |H| while staying always between MA and MD. Of the three shape modifiers (Hc, k, and f), only Hc has physical dimensions (the same as those of H).From the above equations, it is possible to determine∇·M(x)=φ for both the forward and reverse directions: φF=(∂Mm+(x))/∂x and φR=(∂Mm-(x))/∂x, from which HF(x)=∇φF(x) and HR(x), and thusD(x), can be calculated.To avoid the need to deal with equations having ever growing numbers of ever more complex terms, and to provide means for graphically following the evolving analysis, we assign arbitrary (but as will be seen, arguably reasonable) values to the material-dependent parametersHc, k, and f, and set HP = 15. In this way, MA, MD, Mm+ (for 0≤H≤HP), and M-m (for 0≥H≥-HP) become numerical functions of H. For similar reasons and to make HLPeak = HP = 15, we assign the value 5/0.8587 = 17.47 to m in (1).As the magnet moves forward fromx=-5 to x=5, HL(0) and M(0) will follow the sequence shown in Figure 4(a), and during reverse motion (from x=5 to x=-5), the sequence shown in Figure 4(b). The corresponding variations of M(0)F and M(0)R with magnet position are shown in Figure 5(a), with the central regions of ∇(∇·M(0)F) and ∇(∇·M(0)R) shown in Figure 5(b). The starting points in Figure 5(a) are not exactly the same as the numbered points in Figure 4. This is because the field at the end of the motion in either direction does not reach zero at any finite values of x. Figure 4(a) presumes that the value of x at the motion extrema is large enough for the field starting and finishing points to be close enough to zero. In contrast, the start (and finish) of the plots (but not the motion) in Figure 5 clearly occur at finite (±8.8 G) distances from the sensor. Thus the end of the plot in one direction does not have the same magnetization as the start of the plot in the opposite direction.(a) Variation of field and resulting magnetization at the point on the SUT, which is directly over the field sensor during forward motion of the magnet. (b) Same for reverse motion. (Hc=2, k = 1, f=0.135).
(a)(b)(a) Variation ofM at the sensor position during forward (F) and reverse (R) motion. (b) Corresponding sensed field = ∇(∇·M(0)F) and ∇(∇·M(0)R).
(a)(b)Major hysteresis loops created from (4) and (5) for the values of Hc and k indicated are shown in Figures 6(a) and 7(a). Resultant values of D(x) are shown in Figures 6(b) and 7(b). The sizes of the central positive peak and the symmetrical negative peaks are seen to vary directly with the loop squareness, peak M, and Hc, as does hysteresis loss.(a) Hysteresis Loops for constantHc, variable k values. (b) Resulting D(x).
(a)(b)(a) Hysteresis Loops for constantk, variable Hc values. (b) ResultingD(x).
(a)(b)


## 3. Experimental

An experimental apparatus was set up on a vertical milling machine, thereby conveniently accommodating a variety of experimental conditions, including SUT size (width, thickness, and length), field sensor position relative to the SUT edge, length of Stroke and center position relative to the field sensor, adjustability of both the space between the magnet face and SUT (dimensionb in Figure 1), and the location of the magnet relative to the field sensor centerline. An Allegro 3515 UA Hall effect field sensor (sensitivity = 0.0628 mV/A/m), with d (Figure 1) = 1.62 mm, was mounted just below (~0.05 mm) the mounting surface of the SUT. Nd-Fe-B magnets of a variety of sizes and energy levels were mounted into aluminum holders having common mounting features (a 9.5 mm diameter × 20 mm long cylindrical portion) to allow easy interchangeability in the machine spindle.Except that the magnet was stationary while the SUT/field sensor combination (being mounted on the milling machine table) were the movable elements, operation of the apparatus followed the description in the Introduction.For each magnet and gap combination,HL(x) was first measured without any SUT. The comparative amplitudes, ±HLPeak for the various magnet/“gap” (G in Figure 1) combinations were used to gain understanding on how these parameters affect the later test results. On the basis of the plot of (1) in Figure 3, together with assumption 2, G was presumed to be the distance between ±HLPeak. Measurements of this distance with two or more different physical gaps allowed both the location of the equivalent dipole (“a” in Figure 1(b)) and the dipole moment to be calculated.Values ofD(x), calculated from measurements of HF(x) and HR(x) for the samples and magnet indicated, are plotted in Figures 8 and 9. The most salient features of the plots are seen to match the prediction of the analysis. Moreover, since losses in electrical steels are in approximate proportion to their grade number, and the drastic reduction in coercivity by annealing of cold worked steels is common knowledge, these plots verify the theoretically expected correlation between the amplitude of the signature features and the factors contributing to hysteresis loss. The nominal coercivity values of the three samples shown in Figure 8 are 100 A/m, 45 A/m, 35 A/m for the 800-50, 350-50, and the 290-50, respectively. Comparison with Figure 8 teaches us that there is indeed a larger difference between the 800-50 and the 350-50 than between the 290-50 and the 350-50.Figure 8
D(x) from measured values of HF(x) and HR(x) for the 3 grades of NO Si steels indicated.Figure 9
Same as Figure8 for 0.61 mm thick cold rolled AISI 1010 in as-received and annealed conditions.The versatility of the described method was shown by its use to measure the relative losses in regions near the cut edges of strip samples of the 350 and 800 grades. The SUT was placed with the edge being examined ~1 mm over the center of the field sensor. A magnet, 3.18 mm square by 12.7 mm long (in the direction ofm) was positioned 1 mm back from the edge and 2.5 mm above the SUT. D(x) measurements on 2 samples each, cut by a fast moving laser, averaged 40% (350 grade) and 14% (800 grade) less than those cut more slowly, in good agreement with results from other measurement methods [4].The more compact spacing of the signature features predicted in the analysis than those experimentally observed is attributed to several critical but enormously complicating factors being ignored in the analysis; namely, the interaction between the magnetization and the generated fields and the gradients in these throughout the SUT thickness. Nevertheless, the model correctly shows that the magnetization gradients are different for the two directions of motion, a difference founded on the double-valuedM-H relationships which manifest “hysteresis.” We also ignored in this paper the effects of sample thickness, distance of the field sensor from the sample surface, possible effects of the circuit permeance on the field from the magnet, and so forth. However, a significant contributor to the less compact detected field difference profile is this interaction of the ∇·M fields with the M fields.This paper proposes an experimental method for locally assessing hysteresis properties, and a theoretical basis is given for explaining the obtained difference curves. The correlation between the experimentally obtained Figure8 and Figures 7(a) and 7(b) shows that indeed the difference curves are correlated with the hysteresis losses. The same correlation was observed with Figure 9. The results in this paper are qualitative, and by refinement of the numerical model for explaining the difference curves we can aim in further work to more quantitative assessment of hysteresis properties and losses. The numerical model can be refined by, for example, taking into account the thickness of the sample under study so to have a more correct explanation of the measurements. The method we describe does give quantitative results. We do not claim, however, that the quantitative description of salient features of D(x) plots have translatable correlations with quantitative descriptions of salient features of B-H plots.A possibility to more quantitative studies is to determine the input model parameters (i.e.,k, Hc, and f) starting from directly measured B-H loops. Also, more advancements are needed in the theory for better understanding the D(x) measurements. For instance, sensitivity analysis can help us in better understanding the measurements by looking at how the variation in one parameter affects the D(x).


## 4. Conclusions

In this paper, we propose a unique methodology together with its modus operandi and important experimental results. Both analytical and experimental results convincingly show the signature features of the difference measurementsD(x) to be reflective of those same properties of test sample material, which are the underlying determinants of the size, shape, and intercept features of conventional hysteresis loops. As such, and considering both the rapidity and sensitivity by which D(x) measurements might be obtained with purposefully constructed apparatus of obvious simplicity, the utility of the described method for comparative evaluations by nondestructive means, of factors affecting hysteresis properties, seems well established. The proposed technique has the possibility to assess in detail the relationships between microstructure (or internal stress) and magnetic performance.



---

*Source: 870916-2012-05-14.xml*