Computer Program
Aids Magnetics Design

of Brush And Brushless DC Motors

George P. Gogue

Dan Jones

*A program that operates on a
personal computer provides design data based on the equivalent magnetic
circuit of a brush or brushless, permanent magnet DC motor.*

A computer program that can run on a personal computer employs simple lumped reluctance simulation to determine the magnetics design. The program has sufficient accuracy for the majority of motor design applications and can be used in place of the slower but more accurate finite element mesh (FEM) simulation programs.

To demonstrate use of this program, we will simulate a set of brush and brushless DC motor designs for a rare earth magnet motor. Our intent is to show that a simpler program can be used as long as the user understands that the resulting boundary conditions will impact simulation accuracies.

Whether you use computer-aided or calculator-aided design techniques, permanent magnet DC motor design involves many interactive design conditions. Therefore, the designer must:

* Determine the flux from the main field (magnets) and the winding

* Know the resultant magnetic field and the armature reaction

* Maintain a continuous check on the iron losses in the different sections of the magnetic circuit, as well as the copper winding configuration and fill.

The result of these optimization activities is a magnetic design that integrates both electrical and mechanical performance. This design would define the geometry of all parts of the magnetic circuit and also the materials involved.

An electrical equivalent of the magnetic circuit (Refs. 1, 2 and 3) enables development of the performance equations.

Due to the complicated nature of the magnetic circuits, various assumptions are made to simplify the equations. These assumptions either result in a very small error or limit the usefulness of the equations to a certain range with some of the parameters. The computer program used in this work is based on the reluctance model of Figure 1:

where l is the length of the magnetic path in meters and A is the cross sectional area in square meters. The reluctance is, therefore, in Henry/m and its reciprocal is the permeance. The permeability of the material of the particular section of the magnetic circuit is . But is a parameter totally dependent on the reluctance of the circuit. Perhaps Figure 2 is a helpful representation of the interaction between those parameters.

The other basic equations are:

where mmf is the magneto-motive
force in Ampere-turns and* *flux, , is in Webers. This equation
is the magnetic analog of the electric Ohm's law equation:

where B is the flux-density is Tesla.

The computer program sums up the reluctances of all the different parts of the circuit and calculates the flux corresponding to a particular value of total reluctance. But because of the interactive nature of flux and reluctance, this calculation is repeated many times in an iterative fashion. The limit on this calculation is a residue of percentage error that is pre-selected by the designer.

The reluctance terms of Figure 1 correspond to different parts of the motor relevant to the magnetic circuit. Tooth reluctance is perhaps the most difficult to determine accurately, because of the wide variation in flux-density of the different parts of each tooth. Tooth tips are the most likely parts of the tooth to reach saturation because of their narrow cross section. Unchecked flux-density in the tooth tips can lead that part into saturation, making the slot opening effectively larger than the physical opening. This has the disadvantage of reducing the effective area of the steel facing the magnets. The reluctance of the path increases with the reduced area of steel, leading to a drop in the air-gap flux. To avoid this, the tooth tips are made as thick as reasonably possible.

Other sections of the tooth also experience different levels of flux-density, although most teeth are parallel sided. This is the result of varying amounts of leakage flux outside the tooth and along its radial surfaces. For added accuracy of the reluctance calculation, the tooth is divided into sections that have differing leakage factors. In addition to the tooth tip section, this division, therefore, defines the conditions in the tooth sufficiently accurately.

Because there are usually several teeth facing each magnet, it is important to determine the minimum number of teeth including fractions of teeth that constitute the flux-carrying steel.

The core is the steel part behind
the teeth that completes the magnetic circuit between adjacent magnets.
The reluctance of this path is , and is represented by one value in Figure 1*.
*In reality, the reluctance of the core depends on the flux-density
(and )
of the particular segment of the core. To simplify the problem, the
computer program divides up the core into a number of segments
depending on the number of teeth chosen in the design. For each segment
the flux-density is assumed constant, allowing calculation of and the
reluctance. The reluctance of the core is then the sum of the
reluctances of all the parts.

Usually a carbon-steel cylinder around the magnets, the back-iron (or housing) forms the return path of flux between adjacent magnets. Flux-density levels are different at the different parts of the iron. Values are highest at interpolar zones and least along the centerline of the magnet.

In a similar fashion to the procedure followed in the stator core, the back-iron is divided into sections that are individually assumed to have a constant flux-density. This enables the calculations of and the individual reluctances leading to the total reluctance of the back-iron.

