Open access peer-reviewed chapter

Mechanical Design of Rotors with Surface Mounted Permanent Magnets

By Simon Barrans and Levi Mallin

Submitted: February 9th 2018Reviewed: May 18th 2018Published: November 5th 2018

DOI: 10.5772/intechopen.78854

Downloaded: 297

Abstract

Electric machines with permanent magnet rotors are becoming increasingly popular due to the high power density that they offer relative to other configurations. Where the speed of rotation is high, the magnets are typically mounted on the surface of the rotor and retained by an outer sleeve. In the literature, a variety of analytical models have been proposed to aid the mechanical design process. Many of these models contain inherent assumptions about the stress and strain field in the rotor which may not always be apparent. In this article the range of rotor stress models are presented and explained and the limitations due to the inherent assumptions are investigated. This will allow the designer of such rotors to assess mechanical performance without introducing unforeseen errors.

Keywords

  • plane stress
  • plane strain
  • generalized plane strain
  • surface mounted permanent magnets
  • rotor stress
  • rotor failure
  • rotor sleeve

1. Introduction

High speed electric machines are increasingly being used in place of mechanical drive systems. As [1, 2, 3] observe, they are more efficient and reliable than their mechanical equivalents and have high power densities. Gerada et al. [4] observe that this makes them an attractive choice across a number of areas including flywheel energy storage and high speed spindles. A particularly demanding application is in electrically assisted turbochargers [5].

In order to achieve the increasingly high speeds being demanded by a number of applications, synchronous machines with permanent magnet rotors are typically used. These machines can be broadly categorized as having interior permanent magnets (IPM’s) embedded within the iron core of the rotor or surface-mounted permanent magnets (SPM) attached to the outer surface of the rotor. When considering the mechanical design of these rotors, SPM rotors allow higher rotational speeds than IPM rotors and will be the focus of this article.

The materials used for the permanent magnets include samarium cobalt (SmCo5 or Sm2Co17) and neodymium (Nd2Fe14B). These materials are formed from sintering or bonding and are typically very weak in tension. If the magnets were simply glued to the outer surface of the rotor, the tensile stresses induced by rotation would cause them to fail. They are therefore typically retained onto the rotor shaft or core with a sleeve, as shown in Figure 1. The interference fit between the sleeve and the magnets places the magnets into compression. A sufficiently large interference will ensure that this initial compressive state will prevent tensile stresses being generated under rotation. However, a large interference can induce large circumferential forces in the sleeve. Increasing the sleeve thickness increases the effective air gap between the rotor and the stator which reduces electromagnetic efficiency. The mechanical design of SPM rotors is therefore challenging and has been approached in a number of ways using both analytical and numerical approaches: The analytical approaches treat the sleeve as either a thin walled or thick walled cylinder. Where a thick walled cylinder approach is used, it may be approximated as being in plane stress, plane strain or generalized plane strain. All of these approaches incorporate some degree of approximation which may limit their applicability to a particular rotor design. This article will explore these approximations and limitations.

Figure 1.

SPM rotor assembly.

2. Classical analytical methods

2.1. Thin shell approach

One of the simplest approaches to SPM rotor mechanical design is to assume that the sleeve acts as a thin walled cylinder, wall thickness tand radius R, as proposed in [6] and illustrated in Figure 2. When the rotor is stationary, the radial interference, δ, between the sleeve and the magnets generates a circumferential strain, εθand stress, σθ, in the rotor and an interference pressure, S, between the sleeve and the magnets, given respectively by Eqs. (1)(3).

εθ=δRE1
σθ=Eεθ=EδRE2
S=σθtR=EδtR2E3

Figure 2.

Thin cylinder theory.

where Eis the Young’s modulus of the sleeve material. It is important to note here that all the deformation required to assemble the sleeve onto the magnets is assumed to take place in the sleeve, i.e. the magnets, back iron (if present) and shaft are assumed to be rigid.

