Internal Demagnetizing Factor in FerrousMetals

Modelling the saturation major loop of a ferrous metal produces the intrinsic magnetization parameters; fitting the measured commutation curve, however, can yield different results. The relation of the intrinsic loci of the vertices of the minor loops (ND = 0) to the experimental curve (ND / = 0) is investigated. The two-way transformation between the two curves is formulated in closed mathematical form with the help of the internal demagnetization factor, ND . The method is applied to four ferrous metals, with widely different intrinsic properties (soft nonoriented Fe-Si steel, normalized low carbon steel, and Finemet in nanocrystalline and amorphous state) supporting the predictions of the proposal. The developed relationship is model independent and it is shown that the ND factor depends linearly on coercivity based on experimental evidence.


Introduction
A large number of the measurements of ferrous substances are aimed at finding the intrinsic material properties [1] of the tested ferrous sample (as defined by Fiorillo).Due to the ever-presence of demagnetization field, various measuring methods have been developed to minimize its effect.The most commonly accepted way is to make the sample turn into a closed magnetic circuit, such as a toroid or an Epstein square [1][2][3].Although these two methods are not completely free from the ever-present internal demagnetization, they suffer the least from it [1].Researchers went into great length to include the internal demagnetization force into current models like Preisach, Jiles, Stoner-Wohlfarth [2,4,5], and so forth, leading to complicated, so called, dynamic versions.
The saturated major hysteresis loop of the sample carries all the intrinsic magnetic parameters directly recoverable from the measured data.Within this loop lie the unhysteretic loci of the vertices of the symmetrical minor loops, the only curve, which belongs to both the ascending and descending branches of the hysteresis loops [5,6].
A proposal is put forward in this paper to show the relationship between the intrinsic curve and the loci of vertices of the measured minor loops.This relationship between the two curves, independent of models, is formulated in closed mathematical form and its prediction is verified by the experimental data obtained from four different ferrous samples.
Once the intrinsic locus (M 01 for N d = 0) is modelled from saturation or minor loop data, by using any of the static models, the measured curve (N d / = 0) can be calculated from the proposed formulation below with optimization of the N d parameter value as specified by Jiles [2] and used by Fiorillo [1].

The Intrinsic Loci of Vertices
The intrinsic commutation curve (N D = 0) is the locus of the return points or the maxima of the set of symmetrical minor loops.It is a single-valued function and in spite of having no hysteresis, it carries all the hysteretic properties of the ferrous material obtained from the saturation data [7].
M 0 is the arithmetic mean of the ascending and descending magnetization functions, M u and M d , with equal application to N D = 0 case as well. (1)

