The Static Characteristic Loop and the External Demagnetizing Factor ()
1. Introduction
A practical way of characterising a magnetic material is to measure its static hysteresis loop in a closed magnetic circuit. This usually takes form in a toroid or an Epstein square [1], which minimises most of the stray and side effects of other measuring methods. This however is difficult in practice, considering the shapes and sizes of the samples of the magnetic materials in industrial applications. Since the shape of the sample changes its measured magnetic properties, it is desirable to use samples from the final product, which normally come in the shape of a rod or a sheet, for industrial purposes. This however creates some difficulties, with making them into a toroid or Epstein square sample. Most test results therefore are obtained by using VSM (Vibrating-sample Magnetometer) [1-3] method. This method, with excellent repeatability, is now a widely accepted for industrial tests. The data however suffers from the effect of demagnetization, due to the open circuit geometry of the measuring arrangement. The need, to obtain the real, static parameters from measured data of the tested materials, has motivated people to work out various methods to reduce the effects, influencing the measured data [4-7].
At the beginning we have assumed that, only one of the calculated parameters is affected by the external demagnetization. Test results showed, that this assumption was right. It is possible to work out the un-sheared parameters (equivalent to closed circuit toroidal measurements) from the open circuit results of a VSM results. It will be shown, that this procedure requires only change in the effective slope of the loop with no change in the other major parameters of the modelled sample [8].
It was initially assumed that, the error coming from the ever present internal demagnetisation and from other sources is smaller than the [9] (usually at least an order of magnitude smaller) effect of shearing. This is however only conditionally true and all major models (like the Moving Preisach and Dinamic Jiles-Atherton etc.) [10,11] are modified to include, the ever present, internal demagnetization. We also assumed that the effect of shearing depends only on the magnitude of the demagnetization factor and the mechanical treatment not introducing changes in the magnetic properties of the materials. Also assumed, that it is independent of the method of measurement.
Cullity, who studied the effect of shearing on the internal demagnetization in detail, published the relation specifically to the shape and size of the specimen in graphical form [11-15]. In spite of all the assumptions made, the agreement between the theoretical and experimental results is remarkable. Further experiments are in progress for further verification this simplified method. Work is also extended to include the so far neglected effects to improve on the accuracy.
2. Measuring Arrangements
Preceding the test, the toroidal sample was carefully demagnetized by applying an alternating field of f = 10 Hz with logarithmically decreasing amplitude in 5000 steps from saturation to zero [1,2]. The hysteresis loop was measured with a triangular excitation at a frequency of f = 0.001 Hz and integrated with a Walker integrator (see Figure 1). The rate of change of the field here is dependent on both the frequency f and field amplitude Hm therefore at Hm = 20 A/m, the rate of change, comes to dH/dt = 4f·Hm = 0.08 A/m·s.
The static hystersis loop can only be obtained at low rate of change, particularly when the extremely soft sample has a characteristic square-like hysteresis loop [1,2].
To avoid phase shift between H and B, the current measuring resistance R (see Figure 1) has always satisfied the 2πfL << R criterion.
Although the same test was repeated on ultrasoft, Finemet, from nanocrystalline to Mn-Zn ferrite, NO Fe-Si and low and high carbon steel cores, due to the limited length of the paper, here only the data of one illustrative experiment are included.
3. Description of the Model
The analytical approach is based on the assumption, that in general, particularly in case of soft irons, there are three parallel processes, dominating the overall magnetization process i.e. the reversible and irreversible domain wall movement (DWM), the reversible and irreversible
Figure 1. Schematic block diagram of an experimental arrangement for measuring the toroid sample.
domain rotation (DR) and the domain wall annihilation and nucleation (DWAN) processes. Although these processes are interlinked, they can be mathematically formulated separately and combined, by using Maxwell’s superposition principle. They individually dominate the low, middle and near saturated region of magnetization and all are supposed to have sigmoid shapes. This model is already used in number of applications [15-18] and it is well documented in the literature therefore, here only a brief summary of the relevant formulation will be given in, canonic form [19,20].
The contribution of the individual processes to the combined hysteresis loop can be described by the following normalized mathematical equations.
(1a)
(1b)
(2a)
(2b)
for h = hm (3)
Here fuk and fdk are the normalized ascending and descending magnetization functions respectively, h is the field excitation and hck is the coercivity of the kth process. ak is the amplitude of the component processes present, αk is the slope, nS is the shearing factor (nS is unity for un-sheared loops, see later) and f0k is the integration constant [7,21], while hm represents the maximum field excitation. The index k refers to the individual component processes and n is the total number of processes involved. For most of the magnetic materials used in practice, n equals 3 (see beginning of this section).
4. Experimental Results
We proposed at the beginning, that the changes in the shape of the hystetesis loop measured on a closed toroid can be described by an appropriate change in the αknS product (see Equation (2)). Here αk represents the slope of the major loop, measured on a toroid sample (for nS = 1, nD = 0) in the model. For demonstration, we used the measurements on a toroidal sample made of NO Fe-Si soft magnetic material. The geometrical details of the toroids used, in the experiment, were: Dex = 25 mm, Din = 15 mm and thickness d = 0.5 mm [1].
In the first part of the experiment, the toroid was magnetised to Bm = 1.68 T by applying a linearly variable field of Hm = 2750 A/m maximum amplitude for measuring the hysteresis curve. Then, two 0.5 mm air gaps were cut in the magnetic circuit at a slow rate and the measurement was repeated. Great care was taken not to introduce additional side effects at cutting. This was for the simulation of an open-loop VSM measurement with the same ND and NS. The two hystetesis loops (before and after cutting) are shown in Figure 2. The measured loop, without the gap, was modelled by using Equations (1)-(3). When the iteration (manual or computerised) produced the best fit to the measured curve, the normalized and the equivalent physical values can be easily be read from the two coordinate systems (normalized and measured), as shown in Figure 3.
The iteration yielded the following normalised (lower case letters) and corresponding physical parameters (capital letters) for the best fit:
a1 = 1.17, a2 = 0.14, a3 = 0.35α1 = 4.5, α2 = 1.2, α3 = 0.143hc1 = 0.115, hc2 = 0.45, hc3 = 0.575hm = 7.5.
Corresponding to:
A1 =1.17 T, A2 = 0.14 T, A3 = 0.35 THc1 = 44 A/m, Hc2 = 164.9 A/m, Hc3 = 210.75 A/m.
Hm = 2750 A/m.
The normalization: 1b = 1T, 1h = 382.6 A/m.
By altering the shearing coefficient ns, in the iteration and using the new values of α1·nS = 0.28, α2·nS = 0.073, and α3·nS = 0.0087, we obtained the closest fit to the sheared hystetesis loop. The best fit appeared to be, when nS was around 0.061.
From the geometry of the magnetic circuit, by using the magnetic Ohm’s law, ND (see Equation 7 for the relation between ND and NS) the demagnetization factor, can be calculated by using Equations (4) and (5). We assumed here, that the amplitude of the magnetization vector is the same in the iron as in the air gap. So:
Figure 2. The major hysteresis loops of a NO Fe-Si sample measured on a toroid with and without 1 mm gap.
Figure 3. The measured and calculated hysteresis loops of NO Fe-Si material with and without external demagnetization.
(4)
when μir relative permeability is used and it is assumed that la li then
(5)
where
(6)
is the effective magnetic field excitation.
(7)
or in normalised form
(8)
Here li and la are the mean paths of iron and air respectively, μi and μa represent the permeability of iron and air (free space) and A is the cross section area of the iron core. It shows that, when NS the shearing coefficient is unity (no shearing), then the ND demagnetization factor is zero. The expression in (9) represents the effective permeability.
(9)
We must remind the reader of the following: ND is given traditionally in a value with unity dimension, (i.e. when both the H and B measured in A/m). When different unitary system is used then ND has a different physical dimension and must be normalised (see nD as normalised ND).
By using the geometrical details of the toroid, given above with 1 mm effective gap, the physical shearing coefficient numerical value will come to
(10)
here the initial permeability, calculated from the measured static hysteresis loop data of the soft iron sample, was around μi = 1000.
Figure 4 depicts the calculated major loops.
In the second experiment the closed toroidal sample was magnetized into deeper saturation by applying a high Hm = 4125 A/m maximum excitation. Following that, two 1.5 mm cut (effective 3 mm gap) were made in the toroid and the measurement was repeated.
This time however the maximum excitation was increased to approach Bm = 1.68 T, achieved without the gap.
The loops were modelled by using the same numerical parameters at a larger hm maximum field excitation. The resulted loops (measured and calculated) are depicted in Figure 3. The agreement between the experimental result and the modelled one is extremely close.
Other additional factors like the internal, shape demagnetization factors and/or stresses were taken as negligibly small.
The separation and calculation of αk and nS from the experimentally obtained αk·nS product, requires the presence of two invariants in the transformation.
It is imperative to look at the two loops (un-sheared and sheared) at the same maximum magnetization (or as near as possible). Under this condition we can claim to have the necessary two invariants i.e. the coercivity and the maximum magnetisation.
It can be shown, that the two smaller processes (DR, DWAM) have often a negligible effect on the overall coercivity and experience shows that in most applications,
Figure 4. The modelled major loops before and after shearing.
particularly for soft irons, it is accurate enough in most cases to use the equations describing the dominant process. The proposal is valid, however, when the DR and DWAM processes are not negligible. All resultant parameters, required for the transformation can be calculated (see Appendix) from the parameters of the processes involved. It is enough therefore to demonstrate the process by using one major component [22]. On this ground the present example is given entirely for illustrative purposes.
By using (2a), we can write for the sheared and the un-sheared dominant loop:
(11)
(nS being unity in the first case)
From this, the shearing coefficient nS may be calculated as:
(12)
where hm1 and hm1s are the amplitudes of field excitations necessary to achieve the same maximum magnetization in the un-sheared and the sheared sample respectively.
In soft steel, when the coercivity hc is small relative to the peak excitation, it is often good enough to use the ratio of the maximum field, from the saturation loop, measured by VSM. In the knowledge of nS the value of ακ could be calculated from the measured product.
The external shearing coefficients, calculated from the measured data was ns = 0.065 for the first experiment (1 mm gap) and ns = 0.02 for the second one (3 mm gap). Some allowances should be made to the inherent internal demagnetization present in all type of measurement [2]. The estimated value of the internal demagnetization field was max 0.2 T [9].
From the second invariant, the ratio between the arctanh values of the remanences (sheared and un-sheared) (see Equation 2(b)) is also equal to nS.
The measured remanence of the 1 mm gap sample was 0.032 T at Bm = 1.41 T, while the corresponding modelled figure was 0.035 T. In the case of the 3 mm gap sample, the measured and the modelled remanence figures were 0.0085 T and 0.0066 T respectively, taken at Bm = 1.5 T. Based on these figures, Br the projected remanence of the loop, transformed back onto the unsheared (nS = 1 and nD = 0 ) plane came to 0.636 T.
It can be shown also, that the transformation does not change the area enclosed by the hysteresis loop, therefore sheared or un-seared, it has the same hyteretic losses.
5. Conclusion
It was proven, that shearing and/or un-shearing will not affect most of the static characteristic, material parameters of the soft sample. It is also shown, that the external demagnetization affects one parameter only, namely the effective slope (i.e. nSα) of the hysteresis loop. From the sample’s geometry the shearing coefficient can often be calculated. With the proposed method and the knowledge of ND [2] one can exploit the advantage of the open geometry measuring methods for unusual samples (i.e. balls, sheets or rods with known ND) and calculate the equivalent closed circuit static parameters.
The transformation is based on the two invariants, the coercivity and the maximum excitation.
This model independent transformation can be an useful tool in practice, in experimental work as well as in theoretical applications. The very close agreement between the theoretical and experimental results, the relative simplicity of the procedure and accuracy demonstrates the applicability of the proposed method in practical cases.
Appendix
The single equivalent tanh function can be calculated from (14), (15), (16), and (17) as shown in (18), (19), (20) and (21).
(13)
where the new integration constant is:
(14)
The remanence of the resultant hysteresis of the two process magnetization is described from (1b) and (2b) (for h = 0) as:
(15)
The maximum magnetization amplitude a0 is:
(16)
From (15) α0hc0 can be calculated. From the firs derivative of (2b) by h comes:
(17)
In the knowledge a0 and α0 the hc0 can also be calculated.
(18)
and
(19)
From (15), (16), (17), (18) and (19) all parameters can be calculated for
(20)
and
(21)
as functions of the field excitation h.
This calculation can be repeated n – 1 times, for the equivalent set of parameters, where n is the number of processes in magnetization.
NOTES