The air gap referred to here
includes the length of both the physical air gap and the actual radial
thickness of the magnet. The implied assumption here is that most
magnet materials used in DC motors have a permeability equal to that of
air. The combined reluctance of the air gap and magnet are represented
by in
Figure* *1. The magnet arc is usually smaller than the full pitch
(unless a ring magnet is used, which results in a near full-pitch arc).
In addition, the axial length of the magnet is usually longer than that
of the stack. The combination of the above two features results in the
cross sectional area of the magnet being different from that of the air
gap. Besides this difference in area, an account must also be taken of
the flux fringing at the edges of the magnet. The most traditional way
is to use Carter's coefficients that approximate the effect of fringing
to an additional air gap length. An alternate technique of accounting
for fringing flux is discussed by Ireland (Ref. 6).

Reluctance of the air gap represents by far the largest portion of the total reluctance of the circuit. Special care, therefore, must be given to the assumptions and approximations that can result in a large percentage error in this calculation. For example, assuming a flat 5 or 10% increase in air gap area is only a first order approximation based on the designer's experience.

**Magnetic
Characteristics**

With the reluctance values of all the different parts of the motor represented as equations, calculations now center around the amount of flux flowing in the circuit, as shown in Figure 1. The mmf referred to in Equation 2 is the sum of Ampere-turns from the main magnet and from the winding current. Magnet mmf is a function of the flux-density level at which the magnet is operating. The B-H demagnetizing curves of the particular magnet material give the value of the magnetizing force H from which the magnet mmf can be calculated.

where lm is the radial thickness of the magnet (in meters if H is AT/m). Constant k depends on the units selected.

Winding ampere-turns are determined from NI where N is the number of turns carrying current I in amperes. However, care must be taken to include only the coils within the magnetic circuit being analyzed and not the complete winding. The magnetic circuit consists of a pair of magnets, a north and a south. In addition, with some of the recently developed winding configurations, some coils could be totally subtended within one magnet, e.g. when a coil is wound around one tooth or when the number of teeth is odd and there are fractional coils facing each magnet.

Reluctance calculation for steel requires the value of permeability, which in turn requires the B-H curve of the material.

Figure 3. B-H Curves for Typical Steel Materials.

Figure 3 shows B-H curves for typical steel materials. There are many mathematical ways of representing these nonlinear curves. The simplest and least accurate is the approximation of the curve to two straight lines joined at the knee of the typical curve. Much more accurate results can be obtained from such mathematical simulation (Ref. 4)

where , and are constants solved from the equations written by inspecting the B-H curve.

There is also the issue of orientation of steel and its effect on the calculation. In most magnetic circuit designs for DC motors, non-oriented silicon steel is used. Although flux simultaneously flows in many directions inside the steel, the same B-H characteristics can be used in all parts of the steel.

**Armature
Reaction**

High armature reaction and the potential for magnet demagnetization are major parameters in optimum DC motor designs. Because most applications require the motor to run in both directions, the approach of shifting the commutation axis cannot be employed.

Figure 4. Armature Reaction.

Armature reaction is illustrated in Figure 4. This is the interaction of the magnet field axis with the winding axis. It unbalances the magnet and winding flux patterns. It also demagnetizes the magnet in the region where the two fields oppose each other. If soft iron magnetic saturation levels were of no concern, the net impact of armature reaction would be minimal. However, that is certainly not the case in high-power density motors.

Figure 5. Second Quadrant B-H Curve for Nedodymium iron.

Figure 5 shows the second quadrant BH curve for Neodymium iron. The high intrinsic coercive force eliminates permanent demagnetization of the magnet. However, partial demagnetization of the magnet can still occur and can significantly lower the magnet flux. Figure 5 also shows the effect of partial permanent magnet demagnetization. It can be computed using Equation 6.

Z = Number of winding conductors

= Motor peak current (amps)

= Magnet radial thickness (cm)

p = Number of magnet poles

n = Number of winding parallel paths

This equation shows that a higher pole count and a thicker magnet reduces the demagnetization value. The term is the traditional NI or ampere-turn parameter that must be maximized for the high power performance needed in this motor.

**Motor
Performance**

After determining the magnetic circuit parameters, flux and flux-densities in all its parts, the next step is to calculate motor performance. Either the back EMF constant or the number of turns can be assumed or it can be calculated from Equation 7. The value of flux is the parameter that links the two together.