The Effective Field and Its Implications
For the description of the effect of the internal demagnetizing field we will use the concept of the effective field, which is analogous to the Weiss mean field as defined and used by other authors [2,10,11].Although in simple cases, proportionality is assumed between H D interaction field and M magnetization vector, this can only be regarded as a linear approximation [12,13] expressed as H D = N D M in scalar form, leading to where N D is the internal demagnetization factor.
For simplicity we will use normalised quantities in further calculations, where the lower case letters will represent the normalised quantities of the physical equivalents, denoted by the same capital letters.m 0 denote the normalised M 0 function in (1).With the normalized effective field, h eff , m 0 can be expressed as: where n D is also normalized.(1 Tesla = 8.10 5 A/m.The internal demagnetization factor has unity dimension only when the magnetization is measured in A/m.).The first derivative of m 0 by h in (3) leads to an expression, which shows a character similar to the feedback in an electrical circuit [14].
This expression describes a well-known relationship between the inherent (μ i ) and the effective (μ eff ) permeabilities [1,2].
For most magnetic substances the value of N D is small in the order of ∼ −10 −5 with unity dimension [1].Consider By using the hysteretic model [6] (see also the appendix), the integral of ( 4) by h leads us to the following expression for m 0 , when . ( After integration: When the integration constants are appropriately chosen, this form of m 0 is equivalent to the one given in (A.5).The intrinsic locus is entirely a theoretical concept.It was introduced for the free (Gibbs) energy calculations [2].It assumes zero internal demagnetization in a system, where a m moment can freely move around under the influence of H external excitation field without any hindrance from the interaction between the magnetic moments.
We must remind the reader that, for various ferromagnetic substances, this internal demagnetization constant N D is given traditionally in a numerical value with unity dimension (i.e., when both the H and B measured in A/m).When different unitary system is used, N D has a different physical dimension and must be normalised (see n D as normalised N D ).

Experimental Verification
To verify the predictions of the proposed method, it was applied initially to two ferrous materials with very different characters.The first was a soft steel NO Fe-Si with 67.5 A/m coercivity [15], shown in Figure 1.The second material was a normalized low carbon steel (AISI 1040) with coercivity of 450 A/m (see Figure 5).Following the excellent results, the same experiment was also repeated later on two other samples; Finemet in nanocrystalline and amorphous (as cast) state.The detailed data of those samples are not included in this paper due to its limited size.All measurements were carried out under identical conditions by using triangular excitation of f = 0.001 Hz on toroid samples with geometrical details as follows: D e = 25 mm, D int = 15 mm and thickness d = 0.5 mm [1].In making the identical toroid samples, extreme care was taken to avoid any changes in magnetic properties due to mechanical handling of the materials.
For numerical calculations we used the hyperbolic model for its simplicity and speed and the Mathematica program interactively.A brief summary of the model is given in the appendix.For further details we refer the reader to the literature [6,16,17].When the iteration produced the best fit to the measured curve, the normalized and the equivalent physical values can be easily read from the two coordinate systems (normalization) as shown in Figures 1, 2, and 3.
The first sample (see Figure 1) was modelled with the following normalized and physical parameters: In order to check the accuracy of the transformation, between the two curves, 25 of the minor loops were measured with maximum field excitation values between h m = 5.01, Hm = 501 A/m and h m = 0.152, Hm = 15.2A/m.For all the minor loops measured, the corresponding loops were calculated for the reduced maximum magnetization at n D = 0 and n D = −0.151.
For clarity, only one of the minor loops is shown in Figure 3 for the peek excitation field value of 111.3 A/m.
All the calculated loops had an excellent fit to the equivalent measured loops for the same N D value.The proposed method is applicable to all measured hysteretic data, where the external demagnetization field is eliminated or reduced to a negligible level.

Low Carbon Steel Toroid (AISI 1040).
For the second sample, we selected a toroid from low carbon steel with coercivity of 450 A/m.The use of the same experimental setup has yielded the following parameters in normalized and physical units: The measured and the modelled hysteresis loops are depicted in Figure 4 with the intrinsic and the measured

N d as a Function of Coercivity
The experimental results have indicated that N D is greatly dependent on the coercivity of the sample.To verify this dependency, four samples, listed under Section 4, were tested for this purpose, with coercivity H c ranging between 1.6 A/m and 450 A/m. Figure 5 depicts the relationship between N D and H c .
The graph, in Figure 5, was constructed by using the parameters of the samples listed below, whose coercivity values change in steps approximately by an order of magnitude: The graph shows N D linear dependency on H c coercivity and based on experimental evidence, it can be approximated as where B m represents the saturation magnetization.

Conclusions
The relation between the intrinsic (N D = 0) and the measured (N D / = 0) loci of the vertices of the minor loops was investigated.By using the effective field, this relationship was formulated.The method was subjected to tests on four magnetic materials with widely ranging magnetic properties.The test results show that the mathematical approach, presented here, describes the relation well.The N D factor linear dependence on coercivity, based on empirical evidence, was also demonstrated.The paper shows that magnetic parameters (H c and B m ) N D can be estimated very close to the real value.
The proposal's aim is to recover the intrinsic magnetization properties from the measured commutation curves.constant, while h m is the maximum field excitation, common to all.The index k refers to the individual processes and n is the number of total processes involved (for most substances n = 3).
The parameters are calculated by changing the model parameters until the best fit to the measured curve is achieved.When the iteration gives the best fit, the normalized and the equivalent physical values can be read from the two coordinate systems (normalized and measured), as shown in Figures 1, 2, and 3.

Figure 1 :
Figure 1: The measured and the modelled major loop of NO Fe-Si soft steel.

4. 1 .
NO Fe-Si.Starting with conveniently chosen parameter values at the beginning with subsequent changes of the parameters, new curves are calculated and compared with the measured one.
n D = 0 Dotted line calculated n D = −0.151Broken line measured

, m d , m 0 Figure 4 :
Figure 4: Measured (broken line) and modelled (solid line) hysteresis and the commutation curves for n D = 0 and n D = −0.08 (dotted line) in the first quadrant for the low carbon steel (AISI 1040).

Figure 5 :
Figure 5: N D as the function of coercivity H c .