# Modified Bertotti model identification tool

## Introduction

This document presents the modified Bertotti model identification tool available in Flux. This tool is provided to the user as a Microsoft Excel file, and uses iron loss measurements to determine the coefficients and exponents in that model that best fit the measurement set. The identification of these parameters is required to create a material containing a Bertotti model for iron loss computations in Flux.

- how to find the modified Bertotti model identification tool in Flux installation directory;

- a reminder of the modified Bertotti model for iron loss computations in Flux;
- a brief discussion on the data fitting approach implemented by the identification tool;
- how to use the identification tool.

## Location of the identification tool

The 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.

## The modified Bertotti model

The modified Bertotti model implemented in Flux separates iron losses in three different parts: hysteresis losses, classical eddy current losses and excess losses accordingly with the equation below:$dP={k}_{1}{B}_{max}^{{\alpha}_{1}}f+{{k}_{2}{(B}_{max}f)}^{{\alpha}_{2}}+{{k}_{3}{(B}_{max}f)}^{{\alpha}_{3}}$

Flux uses the expression above to evaluate the power density in AC steady state applications. In that equation:

*dP*represents the total iron losses per unit volume (in W/m^{3});*k*is the hysteresis losses coefficient;_{1}*k*is the classical Foucault losses coefficient;_{2}*k*is the excess losses coefficient;_{3}*α*is the hysteresis losses exponent;_{1}*α*is the classical Foucault losses exponent;_{2}*α*is excess (or supplementary) losses exponent;_{3}*f*is the frequency (in Hz) and*B*is the maximum magnetic flux density in an electrical period (in T)._{max}

## Input data

Table 1 contains an example of iron losses measurements provided by an electrical steel supplier. Such measurements are usually performed on an Epstein frame, by imposing a sinusoidal magnetic flux density to the material.

B_{max }(T) |
Losses at 50 Hz (W/kg) | Losses at 100 Hz (W/kg) | Losses at 200 Hz (W/kg) | Losses at 2500 Hz (W/kg) |
---|---|---|---|---|

0.1 | 0.05 | 0.04 | 0.08 | 3.89 |

0.2 | 0.06 | 0.14 | 0.32 | 14.3 |

0.3 | 0.11 | 0.30 | 0.73 | 29.6 |

0.4 | 0.20 | 0.49 | 1.21 | 50.2 |

0.5 | 0.23 | 0.71 | 1.78 | 76.7 |

0.6 | 0.38 | 0.97 | 2.44 | 110 |

0.7 | 0.50 | 1.25 | 3.19 | 153 |

0.8 | 0.62 | 1.57 | 4.03 | 205 |

0.9 | 0.77 | 1.92 | 4.97 | 270 |

1.0 | 0.92 | 2.31 | 6.01 | 349 |

1.1 | 1.10 | 2.75 | 7.19 | |

1.2 | 1.31 | 3.26 | 8.54 | |

1.3 | 1.56 | 3.88 | 10.1 | |

1.4 | 1.92 | 4.67 | 12.2 | |

1.5 | 2.25 | 5.54 | 14.4 | |

1.6 | 2.53 | |||

1.7 | 2.75 | |||

1.8 | 2.94 |

The input data required by the identification tool is similar in format to Table 1. It
consists of a list of triplets (*B _{max}*,

*L*,

*f*) relating the peak magnetic flux density

*B*(in teslas) in the material to both specific iron losses

_{max}*L*(in W/kg) and frequency

*f*(in Hz). The user also needs to provide the density or specific mass of the material

*ρ*(in kg/m

^{3}).

The goal of the identification tool is to find exponents (*α _{1}*,

*α*,

_{2}*α*) and coefficients (

_{3}*k*,

_{1}*k*,

_{2}*k*) that best fit the input data, as discussed in the next section.

_{3}*k*,

_{1}*k*,

_{2}*k*) and (

_{3}*α*,

_{1}*α*,

_{2}*α*) identified by the tool are valid for both AC steady state and Transient applications in Flux.

_{3}## Data fitting approach

To determine coefficients (*k _{1}*,

*k*,

_{2}*k*) and exponents (

_{3}*α*,

_{1}*α*,

_{2}*α*) of the modified Bertotti model for a given set of measurements, the identification tool employs the least squares minimization approach presented below.

_{3}Let *r _{ij}* be the squared residual between specific loss measurement

*m*and the specific loss predicted by the modified Bertotti model

_{ij}*b*=

_{ij}*dP*(

*k*)

_{1}, k_{2}, k_{3}, α_{1}, α_{2}, α_{3}, Bmax_{i}, f_{j}*/ ρ*, both evaluated at frequency

*f*. This residual is computed by the following expression:

_{j}*r _{ij} = *(

*m*)

_{ij}- b_{ij }

^{2}.The identification tool uses an initial guess of parameters (*k _{1},
k_{2}, k_{3}, α_{1}, α_{2},
α_{3}*) and Excel's GRG Nonlinear solver (available in the Solver
add-in) to minimize the total weighted residual:

*R = *Σ* _{j }* {

*w*Σ

_{j}*[*

_{i}*r*] },

_{ij}in which *w _{j}* are weights attributed to each frequency

*f*in the measurement set. These weights and parameters (

_{j}*k*) are constrained to real values greater or equal than zero.

_{1}, k_{2}, k_{3}, α_{1}, α_{2}, α_{3}## How to use the modified Bertotti model identification tool

As already mentioned, the modified Bertotti model identification tool is provided in the form of a Microsoft Excel file. This file contains four spreadsheets, namely:

**Identification Tool**: a spreadsheet containing empty input cells, and ready to be filled by the user.**Single Frequency Example**: an example of identification using measurements performed at a single frequency.**Multi-Frequency Example**: an example of identification using measurements performed at several frequencies.

*k*and the exponents

_{1}, k_{2}, k_{3}*α*using the

_{1}, α_{2}, α_{3}**Identification Tool**spreadsheet:

- Fill the "Measured Losses" table with measurements performed at one or several frequencies.
- Adjust the weights of the measurements for each frequency in the cell range C48:H48;
- Provide the density of the material in cell M5.
- Set an initial guess of the Bertotti coefficients and exponents in the cell range M11:M16.
- Open the Excel solver menu: [Data] → [Solver] (requires enabling the Solver Add-in in Excel).
- Configure and run the Excel solver as follows:
- "Set objective" field: set to $L$19
- Choose "To Min" option
- "By changing variable cells" field: set to $M$11:$M$16
- "Subject to the Constrains" field: set to $M$11:$M$16 >= 0
- "Select a Solving Method" field : choose the "GRG Nonlinear" solver.
- Click on "Solve"

- Find the fitted/updated coefficients in the "Modified Bertotti Model Parameters" table.
- Scroll down to verify the plots comparing the fitted Bertotti model to the measurements at each frequency.

*n*as an initial guess for the trial

*n+1*. The user may also find useful to adjust the weights iteratively between trials, by analyzing the partial residuals obtained for each frequency in the "Least Squares Fitting" table.