Three dimensional optimization of electromagnetic devicesbased on

3D Optimization of Ferrite Inductor
Considering Hysteresis Loss
Hokkaido University:
Muroran Institute of Technology:
Taiyo Yuden Co.:
Kyoto University:
Fujitsu Ltd.:
○ T. Sato, H. Igarashi
K. Watanabe
K. Kawano, M. Suzuki
T. Matsuo, K. Mifune
U. Uehara, A. Furuya
1. Introduction
A)
B)
Background
Purpose
2. Coupled analysis
3. 3D shape optimization of Inductor
4. Conclusion
2
Background

1.
2.
3.
Introductions
Coupled analysis
3D optimizations
Dc-dc converters
They are used for electric and electronic devices,
Their efficiency must be improved.

The energy losses in dc-dc converters are caused in
 control ICs
 FETs and diodes
 Inductors
The purpose of this work is that the efficiency of
the dc-dc converters is improved by reducing
the inductor losses.
3
1.
2.
3.
Purpose

Introductions
Coupled analysis
3D optimizations
The losses in the inductors are




Copper loss due to the winding resistance,
Hysteresis loss due to the magnetic cores.
The inductors must be designed to reduce them simultaneously
under the condition that the inductance satisfies the specifications
for inductance.
Design optimization based on the finite element analysis is
effective to solve this problems
However, optimizations of 3D inductor models are
computationally prohibitive mainly due to time
consuming mesh generations.
 We present the 3D shape optimization of the
inductor on the basis of nonconforming voxel FEM.

4
(Introduction. B): Purpose

1.
2.
3.
Introductions
Coupled analysis
3D optimizations
In the present method, the voxels at the material
boundaries are subdivided into fine voxels.


The fine voxels at the boundaries improve the
accuracy in the voxel FEM.
This method realizes fast generation of FE
meshes.
The multiobjective optimization of the inductors is
performed to minimize the copper and hysteresis.
5
1. Introduction
2. Coupled analysis
A)
B)
C)
D)
E)
Dc-dc converter model
Inductance under dc bias
Hysteresis loss
Inductor model for test
Analysis results
3. 3D shape optimization of Inductor
4. Conclusion
6
1.
2.
3.
Dc-dc converter model

Introductions
Coupled analysis
3D optimizations
Before the optimizations, we evaluate the inductance
and hysteresis loss under the dc bias condition by
performing the coupled analysis of inductors and dc-dc
converters.

A boost circuit shown in Fig 1 is analyzed.



The steady-state mode is evaluated.
The FET and diode are treated as ideal devices,
except for their on-resistance.
The hysteresis loss is evaluated using the Steinmetz
formula in the post processing.
L=1μH
Vout=5 V
Vin=3.3V
RIC
Cin
=2.2 μF
Control logic
RFB
Cout
=1MΩ =4.7 μF
Ro
PWM signal,
3.5MHz
Fig. 1 Boost circuit.
7
1.
2.
3.
Inductance under dc bias
Introductions
Coupled analysis
3D optimizations
1.
(Circuit analysis):Assuming the value of inductance and winding resistance
in the inductor, the circuit analysis of the boost circuit is carried out.
2.
(Nonlinear FEM): The operating points in FEs corresponding to the dc bias
are computed based on the nonlinear analysis using the dc property.
3.
(Coupled analysis):The reluctivity for the operating points is determined
from the ac property. The coupled analysis is performed over one cycle
around the operating point, in which magnetic field is analyzed by the
linear FE equation using the determined reluctivity.
L=1μH
Vout=5 V
Cin
=2.2 μF
RFB
Cout
=1MΩ =4.7 μF
PWM signal,
3.5MHz
0.4
relative permeability
RIC
Control logic
Ro
0.3
B (T)
Vin=3.3V
0.2
0.1
0
0
200
400
H (A/m)
600
800
250
Reluctivity for
operating points
200
150
100
50
0
0
dc property
20000
40000
H (A/m)
60000
ac property
 K
b T
 u
bu   A  0 

,
Z   i  V 
8
1.
2.
3.
Hysteresis loss

Introductions
Coupled analysis
3D optimizations
The hysteresis loss in the inductor, Wh (W), is estimated
based on the Steinmetz formula using the coupled analysis
results.
Wh  K h f
V
e e
Be

