PLECS 3.4 Online Help

Magnetic Modeling

Inductors and transformers are key components in modern power electronic circuits. Compared to other passive components they are rather difficult to model for the following reasons:

In PLECS, the user can build complex magnetic components in a special magnetic circuit domain. Primitives such as windings, cores and air gaps are provided in the Magnetics Library. The available core models include saturation and hysteresis. Frequency dependent losses can be modeled with magnetic resistances. Windings form the interface between the electrical and the magnetic domain.

Alternatively, less complex magnetic components such as saturable inductors and single-phase transformers can be modeled directly in the electrical domain.

Equivalent circuits for magnetic components

To model complex magnetic structures with equivalent circuits, three different approaches exist: Coupled-inductors, the resistance-reluctance analogy and the capacitance-permeance analogy.

Coupled inductors

In the coupled inductor approach, the magnetic component is modeled directly in the electrical domain as an equivalent circuit, in which inductances represent magnetic flux paths and losses incur at resistors. Magnetic coupling between windings is realized either with mutual inductances or with ideal transformers.

Using coupled inductors, magnetic components can be implemented in any circuit simulator since only electrical components are required. This approach is most commonly used for representing standard magnetic components such as transformers. The figure below shows an example for a two-winding transformer, where Lσ1   and L σ2   represent the leakage inductances, Lm   the non-linear magnetization inductance and Rfe   the iron losses. The copper resistances of the windings are modeled with R1   and R2  .

pict
Transformer implementation with coupled inductors

However, the equivalent circuit bears little resemblance to the physical structure of the magnetic component. For example, parallel flux paths in the magnetic structure are modeled with series inductances in the equivalent circuit. For non-trivial magnetic components such as multiple-winding transformers or integrated magnetic components, the equivalent circuit can be difficult to derive and understand. In addition, equivalent circuits based on inductors are impossible to derive for non-planar magnetic components.

Reluctance-resistance analogy

The traditional approach to model magnetic structures with equivalent electrical circuits is the reluctance-resistance analogy. The magnetomotive force (MMF) F   is regarded as analogous to voltage and the magnetic flux P   as analogous to current. As a consequence, magnetic reluctance of the flux path R  corresponds to electrical resistance:

    F
R = P-

The magnetic circuit is simple to derive from the core geometry: Each section of the flux path is represented by a reluctance and each winding becomes an MMF source.

To link the external electrical circuit with the magnetic circuit, a magnetic interface is required. The magnetic interface represents a winding and establishes a relationship between flux and MMF in the magnetic circuit and voltage v   and current i   at the electrical ports:

v = N dP
     dt

   F-
i = N

where N   is the number of turns. If the magnetic interface is implemented with an integrator it can be solved by an ODE solver for ordinary differential equations:

      ∫
P = 1-  v dt
    N

The schematic below outlines a possible implementation of the magnetic interface in PLECS.

pict
Implementation of magnetic interface

Although the reluctance-resistance duality may appear natural and is widely accepted, it is an awkward choice for multiple reasons:

Permeance-capacitance analogy

To avoid the drawbacks of the reluctance-resistance analogy the alternative permeance-capacitance analogy is most appropriate. Here, the MMF F   is again the across-quantity (analogous to voltage), while the rate-of-change of magnetic flux P   is the through-quantity (analogous to current). With this choice of system variables, magnetic permeance P  corresponds to capacitance:

      dF-
P = P  dt

Hence it is convenient to use permeance P  instead of the reciprocal reluctance R  to model flux path elements. Because permeance is modeled with storage components, the energy relationship between the actual and equivalent magnetic circuit is preserved. The permeance value of a volume element is given by:

P =  1-= μ0μrA-
     R      l

where            -7    2
μ0 = 4π × 10 N∕A   is the magnetic constant, μr   is the relative permeability of the material, A   is the cross-sectional area and l   the length of the flux path.

Magnetic resistors (analogous to electrical resistors) can be used in the magnetic circuit to model losses. They can be connected in series or in parallel to a permeance component, depending on the nature of the specific loss. The energy relationship is maintained as the power

Ploss = F P

converted into heat in a magnetic resistor corresponds to the power loss in the electrical circuit.

Windings form the interface between the electrical and the magnetic domain. A winding of N   turns is described with the equations below. The left-hand side of the equations refers to the electrical domain, the right-hand side to the magnetic domain.

v = N P

i = F
   N

Because a winding converts through-quantities (P   resp. i  ) in one domain into across-quantities (v   resp. F  ) in the other domain, it can be implemented with a gyrator, in which N   is the gyrator resistance R  . The figure below shows the symbol for a gyrator and a possible implementation in PLECS.

pict pict
Gyrator symbol and implementation

