Professional Documents
Culture Documents
GangLei
JianguoZhu
YouguangGuo
Multidisciplinary Design
Optimization Methods for
Electrical Machines and
Drive Systems
Power Systems
More information about this series at http://www.springer.com/series/4622
Gang Lei Jianguo Zhu Youguang Guo
Multidisciplinary Design
Optimization Methods
for Electrical Machines
and Drive Systems
123
Gang Lei Youguang Guo
University of Technology Sydney University of Technology Sydney
Sydney, NSW Sydney, NSW
Australia Australia
Jianguo Zhu
University of Technology Sydney
Sydney, NSW
Australia
Electrical machines and drive systems account for about 46 % of all global elec-
tricity consumption, resulting in about 6,040 Mt of CO2 emissions. This is by far
the largest portion of electricity use, easily outstripping lighting, which takes up to
19 % of the worlds demand. Therefore, the energy efciency of electrical drive
systems is very important for the energy conservation, environment and sustainable
development of the world.
Electrical drive systems are key components in many modern appliances, as well
as industry equipment and systems. In order to achieve the best design objectives,
such as high performance and low cost, various optimization methods have been
developed for design optimization of electrical machines and drive systems. The
traditional design optimization is at the component level, e.g. optimization of a
motor design or the parameters of a control algorithm. However, modern appliances
or systems demand that the drive systems be specically designed and optimized to
provide full support to their best functionalities with multiple performance indi-
cators. For such applications, the authors developed an application-oriented
multi-objective system-level design optimization method. Because of the com-
plexity of drive system design that involves many disciplines, such as electro-
magnetics, materials, mechanical dynamics including structural, thermal, and
vibrational analyses, power electronic convertors, and control algorithms, a
multi-level optimization method was developed by the authors to improve the
effectiveness of the optimization of electrical machines as well as drive systems.
On the other hand, the real quality of motors and drives in mass production
highly depends on the available machinery technology and those unavoidable
variations or uncertainties in the manufacturing process, assembly process and
operation environment. The manufacturing precision and tolerances are two main
issues in the manufacturing process, including mainly the variations of material
characteristics, such as magnetization faults in terms of magnitude and magneti-
zation direction for permanent magnets (PMs), and density and permeability of
soft-magnetic-composite (SMC) stator cores manufactured by powder metallic
moulding technology, and dimensional variations of parts of drive systems, such as
v
vi Preface
the rotor, stator, winding and PMs. The assembly process variations mainly include
the lamination of silicon steel sheets and misalignments of stator, rotor and PMs.
The operating uncertainties mainly include the load variations, changes of electrical
and mechanical parameters, such as the changes of resistance and inductance due to
the operational temperature rise, and fluctuations of drive voltage.
Limited by these variations in the practical machinery technology, an aggres-
sively optimized design may be difcult for high-quality batch production and end
up with high rejection rates. Similarly, variations in system parameters and oper-
ational conditions may also lead to sub-optimal performance, and in a severe case,
even instability. To solve this type of problems, the methodology of Six-Sigma
quality control can be adopted to develop a robust design optimization method to
guarantee the high-quality batch production of drive systems.
Based on many years of research experience of the authors, this book aims to
present efcient application-oriented, multi-disciplinary, multi-objective, and
multi-level design optimization methods for advanced high-quality electrical drive
systems. The multi-disciplinary analysis includes materials, electromagnetics,
thermotics, mechanics, power electronics, applied mathematics, machinery tech-
nology, and quality control and management.
This book will benet both researchers and engineers in the eld of motor and
drive design and manufacturing, thus enabling the effective development of the
high-quality production of innovative, high-performance drive systems for chal-
lenging applications, such as green energy systems and electric vehicles.
This book consists of eight chapters, based on our several research projects, and
covering the aspects of electrical machines, drive systems, high-quality mass pro-
duction and application-oriented design optimization methods.
Like most books, this book starts with an introduction in Chap. 1 to provide an
overview of application elds of electrical machines and drives as well as the
state-of-art design optimization methods for electrical machines, drive systems and
high-quality mass production.
Chapter 2 presents an overview of the design fundamentals of electrical
machines and drive systems. Design analysis models in terms of different disci-
plines (domains) are investigated in this chapter, such as the analytical models or
methods for electromagnetic and thermal analyses, magnetic circuit model for
electromagnetic analysis, nite element model (FEM) for all electromagnetic,
thermal and mechanical analyses, and eld-oriented control and direct torque
control algorithms for the control systems.
Chapter 3 reviews the popular optimization algorithms and approximate models
used in the optimization of electrical machines as well as electromagnetic devices.
Optimization algorithms include classical gradient-based algorithms and modern
intelligent algorithms, such as genetic algorithms, differential evolution algorithm
and multi-objective genetic algorithms. Approximate models include response
surface model, radial basis function model and Kriging model.
Chapter 4 presents the design optimization methods for electrical machines in
terms of different optimization situations, including low- and high-dimensional,
single and multi-objectives and disciplines. Five new types of design optimization
Preface vii
1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..... 1
1.1 Energy and Environment Challenges . . . . . . . . . . . . . . . ..... 1
1.2 Introduction of Electrical Machines, Drive Systems,
and Their Applications . . . . . . . . . . . . . . . . . . . . . . . . . ..... 3
1.2.1 General Classication of Electrical Machines. . . . ..... 3
1.2.2 Electrical Machines and Applications . . . . . . . . . ..... 4
1.3 The State-of-Art Design Optimization Methods
for Electrical Machines and Drive Systems . . . . . . . . . . . ..... 8
1.3.1 Design Optimization of Electrical Machines . . . . ..... 8
1.3.2 Design Optimization of Electrical Drive Systems . ..... 11
1.3.3 Design Optimization for High Quality
Mass Production . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
1.4 Major Objectives of the Book . . . . . . . . . . . . . . . . . . . . . . . . . 18
1.5 Organization of the Book . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2 Design Fundamentals of Electrical Machines and Drive Systems . . . 25
2.1 Introduction. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
2.1.1 Framework of Multi-disciplinary Design . . . . . . . . . . . . 25
2.1.2 Power Losses and Efciency . . . . . . . . . . . . . . . . . . . . 26
2.2 Electromagnetic Design . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
2.2.1 Analytical Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
2.2.2 Magnetic Circuit Model . . . . . . . . . . . . . . . . . . . . . . . 30
2.2.3 Finite Element Model . . . . . . . . . . . . . . . . . . . . . . . . . 33
2.3 Thermal Design . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
2.3.1 Thermal Limits in Electrical Machines . . . . . . . . . . . . . 35
2.3.2 Thermal Network Model . . . . . . . . . . . . . . . . . . . . . . . 36
2.3.3 Finite Element Model . . . . . . . . . . . . . . . . . . . . . . . . . 41
2.4 Mechanical Design. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
2.5 Power Electronics Design . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
ix
x Contents
AC Alternating current
ANN Articial neural network
ANOVA Analysis of variance
BLDC Brushless direct current motor
CCD Central composite design
CEMPE Centre for Electrical Machines and Power Electronics
COP Coarse optimization process
CSIRO Commonwealth Scientic and Industrial Research Organization
CSRBF Compactly supported radial basis function
DC Direct current
DEA Differential evolution algorithm
DFSS Design for Six-Sigma
DOE Design of experiments
DPMO Defects per million opportunities
DTC Direct torque control
EAs Evolutionary algorithms
EDA Estimation of distribution algorithms
EMF Electromotive force
FEA Finite element analysis
FEM Finite element model
FOC Field-oriented control
FOP Fine optimization process
FSPMM Flux-switching permanent magnet machine
GA Genetic algorithm
GCD Greatest common divisor
GEVIC Green Energy and Vehicle Innovations Centre
HEVs Hybrid electric vehicles
HTS High-temperature superconductor
LCM Least common multiple
LSA Local sensitivity analysis
LSM Least square method
MCA Monte Carlo analysis
MDO Multi-disciplinary design optimization
xiii
xiv Abbreviations
In an electrical drive system, the role of electric motor is to convert electrical power
into mechanical power. It is found that electric motors account for about 46.2 % of
all global electricity consumption, leading to about 6,040 Megatons (Mt) of CO2
emissions. This is by far the largest portion of electricity use, as shown in Fig. 1.1.
Around the world, over 300 million motors are being used in industry, large
buildings and infrastructure, and about 30 million new electric motors are sold each
year for industrial purposes alone [1, 2].
Figure 1.2 illustrates the estimated electricity demand for all motors by sector. It
can be seen that the motors used in industry consume about 63.1 % of the total
energy consumption. The corresponding energy costs are estimated to be USD 362
billion per year. In the industrial sector, motors are used primarily for four areas of
applications, namely pumps, fans, compressors, and mechanical movement. These
applications and their respective shares are illustrated in Fig. 1.3 [1].
3.2%
18.6%
3.2%
19.9%
16%
25%
In terms of the life cycle cost of motors, the electric energy cost accounts for
more than 90 % of all cost in general, which are much larger than the other two
parts, purchase-price and repair or maintenance cost. As an example, Fig. 1.4
depicts the breakdown of the life cycle cost for an 11 kW motor operated 4000 h per
year, where the electric energy cost accounts for 96.7 % of the total cost [2].
1.1 Energy and Environment Challenges 3
96.7%
Therefore, the energy efciency of motors is a crucial issue for the energy
conservation, environment, and sustainable development of the world, and this is
also the main reason that high efciency motors have attracted so much attention all
over the world. Even 1 % increase in motor efciency would save about 20 billion
kWh per year or USD 1.4 billion in electricity and 3.5 million barrels of oil in the
U.S. alone. These savings would be multiplied by about a factor of four on a
worldwide basis [3, 4].
In general, there are three main kinds of applications for electrical machines, which
are electricity generation, electricity transformation, and electrical drives. The fol-
lowing are some examples.
A. Electricity generation
In most electricity generation systems, except the photovoltaic systems in which the
solar energy is converted directly into DC currents by static solar cells, various
types of energy resources, such as fossil fuels (coal, diesel, and natural gas, etc.),
water, wind, sun light, and atomic energy, etc., are rstly converted into mechanical
energy by rotating turbines and then electricity by AC rotating electrical generators.
Figure 1.5 shows the working principle of a hydroelectric generation with a
synchronous generator. Hydroelectric plants use the energy from water to power a
process that turns water potential energy into electricity. This process involves the
water flowing from the dam, through a tunnel which leads to a turbine. Once the
water reaches the turbine, the force from the water spins a generator to generate
electricity. The generator terminal is connected to a transformer, which is where the
electricity generated is transformed, e.g. to high voltage, and for long distance
transmission [5].
Wind turbines as another application of AC generators have been employed
worldwide. In general, both synchronous and induction machines are commonly
employed for wind power generation. Moreover, various rotating AC generators are
also commonly used for electricity generation from other energy sources, such as
nuclear and solar thermal power plants, gas red turbines, and diesel/petrol engines,
to meet different needs, such as power supply for grid, air planes, trains, and ships.
B. Electricity transformation
Electricity transformers can be dened as a type of static electrical device which
transfers power from one circuit to another by means of electromagnetic induction.
1.2 Introduction of Electrical Machines, Drive Systems, and Their Applications 5
(b) (c)
(a) (b)
Fig. 1.7 a Aurora solar car and b its in-wheel permanent magnet motor
1.2 Introduction of Electrical Machines, Drive Systems, and Their Applications 7
Figure 1.7 shows a picture of the Aurora solar car and a photo of the in-wheel
permanent magnet motor developed by UTS CEMPE and CSIRO. The ratings of
the motor are 5.5 kW, 50 Nm peak for 72 s, with the maximum efciency of
98.5 %.
In-wheel motors and other PM motors are widely used as drive machines in
(plug-in) hybrid electric vehicles (HEVs). In a plug-in hybrid car, the battery bank
is charged by the grid power supply when the car is not in use, and the electrical
motor plays the major role of drive. A small internal combustion engine is
employed to provide extra torque when the car accelerates, or to charge the battery
when the state of charge is low. Since the motor controlled by a power electronic
inverter can operate in all four quadrants of the torque-speed plane, the car is able to
retrieve the kinetic energy by regenerative braking when it is decelerating.
Therefore, hybrid electrical cars have much higher energy efciency than the tra-
ditional internal combustion engine drive cars.
Currently, GEVIC researchers are designing several PM machines including
flux-switching machine for plug-in HEVs. The designed and fabricated machines
will be tested in the powertrain testing facility at UTS Automotive Laboratory, as
depicted in Fig. 1.8. This facility can simulate urban and highway drive cycles and
measure the torque and speed along with the powertrain performance in various
operation modes, including regenerative braking [10].
Figure 1.9 shows the photos of the SolarSailor boat powered by 2 40 kW, 400
Nm, 3 phase, 16 pole, 950 rev/min, 100 V, 250 A direct drive high efciency PM
brushless DC motors and the power electronic controller developed by
UTS CEMPE in 2000.
For propulsion and power supply of large ships, multi-MW electrical machines
are used. In such cases, the efciency, volume and weight of the electrical machines
become a serious concern. With the technological breakthrough, large capacity
HTS generators and motors are being built around the world for application in large
cargo ships and gun boats.
8 1 Introduction
(a)
(b) (c)
Fig. 1.9 a The SolarSailor boat, b power electronic converter, and c PM drive motor
which is a challenge to both research and industry communities. This work consists
of two parts: design analysis and performance optimization.
A. Design analysis
For design analysis, electromagnetic, thermal, mechanical design analyses and
manufacturing design are the main concerns. The electromagnetic design is con-
ducted by using mainly the analytical model, magnetic circuit model, and nite
element model (FEM) to calculate the electromagnetic parameters of the machines
being designed, such as flux linkage, back electromotive force (EMF), inductance
and core losses. These parameters are then used to evaluate the machine perfor-
mance indicators, such as output power and efciency. The thermal design is
conducted by using mainly the FEM and thermal network model to compute the
temperature-rises in the machine. The mechanical design is often done by stress
and/or modal analyses using the FEM to calculate the maximum stress, deforma-
tion, and resonant frequency of the machine structure under various operating
conditions. In general, the thermal and mechanical analyses are usually conducted
to check the insulation, magnetic strength (e.g. the Curie temperature of PMs), and
mechanical material and structural strengths to ensure the machines safety. In
electrical machine design optimization, they are often used as design constraints.
B. Performance optimization
Performance optimization includes two aspects as well, namely optimization
models and optimization methods. There are several popular optimization models.
For example, from the perspective of objective numbers, optimization models can
be classied as single-objective or multi-objective models. Generally, the cogging
torque, torque ripples, material and manufacturing as well as sometimes operating
costs, weight, and energy consumption are the main concerns in the design and
optimization process. From the industrial perspective, the optimization models can
consist of three main types: the deterministic, reliability and robust models.
Regarding the optimization methods, despite many kinds of optimization methods
have been developed, their effective application for design optimization of electrical
machines and systems has always been a research focus in electrical engineering.
Since 1987, it has been selected as one of the most important development directions
in computational electromagnetics by the premier International Conferences on
Magnetics, such as Intermag (International Magnetics Conferences), CEFC
(Conference on Electromagnetic Field Computation) and Compumag (Conference on
the Computation of Electromagnetic Fields). In CEFC 2000, a special academic
lecture about these problems was organized [11]. Electrical machines and drives are
also an important section in several international conferences on electrical machines
and energy systems, such as ICEMS (International Conference on Electrical
Machines and Systems), ECCE (IEEE Energy Conversion Congress and Exposition)
and IECON (Conference of the IEEE Industrial Electronics Society). Recently, the
research interest in electromagnetic optimization design problems, particularly for
electrical machines, has increased signicantly. A special section on optimal design of
electrical machines was presented in the IEEE Transactions on Energy Conversion in
Sep. 2015 [12]. The main driving forces behind this interest are the rapid development
10 1 Introduction
Electrical machines and the corresponding drive systems have a history of over a
century and the design procedure has become almost standard. When designing an
appliance that needs an electrical drive system, the designer rstly selects the motor,
inverter/converter and controller from the existing products. The appliance
designer, on one hand, has to deliver the functions that the appliance is supposed to
have, and on the other hand, has to take into account the availability and perfor-
mance that the existing motor drive can provide. This motor manufacturer-oriented
approach has been the dominant design concept for drive systems for a long time.
However, this approach would apply many constraints to the design and therefore
limit the functions of the appliance [39].
With the fast development of numerical eld analysis, CAD software, and
flexible mechanical manufacturing technology, it is possible to design and manu-
facture a motor to meet the special requirements of a particular application such that
the designer can concentrate on pursuing the best appliance functions. Since early
1990s, this application-oriented approach has been gradually becoming a common
practice [4048]. In many cases, the motor and control electronics are closely
integrated into the appliances.
For example, the solar powered deep well submersible pump drive as shown in
Fig. 1.6 that UTS CEMPE developed in 1991 has an integrated structure. The high
efciency motor and its electronic controller are packed into the pump and installed
down the deep well. The concept of integrated design can be better illustrated by the
12 1 Introduction
in-wheel motor for the Aurora Solar Car drive (Fig. 1.7) developed jointly by
UTS CEMPE and CSIRO in 1997. In this drive, to meet the special requirement of
extremely low weight and high efciency, a core-less in wheel motor topology was
employed.
No doubt, the application-oriented integrated design concept is very advanced,
but the design methodology used in these examples and also by all other motor
designers is still very traditional. As illustrated in Fig. 1.10, the traditional design is
conducted on the component level, i.e. the motor and the electronic controller are
separately designed by the standard procedure. The major part of the motor design
is the electromagnetic design whereas the thermal and sometimes mechanical
analyses are carried out as verication only but not coupled to the electromagnetic
analysis. In the case of design optimization, inaccurate circuit, eld, and material
models, and even empirical formulae are commonly used in order to reduce the
computing time. This approach has two problems: (a) the real optimum design is
not possible because of the inaccurate models used, and (b) the system performance
cannot be optimized because the design is on the component level.
Meanwhile, the system performance optimization is becoming essential, espe-
cially when new materials, for example, the SMC, and new topologies are
employed. SMC is a relatively new soft magnetic material developed for extremely
low cost motor manufacturing using the highly matured powder metallurgical
technology [4953]. However, the magnetic properties of the SMC material are
much poorer than those of the traditionally used silicon sheet steels. In order to
develop low cost high performance SMC motor drive systems, we must explore
new motor topologies of 3D magnetic flux and new drive schemes, and optimize
the design at the system level.
The electric vehicles and HEVs are attracting great attentions and funding from
the governments and general public around the world because of the worldwide
fossil fuel energy crisis and severe greenhouse gas emissions of the conventional
vehicles powered by the internal combustion engines. To improve the efciency
and drive performance with reduced volume, weight, and cost of novel drive sys-
tems to meet the challenging requirements of hybrid electric vehicles, a great
amount of recent efforts are being directed towards the development and optimum
design of high performance drive systems for (plug-in) HEVs [54, 55].
Through the extensive research practice, it is recognized that when designing
such an electrical drive system, it is important to pursue the optimal system per-
formance rather than the optimal components like motors, because assembling
individually optimized components into a system would not necessarily guarantee
an optimal system performance. The optimal system performance can only be
achieved through a holistic approach of integrated simultaneous optimization of all
components at the system level [39].
Figure 1.11 shows a brief design framework and the coupled relations of an
electrical drive system. As shown, design optimization of such a system is a
multi-disciplinary, multi-level, multi-objective, and high dimensional problem. It
mainly includes electromagnetic, material, mechanical, thermal, and power elec-
tronic designs, which are strongly coupled [38, 39, 56, 57].
1.3 The State-of-Art Design Optimization Methods 13
Fig. 1.11 Multi-disciplinary and multi-level design framework of electrical drive systems
The design optimization method mentioned above are all deterministic design
optimization approaches which do not take into account the unavoidable variations
(noise factors) in the engineering manufacturing, including mainly the material
diversity, manufacturing error and assembly inaccuracy, and system parameter
variations in practical operation environment [6668]. The motor and drive system
1.3 The State-of-Art Design Optimization Methods 15
Fig. 1.12 Low cost mold, press, and molded core before and after curing
16 1 Introduction
both mechanically and magnetically, so that they are natural choices for the design
of motors requiring 3D magnetic flux paths. Unlike the laminated cores made of the
traditional silicon steel sheets, SMC cores can be manufactured by compacting
SMC powders in a mold, and thus suitable for constructing motors of complex
structures.
Figure 1.12 also shows a photo of the mold, press in the lab, and molded core of
a three phase PM claw pole motor before (white) and after (black) thermal curing.
On the other hand, because of specic nature of iron powders, the pressing must be
done in multiple steps in order to obtain uniform powder distribution.
The manufacturing cost of SMC cores is directly related to the size of press used
for the molding, while the productivity is inversely proportional to the press size.
For a given SMC core, it is desirable to choose a smaller press in order to keep the
manufacturing cost low. In the case where the volume of the motor is not a big
problem, in order to reduce further the manufacturing cost, it is possible to use a
low density SMC core. For example, for a 100 ton press, it can produce 500 pieces
per hour with a cost of $ 100, resulting in the manufacturing cost of $ 0.2 per piece.
For a 500 ton press, on the other hand, it can produce only 100 pieces per hour with
a cost of $ 500, and thus the manufacturing cost is $ 5 per piece [34, 6971].
Therefore, the manufacturing process is a major parameter in design optimization of
SMC motors.
B. Assembly process variations
The assembly process variations mainly include the lamination of silicon steel
sheets and misplacements of stator, rotor and PMs. Severe misplacements can result
in big variations of the motor quality and cause large vibration and excessive
resistive torque and mechanical power loss. Table 1.1 lists several manufacturing
tolerances (such as magnet strength and skew error) and assembly variations (such
as magnet disposition and rotor/stator eccentricity) obtained from industrial motor
manufacturing experiences [66].
In the assembling process, due to the structural asymmetry (such as keyway and
tag hole), non-uniform material quality (such as the thickness or sand hole) and
manufacturing error (such as drill hole and others), and big mechanical disequi-
librium in the rotating parts (such as rotors and fans) may appear, and the rotating
parts will displace from their gravity centre, resulting in unbalanced centrifugal
force and causing the motor to vibrate. Vibration has large negative effects to the
motor, such as extra energy consumption, efciency reduction, direct damage to the
shaft, acceleration of abrasion, which shortens signicantly the lifetime of the motor
and drive system.
The assembly of stator cores can also be a big challenge for SMC motors. For
manufacturing convenience, the stator core of an SMC motor is often molded in
separate pieces. Any extra air gap between two pieces of the stator core caused by
poor assembly process will result in large reduction of air gap flux density and in
turn the motor efciency. Good assembly structures should be investigated in the
design stage of SMC motors.
Figure 1.13 illustrates a multi-disciplinary design framework of electrical motors
and drive systems that we propose to take into account the manufacturing quality in
mass production in the design stage. The rst step is to dene the acceptable
maximal defect-ratio and system performance. Under these specications, the
motor types, topologies, materials, inverters, and controllers will be designed and
optimized under the multi-disciplinary design and machinery technology design.
Due to these manufacturing tolerances, the design optimization of electrical
drive systems for mass production is really a challenge in both research and
industrial communities as it includes not only the theoretical multi-disciplinary
design and analysis but also the practical engineering manufacturing of electrical
machines and drive systems. Meanwhile, many new control algorithms, e.g. MPC,
have been proposed for motor control, and in the design optimization stage, various
many algorithm parameters should be optimized for the best drive system dynamic
performance. From the industrial application perspective, it is a natural requirement
that the obtained optimal control algorithm parameters are robust against the
variations of motor parameters. This is a crucial issue for the batch production of
novel drive systems [6971].
To nd effective ways to deal with this problem, several robust design opti-
mization methods have been investigated, such as Taguchi method [7277] and
Six-Sigma robust optimization method [6971, 78]. These two methods have been
found useful to optimize motor performances (including torque ripples, cost, and
output power) and quality against the manufacturing tolerances.
Fig. 1.13 Proposed design framework of drive systems for mass production
18 1 Introduction
References
12. Fahimi B, Mohammed O (2015) Optimal design of electrical machines. IEEE Trans Energy
Convers 30(3):1143
13. Hasanien HM, Abd-Rabou AS, Sakr SM (2010) Design optimization of transverse flux linear
motor for weight reduction and performance improvement using RSM and GA. IEEE Trans
Energy Conver 25(3):598605
14. Hasanien HM (2011) Particle swarm design optimization of transverse flux linear motor for
weight reduction and improvement of thrust force. IEEE Trans Ind. Electron 58(9):40484056
15. Storn R, Price K (1997) Differential evolution- a simple and efcient heuristic for global
optimization over continuous spaces. J Global Optim 11:341359
16. Lei G, Zhu JG, Guo YG, Shao KR, Xu W (2014) Multiobjective sequential design
optimization of PM-SMC motors for six sigma quality manufacturing. IEEE Trans Magn 50
(2), Article 7017704
17. Yao D, Ionel DM (2013) A review of recent developments in electrical machine design
optimization methods with a permanent magnet synchronous motor benchmark study. IEEE
Trans Ind Appl 49(3):12681275
18. Ishikawa T, Tsukui Y, Matsunami M (1999) A combined method for the global optimization
using radial basis function and deterministic approach. IEEE Trans Magn 35(3):17301733
19. Lei G, Yang GY, Shao KR, Guo YG, Zhu JG, Lavers JD (2010) Electromagnetic device
design based on RBF models and two new sequential optimization strategies. IEEE Trans
Magn 46(8):31813184
20. Wang LD, Lowther DA (2006) Selection of approximation models for electromagnetic device
optimization. IEEE Trans Magn 42(2):12271230
21. Lebensztajn L, Marretto CAR, Costa MC, Coulomb J-L (2004) Kriging: a useful tool for
electromagnetic device optimization. IEEE Trans Magn 40(2):11961199
22. Mendes MHS, Soares GL, Coulomb J-L, Vasconcelos JA (2013) Appraisal of surrogate
modeling techniques: a case study of electromagnetic device. IEEE Trans Magn 49(5):1993
1996
23. Alotto P, Baumgartner U, Freschi F, Kstinger A, Magele Ch, Renhart W, Repetto M (2008)
SMES optimization benchmark extended: introducing Pareto optimal solutions into TEAM22.
IEEE Trans Magn 44(6):10661069
24. Guimaraes FG, Campelo F, Saldanha RR, Igarashi H, Takahashi RHC, Ramirez JA (2006) A
multiobjective proposal for the TEAM benchmark problem 22. IEEE Trans Magn 42(4):1471
1474
25. Lebensztajn L, Coulomb JL (2004) TEAM workshop problem 25: A multiobjective analysis.
IEEE Trans Magn 40(2):14021405
26. Deb K, Pratap A, Agarwal S, Meyarivan T (2002) A fast and elitist multi-objective genetic
algorithm: NSGA-II. IEEE Trans Evol Comput 6(2):182197
27. Reyes-Sierra M, Coello CAC (2006) Multi-objective particle swarm optimizers: A survey of
the state-of-the-art. Int J Comput Intell Res 2(3):287308
28. Dias AHF, Vasconcelos JA (2002) Multiobjective genetic algorithms applied to solve
optimization problems. IEEE Trans Magn 38(2):11331136
29. dos Santos Coelho L, Alotto P (2008) Multiobjective electromagnetic optimization based on a
nondominated sorting genetic approach with a chaotic crossover operator. IEEE Trans Magn
44(6):10781081
30. Ho SL, Yang SY, Ni GZ, Lo EWC, Wong HC (2005) A particle swarm optimization-based
method for multiobjective design optimizations. IEEE Trans Magn 41(5):17561759
31. Xie DX, Sun XW, Bai BD, Yang SY (2008) Multiobjective optimization based on response
surface model and its application to engineering shape design. IEEE Trans Magn 44(6):1006
1009
32. Di Barba P (2010) Multiobjective shape design in electricity and magnetism. Lecture notes in
electrical engineering, vol. 47
33. Lei G, Shao KR, Guo YG, Zhu JG (2012) Multi-objective sequential optimization method for
the design of industrial electromagnetic devices. IEEE Trans Magn 48(11):45384541
22 1 Introduction
34. Lei G, Guo YG, Zhu JG et al (2015) Techniques for multilevel design optimization of
permanent magnet motors. IEEE Trans Energy Conver 30(4):15741584
35. Lei G, Guo YG, Zhu JG, Chen XM, Xu W (2012) Sequential subspace optimization method
for electromagnetic devices design with orthogonal design technique. IEEE Trans Magn 48
(2):479482
36. Lei G, Xu W, Hu JF, Zhu JG, Guo YG. Shao KR (2014) Multilevel design optimization of a
FSPMM drive system by using sequential subspace optimization method. IEEE Trans Magn
50(2), Article 7016904
37. Wang SH, Meng XJ, Guo NN, Li HB, Qiu J, Zhu JG et al (2009) Multilevel optimization for
surface mounted PM machine incorporating with FEM. IEEE Trans Magn 45(10):47004703
38. Lei G, Liu CC, Guo YG, Zhu JG (2015) Multidisciplinary design analysis for PM motors with
soft magnetic composite cores. IEEE Trans Magn 51(11), Article 8109704
39. Lei G, Wang TS, Guo YG, Zhu JG, Wang SH (2014) System level design optimization
methods for electrical drive systems: deterministic approach. IEEE Trans Ind Electron 61
(12):65916602
40. Ramsden VS, Dunlop JB, Holliday WM (1992) Design of a hand-held motor using a rare earth
permanent magnet rotor and glassy metal stator. In: Proceedings of international conference on
electrical machines, pp 376380
41. Ramsden VS, Watterson PA, Dunlop JB (1992) Optimization of REPM motors for specic
applications (invited). In: Proceedings of the international workshop on rare earth materials
and applications, Canberra, pp 296307
42. Ramsden VS et al (1992) Optimization of a submersible brushless rare-earth
permanent-magnet motor for solar power. In: Proceedings of internationl workshop on
electric and magnetic elds, Liege, pp 461464
43. Lovatt HC, Ramsden VS, Mecrow BC (1997) Design of an in-wheel motor for a
solar-powered electric vehicle. In: Proceedings of international conference on electrical
machines and drives, pp 234238
44. Ramsden VS, Zhu JG et al (2001) High performance electric machines for renewable energy
generation and efcient drives. J Renew Energy 22(13):159167
45. Widdowson GP, Howe D, Evison PR (1991) Computer-aided optimization of rare-earth
permanent magnet actuators. In: Proceedings of international conference on computation in
electromagnetics, 2527 Nov 1991, pp 9396
46. Chhaya SM, Bose BK (1993) Expert system based automated simulation and design
optimization of a voltage-fed inverter for induction motor drive. Int Conf IECI 2:10651070,
1519 Nov 1993
47. Guo YG, Zhu JG, Liu DK, Wang SH (2007) Application of multi-level mult-domain
modelling in the design and analysis of a PM transverse flux motor with SMC core. In:
Proceedings of the 7th international conference on power electronics and drive systems
(PEDS07), Bangkok, Thailand, pp 2731
48. Lei G, Liu CC, Zhu JG, Guo YG (in press) Robust multidisciplinary design optimization of
PM machines with soft magnetic composite cores for batch production. IEEE Trans Magn (in
press)
49. Guo YG, Zhu JG, Dorrell D (2009) Design and analysis of a claw pole PM motor with molded
SMC core. IEEE Trans Magn 45(10):5824585
50. Zhu JG, Guo YG, Lin ZW, Li YJ, Huang YK (2011) Development of PM transverse flux
motors with soft magnetic composite cores. IEEE Trans Magn 47(10):43764383
51. Guo YG, Zhu JG, Watterson PA, Wei Wu (2006) Development of a PM transverse flux motor
with soft magnetic composite core. IEEE Trans Energy Conver 21(2):426434
52. Huang YK, Zhu JG et al (2009) Thermal analysis of high-speed SMC motor based on thermal
network and 3D FEA with rotational core loss included. IEEE Trans Magn 45(106):4680
4683
53. Liu CC, Zhu JG, Wang YH, Guo YG, Lei G, Liu XJ (2015) Development of a low-cost double
rotor axial flux motor with soft magnetic composite and ferrite permanent magnet materials.
J. Appl Phys 117(17), Article# 17B507
References 23
54. Zhu ZQ, Howe D (2007) Electrical machines and drives for electric, hybrid, and fuel cell
vehicles. Proc IEEE 95(4):746765
55. Emadi A, Lee YJ, Rajashekara K (2008) Power electronics and motor drives in electric, hybrid
electric, and plug-in hybrid electric vehicles. IEEE Trans Ind Electron 55(6):22372245
56. Kreuawan S, Gillon F, Brochet P (2008) Optimal design of permanent magnet motor using
multidisciplinary design optimization. In: Proceeings of 18th international conference on
electrical machines, Vilamoura, 69 Sept 2008, pp 16
57. Vese I, Marignetti F, Radulescu MM (2010) Multiphysics approach to numerical modeling of
a permanent-magnet tubular linear motor. IEEE Trans Ind Electron 57(1):320326
58. Buja GS, Kazmierkowski MP (2004) Direct torque control of PWM inverter-fed AC motors
a survey. IEEE Trans Ind Electron 51(4):744757
59. Kouro S, Cortes P, Vargas R, Ammann U, Rodriguez J (2009) Model predictive controla
simple and powerful method to control power converters. IEEE Trans Ind Electron 56
(6):18261838
60. Bolognani S, Bolognani S, Peretti L, Zigliotto M (2009) Design and implementation of model
predictive control for electrical motor drives. IEEE Trans Ind Electron 56(6):19251936
61. Xia CL, Wang YF, Shi TN (2013) Implementation of nite-state model predictive control for
commutation torque ripple minimization of permanent-magnet brushless DC motor. IEEE
Trans Ind Electron 60(3):896905
62. Morel F, Lin-Shi XF, Retif J-M, Allard B, Buttay C (2009) A comparative study of predictive
current control schemes for a permanent-magnet synchronous machine drive. IEEE Trans Ind
Electron 56(7):27152728
63. Zhang YC, Zhu JG, Xu W, Guo YG (2011) A simple method to reduce torque ripple in direct
torque-controlled permanent-magnet synchronous motor by using vectors with variable
amplitude and angle. IEEE Trans Ind Electron 58(7):28482859
64. Wang TS, Zhu JG, Zhang YC (2011) Model predictive torque control for PMSM with duty
ratio optimization. In: Proceedings of ICEMS, 2023 Aug 2011, pp 1-5
65. Liu C-H, Hsu Y-Y (2010) Design of a self-tuning PI controller for a STATCOM using particle
swarm optimization. IEEE Trans Ind Electron 57(2):702715
66. Khan MA, Husain I, Islam MR, Klass JT (2014) Design of experiments to address
manufacturing tolerances and process variations influencing cogging torque and back EMF in
the mass production of the permanent-magnet synchronous motors. IEEE Trans Ind Appl 50
(1):346355
67. Coenen I, Giet M, Hameyer K (2011) Manufacturing tolerances: Estimation and prediction of
cogging torque influenced by magnetization faults. In: Proceedings of 14th European
conference on power electronics and applications, pp 19
68. Gasparin L, Fiser R (2011) Impact of manufacturing imperfections on cogging torque level in
PMSM. In: Proceedings of 2011 IEEE ninth international conference on power electronics and
drive systems (PEDS), pp 10551060
69. Lei G, Wang TS, Zhu JG, Guo YG, Wang SH (2015) System level design optimization
method for electrical drive system: robust approach. IEEE Trans Ind Electron 62(8):4702
4713
70. Lei G, Guo YG, Zhu JG et al (2012) System level six sigma robust optimization of a drive
system with PM transverse flux machine. IEEE Trans Magn 48(2):923926
71. Lei G, Zhu JG, Guo YG, Hu JF, Xu W, Shao KR (2013) Robust design optimization of
PM-SMC motors for Six Sigma quality manufacturing. IEEE Trans Magn 49(7):39533956
72. Taguchi G, Chowdhury S, Wu Y (2004) Taguchis quality engineering handbook, Wiley
73. Omekanda AM (2006) Robust torque and torque-per-inertia optimization of a switched
reluctance motor using the Taguchi methods. IEEE Trans Ind Appl 42(2):473478
74. Kim S-I, Lee J-Y, Kim Y-K et al (2005) Optimization for reduction of torque ripple in interior
permanent magnet motor by using the Taguchi method. IEEE Trans Magn 41(5):17961799
75. Ashabani M, Mohamed YA, Milimonfared J (2010) Optimum design of tubular
permanent-magnet motors for thrust characteristics improvement by combined
Taguchi-neural network approach. IEEE Trans Magn 46(12):40924100
24 1 Introduction
76. Hwang C-C, Lyu L-Y, Liu C-T, Li P-L (2008) Optimal design of an SPM motor using genetic
algorithms and Taguchi method. IEEE Trans Magn 44(11):43254328
77. Lee S, Kim K, Cho S et al (2014) Optimal design of interior permanent magnet synchronous
motor considering the manufacturing tolerances using Taguchi robust design. IET Electr
Power Appl 8(1):2328
78. Koch PN, Yang RJ, Gu L (2004) Design for six sigma through robust optimization. Struct
Multidiscipl Optim 26(34):235248
79. Meng XJ, Wang SH, Qiu J, Zhang QH, Zhu JG, Guo YG, Liu DK (2011) Robust multilevel
optimization of PMSM using design for six sigma. IEEE Trans Magn 47(10):32483251
Chapter 2
Design Fundamentals of Electrical
Machines and Drive Systems
Keywords Electrical drive systems Electromagnetic design Thermal design
Mechanical design Power electronics design Control algorithms Finite element
model Model predictive control
2.1 Introduction
Fig. 2.1 Multi-disciplinary design framework of electrical machines and drive systems
Secondly, select a type of the motor, power electronic converter, and control
algorithm from possible options. The motor options include permanent magnet
(PM) motors, induction machines, synchronous machines, DC machines, and swit-
ched reluctance machines. For servo drives, stepping motors and other types of servo
motors can be considered. In this step, different motor topologies have to be inves-
tigated as well. The power electronic converter options mainly include the different
topologies of AC/DC, DC/DC, and DC/AC converters. The controller design mainly
investigates the control strategies and algorithms, such as eld oriented control
(FOC), direct torque control (DTC), and model predictive control (MPC).
Thirdly, based on the selected motor type, converter circuit, and control scheme,
various disciplinary-level analyses should be conducted to evaluate the performance
of the drive system. For example, the motor design analysis consists of mainly the
electromagnetic, thermal and mechanical analyses (the shaded boxes in the gure).
Coupled-eld analyses may be required in the design process, such as
electromagnetic-thermal and electromagnetic-mechanical stress analyses.
In summary, the design of electrical machines and drive systems mainly consists
of the analyses of ve coupled disciplines or domains: electromagnetic, thermal,
mechanical, power electronics, and controller designs. The following sections will
present the popular design analysis models and methods for each discipline.
Power losses and efciency are two main issues in the design analysis of electrical
machines and drive systems. The power losses are mainly composed of the copper
loss, core loss, mechanical loss, and stray loss.
2.1 Introduction 27
(1) The copper loss or Ohmic loss: PCu I 2 R is the power dissipated in stator and
rotor windings due to the resistance of copper wire, where I is the winding
current and R the winding resistance. Normally the DC resistance is used in
the calculation. However, it should be noted that the winding resistance
depends on the operating conditions, i.e., temperature and frequency (due to
the skin effects). In case of the brushes and slip rings/commutator, the effect of
contact resistance is often accounted for by assuming a voltage drop of 2 V.
(2) The core loss is the power dissipated in a magnetic core due to the variation of
magnetic eld. This occurs in the stator and/or rotor iron core of an electrical
machine subject to AC excitations. Practically, it can be measured by
open-circuit or no-load tests. When the magnetic material is under an alter-
nating sinusoidal flux excitation, the alternating core loss can be calculated by
where
" #
Phr 1=s 1=2 s
a1
f a2 1=s2 a23 a2 1=2 s2 a23
and
q
B
s1 1 1=a22 a23
Bs
Bs is the saturation flux density, and Cer, Car, a1, a2 and a3 are the rotational
core loss coefcients.
When the material is under a two dimensional elliptically rotating B excitation,
the core loss can be computed by
Per RB Pr 1 RB 2 Pa 2:3
28 2 Design Fundamentals of Electrical Machines and Drive Systems
where RB Bmin =Bmaj is the axis ratio, Bmin and Bmaj are the magnitudes of
the minor and major axes of the ellipse, respectively, and Pr and Pa the
corresponding rotational and alternating core losses when B = Bmaj. More
details about the rotational core losses can be found in [39].
(3) The mechanical losses are the power losses caused by the friction (brushes,
slip rings/commutator, shaft and bearing), damping, windage, and cooling fan.
It can be approximately determined by no-load test. In design, empirical data
are used.
(4) The stray loss is the power loss caused by stray factors that are hard to
determine separately, such as the non-uniform current distribution in con-
ductors and additional core loss due to distorted magnetic flux distribution for
various reasons. Because it is usually difcult to determine accurately the stray
loss, estimations based on experimental tests and empirical judgment are
2.1 Introduction 29
ZT ZT
m m 2p 2p m
Pin etitdt Em sin t Im sin t dt Em Im 2:5
T T T T 2
0 0
where m is the number of phases, Em the peak value of back electromotive force
(EMF), Im the peak value of phase current, and T the electrical time period. The
output torque can be calculated by
Pout m
Tout g pkp Ksf As Jm 2:6
xr 2
where is the efciency, p the number of pole pairs, p the peak value of PM flux
linkage, r the mechanical rotary speed, Ksf the slot ll factor, As the slot area, and
Jm the peak of current density. For different kinds of PM motors, p and As are
related differently to their dimensions [2123].
(a)
(b)
Fig. 2.3 Photo of the PM-SMC TFM prototype, a PM rotor, and b 3 stack SMC stator
magnetic reluctance of PM, Rg the magnetic reluctance of the air gap, Rst1 the
magnetic reluctance of the stator teeth, Rst2 and Rsy stand for the magnetic reluctance
of the stator yoke. By analyzing this model, the main magnetic flux can be calculated.
Meanwhile, the magnetic flux leakage is a serious problem in this TFM, thus it
should be considered in the magnetic circuit model. Several flux leakage models
can be constructed for this TFM. Figure 2.5 illustrates the main flux leakage model.
In this model, the adjacent PM in the one side of the machine is modeled, where
Rry1 represents the magnetic reluctance of rotor, Rg1 and Rg2 represent the magnetic
reluctance of the air gap, Rs1 stands for the magnetic reluctance of the stator.
With the computed flux linkage, the resultant magnetic flux density in the air gap
and the flux per turn of coil can be estimated. After calculation, the obtained flux
per turn of this PM-SMC TFM is 0.32 mWb, which is higher than the calculated
result (0.28 mWb) by using the FEM [24].
This model can be also used to evaluate the performance of the motor. Based on
the calculated magnetic flux of the motor, the flux linkage per phase equals the
32 2 Design Fundamentals of Electrical Machines and Drive Systems
(a)
(b)
Fig. 2.4 Main flux circuit and flux path of the PM-SMC TFM, a magnetic circuit model, b flux
path in 2D plane
(a)
(b)
Fig. 2.5 Flux leakage circuit and path of the PM-SMC TFM, a magnetic circuit model, b leakage
path in 2D plane
number of coil turns multiplied by the magnetic flux of each coil turn, and it can be
computed as
where PM is the PM flux linkage per phase, kl the leakage coefcient, Ncoil the
number of turns of the phase winding, p the number of pole pairs, and gap is flux
per coil turn. The back EMF can be expressed as
2.2 Electromagnetic Design 33
After the calculation, the no-load back EMF is 53.26 V at the rated speed of
1800 rev/min. According to (2.9), the electromagnetic torque is 4.66 Nm at the
rated current of 5.5 A (RMS value). Compared to the electromagnetic torque
obtained from FEM, i.e., 4.08 Nm, the relative error is about 0.58/4.08 = 14.21 %.
FEM is a widely used analysis model for eld analysis in electrical machines as
well as other electromagnetic devices. The theory of FEM can be found in many
books and research papers. The PM-SMC TFM investigated above will be
employed as an example to show the application of FEM for designing electrical
machines.
When analyzing the magnetic eld distribution, we used eld analysis software
package ANSYS, and taking advantage of the periodical symmetry, we only need
to analyze one pole-pair region of the machine, as shown in Fig. 2.6a. At the two
radial boundary planes, the magnetic scalar potential obeys the periodical boundary
conditions:
Fig. 2.6 a One pole pitch of FEM solution region for one phase (stack), and b magnetic eld
distribution under no-load
34 2 Design Fundamentals of Electrical Machines and Drive Systems
where Dh 18 mechanical is the angle of one pole pitch. The origin of the
cylindrical coordinate is located at the center of the stack.
Figure 2.6b illustrates the magnetic eld distribution under no-load. Based on
the FEM analysis, the calculated key motor parameters for this machine are listed in
Table 2.2. The measured parameters are also listed in the table to show the effec-
tiveness of the FEM method. As shown, the measured motor back EMF constant is
0.244 Vs, 1 % lower than the calculated value of 0.247 Vs. The calculated phase
resistance and inductance, and maximal cogging torque are 0.310 , 6.68 mH and
0.339 Nm, respectively, which are very close to the measured values (0.305 ,
6.53 mH and 0.320 Nm). In summary, the estimated parameters calculated by the
FEM-based method are well aligned with the experimental results. Therefore, FEM
is better than magnetic circuit model, and it is reliable to be used for optimization of
the electromagnetic design of electrical machines.
Moreover, the output performance parameters, such as output power, torque and
efciency, can be estimated with the calculated electromagnetic parameters men-
tioned above. In the estimation, the control method is assumed to maintain that the
d-axis component of current equals zero. Figure 2.7 shows the per phase equivalent
electric circuit of this motor under the assumed control method.
Based on this per phase equivalent electrical circuit, the main relationships of the
motor can be predicted by
q
Vin Ea Ia Ra 2 xe La Ia 2 2:11
Pout
Tout 2:14
xr
where Vin is the input voltage, Ea the back EMF, Ia the armature current, e the
electric angular frequency, La the inductance, Ra the resistance, the angle between
Vin and Ea, Pin the input power, Pout the output power, Pcore the core loss, Pcopper
the copper loss, Pmech the mechanical loss, Tout the output torque, and r the
mechanical angular speed.
In motor with SMC cores, unlike the conventional motors made of silicon sheet
steels, the core loss can be a major part among all power losses, and the mechanical
2.2 Electromagnetic Design 35
loss is generally considered as 11.5 % of the output power. In general, the core
loss prediction in the TFM should be calculated by using the FEM based on the
multi-frequency core loss characteristic of the material. More comparison results
can be seen in [24, 25].
The rating of an electrical machine gives its working capability under the specied
electrical and environmental conditions. Major factors that determine the ratings are
thermal and mechanical considerations. To obtain an economic utilization of the
materials and safe operation of the motor, it is necessary to predict with reasonable
accuracy the temperature rise of the internal parts, especially in the coils and
magnets.
The temperature rise resulted from the power losses in an electrical machine
plays a key role in rating the power capacity of the machine, i.e., the amount of
power it can convert without being burnt for a specied length of life time. The life
expectancy of a large industrial electrical machine ranges from 10 to 50 years or
more. In an aircraft or electronic equipment, it can be of the order of a few thousand
hours, whereas in a military application, e.g. missile, it can be only a few minutes.
36 2 Design Fundamentals of Electrical Machines and Drive Systems
where A and B are constants and T is the absolute temperature. Roughly, it says that
for each 10 C temperature rise exceeding the maximum allowable temperature rise,
the life time of insulation is halved.
Insulation materials used in electrical machines are classied by the maximum
allowable temperature rise that can be safely withstood. Table 2.3 lists the classi-
cation of electrical insulation materials by the IEC (International Electrotechnical
Commission).
Generally, there are two kinds of analysis models for thermal analysis in elec-
trical machines, namely the thermal network model and the FEM [14, 15, 20, 25].
The following sections will present examples for the two methods.
Two design examples will be illustrated to show the usage of thermal network
model for the thermal analysis of PM-SMC motors. The rst one is a TFM, and the
second is a high speed claw pole motor.
A. Transverse flux machine
In this study, the temperature rise was calculated by using a hybrid thermal network
model with distributed heat sources, as shown in Fig. 2.8.
For high computation accuracy, every part, e.g. the air gap, is divided into two or
more segments. The thermal resistances to heat conduction in the following sections
are calculated: rotor yoke (Rry), magnets (Rm), glue between magnets and rotor yoke
2.3 Thermal Design 37
(Rmg), air gap (Rag), stator yoke (RFe1), stator side discs (RFe2), stator teeth (RFe3),
varnished copper wire (Rcu), and insulations (RI1, RI2, RI3) between the winding and
the stator yoke, the stator wall disc, and the air gap, respectively. In addition, the
thermal resistances of the stator shaft (Rss), the aluminum end plates (Ral), and the
stationary air (Rsa) between the side stator discs and the end plates are calculated
separately [25].
The equivalent thermal resistances to the heat convection of the following
sections are calculated: that between the stator tooth surface and the inner air in the
air gap (RFeA), that between the winding and the inner air (RWA), that between the
magnet and the inner air (RmA), that between the rotor yoke and the inner air
(RryA1), and that between the rotor yoke and the outer air (RryA2).
The heat sources include the stator winding copper losses (Pcu), the stator and
rotor core losses (PFes, PFer), and the mechanical losses due to windage and friction
(Pmec). The improved method for core loss calculation can obtain the loss distri-
bution, which is a great advantage for thermal calculation by the hybrid thermal
model.
The temperature rises in the middle of several parts are calculated as 64.9 C in
the stator winding, 78.6 C in the stator core, 59.3 C in the air gap, 36.1 C in the
magnets, and 25.3 C in the rotor yoke outer surface. The experimentally measured
results are 66 C in the stator winding and 27 C in the rotor yoke, and it can be
seen that the maximum relative error between the calculated and measured results is
only 3 %. Thus, it is reliable to use the thermal network method for design of this
TFM.
38 2 Design Fundamentals of Electrical Machines and Drive Systems
such as the 2D and 3D heat flow, convection, internal heat generation, and varia-
tions in material properties. To account for the three dimensional heat flows at a
node, the thermal structure shown in Fig. 2.13 can be employed.
As shown in Fig. 2.13, the thermal resistances of an element are built in three
directions, and the heat source if any can be placed at the center point. In this
model, the thermal conduction equation can be expressed as
Tb Ta Tc Ta Td Ta Te Ta Tf Ta Tg Ta
qa
Rab Rac Rad Rae Raf Rag
0
@Ta Ta
Ca 2:16
@t
where Ta, Tb, Tc, Td, Te, Tf, and Tg are the temperatures at nodes a, b, c, d, e, f, and
g, Rab, Rac, Rad, Rae, Raf, and Rag the thermal resistance between nodes a-b, a-c, a-d,
0
a-e, a-f, and a-g, respectively, qa is the heat source, Ca the heat specic, and Ta the
2.3 Thermal Design 41
temperature of node a at the next time instant. The thermal resistance in Fig. 2.13
can be calculated by
DX
Rab Rac 2:17
2kx DYDZ
DY
Rad Rae 2:18
2ky DXDZ
DZ
Raf Rae 2:19
2kz DXDY
In the thermal network, the core loss at each node cannot be obtained easily from
the magnetic eld calculation. In most cases, the average value is used. Since the
core loss distribution is quite different in different positions of the stator core, 3D
FEM is used to analyze the temperature distribution in this section. Two design
examples investigated in the previous section will be illustrated to show the usage
of FEM for the thermal analysis of PM-SMC motors.
42 2 Design Fundamentals of Electrical Machines and Drive Systems
Figure 2.14 illustrates the temperature distribution of the PM-SMC TFM based
on FEM. As shown, the average temperature rises in the winding is 62.5 C, which
is close to the measured value 65 C.
Figure 2.15 depicts the distributions of core loss and temperature at full load in
the SMC core of the high speed claw pole motor. The temperature is measured by
an infrared temperature probe. At 20,000 rev/min and no load, the frame temper-
ature is 331.4 K and the stator yoke temperature is 333.5 K, respectively. The
Fig. 2.15 a Distributions of core loss, and b temperature in SMC core of the high speed claw pole
motor
2.3 Thermal Design 43
measured temperatures are slightly higher than the FEM results, because the actual
loss is greater than the calculation. The FEM method is more accurate than the
thermal network method because there are only ten nodes in the network. The
advantage of the thermal network is the calculation speed, which is much faster than
the FEM method [8, 9].
Fig. 2.16 Illustration of rst order modal analysis for PM-SMC TFM
Fig. 2.17 Vibration patterns at a 4,102 Hz (Y axis), b 4,102 Hz (X axis), c 9,562 Hz (Z axis), and
d 10,321 Hz (Y axis)
2.5 Power Electronics Design 45
The design of power electronics for electrical machines and drive systems is also an
important and complex stage. Among many aspects in power electronics, the
converter/inverter and switching scheme are two main concerns in the design of
electrical machines and drive systems.
The converter/inverter is an important component to drive an electrical machine.
An inverter, for example, is an electronic apparatus that can convert a DC voltage to
an AC voltage of specied waveform, frequency, magnitude, and phase angle.
Among many different topologies, the three phase bridge power circuit as shown in
Fig. 2.18 has become favorite and standard for use in the control systems of
electrical machines. Many different topologies can be obtained from this structure
for different applications. For example, two extra switches can be added to establish
two bridges for the fault tolerant control scheme [33, 34].
For controlling the waveform, frequency, magnitude, and phase angle of the AC
voltage, many switching schemes can be used, such as square wave and sine wave
pulse width modulations (PWMs) and space vector modulation (SVM), as well as
hard and soft switching.
The six-step control method was oriented to drive brushless DC (BLDC) motors
with trapezoidal back EMF waveforms. In many applications, however, the
trapezoidal excitation is also used to drive PMSMs with sinusoidal back EMF
waveforms because the trapezoidal excitation or six-step method based drive is
robust and low cost [35].
In the six-step control scheme, the stationary reference frame is always used to
model the PMSM. The phase variables are used to express the machine equations as
they can account for the real waveforms of the back EMF and phase current.
Assuming that the resistances of three phase stator windings are equal, the three
phase voltage equations of the motor can be written as
2.6 Control Algorithms Design 47
2 3 2 32 3 82 32 39 2 3
d <4 aa
va Rs 0 0 ia L Lba Lca ia = ea
4 vb 5 4 0 Rs 0 54 ib 5 Lba Lbb Lcb 54 ib 5 4 eb 5
dt : ;
vc 0 0 Rs ic Lca Lcb Lcc ic ec
2:20
where va, vb, and vc are the phase voltages, ia, ib, and ic the 2phase currents, e3a, eb,
Laa Lba Lca
and ec the phase back EMF, Rs is the phase resistance, and 4 Lba Lbb Lcb 5 the
Lca Lcb Lcc
inductance matrix, including both the self-and mutual-inductances.
Assuming further that the reluctance is independent of the rotor position, one can
obtain
La Lb Lc Ls
2:21
Lab Lca Lbc M
Assuming linear system, the machine model in state space form can be expressed
as
2 3 2 382 3 2 32 3 2 39
i 1=L M 0 0 < va Rs 0 0 ia ea =
d 4 a5 4 5 4 vb 5 4 0
ib 0 1=L M 0 Rs 0 54 i b 5 4 eb 5
dt : ;
ic 0 0 1=L M vc 0 0 Rs ic ec
2:23
Te ea ia eb ib ec ic =xm 2:24
dxm
Te J Fxm TL 2:25
dt
where J is the inertia of the machine rotating parts, F the friction coefcient, and TL
the load torque on the rotor shaft.
48 2 Design Fundamentals of Electrical Machines and Drive Systems
DC voltage source or
rectified from AC power
iDC
-
ref + Speed + Current PWM & Inverter
Controller Controller Commutation
r -
Speed Hall
Sensor PMSM
Calculation
Figure 2.19 shows the block diagram of six-step drive scheme. The drive system
is operated with the feedback information of rotor position, which is obtained at
xed points, typically every 60 electrical degrees for commutation of the phase
currents.
The 120 conduction mode is applied to drive the PMSM. The voltage may be
applied to the motor every 120 (electrical), with a current limit to hold the phase
currents within the motors capabilities. Because the phase currents are excited in
synchronism with the back EMF, a constant torque is generated. A simulation
model is built in MATLAB/SIMULINK as shown in Fig. 2.20.
Fig. 2.20 Simulation block diagram of six-step controlled PMSM drive system
2.6 Control Algorithms Design 49
As shown in Fig. 2.20, the rotor position information comes from the Hall effect
sensors, which are integrated in the machine model in MATLAB/SIMULINK. The
resolution of the feedback signals is only 60 (electrical). Since most applications
require a stable speed, a speed feedback loop is employed. The rotor speed infor-
mation can be deduced from the low resolution Hall signals, which is marked as
Speed Calculation in Fig. 2.20. Typically, the average speed in one 60 section is
used as the speed feedback.
However, by using the average speed, there is always a lag when the motor
speed is not constant in accelerating or other dynamic state. To overcome this, the
rotor position can be expressed in Taylors series as the following:
2
1 h2k
ht hk t h1k t tk t t k 2 2:26
2!
1
p=3
where tk is the last commutation time, h1k tk tk1
the average speed of last section,
1 1
2 h h
and h2k 1ktk tk1
1k1
the average acceleration of last section.
As shown above, with the higher order calculation, more accurate speed and
position information can be deduced, whereas the computing cost rises. As a
compromise, in some situations, the following equations are used to estimate the
rotor position and speed:
( 2
1 h2k
ht hk t h1k t tk 2! t tk 2 2:27
1 2
xt h1k t h2k t tk
For a PMSM under sinusoidal excitations, the original voltage equations can be
expressed in the stationary reference frame as the following
2 3 2 3 2 3
va ia ka
4 vb 5 Rs 4 ib 5 d 4 kb 5 2:28
dt
vc ic kc
2 3 r2 32 3
rd cos h cos h 23 p cos h 23 p ra
4 rq 5 24 sin h sin h 23 p sin h 23 p 54 rb 5 2:29
3 p1 p1 p1
r0 2 2 2
rc
where vd, vq, and v0 are the phase voltages, id, iq, and i0 the phase currents, and d,
q, and 0 the phase flux linkages.
For the linear PMSM model, the magnetic saturation saliency is not considered.
The flux linkages of the d- and q-axes can be further expressed as
kd Ld id km
2:31
kq Lq iq
where Ld and Lq are the constant d- and q-axes inductances, respectively, and m is
the flux linkage generated by the rotor PMs.
On the other hand, the voltage equation of the 0 axis in (2.30) is usually ignored
by assuming well-balanced three-phase windings for the controller design.
2.6 Control Algorithms Design 51
Therefore, the electrical voltage equations in the rotor reference frame can be
rewritten as
(
vd Rs id Ld didtd Lq iq dh
di
dt
2:32
vq Rs iq Lq dtq Ld id km dh
dt
3
Te p kd iq kq id
2 2:33
3
p km iq Ld Lq id iq
2
DC voltage source or
rectified from AC power
ref + + vd va
PI PI va
r - - 1/ Park
v
1/ Clark vb
PWM VSI
+ vq Transform Transform vc
PI
id_ref -
id ia ia
iq Park i Clarke
Transform ib
Transform
Position and Encoder
Speed PMSM
Estimation
Fig. 2.23 Simulation block diagram of typical FOC based PMSM drive system
Similar to the six-step method, a simulation model of the FOC scheme based
PMSM drive is built in MATLAB/SIMULINK. The sinusoidal back EMF machine
model is selected from the SimPowerSystem tool box, in which the current sensors
and rotor position sensor are integrated. The Park and Clark transformations are
synthesized as one abc_to_dq block to transfer the variables between the sta-
tionary and rotating reference frames, as shown in Fig. 2.23. Two discrete PI
controllers are used for the speed and current feedback loops.
The traditional triangulation PWM generation technique is applied. A triangular
carrier wave sampling signal is compared directly with a sinusoidal modulating
wave to determine the switching instants, and therefore the resultant pulse widths.
In the DTC strategy, the flux linkage and torque are calculated in the two-phase
stator reference frame, i.e., the - frame, which is transformed from the
three-phase a-b-c reference frame by using the Clark transformation. The Clark
transformation can be expressed in the matrix form as
r" #2 3
1
p12 ra
ra 2 1 p2 4 rb 5
2:34
rb 3 0 3
23
2 rc
After the measured phase voltages and currents are transformed to the - frame,
the flux linkage components of the - and -axes can be calculated as
R
ka R va Rs ia dt
2:35
kb vb Rs ib dt
2.6 Control Algorithms Design 53
DC voltage source or
rectified from AC power
ref + e Hysteresis
- controller
Switch
Inverter
table
ref + + eT Hysteresis
Controller
- - controller
r s
vdc
Flux & Torque ia
calculation ib
Encoder
1/s PMSM
Fig. 2.24 Block diagram of typical DTC scheme based PMSM drive
3 pm
Te ka ib kb ia 2:36
2 2
Figure 2.24 shows the block diagram of a typical DTC scheme for PMSM drive.
Two hysteresis controllers are applied to the flux linkage and torque control loops.
The calculated flux linkage is also sent to the switching table to identify the current
flux vector position.
From (2.35), the stator flux linkage is
Z
ks vs Rs is dt 2:37
where vs and is are the stator voltage and current spatial vectors, respectively.
In the case of a PMSM, s always varies even when the zero voltage vectors are
applied because of the rotating rotor magnets, and thus, zero voltage vectors are not
used for DTC driven PMSM. s should always be in motion with respect to the rotor
flux.
According to (2.36), the electromagnetic torque can be controlled effectively by
controlling the amplitude and rotating speed of s. For counter-clockwise operation,
if the actual torque is smaller than the reference, the voltage vectors that keep s
rotating in the same direction are selected. The angle increases as fast as it can, and
the actual torque increases as well. Once the actual torque is greater than the
reference, the voltage vectors that keep s rotating in the reverse direction are
selected instead of the zero voltage vectors. The angle decreases, so does the torque.
By selecting the voltage vectors in this way, s will rotate all the time in the
54 2 Design Fundamentals of Electrical Machines and Drive Systems
Table 2.5 Switching table of typical DTC scheme for PMSM drive
e eT
1 2 3 4 5 6
1 1 V2(110) V3(010) V4(011) V5(001) V6(101) V1(100)
0 V6(101) V1(100) V2(110) V3(010) V4(011) V5(001)
0 1 V3(010) V4(011) V5(001) V6(101) V1(100) V2(110)
0 V5(001) V6(101) V1(100) V2(110) V3(010) V4(011)
direction determined by the output of the hysteresis controller for the torque. The
switching table for controlling both the amplitude and rotating direction is shown in
Table 2.5, in which the inverter voltage vector and spatial sector denitions are
illustrated in Fig. 2.25.
Figure 2.26 shows the simulation model built based on the typical DTC scheme.
The inverter switching status and DC bus voltage are utilized to calculate the stator
voltage. The stator flux linkage is obtained in the observer. The traditional two-level
hysteresis controllers are applied and the switching table is designed based on
Table 2.5.
The principle of MPC was introduced for industrial control applications in the
1970s after the publication of this strategy in the 1960s. The MPC requires great
computational effort and it has been formerly limited to slowly varying systems,
such as chemical processes. With the availability of inexpensive high computing
2.6 Control Algorithms Design 55
Fig. 2.26 Simulation block diagram of typical DTC based PMSM drive system
where Te and ws are the reference torque and flux, Tek 1 and wsk 1 the predicted
values of torque and flux, respectively, and k1 is the weighting factor. Because the
physical natures of electromagnetic torque and stator flux are different, the
weighting factor k1 is introduced to unify these terms. In this work, k1 is selected to
be Tn =wn , where Tn and wn are the rated values of torque and stator flux, respec-
tively. It should be noted that when a null vector is selected, the specic state (V0 or
V7) will be determined based on the principle of minimal switching commutations,
which is related to the switching states of the previous voltage vector.
The voltage equations in the d-q reference frame are as follows:
56 2 Design Fundamentals of Electrical Machines and Drive Systems
did
ud Rs id Ld xLq iq 2:39
dt
diq
u q R s i q Lq xLd id xwf 2:40
dt
Given the voltage and current values at sampling instant k, the predicted current,
torque and flux at instant k + 1 can be expressed as follows:
1
ikd 1 ikd Rs ikd xk Lq ikq ukd Ts 2:41
Ld
1
k k
ikq 1 ikq x Ld id Rs ikq ukq xk wf Ts 2:42
Lq
wks 1 Ld ikd 1 wf jLq ikq 1 2:43
3
Tek 1 pwks 1 iks 1 2:44
2
where ikd 1 and ikq 1 are the predicted values of stator current for the sampling
instant k + 1, Ts is the sampling period, Tek 1 and wks 1 are the predicted values of
torque and flux, respectively, which are also the main concerns for the cost function
in the following MPC control scheme [1, 36, 49].
The block diagram of MPC is shown in Fig. 2.27. The inputs of the system are
the reference and estimated values of torque and flux. By evaluating the effects of
each voltage vector when applied to the machine, the voltage vector which mini-
mizes the difference between the reference and predicted values is rst selected, and
then it is generated by the inverter.
The cost function in (2.38) assumes that all calculations and judgments are
implemented at the kth instant and the selected vector will be applied immediately.
However, in practical digital implementation, this assumption is not true and the
applied voltage vector is not applied until the (k + 1)th instant.
In other words, for the duration between the kth and (k + 1)th instants, the
applied rotor voltage vector uks has been decided by the value in the (k-1)th instant
and the evolutions of ws and Te for this duration are uncontrollable. What is left to
be decided is actually the stator voltage vector uks 1 , which is applied at the
beginning of the (k + 1)th instant. To eliminate this one step delay, the variables of
wks 2 and Tek 2 should be used rather than wks 1 and Tek 1 for the evaluation of the
cost function in (2.38). This fact is clearly illustrated in Fig. 2.28, where x indicates
the state variables of a dynamic system and u is the input to be decided. For PMSM,
x represents torque or stator flux value.
To eliminate the one-step delay in digital implementation, the cost function in
(2.38) should be changed to (2.45) as shown below
min: G jTe Tek 2 j k1 jws j jwks 2 j
2:45
s.t. uks 2 fV0 ; V1 ; . . .; V7 g
Obtaining wks 2 and Tek 2 in (2.45) requires a two-step prediction. To obtain the
best voltage vector minimizing the cost function in (2.45), each possible congu-
ration for uke 1 will be evaluated to obtain the value at the (k + 2)th instant.
In this section, the simulation tests of DTC and MPC are carried out by using
Matlab/Simulink. The parameters of the motor are listed in Table 2.6. The sampling
frequency of both methods is set to 5 kHz. The values of control parameters are
k1 25:4; A 0:1; and N 10 [36].
This simulation test combines start-up, steady-state and external load tests. The
motor starts up from 0 s with several reference speeds (500 rev/min, 1000 rev/min,
1500 rev/min and 2000 rev/min). After reaching the reference speed, the motor
maintains the speed for at least 0.2 s and an external load is applied at 0.3 s.
Figures 2.29, 2.30, 2.31 and 2.32 show the combined load test for four control
strategies for one reference speed, 1000 rev/min. From top to bottom, the curves are
the stator current, stator flux, torque, motor speed, and switching frequency,
respectively. The test results for other speed situations can be found in [36].
ia/A
0
-10
0.2
Phi/Wb
0.15
0.1
0.05
6
4
Te/Nm
2
0
2000
Speed/rpm
1500
1000
500
0
4000
fs/Hz
2000
0
0 0.1 0.2 0.3 0.4 0.5
t/s
By comparing Fig. 2.30 with Fig. 2.29, it is shown that the torque and flux
ripples of MPC are lower than that of DTC. In Fig. 2.31, MPC with one-step delay
compensation (indicated as MPC + comp) presents torque and flux ripples even
lower than MPC along with an increase in switching frequency. Figure 2.32
illustrates the responses by using cost function (2.48), where factor A is included in
the simulation. As shown, the introduction of linear multiple horizon prediction
(factor A, and indicated as MPC + A) can greatly reduce the switching frequency
only with a quite limited degradation of torque and flux ripples. As shown, all these
methods present similar dynamic performance and the motor can reach the refer-
ence speed rapidly. When the load was applied, the motor speed returned to its
original value in a very short time period.
The recorded data from 0.1 to 0.3 s are picked to calculate the torque and flux
ripples (obtained by standard deviations). The torque and flux ripples of these
control methods are summarized in Table 2.7. A segment (three periods) of the
stator current of phase A is used to calculate the total harmonic distortion
(THD) and current harmonic spectrum.
As shown, MPC can achieve lower torque ripple than that of DTC as proven.
However, MPCs characteristic in flux ripple reduction is quite unstable. With the
help of one-step delay compensation, the steady-state performance of MPC is
60 2 Design Fundamentals of Electrical Machines and Drive Systems
ia/A
0
-10
0.2
Phi/Wb
0.15
0.1
0.05
6
4
Te/Nm 2
0
2000
Speed/rpm
1500
1000
500
0
4000
fs/Hz
2000
0
0 0.1 0.2 0.3 0.4 0.5
t/s
In addition to the simulation study, the control methods mentioned above are
further experimentally tested on a two-level inverter-fed PMSM motor drive. The
experimental setup is illustrated in Fig. 2.33. A dSPACE DS1104 PPC/DSP control
board is employed to implement the real-time algorithm coding using C language.
A three phase intelligent power module equipped with an insulated-gate bipolar
transistor (IGBT) is used as an inverter. The gating pulses are generated in the
DS1104 board and then sent to the inverter. The load is applied using a pro-
grammable dynamo-meter controller DSP6000 (Fig. 2.34). A 2500-pulse incre-
mental encoder is equipped to obtain the rotor speed of PMSM. All experimental
results are recorded by the ControlDesk interfaced with DS1104 and PC at 5 kHz
sampling frequency [36].
2.6 Control Algorithms Design 61
ia/A
0
-10
0.2
Phi/Wb
0.15
0.1
0.05
6
4
Te/Nm 2
0
2000
Speed/rpm
1500
1000
500
0
4000
fs/Hz
2000
0
0 0.1 0.2 0.3 0.4 0.5
t/s
The steady-state responses at 1000 rpm are presented in this section. From top to
bottom, the curves shown are torque, stator flux and switching frequency,
respectively.
Figures 2.35 and 2.36 show the measured steady-state performance at 1000 rpm.
It is seen that the implementation of MPC can reduce the torque ripple, but does not
reduce the flux ripple. When the one-step delay is compensated, a signicant
decrease of torque and flux ripples can be found as well as an obvious increase of
switching frequency. When the linear multiple horizon prediction is added to MPC,
it can be seen that the torque and flux ripples are slightly decreased along with a
limited reduction of the switching frequency.
Table 2.8 lists the torque and flux ripples of these control methods in experiment.
As shown, similar conclusions can be obtained as those from Table 2.7. According
to the analysis above, it can be concluded that:
(1) MPC can achieve lower torque ripple than that of DTC whilst
maintaining/reducing the switching frequency as proven in both simulation
and experimental tests. However, MPCs ability in flux ripple reduction is
insignicant and even unstable.
62 2 Design Fundamentals of Electrical Machines and Drive Systems
ia/A
0
-10
0.2
Phi/Wb
0.15
0.1
0.05
6
4
Te/Nm 2
0
2000
Speed/rpm
1500
1000
500
0
4000
fs/Hz
2000
0
0 0.1 0.2 0.3 0.4 0.5
t/s
Fig. 2.33 Experimental setup of testing system: a overview of the testing platform and b front
view of the PMSM and inverter control board
There are many improvements for these control algorithms. One of them known as
MPC with duty ratio optimization will be selected for the control of PM-SMC
TFM. As a general algorithm, the theory and test results will be presented in this
section.
64 2 Design Fundamentals of Electrical Machines and Drive Systems
1 1
Te/Nm
Te/Nm
0 0
-1 -1
-2 -2
0.16 0.16
s/Wb
s/Wb
0.10 0.10
0.04 0.04
4000 4000
fs/Hz
fs/Hz
2000 2000
0 0
0 1 2 0 1 2
t/s t/s
In the conventional MPC, the selected voltage vector works during the whole
sampling period. In many cases, it is not necessary to work for the entire period to
meet the performance requirement of torque and flux. This is one of the main
reasons for the torque and flux ripples. By introducing a null vector to each sam-
pling period, the effects of voltage on torque can be adjusted to be more moderate,
in order to diminish the ripples of torque and flux.
2.6 Control Algorithms Design 65
Te/Nm
0 0
-1 -1
-2 -2
0.16 0.16
s/Wb
s/Wb
0.10 0.10
0.04 0.04
4000 4000
fs/Hz
fs/Hz
2000 2000
0 0
0 1 2 0 1 2
t/s t/s
Fig. 2.36 Steady-state response for: a MPC with one-step delay compensation, and b MPC with
linear multiple horizon prediction
Actually, the torque can be changed by adjusting the amplitude and time
duration of us . The amplitude is decided by the DC bus voltage and is usually xed,
while the time duration of us can be varied from zero to the whole period, which is
equivalent to changing the voltage vector length. The null vector only decreases the
torque, while appropriate non-zero vectors can increase the torque, and it is possible
to employ both null and nonzero vectors during one cycle to reduce the torque
ripple. The appropriate non-zero vectors are also referred as active vector. The key
issue is how to determine the time duration of the two vectors, or the duty ratio of
the active vector.
The expression of duty ratio for MPC is shown as follows
T Tek 1 ws wek 1
d e
;
2:49
C T C w
66 2 Design Fundamentals of Electrical Machines and Drive Systems
Fig. 2.37 Diagram of an improved MPC with duty ratio optimization in MATLAB/SIMULINK
where d is the duty ratio of the active voltage vector, and CT and C are two
positive parameters. The idea of this method is that the larger the difference between
the reference and predicted torque values, the larger is the duty ratio value [36]. On
the other hand, the lower the CT and C values, the quicker is the dynamic response
(e.g. take less time to reach the given speed), but the poorer will be the steady-state
response (e.g. higher torque and flux ripples). Higher values of CT and C could
lead to better steady-state responses, but slower dynamic responses. Therefore, the
determination of these values is a compromise between the steady-state and
dynamic performances. Extensive simulation and experimental results have proven
that the PM flux value and half-rated torque value for CT and C can provide a good
compromise between the steady state performance and the dynamic response. The
block diagram of the proposed improved MPC is shown in Fig. 2.37.
The parameters of the motor and control system simulated in this section are listed
in Table 2.9. Similar to the previous test example, this simulation test combines the
start-up, steady-state and external load tests. The motor starts up from 0 s with
several reference speeds (500 rev/min, 1000 rev/min, 1500 rev/min and 2000
rev/min). After reaching the reference speed, the motor maintains the speed for at
least 0.2 s and an external load is applied at 0.3 s. Figure 2.38 shows the combined
load test for one reference speed, 1000 rev/min. From top to bottom, the curves are
2.6 Control Algorithms Design 67
-10
0.2
Phi/Wb
0.15
0.1
0.05
6
4
Te/Nm
2
0
2000
Speed/rpm
1500
1000
500
0
1
Duty
0.5
6000
fs/Hz
4000
2000
0
0 0.1 0.2 0.3 0.4 0.5
t/s
the stator current, stator flux, torque, motor speed, and switching frequency,
respectively. The test results for other speed situations can be found in [36].
It can be found that the proposed MPC scheme present very low torque and flux
ripples and excellent dynamic response. The proposed MPC scheme also presents
very low stator current THDs and narrow harmonic spectrums with the dominant
harmonics of around 5 kHz.
68 2 Design Fundamentals of Electrical Machines and Drive Systems
The experimental tests are performed on the same testing platform introduced in the
last section. Figure 2.39 shows the steady-state responses at 1000 rpm for three
control strategies, namely (a) DTC, (b) MPC, and (c) MPC with duty ratio
optimization.
Te/Nm
0 0
-1 -1
-2 -2
0.16 0.16
Phis/Wb
Phis/Wb
0.10 0.10
0.04 0.04
4000 4000
fs/Hz
fs/Hz
2000 2000
0 0
0 1 2 0 1 2
t/s t/s
MPC Duty
(c) 2
1
Te/Nm
0
-1
-2
0.16
Phis/Wb
0.10
0.04
4000
fs/Hz
2000
0
1
duty
0.5
0
0 1 2
t/s
Fig. 2.39 Steady-state response at 1000 rpm for: a DTC, b MPC and c MPC with duty ratio
optimization
2.6 Control Algorithms Design 69
It can be seen that in MPC with duty ratio optimization, the torque and flux
ripples are reduced signicantly compared to other methods. The duty ratio
increases along with the increase in motor speed.
According to the analysis above, it can be concluded that:
(1) MPC with duty ratio optimization can achieve a better performance than DTC
and original MPC in terms of torque and flux ripples reduction;
(2) Under the same system sampling frequency (5 kHz), the switching frequency
of the improved method is much higher than other methods; and
(3) In DTC and MPC, the switching frequency slightly decreases along with the
increase of motor speed. However, the switching frequency is almost stable in
the proposed method.
More experimental results including different speed, dynamic response and data
analysis can be found in [36].
2.7 Summary
This chapter presents the multi-disciplinary design analysis models and methods for
electrical machines and drive systems. All the models and methods are discussed in
terms of the three major parts of electrical drive systems, namely electrical
machines, power electronic converters and controllers. Electromagnetic, thermal
and mechanical analyses based on different models, e.g. FEM, have been investi-
gated for the design of electrical machines with several prototypes developed in our
research center. Various kinds of popular control algorithms have been described
for the controller design. Several examples investigated in our previous work have
been presented to show the effectiveness of the proposed models and analysis
methods.
References
1. Lei G, Wang TS, Guo YG, Zhu JG, Wang SH (2014) System level design optimization
methods for electrical drive systems: deterministic approach. IEEE Trans Ind Electron 61
(12):65916602
2. Lei G, Wang TS, Zhu JG, Guo YG, Wang SH (2015) System level design optimization
method for electrical drive system: robust approach. IEEE Trans Ind Electron 62(8):4702
4713
3. Zhu JG, Guo YG, Lin ZW, Li YJ, Huang YK (2011) Development of PM transverse flux
motors with soft magnetic composite cores. IEEE Trans Magn 47(10):43764383
4. Zhu JG, Ramsden VS (1998) Improved formulations for rotational core losses in rotating
electrical machines. IEEE Trans Magn 34(4):22342242
5. Guo YG, Zhu JG, Lu HY, Lin ZW, Li YJ (2012) Core loss calculation for soft magnetic
composite electrical machines. IEEE Trans Magn 48(11):31123115
70 2 Design Fundamentals of Electrical Machines and Drive Systems
6. Guo YG, Zhu JG, Lu HY, Li YJ, Jin JX (2014) Core loss computation in a permanent magnet
transverse flux motor with rotating fluxes. IEEE Trans Magn 50(11). Article#: 6301004
7. Guo YG, Zhu JG, Zhong JJ, Wu W (2003) Core losses in claw pole permanent magnet
machines with soft magnetic composite stators. IEEE Trans Magn 39(5):31993201
8. Huang YK, Zhu JG, Guo YG, Lin ZW, Hu Q (2007) Design and analysis of a high speed claw
pole motor with soft magnetic composite core. IEEE Trans Magn 43(6):24922494
9. Huang YK, Zhu JG et al (2009) Thermal analysis of high-speed SMC motor based on thermal
network and 3-D FEA with rotational core loss included. IEEE Trans Magn 45(10):46804683
10. Pster P-D, Perriard Y (2010) Very-high-speed slotless permanent-magnet motors: analytical
modeling, optimization, design, and torque measurement methods. IEEE Trans Ind Electron
57(1):296303
11. Komeza K, Dems M (2012) Finite-element and analytical calculations of no-load core losses
in energy-saving induction motors. IEEE Trans Ind Electron 59(7):29342946
12. Wang SH, Meng XJ, Guo NN, Li HB, Qiu J, Zhu JG et al (2009) Multilevel optimization for
surface mounted PM machine incorporating with FEM. IEEE Trans Magn 45(10):47004703
13. Barcaro M, Bianchi N, Magnussen F (2012) Permanent-magnet optimization in
permanent-magnet-assisted synchronous reluctance motor for a wide constant-power speed
range. IEEE Trans Ind Electron 59(6):24952502
14. Vese I, Marignetti F, Radulescu MM (2010) Multiphysics approach to numerical modeling of
a permanent-magnet tubular linear motor. IEEE Trans Ind Electron 57(1):320326
15. Bornschlegell AS, Pelle J, Harmand S, Fasquelle A, Corriou J-P (2013) Thermal optimization
of a high-power salient-pole electrical machine. IEEE Trans Ind Electron 60(5):17341746
16. Lee D-H, Pham TH, Ahn J-W (2013) Design and operation characteristics of four-two pole
high-speed SRM for torque ripple reduction. IEEE Trans Ind Electron 60(9):36373643
17. Flieller D, Nguyen NK, Wira P, Sturtzer G, Abdeslam DO, Merckle J (2014) A self-learning
solution for torque ripple reduction for nonsinusoidal permanent-magnet motor drives based
on articial neural networks. IEEE Trans Ind Electron 61(2):655666
18. Hasanien HM, Abd-Rabou AS, Sakr SM (2010) Design optimization of transverse flux linear
motor for weight reduction and performance improvement using response surface
methodology and genetic algorithms. IEEE Trans Energy Convers 25(3):598605
19. Hasanien HM (2011) Particle swarm design optimization of transverse flux linear motor for
weight reduction and improvement of thrust force. IEEE Trans Ind Electron 58(9):40484056
20. Lei G, Liu CC, Guo YG, Zhu JG (2015) Multidisciplinary design analysis for PM motors with
soft magnetic composite cores. IEEE Trans Magn 51(11). Article 8109704
21. Hua W, Cheng M, Zhu ZQ, Howe D (2006) Design of flux-switching permanent magnet
machine considering the limitation of inverter and flux-weakening capability. In: Proceedings
of 41st IAS annual meeting-industry applications conference, vol 5, pp 24032410
22. Liu CC, Zhu JG, Wang YH, Lei G, Guo YG, Liu XY (2014) A low-cost permanent magnet
synchronous motor with SMC and ferrite PM. In: Proceedings of 17th international conference
on electrical machines and systems (ICEMS), pp 397400
23. Fei W, Luk PCK, Shen JX, Wang Y, Jin M (2012) A novel permanent-magnet flux switching
machine with an outer-rotor conguration for in-wheel light traction applications. IEEE Trans
Ind Appl 48(5):14961506
24. Guo YG (2003) Development of low cost high performance permanent magnet motors using
new soft magnetic composite materials, UTS thesis (PhD)
25. Guo YG, Zhu JG, Watterson PA, Wei Wu (2006) Development of a PM transverse flux motor
with soft magnetic composite core. IEEE Trans Energy Conver 21(2):426434
26. Guo YG, Zhu JG, Watterson PA, Wei Wu (2003) Comparative study of 3-D flux electrical
machines with soft magnetic composite cores. IEEE Trans Ind Appl 39(6):16961703
27. Guo YG, Zhu JG, Dorrell D (2009) Design and analysis of a claw pole PM motor with molded
SMC core. IEEE Trans Magn 45(10):5824585
28. Lei G, Shao KR, Guo YG, Zhu JG (2012) Multi-objective sequential optimization method for
the design of industrial electromagnetic devices. IEEE Trans Magn 48(11):45384541
References 71
29. Lei G, Guo YG, Zhu JG et al (2012) System level six sigma robust optimization of a drive
system with PM transverse flux machine. IEEE Trans Magn 48(2):923926
30. Lei G, Zhu JG, Guo YG, Hu JF, Xu W, Shao KR (2013) Robust design optimization of
PM-SMC motors for Six Sigma quality manufacturing. IEEE Trans Magn 49(7):39533956
31. Lei G, Zhu JG, Guo YG, Shao KR, Xu W (2014) Multiobjective sequential design
optimization of PM-SMC motors for six sigma quality manufacturing. IEEE Trans Magn, 50
(2). Article 7017704
32. Liu CC, Zhu JG, Wang YH, Guo YG, Lei G, Liu XY (2015) Development of a low-cost
double rotor axial flux motor with soft magnetic composite and ferrite permanent magnet
materials. J Appl Phys, 117(17). Article # 17B507
33. Teng QF, Zhu JG, Wang TS, Lei G (2012) Fault tolerant direct torque control of three-phase
permanent magnet synchronous motors. WSEAS Trans Syst 8(11):465476
34. Teng QF, Bai J, Zhu JG, Sun Y (2013) Fault tolerant model predictive control of three-phase
permanent magnet synchronous motors. WSEAS Trans Syst 12(8):385397
35. Wang Y (2011) Investigation of rotor position detection schemes for PMSM drives based on
analytical machine model incorporating nonlinear saliencies, UTS thesis (PhD)
36. Wang TS (2013) Model predictive torque control of PMSM with duty ratio optimization for
torque ripple reduction, UTS thesis (Master degree)
37. Kim SY, Lee W, Rho MS, Park SY (2010) Effective dead-time compensation using a simple
vectorial disturbance estimator in PMSM drives. IEEE Trans Ind Electron 57(5):16091614
38. Lee J, Hong J, Nam K, Ortega R, Praly L, Astol A (2010) Sensorless control of
surface-mount permanent-magnet synchronous motors based on a nonlinear observer. IEEE
Trans Power Electron 25(2):290297
39. Genduso F, Miceli R, Rando C, Galluzzo GR (2010) Back EMF sensorless-control algorithm
for high-dynamic performance PMSM. IEEE Trans Ind Electron 57(6):20922100
40. Takahashi I, Noguchi T (1986) A new quick-response and high-efciency control strategy of
an induction motor. IEEE Trans Ind Appl 22(5):820827
41. Depenbrock M (1988) Direct self-control (DSC) of inverter-fed induction machine. IEEE
Trans Power Electron 3(4):420429
42. Buja GS, Kazmierkowski MP (2004) Direct torque control of PWM inverter-fed AC motors-A
survey. IEEE Trans Ind Electron 51(4):744757
43. Lai YS, Chen JH (2001) A new approach to direct torque control of induction motor drives for
constant inverter switching frequency and torque ripple reduction. IEEE Trans Energy
Convers 16(3):220227
44. Lascu C, Trzynadlowski A (2004) A sensorless hybrid DTC drive for high-volume low-cost
applications. IEEE Trans Ind Electron 51(5):10481055
45. Zhang Y, Zhu J, Xu W, Hu J, Dorrell DG, Zhao Z (2010) Speed sensorless stator flux oriented
control of three-level inverter-fed induction motor drive based on fuzzy logic and sliding mode
control. In: Proceedings of 36th IEEE IECON, pp 2926293
46. Zhang Y, Zhu J (2011) Direct torque control of permanent magnet synchronous motor with
reduced torque ripple and commutation frequency. IEEE Trans Power Electron 26(1):235248
47. Zhang Y, Zhu J (2011) A novel duty cycle control strategy to reduce both torque and flux
ripples for DTC of permanent magnet synchronous motor drives with switching frequency
reduction. IEEE Trans Power Electron 26(10):30553067
48. Zhang Y, Zhu J, Xu W, Guo Y (2011) A simple method to reduce torque ripple in direct
torque-controlled permanent-magnet synchronous motor by using vectors with variable
amplitude and angle. IEEE Trans Ind Electron 58(7):28482859
49. Wang TS, Zhu JG, Zhang YC (2011) Model predictive torque control for PMSM with duty
ratio optimization. In Proceedings of 2011 international conference on electrical machines and
systems (ICEMS), pp 15, 2023 August 2011
50. Miranda H, Cortes P, Yuz J, Rodriguez J (2009) Predictive torque control of induction
machines based on state-space models. IEEE Trans Ind Electron 56(6):19161924
51. Geyer T, Papafotiou G, Morari M (2009) Model predictive direct torque controlPart I:
Concept, algorithm, and analysis. IEEE Trans Ind Electron 56(6):18941905
72 2 Design Fundamentals of Electrical Machines and Drive Systems
52. Kouro S, Cortes P, Vargas R, Ammann U, Rodriguez J (2009) Model predictive controla
simple and powerful method to control power converters. IEEE Trans Ind Electron 56
(6):18261838
53. Morel F, Retif J-M, Lin-Shi X, Valentin C (2008) Permanent magnet synchronous machine
hybrid torque control. IEEE Trans Ind Electron 55(2):501511
54. Drobnic K, Nemec M, Nedeljkovic D, Ambrozic V (2009) Predictive direct control applied to
AC drives and active power lter. IEEE Trans Ind Electron 56(6):18841893
55. Zhang Y, Xie W (2014) Low complexity model predictive control-single vector-based
approach. IEEE Trans Power Electron 29(10):55325541
56. Zhang Y, Qu C (2015) Model predictive direct power control of PWM rectiers under
unbalanced network conditions. IEEE Trans Ind Electron 62(7):40114022
Chapter 3
Optimization Methods
3.1 Introduction
In Chap. 2, the design fundamentals, and various design analysis models for
electrical machines and drive systems have been investigated in terms of different
disciplines or subject domains, such as the analytical models or methods for
electromagnetic and thermal analyses, magnetic circuit model for electromagnetic
analysis, nite element model (FEM) for all electromagnetic, thermal and
mechanical analyses, and eld oriented control (FOC), direct torque control (DTC),
and model predictive control (MPC) algorithms for the control systems. While
some of them are physical analysis models which can reveal the basic operational
principles of the electrical machines and drive systems, the FEM is a kind of
numerical analysis model, which is widely used in the design optimization of
electrical machines to get a further understanding and illustrations for the eld
analysis. All these design analysis models can be employed for the performance
evaluation of electrical machines and drive systems.
On the other hand, solving an optimization problem consists of two main issues:
denition of optimization models and selection/development of (new) optimization
methods. Generally, an optimization model for a single-objective with m constraints
can be dened as
min : f x
s:t: gi x 0; i 1; . . .; m ; 3:1
x l x xu
where x, f and g are the design parameter vector, objectives and constraints,
respectively, xl and xu the lower and upper boundaries of x, respectively, and m is
the number of constraints. Typical optimization objectives for electrical machines
and drive systems are the minimization of cost, cogging torque, torque ripples,
overshoot, and maximization of output power and efciency. Popular constraints
are volume, mass, current density and temperature rises [1].
Theoretically, the optimization model (3.1) is always a strongly-constrained,
highly- nonlinear and high-dimensional problem for electrical machines and drive
systems. Many kinds of optimization algorithms have been employed to nd the
optimum for the above equation, such as the sequential quadratic programming
algorithm, genetic algorithm (GA), differential evolution algorithm (DEA), and
particle swarm optimization (PSO) algorithm [25]. Section 3.2 presents an over-
view for the classical and modern optimization algorithms for solving (3.1).
On the other hand, from the perspective of practical engineering applications, the
design optimization of electrical machines is actually a multi-objective problem as
there are many objectives which can be dened and different objectives can be
selected for different applications. For example, for applications in home appli-
ances, such as washing machines and refrigerators, the motor price and output
power may be the two most important issues; while for applications in hybrid
electric vehicles, the volume, power density and torque ripple are very important.
Therefore, multi-objective optimization design problems of electrical machines as
well as other electromagnetic devices have become a topic of much interest recently
[68].
Generally, a multi-objective optimization model with p objectives and m con-
straints can be dened as
min : ff1 x; f2 x; . . . fp xg
s:t: gi x 0; i 1; . . .; m ; 3:2
x l x xu
The solutions of (3.2) are often illustrated by Pareto optimal gure, which can be
obtained by using multi-objective optimization algorithms. Many multi-objective
3.1 Introduction 75
Many kinds of classic optimization algorithms have been introduced to solve the
constrained and nonlinear optimization problem (3.1). Some are gradient-based
algorithms, such as conjugate gradient algorithm, sequential quadratic program-
ming algorithm and augmented Lagrange multiplier method [5, 912]. Generally,
the rst or second order derivative or Hessians matrix is required in the imple-
mentation. To use these algorithms efciently, there are several constraints, such as
(a) The objective functions should be continuous and derivable;
(b) The objective functions and constraints can be expressed analytically; and
(c) The constrained optimization models have to be converted to unconstrained
forms for some initial gradient-based algorithms, e.g. the conjugate gradient
algorithm.
Analytical models or methods for electromagnetic, thermal and other disci-
plinary analyses should be constructed before the optimization. However, many
analysis models for electrical machines are based on FEM, and there is no analytical
expression for the optimization model. Therefore, various intelligent optimization
algorithms using non-analytical machine models have been employed, such as
those based on the GA and PSO algorithms.
76 3 Optimization Methods
3.2.2.1 GAs
GAs have been employed in science and engineering as adaptive intelligent algo-
rithms for solving practical problems. They are inspired by Darwins theory about
evolution. Solution to a problem provided by GAs is evolved. Figure 3.2 illustrates
a general optimization flowchart of GAs. As shown, the algorithm is started with a
set of population (represented by chromosomes). Solutions from one population
(known as parent) are taken and used to form a new population (known as offspring
or children) by three genetic operations, crossover, mutation and selection.
Solutions which are used to generate new solutions (offspring or children) are
selected in terms of their tness, which means that the more suitable they are, the
more opportunities they have to reproduce in the evolution process [14, 19, 20].
This is repeated until some conditions or criteria are satised, for example the
maximal iteration number. The outline of the basic GA is listed as follows.
(a) StartGenerate initial population of NP chromosomes;
(b) FitnessEvaluate the tness f(x) of each chromosome x in the initial
population;
(c) New populationCreate a new population by repeating following steps;
SelectionSelect two parent chromosomes from a population in terms of
their tness (the better tness, the bigger chance to be selected);
CrossoverForm a new offspring with a crossover probability over the
parents;
MutationMutate new offspring at each locus (position in chromosome)
with a mutation probability;
78 3 Optimization Methods
(d) ReplaceUse new generated population for a further run of algorithm; and
(e) TestIf the end condition is satised, stop, and return the best solution in
current population. Otherwise, go to step (b) and do the optimization loop till
convergence.
In the implementation, there are several algorithm parameters needed to dene.
They are the population size (NP), crossover probability (Pc), mutation probability
(Pm) and determination of the type of selection strategies. Generally, the population
size can be dened as 510 times of the problem dimension, Pc = 0.61.0, and
Pm = 0.0050.05. For the selection operation (regarding the problem of how to
select parents for crossover), this can be done in many ways. A popular one is to
select the better parents (assuming that the better parents will produce better off-
spring), and this is generally called elitism select strategy. It means that at least one
best solution is copied without changes to a new population, so that the best
solution found can survive to the end of run. Other selection strategies are roulette
wheel selection and rank selection methods [19, 20].
GAs have many advantages. For example, GAs work on the chromosome, which
is an encoded version of potential solutions parameters, rather than the parameters
themselves. On the other hand, they use tness score, which is obtained from
3.2 Optimization Algorithms 79
3.2.2.2 DEA
L U L
x0ji xj randji 0; 1 xj xj ; j 1; 2; . . .; D 3:3
L U
where xj and xj are the lower and upper boundaries of x. In details, it
can be expressed as
8
> L U L
>
> x01i x1 rand1i 0; 1 x1 x1
>
< x0 xL rand 0; 1 xU xL
2i
2i 2
..
2 2
3:4
>
>
>
> .
: 0 L U L
xDi xD randDi 0; 1 xD xD
where r1, r2, and r3 are three different numbers in [1, NP], that are
different from i, and F 2 0; 2 is the mutation factor. This mutation
method is called as DE/rand/1. There are several other situations:
where r4 and r5 are two different numbers in [1, NP], and subscript
best means the best one in the iteration t.
3.2 Optimization Algorithms 81
3.2.2.3 EDA
3.2.2.4 PSO
where w is the inertial weight factor, subscript D the dimension of parameter, pi the
local best vector of the tth particle, pg the global best vector, c1 and c2 are adjustable
social factors, r1 and r2 random numbers between 0 and 1, respectively, and is the
time step.
Figure 3.5 shows a flowchart of the PSO algorithm. The PSO algorithm has been
used in many applications and has had many improvements. Compared with GAs,
PSO is very much similar in many aspects. It is also a kind of evolutionary tech-
nique with its algorithm starting with a group of a randomly generated population,
using tness value to evaluate the population, updating the population, and
searching for the optimum with random techniques.
3.2 Optimization Algorithms 83
However, PSO does not have genetic operators like crossover and mutation.
Particles update themselves with the internal velocity. They also have memory,
which is important to the algorithm.
The information sharing mechanism in PSO is also signicantly different from
that of GAs. In GAs, the chromosomes share information with each other, so that
the whole population moves like a group towards an optimal area. In PSO, only pg
and pi give out the information to others, which is a one-way information sharing
mechanism. The evolution only looks for the best solution. Unlike GAs, in PSO, all
the particles tend to converge to the best solution quickly even in the local version
in most cases [11, 2830].
84 3 Optimization Methods
There are many kinds of multi-objective design optimization problems in the design
of electrical machines and other electromagnetic devices [7, 8, 3035].
Theoretically, the objectives in multi-objective optimization problems are always
conflicting. The improvement of an objective may result in performance decrease of
the other objectives. For example, the material cost and output power are two
important issues for designing the transverse flux machines. The improvement of
output power is often accompanied by the increase of material cost [7]. Therefore, it
is always impossible to achieve the optimum for each of these objectives, and the
corresponding optimal solutions are actually a compromise between these objec-
tives by making the objectives close to their optimums as much as possible. The
corresponding optimal solutions are called the Pareto optimal solutions.
Theoretically, the Pareto solutions are only acceptable solutions or non-inferior
solutions. The number of these solutions may be very large or even innite.
There are several conceptions which are widely mentioned in the multi-objective
optimization to dene the Pareto optimal solutions [8].
Denition 1 Given two vectors, x; y 2 Rq , we say that x y if xi yi for i = 1,,
q, and that x dominates y (denoted by x y if x y and x 6 y:
Figure 3.6 shows a particular case of the dominance relation in the presence of
two objective functions for a minimization situation.
Denition 2 We say that a vector of decision variables x 2 X is non-dominated in
X, if there does not exist another x0 2 X such that f x0 f x:
Denition 3 We say that a vector of decision variables x 2 F (F is the feasible
region) is Pareto-optimal if it is non-dominated in terms of F:
Denition 4 The Pareto optimal set is dened by
P fx 2 Fjx is Pareto-optimalg
PF ff x 2 Rq jx 2 P g
Figure 3.7 shows a particular case of the Pareto front in the presence of two
objective functions [8].
Different from the single objective optimization algorithms, the multi-objective
optimization algorithms have to provide a set of non-inferior solutions with large
population, and this set approaches the front of the global Pareto optimal solutions.
Those solutions should be uniformly distributed at the front of Pareto solutions as
much as possible.
Based on these basic principles, a number of multi-objective optimization
algorithms have been developed in the eld of evolutionary computation and have
been employed for the design optimization of multi-objective problems, such as
MOGA, NSGA and its improvement NSGA II, MPSO algorithm [6, 3640]. Three
of them, MOGA, NSGA II and MPSO will be introduced in the following sections.
3.3.2 MOGA
20
f2
15
10
-5
0 2 4 6 8 10 12 14 16 18
f1
NSGA stands for non-dominated sorting genetic algorithm, which was rst pre-
sented by Srinivas and Deb in 1994 [42, 43]. In this NSGA, a new method was
presented to classify individuals in layers before the selection is performed.
Individuals of the rst layer have the highest tness while the members of the last
layer have the smallest tness. Individuals from the rst layer produce more copies
than other layers in the next generation.
88 3 Optimization Methods
The NSGA II is an improved version of NSGA. It is one of the most efcient and
famous multi-objective evolutionary algorithms and has been widely applied in
many kinds of engineering multi-objective optimization problems. Figure 3.11
shows a flow chart of the algorithm. The method includes two important compo-
nents: the non-dominated sorting approach and the crowd comparison operator.
A detailed description can be found in Ref. [41].
3.3.4 MPSO
Figure 3.12 shows a general MPSO framework. MPSO has similar optimization
procedures as multiobjective evolution algorithms, such as MOGA and NSGA.
Many successful strategies, for example, external archive, have been introduced to
MPSO. On the other hand, tness evaluation is not a necessary step in MPSO, so
that the algorithm design can be simplied. However, a global optimal position
should be selected from the external archive for each particle, and this step is not
required for multi-objective evolution algorithms. There are many kinds of
improvements for MPSO, and many of them have been widely employed in the
design optimization of electrical machines. Detail descriptions of them can be found
in references [39, 40, 43].
3.4.1 Introduction
The above contents are about the single- and multi- optimization algorithms. As we
know, they are only one important component in the optimization as well as a direct
factor for the optimization efciency. Another important issue for the optimization
efciency is related to the optimization models. FEM has been widely used in the
design optimization of electrical machines. However, as mentioned previously, the
computational cost of FEM is always high, especially for complex structured
electrical machines, such as permanent magnetic (PM) soft magnetic composite
(SMC) machines [1, 7]. As an alternative, some approximate models (or surrogate
models) are employed in the practical engineering design problems to ease the
computational burden of optimization process, such as RSM and Kriging model
[44, 45].
Many research works have found that the optimization design based on
approximate models presents an effective way to solve the aforementioned prob-
lems. By using the design of experiments (DOE) technique and statistical analysis
methods, approximate models can be established as surrogate models for those
physical models, such as FEMs and circuit models, so as to reduce the high sim-
ulation cost in the iterative process of optimization. Meanwhile, approximate
models can degrade the nonlinearity of the practical problems, which can benet the
nding of the global optimal solution.
Generally, constructing an approximate model consists of the following two
steps.
Sampling: determine the required samples for constructing the approximate
models by using the DOE technique, X fx1 ; x2 ; . . .; xn jxi 2 RD g and their
responses fy1 ; y2 ; . . .; yn jyi 2 Rq g ; and
Modelling construction or tting: Fit the samples X, y with a suitable
approximate model, and test the model accuracy with some new samples.
There are four kinds of approximate models which have been widely used in
optimization design of electromagnetic devices, namely RSM, RBF model, Kriging
model, and ANN model. While RSM and RBF models are parametric models,
Kriging is a semi-parametric model, and ANN model a non-parametric model.
3.4.2 RSM
out. RSM can effectively replace the simulation and rapidly investigate tradeoffs
between various optimization tasks and conditions [46].
Generally, RSM ts the response data to lower-order quadratic polynomials
using the least-square method, and the sample data are usually obtained from the
uniform sampling method. Quadratic polynomials are generally used in electro-
magnetic problems and have been successfully developed for RSM in the form as
X
D X
D X
D
^yx b0 bi xi bij xi xj 3:14
i1 i1 j1;i j
where s are the regression parameters. To estimate the regression parameters, the
least square method (LSM) is commonly used to minimize the quadratic sum of
residual errors
X
n
Se yi ^yi 2 3:15
i1
Let the partial derivatives for all parameters be zero, and one can obtain a linear
equation system in the matrix form as the following
XT Xb XT y 3:16
b XT X1 XT y 3:17
However, there is a conflict between the accuracy of tting value and the step
size of sample data. A smaller step size may induce higher accuracy, but the
structural matrix may be poorly conditioned because the quadratic polynomials are
in general globally supported in LSM. The moving least square method (MLSM)
can overcome this defect. The centers of MLSM are randomly chosen and the
randomness is controlled by the point density and surface geometry [47].
There are two main reasons for why the RSM based on MLSM is better than
LSM. The rst one is the construction of tting functions that consist of a vector
92 3 Optimization Methods
function with coefcients and a basis function instead of the traditional polyno-
mials. Both are functions of sample points. The second is the introduction of
compactly supported domain in MLSM. Within MLSM, the response of y f x is
only determined by those samples in a small sub-domain around the point x, and
this sub-domain is called the compactly supported domain. Then the samples
outside this sub-domain do not have any effects on the response. Therefore, a
weighting function is dened in this compactly supported domain. If the weight
function is constant in the whole design space, MLSM is the traditional LSM. In
other words, LSM is a special case of MLSM.
In the MLSM, tting function ^yx can be approximated as a sum of linearly
independent function as the following
X
m
^yx pi xai x pT xax 3:19
i1
X
np
Jx wx xi ^yxi yxi 2
i1
3:21
X
np
2
wx xi p xi ax yxi
T
i1
It should be noted that the tting accuracy highly depends on the selection of
weight functions, which should be equal to 0 outside the compactly supported
3.4 Approximate Models 93
domain. The cube spline function is widely used to get the coefcient matrix by
minimizing the weighted square sum of error, which has the form as
8
< 2=3 4s2 4s3 s 1=2
Ws 4=3 4s 4s2 4s3 =3 1=2\s 1 ; 3:23
:
0 s[1
where s jjx xi jj=s0 ; and s0 is the radius of the compactly supported domain.
As a summary, the following lists the computational procedure of the RSM
based on MLSM:
(a) Generate samples by DOE technique, and calculate the responses (objectives)
of those samples;
(b) For each new sample x needing evaluation, implement the following steps:
determine the compactly supported domain for x;
count the sample size inside the supported domain;
calculate the model coefcient vector; and
compute the tting value of x;
(c) Draw the response surface by connecting all samples.
Compared with the RSM, RBF is also an empirical modeling approach for deter-
mining the relationship between various process parameters and responses with the
various desired criteria. RBF can effectively replace the time consuming simulation
(or measurement) and investigate very rapidly the tradeoffs between conflicting
performance criteria for optimization tasks. RBFs are commonly used in electro-
magnetic problems and have been successfully developed for constructing the
response surface [4850].
In general, the multivariate functions H : Rd ! R can be efciently evaluated
if they are expressible as univariate functions H Hjj jj of the Euclidean
norm jj jj; and such functions are called RBFs. With a set of scattered points
xi 1 i n ,the analytical expression of RBF can be given by
X
n
f x b j H x x j 3:24
j1
where Hr H x xj is the RBF, r x xj the Euclidian norm, and bj
(j = 1,2,,n) are the unknown parameters. Three most widely used RBFs in the
94 3 Optimization Methods
MQ : Hr r 2 c2 1=2 3:26
f xi yi ; i 1; 2; . . .; n 3:28
Then
8
>
> b H k x 1 x 1 k b2 H k x 1 x 2 k bn H k x 1 x n k y 1
> 1
< b1 H kx2 x1 k b2 H kx2 x2 k bn H kx2 xn k y2
.. 3:29
>
> .
>
:
b1 H kxn x1 k b2 H kxn x2 k bn H kxn xn k yn
Hb Y 3:30
2 3
H kx1 x1 k; H kx1 x2 k; . . .; H kx1 xn k
6 H kx2 x1 k; H kx2 x2 k; . . .; H kx2 xn k 7
6 7
where H 6 .. 7:
4 . 5
H kxn x1 k; H kxn x2 k; . . .; H kxn xn k
When there are no superposition points, H is a positive denite matrix and (3.30)
has an unique solution as
b H1 Y 3:31
However, RBFs are generally globally supported and poorly conditioned (sim-
ilar to the LSM). Although there are several remedies for these problems, such as
domain decomposition, preconditioning, and ne tuning of the variable parameter
of RBF, the compactly supported radial basis function (CSRBF) provides a
promising approach.
The centers of CSRBF are randomly chosen from the points and the randomness
is controlled by the point density and surface geometry. When the CSRBFs are
used, the evaluation of (3.31) will not run over the whole set of n summands and the
coefcient matrix will be sparse. The following two classes of CSRBF will be
studied in this book [51, 52]
3.4 Approximate Models 95
and
where r jjx xj jj=r0 , r0 is the radius of the compactly supported domain, and
1 r is given by
1r if 0 r 1
1 r : 3:34
0 otherwise
Given n sample points fx1 ; x2 ; . . .; xn g and their responses fyx1 ; yx2 ; . . .; yxn g;
for an input x, the response value y(x) of the Kriging model can be expressed as
yx f xT b zx 3:35
where
X
D
rxi ; xj expf ak jxik xjk j2 g 3:38
k1
and
respectively.
By using the maximum-likelihood estimation (MLE) method, the estimation of
r2 can be obtained as
1
^2 y F^bT R1 y F^b
r 3:44
n
Substituting the predicted b and r2 , i.e. (3.40) and (3.44), into the above
equation, the only unknown parameter ak can be estimated.
3.4 Approximate Models 97
y f Wx b 3:46
where jj jj is the Euclid norm. Gaussian function is always used as the transfer
function in this network.
Figure 3.13 shows the main three design steps for the approximate models. Firstly,
the samples are generated by using the DOE techniques. Then, the approximate
models are constructed, including the selections of model basis functions and tting
methods. Finally, the effectiveness/accuracy of the constructed approximate models
is veried [63, 64]. While the second step has been discussed previously, this
section presents the rst and the last steps.
98 3 Optimization Methods
Many DOE techniques have been developed, such as the full-factor design, Latin
hypercube design, orthogonal design, and central composite design. To improve the
tting accuracy, the full-factor design has been employed in many situations. For
example, consider an optimization problem with three parameters, and 5 levels are
dened for each parameter. In total, 53 = 125 samples are needed to construct an
approximate model. The Latin hypercube design is a sampling method with con-
straints. It rstly divides the initial design space into a number of non-overlapping
sub-spaces, and then sampling with equal probability is implemented in each
sub-space [57, 63, 64].
3.5 Construction and Verication of Approximate Models 99
The smaller the RMSE value is, the better the accuracy of tting model.
(b) Multiple correlation coefcient
The multiple correlation coefcient, R, or the coefcient of multiple correlation
is used to reveal the correlation between one factor/parameter with another one
or several other factors/parameters. It is dened as
It is a number that indicates how well the data t a kind of model sometimes
simply a line or curve. It is one of the indexes for evaluating the modelling accuracy
of approximate models. A high modelling accuracy is achieved when R2 approa-
ches to 1. However, this R2 coefcient can be easily affected by the number of
parameters. To avoid this undesired effect, the adjusted coefcient of determination
dened as
SSe =Ne D 1
R2adj 1 3:51
SST =Ne 1
can be employed.
100 3 Optimization Methods
Two examples of classical test functions with multiple local minimums will be
presented as follows to show the verications for the RSM, RBF (Gauss type) and
Kriging models.
The rst one is the Rastrigin function expressed in the form of
X
2
f1 x1 ; x2 x2i 10 cos2pxi 10; xi 2 5:12; 5:12 3:52
i1
The global minimum of this function is (0, 0), and the objective is 0. There are
99 local minimums around the global minimal point (0, 0). To have a clear illus-
tration, Fig. 3.14 shows the 3D surface prole of this function in a smaller region of
xi 2 2; 2 instead of the whole region [17, 41, 48, 63].
Figure 3.15 illustrates the RMSE curves for this Rastrigin function with the
RSM, RBF, and Kriging models, where the horizontal axis is the sample size for
each parameter. The minimal and maximal sample sizes are 5 and 33, respectively.
For the situation of minimal sample 5, the whole optimization region [5.12 5.12]
is rstly divided into 4 parts uniformly, and then ve points at [5.12, 2.56, 0,
2.56, 5.12] are sampled with a step size of 2.56. For the situation of maximal
sample 33, the whole optimization region [5.12 5.12] is divided into 32 parts with
equal length, and 33 points are sampled with a step size of 0.32. Thus, we have in
total 29 situations (from 5 to 33) for each approximate model. To verify the
accuracy of the constructed approximate model, 50 sample points are generated for
each parameter in [5.12, 5.12], with a step size of 0.21.
As shown in Fig. 3.15, all the three curves oscillate when the sample size is smaller
than 12. When the sample size is larger than 19, the RMSE values of the RBF and
Kriging models are smaller than those of the RSM, meaning that the RBF and Kriging
models are better than the RSM model in terms of the modelling accuracy.
X
2
f2 x1 ; x2 0:01 xi 0:54 30x2i 20xi ; xi 2 5:12; 5:12 3:53
i1
15
RMSE
10
0
5 10 15 20 25 30 35
Sample size for each parameter
0.6
R2adj
0.4
0.2 RSM
RBF
Kriging
0
5 10 15 20 25 30 35
Sample size for each parameter
102 3 Optimization Methods
This function has four minimal points located in each of the four quadrants,
respectively. The global minimum is 5.23 and is located at [4.45, 4.45] in the
third quadrant [63, 6567]. The other minimal points are 3.68 at [3.29, 3.29] (the
rst quadrant), 4.46 at [4.45, 3.29] (the second quadrant), and 4.46 at [3.29,
4.45] (the fourth quadrant).
Figure 3.18 illustrates the RMSE curves of the RSM, RBF, and Kriging models
for this test function. The horizontal axis is the sample size for each parameter.
Similar to the last example, the minimal and maximal sample sizes are chosen as 5
and 33, respectively. To verify the accuracy of the constructed approximate models,
100 sample points are generated for each parameter. As shown, the RMSE values of
the RBF and Kriging models are signicantly smaller than those of the RSM model.
It can also be seen that the Kriging model has fastest convergence and the smallest
RMSE values for almost all sample sizes.
Figure 3.19 illustrates the curves of the adjusted coefcient of these approximate
models for the second test function. As shown, when the sample size is larger than
0.8
RMSE
RSM
0.6
RBF
0.4 Kriging
0.2
-0.2
5 10 15 20 25 30 35
Sample size for each parameter
3.5 Construction and Verication of Approximate Models 103
0.6
R2adj
0.4
0.2 RSM
RBF
0 Kriging
-0.2
5 10 15 20 25 30 35
Sample size for each parameter
19, the R2adj values of the RBF and Kriging models approach 1 quickly, whereas
those of the RSM model remain around 0.4, which means the RBF and Kriging
models are better than the RSM model in terms of the modelling accuracy.
As shown by these two examples, the RBF and Kriging models are better than the
RSM model in terms of the modelling accuracy, and in both cases, the Kriging model
is better than the RBF model. It can also be seen that the modelling accuracy of these
models for the second test function is higher than that for the rst one because of the
stronger nonlinearity and much more number of minimums of the rst function.
3.6 Summary
This chapter presents a brief summary of the most commonly used numerical opti-
mization algorithms for electrical machines and drive systems, including the classic
gradient-based algorithms, modern intelligent algorithms, and multi-objective opti-
mization algorithms. Four kinds of approximate/surrogate models, namely RSM,
RBF, Kriging, and ANN models, have been presented with details. Two examples are
studied to illustrate the procedure of construction and verication of approximate
models.
References
1. Lei G, Wang TS, Guo YG, Zhu JG, Wang SH (2014) System level design optimization
method for electrical drive system: deterministic approach. IEEE Trans Ind Electron 61
(12):65916602
2. Yao D, Ionel DM (2013) A review of recent developments in electrical machine design
optimization methods with a permanent-magnet synchronous motor benchmark study. IEEE
Trans Ind Appl 49(3):12681275
104 3 Optimization Methods
3. Chechurin VL, Korovkin NV, Hayakawa M (2007) Inverse problems in electric circuits and
electromagnetics. Springer, New York
4. Coelho L, Alotto P (2007) Electromagnetic device optimization by hybrid evolution strategy
approaches. Int J Comput Math Electr Electron Eng (COMPEL) 26(2):269279
5. Nguyen TD, Lanfranchi V, Doc C, Vilain JP (2009) Comparison of optimization algorithms
for the design of a brushless DC electric machine with travel time minimization. In:
Proceedings of the 8th international symposium on advanced electromechanical motion
systems & electric drives (ELECTROMOTION), pp 16
6. Di Barba P (2010) Multiobjective shape design in electricity and magnetism. Lect Notes Elect
Eng 47
7. Lei G, Shao KR, Guo YG, Zhu JG (2012) Multiobjective sequential optimization method for
the design of industrial electromagnetic devices. IEEE Trans Magn 48(11):45384541
8. Reyes-Sierra M, Coello CAC (2006) Multi-objective particle swarm optimizers: A survey of
the state-of-the-art. Int J Comput Intell Res 2(3):287308
9. Dubas F, Sari A, Hissel D, Espanet C (2008) A comparison between CG and PSO algorithms
for the design of a PM motor for fuel cell ancillaries. In: Proceedings of the IEEE Vehicle
Power and Propulsion Conference, pp 17
10. Razik H, Defranoux C, Rezzoug A (2000) Identication of induction motor using a genetic
algorithm and a quasi-Newton algorithm. In: Proceedings of the IEEE International Power
Electronics Congress, pp 6570
11. Qteish A, Hamdan M (2010) Hybrid particle swarm and conjugate gradient optimization
algorithm. Proc Adv Swarm Intell 6145:582588
12. Popa D-C, Micu DD, Miron O-R, Szabo L (2013) Optimized design of a novel modular
tubular transverse flux reluctance machine. IEEE Trans Magn 49(11):55335542
13. Whitley D (1994) A genetic algorithm tutorial. Stat Comput 4(2):6585
14. Hasanien HM, Abd-Rabou AS, Sakr SM (2010) Design optimization of transverse flux linear
motor for weight reduction and performance improvement using response surface
methodology and genetic algorithms. IEEE Trans Energy Convers 25(3):598605
15. Hasanien HM (2011) Particle swarm design optimization of transverse flux linear motor for
weight reduction and improvement of thrust force. IEEE Trans Ind Electron 58(9):40484056
16. Chun JS, Kim MK, Jung HK (1997) Shape optimization of electromagnetic devices using
immune algorithm. IEEE Trans Magn 33(2):18761879
17. Campelo F, Guimaraes FG, Igarashi H et al (2005) A clonal selection algorithm for
optimization in electromagnetics. IEEE Trans Magn 41(5):17361739
18. Ho SL, Yang SY, Ni G et al (2005) A particle swarm optimization-based method for
multiobjective design optimizations. IEEE Trans Magn 41(5):17561759
19. http://www.obitko.com/tutorials/genetic-algorithms/ga-basic-description.php
20. Srinivas M (1994) Genetic algorithms: a survey. Computer 6:1726
21. Storn R, Price K (1997) Differential evolution - a simple and efcient heuristic for global
optimization over continuous spaces. J Global Optim 11:341359
22. Lei G, Shao KR, Guo Y, Zhu J, Lavers JD (2008) Sequential optimization method for the
design of electromagnetic device. IEEE Trans Magn 44(11):32173220
23. Lei G, Shao KR, Guo Y, Zhu J, Lavers JD (2009) Improved SOM for high dimensional
electromagnetic optimization problems. IEEE Trans Magn 45(10):39933996
24. Larranaga P, Lozano JA (2002) Esitimation of distribution algorithm: a new tool for
evolutionary computation. Kluwer Acadmic Publisher, Boston
25. Pelikan M (2005) Hierarchical Bayesian optimization algorithm: toward a new generation of
evolutionary algorithms. Springer, New York
26. Hauschild M, Pelikan M (2011) An introduction and survey of estimation of distribution
algorithms. Swarm and Evol Comput 1(3):111128
27. Campelo F, Guimaraes FG, Ramirez JA, Igarashi H (2009) Hybrid estimation of distribution
algorithm using local function approximations. IEEE Trans Magn 45:15581561
References 105
28. Ho SL, Yang SY, Ni GZ, Wong HC (2006) A particle swarm optimization method with
enhanced global search ability for design optimizations of electromagnetic devices. IEEE
Trans Magn 42(4):11071110
29. Al-Aawar N, Hijazi TM, Arkadan AA (2011) Particle swarm optimization of coupled
electromechanical systems. IEEE Trans Magn 47(5):13141317
30. Elhossini A, Areibi S, Dony R (2010) Strength Pareto particle swarm optimization and hybrid
EA-PSO for multi-objective optimization. Evol Comput 18(1):127156
31. Beniakar ME, Sarigiannidis AG, Kakosimos PE, Kladas AG (2014) Multiobjective
evolutionary optimization of a surface mounted PM actuator with fractional slot winding
for aerospace applications. IEEE Trans Magn 50(2), Article 7016404
32. dos Santos Coelho L, Barbosa LZ, Lebensztajn L (2010) Multiobjective particle swarm
approach for the design of a brushless DC wheel motor. IEEE Trans Magn 46(8):29942997
33. Alotto P, Baumgartner U, Freschi F, Kstinger A, Magele Ch, Renhart W, Repetto M (2008)
SMES optimization benchmark extended: introducing Pareto optimal solutions into TEAM22.
IEEE Trans Magn 44(6):10661069
34. Guimaraes FG, Campelo F, Saldanha RR, Igarashi H, Takahashi RHC, Ramirez JA (2006) A
multiobjective proposal for the TEAM benchmark problem 22. IEEE Trans Magn 42(4):1471
1474
35. Lebensztajn L, Coulomb JL (2004) TEAM workshop problem 25: a multiobjective analysis.
IEEE Trans Magn 40(2):14021405
36. Dias AHF, Vasconcelos JA (2002) Multiobjective genetic algorithms applied to solve
optimization problems. IEEE Trans Magn 38(2):11331136
37. dos Santos Coelho L, Alotto P (2008) Multiobjective electromagnetic optimization based on a
nondominated sorting genetic approach with a chaotic crossover operator. IEEE Trans Magn
44(6):10781081
38. Caravaggi Tenaglia G, Lebensztajn L (2014) A multiobjective approach of differential
evolution optimization applied to electromagnetic problems. IEEE Trans Magn 50(2), Article
7015404
39. Ho SL, Yang JQ, Yang SY, Bai YN (2015) Integration of directed searches in particle swarm
optimization for multiobjective optimization. IEEE Trans Magn 51(3) Article: 7000804
40. Al-Aawar N, Hijazi TM, Arkadan AA (2011) Particle swarm optimization of coupled
electromechanical systems. IEEE Trans Magn 47(5):13141317
41. Deb K, Pratap A, Agarwal S, Meyarivan T (2002) A fast and elitist multi-objective genetic
algorithm: NSGA-II. IEEE Trans Evol Comput 6(2):182197
42. Srinivas N, Deb K (1994) Multiobjective optimization using nondominated sorting in genetic
algorithms. Evol Comput 2(3):221248
43. Sooksaksun N (2012) Pareto-based multi-objective optimization for two-block class-based
storage warehouse design. Ind Eng Manage Syst 11(4):331338
44. Wang LD, Lowther DA (2006) Selection of approximation models for electromagnetic device
optimization. IEEE Trans Magn 42(2):12271230
45. Mendes MHS, Soares GH, Coulomb J-L, Vasconcelos JA (2013) Appraisal of surrogate
modeling techniques: a case study of electromagnetic device. IEEE Trans Magn 49(5):1993
1996
46. Dyck D, Lowther DA, Malik Z et al (1998) Response surface models of electromagnetic
devices and their application to design. IEEE Trans Magn 35(3):18211824
47. Liu WK, Li S, Belytschko T (1997) Moving least-square reproducing kernel methods
(I) Methodology and convergence. Comput Methods Appl Mech Eng 143:113154
48. Ishikawa T, Tsukui Y, Matsunami M (1999) A combined method for the global optimization
using radial basis function and deterministic approach. IEEE Trans Magn 35(3):17301733
49. Gutmann HM (2001) A radial basis function method for global optimization. J Global Opt 19
(3):201227
50. Lei G, Yang GY, Shao KR, Guo YG, Zhu JG, Lavers JD (2010) Electromagnetic device
design based on RBF models and two new sequential optimization strategies. IEEE Trans
Magn 46(8):31813184
106 3 Optimization Methods
51. Wu ZM (1995) Compactly supported positive denite radial basis functions. Adv Comput
Math 4:283292
52. Song KZ, Zhang X, Lu MW (2004) Meshless method based on collocation with consistent
compactly supported radial basis functions. Acta Mech Sin 20(5):551557
53. Lebensztajn L, Marretto CAR, Costa MC et al (2004) Kriging: a useful tool for
electromagnetic device optimization. IEEE Trans Magn 40(2):11961199
54. Xia B, Pham M-T, Zhang YL, Koh C-S (2013) A global optimization algorithm for
electromagnetic devices by combining adaptive Taylor Kriging and particle swarm
optimization. IEEE Trans Magn 49(5):20612064
55. Lim D-K, Woo D-K, Kim I-W, Ro J-S, Jung H-K (2013) Cogging torque minimization of a
dual-type axial-flux permanent magnet motor using a novel optimization algorithm. IEEE
Trans Magn 49(9):51065111
56. Kim JB, Hwang KY, Kwon BI (2013) Optimization of two-phase in-wheel IPMSM for wide
speed range by using the Kriging model based on Latin hypercube sampling. IEEE Trans
Magn 47(5):10781081
57. Lei G, Guo YG, Zhu JG, Wang TS, Chen XM, Shao KR (2012) System level six sigma robust
optimization of a drive system with PM transverse flux machine. IEEE Trans Magn 48
(2):923926
58. Lei G, Guo YG, Zhu JG, Chen XM, Xu W (2012) Sequential subspace optimization method
for electromagnetic devices design with orthogonal design technique. IEEE Trans Magn 48
(2):479482
59. Lophaven SN, Nielsen HB, Sondergaard J (2002) DACE: A MATLAB Kriging toolbox
version 2.0. Technical Report IMM-TR-200212. Technical University of Denmark,
Copenhagen
60. Dorica M, Giannacopoulos DD (2006) Response clustering for electromagnetic modeling and
optimization. IEEE Trans Magn 42(4):11271130
61. Marinova I, Panchev C, Katsakos D (2000) A neural network inversion approach to
electromagnetic device design. IEEE Trans Magn 36(4):10801084
62. Burrascano P, Fiori S, Mongiardo M (1999) Review of articial neural networks applications
in microwave computer-aided design. Int J RF Microw CAE 9(3):158174
63. Lei G (2009) Statistical analysis method for electromagnetic inverse problems . Ph D thesis.
Huazhong University of Science and Technology, Wuhan, China (in Chinese)
64. Wu CFJ, Hamada MS (2000) Experiments: planning, analysis and parameter design
optimization. Wiley, New York
65. Rashid K, Farina M, Ramirez JA et al (2001) A comparison of two generalized response
surface methods for optimization in electromagnetics. Int J Comput Math Electr Electron Eng
(COMPEL) 20(3):740752
66. Coulomb JL, Kobetski A, Caldora M, Costa et al (2003) Comparison of radial basis function
approximation techniques. Int J Comput Math Electr Electron Eng (COMPEL) 22(3):616629
67. Pahner U, Hameyer K (2000) Adaptive coupling of differential evolution and multiquadrics
approximation for the tuning of the optimization process. IEEE Trans Magn 36(4):10471050
Chapter 4
Design Optimization Methods
for Electrical Machines
Abstract This chapter presents the design optimization methods for electrical
machines in terms of different optimization situations, including low- and
high-dimensional, single- and multi- objectives and disciplines. Firstly, the tradi-
tional design optimization methods are briefly reviewed, and the challenges pre-
sented. Then, ve new types of design optimization methods are presented to improve
the optimization efciency of electrical machines, particularly those complex struc-
tured permanent magnet machines, in terms of different optimization situations. They
are (a) a sequential optimization method for design optimization of low-dimensional
problems of electromagnetic devices including electrical machines, (b) a
multi-objective sequential optimization method for engineering multi-objective
problems, (c) a multi-level design optimization method (or sequential subspace
optimization method) for high dimensional problems, (d) a multi-level genetic
algorithm for high dimensional optimization problems as well, and (e) the
multi-disciplinary design optimization method. Design examples with detailed
experimental and optimization results are illustrated for each optimization method.
Keywords Design optimization Sequential optimization method Multi-objective
sequential optimization method Multi-level optimization method Multi-level
genetic algorithm Multi-disciplinary optimization method Permanent magnet
motors
4.1 Introduction
Design optimization methods actually consist of two parts, design methods with
analysis models, and optimization methods with algorithms. In Chap. 2, the popular
design methods and analysis models for electrical machines and drive systems are
reviewed. As shown, the design of electrical machines is a complex
multi-disciplinary or multi-physics problem, including electromagnetics, thermotics,
mechanics and control, and each discipline has its own design methods and analysis
models [13]. For example, the electromagnetic design is mainly based on analytical
model, magnetic circuit model and nite element model (FEM). The thermal design
is mainly based on FEM and thermal network model. The control system design is
mainly based on the topologies of power electronic circuits and control algorithms,
such as the eld oriented control algorithm and the direct torque control algorithm.
In Chap. 3, the popular optimization algorithms and approximate models used in
optimization of electrical machines as well as other electromagnetic devices are
discussed. The optimization algorithms include the classical gradient-based algo-
rithms and modern intelligent algorithms, such as genetic algorithms (GA), dif-
ferential evolution algorithm (DEA), and multi-objective genetic algorithms
(MOGA). The approximate models include the response surface model (RSM),
radial basis function (RBF) model, compactly supported radial basis function
(CSRBF) model, and Kriging model [46].
The general procedures for design optimization of electrical machines are listed
as follows:
(1) Determine the design analysis model for the investigated machine, including
selection of material, motor type and topology, and develop the
multi-disciplinary analysis model, for example, electromagnetic-thermal cou-
pled model based on FEM.
(2) Establish the optimization model, including the denition of objectives (such
as maximization of output power and efciency, and minimization of cost),
constraints (such as volume, mass and temperature rises) and design param-
eters (such as material and structural parameters). The optimization regions of
these parameters and their types, such as discrete and continuous, have to be
dened in this step as well.
(3) Select an optimization method and optimize the established optimization
model. There are many kinds of available optimization methods. The most
popular one is the direct optimization method, which uses an optimization
algorithm to optimize the optimization model established on FEM.
(4) Validate the effectiveness of the obtained optimal solutions by experiments or
other ways. If necessary, update the optimization model in terms of the
experimental results and do the optimization again.
In general, there are two kinds of classical optimization methods, and they are
(1) The direct optimization methods, and
(2) The optimization methods based on approximate models.
The direct optimization methods use the optimization algorithms to directly
optimize the physical models, such as the analytical model, magnetic circuit model,
and FEM for the design of electrical machines. For example, the conjugate gradient
algorithm and sequential quadratic programing algorithm have been employed to
4.2 Classical Optimization Methods 109
optimize several kinds of motors based on the analytical models for electromagnetic
analysis. The most popular formula of this type of methods is the combination of
intelligent algorithms and FEM, such as GA&FEM, DEA&FEM and
MOGA&FEM [2, 711].
This approach can present global optimal design schemes for electrical machines.
However, the computational cost of this kind of optimization methods is always
huge due to the extensive computational burdens of FEM, practically for some
complex structured electrical machines requiring 3D FEM, such as the transverse
flux machines (TFM) and claw pole motors. Most importantly, the computational
cost will increase greatly with the increase of problem dimension. For example, for a
motor design problem with four parameters, about 4,000 (4 5 200, where 4 5 is
the population size and 200 the general iteration number of GAs) FEM samples are
needed if the GA&FEM method is applied. However, if a motor has 10 design
parameters, about 20,000 (10 5 200) FEM samples are required, which is a huge
computational burden for many situations.
An effective way to solve this problem is to use the second kind of optimization
method: the optimization based on approximate models. This method replaces the
FEM with a kind of approximate model, such as RSM and Kriging, so as to form
the optimization method, GA&RSM and GA&Kriging. As mentioned in Chap. 3,
the approximate models can degrade the nonlinearity of the optimization problem.
Therefore, the optimization efciency can be improved signicantly.
This method is relatively simple to implement. However, the optimization
accuracy is an important problem for this kind of methods, particularly for high
dimensional problems. Actually, it is very hard or impossible to replace the FEM
with approximate models for high dimensional problems, because they cannot
approximate high dimensional problems with sufcient accuracy by using limited
number of samples. As we know, the rst step in the construction of approximation
models is to obtain the initial samples by using the design of experiments
(DOE) technique. If 5 samples are required for each parameter in a design problem
with ten parameters, 510 FEM samples are required in total, which is greater than
that required by direct optimization method of GA and FEM [1].
Therefore, the classical or traditional direct optimization methods based on both
FEM and the approximation models have challenges for the design optimization of
electrical machines, especially for the high dimensional design optimization
problems. The following sections will present several new optimization methods.
simultaneous updating process about them. This is one of the main reasons why the
traditional optimization methods cause huge computational costs. To overcome this
problem, a sequential optimization method (SOM) was developed to reduce the
computational cost for single-objective low-dimensional optimization problems of
electromagnetic devices in 2008 [12, 13].
Figure 4.1 illustrates the optimization strategy of SOM. Assume that square
ABCD is the initial design space and point S is the optimal point. The traditional
direct optimization method, e.g. GA&FEM, searches the whole design space for the
optimal point by an iterative process, and many new samples in the population are
required in the whole design space. However, the optimal point is located in a small
subspace (shown as the shaded rectangle) around point S. If we can nd this
interested subspace and sample more points in it instead of the whole design space,
the optimization efciency would be improved greatly. Based on this idea, the main
question is how to nd this subspace efciently. SOM is a method to deal with this
problem. It can reduce the design space step by step. As shown in Fig. 4.1, it
reduces the design space from ABCD to ABCD in the rst step, then to ABCD
in the second step, and so on until the shaded rectangle is reached in the last step.
Figure 4.2 depicts a brief flowchart of SOM. Basically, SOM can be regarded as
a space-to-space optimization strategy compared with the point-to-point optimiza-
tion strategy of the traditional intelligent algorithms, such as GA and DEA. SOM
consists of two optimization processes, coarse and ne optimization processes. The
main aim of the former is to reduce the design space to a small space (the shaded
rectangle as shown in Fig. 4.1). The purpose of the latter is to update the model in
the local space and nd the optimal solutions [12, 13].
SOM is conducted in the following six main steps:
Step 1: Dene optimization model, such as objectives, constraints and design
parameters. Select the approximate models and optimization algorithms
that will be used in the SOM, and dene the algorithm parameters, such as
population, genetic operators and maximal iteration number.
4.3 Sequential Optimization Method 111
Step 2: Use the coarse optimization process (COP) to reduce the design space.
The traditional optimization methods based on approximate models will
be employed in this step. The DOE is required in this method to generate
samples for the construction of approximate models.
k k
In general, one can assume that the design space of the kth COP is xli ; xui ,
k
i 1; 2; . . .; D, where D is the dimension, li the interval of the ith variable,
k k1 k
hi hi =2 the step size, and Ni the sample size, respectively. The sample set
S(k) can be obtained by using a DOE technique, e.g. the full-factor design. Based on
these samples, the approximate model can be constructed. Finally, through the
k k
optimization of the model, the current optimal point xo fxoi ji 1; 2; . . .; Dg
and objective f k can be obtained.
Note that the step size must ensure that the minimum number of sample points is
no less than 3. Otherwise, a singular matrix may appear in the matrix inversion
process of model construction.
112 4 Design Optimization Methods for Electrical Machines
k 1
Step 3: Terminate COP by the step size. If hi =li \d, where d is a positive
constant and can be a value in [1 %, 5 %], stop COP and go to step 5.
Otherwise, go to step 4.
Step 4: Reduce the current design space with the obtained optimal value. Under
the boundary condition of design space, the design space of next step can
be updated as follows:
8 k
>
> k 1 k x Dl
>
< xli max xli ; round oiDh Dh
k ; 4:1
>
> x Dl
> k 1
: xui
k
min xui ; round oi Dh Dh
k k
where function round(x) rounds x to the nearest integer, Dl li =nl, Dh hi =nh,
and nl and nh are the reduction factors.
For a practical problem, three parameters, nl, nh and N, can be used to determine
the construction of an approximate model. nl can be 4, 6 or 8, meaning that the
corresponding intervals of reduced space are 1/2, 1/3 and 1/4 of the former space,
respectively. nh can be 2, 4, or 8, which let function round(x) round x with 1/2, 1/4
and 1/8 of the current step size, respectively, and N can be 3, 4 or 5 for a stan-
dardization space of [0, 1].
To nd the best values of these parameters, the Monte Carlo analysis
(MCA) method is employed. Firstly, assume that the initial design space is [0, 1]
and nl is 4. Then, generate 106 random numbers by using the Monte Carlo method
as the optimal results of an approximate model. Thereafter, use these numbers to
reduce the design space under N = {3, 4, 5} and nh = {2, 4, 8}. The target of this
analysis is to compare the mean of the errors between . the current optimal results
k k
and the average of the reduced space xli xui 2. Table 4.1 lists the MCA
results for the case of nl = 4. Similarly, Tables 4.2 and 4.3 list the results for the
cases of nl = 6 and nl = 8, respectively.
As shown, the mean error decreases with the increase of N and nh for all three
cases. Thus, in the later implementation, the default value of nh is 8 and N is 5. The
A test function shown in Chap. 3 is used here to verify the efciency of the
proposed SOM [12], which has the form as
X
2
f2 x1 ; x2 0:01 xi 0:54 30x2i 20xi ; xi 2 5:12; 5:12 4:3
i1
The parameters in SOM are = = 5 % and N(1) = [5, 5]. The DEA is used as
the optimization algorithm in this example [14]. The algorithm parameters are
chosen as 0.8 for the mutation scaling factor, 0.8 for the crossover factor, 1000 for
the maximum number of iteration, and 100 for the maximum stall generation, which
is selected for the stop criterion [12, 13].
Table 4.4 tabulates the optimization results obtained by SOM with ve different
approximate models. They are the general RSM model based on least square
method, improved RSM model based on moving least square (MLS) method,
Gauss RBF model, CSRBF 1 model, and Kriging model. As shown, the obtained
optimal solutions by using the SOM are of good accuracy compared with the exact
values.
the efciencies of different optimization methods [12, 13, 1519]. This problem
consists of several cases, including low and high dimensional, single- and
multi-objective, and discrete as well as continuous parameter optimization cases.
Figure 4.3 illustrates the optimization structure of this SMES. As shown, it consists
of two solenoids and there are eight design parameters.
For this benchmark problem, there are three optimization objectives as listed
below:
(1) The expect value of the stored energy, E, in this SMES is 180 MJ;
(2) The mean stray elds, Bstray, should be as small as possible; The value of Bstray
21
P jBistray j2 1=2
can be calculated by Bstray 21 , where Bistray is the magnetic flux
i1
density evaluated along 21 equidistant points on lines a and b, and
(3) The magnetic eld should maintain the superconducting condition of the
storage. The superconducting material employed in this SMES is NbTi.
Figure 4.4 illustrates the critical curve of this material. To ensure the super-
conducting condition, the current density, Ji, in the solenoids and magnetic
eld density in the storage must follow the following constraint
30
20
10
0
0 1 2 3 4 5 6
Magnetic field density |B| [T]
Table 4.5 Design parameters Par. Unit Min Max Step size xed
and values of SMES for
discrete case R1 m 2.0
R2 m 2.6 3.4 0.01
h1/2 m 0.8
h2/2 m 0.204 1.1 0.007
d1 m 0.27
d2 m 0.1 0.4 0.003
J1 A/mm2 22.5
J2 A/mm2 22.5
Table 4.6 lists the optimization results for this SMES by using different opti-
mization methods, including direct optimization method and SOM based on dif-
ferent approximate models. The parameters used in SOM are d e 2:5 % and
4.3 Sequential Optimization Method 117
0.45
0.4
0.35
0.3
0.25
1 2 3 4 5
Iterative number of SOM
N 1 5; 5; 3. The best solution of this problem given by the TEAM workshop is
also listed in the table. Figure 4.5 illustrates the convergence rates of SOM by using
different approximate models. The rst three iteration processes are for the COPs.
The last one or two iteration processes are for the FOPs. The following conclusions
can be drawn from the table and gure:
(1) By using the direct optimization method based on FEM and DEA, the obtained
optimal result is 180.00 MJ for the stored energy (slightly better than the TEAM
value, 179.86 MJ), 3.83 T for the maximum magnetic flux density, and
1.0323 mT for the mean stray flux density (slightly higher than the TEAM value,
0.9084 mT). The objective is 0.3441, which is much higher than the TEAM
objective, 0.3034. To obtain this optimal result, 2310 FEM samples are used.
(2) By using the SOM based on RSM, the obtained optimal result is 179.68 MJ
for the stored energy, 0.9051 mT for the mean stray flux density, and 0.3049
for the objective. The obtained objective is lower than that obtained by the
118 4 Design Optimization Methods for Electrical Machines
0.4
0.35
0.3
0.25
d2
0.2
0.15
0.1
1
0.8
0.6 3.4
3.2
0.4 3
2.8
h2/2 0.2 2.6 R2
direct optimization method and close to the TEAM value. Most importantly,
this method requires only 164 FEM samples to obtain the optimal result.
Figure 4.6 depicts these samples using +, and to denote the points
sampled from three COPs, and o the points sampled in the FOP. As shown,
the samples are distributed non-uniformly in the whole design space, and the
sampling processes approach the nal optimal solutions step by step.
Figure 4.7 illustrates the magnetic flux density distribution in the SMES for this
optimal solution. As shown, the maximal flux density is 4.702 T, which is very
close to 4.70 T, the value listed in Table 4.6.
Therefore, SOM based on RSM has high optimization efciency while main-
taining the accuracy of optimum.
(3) The detailed optimization results of SOM based on other models are shown in
Table 4.6. As shown, all optimal objectives obtained from SOM with different
models are slightly higher than the TEAM value and smaller than that of DEA.
Both the SOM based on CSRBF1 model and the SOM based on Kriging
model used only 157 FEM samples, which is 6.8 % (=157/2310) of that used
by the DEA method.
4.3 Sequential Optimization Method 119
In summary, all solutions obtained by the SOM based on various models are
better than that of the direct optimization or DEA method, and there is not much
difference between the approximate models. In other words, SOM has high opti-
mization efciency whereas the model type does not affect its optimization accuracy
very much. Unlike the case presented in Chap. 3 that MLS is better than LSE for
RSM and CSRBF is better than RBF, they do not have much difference within the
framework of SOM. On the other hand, the MLS and CSRBF require extra
parameters, which make them more complex compared with RSM, RBF and
Kriging. Therefore, RBF and Kriging are two good models for SOM.
As mentioned in the above section, space reduction strategy plays an important role
in SOM. The design purpose of the former space reduction strategy is to minimize
the distance between the mean of next design range and the optimal result [12]. It is
accurate from the point of view of distance minimization, but it has not considered
the issue that reduces the number of FEM samples effectively. To make full use of
most points sampled in the last set, we present a new space reduction strategy in this
section.
k k
Assume xk xli ; xui is the boundary of the ith variable in the kth opti-
mization process, lk the interval, hk the step size, N k the number of sample
120 4 Design Optimization Methods for Electrical Machines
k
points, Sk the sample set, and xo and f k are the optimal result and corresponding
function value, respectively. The new space reduction strategy consists of the fol-
lowing two steps.
Reduction step:
( ! )
k
k 1 k xoi Dl
^xli max xli ; round Dh ; 4:7
Dh
( ! )
k
k 1 k xoi Dl
^xui min xui ; round Dh : 4:8
Dh
Correction step:
!
k 1 k k
k 1 k 2^xli xli hi
xli xli round k
; 4:9
h 2
i
!
k 1 k k
k 1 k 2^xui xli hi
xui xli round k
; 4:10
h 2
i
k
where round(x) is a function to round x to its nearest integer, Dl li =nl ,
k
Dh hi =nh , and nl and nh are the two reduction factors with defaults 4 and 8,
respectively [20].
To check the efciency of the new method, a comparison with the former space
reduction strategy is conducted.
As an example, assume that the initial design space is [0, 1], N = 6, and the
uniform sampling method is used. The rst sample data is then S(1) = {0.0, 0.2, 0.4,
0.6, 0.8, 1.0}, and the optimal value is supposed as 0.35. By the former space
reduction strategy, the next design space is [0.1, 0.6], and the next sample set
S(2) = {0.1, 0.2, 0.3, 0.4, 0.5, 0.6}. Thus, 3 sample points have been sampled in S(1),
or 50 % computation cost is saved. On the other hand, by the new space reduction
strategy, the next sample space is [0.2, 0.6], and S(2) = {0.2, 0.3, 0.4, 0.5, 0.6},
resulting in a save of 60 % computation cost, which is better than the former
strategy.
As another example, if the optimal value is 0.3, the next sample space is [0.05,
0.55] by the former strategy, and S(2) = {0.05, 0.15, 0.25, 0.35, 0.45, 0.55}. Thus,
no sample point has been sampled. By the new strategy, however, the next sample
space is [0.0, 0.6], and S(2) = {0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6}. Thus, 4 sample
points have been sampled in the S(1), resulting in a reduction of 57.14 % compu-
tation cost.
Therefore, the new strategy is more effective than the former one.
Table 4.7 shows mean saving rates of sample points for the former and new
space reduction strategies by using MCA. For each strategy and every sample
4.3 Sequential Optimization Method 121
number N, 106 random numbers are generated as the current optimal points by
Monte Carlo method. Then, the mean saving rate is obtained for each case. As
shown, all saving rates by the new strategy are more than 50 %, which are higher
than those obtained by the former strategy [20].
Table 4.8 lists the optimization results of the SMES by using two kinds of RBF
models. As shown, the Gauss RBF model gives the best result, and the
multi-quadrics (MQ) RBF model requires the least number of FEM samples.
Furthermore, either RBF model requires less than 1/10 of the FEM samples needed
by DEA.
(a) (b)
Fig. 4.8 a Molded stator of a PM claw-pole motor, and b FEM region of a pole of a stack
Regarding the optimization of this motor, the objective is chosen to minimize the
material cost while maximizing the output power or torque at 3000 rev/min. The
outer radius and axial length, or volume, of the motor is xed in the optimization.
The material cost mainly includes the costs of PM, copper, SMC core, and steel.
4.3 Sequential Optimization Method 123
Three constraints are also considered. The optimization model can be dened as the
following
Cost P0
min: f x
C0 Pout
s.t: g1 x 0:78 g 0 4:11
g2 x 60 Pout 0
g3 x Jc 4:5 0
where C0 and P0 are the cost and output power of the initial prototype, 0.78 and 60
the rated values of efciency () and output power (Pout), respectively for the initial
design, and the last constraint is the current density, Jc, of winding, which should be
no more than 4.5 A/mm2 in terms of its specications [22].
From previous design experience, three parameters Rsi, bs and hrm, are important
to the motor performance. Therefore, they are selected as the optimization variables.
Table 4.10 lists the optimization results by using improved SOM and Kriging
model. As shown, the obtained output power and material cost are 98 W and $8.99,
which are better than those of the initial design (60 W and $14.18). Meanwhile,
only 197 FEM samples are reuqired by improved SOM to obtain the optimal
124 4 Design Optimization Methods for Electrical Machines
solution. Compared with the samples required by direct optimization method, for
example, DEA requiring about 3 5 170 or 2550 FEM samples to obtain the
optimization results, where 3 5 is the population size and 170 the average iteration
number, 92.27 % of FEM computation cost has been saved. Therefore, the
improved SOM is efcient for low-dimensional electromagnetic design problems.
min : ff1 x; f2 x; . . .; fp xg
s.t: gi x 0; i 1; 2; . . .; m ; 4:12
xl x xu ; x x1 ; . . .; xD T
where p, m and D are the numbers of objectives, inequality constraints and vari-
ables, respectively.
In general, the solutions of a multi-objective optimization problem can be
illustrated as a Pareto optimal set, and its front is not a point, but a continuous or
non-continuous curve or surface. Thus, we should pay more attention to the sub-
space as shown in Fig. 4.12 which includes all these Pareto points and seek for a
new method to sample more points in this subspace rather than the total design
space. This approach may improve the modeling efciency because it includes the
investigation of model characteristic.
Figure 4.13 shows the flowchart of MSOM, which mainly includes the following
three steps [23]:
(1) Generate an initial sample set S(0), and construct Kriging models for all FEMs
in the design optimization problem (4.12) to get a Kriging multi-objective
optimization model.
(2) Optimize the Kriging multi-objective optimization model with a
multi-objective optimization algorithm, for example, NSGA II, and get the
Pareto optimal solutions P(k). Then compute the root mean square error
(RMSE) of the obtained Pareto points for each model. If all the RMSEs are no
more than a constant , output the solutions; otherwise go to the next step.
(3) Update the sample set S(k) and Kriging model. As the constructed models are
getting more and more accurate through the optimization process, the true
Pareto solutions are probably located around the current P(k). To improve the
modeling efciency, a modied central composite design (CCD) sampling
method is presented to update the sample sets.
The CCD is a classic sampling method for the construction of RSM. It divides the
samples into two subsets, one for the property estimation of the linear term, and the
4.4 Multi-objective Sequential Optimization Method 127
other for the curved surface. It is claimed to be superior in the modeling of RSM [27].
Considering that RSM was used as the determined term in the Kriging model, a
modied CCD sampling method has been presented in our previous work.
Figure 4.14 shows two illustrations of the proposed method for (a) a two-variable
case, and (b) a three-variable case. The circle points in the gure are the Pareto points.
The triangle points (with number of 2D) are sampled by the two levels full factor
design method. The square points (with number of 2D) in the axial direction are the
peaks of circumscribed circle or sphere. The relationship between the diameter (d) of
p
the circle (sphere) and the side length (l) of the square (cube) is d l D. For the kth
optimization step, the side length is dened as half of the step size in the last step.
(a) (b)
Figure 4.16 illustrates the Pareto optimal solutions obtained by two methods, the
direct function optimization with NSGA II and the proposed MSOM. For the latter
case, only one model updating process is needed to get the nal Pareto solutions. As
shown, the Pareto solutions of this function are separated to two parts. It is
non-continuous on the whole. The Pareto front of the 2nd Kriging model (k = 2 in
MSOM) ts that from the true function very well. Figure 4.17 illustrates the total
sample points sampled by MSOM which includes only 109 points. As shown, the
obtained samples are non-uniformly distributed in the whole space and include
more points in the subspace to which the Pareto points belong, so that the sampling
efciency can be improved by the proposed method.
10
0 2 4 6 8 10 12 14 16 18
f1
4.4 Multi-objective Sequential Optimization Method 129
x2
0
-1
initial samples
-2 1st Pareto points
updated samples
-3
-3 -2 -1 0 1 2 3
x1
In our previous work, a PM TFM with soft magnetic composite (SMC) stator core
was developed [1, 2932]. A multi-objective design optimization is conducted here
for the broad industrial applications of this machine. A PM-SMC TFM prototype
was illustrated previously in Fig. 2.3, Chap. 2 This machine was initially designed
to deliver an output power of 640 W at 1800 rev/min. Table 2.1 tabulates the main
dimensions. Figure 2.6 illustrates the FEM model used in ANSYS. The computa-
tion of FEM of this machine is very time-consuming as 3D FEM is required for the
performance evaluation.
To get reliable analysis and optimization results, the analysis model based on
FEM (as shown in Fig. 2.6) should be veried by experimental results. Figures 4.18,
4.19, 4.20 and Table 2.2 show the calculated and measured key parameters for this
machine. Figure 4.18 shows the measured motor speed against output torque with
different DC link voltages. Figure 4.19 illustrates the measured electromotive force
(EMF) waveforms at 1800 rev/min. The measured motor back EMF constant is 0.
244 Vs, which is 1 % lower than the calculated value of 0.247 Vs. The calculated
-40
Fig. 4.20 Measured curves 900 Pin & Pout (W) 0.9
for input, output powers and
motor efciency in terms of Efficiency
torque 700 0.7
Pin
500 Pout 0.5
300 0.3
100 0.1
phase resistance and inductance, and maximal cogging torque are 0.310 , 6.68 mH
and 0.339 Nm, respectively, which are very close to the measured values of 0.305 ,
6.53 mH and 0.320 Nm. Figure 4.20 shows the measured curves of the input power,
output power, and efciency against the output torque. It is found that the estimated
performance parameters calculated from FEM-based method are well aligned with
the experimental results, such as the inductance and cogging torque. Therefore, all
the experimental results have veried the effectiveness of this FEM-based analysis
method, and it is reliable to be used for optimizing the electrical machine under
investigation.
The multi-objective optimization model of this machine can be dened as
f1 x = Cost(PM) + Cost(Cu)
min :
8 f2 x 640 Pout
> g1 x 0:795 g 0;
> 4:14
<
g2 x 640 Pout 0;
s.t:
>
> g x sf 0:8 0;
: 3
g4 x Jc 6 0:
4.4 Multi-objective Sequential Optimization Method 131
790
730
700
FEM
670 MSOM-1st model
MSOM-3rd model
640
28 30 32 34 36 38 40
Cost [$]
where and Pout are the efciency and output power of the machine, respectively, sf
is the winding ll factor, and Jc in A/mm2 is the current density of the copper wire.
The rst objective cost mainly consists of the costs of PM and copper winding.
Four parameters are selected as the optimization variables in this work. They are
circumferential angle and axial width of PM, and the number of turns and diameter
of copper wire winding. These are the signicant parameters for the objectives from
our previous design experience [30].
For this problem, three model updating processes (k = 3 in MSOM) are needed
for the MSOM to get the nal Pareto solutions as shown in Fig. 4.21, which
includes the initial and the last Pareto points of MSOM. Figure 4.21 also illustrates
the Pareto points obtained by the direct optimization of FEM with NSGA II. As
shown, the Pareto front from the MSOM ts that from FEM very well. Moreover,
the needed FEM sample points of MSOM are only 556, which is about 6 % that by
the direct multi-objective optimization of FEM, in which about 10,000 FEM
samples are needed [23].
On the other hand, it is impossible to replace the FEM with approximate models,
such as RSM and Kriging model, because they cannot approximate highly
dimensional problems with sufcient accuracy by using reasonably small number of
samples. For example, the rst step in the construction of approximation models is
to use the DOE technique to obtain the initial samples. If 5 samples are required for
each parameter, in total, 510 FEM samples are required, which is greater than those
required by direct optimization method of GA and FEM.
Therefore, the traditional direct optimization method based on FEM and the
approximation models cannot solve the highly dimensional design optimization
problems. To solve these problems, a multi-level optimization method was pre-
sented for electrical machines and drive systems and other electromagnetic devices
in our previous work [22, 33, 34]. The main idea of the multi-level optimization
method is that the high dimensional design space can be divided into two or several
low dimensional design spaces in terms of the order of their sensitivities. The
detailed discussion of multi-level optimization method will be presented in the next
section and Chaps. 5 and 6. This section presents a brief investigation for sensitivity
analysis (SA) techniques. In general, there are four types of techniques for the
signicance analysis of parameters in the design of PM motors. They are the sizing
equation [3537], local sensitivity analysis (LSA) [38], global sensitivity analysis
[3942], and analysis of variance (ANOVA) [22, 27] techniques, respectively. The
last one, ANOVA, is based on the DOE technique. Two of them, LSA and DOE
will be introduced in the following sections.
Assume that f(x) is the objective function (such as output power, torque and cost) to
be optimized. Mathematically, the sensitivity of the ith parameter, xi, at the point x0
can be dened as
@f x
Si 4:15
@xi xx0
where Si is the sensitivity. The larger the jSi j, the more sensitive the objective
function f(x) is to the parameter xi [22, 38].
It should be noted that an analytical expression of objective function is required
in (4.15). However, there is no analytic form of objective function if motors
performance is calculated by using FEM. In this case, a differential form of (4.15)
should be used to calculate the sensitivity as the following
f x0 Dxi f x0
Si 4:16
Dxi
4.5 Sensitivity Analysis Techniques 133
where xi is the increment of parameter xi. In general, there are two methods to
determine this increment. The rst one is known as the parameter variation method,
in which xi is usually dened as 10, 20 % or both of its initial value. The other
method is known as the deviation variation method, in which xi is usually dened
as the standard deviation of xi . In this work, as the parameters deviation is not
given, the rst method will be used to calculate the sensitivity.
As an example, let us consider the output power as the objective and the
dimensions of PMs, such as width, x1, and height, x2, as the design parameters. The
sensitivity of PM width can be calculated as the following. Firstly, apply 10 %
perturbations to the PM width, and calculate the objective function (the output
power) corresponding to the two samples, i.e. (1 10 %)x1 and (1 + 10 %)x1.
Secondly, calculate the relative errors of these two samples by comparing the
objective function of these two samples to that obtained from the initial reference
point (x0). Finally, the average of absolute values of these two relative errors is
taken as the sensitivity of PM width on the output power.
It should be noted that different parameters have different units. To make the
obtained sensitivity values comparable, a normalization step below is needed
@f x=f x @f x 1 Df x 1
SSi 4:17
@xi =xi f x d f x d
where is the ratio of changed amplitude of parameter xi. By taking this normal-
ization, only the ratios of Dff x
x
are compared to acquire the sensitivity value for each
of the design parameters.
Basically, DOE is a kind of statistical method which has been widely used in the
design and data analysis of experiments in many areas, such as experiments in
agriculture, chemistry, and industrial design. The main aim of DOE is to arrange an
efcient experiment with smaller number of experiments, shorter experimental
cycle, and lower experimental cost, so as to obtain good experimental results and
scientic analysis conclusions. There are two types of DOE techniques, the full
factor design and partial factor design. The latter includes many further types, such
as orthogonal design and Latin hypercube design [22, 27].
ANOVA is a technique based on DOE, which can be used to determine the
signicant factors from all the design parameters. To implement the ANOVA, an
experiment table should be designed rstly by using DOE. The full-factor design
and orthogonal design are the two most popular DOE techniques. Because of the
high dimensional feature of motor design, it is time consuming to use the full-factor
method. For example, a motor design optimization problem of 8 parameters would
need 28 or 256 samples if a two-level full-factor design scheme is used. However,
134 4 Design Optimization Methods for Electrical Machines
only 12 samples are needed if the orthogonal design is used. Similarly, 38 or 6,561
samples are needed if a three-level full-factor design scheme is used, whereas if we
use the orthogonal design, only 27 points are needed. Therefore, orthogonal design
should be a good choice for most motor design optimization problems.
To implement the orthogonal design, the rst step is to select an orthogonal
design table from available tables. Table 4.11 shows an orthogonal design table of
L9(34), where subscript 9 indicates the number of experiments, and 3 the levels for
each parameter, while superscript 4 indicates that this table can be used for a
problem with no more than 4 design parameters and no interactions between them.
Then, the numbers 1, 2, and 3 in the table are the corresponding levels for each
parameter. Table 4.12 illustrates an orthogonal design table of L12(211), which can
be used for a problem with no more than 11 design parameters and no interactions
between them. Similarly, there are many available orthogonal design tables, such as
L12(241), L16(45), and L27(313).
After the DOE step, ANOVA can be used to determine the signicance of each
parameter. ANOVA is a family of multivariate statistical technique for helping infer
whether there are real differences between the means of three or more groups or
variables in a population, based on the sample data. In order to determine whether
the differences are signicant, ANOVA is concerned with differences between the
samples, known as the variance. By comparing the variance among sample mem-
bers, the differences are considered to be signicant if the variance is larger between
samples. Therefore, ANOVA can be regarded as a statistical test that looks for
signicant differences between means.
The understanding of ANOVA requires the background of multivariate statistics.
Fortunately, its implementation is very simple and can be realized by various
software packages, such as SPSS, Minitab, Matlab, and Excel.
This section presents an example for the sensitivity analysis of a PM claw pole
motor investigated in SOM. Six design parameters listed in Table 4.9 (Rsi, bs, hrm,
g1, hp and hsy) will be investigated for the sensitivity analysis of this motor by using
LSA. Table 4.13 lists the samples needed for the data analysis of LSA. For each
parameter, its initial value and four variation amplitudes, 20, 10, 10, and 20 %
will be considered. Totally, 25 samples are needed for the calculation, which
includes 24 points for those four variations (20, 10, 10, and 20 %) of six
parameters and 1 initial sample (0 column in the table).
Table 4.14 tabulates the analysis data obtained from the LSA technique. For the
sake of comparison, an average column, i.e. mean sensitivity, is listed in this table.
The sensitivity order can be obtained from the data in this column as
jbs j [ jhrm j [ jRsi j [ jg1 j [ hp [ hsy 4:18
To balance the optimization framework, we can take three of them as the sig-
nicant factors, which are Rsi, bs, and hrm. Actually, they are the parameters used
for the analysis of improved SOM in Sect. 4.3.5 [22].
It should be noted that two subspaces maybe reasonable for some problems. In
that case, combine X1 and X2 into one subspace, which may be known as the
signicant factor subspace. Meanwhile, the dimension of subspace is much smaller
than that in the initial space. Thus, the traditional direct optimization methods and
approximate models can be used in each of them.
138 4 Design Optimization Methods for Electrical Machines
Figure 4.3 illustrated a benchmark design example about the SMES. A three
parameter discrete optimization case has been investigated in Sect. 4.3.3. This
section will investigate an eight parameter continuous case. Table 4.15 lists the
scope of these eight parameters. From the previous discussion, it can be found that
those eight parameters can be divided into two subspaces in terms of their signi-
cances. One subspace involves the four parameters of inner coil {R1, h1, d1, J1}, and
the other one includes the other four parameters of outer coil {R2, h2, d2, J2} [33].
Table 4.16 shows obtained optimization solutions for this continuous case of
SMES. As shown, 4720 FEM samples are needed for DEA to get the optimal
solution under the direct optimization framework. The optimal stored energy in
SMES is 178.75 MJ (the error is 1.25 MJ), and mean stray eld is 2.27 mT.
By the multi-level optimization method, only 1078 FEM samples are needed for
the optimization, which is less than 1/4 of that by DEA. The resultant optimal mean
stray eld is 3.23 mT, slightly higher than that by DEA, and the error of energy is
1.01 MJ which is smaller than that given by DEA. Therefore, the proposed
multi-level optimization method is more efcient than DEA for high dimensional
design problems of electromagnetic devices. It will be employed to optimize
electrical machines and drive systems in the following two chapters.
where IP, PS, EL and FIT represent the initial population, population size,
encoding length and tness value, respectively, SO, CO and MO are the genetic
operations, namely, selection, crossover and mutation operations.
The MLGA can be described as follows.
where GAij stands for applying the independent GA to the ith level and the jth
module. With the reaction between different levels and adjacent sub-modules on the
same level, GAij can be described as the following:
GAij IPi j GAi ;j1 ; GAi1 ;j ; GAi ;j 1
PSi j GAi ;j1 ; GAi1 ;j ; GAi ;j 1
ELi j GAi ;j1 ; GAi1 ;j ; GAi ;j 1
FITi j GAi ;j1 ; GAi1 ;j ; GAi ;j 1 4:21
SOi j GAi ;j1 ; GAi1 ;j ; GAi ;j 1
COi j GAi ;j1 ; GAi1 ;j ; GAi ;j 1
MOi j GAi ;j1 ; GAi1 ;j ; GAi ;j 1
4.7 Multi-level Genetic Algorithm 141
Step 3: Determine the architecture of MLGA, including the number of levels and
the number of modules in each level.
Step 4: Allocate the design variables, objective functions and constraints on dif-
ferent levels according to the problem matrix, and build up the relation-
ships among different levels and different modules on each level. Each
module corresponds to a genetic algorithm module.
Step 5: Implement the MLGA modules of each level from the top to the lowest
level. The upper level module sends control messages and parameter
values to the lower level module. Feedback messages from the lower level
are used as the evaluation function by the upper level.
Step 6: The total solving process ends when the termination criterion of the top
level has been reached. Otherwise, Step 5 will be repeated.
max f x w 100g
K
1 100 w2 Costmax
Cost
where the design variable x hm; bm; Ns; WindD, Ns and WindD are discrete
variables, is the efciency of the SPMSM, P2 the output power, sf the ll factor,
and K, w1 and w2 are the weighting factors specied by the designer [44].
4.7 Multi-level Genetic Algorithm 143
(a)
(b)
(b)
parameters when the magnet thickness and width are changed. Table 4.18 lists the
main motor parameters obtained from the experimental results, and the calculated
results based on FEM. Figure 4.28 illustrates the magnetic eld distribution when
Lad and Laq are calculated. Figure 4.29 shows the bi-level architecture of opti-
mization for SPMSM [44].
146 4 Design Optimization Methods for Electrical Machines
Figure 4.31 shows a classic design framework and the coupled relations for elec-
trical machines. As shown, the design optimization is really a multi-domain
problem which includes electromagnetic, material, mechanical, and thermal aspects,
and they can be strongly coupled. In order to achieve high performance, the
multi-disciplinary design optimization (MDO) of electrical machines must be
investigated. It includes the following ve main steps [47]:
Step 1: Denition of the motor specications. It mainly consists of cost, such as
material cost and manufacturing cost, output performance, such as power,
torque, efciency and speed, and other constraints, such as volume,
weight, temperature rise, mechanical strength and resonance frequency,
etc.
Step 2: Selection of motor type and its topology. For example, for the motor types
of the PM-SMC motor investigated above, there are several options, such
as TFM, claw pole and flux switching motors. The topology options may
include outer rotor, inner rotor and numbers of poles for SMC motors.
Step 3: Initial design. To acquire a possible design scheme, three main
designs/selections in terms of dimensions, materials and manufacturing
methods are required to investigate in this step. For example, for the
materials of PM-SMC motors, PM, steel, SMC and ferrite can be the
options. For the manufacturing method, the press method (the moulding
method) is recommended based on our design and prototyping experi-
ences as it is widely used for batch production of PM-SMC motors.
148 4 Design Optimization Methods for Electrical Machines
machine has been investigated in Sect. 2.4, the electromagnetic and thermal anal-
yses will be the main contents in this section. Meanwhile, a lumped 3D thermal
network model is developed for the thermal analysis. Then, a multi-disciplinary
optimization model is proposed to minimize the material cost and maximize the
output power based on the proposed thermal network model. Finally, the FEM is
employed to verify the performances obtained from optimal results in terms of
thermal analysis and modal analysis.
1.5
B (T)
1
7.32 g/cm3
Cost ($/hour)
300
200
100
100 200 300 400 500
(b) Press size (ton)
500
Productivity (pieces/hour)
400
300
200
100
100 200 300 400 500
Press size (ton)
Thermal analysis is used to calculate the temperature rises in winding and PM rotor
for this machine, so as to ensure that the motor works safely [47, 48]. For this
machine, a 2D thermal network model as shown in Fig. 2.8 was developed in our
previous work to simulate the thermal analysis. Considering the thermal isotropy of
SMC material, we developed a 3D thermal network of lumped parameters as it can
provide more accurate results than the 2D network model.
4.8 Multi-disciplinary Optimization Method 151
Figure 4.34 describes the 3D thermal network model for the middle stack of the
studied PM TFM. The two side stacks can be neglected as the temperature of the
middle stack is higher than those of two side stacks. Therefore, only the middle
stack is investigated in this model. Meanwhile, the major heat dissipates from the
rotor. The resistances to the thermal conduction of the following sections are cal-
culated: two segments of stator yoke (Rsy1and Rsy2), stator side disk (Rsd1 and Rsd2),
coils (Rcu1, Rcu2, and Rcu3), stator teeth (Rst1), air gaps (Rg1, Rg2, Rg3, Rg4 and Rg5),
PMs (Rpm1, Rpm2), rotor in radial direction (Rrt1, Rrt2), rotor in axial direction
(Rry1, Rry2), and shaft (Rsf1 and Rsf2). The equivalent resistances are calculated for
152 4 Design Optimization Methods for Electrical Machines
thermal convection between the air and stator teeth (Rsta), air and stator disk (Rsda),
air and coil (Rcua), air and PM (Rpm1, Rpm2), air and rotor (Rrta), rotor and outer
ambient air (Rrya), and shaft and air (Rsfa), respectively. The thermal resistances in
the circumferential direction in this motor are calculated for thermal conduction in
the stator yoke (Rsya1 and Rsya2), coil (Rcu3), stator teeth (Rsta1), and rotor (Rra1 and
Rra2).
The heat sources in this model include the stator and rotor core losses, copper
loss, and mechanical loss. In order to gain a relatively high accuracy, each loss is
divided into several parts. The stator core loss is divided into six parts (PFes1), the
copper loss two parts (Pcu1), the rotor core loss four parts (PFer1), and mechanical
loss six parts (Pmech1), respectively.
Based on this 3D thermal network model, it is found that there is 68 C tem-
perature rise in the coil and 27 C in the rotor yoke surface for this PM-SMC TFM
prototype. Comparison with experimental results will be shown in Sect. 4.8.5.
min: C0 Pout
f x Cost P0
s.t: g1 x 0:795 g 0;
g2 x 640 Pout 0; 4:23
g3 x sf 0:7 0;
g4 x TPM 65 0;
g5 x TCoil 65 0;
where x is a vector of design parameters, C0 and P0 are the cost and output power of
the initial design scheme, and Pout in g1 and g2 the motors efciency and output
power, respectively, sf in g3 is the ll factor, and TPM and TCoil in g4 and g5 are the
temperature rises in the PM and windings, respectively. From our design experi-
ence, six parameters as shown in Table 4.20 are signicant to the performance of
this machine. The cost in the objective function mainly includes the material and
manufacturing costs of the SMC core.
Modal analysis is not included in this optimization model. However, to ensure
that the optimized motor has good mechanical performance, modal analysis will be
presented in the next section to verify the performance of the optimized motor. For
the thermal analysis in optimization, the lump 3D thermal network model is used to
replace the FEM analysis to improve the optimization efciency. Then, a FEM
method will be presented for the thermal analysis of the nal optimal scheme in the
next section.
4.8 Multi-disciplinary Optimization Method 153
Fig. 4.36 Distribution of temperature rise in the coil for MDO optimum
Fig. 4.37 Illustration of rst order modal analysis for MDO optimum
156 4 Design Optimization Methods for Electrical Machines
4.9 Summary
References
1. Lei G, Wang TS, Guo YG, Zhu JG, Wang SH (2014) System level design optimization
method for electrical drive system: deterministic approach. IEEE Trans Ind Electron 61
(12):65916602
2. Vese I, Marignetti F, Radulescu MM (2010) Multiphysics approach to numerical modeling of
a permanent-magnet tubular linear motor. IEEE Trans Ind Electron 57(1):320326
3. Kreuawan S, Gillon F, Brochet P (2008) Optimal design of permanent magnet motor using
multi-disciplinary design optimisation, In: proceedings of 18th international conference on
electrical machines, Vilamoura. pp 16, 69 Sep 2008
4. Yao D, Ionel DM (2013) A review of recent developments in electrical machine design
optimization methods with a permanent-magnet synchronous motor benchmark study. IEEE
Trans Ind Appl 49(3):12681275
5. Di Barba P (2010) Multi-objective shape design in electricity and magnetism, Lecture Notes
Elect Eng, vol 47
6. Reyes-Sierra M, Coello CAC (2006) Multi-objective particle swarm optimizers: A survey of
the state-of-the-art. Int J Computat Intell Res 2(3):287308
7. Yamazaki K, Ishigami H (2010) Rotor-shape optimization of interior-permanent-magnet
motors to reduce harmonic iron losses. IEEE Trans Ind Electron 57(1):6169
8. Barcaro M, Bianchi N, Magnussen F (2012) Permanent-magnet optimization in
permanent-magnet-assisted synchronous reluctance motor for a wide constant-power speed
range. IEEE Trans Ind Electron 59(6):24952502
9. Komeza K, Dems M (2012) Finite-element and analytical calculations of no-load core losses
in energy-saving induction motors. IEEE Trans Ind Electron 59(7):29342946
10. Zhao W L, Lipo T A, Kwon B-I (2015) Optimal design of a novel asymmetrical rotor structure
to obtain torque and efciency improvement in surface inset PM motors, IEEE Trans Magn, 51
(3), Article no 8100704
11. Kim H-W, Kim K-T, Jo Y-S, Hur J (2013) Optimization methods of torque density for
developing the neodymium free SPOKE-type BLDC motor. IEEE Trans Magn 49(5):2173
2176
12. Lei G, Shao KR, Guo Y, Zhu J, Lavers JD (2008) Sequential optimization method for the
design of electromagnetic device. IEEE Trans Magn 44(11):32173220
13. Lei G, Zhu JG, Guo YG, Zou Y (2014) State of art of sequential optimization strategies for the
design of electromagnetic devices. In Proceedings of the 17th international conference on
electrical machines and systems (ICEMS), pp 706709, Oct. 2225, 2014
158 4 Design Optimization Methods for Electrical Machines
14. Storn R, Price K (1997) Differential evolution - a simple and efcient heuristic for global
optimization over continuous spaces. J Global Optim 11:341359
15. Takahashi RHC, Vasconcelos JA, Ramrez JA, Krahenbuhl L (2003) A multi-objective
methodology for evaluating genetic operators. IEEE Trans Magn 39(3):13211324
16. Campelo F, Guimaraes FG, Igarashi H, Ramirez JA (2005) A clonal selection algorithm for
optimization in electromagnetics. IEEE Trans Magn 41(5):17361739
17. Alotto P, Baumgartner U, Freschi F, Kstinger A, Magele Ch, Renhart W, Repetto M (2008)
SMES optimization benchmark extended: introducing Pareto optimal solutions into TEAM22.
IEEE Trans Magn 44(6):10661069
18. Wang LD, Lowther DA (2006) Selection of approximation models for electromagnetic device
optimization. IEEE Trans Magn 42(2):12271230
19. Lei G, Shao KR, Guo Y, Zhu J, Lavers JD (2009) Improved SOM for high dimensional
electromagnetic optimization problems. IEEE Trans Magn 45(10):39933996
20. Lei G, Yang GY, Shao KR, Guo YG, Zhu JG, Lavers JD (2010) Electromagnetic device
design based on RBF models and two new sequential optimization strategies. IEEE Trans
Magn 46(8):31813184
21. Guo YG, Zhu JG, Dorrell D (2009) Design and analysis of a claw pole PM motor with molded
SMC core. IEEE Trans Magn 45(10):5824585
22. Lei G, Liu CC, Zhu JG, Guo YG (2015) Techniques for multi-level design optimization of
permanent magnet motors. IEEE Trans Energy Conver 30(4):15741584
23. Lei G, Shao KR, Guo YG, Zhu JG (2012) Multi-objective sequential optimization method for
the design of industrial electromagnetic devices. IEEE Trans Magn 48(11):45384541
24. Lebensztajn L, Coulomb JL (2004) TEAM workshop problem 25: A multi-objective analysis.
IEEE Trans Magn 40(2):14021405
25. Xie DX, Sun XW, Bai BD, Yang SY (2008) Multi-objective optimization based on response
surface model and its application to engineering shape design. IEEE Trans Magn 44(6):1006
1009
26. Di Barba P (2009) Evolutionary multi-objective optimization methods for the shape design of
industrial electromagnetic devices. IEEE Trans Magn 45(3):14361441
27. Wu CFJ, Hamada MS (2000) Experiments: Planning, Analysis Parameter Design
Optimization. Wiley, New York
28. Deb K, Pratap A, Agarwal S, Meyarivan T (2002) A fast and elitist multi-objective genetic
algorithm: NSGA-II. IEEE Trans Evol Comput 6(2):182197
29. Zhu JG, Guo YG, Lin ZW, Li YJ, Huang YK (2011) Development of PM transverse flux
motors with soft magnetic composite cores. IEEE Trans Magn 47(10):43764383
30. Guo YG, Zhu JG, Watterson PA, Wei Wu (2006) Development of a PM transverse flux motor
with soft magnetic composite core. IEEE Trans Energy Conver 21(2):426434
31. Lei G, Guo YG, Zhu JG et al (2012) System level six sigma robust optimization of a drive
system with PM transverse flux machine. IEEE Trans Magn 48(2):923926
32. Lei G, Zhu JG, Guo YG, Hu JF, Xu W, Shao KR (2013) Robust design optimization of
PM-SMC motors for Six Sigma quality manufacturing. IEEE Trans Magn 49(7):39533956
33. Lei G, Guo YG, Zhu JG, Chen XM, Xu W (2012) Sequential subspace optimization method
for electromagnetic devices design with orthogonal design technique. IEEE Trans Magn 48
(2):479482
34. Lei G, Xu W, Hu JF, Zhu JG, Guo YG, Shao KR (2014) Multi-level design optimization of a
FSPMM drive system by using sequential subspace optimization method, IEEE Trans Magn,
50(2), Article no 7016904
35. Liu CC, Zhu JG, Wang YH, Lei G, Guo YG, Liu XJ (2014) A low-cost permanent magnet
synchronous motor with SMC and ferrite PM, In: Proceedings of 17th international conference
on electrical Machines and Systems (ICEMS), pp 397400
36. Fei W, Luk PCK, Shen JX, Wang Y, Jin M (2012) A novel permanent-magnet flux switching
machine with an outer-rotor conguration for in-wheel light traction applications. IEEE Trans
Ind Appl 48(5):14961506
References 159
37. Hua W, Cheng M, Zhu ZQ, Howe D (2006) Design of flux-switching permanent magnet
machine considering the limitation of inverter and flux-weakening capability. In: proceedings
of 41st IAS annual meeting-industry applications conference, vol 5, pp 24032410
38. Morio J (2011) Global and local sensitivity analysis methods for a physical system. Eur J Phys
32(6):15771583
39. Rodriguez-Fernandez M, Banga JR, Doyle FJ (2012) Novel global sensitivity analysis
methodology accounting for the crucial role of the distribution of input parameters: application
to systems biology models. Int J Robust Nonlin Control 22:10821102
40. Morris MD (1991) Factorial sampling plans for preliminary computational experiments.
Technometrics 33(2):161174
41. Sobol IM (2011) Global sensitivity indices for nonlinear mathematical models and their Monte
Carlo estimates. Math Comput Simul 55:271280
42. Herman JD, Kollat JB, Reed PM, Wagener T (2013) Technical note: method of morris
effectively reduces the computational demands of global sensitivity analysis for distributed
watershed models. Hydrol Earth Syst Sci 17:28932903
43. Li QS, Liu DK, Leung AYT, Zhang N, Luo QZ (2002) A multi-level genetic algorithm for the
optimum design of structural control systems. Int J Numer Meth. Engng 55:817834
44. Wang SH, Meng XJ, Guo NN, Li HB, Qiu J, Zhu JG et al (2009) Multi-level optimization for
surface mounted PM machine incorporating with FEM. IEEE Trans Magn 45(10):47004703
45. Meng XJ, Wang SH, Qiu J, Zhu JG, Wang Y, Guo YG et al (2010) Dynamic multi-level
optimization of machine design and control parameters based on correlation analysis. IEEE
Trans Magn 46(8):27792782
46. Meng XJ, Wang SH, Qiu J, Zhang QH, Zhu JG, Guo YG, Liu DK (2011) Robust multi-level
optimization of PMSM using design for six sigma. IEEE Trans Magn 47(10):32483251
47. Lei G,Liu CC, Guo YG, Zhu JG (2015) Multi-disciplinary design analysis for PM motors with
soft magnetic composite cores, IEEE Trans Magn, 51(11), Article no 8109704
48. Huang YK, Zhu JG et al (2009) Thermal analysis of high-speed SMC motor based on thermal
network and 3D FEA with rotational core loss included. IEEE Trans Magn 45(106):4680
4683
Chapter 5
Design Optimization Methods
for Electrical Drive Systems
Keywords Electrical drive systems System-level design optimization Multi-level
design optimization Field oriented control Finite element methods Model pre-
dictive control Transverse flux machine Permanent magnet synchronous machine
5.1 Introduction
Electrical machines and the corresponding drive systems have a history of over a
century and the design procedure has become almost standard. When designing
an appliance that needs an electrical drive system, the designer rstly selects the
motor, inverter/converter and controller from the existing products. The appliance
designer, on one hand, has to deliver the functions that the appliance is supposed to
have, and on the other hand, has to take into account the availability and perfor-
mance that the existing motor drive can provide. This motor manufacturer-oriented
approach has been the dominant design concept for drive systems for a long time.
However, this approach would apply many constraints to the design and therefore
limit the functions of the appliance [1].
With the fast development of CAD/CAE software, new material, flexible
mechanical manufacturing technology, advanced optimization and control algo-
rithms, it is possible to design a motor to meet the special requirements of a
particular application. Since early 1990s, this application oriented approach has
become a common practice. Nowadays, the motors and their control systems are
generally closely integrated into the appliances. Therefore, more and more holistic
integrated design problems of the electrical drive systems have boomed in industry,
for example, the drive systems for hybrid electric vehicles (HEVs) [24].
Through the extensive research practice, it is recognized that when designing
such an electrical drive system, it is important to pursue the optimal system per-
formance rather than the optimal components like motor, because assembling
individually optimized components into a system cannot necessarily guarantee an
optimal system performance. The optimal system performance can only be
achieved through a holistic approach of integrated simultaneous optimization of all
components at the system level [1, 5, 6].
Figure 2.1 illustrated a general design framework and the interactions between
different disciplines/domains for electrical drive systems. As shown, the design
optimization of such a drive system is really a multi-disciplinary problem. It mainly
includes electromagnetic, material, mechanical, thermal and power electronic
designs, which are strongly coupled [79]. In order to achieve high system-level
performance, the perfect cooperation of motor and its drive and control systems
must be designed and optimized simultaneously.
Although the importance of system-level design optimization of electrical drive
systems is noted, not much work has been reported in the literature [1, 5, 6].
Traditional design and optimization methods are mostly on the component level of
different kinds of motors [1016]. Generally, cogging torque, torque ripple, cost,
weight and energy consumption are the main concerns for motors performance
parameters in the design and optimization process [1719]. For the design opti-
mization of these motors, Chap. 4 presented several kinds of optimization methods,
including the combinations of intelligent optimization algorithms and nite element
model (FEM) or approximate models, for example, the differential evolution
algorithm (DEA) plus Kriging model. Approximate models are generally used to
replace the FEM in the performance evaluation of motors so as to reduce the FEM
computational costs [1, 5, 10, 20, 21].
On the other hand, for the controller part, though a lot of control algorithms have
been developed, such as eld oriented control (FOC), direct torque control
(DTC) and model predictive control (MPC) [2226], they are also generally
designed and optimized on the controller level, and have not been combined with
the design optimization of motors [27].
This component-level-based method may be reasonable for some traditional
motors and their drive systems where there is much design experience that can be
used. However, there is not much design experience for novel electrical drive
systems. Furthermore, these methods are not system-level holistic design basically.
5.1 Introduction 163
As previously discussed, by this component level approach, one can hardly achieve
the optimal system performance. Therefore, how to design and optimize novel high
performance drive systems is an important problem in both research community
and industrial applications [1, 5].
In order to deal with the above problems, this chapter presents three types of
system-level design optimization methods for electrical drive systems. This chapter
is structured as follows. Section 5.2 presents the system-level design optimization
framework and models for electrical drive systems. Section 5.3 presents a
single-level (only at the system level) optimization method for the design of drive
systems. Section 5.4 presents a multi-level design optimization method for drive
systems, including the investigation of the rst design example: a drive system
consists of a permanent magnet (PM) transverse flux machine (TFM) with soft
magnetic composite (SMC) core and an improved MPC system. Section 5.5
introduces the multi-level Genetic Algorithm (MLGA) for the design optimization
of drive systems, including the second drive system example which is composed of
a surface-mounted permanent magnet synchronous machine (SPMSM) and a
classical FOC system, followed by the summary section.
min : fm xm
s:t: gmi xm 0; i 1; . . .; Nm ; 5:1
xml xm xmu
where xm, fm and gm are the motor design parameter vector, objectives and con-
straints, xml and xmu the lower and upper boundaries of xm, respectively, and Nm is
the number of the constraints. It should be noted that the objectives and constraints
5.2 System-Level Design Optimization Framework 165
in (5.1) must be dened in terms of the required system objectives and constraints
in Step 1.
The control design optimization model can be dened as
min : fc xc
s:t: gci xc 0; i 1; . . .; Nc ; 5:2
xcl xc xcu
where xc, fc and gc are the control design parameter vector, objectives and con-
straints, xcl and xcu the lower and upper boundaries of xc, respectively, and Nc is the
number of controller constraints. Similarly, the objectives and constraints in (5.2)
must also be dened in terms of the required system objectives and constraints.
Combining the motor and controller design optimization models, (5.1) and (5.2),
one obtains the system-level design optimization model as the following:
min : fs xs Ffm ; fc
s:t: gmi xs 0; i 1; 2; . . .; Nm
; 5:3
gci xs 0; i 1; 2; . . .; Nc
xsl xs xsu
where xs xm ; xc , xsl and xsu are the lower and upper boundaries of xs, respec-
tively, and fs is the system objective which is generally a function of fm and fc.
Step 5 Evaluation of the system performance. This step consists of two parts. One
is the evaluation of steady performance of motor, such as cost and ef-
ciency. The other is the evaluation of dynamic performance of controller or
the whole drive system, such as overshoot, settling time, torque ripple, and
speed ripple.
Figure 5.2 illustrates the rst type of optimization method for electrical drive
systems. It can be seen that the optimization process is implemented at a single level
for the whole system, which is thus known as the single-level design optimization
method [1].
This method mainly includes the following three steps:
Step 1 Determination of system level optimization model (5.3). It includes the
selection of motor and controller for the specic drive system.
Step 2 Selection of an optimization method. As drive systems are always high
dimensional and non-linear design problems, intelligent algorithms, such
as genetic algorithm (GA) and DEA, can be good choices in many situ-
ations. Therefore, the algorithm parameters should be determined in this
step, such as genetic operators in GA and mutation operator in DEA.
166 5 Design Optimization Methods for Electrical Drive Systems
Figure 5.3 depicts a multi-level design optimization framework for electrical drive
systems. Three levels are considered in this framework, namely the motor, control,
and system levels [1].
5.4 Multi-level Design Optimization Method 167
The control levelThe aim of this level is to optimize the control model (5.2) and
evaluate the system dynamic performance, such as the overshoot and settling time.
Step 3 Verication of system level performance (5.3). The aim of this step is to
evaluate the system performance and output the optimization results.
As outlined above, the drive system design optimization is always high
dimensional and nonlinear. The optimization efciency of multi-level optimization
method mainly depends on two issues. The rst one is how to construct an efcient
multi-level optimization framework, especially for high dimensional problems. The
second one is how to reduce the computational cost of optimization models, which
are the main contents of the previous chapter. Therefore, the proposed new opti-
mization methods for electrical machines can be employed for the optimization of
motor and control levels, respectively. For example, if there are seven optimization
parameters for the motor, then the multi-level optimization method presented in
Sects. 4.6 and 4.7 can be introduced to improve the optimization efciency, and the
sensitivity analysis techniques presented in Sect. 4.5 can be employed to determine
the multi-level optimization framework.
min : fm xm w1 Cost
C0 w2 Pout
P0
Fig. 5.4 Structure of the PM-SMC TFM with SMC stator (one phase)
where w1 and w2 are weighting factors, C0 and P0 the cost and output power (Pout)
of the initial PM TFM, , sf and Jc the motor efciency, winding ll factor, and
current density, respectively [1, 5, 28, 29].
In Chap. 2, an improved MPC scheme with a duty ratio optimization module was
presented to drive PM motors [30]. It will be used in this chapter as the control
method for the PM-SMC TFM. This control scheme was shown in Fig. 2.37.
170 5 Design Optimization Methods for Electrical Drive Systems
As mentioned in Chap. 2, the key issue of MPC is the denition of cost function.
Since the two greatest concerns of a PM motor are the torque and stator flux, the
cost function is dened to ensure that both the torque and stator flux at the end of
control period are as close as possible to the reference values. To illustrate the
optimization parameters of the control level, the cost function is rewritten as
G jTe Tek 1 j k1 jws j jwks 1 j
5:5
A jT T k N j k1 jw j jwk N j
e e s s
where Te , ws , Tek 1 and wsk 1 are the reference torque and flux, predicted torque
and flux, Tek N and wks N the linear predictions of torque and flux at the (k + N)-th
instant, and k1 and A weighting factors, respectively.
The other important part of this improved MPC is the duty ratio module. The
expression of duty ratio optimization module has the form as
Te Tek 1 ws wek 1
d ;
C C 5:6
T w
where CT and C are two positive parameters. The idea of this method is that a
larger difference between the reference and predicted torque values would lead to a
larger duty ratio value [30].
Six parameters should be optimized in the control level. They are A, N, CT, C,
Kp and Ki, where Kp and Ki are the PI controller parameters as shown in Fig. 2.37.
One objective and four constraints are considered for this level. The objective is to
minimize the sum of root mean square errors (RMSE) of torque (T) and speed (n) in
the steady state operation period. At the same time, the speed overshoot should be
minimized for this control system. The optimization model of the control system
can be expressed as
5.4 Multi-level Design Optimization Method 171
min : fc xc w3 RMSE
Trated
T
w4 RMSE
nrated
n
w5 nos
s:t: gc1 xc RMSET=Trated 0:06 0;
gc2 xc RMSEn=nrated 0:05 0; 5:7
gc3 xc nos 0:02 0;
gc4 xc ts 0:02 0;
xcl xc xcu
where w3 to w5 are weighting factors, the subscript rated indicates that the values
are obtained from the motor optimization model (5.4), nos is the speed overshoot,
which should be no larger than 2 % of the rated speed, 1800 rev/min, and ts the
settling time, which should be no larger than 0.02 s after the load is applied to the
control system [1].
However, since the signicant space only has 1 factor, we just select all six
parameters as the same level.
In summary, the total optimization framework of this drive system has three
levels as shown in Fig. 5.6. The rst and second levels are the subspaces X1 and
X2, respectively, which should be optimized with model (5.4). After the opti-
mization of motor level, the motor characteristic parameters, such as R, L and flux
can be obtained. They will be used as the input parameters of MPC control system
in the control level. The third level is the subspace of all the control parameters in
model (5.7).
B.Optimization Results
First of all, DEA is selected as the optimization algorithm in the multi-level opti-
mization of this drive system. The algorithm parameters include the mutation
scaling factor of 0.8, crossover factor of 0.8, and the maximum number of iterations
of 1000. Then in the multi-level optimization method is dened as 1 %. All
weighting factors are assumed to be 1 in this work. Tables 5.3 and 5.4 list the
optimization results of motor and control levels obtained by the single-level and
multi-level optimization methods, respectively. From these tables, the following
conclusions can be drawn:
(1) The motor level. For the initial design scheme, the motor efciency is 79.5 %,
the output power 640 W, average torque 3.40 Nm, and material cost $35.8.
By the single-level optimization method with DEA and FEM, all the 14
parameters (8 motor and 6 control parameters) are optimized as shown in
Fig. 5.2. The obtained motor efciency is 81.5 %, output power 658 W,
average torque 3.49 Nm, and material cost $28.3. They are better than those of
the initial design.
For the multi-level optimization method, 3 iteration processes are required to
get the optimal results. Figure 5.7 shows the iteration process of multi-level
optimization for the motor level. As shown, level 1 is optimized twice while
level 2 is optimized only once. After the optimization, the optimal output
power reaches 670 W and the average torque 3.55 Nm while the motor ef-
ciency decreases to 81.3 %. The cost is only $26.9, which is the smallest
among these three designs, and the output power is increased by 4.7 %
(30/640) by using the proposed multi-level optimization method compared
with the initial design.
(2) The control level. As shown in Table 5.4, the relative RMSEs of torque and
speed are 4.17 and 0.10 %, respectively, the speed overshoot is 1.03 % and the
objective 5.30 % for control level by using the single-level optimization
method. After the multi-level optimization, all these objectives have been
increased signicantly. For example, the relative RMSE of speed has been
5.4 Multi-level Design Optimization Method 173
decreased from 0.10 to 0.02 %, and the objective of control level from 5.30 to
4.87 %. Therefore, the dynamic performances have been greatly improved by
using the multi-level optimization method. Figure 5.8 illustrates the dynamic
performance of the drive system obtained by the multi-level method. As
shown, the dynamic performances of speed and torque are very good.
(3) For the computational cost, the cost of FEM analysis at the motor level and the
cost of Simulink simulation calls at control level are the largest computational
burden for the whole optimization process. For the single-level optimization
method with DEA, about 14,000 FEM samples and 14,000 Simulink simu-
lation calls (5 14 200, where 5 14 is the population size and 200 the
5.4 Multi-level Design Optimization Method 175
1.6
1.4
1.2
1
Initial Level 1 - 1st Level 2 - 1st Level 1 - 2nd
Iteration number of multilevel optimization
Fig. 5.7 Iteration process of multi-level optimization for the motor level
6
Torque [N*m]
-2
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14
Time [s]
2000
Speed [r/min]
1500
1000
500
0
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14
Time [s]
average iteration number of DEA) are needed to achieve the optimal results.
Actually, it is hard for an intelligent algorithm to deal with this kind of
high-dimensional and highly nonlinear optimization problem. It is time con-
suming and tends to nd a local optimal point.
On the other hand, only 1625 (see Table 5.3) FEM samples are needed for the
motor level optimization by using the proposed multi-level optimization method,
which is about 11.6 % of the direct single-level optimization method. Moreover,
about 6,000 Simulink simulation calls are needed for the control level by using the
multi-level optimization method. This is less than half of the simulation cost of
single level optimization method. Therefore, the proposed multi-level design
optimization method can signicantly reduce the computational cost, and produce
better solutions than the single-level optimization method.
In summary, compared with the schemes gained from the initial design and
single-level method, the solutions obtained by the multi-level optimization method
have many improvements, such as larger output power, less material cost, and less
overshoot. As a matter of fact, both the steady state and dynamic performances of
this drive system have been improved by using the multi-level method.
MLGA was introduced as a kind of optimization method for the high dimensional
optimization problems of electrical machines in Chap. 4. Its efciency has been
veried by a design example of an SPMSM [31, 32]. MLGA can be applied to
design optimization of electrical drive systems with a similar optimization frame-
work. In this section, MLGA is presented for design optimization of a motor drive
system consisting of a SPMSM and FOC control scheme. The main design
parameters of this SPMSM can be seen in Sect. 4.7, and the FOC control scheme in
Fig. 2.22. As shown, three PI controllers are used for the d- and q-axis components
of stator current, and speed control, respectively.
For the motor level optimization, the objective is to minimize the cost of copper
and permanent magnets, and to maximize the motor efciency, . The optimization
model can be expressed as
K
max : fm xm
w1 100g w2 Cost
Cost
100 max 5:8
s.t. P2 [ 945W
sf \78 %
5.5 MLGA for a SPMSM Drive System with FOC 177
where the design variable xm hm; bm; Ns; WindD, Ns and WindD are discrete
variables, K, w1 and w2 the weighting factors dened by the designer, P2 is the
output power, and sf the ll factor [6].
For the control level optimization, the six integral and proportional gain factors
in the three PI controllers as shown in Fig. 2.22 are chosen as the design
optimization variables. The optimization model can be dened as
where xc stands for the vector of six PI variables, Trip is the torque ripple, nos the
overshoot of speed, Id the d-axis component of stator current, and i (i = 1, 2, 3) are
weighting factors.
In the implementation of MLGA, all weighting factors are dened as 1 in this work.
Table 5.5 lists the optimization results obtained by MLGA, and Table 5.6 the
proportional and integral gains calculated on the third level. It can be seen that the
motor efciency and output power have been increased greatly by using MLGA
compared with the initial design. For example, the motor efciency after MLGA
optimization reaches 86.4 %, an increase by 2.7 % compared with the initial
83.7 %. Figure 5.10 illustrates the speed responses of SPMSM before and after the
178 5 Design Optimization Methods for Electrical Drive Systems
5.6 Summary
From the above discussions, it can be seen that the system-level design optimization
method is necessary for electrical drive systems in order to achieve high steady state
and dynamic performances at the system level. The proposed multi-level method is
180 5 Design Optimization Methods for Electrical Drive Systems
efcient for the design optimization of high dimensional drive systems. This
method will build a solid foundation to enable the effective development of novel
high performance drive systems with new materials, low cost and high efciency
for industrial applications. This method can be also applied to other high dimen-
sional design optimization problems in industrial applications. It will shorten the
design cycle, reduce the design cost and improve the design efciency for the
industrial products in the early stage of product development.
References
1. Lei G, Wang TS, Guo YG, Zhu JG, Wang SH (2014) System level design optimization
method for electrical drive system: deterministic approach. IEEE Trans Ind Electron 61
(12):65916602
2. Zhu ZQ, Howe D (2007) Electrical machines and drives for electric, hybrid, and fuel cell
vehicles. Proc IEEE 95(4):746765
3. Chan CC (2007) The state of the art of electric, hybrid, and fuel cell vehicles. Proc IEEE 95
(4):704718
4. Emadi A, Lee YJ, Rajashekara K (2008) Power electronics and motor drives in electric, hybrid
electric, and plug-in hybrid electric vehicles. IEEE Trans Ind Electron 55(6):22372245
5. Lei G, Guo YG, Zhu JG et al (2012) System level six sigma robust optimization of a drive
system with PM transverse flux machine. IEEE Trans Magn 48(2):923926
6. Meng XJ, Wang SH, Qiu J, Zhu JG, Wang Y, Guo YG et al (2010) Dynamic multi-level
optimization of machine design and control parameters based on correlation analysis. IEEE
Trans Magn 46(8):27792782
7. Lei G, Liu CC, Guo YG, Zhu JG (2015) Multi-disciplinary design analysis and optimization of
a PM transverse flux machine with soft magnetic composite core, IEEE Trans Magn, 2015,
vol 51(11). Article 8109704
8. Vese I, Marignetti F, Radulescu MM (2010) Multiphysics approach to numerical modeling of
a permanent-magnet tubular linear motor. IEEE Trans Ind Electron 57(1):320326
9. Kreuawan S, Gillon F, Brochet P (2008) Optimal design of permanent magnet motor using
multi-disciplinary design optimisation. In: Proceedings of the 18th International Conference
on Electrical Machines, 69 Sept., Vilamoura, pp 16
10. Lei G, Shao KR, Guo YG, Zhu JG (2012) Multi-objective sequential optimization method for
the design of industrial electromagnetic devices. IEEE Trans Magn 48(11):45384541
11. Li WL, Zhang XC, Cheng SK, Cao JC (2013) Thermal optimization for a HSPMG used for
distributed generation systems. IEEE Trans Ind Electron 60(2):474482
12. Pster P-D, Perriard Y (2010) Very-high-speed slotless permanent-magnet motors: analytical
modeling, optimization, design, and torque measurement methods. IEEE Trans Ind Electron
57(1):296303
13. Laskaris KI, Kladas AG (2011) Permanent-magnet shape optimization effects on synchronous
motor performance. IEEE Trans Ind Electron 58(9):37763783
14. Yamazaki K, Ishigami H (2010) Rotor-shape optimization of interior-permanent-magnet
motors to reduce harmonic iron losses. IEEE Trans Ind Electron 57(1):6169
15. Ma C, Qiu L (2015) Multi-objective optimization of switched reluctance motors based on design
of experiments and particle swarm optimization. IEEE Trans Energy Convers 30(3):11441153
16. Barcaro M, Bianchi N, Magnussen F (2012) Permanent-magnet optimization in
permanent-magnet-assisted synchronous reluctance motor for a wide constant-power speed
range. IEEE Trans Ind Electron 59(6):24952502
References 181
17. Lee D-H, Pham TH, Ahn J-W (2013) Design and operation characteristics of four-two pole
high-speed SRM for torque ripple reduction. IEEE Trans Ind Electron 60(9):36373643
18. Flieller D, Nguyen NK, Wira P, Sturtzer G, Abdeslam DO, Merckle J (2014) A self-learning
solution for torque ripple reduction for nonsinusoidal permanent-magnet motor drives based
on articial neural networks. IEEE Trans Ind Electron 61(2):655666
19. Lei G, Xu W, Hu JF, Zhu JG, Guo YG, Shao KR (2014) Multi-level design optimization of a
FSPMM drive system by using sequential subspace optimization method. IEEE Trans Magn
50(2). Article 7016904
20. Yao D, Ionel DM (2013) A review of recent developments in electrical machine design
optimization methods with a permanent-magnet synchronous motor benchmark study. IEEE
Trans Ind Appl 49(3):12681275
21. Hasanien HM, Abd-Rabou AS, Sakr SM (2010) Design optimization of transverse flux linear
motor for weight reduction and performance improvement using response surface
methodology and genetic algorithms. IEEE Trans Energy Convers 25(3):598605
22. Buja GS, Kazmierkowski MP (2004) Direct torque control of PWM inverter-fed AC motors
a survey. IEEE Trans Ind Electron 51(4):744757
23. Kouro S, Cortes P, Vargas R, Ammann U, Rodriguez J (2009) Model predictive controla
simple and powerful method to control power converters. IEEE Trans Ind Electron 56
(6):18261838
24. Bolognani S, Bolognani S, Peretti L, Zigliotto M (2009) Design and implementation of model
predictive control for electrical motor drives. IEEE Trans Ind Electron 56(6):19251936
25. Xia CL, Wang YF, Shi TN (2013) Implementation of nite-state model predictive control for
commutation torque ripple minimization of permanent-magnet brushless DC motor. IEEE
Trans Ind Electron 60(3):896905
26. Morel F, Lin-Shi XF, Retif J-M, Allard B, Buttay C (2009) A comparative study of predictive
current control schemes for a permanent-magnet synchronous machine drive. IEEE Trans Ind
Electron 56(7):27152728
27. Liu C-H, Hsu Y-Y (2010) Design of a self-tuning PI controller for a STATCOM using particle
swarm optimization. IEEE Trans Ind Electron 57(2):702715
28. Zhu JG, Guo YG, Lin ZW, Li YJ, Huang YK (2011) Development of PM transverse flux
motors with soft magnetic composite cores. IEEE Trans Magn 47(10):43764383
29. Guo YG, Zhu JG, Watterson PA, Wei Wu (2006) Development of a PM transverse flux motor
with soft magnetic composite core. IEEE Trans Energy Convers 21(2):426434
30. Wang TS, Zhu JG, Zhang YC (2011) Model predictive torque control for PMSM with duty
ratio optimization. In: Proceedings of the 2011 International Conference on Electrical
Machines and Systems (ICEMS), 2023 Aug 2011, pp 15
31. Wang SH, Meng XJ, Guo NN, Li HB, Qiu J, Zhu JG et al (2009) Multi-level optimization for
surface mounted PM machine incorporating with FEM. IEEE Trans Magn 45(10):47004703
32. Meng XJ, Wang SH, Qiu J, Zhang QH, Zhu JG, Guo YG, Liu DK (2011) Robust multi-level
optimization of PMSM using design for six sigma. IEEE Trans Magn 47(10):32483251
Chapter 6
Design Optimization for High Quality
Mass Production
Abstract In the last two chapters, the design optimization methods under the
framework of deterministic approach were presented for electrical machines and
drive systems. By the deterministic approach, all material and structural parameters
in the manufacturing process are exact values that do not have any variations from
their nominal values. However, there are many unavoidable uncertainties or vari-
ations in the industrial manufacturing process of electrical machines and drive
systems, including mainly material diversity, manufacturing errors and assembly
inaccuracy. These will result in big variations affecting the reliability and quality of
electrical machines and drive systems in mass production. These variations are not
investigated in the deterministic approach. The main aim of this chapter is to
present a robust approach based on the technique of design for six-sigma (DFSS)
for the design optimization of high-performance and high-quality electrical
machines and drive systems in mass production. Meanwhile, two multi-level
optimization strategies are presented to improve the optimization efciency for high
dimensional problems. Through the investigation of several design examples, it is
shown that the reliability and quality of the investigated electrical machines and
drive system can be increased greatly by using the proposed robust approach.
Keywords Design for six-sigma Electrical drive systems Model predictive
control Robust design optimization
System-level design optimization
Multi-objective optimization
Transverse flux machine
Permanent magnet
motors Monte carlo analysis
6.1 Introduction
In general, there are three kinds of design approaches from industry perspective,
namely deterministic, reliability and robust approaches. Figure 6.2 illustrates a
general optimization flowchart and features for them.
Figure 6.2a illustrates the optimization solutions obtained from the robust and
deterministic approaches, respectively. As shown, the optimal solution obtained by
deterministic approach (namely function minimum) is smaller than the robust
optimum. However, when a variation or noise x happens, the objectives variation
f(x) of the deterministic approach is obviously larger than that of the robust
approach. Most importantly, some variations from deterministic approach violate
the basic constraints of design problem, e.g. lying inside the infeasible domain, and
this is prohibited in engineering design [2].
Figure 6.2b illustrates the optimizing flowcharts for these three design approa-
ches. As shown, the distance between constraint boundary and the solution obtained
from robust approach is the furthest one, which means the reliability of the product
is the highest. Meanwhile, the objectives variance of robust approach is the
smallest, which means that the quality variance of robust approach is also the
smallest and the products quality is the highest. Generally, if a design scheme is
not robust, it may be very difcult or even impossible to manufacture (e.g. requiring
extreme material characteristics or unrealistically high manufacturing precision) due
to current manufacturing technology or to operate (e.g. unstable system perfor-
mance in the application environment).
As shown in Fig. 6.2b, the deterministic approach tends to push a design toward
one or more constraint boundaries until those constraints are reached, which pro-
vides a high-risk design to the designer. Furthermore, the deterministic approach
tends to search for the valley solutions or global minimal values from the point of
view of mathematics. However, the valley point is highly sensitive to design
parameter variations, i.e. the product performance will be degraded signicantly in
practical industry manufacturing [2]. Therefore, the robust approach is very
important for modern quality control and design and should be taken into account in
the system-level design optimization of drive systems.
Generally, a deterministic design with respect to an objective f(x) and m con-
straints g(x) has the form as
min : f x
s:t: gi x 0; i 1; . . .; m ; 6:1
xl x x u
where xl and xu are the boundaries of design parameter x which is deterministic and
does not cover any uncertain information. As mentioned above, there are many
unavoidable noise factors in the industrial design and manufacturing process, such
as the assembly tolerances and manufacturing imperfections in mass-production
[46, 2023]. Therefore, reliability design is developed to include the noise factors
in the constraints to improve the reliability of products, in which g(x) is converted
to a probability function as
Pf Pgx [ 0 PU ; 6:2
where x is a vector of random variables, Pf the failure probability, and PU its upper
bound [2].
However, reliability design just focuses on the constraint boundary, and does not
consider the variations of objectives and constraints in terms of those noise factors.
Therefore, the quality distribution and average product performance cannot be
evaluated. Fortunately, DFSS technique can deal with these problems very well.
Actually, DFSS is a kind of robust design approach which was originated from the
Six-Sigma Methodology developed by MOTOROLA and GE [16]. It is generally
188 6 Design Optimization for High Quality Mass Production
used to develop products to meet customer needs with very low defect levels. It has
the design form as
min : Flf x; rf x
s:t: : gi lf x; rf x 0; i 1; . . .; m
; 6:3
xl nrx lx xu nrx
LSL lf nrf USL
where and are the mean and standard deviation of the corresponding terms
which are generally estimated by MCA method, LSL and USL the lower and supper
specication limits, n is the sigma level, which is generally equivalent to a prob-
ability of a standard normal distribution as shown in Fig. 6.3.
Table 6.1 tabulates the equivalent percentage/probability for each sigma level.
For example, 3 means that the probability of P3r x 3r is 99.73 %
assuming that x follows a standard normal distribution with mean 0 and variation 1
(r 1). In other words, 3 is equivalent to a probability of 99.73 % or the POF is
0.27 %. This probability was deemed acceptable in statistical terms, and this value
can be regarded as the quality control of short term, which means that there are
2,700 defects per million products.
However, with the development of long term quality control and management, this
3 quality level is insufcient from the manufacturing perspective. From the rich
experience of MOTOROLA, GE and others, an approximate 1.5 sigma shift in the mean
(as shown in Fig. 6.4) was observed and this has been used to dene the long term sigma
quality as opposed to the above short term sigma quality [2, 1618, 2427].
For example, if there is a 1.5 shift for 3 quality control, the equivalent
probability is
0.2
0.1
1
-6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6
68.27%
99.73% (3 )
99.9999998% (6 )
6.2 Design for Six-Sigma 189
Table 6.1 Percentages and defects per million in terms of sigma level
Sigma level () Percentage Defects per million (short term) Defects per million (long term)
1 68.26 317,400 697,700
2 95.46 45,400 308,733
3 99.73 2,700 66,803
4 99.9937 63 6,200
5 99.999943 0.57 233
6 99.9999998 0.002 3.4
0.3
0.1
-6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6
Therefore, the failure rate is 6.6803 %, which means that the 3 quality control is
actually equivalent to 66,803 defects per million products being manufactured in
the long term quality control technique. Obviously, this quality is not acceptable for
the mass production of a product in industry. Similarly, 4 and 5 are equivalent to
6,200 and 233 defects per million, respectively, and they are not good choices for
quality control either. To achieve the highest prot, 6 level should be selected as it
is equivalent to only 3.4 defects per million products, and it has been adopted in
many companies worldwide nowadays.
To compare the product reliability by using different design approaches, a criterion
called product probability of failure (POF) was used in many works. Assuming that all
constraints in (6.3) are independent events, and then according to the Multiplication
Theorem of Probability, the POF of the system described by (6.3) has the form as
Y
m
POF 1 Pgi 0 6:5
i1
where Pgi 0 means that the probability of event constraint gi is correct for all
samples in the MCA [2].
190 6 Design Optimization for High Quality Mass Production
With the robust optimization framework of (6.3), the robust optimization model
of (6.6) can be expressed as
min : lf x
s:t: : lgi x 6rgi x 0; i 1; . . .; 4
xl 6rx lx xu 6rx 6:7
lf 6rf LSL
lf 6rf USL
MCA is used to estimate the mean and standard deviation terms in (6.7), and the
sample size is 104. It should be noted that the optimization parameters in (6.6) and
(6.7) are discrete values, and their step sizes are shown in Table 6.2 as well.
Table 6.4 Sigma levels for constraints and POF for TFM after robust optimization
Par. g1 g2 g3 g4 POF (%)
Deter. 6 6 6 50.37 % 49.63
Robust 6 6 6 6 *0
less than 1. Actually, the current density is 6.00 A/mm2, which is the same as
the limit of this constraint. As a result, the POF of motor is 49.63 %. For the
robust scheme, the sigma levels for all constraints are larger than 6 and the
POF is almost 0.
(4) Figures 6.5, 6.6, 6.7, 6.8 show the distributions of cost, output power, ef-
ciency and current density respectively for both methods. As shown, the
standard deviations of cost and output power of robust optimal scheme are
smaller than those of the deterministic optimization scheme. As shown in
Fig. 6.7, the robust optimal scheme can produce a larger mean and a smaller
standard deviation for the efciency of this TFM compared with the deter-
ministic scheme.
As shown in Fig. 6.8, all current distribution points of robust design scheme
are satised with the condition g4 of no larger than 6.0 A/mm2. It can also be
seen that many points of deterministic design scheme are not satised with
this condition. Therefore, the reliability and sigma level of this constraint of
deterministic method is very low, and the POF of motor is high. Actually, the
lower cost of deterministic optimization is obtained at the cost of low relia-
bility and robustness. This is not acceptable in engineering design.
6.3 Robust Design Optimization of Electrical Machines 193
0.2
Probability
0.15
0.1
0.05
0
27 27.5 28 28.5 29 29.5 30
Cost [$]
0.2
Probability
0.15
0.1
0.05
0
680 690 700 710 720 730 740 750 760
Pout [W]
0.2
Probability
0.15
0.1
0.05
0
0.815 0.82 0.825 0.83 0.835 0.84
Efficiency
194 6 Design Optimization for High Quality Mass Production
0.2
Probability
0.15
0.1
0.05
0
5.6 5.7 5.8 5.9 6 6.1 6.2 6.3
Jc [A/mm2]
Section 6.3.1 presents a robust design optimization for a PM-SMC TFM under
single-objective optimization situation. From the discussion, it can be found that the
manufacturing quality of the motor has been increased signicantly. However, two
other issues should be investigated for the industrial applications of PM-SMC
motors besides the robust analysis. Firstly, multi-objective design schemes are
necessary as it is hard to determine the weights for different objectives without
detailed information of industrial applications. Secondly, high computational cost is
also an important issue as this is a high dimensional optimization problem and 3D
nite element model (FEM) is involved. Therefore, this section presents an
improved multi-objective sequential optimization method (MSOM) for the robust
multi-objective design optimization of these PM-SMC motors to improve their
industrial applications. The PM-SMC TFM discussed in the last section will be
investigated here to illustrate the proposed new method.
Considering the manufacturing condition and material characteristic of SMC
cores, the multi-objective design optimization model of PM-SMC motors can be
dened as
f1 x Cost
min :
8 2 x Pout
f
>
> g1 x 0:795 g 0;
< ; 6:8
g2 x 640 Pout 0;
s:t:
>
> g x sf 0:8 0;
: 3
g4 x Jc 6 0:
6.3 Robust Design Optimization of Electrical Machines 195
where xs, xmt, and xmf are the structure, material and manufacturing parameters,
respectively. To achieve the six sigma quality manufacturing, the design model can
be converted into (6.9) within the framework of DFSS [2, 18].
min : Fk lfk x; k 1; 2
: 6:9
s:t: lgi x nrgi x 0; i 1; . . .; 4
where and stand for the mean and standard deviation of the corresponding
terms.
Figures 6.9, 6.10, 6.11, 6.12 illustrate the optimization results for both deterministic
and robust design approaches. The following conclusions can be drawn:
(1) Figure 6.9 illustrates the Pareto optimal solution obtained from deterministic
model (6.8) and robust model (6.9), respectively. It can be found that the
output power increases with the increase of cost and vice versa. The front of
196 6 Design Optimization for High Quality Mass Production
760
Pout [W]
720
680
640 Deterministic
Robust
25 27 29 31 33 35
Cost [$]
0.6
0.4
0.2
0 10 20 30 40 50
Number of Pareto Points
6.1
Mean of current density (A/mm2)
5.9
5.8
5.7
5.6
5.5
0 10 20 30 40 50
Number of Pareto Points
6.3 Robust Design Optimization of Electrical Machines 197
6.8
6.6
0 10 20 30 40 50
Number of Pareto Points
Pareto solutions obtained by the robust approach is lower than that by the
deterministic approach, meaning that to achieve the same output power, the
needed cost of robust design approach is higher than that of deterministic one.
(2) Figure 6.10 illustrates the POF values of all Pareto points for both approaches.
It can be found that the POF values of deterministic design approach (or Pareto
points) are unstable and higher than those of robust approach. Some of them
are even over 50 % higher. These are bad designs from the point of view of
high quality industrial design.
For the robust multi-objective design approach, the POF values are almost 0.
Therefore, although the needed cost for the same output power of the deter-
ministic scheme is less than that of the robust approach, its lower cost is
achieved at the cost of lower POF.
(3) Figure 6.11 shows the means of current density (Jc) for all Pareto points. It can
be seen that deterministic designs have higher means of Jc. The means of Jc of
the robust designs are obviously smaller than the limit of 6 A/mm2, and the
average of these means is 5.67 A/mm2. However, many of the deterministic
designs are beyond the limit, and the average is 5.99 A/mm2. Thus, the POF
values of g4 of the deterministic approach are higher than those of the robust
approach. For other constraints, the POFs and means for all Pareto points are
also obtained, but not much difference has been founded. Therefore, the
current density issue is the main reason why the deterministic approach has
higher POFs than the robust approach as shown in Fig. 6.10.
(4) Figure 6.12 illustrates the core density for all Pareto points. It can be found
that the core densities of deterministic designs are around 7.2 g/cm3, which
means that a 200-ton press is needed to compact these cores. On the other
hand, the core densities of all robust designs are around 6.6 g/cm3, and only a
100-ton press is required. Therefore, the robust approach needs lower manu-
facturing condition and cost than the deterministic approach.
198 6 Design Optimization for High Quality Mass Production
(5) In terms of the computational cost, the direct optimization method (NSGA II
with FEM) requires about 12,000 FEM points, in which half points were
sampled for the no-load analysis and others were used for the operation
analysis of this machine. The proposed method requires only about 3,800
FEM points, which is much less than that required by the direct optimization
method.
In summary, the proposed robust multi-objective optimization method can sig-
nicantly improve the reliability and manufacturing quality of the motor with lower
manufacturing condition and cost. Consequentially, it will promote signicantly the
industrial applications of PM-SMC motors.
Fig. 6.13 Single level robust optimization method for drive systems
6.4 Robust Design Optimization of Electrical Drive Systems 199
To estimate the mean and standard deviation terms in (6.10), each design
parameter xi in xs is assumed normally distributed as Nxi ; r2xi with rxi Dxi =3,
where Dxi is the manufacture tolerance of xi.
Step 2: Selection of an optimization method for model (6.10)
As drive systems are always high dimensional and non-linear design problems,
intelligent algorithms can be good choices in many situations, such as the genetic
algorithm (GA) and DEA [13, 3537].
Step 3: Implementation of optimization process
Firstly, determine the manufacture tolerance for xs and obtain the distribution
parameters for each parameter. Secondly, generate an initial population of xs and its
noise population. Thirdly, evaluate the steady-state and dynamic performance
parameters of the drive system. Meanwhile, the objectives, constraints in (6.10),
and their means and standard deviations can be gained by using the MCA method.
Finally, apply the optimization algorithm until the convergence criterion is met.
Because, as mentioned above, the computational cost of this single-level method
is always huge as these design problems are generally high dimensional and non-
linear, the computational cost of whole system optimization is very expensive, in
which the major part is the computational cost of FEM for the motors. To solve this
problem, a multi-level robust optimization method is proposed.
Figure 6.14 illustrates the framework of multi-level robust optimization method for
drive systems, which consists of three levels, namely the motor, control, and system
levels. In the implementation, the rst step is to dene the deterministic and robust
optimization models for the motor and control levels, respectively. For the motor
level, its deterministic optimization model has been dened in Chap. 5, and has the
form as
200 6 Design Optimization for High Quality Mass Production
min : fm xm
s:t: gmi xm 0; i 1; . . .; Nm ; 6:11
xml xm xmu
where xm, fm and gm are the design parameter vector, objective and constraint,
and xml and xmu the lower and upper boundaries of xm, respectively.
With the DFSS technique, the corresponding robust model can be expressed as
where nm is the sigma level specied for the motor level. It should be noted that the
required sigma levels for motor level and control level may be different from that in
system level model, and different symbols are specied for them.
For the control level, its deterministic design optimization model has the form as
min : fc xc
s:t: gcj xc 0; j 1; . . .; Nc ; 6:13
xcl xc xcu
where xc, fc and gc are the design parameter vector, objective and constraint, and xcl
and xcu the lower and upper boundaries, respectively.
Similarly, its robust model can be expressed as
inductance and PM flux, should be calculated at the same time, because they will be
considered as the noise factors in the optimization of control level.
Step 3: Optimization of control level with model (6.14)
The main aim of this step is to evaluate the fluctuations of dynamic performance
with respect to the noise factors from the motor characteristic parameters.
Step 4: Performance evaluation of system level
Here are two remarks for this multi-level robust optimization method. Firstly, if
the dimension of motor level or control level is large, the optimization of corre-
sponding level can be divided into two or three sub-levels. Secondly, as mentioned
before, the optimization process is usually quite time-consuming because of the
huge computational costs of FEM for electromagnetic analysis of the motor and the
MCA process in robust optimization. Approximate models, such as the response
surface model and Kriging model, can be used to replace the FEM [3234]. The
Kriging model will be employed for the design example in this work.
min : fm xm w1 Cost
C0 w2 Pout
P0
where C0 and P0 are the cost and output power (Pout) of the initial prototype, is
motor efciency, sf the winding ll factor, and Jc the current density of copper wire
winding. The optimization parameters are illustrated in Table 5.1 and Fig. 5.5.
B. Optimization model for control level
Figure 2.37 showed a diagram of this improved MPC control scheme. There are
two important modules needed to be designed and optimized in this improved
MPC, namely cost function and duty ratio module [3840]. Six parameters pre-
sented in Sect. 5.4.2 should be optimized in the control level. They are A, N, CT, C,
6.4 Robust Design Optimization of Electrical Drive Systems 203
Kp and Ki, where Kp and Ki are the parameters of PI controller. One objective and
four constraints are considered for this level as the following
min : fc xc w3 RMSE(
Trated
T
w4 RMSE(
nrated
n
w5 nos
s:t: : gc1 xc RMSE(T=Trated 0:06 0;
gc2 xc RMSE(n=nrated 0:05 0; ; 6:16
gc3 xc nos 0:02 0;
gc4 xc ts 0:02 0;
xcl xc xcu
where RMSE is the root mean square error of an item in the steady operation
period, the motor speed, nos the overshoot of speed, and ts the settling time.
C. Robust optimization models
Firstly, based on (6.15), the robust optimization model for the motor level has the
form as
min : 8
lfm xm
< lgmi xm nm rgmi xm 0; i 1; . . .; 4
>
xmlj nm rxmj lxmj xmuj nm rxmj ; : 6:17
s:t:
>
:
j 1; . . .; 8
min : lfc xc
lgci xc nc rgci xc 0; i 1; . . .; 4 : 6:18
s:t:
xclj xcj xcuj ; j 1; . . .; 6
It should be noted that the parameters of control level (xc1 to xc6) are digital
parameters and do not have disturbances for the MCA. There are only four robust
constraints in (6.18).
However, there are four random variables in this model, namely the resistance
(R), inductance (L), torque (T), and PM-flux corresponding to the random variables
of all motor parameters xm1 to xm8. After the MCA, the output parameters are the
means and deviations of objectives and constraints of control level. Then, model
(6.18) can be calculated and the POF of control level can be obtained.
204 6 Design Optimization for High Quality Mass Production
Probability
0.2
0.15
0.1
0.05
0
5.6 5.8 6 6.2
2
Jc [A/mm ]
(b)
0.3
Robust
0.25
Probability
0.2
0.15
0.1
0.05
0
5.5 5.6 5.7 5.8 5.9 6 6.1 6.2
Jc [A/mm2]
83.45 %, and the corresponding sigma level is 1.4. As a result, the POF of
motor is 16.55 %. As an example, Fig. 6.16 illustrates the current density
distribution of copper wire corresponding to constraint gm4. As shown, many
points in Fig. 6.16a are beyond the limit of current density, 6 A/mm2.
Meanwhile, the reliabilities of constraints gc2 and gc4 are 99.90 % and 19.40 %, and
the corresponding sigma levels are 3.3 and 0.2, respectively. As a result, the POF of
control level is 80.62 %. As another example, Fig. 6.17 illustrates the distribution of
the settling time corresponding to constraint gc4. As shown, most of the points
violate constraint gc4 of no more than 0.02 s after the load is applied. This results
in a high POF for the control level. It should be noted that the last column
(t = 0.14 s) in deterministic gure shows the probability for the issue that settling
times are no less than 0.14 s instead of exactly equaling 0.14 s. Figure 6.18
illustrates the distribution of g2, which is related to speed. As shown, there are also
several points violating the limit no more than 0.05.
208 6 Design Optimization for High Quality Mass Production
0.8
0.6
Probability Deterministic
0.4
0.2
0
0.06 0.08 0.1 0.12 0.14
Robust
0.8
0.6
Probability
0.4
0.2
0
0.06 0.08 0.1 0.12 0.14
Settling time [s]
6.4 Robust Design Optimization of Electrical Drive Systems 209
0.2
Probability
0.15
0.1
0.05
0
0 0.01 0.02 0.03 0.04 0.05 0.06
RMSE(n) / nrated
(b) 1
0.8
Probability
0.6
Robust
0.4
0.2
0
0 0.01 0.02 0.03 0.04 0.05 0.06
RMSE(n) / nrated
Finally, the POF of the whole drive system is 83.83 %. This is absolutely not an
acceptable system design for engineering applications.
(2) By the robust approach, three iteration processes are required to get the
optimal results of the multi-level optimization method. Figure 6.19 shows the
iteration process of multi-level optimization for the motor level. As shown,
level 1 is optimized twice while level 2 is optimized only once. The optimum
motor efciency is 81.1 %, and output power 671 W, which are close to those
obtained by the deterministic approach, whereas the material cost is $28,
which is bigger than that of deterministic design. Regarding the control level,
as shown in Table 6.6, the dynamic performance of this drive system after
robust optimization is slightly worse than that of the deterministic optimiza-
tion. However, all comply well with those constrains in control level.
210 6 Design Optimization for High Quality Mass Production
Motor objective
1.6
1.4
1.2
1
Initial Level 1 - 1st Level 2 - 1st Level 1 - 2nd
Iteration number of multilevel optimization
optimization
0
-5
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14
Time [s]
2000
Speed [r/min]
1000
0
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14
Time [s]
Figure 6.20 illustrates the dynamic performance of drive system by using the
algorithm parameters obtained from the robust optimization approach.
After MCA, the reliabilities of all constraints are 1, and the corresponding sigma
levels are all more than 6. As a result, the POFs for the motor, controller, and whole
system are almost zero. This is much better than those obtained by the deterministic
design approach and satises the initial reliability requirements. Therefore, the POF
of deterministic optimization is bigger than the robust optimization. As a matter of
fact, the lower cost of deterministic optimization is obtained at the expense of low
reliability and robustness. This is not acceptable in engineering design.
In summary, compared with the initial design, the solutions of robust opti-
mization have many improvements, such as higher output power and efciency,
lower cost, and better reliability and robustness. Meanwhile, the objectives obtained
by the deterministic design are smaller than those by the robust design optimization,
6.4 Robust Design Optimization of Electrical Drive Systems 211
but the reliabilities, robust levels and POFs of the motor, control, and system levels
are obviously worse than the latter. Finally, the obtained control parameters by the
robust multi-level optimization are not sensitive to the disturbances of motor output
parameters. Therefore, the system dynamic performances can be ensured by the
proposed method. This is very valuable for engineering batch production.
6.5 Summary
The manufacturing quality and reliability are two terms used to evaluate the vari-
ations of motor performance against the material variations, manufacturing toler-
ances, and assembling errors, which is very important for the industry. The robust
approach based on DFSS has been presented for the design optimization of elec-
trical machines and drive systems so as to improve their manufacturing quality and
reliability in mass production in this chapter. Several cases of single- and
multi-objective optimization are investigated. The multi-level robust design opti-
mization method has been presented for the system-level design optimization of
electrical drive systems. From the investigation of the design example, it can be
seen that this approach can signicantly improve the drive system reliability, which
will benet the manufacturing of those devices and extent their applications in
industry. Therefore, the robust design optimization is a valuable and necessary step
for design optimization of electromagnetic devices and systems from the perspec-
tive of engineering design, which can improve the product reliability and quality,
and save the design cost and cycle.
References
1. Lei G, Wang TS, Guo YG, Zhu JG, Wang SH (2014) System level design optimization
method for electrical drive system: deterministic approach. IEEE Trans. Ind. Electron. 61
(12):65916602
2. Lei G, Wang TS, Zhu JG, Guo YG, Wang SH (2015) System level design optimization
method for electrical drive system: robust approach. IEEE Trans Ind Electron 62(8):4702
4713
3. Lei G, Guo YG, Zhu JG et al (2012) System level six sigma robust optimization of a drive
system with PM transverse flux machine. IEEE Trans Magn 48(2):923926
4. Khan MA, Husain I, Islam MR, Klass JT (2014) Design of experiments to address
manufacturing tolerances and process variations influencing cogging torque and back EMF in
the mass production of the permanent-magnet synchronous motors. IEEE Trans Ind Appl 50
(1):346355
5. Coenen I, Giet M, Hameyer K (2011) Manufacturing tolerances: estimation and prediction of
cogging torque influenced by magnetization faults. In: Proceedings of 14th european
conference on power electronics and applications, pp 19, Aug./Sep. 2011
6. Gasparin L, Fiser R (2011) Impact of manufacturing imperfections on cogging torque level in
PMSM. In: Proceedings of 2011 IEEE ninth international conference on power electronics and
drive systems (PEDS), pp 10551060
212 6 Design Optimization for High Quality Mass Production
7. Lei G, Liu CC, Guo YG, Zhu JG (2015) Multi-disciplinary design analysis and optimization of
a PM transverse flux machine with soft magnetic composite core. IEEE Trans Magn 51(11),
Article 8109704
8. Vese I, Marignetti F, Radulescu MM (2010) Multiphysics approach to numerical modeling of
a permanent-magnet tubular linear motor. IEEE Trans Ind Electron 57(1):320326
9. Kreuawan S, Gillon F, Brochet P (2008) Optimal design of permanent magnet motor using
multi-disciplinary design optimisation In: Proceedings of 18th international conference on
electrical machines, Vilamoura, pp 1-6, 69 Sept 2008
10. Taguchi G, Chowdhury S, Wu Y (2004) Taguchis quality engineering handbook. Wiley, New
Jersey
11. Omekanda AM (2006) Robust torque and torque-per-inertia optimization of a switched
reluctance motor using the Taguchi methods, IEEE Trans Ind Appl, 42(2):473478
12. Kim S-I, Lee J-Y, Kim Y-K et al (2005) Optimization for reduction of torque ripple in interior
permanent magnet motor by using the Taguchi method. IEEE Trans Magn 41(5):17961799
13. Ashabani M, Mohamed YA, Milimonfared J (2010) Optimum design of tubular
permanent-magnet motors for thrust characteristics improvement by combined
Taguchi-neural network approach. IEEE Trans Magn 46(12):40924100
14. Hwang C-C, Lyu L-Y, Liu C-T, Li P-L (2008) Optimal design of an SPM motor using genetic
algorithms and Taguchi method. IEEE Trans Magn 44(11):43254328
15. Lee S, Kim K, Cho S et al (2014) Optimal design of interior permanent magnet synchronous
motor considering the manufacturing tolerances using Taguchi robust design. IET Electr
Power Appl 8(1):2328
16. Koch PN, Yang RJ, Gu L (2004) Design for six sigma through robust optimization. Struct.
Multidiscip. Optim. 26(34):235248
17. Lei G, Zhu JG, Guo YG, Hu JF, Xu W, Shao KR (2013) Robust design optimization of
PM-SMC motors for Six Sigma quality manufacturing. IEEE Trans Magn 49(7):39533956
18. Lei G, Zhu JG, Guo YG, Shao KR, Xu W (2014) Multi-objective sequential design
optimization of PM-SMC motors for Six Sigma quality manufacturing. IEEE Trans Magn 50
(2):7017704
19. Meng XJ, Wang SH, Qiu J, Zhang QH, Zhu JG, Guo YG, Liu DK (2011) Robust multi-level
optimization of PMSM using design for six sigma. IEEE Trans Magn 47(10):32483251
20. Gasparin L, Cernigoj A, Markic S, Fiser R (2008) Prediction of cogging torque level in PM
motors due to assembly tolerances in mass-production. Int J. Comput Math Elect Electron Eng
(COMPEL) 27(4):911918
21. Gasparin L, Cernigoj A, Markic S, Fiser R (2009) Additional cogging torque components in
permanent-magnet motors due to manufacturing imperfections. IEEE Trans Magn 45
(3):12101213
22. Islam MS, Islam R, Sebastian T, Ozsoylu S, Chandy A (2010) Cogging torque minimization in
PM motors using robust design approach. In: Proceedings of IEEE conference on energy
conversion congress and exposition (ECCE), pp 38033810
23. Islam MS, Mir S, Sebastian T (2004) Issues in reducing the cogging torque of mass-produced
permanent-magnet brushless DC motor. IEEE Trans Ind Appl 40(3):813820
24. Tjahjono B, Ball P, Vitanov VI et al (2010) Six Sigma: a literature review. Int J Lean Six
Sigma 1(3):216233
25. Reosekar RS, Pohekar SD (2014) Six Sigma methodology: a structured review. Int J Lean Six
Sigma 5(4):392422
26. Aboelmaged MG (2010) Six Sigma quality: a structured review and implications for future
research. Int J Quality Reliab Manag 27(3):268317
27. Wang SH, Meng XJ, Guo NN, Li HB, Qiu J, Zhu JG et al (2009) Multi-level optimization for
surface mounted PM machine incorporating with FEM. IEEE Trans Magn 45(10):47004703
28. Zhu JG, Guo YG, Lin ZW, Li YJ, Huang YK (2011) Development of PM transverse flux
motors with soft magnetic composite cores. IEEE Trans Magn 47(10):43764383
29. Guo YG, Zhu JG, Watterson PA, Wei Wu (2006) Development of a PM transverse flux motor
with soft magnetic composite core. IEEE Trans Energy Convers 21(2):426434
References 213
30. Guo YG, Zhu JG, Dorrell D (2009) Design and analysis of a claw pole PM motor with molded
SMC core. IEEE Trans Magn 45(10):45824585
31. Deb K, Pratap A, Agarwal S, Meyarivan T (2002) A fast and elitist multi-objective genetic
algorithm: NSGA-II. IEEE Trans Evol Comput 6(2):182197
32. Yao D, Ionel DM (2013) A review of recent developments in electrical machine design
optimization methods with a permanent-magnet synchronous motor benchmark study. IEEE
Trans Ind Appl 49(3):12681275
33. Mendes MHS, Soares GL, Coulomb J-L, Vasconcelos JA (2013) Appraisal of surrogate
modeling techniques: a case study of electromagnetic device. IEEE Trans Magn 49(5):1993
1996
34. Wang LD, Lowther DA (2006) Selection of approximation models for electromagnetic device
optimization. IEEE Trans Magn 42(2):12271230
35. Meng XJ, Wang SH, Qiu J, Zhu JG, Wang Y, Guo YG et al (2010) Dynamic multi-level
optimization of machine design and control parameters based on correlation analysis. IEEE
Trans Magn 46(8):27792782
36. Hasanien HM, Abd-Rabou AS, Sakr SM (2010) Design optimiz-ation of transverse flux linear
motor for weight reduction and perfo-rmance improvement using response surface
methodology and GAs. IEEE Trans Energy Convers 25(3):598605
37. Hasanien HM (2011) Particle swarm design optimization of transverse flux linear motor for
weight reduction and improvement of thrust force. IEEE Trans Ind Electron 58(9):40484056
38. Bolognani S, Bolognani S, Peretti L, Zigliotto M (2009) Design and implementation of model
predictive control for electrical motor drives. IEEE Trans Ind Electron 56(6):19251936
39. Morel F, Lin-Shi XF, Retif J-M, Allard B, Buttay C (2009) A comparative study of predictive
current control schemes for a permanent-magnet synchronous machine drive. IEEE Trans Ind
Electron 56(7):27152728
40. Wang TS, Zhu JG, Zhang YC (2011) Model predictive torque control for PMSM with duty
ratio optimization. In Proceedings of 2011 international conference on electrical machines and
systems (ICEMS), 20-23 Aug. 2011, pp. 15
Chapter 7
Application-Oriented Design Optimization
Methods for Electrical Machines
Keywords Application-oriented design optimization method Motor topologies
Hybrid electric vehicles Transverse flux machine Flux-switching permanent
magnet machine
7.1 Introduction
Since each vehicle company has its own design conception, the requirements for
drive motors of different companies are different. On the other hand, there are many
motor types and design schemes for each application. Therefore, for a given
application, all possible motor topologies and structural parameters should be
investigated to get a globally optimum among different options.
This chapter presents an application-oriented design optimization method for
novel electrical machines for domestic appliances and HEVs, respectively.
Section 7.2 presents the proposed optimization method for electrical machines by
the deterministic approach. Section 7.3 presents a robust approach for the
application-oriented design optimization method with a design example of the
plug-in HEV (PHEV) drive, followed by the remarks and a summary.
Fig. 7.1 Framework of application-oriented design optimization method for electrical machines
218 7 Application-Oriented Design Optimization
algorithm (GA) and particle swarm optimization (PSO) algorithm [8, 9],
and response surface model (RSM) and Kriging model [1012]. It should
be noted that, sequential optimization method can be employed here to
improve the optimization efciency for the low-dimensional design
problems, and multi-level optimization method can be employed for
high-dimensional design problems [1217].
Step 5: Compare the optimal results of all options, and output the best one as the
nal optimal solution for that specic application.
Refrigerators are commonly used domestic appliances, and each has different
specications for its drive machine. Several popular ones are the cost, output power
and efciency. Meanwhile, many types of electrical machines have been designed
for driving refrigerators so far. Therefore, in order to achieve a best design, all
possible motor types and structures should be investigated, and each of them should
be optimized for its best performance. The best ones of these motors will be
compared to nd the nal optimal design.
In our previous work, a PM transverse flux machine (TFM) with SMC core as
shown in Fig. 2.3 has been developed for the drive machine used in a kind of
refrigerator [18, 19]. This refrigerator will be considered as the rst application for
the proposed application-oriented design optimization method. Table 7.1 lists some
of the specications for this application.
This PM-SMC TFM was designed to deliver an output power of 640 W (or
torque of 3.4 Nm) at the rated speed of 1800 rev/min. From the experimental
results, it can be seen that this prototype can present good performance [18, 19].
Meanwhile, this motor was compared with two other commercial motors of lami-
nated cores. One is a high efciency induction motor with the rated torque of
3.72 Nm at 1410 rev/min, 75 % efciency and 80C temperature rise in the coil. Its
outer diameter is 160 mm and total length is 234 mm. The second one is a radial
eld NdFeB brushless DC servo motor with rated torque of 3.45 Nm at the rated
speed of 3000 rev/min. Its outer diameter is 100 mm and axial length 217 mm.
Through the comparison, it can be seen that the proposed SMC motor features a
Fig. 7.2 Analysis model in Ansoft, a one phase motor stack, and b one phase SMC stator stack
Fig. 7.3 Stator claw pole topologies: a rectangular-shoe pole structure, b and c arc pole structures
torque per unit volume of 4.5 times that of the laminated induction motor.
Meanwhile, it delivers 2.25 times the torque per unit volume of the second
brushless DC servo motor [18]. Therefore, the proposed PM-SMC TFM has better
performances than those two commercial motors.
As shown in Fig. 7.2, only one stator phase stack has been investigated for that
PM-SMC TFM in our previous work. In practice, the PM claw pole motor with
SMC core may have various other kinds of stator phase stack structures besides the
original one shown in Fig. 7.2b. Figure 7.3 shows several possible SMC claw pole
stator topologies. The rst one has rectangular-shoe poles, which is derived from
the design of SMC stator for claw pole motors [20, 21]. The second and third are
arc tooth stator. The following work will investigate several stator topologies for
this TFM and present an optimal design by using the proposed application-oriented
design optimization method.
To obtain the optimal SMC core for this PM TFM, all three topology structures
of phase stacks and the dimensions of SMC stator have to be investigated.
Meanwhile, in the manufacturing process, besides the density and dimensions of
SMC stator, the dimensions of PMs are also important to the machine performance.
Therefore, all the factors mentioned above are considered to minimize the material
220 7 Application-Oriented Design Optimization
cost and maximize the output power of this TFM. Based on the specications listed
in Table 7.1, the uniform optimization model for these motors can be dened as
Cost P0
min : f x w1 w2
C0 Pout
s:t: g1 x 0:795 g 0
; 7:1
g2 x 640 Pout 0
g3 x sf 0:8 0
g4 x Jc 6 0
where w1 and w2 are the weighting factors, C0 and P0 the material cost and output
power of the initial prototype, , Pout, sf and Jc the motor efciency, output power,
ll factor and current density, respectively.
Five parameters are selected as the optimization parameters for these three
topology structures of stator phase stacks. Three of them are the circumferential
angle, axial width of PMs, and core density. In the case of rectangular-shoe poles,
the other two parameters are the axial and circumferential widthes of SMC teeth
while the tooth height is xed as 3 mm. In the case of arc poles, the other two
parameters are the inner and outer circumferential angles of SMC teeth. Excluding
the claw pole part, the other parameters of the three structures are the same.
Tables 7.2,7.3,7.4,7.5 and 7.6 list the optimization results for this motor in terms
of different weighting factors. From the tables, the following conclusions can be
drawn:
Table 7.2 Optimization Parameter Unit Original design Arc Rectangular shoe
results for w1 = 0.3 and
Cost $ 35.6 40.5 43.6
w2 = 0.7
Pout W 768 870 878
Efciency % 82.2 80.3 81.6
Density g/cm3 5.84 5.81 5.80
Objective 0.859 0.829 0.849
Table 7.3 Optimization Parameter Unit Original design Arc Rectangular shoe
results for w1 = 0.4 and
Cost $ 32.5 35.1 38.4
w2 = 0.6
Pout W 730 781 802
Efciency % 82.2 79.5 80.6
Density g/cm3 6.15 6.02 5.91
Objective 0.861 0.855 0.876
7.2 Application-Oriented Design Optimization Method 221
Table 7.4 Optimization Parameter Unit Original design Arc Rectangular shoe
results for w1 = 0.5 and
Cost $ 28.7 32.6 32.4
w2 = 0.5
Pout W 672 734 690
Efciency % 83.0 79.5 80.0
Density g/cm3 6.27 6.13 6.35
Objective 0.847 0.858 0.883
Table 7.5 Optimization Parameter Unit Original design Arc Rectangular shoe
results for w1 = 0.6 and
Cost $ 27.1 32.3 30.1
w2 = 0.4
Pout W 640 726 642
Efciency % 83.3 79.5 79.5
Density g/cm3 6.33 6.15 6.55
Objective 0.789 0.848 0.844
Table 7.6 Optimization Parameter Unit Original design Arc Rectangular shoe
results for w1 = 0.7 and
Cost $ 27.1 32.3 30.2
w2 = 0.3
Pout W 642 726 645
Efciency % 83.3 79.5 79.5
Density g/cm3 6.31 6.15 6.54
Objective 0.820 0.853 0.865
(1) For the situation that w1 = 0.3 and w2 = 0.7 (Table 7.2), TFM with original
design stator has the least material cost of $ 35.6 and the highest efciency of
82.2 %, TFM with rectangular-shoe stator has the highest output power of
878 W. However, the best topology structure is the TFM with arc-teeth stator.
The objective is 0.829, which is the minimal one among them and the cor-
responding material cost and output power are $ 40.5 and 870 W, respectively.
(2) For the situation that w1 = 0.4 and w2 = 0.6 (Table 7.3), similarly, the best
topology structure is still the arc teeth stator. The objective is 0.855, which is
the minimal one among them and the corresponding material cost and output
power are $ 35.1 and 781 W, respectively.
(3) For the situation that w1 = 0.5 and w2 = 0.5 (Table 7.4), TFM with original
design stator has the least material cost of $ 28.7 and the highest efciency of
83.0 %, TFM with arc teeth stator has the highest output power of 734 W. The
best topology structure is the TFM with original design stator. The objective is
0.847.
(4) For the situations that w1 = 0.6 (Table 7.5) and w1 = 0.7 (Table 7.6), similarly,
the best topology structure is the TFM with original design stator. The optimal
objectives and other motor performances can be seen in those two tables.
Therefore, the TFM with arc-teeth stator is the best one for the rst two situa-
tions of w1 = 0.3 and 0.4, and the TFM of the original stator is the best one for the
222 7 Application-Oriented Design Optimization
latter three situations of w1 = 0.5, 0.6 and 0.7. The best topology structure and
dimension parameters are related to the weighting factors in the optimization model
chosen according to the applications. In a situation that the cost is more important
than the output power, the TFM with arc teeth stator should be taken as the best
design schemes.
Figure 7.4 shows the framework for the proposed application-oriented design
optimization method based on the robust approach for PM machines. Compared to
the framework under deterministic approach, there are two main differences.
The rst one is the development of robust analysis model for each option. The
robust analysis model mainly includes the determination of variations or noise
factors and the manufacturing tolerances or distribution parameters. The robust
optimization model can be constructed by using the design for Six-Sigma (DFSS)
technique. The second one is that the Monte Carlo analysis (MCA) method is
required in this method to evaluate the manufacturing quality of the motors in mass
production, e.g. one million products per batch, which mainly includes the mean
and standard deviation of variations of motor performance and reliability.
PHEVs have been developed in many countries due to the shortage of fossil fuels.
The electric drive system as one of key units in the PHEV plays crucial role for its
widely successful commercialization. Figure 7.5 shows a novel PHEV powertrain
dependent on one electric machine has been proposed by the University of
Technology Sydney (UTS) since 2009 [4, 22]. It consists of an energy storage unit
comprising of batteries and super-capacitors, a power control unit including the DC
link, DC/DC converters and a back to back inverter/rectier, an electric machine,
functioning as either a motor or a generator (M/G), and an internal combustion
engine (ICE) working mainly during fast acceleration to provide the extra torque.
The system operation is governed by a special energy management strategy as
illustrated in Fig. 7.6 [22], where SOC stands for the state of charge of the energy
storage unit, and EM the electric machine. Initially, it is assumed that the battery
and super-capacitor banks are fully charged from the grid, and the capacity of the
7.3 Robust Approach for the Application-Oriented Design Optimization Method 223
Inverter/
Rectifier
Charger Battery DC/DC
M/G
Wheel
energy storage is designed such that the car could cover a reasonable long range. In
the normal operation mode (high SOC and moderate load), the EM works alone as
the prime mover of the car. When it needs extra torque for fast acceleration, the
internal combustion engine (ICE) will provide the assistance. When SOC drops, the
ICE will recharge the battery through EM while the system is idle. If the load is
high and SOC very low, the ICE will work alone to drive the car and recharge
battery through EM. It can be seen that the EM in different working conditions has
to work continuously. Hence, it must have good attributes of high torque density,
high efciency, strong robustness, and convenience of cooling, etc.
The PHEVs have strict requirements on the drive machine, mainly including
high torque/power density, strong flux weakening ability (wide speed range for
cruising), good mechanical robustness, strong thermal dissipation capability, etc.
[2326]. Due to these requirements, several FSPMMs with different topology
structures have been investigated for a PHEV system in our previous work.
Figure 7.7a illustrates the structure of an FSPMM with 12 stator poles and 10 rotor
7.3 Robust Approach for the Application-Oriented Design Optimization Method 225
Fig. 7.7 Structures for FSPMM with 12/10 poles, a whole motor, and b rotor topology
poles (abbreviated as 12/10 FSPMM). Figure 7.7b shows the structure of the rotor.
Compared with those traditional PM machines, such as rotor surface mounted or
rotor inserted PM machines, FSPMMs have the following main advantages.
(a) Strong thermal dissipation capabilityAs the PMs are inserted in the stator,
they can have greater cross sectional area and are less likely to suffer the
demagnetization problem. The winding current density can reach 78 A/mm2
or even larger. In continuous operation, the stator temperature can be main-
tained well below 125 C, which is in the range of H-class insulation by water
cooling.
(b) Strong structure robustnessSimilar to switched reluctance motors, the rotor
of FSPMM has no PMs or brushes as shown in Fig. 7.7, and therefore is
suitable for high speed operation, e.g. above 20,000 rev/min. For a given
power rating, as the rated speed increases, the machine volume can be reduced
gradually.
(c) Concentrated windingThe edge connection of stator winding is shorter than
distributed ones, which means less copper loss with the same amplitude of
stator current.
(d) High power or torque densitySame as the traditional PM machines, PMs in
the FSPMM are employed to generate the major air gap flux linkage, and the
merit of high power or torque density is retained without extra excited loss
[23, 27, 28].
From the extensive research work, it is found that the combination of stator/rotor
poles is a very important topology issue for the motor performance [5, 6, 2931].
Generally, there are many feasible options for the combinations of stator and rotor
poles of FSPMMs by the following equation.
226 7 Application-Oriented Design Optimization
Fig. 7.8 Rotor structures with different poles: a 11 poles, b 13 poles, and c 14 poles
N s k1 m
k1 ; k2 1; 2; 3. . .; 7:2
Nr Ns k2
where Ns and Nr are the numbers of stator and rotor poles, respectively, and m is the
number of phases. For example, for the FSPMM with 12 stator poles as shown in
Fig. 7.7a, there are several promising numbers of rotor poles, such as 10, 11, 13, 14
and 16 poles, respectively, which have been widely investigated in many research
works [5, 6, 17, 30]. Figure 7.8 shows three other rotor topology structures for this
case. For another example, for the FSPMM with 6 stator poles as shown in Fig. 7.9,
several promising numbers of rotor poles are 5, 7 and 8.
On the other hand, there are also some disadvantages for this kind of machine
compared with the traditional doubly salient structure machines, such as high
cogging torque and odd harmonics in the back electromotive force (EMF). These
will reduce the efciency and increase the torque ripple [7, 22, 2729]. To over-
come these disadvantages, many new topology structures including different
combinations of stator/rotor poles, different laminated structures for sheet steels,
such as the traditional radially-laminated and new axially-laminated structures [7],
and different rotor structures, such as the pole-pairing and pole-skewing [30], have
been investigated for the FSPMMs.
Take the FSPMM with 6/7 poles for example. In order to operate at
2000 rev/min, it should be excited by 233 Hz current. As the speed or frequency
goes up, the core loss will increase greatly. To reduce this negative influence, a new
laminated-structure FSPMM (LSFSPMM) as illustrated in Fig. 7.9b has been
proposed in our previous work [7]. Different from the traditional FSPMM,
LSFSPMM is laminated axially in parallel to the shaft. As shown, the stator
includes 6 respective lamination modules, while the rotor involves 7 modules,
which are all made of 0.3 mm high grain oriented silicon steel sheet (HiB). For high
magnetic permeability of HiB, each phase flux linkage loops along the lamination
and make full use of PMs. Figure 7.10 shows the manufacturing modules for rotor
and stator, respectively. Table 7.7 lists the main dimensions for both machines [7].
7.3 Robust Approach for the Application-Oriented Design Optimization Method 227
Fig. 7.9 Structure diagrams of two machines with 6/7 poles, a traditional FSPMM, and
b LSFSPMM
Fig. 7.10 Manufactured modules for the 6/7 LSFSPMM, a rotor, and b stator
(a) (b)
4
4
Torque T (Nm)
Torque T (Nm)
3
3
2 12/10 poles
2 6/5 poles
12/14 poles
6/7 poles
12/16 poles
1 6/8 poles 1
0 0
0 90 180 270 360 0 90 180 270 360
Rotor position r (elec.deg) Rotor position r (elec.deg)
Fig. 7.11 Torque curves versus different combinations of stator/rotor poles, a 6-pole stators, and
b 12-pole stators
Size equation is a good method to present a qualitative analysis and comparison for
FSPMMs with different rotor topology structures. The effect of the combinations of
stator/rotor poles on electromagnetic torque can be expressed by the sizing equation
as,
p
2p Nr
Tem ks kd As Bg D2si la cs g 7:3
4 Ns
where ks and kd are the skew and leakage factors, respectively, As is the armature
current electrical loading, Bg the peak flux density at no-load situation, Dsi the inner
diameter of stator, la the active stack length, cs the stator tooth arc factor, and the
efciency of the machine [32].
As shown, the electromagnetic torque is directly proportional to the number of
rotor poles and inversely proportional the number of stator poles. The combinations
of stator/rotor poles have also great impacts on the cogging torque. As stated in
[32], the period number and magnitude of cogging torque in the FSPMMs can be
briefly evaluated by the least common multiple (LCM) and greatest common
divisor (GCD) of Nr and Ns, respectively. A higher LCM and a lower GCD have
been suggested for effective decrease of the cogging torque.
Table 7.8 tabulates the comparison of FSPMMs with two different combinations
of stator/rotor poles based on the size equation. It can be seen that the FSPMM with
12/14 structure has the higher LCM and ratio of Nr over Ns, which can be regarded as
the index of smaller cogging torque and bigger magnitude of electromagnetic torque.
As aforementioned, Table 7.8 only shows a brief comparison for FSPMMs with
different rotor topology structures. Since the efciency and structural parameters are
involved in (7.3) but not included in the performance evaluation in this table, an
accurate quantitative analysis method involving both topology structures and
structural parameters are required. This is the main aim of the next sub-section.
performances for all topology structure of FSPMMs for the same objectives and
constraints. Therefore, the nal optimal FSPMM will possess not only the best
topology but also the optimal structural parameters for this topology structure.
A. Uniform Optimization Model
As an example, a FSPMM with two topology structures will be investigated to
illustrate the efciency of the proposed method. The FSPMM is designed to deliver
75 kW output power at a rated speed of 3000 rev/min for a drive machine proposed
in the UTS PHEV [23]. For the optimization, objectives are minimizing cogging
torque (Tcog) and torque ripple (Trip) and maximizing average torque (Tave); and
objective function has the form as
where the subscript initial means the values calculated from the initial design
scheme as shown in Table 7.9, is the efciency, and sf the slot lling factor.
Considering the manufacturing tolerances, the robust optimization model of
(7.4) can be obtained within the framework of design for Six-Sigma (DFSS) shown
in the last chapter, and it has the form as
min : lf x
7:5
s:t: lgi x nrgi x 0; i 1; . . .; 4
where and are the mean and standard deviation of the corresponding terms
estimated by MCA, n is the sigma level and is selected as 6 in this work to
guarantee that the obtained optimal design can achieve the Six-Sigma quality in
industry manufacturing, namely 3.4 defects per million products. To obtain the
mean and standard deviation in this equation, the sample size in the MCA is dened
as 10,000 [3336].
From previous experience, it is found that four parameters including depth of
stator pole, width of rotor pole, height of rotor yoke and width of permanent magnet
are signicant to the motor performance. Therefore, they are selected as the opti-
mization parameters in the optimization. For the robust optimization, 0.05 mm was
chosen to be the manufacturing tolerances for all these parameters, which means
that the standard deviations are 0.05/3 mm. To calculate the reliability of different
motors, a term known as probability of failure (POF) will be calculated as follows
Y
4
POF Probgi 0 7:6
i1
This POF will be used to calculate the defect rate in the 10,000 products used in
the MCA.
B. Optimization Results and Comparison
Figure 7.12 illustrates the comparison of cogging torque curves of 12/10 and 12/14
FSPMMs after robust optimization of (7.5). As shown, the minimal and maximal
cogging torques for 12/10 FSPMM are 9.27 and 8.95 Nm, respectively, while they
are 6.95 and 8.52 Nm, respectively for 12/14 FSPMM. Therefore, 12/14 has
smaller amplitudes of cogging torques.
10
Cogging torque(N*m)
-5
-10
-15
0 90 180 270 360
Rotor position (electrical degree)
232 7 Application-Oriented Design Optimization
Moving torque(N*m)
240
230
220
210
200
0 90 180 270 360
Rotor position (electrical degree)
Figure 7.13 shows the comparison of torque curves of 12/10 and 12/14 FSPMMs
at rated speed after robust optimization. As shown, the minimal and maximal
torques for the 12/10 FSPMM are 215.09 and 236.71 Nm, respectively, and
therefore the average torque is 225.90 Nm, the torque ripple is 21 62 Nm (or
9.57 %). The minimal and maximal torques for the 12/14 FSPMM are 230.18 and
244.89 Nm, respectively. Thus, the average torque is 237.54 Nm, and the torque
ripple is 14.71 Nm (or 6.19 %). It can be seen that the average torque of the optimal
12/14 motor is higher than that of the optimal 12/10 motor, while the torque ripple
is smaller than that of 12/10 motor.
Meanwhile, the reliability and POF of both optimal 12/10 and 12/14 FSPMMs
can be calculated by using (7.6) and MCA. It can be found both motors reliabilities
are almost 100 % after calculation.
In conclusion, 12/14 FSPMM has better performance parameters compared with
those of 12/10 motor in terms of average torque, torque ripple and cogging torque.
It should be noted that only two rotor topologies were investigated in that example.
However, other topology structures can be easily investigated and compared by
using this method, including different stator/rotor pole combinations and different
laminated structures of the steel sheets.
dimension parameters are investigated to acquire the best performance for each
design option.
On the other hand, several remarks could be presented here for this
application-oriented design optimization method. The rst also the most important
one is that only one discipline (electromagnetic analysis) has been investigated for
these examples. However, system-level and integrated design optimizations should
be the starting points of this application-oriented design optimization method,
which means that besides the performance of motor itself, the integrated perfor-
mance of the whole drive system as well as the whole appliance should be
investigated. Therefore, multi-disciplinary design analysis is more important from
the perspective of industry applications, particularly the integration of motor and
control systems, and should be involved for the application-oriented design opti-
mization method.
The second one is that only ve structural parameters (TFM example) and four
parameters (FSPMM example) are investigated in this work. For the
high-dimensional optimization situation, the multi-level (robust) optimization
method presented in Chaps. 46 should be employed to improve the optimization
efciency. Due to the high efciency of the multi-level optimization method, it is
able to efciently optimize and accurately compare the optimal performances
among several electrical machines and drive systems with different topologies and
parameters.
The third one is that only single objective has been considered in these exam-
ples. For the multi-objective optimization situation, the multi-objective sequential
optimization strategy presented in Chap. 4 can be taken to improve the optimization
efciency. Therefore, all the developed design optimization methods as shown in
Chaps. 46 can be employed for the proposed application-oriented design opti-
mization method.
References
1. Zhu ZQ, Howe D (2007) Electrical machines and drives for electric, hybrid, and fuel cell
vehicles. Proc IEEE 95(4):746765
2. Emadi A, Lee YJ, Rajashekara K (2008) Power electronics and motor drives in electric, hybrid
electric, and plug-in hybrid electric vehicles. IEEE Trans Ind Electron 55(6):22372245
3. Lei G, Wang TS, Guo YG, Zhu JG, Wang SH (2014) System level design optimization
method for electrical drive system: deterministic approach. IEEE Trans Ind Electron 61
(12):65916602
4. Salisa AR, Zhang N, Zhu JG (2011) A comparative analysis of fuel economy and emissions
between a conventional HEV and the UTS PHEV. IEEE Trans Veh Techno 60(1):4454
5. Chen JT, Zhu ZQ (2010) Winding congurations and optimal stator and rotor pole
combination of flux-switching PM brushless AC machines. IEEE Trans Energy Convers 25
(2):293302
6. Chen JT, Zhu ZQ, Iwasaki S, Deodhar RP (2011) Influence of slot opening on optimal stator
and rotor pole combination and electromagnetic performance of switched-flux PM brushless
AC machines. IEEE Trans Ind Appl 47(4):16811691
234 7 Application-Oriented Design Optimization
7. Xu W, Zhu JG, Guo YG et al (2011) New axial laminated-structure flux switching permanent
magnet machine with 6/7 poles. IEEE Trans Magn 47(10):28232826
8. Hasanien HM, Abd-Rabou AS, Sakr SM (2010) Design optimization of transverse flux linear
motor for weight reduction and performance improvement using response surface
methodology and genetic algorithms. IEEE Trans Energy Convers 25(3):598605
9. Hasanien HM (2011) Particle swarm design optimization of transverse flux linear motor for
weight reduction and improvement of thrust force. IEEE Trans Ind Electron 58(9):40484056
10. Yao D, Ionel DM (2013) A review of recent developments in electrical machine design
optimization methods with a permanent magnet synchronous motor benchmark study. IEEE
Tran Ind Appl 49(3):12681275
11. Lebensztajn L, Marretto CAR, Costa MC, Coulomb J-L (2004) Kriging: a useful tool for
electromagnetic device optimization. IEEE Trans Magn 40(2):11961199
12. Lei G, Shao KR, Guo YG, Zhu JG, Lavers JD (2008) Sequential optimization method for the
design of electromagnetic device. IEEE Trans Magn 44(11):32173220
13. Lei G, Shao KR, Guo YG, Zhu JG, Lavers JD (2009) Improved Sequential optimization
method for high dimensional electromagnetic optimization problems. IEEE Trans Magn 45
(10):39933996
14. Lei G, Yang GY, Shao KR, Guo YG, Zhu JG, Lavers JD (2010) Electromagnetic device
design based on RBF models and two new sequential optimization strategies. IEEE Trans
Magn 46(8):31813184
15. Lei G, Shao KR, Guo YG, Zhu JG (2012) Multi-objective sequential optimization method for
the design of industrial electromagnetic devices. IEEE Trans Magn 48(11):45384541
16. Lei G, Guo YG, Zhu JG, Chen XM, Xu W (2012) Sequential subspace optimization method
for electromagnetic devices design with orthogonal design technique. IEEE Trans Magn 48
(2):479482
17. Lei G, Xu W, Hu JF, Zhu JG, Guo YG, Shao KR (2014) Multi-level design optimization of a
FSPMM drive system by using SSOM. IEEE Trans Magn 50(2), Article no. 7016904
18. Guo YG, Zhu JG, Watterson PA, Wu W (2006) Development of a PM transverse flux motor
with soft magnetic composite core. IEEE Trans Energy Conver 21(2):426434
19. Zhu JG, Guo YG, Lin ZW et al (2011) Development of PM transverse flux motors with soft
magnetic composite cores. IEEE Trans Magn 47(10):43764383
20. Guo YG, Zhu JG, Dorrell D (2009) Design and analysis of a claw pole PM motor with molded
SMC core. IEEE Trans Magn 45(10):45824585
21. Huang YK, Zhu JG et al (2009) Thermal analysis of high-speed SMC motor based on thermal
network and 3D FEA with rotational core loss included. IEEE Trans Magn 45(106):
46804683
22. Xu W, Lei G, Zhu JG, Guo YG (2012) Theoretical research on new laminated structure flux
switching permanent magnet machine for novel topologic plug-in HEV. IEEE Trans Magn 48
(11):40504053
23. Xu W, Zhu JG, Zhang YC, Wang TS (2011) Electromagnetic design and performance
evaluation on 75 kW axially laminated flux switching permanent magnet machine.
Proc ICEMS, pp 16
24. Cao R, Mi C, Cheng M (2012) Quantitative comparison of flux-switching permanent-magnet
motors with interior permanent magnet motor for EV, HEV and PHEV applications. IEEE
Trans Magn 48(8):23742384
25. Dorrell DG, Knight AM, Evans L, Popescu M (2012) Analysis and design techniques applied
to hybrid vehicle drive machinesassessment of alternative IPM and induction motor
topologies. IEEE Trans Ind Electron 59(10):36903699
26. Takeno M, Chiba A, Hoshi N, Ogasawara S, Takemoto M, Rahman MA (2012) Test results
and torque improvement of the 50-kW switched reluctance motor designed for hybrid electric
vehicles. IEEE Trans Ind Appl 48(4):13271334
27. Zhu ZQ, Chen JT (2010) Advanced flux-switching permanent magnet brushless machines.
IEEE Trans Magn 46(6):14471453
References 235
28. Hua W, Cheng M, Zhu ZQ, Howe D (2008) Analysis and optimization of back EMF
waveform of a flux-switching permanent magnet motor. IEEE Trans Energy Conver 23
(3):727733
29. Fei W, Luk PCK, Shen JX, Wang Y, Jin M (2012) A novel permanent-magnet flux switching
machine with an outer-rotor conguration for in-wheel light traction applications. IEEE Trans
Ind Appl 48(5):14961506
30. Xu W, Zhu JG, Zhang YC, Hu JF (2011) Cogging torque reduction for radially laminated
flux-switching permanent magnet machine with 12/14 poles. In: Proceedings of the 37th
Annual Conference on IEEE Industrial Electronics Society (IECON), pp 35903595
31. Xu W, Lei G, Zhang Y, Zhu JG et al (2011) Development of electrical drive system for the
UTS PHEV. In: Proceedings of the IEEE Energy Conversion Congress and Exposition
(ECCE), pp 18861893
32. Hua W, Cheng M, Zhu ZQ, Howe D (2006) Design of flux-switching permanent magnet
machine considering the limitation of inverter and flux-weakening capability. In: Proceedings
of the 41st IAS Annual MeetingIndustry Applications Conference, vol. 5, pp 24032410
33. Lei G, Guo YG, Zhu JG et al (2012) System level six sigma robust optimisation of a drive
system with PM transverse flux machine. IEEE Trans Magn 48(2):923926
34. Lei G, Zhu JG, Guo YG, Hu JF, Xu W, Shao KR (2013) Robust design optimization of
PM-SMC motors for Six Sigma quality manufacturing. IEEE Trans Magn 49(7):39533956
35. Lei G, Zhu JG, Guo YG, Shao KR, Xu W (2014) Multi-objective sequential design
optimization of PM-SMC motors for six sigma quality manufacturing. IEEE Trans Magn 50(
2):701720
36. Lei G, Wang TS, Guo YG, Zhu JG, Wang SH (2015) System level design optimization
method for electrical drive system: robust approach. IEEE Trans Ind Electron 62(8):4702
4713
Chapter 8
Conclusions and Future Works
Abstract This chapter concludes the book. In summary, this book has focused on
the development of new efcient design optimization methods for novel
high-performance electrical machines and drive systems under two design
approaches, namely the deterministic and robust approaches. These new methods
include sequential optimization method and its multi-objective form, multi-level
optimization method, multi-level Genetic Algorithm, multi-disciplinary design
optimization method and application-oriented system-level design optimization
method. To illustrate the efciency of those proposed methods, several classical test
functions, a TEAM benchmark problem, and four kinds of motors have been
investigated. As shown, the proposed new design optimization methods can achieve
better design objectives for electrical motors and drive systems, such as higher
output power and lower material cost, with much smaller computational cost than
the traditional methods. The proposed robust design optimization approach can
yield optimal designs of electrical drive systems for high quality and high reliability
mass production. Based on these investigations and outcomes, several directions
have been recommended for the future research.
8.1 Conclusions
After a review of the design fundamentals for electrical machines and drive sys-
tems, this book presents several novel efcient design optimization methods in
terms of different optimization situations including:
(1) Low-dimensional situation: The sequential optimization method (SOM) was
presented for this situation. It consists of two processes, the coarse and ne
optimization processes. The main aim of the rst process is to reduce the initial
big design space to a small one by using the space reduction technique. The
main aim of the second process is to nd the optimal solution by using local
sample updating method. From the investigation of a test function, a TEAM
Workshop problem (superconducting magnetic energy storage: SMES) and a
permanent magnet (PM) claw pole motor with soft magnetic composite
(SMC) core, it can be found that SOM and its improved form can present
better optimal solutions while the required computational cost of nite element
analysis (FEA) can be reduced by about 90 %.
(2) Multi-objective situation: The multi-objective sequential optimization method
(MSOM) was proposed for this situation. It uses the strategy of SOM and
improved central composite design technique to reduce the FEA computa-
tional cost. From the investigation of a classic test function and a PM-SMC
transverse flux machine (TFM), it can be seen that the obtained Pareto fronts
are very close to the exact ones.
(3) High-dimensional situation: The multi-level optimization method based on the
sequential subspace optimization strategy and multi-level genetic algorithm
(MLGA) was proposed for this situation. Two popular techniques, the local
sensitivity analysis and the design of experiments (DOE) techniques have been
presented for the sensitivity analysis of design parameters so as to establish the
framework of multi-level optimization method. From the investigation of a
SMES and a surface-mounted permanent magnetic synchronous machine
(SPMSM), it can be found that the multi-level optimization method is efcient
and the obtained optimal solutions are better than those obtained by the tra-
ditional single-level optimization methods.
(4) Multi-disciplinary situation: A multi-disciplinary analysis, design, and opti-
mization framework was presented for PM motors. The multi-disciplinary
analysis includes electromagnetic, thermal, modal and manufacturing analy-
ses. From the investigation of a PM-SMC TFM, it is shown that the obtained
optimal solutions are better than the initial design and those obtained by
non-multi-disciplinary design optimization methods.
(5) Electrical drive systems situation: The system-level design optimization
method based on multi-level optimization strategy and multi-level genetic
algorithm is proposed for this situation. The design parameters in both motor
and control parts have been optimized at the system level rather than the
component level to achieve good steady-state and dynamic performances for
the whole drive systems, including larger output power and higher efciency,
and lower material cost and dynamic overshoot. Two examples have been
investigated to show the efciency of the proposed method. The rst one is a
drive system consisting of a PM-SMC TFM and an improved model predictive
control (MPC) system. The second example is composed of an SPMSM and a
eld oriented control (FOC) system. Through the investigation of these two
drive systems, it is found that both the steady-state and dynamic performances
of the whole drive systems have been greatly improved, and the computational
costs required to obtain the solutions have been reduced signicantly.
(6) Mass production situation: The main concerns in mass production of electrical
machines and drive systems are the product quality and reliability. There are
many unavoidable uncertainties or variations in the industrial manufacturing
process of electrical machines and drive systems, including mainly the
material diversity, manufacturing error, and assembly inaccuracy, which can
result in big variations for the reliability and quality of electrical machines and
8.1 Conclusions 239
(6) Application
While there are many sophisticated new optimization methods, each with its
own advantages and disadvantages, a problem we are facing is to decide
which one would be the best for a particular engineering problem. Therefore, a
selection strategy should be established and a design platform is required to
include all these information. The platform will be a powerful tool for
designing and testing novel high performance electrical drive systems with
new materials, novel topologies, low cost, high efciency, reliability and
robust system performance, and can greatly shorten the design cycles.
Meanwhile, the platform will enable designers to focus their attention on the
system performance which they have the best expertise, rather than the
mathematical algorithms. This will be very valuable for industry mass
production.