The ‘thin walled’ cylinder approximation used here assumes that there is no variation in circumferential stress and strain through the thickness of the sleeve. Hence, in Eq. (1)(3) only a single value for these factors is generated. The rather misleadingly termed ‘thick walled’ cylinder approach correctly allows for variation in stress and strain through the sleeve wall. The ‘thin walled’ model is therefore an approximation to the ‘thick walled’ model and for any sleeve thickness, will introduce an error in the calculation of interference pressure. The various forms of the thick walled cylinder approach will be discussed in the following sections. Taking either the plane stress approach or the generalized plane strain approach gives the relationship between interference pressure and interference as:

S=δERo2RC2RCRO21+ν+RC21νE4

where ROand RCare the outer and inner radii of the sleeve respectively and νis the Poisson’s ratio of the sleeve material and the magnets, back iron and shaft are assumed to be rigid.

Assuming that the radius of the ‘thin walled’ sleeve, R, is the mean of the outer and inner radii:

RC=Rt2andRO=R+t2E5

Taking typical properties for an Inconel sleeve (E=205GPa, ν=0.294) and δ=0.1mm, the error in the shrinkage pressure introduced by the ‘thin walled’ cylinder approximation is as shown in Figure 3. For small sleeve thicknesses, the errors are small. However, as sleeve thickness increases and sleeve diameter reduces, the thin cylinder approach significantly underestimates the circumferential stress.

Figure 3.

Thin cylinder approximation error.

Having established the interference pressure acting between the sleeve and the magnets when the rotor is stationary, the effect at maximum rotor speed must be determined. Various approaches can be taken with different criteria being used to assess the mechanical integrity of the rotor. The method proposed in [6] assumes that the interference pressure between the magnets and the back iron is the same as that between the sleeve and the magnets and then determines the centripetal force required to reduce this pressure to zero (indicating that the magnets may then slip on the back iron). This approach introduces two further approximations: Firstly, to maintain force equilibrium in the radial direction, the inner surface area of the magnet must be identical to the outer surface area. Secondly, the magnets must have a sufficiently small circumferential dimension (i.e. be segmented) for no significant circumferential stress and no bending to be generated during rotation. The impact of the first effect can be readily determined and indeed can be corrected for. The impact of the second effect is more difficult to determine. If the permanent magnet is a single cylinder, it will resist circumferential expansion during rotation and there will be a smaller reduction in the interference pressure between the magnets and the back iron. However, if the magnet is formed from segments which abut in the circumferential direction, the centripetal forces will tend to separate them. This will lead to a substantial increase in the stresses in the sleeve at the point of separation and the generation of tensile stresses due to bending in the magnets. Whilst the magnet segments remain in full contact in the circumferential direction, their behavior will be identical to that of a single cylinder. This article will be restricted to rotors with magnets which are either single cylinders or are formed from segments which abut and remain in full contact during rotation.

2.2. Thick walled cylinder or disk analysis

All ‘thick walled’ approaches to the mechanical analysis of rotors are founded on the fundamental principles of force equilibrium, compatibility and the constitutive relationship between stress and strain within the material. These relationships are initially applied to a two dimensional cross section of the rotor, perpendicular to the axis and are presented in many text books (see for example [7]). This theoretical material is repeated here to allow the reader to appreciate the differences between the different analysis approaches.

Figure 4 shows a segment, angle and thickness dr, of an annular ring, radius r, taken from a cylinder or disk which may be rotating at an angular velocity, ω, and may be subjected to either internal or external pressures. The disk material has a density, ρ. It is assumed that the radial stress, σr, and circumferential stress vary with radius.

Figure 4.

Internal stresses in the cylinder/disc.

Considering equilibrium of forces on this segment in the radial direction, the following relationship is obtained:

σθσrrrdr=ρω2r2E6

The constitutive relationship can then be expressed as Eqs. (7)(9) where εris the radial strain, εzthe axial strain, αthe coefficient of thermal expansion of the material and ΔTthe change in temperature imposed on the material.

εr=dsdr=1Eσrνσθνσz+ΔTαE7
εθ=sr=1Eσθνσrνσz+ΔTαE8
εz=1Eσzνσrνσθ+ΔTαE9