In principle, the gyrator component could be used with regular capacitors to build magnetic circuits. However, neither the gyrator symbol nor the capacitor adequately resemble a winding respectively a flux path. Moreover, any direct connection between the electrical and magnetic domain made by mistake would lead to non-causal systems that are very difficult to debug. Therefore, dedicated magnetic components should be used when modeling magnetic circuits.

Magnetic Circuit Domain in PLECS

The magnetic domain provided in PLECS is based on the permeance-capacitance analogy. The magnetic library comprises windings, constant and variable permeances as well as magnetic resistors. By connecting them according to the physical structure the user can create equivalent circuits for arbitrary magnetic components. The two-winding transformer from above will look like the schematic below when modeled in the magnetic domain.

pict
Transformer implementation in the magnetic domain

Pσ1   and Pσ2   represent the permeances of the leakage flux path, Pm   the non-linear permeance of the core, and Gfe   dissipates the iron losses. The winding resistances R1   and R2   are modeled in the electrical domain.

Modeling Non-Linear Magnetic Material

Non-linear magnetic material properties such as saturation and hysteresis can be modeled using the variable permeance component. The permeance is determined by the signal fed into the input of the component. The flux-rate through a variable permeance P(t)   is governed by the equation:

    d             dF   d
P = dt (P ⋅F) = P ⋅dt-+ dtP ⋅F

Since F   is the state variable the equation must be solved for dFdt-  . Therefore, the control signal must provide the values of both P(t)   and ddtP(t)  .

The control signals must also provide the flux P(t)   through the permeance. This enables the solver to enforce Kirchhoff's current law for all branches k   of a node:

 n
∑
   Pk = 0
k=1

When specifying the characteristic of a non-linear permeance, we need to distinguish carefully between the total permeance Ptot(F) = P ∕F   and the differential permeance Pdiff(F) = dP∕dF  .

If the total permeance Ptot(F )   is known the flux-rate P   through a time-varying permeance is calculated as:

P  =   dP-
       dt
   =   d-(P   ⋅F)
       dt  tot
   =   Ptot ⋅ dF-+ dPtot⋅F
            dt     dt
   =   Ptot ⋅ dF-+ dPtot⋅ dF ⋅F
       (    dt    dF  ) dt
   =    Ptot + dPtot⋅F  ⋅ dF
               dF        dt

In this case, the control signal for the variable permeance component is:

⌊  P(t)  ⌋  ⌊        d--      ⌋
⌈ -d     ⌉ = ⌈ Ptot + dF Ptot ⋅F ⌉
  dtP(t)             0
   P(t)            Ptot ⋅F

In most cases, however, the differential permeance P   (F )
 diff   is provided to characterize magnetic saturation and hysteresis. With

       dP-
P  =   dt
       dP- dF-
   =   dF ⋅ dt
   =   P   ⋅ dF  ,
        diff  dt

the control signal is

⌊  P(t)  ⌋   ⌊  P     ⌋
⌈ -dP(t) ⌉ = ⌈   d0iff  ⌉
  dtP(t)        P   ⋅F
                tot

Saturation Curves for Soft-Magnetic Material

Curve fitting techniques can be employed to model the properties of ferromagnetic material. As an example, a saturation curve adapted from the modified Langevian equation for bulk magnetization without interdomain coupling is used, which is referred to as the coth function:

        (             )
              3H-  -a-
B = Bsat  coth  a - 3H   + μsatH

The coth function has three degrees of freedom which are set by the coefficients Bsat  , a   and μsat  . These coefficients can by found e.g. using a least-squares fitting procedure. Calculating the derivate of B   with respect to H   yields

          (    2              )
dB- = Bsat tanh-(H-∕a)--1- -a-  + μsat
dH          a tanh2 (H∕a)   H2

With the relationships P = B ⋅A   and F = H  ⋅l   the control signal Pdiff   for the variable permeance is easily derived from the equation above.

References

S. El-Hamamsy and E. Chang, "Magnetics modeling for computer-aided design of power electronics circuits," in Power Electronics Specialists Conference, vol. 2, pp. 635-645, 1989.
R. W. Buntenbach, "Improved circuit models for inductors wound on dissipative magnetic cores," in Proc. 2nd Asilomar Conf. Circuits Syst., Pacific Grove, CA, Oct. 1968, pp. 229-236 (IEEE Publ. No. 68C64-ASIL).
R. W. Buntenbach, "Analogs between magnetic and electrical circuits," in Electron. Products, vol. 12, pp. 108-113, 1969.
D. Hamill, "Lumped equivalent circuits of magnetic components: the gyrator-capacitor approach," in IEEE Transactions on Power Electronics, vol. 8, pp. 97-103, 1993.
D. Hamill, "Gyrator-capacitor modeling: A better way of understanding magnetic components," in APEC Conference Proceedings pp. 326-332, 1994.