Custom Transformer (B-H Curve)

Open Model

This example shows calculation and confirmation of a nonlinear transformer core magnetization characteristic. Starting with fundamental parameter values, the core characteristic is derived. This is then used in a Simscape™ model of an example test circuit which can be used to plot the core magnetization characteristic on an oscilloscope. Model outputs are then compared to the known values.

Open Model

Specification of Parameters

Fundamental parameter values used as the basis for subsequent calculations:

  • Permeability of free space, $\mu_0, \rm{H/m}$

  • Relative permeability of core, $\mu_r$

  • Number of primary turns, $N_1$

  • Number of secondary turns, $N_2$

  • Effective magnetic core length, $l_e, \rm{m}$

  • Effective magnetic core cross-sectional area, $A_e, \rm{m^2}$

  • Core saturation begins, $B_{sat_{begin}}, \rm{T}$

  • Core fully saturated, $B_{sat}, \rm{T}$

Calculate Magnetic Flux Density Versus Magnetic Field Strength Characteristic

Where:

  • Magnetic flux density, $B, \rm{T}$

  • Magnetic field strength, $H, \rm{A/m}$

Linear representation:

  • $B = \mu_0 \mu_r H$

Nonlinear representation (including coefficient, a):

  • $B = B_{sat} \tanh(a.H)$

Use Parameters in Simscape Model

The parameters calculated can now be used in a Simscape model. Once simulated, the model is set to output a Simscape logging variable, simlog_CustomTransformerBHCurve, and some signals using output ports, yout. Circuit parameters are:

  • Voltage source magnitude, $V_s=10 \rm{V}$

  • Voltage source frequency, $Freq_{Hz}=60 \rm{Hz}$

  • Voltage source resistance, $R_{Vs}=10 \rm{\Omega}$

  • Operational amplifier input resistance, $R_1=1 \rm{k\Omega}$

  • Operational amplifier feedback resistance, $R_2=1 \rm{M\Omega}$

  • Operational amplifier feedback capacitance, $C_2=1 \rm{\mu F}$

Calculations on Logging and Output Data

The data needs to be processed to provide the Magnetic field strength and Magnetic flux density data for comparison. Where:

  • Magnetomotive force, $F_m, \rm{A}$

  • Magnetic flux, $\phi, \rm{Wb}$

  • Operational amplifier input voltage, $V_{in}, \rm{V}$

  • Operational amplifier output voltage, $V_{out}, \rm{V}$

Equations to be used are as follows:

  • $F_m = I.N_1$

  • $H = F_m/l_e$

  • $B = \phi/A_e$

  • $V_{out}=\frac{1}{R_1C_2}\int V_{in}.dt+c$

  • $\phi =\frac{1}{N_2}\int v.dt$

Conclusion

The three characteristics can now be overlaid:

  • Defined characteristic: calculated from fundamental parameters

  • Characteristic from logging: calculated from internal Simscape logging data

  • Characteristic from measurement: obtained by measurement and calculation using electronic test circuit

Due to leakage and parasitic parameters, the characteristic obtained from the electronic test circuit differs to the defined characteristic. However, the test circuit and its parameterization is shown to find the characteristic for the given transformer within suitable tolerances.

Results from Real-Time Simulation

This example has been tested on a Speedgoat Performance real-time target machine with an Intel® 3.5 GHz i7 multi-core CPU. This model can run in real time with a step size of 50 microseconds.