Iron losses model

Introduction

Iron losses model can be defined:
  • In the material
  • In the post-processor, in Iron losses computation box

In both cases, the defined model is only used for the computation in the post-processor. Iron losses are not taken into account in the solving process.

In Flux, the iron losses models are only available on the laminated regions, the magnetic flux density used in these models is the homogenized magnetic flux density in a block of iron sheets defined by:

B h = K f B f + 1 - K f B a

with:
  • K f : stacking factor
  • B f : material magnetic flux density (iron, steel...
  • B a : air magnetic flux density
  • B h : homogenized magnetic flux density

In the presence of rotating magnetic fields the curve By = f (Bx) can have an elliptical shape, thus it is more judicious to decompose homogenized induction vector of induction homogenize as follows:

B h = B \ \ B

with:
  • B \ \ : magnetic induction in the longitudinal direction
  • B : magnetic induction in the transversal direction

Those both axis correspond to the minor and major axes of the elliptical curve described by the Figure 1.



Figure 1. Curve By = f(Bx)

Modified Bertotti model

Presentation of the model:

The theory of Bertotti gives us the expression of the magnetic losses depending on the frequency and the magnetic flux density.

The power density is expressed with this transient magnetic relation:

d P ( t ) = K f k 1   f B \ \ +   B K f α 1 + k 2     g α 2 1 K f α 2 d B h d t α 2 + k 3   g α 3 1 K f α 3 d B h d t α 3

g α i =   2 ( 2 π ) α i - 1   4 c o s α i θ   d θ

where:
  • k 1 : is the coefficient of losses by hysteresis
  • k 2 : is the coefficient of classical Foucault currents losses
  • k 3 : is the coefficient of supplementary losses or in excess*
  • α 1 : is the exponent of losses by hysteresis
  • α 2 : is the exponent of classical Foucault currents losses
  • α 3 : is the exponent of supplementary losses or in excess*
  • f : frequency
  • B h : homogenized magnetic induction
  • B \ \ : magnetic induction in the longitudinal direction
  • B : magnetic induction in the transversal direction
  • K f : fill factor of the region

* the distinction between supplementary losses / in excess losses and classical losses is artificial. They can be grouped in one term and they therefore correspond to real induced currents flowing in the sheet.

An Excel sheet permitting to determine, coefficients and exponents (Modified Bertotti model identification tool) is available in the following path:

INSTALLATION_FOLDER\Flux\DocExamples\Tools\BertottiLossesCoefficients

in which INSTALLATION_FOLDER represents the Flux installation directory in your system.

To run the modified Bertotti model identification tool, open the file BertottiLossesCoefficients.xls available in that location in Microsoft Excel.

Limits of validity

In a Steady state AC Magnetic application, power density is expressed with the following equation:

d P = k 1   B m a x α 1   f + k 2   ( B m a x   f ) α 2 + k 3   ( B m a x   f ) α 3   K f

The software use the value of the magnetic flux density in each point. Consequently, it is convenient to be very careful with the results concerning the problems represented by the rotating machines with the Steady state AC Magnetic simulation. Indeed, for this type of simulation the rotor has a fixed position with respect to the stator, and the real rotor movement is modeled by changing the resistivity of the conductors of the rotor electric circuit. Thus, the calculated magnetic flux density is maximum in a point on the given position of the rotor in relationship with the stator, due to space harmonics. It follows that the calculated magnetic flux density does not correspond to the peak value of the magnetic flux density over a period in the time domain if the rotor were turning. Consequently, the computation of the magnetic losses must be utilized in this case with much caution.

Moreover, in the case of a non-linear approximation for the magnetic behavior law B(H), the saturation phenomenon, introduced by means of an equivalent model of magnetization, can alter the local values of the magnetic flux density.

Steps of the modified Bertotti model definition in the material

The steps of the modified Bertotti model definition in the material are as follows
  1. Create material
  2. In Iron losses tab, check Model for iron losses computation and choose Modified Bertotti model
  3. Enter coefficients and exponents

LS model

Presentation:

The LS (Loss Surface) model is a method of estimation of the magnetic losses a posteriori, based on a model of dynamic hysteresis associated to a finite elements simulation. The LS model requires that the magnetic behavior of a material be perfectly well defined, having knowledge of a characteristic surface H(B,dB/dt) (determined experimentally). Thus, for a B(t) signal of a certain shape and frequency, we can go up via the H(B,dB/dt) surface to the H(t) field, and thus reconstruct the dynamic cycle of hysteresis corresponding to it. This principle is represented below.
  • B(t) signal with common form and frequency
  • Creation of a characteristic surface H(B, dB/dt) of the material measured experimentally.
  • Reconstitution of H(t) signal (i.e. of the hysteresis cycle)
  • Calculus of losses

Characteristic surface H(B, dB/dt):

For each of the materials, the characteristic surface H(B,dB/dt) can be obtained by using a Epstein type device for magnetic measurements in medium frequency.

An example of this type of surface is represented in the figure below.


Figure 2. H(B, dB/dt) surface

Reconstruction of the cycle:

An analytical model permits the reconstruction of H field from values of B:

H(B,dB/dt) = Hstatic(B) +Hdynamic (B,dB/dt)

B(H) curves like those in the figure below can therefore be obtained, permitting the computations of iron losses quite accurately.


Figure 3. Sinusoidal magnetic flux density & 5th harmonic at 200Hz

The LS model is assigned to different nuances of laminations, which have been especially described to this purpose.

Available materials Flux:
  • M270-35A
  • M330-35A
  • M330-65A
  • M400-50A
  • M600-50A
  • M600-65A
  • M800-50A
  • M800-65A

For more informations about electrical steel sheets, see this page: Standards for electrical steel sheets

With this model, the user can also import his own materials in Flux (see LS model identification with MILS).

Definiton of the LS model in the material

See the section below: Isotropic soft material : Iron sheets described by LS model

Note:
  • Additional information on iron loss with LS model is available in the following document: « Module LS pour l’estimation des pertes fer dans les machines électriques » - A. Lebouc – Rapport d’étude mai 2004