Kh [W/(m3∙s)]=9945.7, α=2.26 : hysteresis coefficients
f : switching frequency, Ve : volume of element e
Be:amplitude of flux density in element e

Kh and α are measured under the condition that f is 3.5MHz.
L=1μH
Vout=5 V
Vin=3.3V
RIC
Cin
=2.2 μF
Control logic
RFB
Cout
=1MΩ =4.7 μF
Ro
PWM signal,
3.5MHz
Fig. 1 Boost circuit.
9
1.
2.
3.
Inductor model
Introductions
Coupled analysis
3D optimizations
To test the validity of the analysis, the efficiency of the
boost circuit in the steady-state mode is computed.
 The inductor shown in Fig. 2 is used in the circuit.

Coil, 6 turn
Ferrite core
unit: (mm)
1.2
0.55
z

Inductance:1.09(μH)

Winding resistance:0.14(Ω)
3.0
1.2
y
0.275
x
x
Fig. 2 Original model
10
Nonconforming connections
1.
2.
3.
Introductions
Coupled analysis
3D optimizations
In the FE analysis, the unknowns assigned to
nonconforming edges (slave edge) are interpolated by
those assigned to the conforming edges (master edge).
 Slave unknown, As, is expressed by the linear combination
of the master unknowns.

As 
 
Aim N im  dl
i 1
Aim:master unknowns
es
N im:interpolated function
Master
l:tangential vector of slave
Slaves
11
1.
2.
3.
Analysis results
The simulation results and measurement data are shown in Fig. 3.


The simulation results are in good agreement with the
measurement results. The efficiency in the heavy load
conditions obtained by the analysis and measurement is almost
identical.
There exist small discrepancies for the light load conditions.
This would be due to the fact that turn-on and off times in FET
are not taken into account.
The inductance and hysteresis loss, as well as the
copper loss, under the bias condition are correctly
considered
Analysis
Measurement
100
80
60
40
20
0
efficiency (%)

Introductions
Coupled analysis
3D optimizations
0
0.1
0.2 0.3 0.4 0.5
load current (A)
Fig. 3 Efficiency.
0.6
12
1.
2.
3.
Analysis results
Introductions
Coupled analysis
3D optimizations
Fig. 4 shows the computed hysteresis and copper losses.
 Hysteresis loss is dominant in the light loads, whereas the
copper loss is dominant in the heavy loads.

The inductor current is shown in Fig. 5.
 The amplitudes of the current are about 320mA being
independent of the load resistance.
30
Copper loss
180
20
120
10
60
0
0
0
0.1 0.2 0.3 0.4 0.5 0.6
load current (A)
Fig. 4 Hysteresis. and Copper losses
inductor current (A)
Hysteresis loss (mW)
Hysteresis loss
Copper loss (mW)

50Ω
1.6
16Ω
10Ω
1.2
320mA
0.8
0.4
0
0
2
4
6
8
time (μs)
10
12
14
Fig. 5 Inductor current
13
1.
2.
3.
Analysis for optimization
If the hysteresis loss is evaluated using the coupled analysis, long
computational time is necessary in the optimizations.


The inductor current is independent
of the load resistance.
We use Steinmetz formula.
inductor current (A)

Introductions
Coupled analysis
3D optimizations
0.5
320mA
0.25
0
0
2
4
6
8
time (μs)
Wh  K h f

e
12
14
Ve Be

The hysteresis loss is estimated only by the
magnetostatic analysis.
1. The operating points in FEs for a bias
current are computed.
2. Linearized FE equation around the
operating point is solved, in which the
inductor current is set to 160mA.
From the resultant flux density, the
hysteresis loss can be evaluated,
in addition to the inductance.
inductor current (A)

10
0.5
160mA
0.25
Bias current
0
0
Wh  K h f
2
4

e
6
8
time (μs)
Ve Be
10
12
14

14
1. Introduction
2. Coupled analysis
3. 3D shape optimization of Inductor
A)
B)
C)
Optimized inductor model
Objectives
Optimization results
4. Conclusion
15
1.
2.
3.
Optimized inductor model