In Eqs. (7) and (8), sis the change in radius of the segment and dsis the change in the radial dimension, as indicated in Figure 5.

Figure 5.

Change in radial dimensions of the segment.

To maintain compatibility:

drεθdr=dsdr=εrE10

or, with reference to Eqs. 7 and 8:

σrσθ1+ν=rdσθdrνdσrdr+dσzdr+EαrdΔTdrE11

As Li et al. [8] determined, the thermal barrier resisting heat transfer from the outer surface of the rotor to the air gap is much larger than the thermal resistance within the rotor components. Hence, temperatures are essentially uniform with respect to radius and vary only slightly along the axis. Assuming that dΔT/dr=0and combining Eq. 11 with the equilibrium equation (6), gives:

dσrdr+dσθdr+dσzdr=1+νρrω2E12

It is worth emphasizing that Eq. 12 is common to all the analytical approaches presented. However, further development of the analysis requires additional information and the solution paths diverge.

2.2.1. Plane stress approach

An approach previously used in some work [9, 10, 11, 12] has been to assume that the axial cross section is in a state of plane stress, where σz=0. Introducing this condition into Eq. (12) gives:

dσrdr+dσθdr=1+νρrω2E13

which can be integrated to give:

σθ+σr=1+νρω2r22+2AE14

where Ais a constant of integration. Reintroducing the equilibrium equation, 6, and further integrating leads to:

σr=ABr23+ν8ρω2r2E15
σθ=A+Br21+3ν8ρω2r2E16

with Bbeing a further constant of integration.

In [10, 11, 12] the plane stress assumption is justified by the fact that there is no significant axial stress in the rotor as there is nothing to constrain axial extension or compression. However, as will be demonstrated later, the lack of an external force does not prevent an axial stress being generated in the rotor under rotation. The only valid justification for the plane stress assumption is that the rotor is in the form of a disk which is of insufficient thickness for axial stress to develop.

Where the rotor is modeled as a pair of interfering cylinders, there are two sources of stress: stresses due to the interference and stresses due to rotation. Assuming linear behavior, these stresses can be determined separately and then simply summed to determine the total stress. There are also two load cases to consider:

  • Stationary rotor with stress due only to interference.

  • Rotating rotor with stresses generated by both rotation and interference.

For the second load case it is important to account for the change in interference due to rotation [13].

Considering Figure 6, the two interfering cylinders must achieve a common radius, Rc, and an interference or shrinkage pressure, S, is generated. Treating the two cylinders separately, the boundary conditions are:

Inner cylinder:Atr=Ri;σrii=0
r=Rc;σric=S
Outer cylinder:Atr=Ro;σroo=0
r=Rc;σroc=S

Figure 6.

Interference condition.

where σriiis the radial stress in the inner cylinder at the inner radius, σricis the radial stress in the inner cylinder at the common radius, σroois the radial stress in the outer cylinder at the outer radius and σrocis the radial stress in the outer cylinder at the common radius.

Substituting these boundary conditions into Eq. (15) allows the constants of integration to be found and Eqs. (15) and (16) become:

Inner cylinder:σri=SRc2Ri2Rc21Ri2r2E17
σθi=SRc2Ri2Rc21+Ri2r2E18
Outer cylinder:σro=SRc2Ro2Rc21Ro2r2E19
σθo=SRc2Ro2Rc21+Ro2r2E20

The compatibility requirement that the two cylinders establish a common radius can be expressed as:

Rci1+εθics=Rco1+εθocsE21

where εθicis the circumferential strain in the inner cylinder at the common radius, εθocis the circumferential strain in the outer cylinder at the common radius, Rciis the outer radius of the inner cylinder prior to shrinkage and Rcois the inner radius of the outer cylinder prior to shrinkage.

Noting that σz=0, Eqs. (17)(20) can be used with Eq. (8) to determine the strain terms in Eq. (21) in terms of the dimensions of the cylinders and the shrinkage pressure. The shrinkage pressure can then be determined from:

S=RciRco+ΔT[RciαiRcoαoRcoRo21+νo+Rc21νoEoRo2Rc2+RciRc21νi+Ri21+νiEiRc2Ri2E22

where the subscripts iand oon the material properties E, νand αindicate different properties for the inner and outer cylinders respectively.

Having determined the shrinkage pressure, the stresses due to shrinkage can be determined from Eqs. (17)(20).

To determine the stresses due to rotation, Eqs. (15) and (16) can again be used, with the boundary conditions:

Inner cylinder:Atr=Ri;σrii=0
r=Rc;σric=0
Outer cylinder:Atr=Ro;σroo=0
r=Rc;σroc=0

being used to determine the constants of integration.

2.2.2. Plane strain approach

Perhaps the first application of plane strain to the analysis of electric machine rotors is in [6] with reference to this paper being made by [14, 15, 16]. In this approach it is assumed that the axial strain is zero and therefore also that z/dr=0. Hence, differentiating Eq. (9) gives:

1Edσzdrνdσrdr+dσθdr+αdΔTdr=0E23

As before, assuming dΔT/dr=0:

dσzdr=νdσrdr+dσθdrE24

and introducing this into Eq. (12) gives:

dσrdr+dσθdr1ν=1+νρrω2E25

As with the plane stress case, this equation is integrated, the equilibrium condition (Eq. (6)) reintroduced and a further integration leads to:

σr=ABr232ν81νρω2r2E26
σθ=A+Br21+2ν81νρω2r2E27

For the case of no rotation, Eqs. (26) and (27) are identical to Eqs. (15) and (16) so Eqs. (17)(20) can be used to determine the stresses. However, with rotation present, the stresses and strains will not be the same. In order to determine the circumferential strain due to rotation (and hence the change in interference), the axial stress, σz, is required in Eq. (8). To determine this stress, the condition of plane strain is imposed on Eq. (9) to give:

σz=νσr+σθEE28

To determine the shrinkage pressure, Eqs. (26)(28) are substituted into Eq. (8) for each cylinder at the common radius and then Eq. (21) is used to enforce compatibility between the cylinders. This leads to:

S=RciRco+ΔTαi1+νiRciαo1+νoRcoRcoRc21νo2νo2+Ro21+νoEoRo2Rc2+RciRc21νi2νi2+Ri21+νiEiRc2Ri2E29

The analysis methodology is otherwise identical to that for the plane stress case with stresses due to interference and rotation being found individually and then summed. Again, as [13] observe, it is important to account for the change in interference when analyzing the combined rotation and interference load case. The impact of this can be seen in Figure 7 for a cylinder pair with the properties given in Table 1 and an initial interference of 0.03 mm. It can be seen that if the interference is not updated to account for rotation, the compressive stress at the outer radius of the magnet, σθic, is overestimated and similarly, the von Mises stress on the inner radius of the sleeve, σVMoc, is overestimated as speed increases.

Figure 7.

Impact of not updating interference.

Inner radius, Ri6.35 mm
Common radius, Rc11.4 mm
Outer radius, Ro13 mm
Inner Poisson’s ratio, νi0.27
Outer Poisson’s ratio, νo0.284
Inner density, ρi8400 kg/m3
Outer density, ρo8220 kg/m3

Table 1.

Properties for axial stress study.

A criticism of the plane strain approach is the requirement for axial strain to be zero. Although a rotor assembly maybe mounted on a solid shaft with end caps to retain the magnet and sleeve, this shaft will be a deformable body. Hence, it will be possible for the rotor components to change in length.

3. Generalized plane strain approach

3.1. Axial stress generation

The generalized plane strain approach is similar to the plain strain approach and Eq. (24) can be applied. Using this in conjunction with Eq. (12), the variation of axial stress with radius is given by:

dσzdr=rνiρiω21νiE30

Integrating gives the axial stress as:

σz=νiρiω221νiCr2E31

where Cis a constant of integration.

Assuming that an external axial force acts on the cylinders, as indicated in Figure 6, and taking the inner cylinder as an example, there must be equilibrium between the internal and external axial forces, expressed as:

Fzi=2πRiRcrσzdrE32

where Fziis the external axial force acting on the inner cylinder. Carrying out the integration allows the constant Cto be determined and the axial stress is then given by:

σzi=FziπRc2Ri2+νiρiω241νiRc2+Ri22r2E33

Similarly, for the outer cylinder:

σzo=FzoπRo2Rc2+νoρoω241νoRc2+Ro22r2E34

where Fzois the external axial force acting on the outer cylinder.

As an aside at this point, it is worth considering the possible impact of axial stress generation. Taking a cylinder pair with the dimensions and material properties shown in Table 1 and assuming that no axial force is applied (Fzi=Fzo=0), the axial stress due to rotation alone, over a range of speeds is shown in Figure 8. With no rotation applied, no axial stress is generated and the generalized plane strain approach is then identical to the plane stress approach. However, once rotation is introduced, axial stress is generated with peak values at the inner and outer radii. These are the positions where the circumferential strain needs to be calculated, according to Eq. (8), in order to determine the degree of interference and shrinkage pressure. The axial stress will therefore impact on these critical values.

Figure 8.

Axial stress generation.

A further point worth noting is that at the end of the cylinder axial stress cannot exist as this would place the material out of equilibrium. On the end face, a state of plane stress must exist. This change in stress state may be responsible for the end failure noted by [17].

3.2. Combining the cylinders

To account for the interaction between the cylinders, two equilibrium and two compatibility conditions need to be considered:

  1. A. In the radial direction there must be force equilibrium at the interface between the cylinders.

  2. B. At the interface there must be circumferential strain compatibility. With reference to Figure 6, this is expressed as Eq. (21).

  3. C. Assuming that there is sufficient friction acting between the two cylinders to prevent a differential change in length, there must be axial strain compatibility. This change in length is shown in Figure 6 as εzLwhere Lis the rotor length. This compatibility condition therefore becomes:

εzi=εzoE35

Where εziand εzoare the axial strains of the inner and outer cylinders respectively.

  1. D. There must be equilibrium of forces in the axial direction between the two cylinders. With reference to Figure 6:

Fzi=FzoE36

To enforce condition A, the constants in Eqs. (26) and (27) are determined with the shrinkage pressure, S, applied at the common radius for the inner and outer cylinders. A pair of simultaneous equations are then determined by considering conditions D and C together and then D and B together. Details of this analysis are given in [18], the result being:

UVXYFzoS=WZE37

where:

U=1πEiRi2Rc2EoRo2Rc2
V=2Rc2νOEiRi2Rc2νiEoRo2Rc2
W=Ro2Rc2Ri2Rc2ω22νOEiρoRc2+Ro2νiEoρiRc2+Ri2+TαiαoEoEi
X=1πνiEoRciRo2Rc2+νoEiRcoRc2Ri2
Y=EoRciRo2Rc2Rc21νi+Ri21+νi+RcoEiRc2Ri2Ro21+νo+Rc21νo}
Z=Rc2Ri2Ro2Rc2{EoEiRcoRci+TαoRcoαiRci
+ρoRcoω2Ei81νo2Rc2+3Ro24νoRo2+Rc22νo2Ro2Rc2
ρiRciω2Eo81νi2Rc2+3Ri24νiRc2+Ri22νi2Ri2Rc2}