Here, p is the total number of poles, N the number of turns/coil, m the number of coils/phase and is the back EMF constant in V/krpm.

Next, , the torque constant, can be determined from using one of the two following equations:

The developed torque can be determined approximately from the following Equation 10.

where I is the winding current. This approach assumes that is a constant value not affected by current because it is based on the no-load value of . Naturally, this is not absolutely true since the load current has a continuously changing effect on and . Alternatively, the torque can be calculated from the stored energy in the magnetic circuit.

where is the magnetic energy and is the angle of rotation. The energy itself can be determined from the parameters already calculated.

The above two methods have different accuracy levels that may still not match the values determined by tests. In these tests, torque and current are measured at various values of load. The ratio of torque to current is then the at that particular amount of load.

**Loss
Calculation**

The most relevant losses are the electrical and magnetic losses. Electrical losses are simply the losses in the copper wire. Once the winding selection of turns and wire resistance can easily be determined. The calculation of the electrical losses is then fairly straightforward even with temperature effect on the resistance taken into account.

Magnetic losses are commonly called iron or core losses. The difficulty in the calculation is similar to that encountered previously in the calculation of reluctance. Both are highly dependent on the flux-density level that the material is working at. Since the flux-density levels are different at different parts of the magnetic circuit and also at different times during one revolution of the rotor, an approach similar to that in the reluctance calculation is taken here.

The magnetic circuit is divided into many sections and each section is small enough to allow the following assumption to hold true. The flux-density in any particular section is assumed to be constant throughout that section. The accuracy of the results, therefore, improves as the magnetic circuit is divided into smaller and smaller sections. However, too many small sections enlarge the program and lengthen the solution time, and this is not intended here.

Iron losses include hysteresis and eddy-current losses. The variables affecting these losses are the flux-density and the frequency of flux reversals. This can be seen from the following two equations.

where:

= Hysteresis loss (W/kg)

= Hysteresis constant

f = Frequency (Hz)

B = Flux-density (T)

= Eddy-current loss (W/kg)

= Eddy-current constant

t = Material thickness (mm)

Steel material manufacturers usually publish the above losses combined. The units are Watts per unit weight and the curves are plotted against the flux-density level of the material. The curves are sometimes made available at various values of frequency. This is the exception, however. In a similar fashion to the B-H curves discussed earlier, the loss curves can also be represented mathematically. Thus, the computer program determines the iron losses for the continuously changing flux-density. In the absence of the loss data at the required frequency, the effect of the higher frequency (since the curves are given at 60 Hz) on the loss is taken into consideration on the basis of the above formulae.

To determine the iron loss of each section of the magnetic circuit, the weight of each section must be determined. This is an exercise in geometric calculation of the cross sectional areas and lengths of these sections leading to the volume and weight by knowing the material density.

In summary, the calculation of the iron losses requires a large number of formulae because of the complicated geometry and large numbers of variables. The following list summarizes the steps taken in those calculations.

1. Divide the magnetic circuit into a reasonable number of segments.

2. Determine the weight of each segment.

3. Determine flux-density in each segment.

4. Determine the iron loss in each segment of the soft iron at the appropriate B level.

5. Determine the effect of the operating frequency on the iron loss (eddy-current and hysteresis).

6. Determine the total losses by multiplying items 2 and 5 for each segment and adding the results for all the segments.

A recent presentation (Ref. 5) was written around the results obtained from using this computer program. In that work the losses were calculated for different types of steel. The material selection was demonstrated to have a significant effect on the total losses, especially since the investigation was aimed at high speed motor performance. Figure 6 is a printout of a typical design performed with the computer program.

Figure 6

References:

1. C.K. Taft**, **H.D. Chai,
"A Parametric Model for Brushless DC Motors," 1987, IMCSD Proceedings.

2. A. Cassat, "Brushless DC Motors - Torque and Inductance Determination," 1987, IMCSD Proceedings.

3. D. Jones, "Motor Magnetics Design Changes in Higher Energy Rare Earth PM DC Motors," 1985, IMCSD Proceedings.

4. J.R. Brauer, "Simple Equations for the Magnetization and Reluctivity Curves of Steel," 1975, IEEE Magnetics.

5. D. Jones, G. Gogue, "The Impact of High Speed Iron Losses in BLDC Motor Operation at Speeds Above 10,000 RPM," 1988, IMCSD Proceedings.

6. J. Ireland, "Ceramic Permanent Magnet Motors," McGraw-Hill Book Company.