Introductions
Coupled analysis
3D optimizations
The inductor shown in Fig. 6 is optimized.

Design variables g:five parameters
g  r1
r2
ch
w cr T
cr is the diameter of the winding wire, which
takes 0.1 or 0.15mm.
y
80
unit: (mm)
z
80
Ferrite core
w
Coil
0.6
r2
ch
r1
r1
r2 w
80
x
r1
r2
w
80
x
Fig. 6 Optimized inductor model (1/8)
16
1.
2.
3.
Objectives

Introductions
Coupled analysis
3D optimizations
Two objective functions and constraints are defined.
Objectives:
constraints:
L  LT 2
LT
2
f1  R,
 0.03,
f 2  Wh ,
Ll  0.7 LT ,
r2  w  1.5mm
R (Ω): winding resistance,Wh (mW): Hysteresis loss.
LT=1μH: specified inductance.
L (μH): inductance value when dc bias is 0.2A.
Ll (μH) : inductance value when dc bias is 2.7A.
y
80
unit: (mm)
z
80
Ferrite core
w
Coil
0.6
r2
ch
r1
r1
r2 w
80
x
r1
r2
w
80
x
17
1.
2.
3.
optimization results

The optimization is performed using the real-ceded Generic
Algorithm and Strength Pareto Evolutionally Algorithm 2.
The optimization is performed over 200 generations.


The Pareto solutions are shown in Fig. 6 in which the red and
blue makers denote the solutions with cr=0.15 and 0.1 mm,
respectively.
All Pareto solutions are superior to the original inductor.
hysretesis loss (mW)

Introductions
Coupled analysis
3D optimizations
Optimized (0.15)
12
Optimized (0.1)
Original
(A)
10
8
6
(B)
4
2
0.1
0.12
0.14
0.16
winding resistance (Ω)
Fig. 8 Pareto solutions
0.18
18
1.
2.
3.
optimization results



The resultant inductor shapes are shown below.
Shape(A) has large cross-section of the coil. The cross-section on the
cylindrical core is so small that the magnetic induction is strong here
solutions (0.15) solutions (0.1) test model
and hysteresis loss is large.
Shape (C) has small and large crosssections of the coil and cylindrical core,
respectively.
Shape (B) is in between (A) and (C). This
seems to be optimum in this case.
(A)
1.45mm
Upper core
volume:0.78mm3
3.75mm
L: 0.976μH, Ll: 0.96μH,
R: 10.5mΩ, Wh: 9.08mW
Num. turn: 3
(B)
1.37mm
hysretesis loss (mW)

Introductions
Coupled analysis
3D optimizations
12
10
8
(A)
6
(C)
(B)
4
2
Upper core
volume:0.38mm3
0.1
0.12
0.14
0.16
winding resistance (Ω)
(C)
1.15mm
Upper core
volume:0.28mm3
3mm
1.5mm
L: 0.97μH, Ll: 0.94μH,
R: 12.07mΩ, Wh: 6.67mW
Num. turn: 4
0.18
L: 0.99μH, Ll: 0.98μH,
R: 16.1mΩ, Wh: 4.8mW
Num. turn: 4
19
1.
2.
3.
optimization results
The circuit efficiency is computed for the shape (B).


Because both losses in the inductor are reduced,
the efficiency in the all loads is improved.
It can be concluded that the present method can
find the inductor to improve the efficiency of the
converter.
Solution (B)
Original model (A)
94
efficiency (%)

Introductions
Coupled analysis
3D optimizations
92
90
88
86
84
0
0.1
0.2
0.3
0.4
load current (A)
0.5
0.6
Fig. 9 Efficiency after optimization.
20
1. Introduction
2. Coupled analysis
3. 3D shape optimization of Inductor
4. Conclusion
21
Conclusions

The efficiency in the dc-dc converters must be improved.

The hysteresis and copper losses in
the inductors should be minimized.
Multiobjective optimization of 3D inductor shape is
performed based on the nonconforming voxel FEM.
 The hysteresis loss and winding resistance can be
reduced under the restriction that the inductance
satisfies the specifications.
 The circuit efficiency is improved using the resultant
inductor.

Future works

Hysteresis modeling is introduced for the accurate analysis.
22