4. Comparison of models

4.1. Basic model parameters

The material properties used in the comparison of the various models are as shown in Table 1, these being typical properties for samarium cobalt magnets and an Inconel sleeve. Analyses were carried out for the parameters shown in Table 2.

ParameterMinimumMaximum (or discrete values)
Magnet inner radius, Ri, (mm)4, 6, 8
Common radius, Rc, (mm)14
Interference, δ, (mm)00.1
Sleeve thickness (mm)010
Rotor speed (ω) (krpm)0100

Table 2.

Geometric parameters for comparative study.

For this study temperature was not treated as a variable so thermal expansion effects were ignored.

A critical result is the circumferential stress on the outer radius of the magnet, σθic, at maximum speed. This is shown in Figure 9 for a magnet inner radius of 6 mm. When the degree of interference and sleeve thickness are low there is insufficient compression generated in the magnets by the sleeve and this stress becomes tensile. This may lead to magnet failure dependent on the precise material used. Or, if magnets are segmented, it would lead to magnet separation and localized over-loading of the sleeve.

Figure 9.

Circumferential stress on the magnet outer radius at speed.

Figure 10 shows the von Mises stress generated at the inner surface of the sleeve, σVMoc, at zero speed, again with a magnet inner radius of 6 mm. This stress or its equivalent at maximum speed will be critical for the sleeve which is normally a ductile material. This stress becomes larger as interference increases (and hence both radial and circumferential stress magnitudes increase) but reduces as sleeve thickness increases.

