EMD504 Assignment 14ETEE034003
EMD504 Assignment 14ETEE034003
i
M. S. Ramaiah University of Applied Sciences
University House, Gnanagangothri Campus, New BEL Road,
M S R Nagar, Bangalore, Karnataka, INDIA - 560 054
Declaration
The assignment submitted herewith is a result of my own investigations and that I have
conformed to the guidelines against plagiarism as laid out in the Student Handbook.
All sections of the text and results, which have been obtained from other sources, are
fully referenced. I understand that cheating and plagiarism constitute a breach of
University regulations and will be dealt with accordingly.
Signature of the
Date
student
Signature of the Module Leader and date Signature of Reviewer and date
iii
Abstract
____________________________________________________________________________
In the recent years the use of the electrical products like the turbines, relays, motors and
generators etc have become quite common in most of the fields, in various industries. The advent of
technological advancements in machine design and operation over the recent years has revolutionized
the industrial sectors starting from small toys and vehicles to bigger electric models in spacecrafts. This
aspect has significantly pushed its performance boundaries beyond the conventional counterparts. In
short based on the wide range of sustained, efficient and cost effective solution offered, they have
garnered the prestigious tag of most reliable machine. Thus first part deals with the debate on the axial
flux permanent magnet machine and its comparison with the conventional machines. Next, the
constructional features of the AFPMSM are dealt and its characteristics are known. Next, the
advantages of the AFPMSMs is learnt and its limitations. Next, the important aspect of the AFPMSMs
that is its latest applications and can it replace the conventional ones with the new application and a
proper justification is done.
The second part deals with the active magnetic bearings. Modeling of the AMB is carried
out in FEM analysis and comparison is made between 3 coils and 4 coils structure such that flux pattern
is analyzed along with tensor forces. Thermal distribution is carried out for both the structures and
finally plots are obtained for various parameters.
The final part deals with the PMSM and its been analyzed in Ansoft Maxwell software.
The analysis is carried out for both the slotless and the slot stator structures. Excitation circuit is drawn
in the Maxwell circuit editor under star and delta connection and flux variations is carried out in both
the cases and finally when the rotor poles are increased, the flux variations are examined in both the
cases and comparison is explained and various parameters are observed.
Declaration Sheet...............................................................................................................................ii
Abstract.............................................................................................................................................iii
Contents............................................................................................................................................iv
List of Tables......................................................................................................................................vi
List of Figures....................................................................................................................................vii
List of Symbols.................................................................................................................................viii
PART-A CHAPTER 1...........................................................................................................................09
Axial Flux Permanent Magnet Synchronous Machines (AFPMSM) in Future.................................09
1.1 Subtitle: Constructional Features of AFPMSM..........................................................................09
2.2 Subtitle: Assign the materials and apply the load conditions...................................................17
3.3 Subtitle: Analyze the flux variations for designed motor at different time instants (star
connected)......................................................................................................................................55
v
3.4 Subtitle: Analyze the flux variations for designed motor at different time instants (delta
connected)......................................................................................................................................57
3.5 Subtitle: Analyze the torque, speed, phase current, and flux linkages.....................................58
3.7 Subtitle: Analyze the flux variations if the poles are increased to six.......................................66
Benefits you have derived by solving this assignment. Whether the assignment was able to assess
module learning outcomes or not?................................................................................................70
vii
List of Symbols
____________________________________________________________________________
viii
PART-A CHAPTER 1
The history of electrical machines shows that the first machines were – more or less – realized in a form
of the axial-flux machine. The first one was invented by Faraday in 1821 and was practically a primitive
permanent-magnet DC machine (Atherton, 1984). Radial-flux machines were invented later and were
patented firstly by Davenport in 1837 (Chan, 1987). Since then radial-flux machines have dominated
excessively the markets of the electrical machines. Manufacturers introduced permanent-magnet
machines successfully during the latest decade. Regardless of the success of radial-flux permanent magnet
machines, axial-flux permanent magnet machines have also been under research interest particularly due
to special-application limited geometrical considerations. A possibility to obtain a very neat axial length for
the machine makes axial-flux machines very attractive into applications in which the axial length of the
machine is a limiting design parameter. Such applications are, for example, electrical vehicles wheel motors
(Profumo et al., 1997) and elevator motors (Hakala et. al., 2000). Axial flux machines have usually been
used in integrated high-torque applications. Axial Flux Permanent Magnet (AFPM) synchronous machines
have been studied with increasing interest since 1990’s, mainly because the AFPM synchronous machine
structure and performance fit very well into the requirements of the electrical vehicles’ wheel motors and
other AFPM motor applications. A permanent magnet design and performance characteristics for axial flux
permanent magnet machine is different when compared to conventional Permanent Magnet Synchronous
Machine (Asko Parviainen 2004).
The use of bearings is essential to all types of machines, in that they provide the function of
supporting another piece or component in a desired position. Two major types include radial and axial
bearings. A further classification can be made into active and passive bearings. Active bearings are
electrically controlled with some type of controller, whereas passive bearings are not electrically powered
and thus have no control mechanism. Magnetic bearings are electromagnetic actuators. These are working
according to electromagnetic principle and it supports moving machinery without physical contact. For
example, they can levitate a rotating shaft and permit relative motion without friction. They are in use
industrial applications such as electric power generation, petroleum refining, machine tool operation and
natural gas pipelines. Active Magnetic Bearing (AMB) is a new kind of bearing. Unlike conventional
bearings, which rely on mechanical forces originating from fluid films or physical contact to support bearing
loads, AMB has advantages of non-contact, friction-less, high speed, low power loss and high accuracy. Due
to these advantages, AMB is used widely in many fields such as transportation, high-speed machine tool,
aerospace and precision instrumentation.
The nonlinear solution of the magnetic vector potential is determined by using the 2-D finite element
method. The force is calculated by Maxwell’s stress tensor method. The electromagnetic field distribution
and density analysis allow verifying the designed AMB and the influence of the shaft and coil current
changes on the bearing parameters (Piłat A 2004).
Finish
Figure 2.3: Flow chart of calculation of flux pattern and force by using ANSYS software for AMB 3 and 4-Coil
structure
To model structure, Select> Peferences> Magnetic-nodal analysis. Since both flux pattern and force
has to be calculated, magnetic-nodal analysis is selected and for thermal distribution thermal analysis is
selected. In preprocessor, Modeling> Create> Lines> Arc –By centre and radius.
Fist, for AMB 3-coil structure, the outer radius is taken as, in global Cartesian> (0,0,0) and (5,0,0) and for
the second circle same procedure in global Cartesian> (0,0,0) and (4.5,0,0), and for the third circle in global
To draw the coils, same procedure is implemented, on stator pole structure for first circle, Create > Lines>
Arc –By centre and radius, in global Cartesian> (0,0,0) and (3.75,0,0) and for the second circle in global
Cartesian> (0,0,0) and (3.25,0,0). And using keypoints on the working plane join the lines such that desired
shapes of coils are obtained. And To draw the circle outside the stator, Create> Lines> Arc –By centre and
radius, in global Cartesian> (0,0,0) and (5.5,0,0). Then complete model is united. Since lines cannot be
overlapped in Modeling> Operate> Boolean> Divide> Area by lines, pick all and pick the lines of the stator
poles and coils. Similarly, same procedure is done for AMB 4-coil structure such that stator poles are 90
degrees apart from each other. And for the different air gaps, the rotor is modeled and in figure 1.2 and
1.3, the rotor is modeled for 0.5cm air gap.
The element type selected is PLANE13. PLANE13 has a 2-D magnetic, thermal, electrical, piezoelectric, and
structural field capability with limited coupling between the fields. PLANE13 is defined by four nodes with
up to four degrees of freedom per node. The element has nonlinear magnetic capability for modeling B-H
curves or permanent magnet demagnetization curves.
For assigning materials,
Materials Permeability Resistivity
(H/m) (Ωm)
1. Air 1 -
2. Coil (Copper) 1 3e-8
3. Stator (Iron) 2000 70e-8
4. Rotor (Aluminum) 4000 -
Table 1.0: selection of Materials (Engineering tool box 2015)
The materials are assigned as per table 1.0 in Preprocessor> Meshing> Mesh tool and meshing is done for
the both 3 and 4-coil structure using smartsize, value to 1 in Element Attributes, set to global and the
whole structure is meshed.
For different air gaps, different excitation current and turns, the following values of flux pattern, flux
density and tensor forces are obtained and are shown in table 1.4.
For AMB 4-coil structure,
Air Flux pattern Flux density Tensor forces
gap(cm),Excitation (Wb) (Tesla) (N)
current(A) and No. of
turns
0.25, 6, 500 0.729e-3 0.001729 0.015526
0.5, 6, 500 0.564e-3 0.001339 0.019438
0.75, 6, 500 0.514e-3 0.001219 0.02377
1, 6, 500 0.496e-3 0.001172 0.03576
Table 1.4: Flux density and forces table for different air gaps
In 3-coil structure, as the area of the air gap increases, flux density decreases and hence there will be few
flux lines pattern. As excitation current increases for example in air gap 0.5 cm, flux will be more and hence
there will be more flux density. As the air gaps are increased, tensor forces will be nearly approximately
constant but not always but depends on the coil structure. As if the coil structures are increased to 4, then
there will be less tensor force compared to 3-coil structure. Thus, a comparative study of three and four
Figure 5.4: Flux density v/s Air gap, 1A, 500 turns
Figure 5.6: Tensor force v/s Air gap, 1A, 500 turns
Figure 6.6: Flux v/s No. of turns, 1A, 0.5 cm air gap
Figure 6.8: Force v/s No. of turns, 1A, 0.5 cm air gap
For the rotor construction, create a circle of radius 2.5 mm and create rectangle of same dimension as
stator pole dimension, with air gap between stator and rotor is of 0.5 mm and rotate the rectangle around
the Z axis with angle of 90 dergrees and are united and the complete rotor is rotated at an angle of 10
degrees.
Stator windings are created in the form of circles on the stator poles and complete model is obtained by
combining both stator and rotor construction.
Next region is created around the model and band is assigned to the rotor model, such that the band
dimension should be more than the rotor model in mm and it is assigned to obtain the motion.
Angle between the two sides of the same coil is 180 degree.
180 degree corresponds to = 180/ 30= 6 slots
The start of the phases must be displaced by 120 degree and so must be finishes. If the start of R or A
phase lies in slot no. 1, the start of Y or B phase must be in slot no. (1 +120/ 30 =) 5 and
that of B or C phase in slot no. (5 +120/ 30 =) 9. This makes the phase sequence RYB or ABC.
Starting of phases: Rs=1 (0 degree) Ys=5 (120 degree) Bs=9 (240 degree).
From the manual calculation mentioned, the following winding parameters are assigned to the command
window and the windings are balanced. Same procedure is implemented when the rotor poles are
increased to 6 and the stator slots are 36 and again with the winding manual calculation, the windings are
balanced.
For the rotor construction, enter the outer and inner dimensions with the number of poles to 4. From the
pole dimension model the thickness of the magnets, magnet duct, width of ribs of the adjacent magnets
and the limited diameter of the four magnets are entered.
4.00
2.00
AirGapFluxDensity [mTesla]
0.00
-2.00
-4.00
-6.00
0.00 125.00 250.00 375.00
ElectricalDegree [deg]
For the slotless model, the flux variation at 5 ms is 2.5616e-2(Wb/m) maximum and at 100ms is 2.3074e-
2(Wb/m) maximum for star connection. This states that as time increases, the flux variation of the pattern
decreases. But in delta connection, the flux variation at 5ms is 2.33e-2(Wb/m) and at 100ms is 2.401e-
2(Wb/m). This states that as time increases, the flux variation also increases for delta connection.
3.4: Analyze the flux variations for designed motor at different time instants (delta connected)
For the slotless model, following results are obtained for the flux variation.
4.00
AirGapFluxDensity [mTesla]
2.00
0.00
-2.00
-4.00
-6.00
0.00 125.00 250.00 375.00
ElectricalDegree [deg]
3.5: Analyze the torque, speed, phase current, and flux linkages
For slotless model, the following results are obtained.
Under star connection:
Torque Maxwell2DDesign1
375.00
Curve Info
Moving1.LoadTorque
Setup1 : Transient
Moving1.Torque
Setup1 : Transient
250.00
125.00
Y1 [NewtonMeter]
0.00
-125.00
-250.00
0.00 100.00 200.00 300.00 400.00 500.00 600.00 700.00 800.00 900.00
Time [ms]
1500.00
1000.00
Moving1.Speed [rpm]
500.00
0.00
-500.00
-1000.00
0.00 100.00 200.00 300.00 400.00 500.00 600.00 700.00 800.00 900.00
Time [ms]
XY Plot 1 Maxwell2DDesign1
1.00 2.00
Curve Info Y Axis
L(L1,L1)
L(L1,L1)
0.90 Setup1 : Transient
FluxLinkage(L1)
Y2
Setup1 : Transient
1.50
0.80 FluxLinkage(L2)
Y2
Setup1 : Transient
0.70 FluxLinkage(L3)
Y2
Setup1 : Transient
1.00
L(L1,L1) [mH]
0.60
Y2 [Wb]
0.50
0.50
0.40
0.30
0.00
0.20
0.10 -0.50
0.00 100.00 200.00 300.00 400.00 500.00 600.00 700.00 800.00 900.00
Time [ms]
8000.00
Moving1.Position [deg]
6000.00
4000.00
2000.00
0.00
0.00 100.00 200.00 300.00 400.00 500.00 600.00 700.00 800.00 900.00
Time [ms]
0.00
-20.00
-40.00
-60.00
For star connection, the maximum torque at 100ms is 260Nm and it remains nearly constant as time
increases. The maximum speed at 100ms is 750 rpm and is nearly remains constant as time increases. The
flux linkage at 100ms is 0.5mH and it nearly remains constant throw ought as the time increases. As the
time increases the position also increases linearly. This shows that the machine is rotating and is changing
the position rotationally. The maximum phase current is 50A and all three phases are 120 degrees apart.
250.00
Moving1.Torque [NewtonMeter]
125.00
0.00
-125.00
-250.00
0.00 100.00 200.00 300.00 400.00 500.00 600.00 700.00 800.00 900.00
Time [ms]
1500.00
1000.00
Moving1.Speed [rpm]
500.00
0.00
-500.00
-1000.00
0.00 100.00 200.00 300.00 400.00 500.00 600.00 700.00 800.00 900.00
Time [ms]
1.00
Y1 [Wb]
0.50
0.00
-0.50
0.00 100.00 200.00 300.00 400.00 500.00 600.00 700.00 800.00 900.00
Time [ms]
8000.00
Moving1.Position [deg]
6000.00
4000.00
2000.00
0.00
0.00 100.00 200.00 300.00 400.00 500.00 600.00 700.00 800.00 900.00
Time [ms]
0.00
-20.00
-40.00
-60.00
40.00
TotalTorque [kNewtonMeter]
30.00
20.00
10.00
0.00
-10.00
0.00 0.20 0.40 0.60 0.80 1.00
RSpeed [rpm]
1.56E+005
1.54E+005
InputLineCurrent [A]
1.52E+005
1.50E+005
1.48E+005
1.46E+005
1.44E+005
-25.00 0.00 25.00 50.00 75.00 100.00 125.00 150.00 175.00 200.00
Ang [deg]
Figure 10.9: Phase current v/s angle plot under star (slot)
XY Plot 1 RMxprtDesign1
6.00
Curve Info
AirGapFluxDensity
Setup1 : Performance
4.00
AirGapFluxDensity [mTesla]
2.00
0.00
-2.00
-4.00
-6.00
0.00 125.00 250.00 375.00
ElectricalDegree [deg]
Figure 11.0: Air gap flux density v/s electrical degree plot under star (slot)
125.00
TotalTorque [kNewtonMeter]
100.00
75.00
50.00
25.00
0.00
-25.00
0.00 0.20 0.40 0.60 0.80 1.00
RSpeed [rpm]
XY Plot 3 RMxprtDesign1
4.80E+005
Curve Info
InputLineCurrent
Setup1 : Performance
4.70E+005
InputLineCurrent [A]
4.60E+005
4.50E+005
4.40E+005
4.30E+005
-25.00 0.00 25.00 50.00 75.00 100.00 125.00 150.00 175.00 200.00
Ang [deg]
Figure 11.2: Phase current v/s angle plot under delta (slot)
4.00
AirGapFluxDensity [mTesla]
2.00
0.00
-2.00
-4.00
-6.00
0.00 125.00 250.00 375.00
ElectricalDegree [deg]
Figure 11.3: Air gap flux density v/s electrical degree plot under delta (slot)
For delta connection, initially at no speed the torque is maximum at 175KNm and as speed increases the
torque decreases and finally becomes zero. The phase current is 4.78e+5 maximum at 180 degrees and is
4.32e+5 minimum at 90 degrees and continues. The air gap flux density is nearly 5mT maximum at 125
electrical degrees and is minimum at 250 electrical degrees. This clearly states that air gap flux density
remains same for both the star and delta connection but torque is maximum in delta connection compared
to star and phase current is maximum in delta compared to star in slot model.
50.00
40.00
StrandedLoss [kW]
30.00
20.00
10.00
0.00
0.00 20.00 40.00 60.00 80.00 100.00
Time [ms]
50.00
40.00
StrandedLoss [kW]
30.00
20.00
10.00
0.00
0.00 100.00 200.00 300.00 400.00 500.00 600.00 700.00 800.00 900.00
Time [ms]
The stranded loss under star connection at 5ms is nearly 55 KW and at 100ms is nearly 15 KW. The torque
created by rotor at 100ms is nearly 260 Nm. The stranded loss under delta connection at 5ms is nearly 50
KW and at 100ms is nearly 18 KW. The torque created by rotor at 100ms is nearly 260 Nm.
In slot model, the total torque created by rotor under star connection is nearly 58 KNm and the total
torque created by rotor under delta connection is nearly 175 KNm.
3.7: Analyze the flux variations if the poles are increased to six
In slotless model, when the rotor poles are increased to six the following results are obtained.
In slotless model, for 4 poles the flux variation at 5 ms is 2.5616e-2(Wb/m) maximum and at 100ms is
2.3074e-2(Wb/m) maximum and when the poles are increased to 6, the flux variation at 5 ms is 2.313e-
3(Wb/m) maximum and at 100ms is 1.6502e-3(Wb/m) maximum. This clearly states that flux variation is
more in 4 poles as compared to 6 poles and as the time increases the flux decreases in both the cases.
25.00
AirGapFluxDensity [uTesla]
0.00
-25.00
-50.00
0.00 125.00 250.00 375.00
ElectricalDegree [deg]
For the slot model, in the 4 poles structure, the maximum air gap density at 125 electrical degrees is nearly
5mT and it remains constant up to 200 degrees and afterwards flux density becomes negative for star
connection and for delta connection also it remains the same. Thus, for both the connection the air gap flux
density remains same. But as compared to 6 poles structure the air gap flux density decreases i.e, the
maximum air gap density at 125 electrical degrees is nearly 45 µT and it remains constant up to 200
degrees and afterwards flux density becomes negative for both star and delta connection. Thus, this clearly
states that 4 poles structure is having higher air gap flux density as compared to 6 poles structure.
Solving the assignment has provided with a lot of beneficiaries and it has helped a lot to access
module learning outcomes and the entire syllabus has been covered including that of the laboratory
aspects. In the part A of this assignment learnt about the features and characteristics of AFPMSMs
compared with the conventional ones. Learnt about the advantages and limitations of the axial machines.
More importantly learnt the various applications of the axial machines and how it has replaced the
conventional ones.
In the second part helped in understanding the applications of the motor used in today’s
technology and learnt about Ansys FEM software tool. In the final part learnt about the Ansys Maxwell
software and came across the application type problem of what kind of motor to be installed and learnt the
importance of both the slotless and slotted structure and its observations and learnt about the magnetic
bearings and its observation is carried out in ANSYS software.
Thus solving this assignment was not only the part of the academics but it really gave a good
knowledge of the industrial electrical products and protocols and the tools in order to create 2D and 3D
models .It was like fun solving the questions, modeling and analysis and it also helped in accessing all the
module syllabus and was of a good benefit.
A. Anonymous, (2015), ‘List of Properties with values of several known metals’, Engineering tool box.
B. Anonymous, (2015), ‘Axial flux permanent magnet machine applications’, Elektromaschinen Und
Atriebe.
C. Anonymous, (2015), ‘Axial flux permanent magnet motors in swimming pool pump’, Fasco
Australia.
D. Asko Parviainen, (2004), ‘Performance of axial flux permanent magnet synchronous low speed
machines’, IEEE Transaction on Industry Applications, Vol. 40, No. 5, 2004, pp. 1333-1340.
E. Adrian Augustin Pop, (2013), ‘Axial-flux vs. radial-flux permanent-magnet synchronous generators
for micro wind turbine application’, EPE'13 ECCE Europe, ISBN: 978-90-75815-17-7.
F. A. Mahmoudi, (2012), ‘Axial flux permanent magnet motor design and applications’, Progress in
Electromagnetic Research, Vol. 122, pp.467–496.
G. Dean J Patterson, (2015), ‘A Comparison of Radial and Axial Flux Structures in Electrical Machines’,
The Timken Company, Mail Code TEC-05, University of Nebraska, Lincoln, NE 68588-0511.
H. Kim J.S. and Sul S.K, (1997), ‘New approach for high-performance PMSM drives without rotational
position sensors’, IEEE Transactions on Power Electronics, Vol. PE-12, Issue 5, pp 904−911.
I. Piłat A, (2004), ‘FEMLab software applied to Active Magnetic Bearing analysis’, International
Journal of Applied Mathematics and Computer Science, Spec. iss.: Issues in modeling, optimization
and controls, Vol. 14, 4: 497-501.
J. Pabitra Kumar Biswas and Subrata Banerjee, (2013), ‘ANSYS Based FEM Analysis for three and four
coil Active Magnetic Bearing-a Comparative Study’, International Journal of Applied Science and
Engineering, Vol.11, 3: 277-292.
K. R.W.Powell, (2001), ‘Thermodynamics and thermo physical properties of materials’, National
standard reference data series.
L. S.Geetha, (1993), ‘Axial flux permanent magnet servo motors with sixteen poles’, University of
Wollongong thesis collection, Vol.12, 1339.
.
Figure 1.0: Calculated and measured cogging torque (Asko Parviainen 2004)
Figure 1.1: Machine efficiency as the function of output power (Asko Parviainen 2004)
Figure 1.8: Power consumption against speed (Elektromaschinen Und Atriebe 2015)
Figure 2.1: Pump assemblies, with the original single phase induction machine version on left, PM axial
flux version on right (Fasco Australia 2015)