Figure 10.

Von Mises stress at the inner radius of the sleeve at speed.

4.2. Comparing plane stress and plane strain models

Figure 11 shows the percentage difference in von Mises stress between the plane stress and plane strain theories at maximum speed for a rotor with a magnet inner radius of 4 mm. For smaller levels of interference, this difference is substantial. The percentage difference in circumferential stress in the magnet for the same case is shown in Figure 12. Where this stress approaches zero, the results approach infinity and have been omitted from the figure for clarity. In order to optimize rotor design it will be desirable to approach zero circumferential stress in the magnet as maximum speed is reached. Hence, whilst the absolute differences between the two theories are small (because the absolute values are small), they are still critical. These results highlight the need for a theoretical approach that does not incorporate the conflicting approximations of the plane stress and plane strain models.

Figure 11.

Percentage difference in σ VMoc at maximum speed.

Figure 12.

Percentage difference in σ θic at maximum speed.

4.3. Generalized plane strain model

The generalized plane strain model does not include any of the approximations of the plane stress or plane strain theories. For the case of relatively long cylinders with high levels of friction between the components, this theory can be regarded as accurate away from the ends of the rotor. Figure 13 shows a selection of results comparing the plane stress and plane strain models to the generalized plane strain model for a rotor with a magnet inner radius of 4 mm. When stationary the error in σθicincreases with decreasing sleeve thickness for both the plane stress and plane strain models. This is much more significant for the plane stress model. Looking at the von Mises stress at maximum speed, the plane strain model shows increasing error as the interference is reduced but relatively little dependence on sleeve thickness. The error in the plane stress result is dependent on both sleeve thickness and interference but is never as large as the peak error value for the plane strain theory.

Figure 13.

Errors in plane stress and plane strain theories.

5. Conclusions and further work

The work presented here has shown that:

  • Using a thin walled cylinder model for the rotor sleeve will introduce increasing errors as the sleeve thickness increases.

  • Neglecting the change in interference due to rotation will introduce significant errors.

  • The plane stress and plane strain models will often give similar results for rotors operating well within their design limits. However, when rotors are optimized to minimize sleeve thickness or increase speed such that the circumferential stress in the magnet becomes a limiting factor, there is significant difference between the theories. Both theories incorporate approximations which are not generally true for this application.

  • The plane stress and plane strain models show significant differences from the generalized plane strain model, particularly at the design limits.

Despite avoiding the limiting approximations of the other theories, the generalized plane strain model still has limitations. Further development is required to consider:

  • The impact of lower levels of friction where the cylinders may slip against each other in the axial direction.

  • Situations where more than two cylinders are interfering such as with magnets affixed to a back iron and retained with a sleeve.

  • Sleeves made from composite materials where the constitutive equations, 7 to 9, are not applicable in the form given here.

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Simon Barrans and Levi Mallin (November 5th 2018). Mechanical Design of Rotors with Surface Mounted Permanent Magnets, Electric Machines for Smart Grids Applications - Design, Simulation and Control, Adel El-Shahat, IntechOpen, DOI: 10.5772/intechopen.78854. Available from:

chapter statistics

297total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Designing a DC Motor Simulator Based on Virtual Instrumentation

By Nicolae Patrascoiu

Related Book

First chapter

Introductory Chapter: Recent Advances

By Adel El-Shahat

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us