## Friday, July 12, 2019

### BLDC (PMSM) end turns and torque production

Jan 2020 update: I have removed the discussion of end turns in iron-core motors from this post as I have since learned that the arguments used were flawed. I will address this topic in my next post.

End turns are not considered when calculating the torque of a hobbyist grade BLDC (PMSM) electric motor. This post explores why.

## 1.End turns and BLDC motors

End turns are the region of a BLDC motor's windings which extend out from the slots at either end of a motor.
Current supplied to a motor passes from one slot to the next via the end turns. However, end turns are not typically considered when calculating the torque produced by a BLDC motor in 2D FEA software packages such as FEMM.

Although end turns can be ignored when it comes to calculating torque, they do contribute to the resistance of the windings. It is therefore ideal to have the end turns be as short as possible when designing a motor so as to improve motor efficiency by minimising ohmic losses in the windings.

Electric motors with concentrated windings, like those found in most hobbyist outrunners, have shorter end turns than those used in motors with distributed windings. Motors that use concentrated windings are therefore more efficient in some use cases, such as in high torque density, low-speed applications. However, as only a limited number of slot and pole combinations are available for concentrated windings they are not suitable for all use cases.

## 2. Torque generation and end turns in core-less motors

To understand why end turns are not considered when calculating the torque output of a BLDC motor let us consider a simplified scenario. In the diagram below a single turn copper wire is placed inside a homogenous magnetic field provided by permanent magnets placed above and below the wire. The magnets are arranged so that the applied field direction is different for either side of the wire.

This situation is an analogous to core-less electric motors which have no ferromagnetic material in the stator. If a current flows through the wire then a Lorentz force is produced on that wire as described by:

Force = Magnetic field \timesCurrent\times Conductor length

where the magnetic field (B) is supplied from the magnets, the current (I) is that flowing through the wire of a given length (L). Looking at this setup from the top down we see the following.

It is clear that a Lorentz force is only generated on the wire which is located between the magnets, which is analogous to the wire located within the slots of a BLDC motor. The parts of the wire which extend outside of the magnets are not exposed to a magnetic field and so the force generated is zero. This situation is the same for end turns in a BLDC motor, and so is ignored by FEA software packages.

However, for the sake of argument, let us also consider a situation where the magnets could be made long enough that they also covered the end turns. In this situation, a force is generated by the end turns. However, this force is tangential to that generated by the wires in the 'slots' and so would not contribute to the production of useful torque within a motor. This force would also be mostly cancelled out as can be seen from below.

Another arrangement I have seen used by DIY motor makers is to take round coils and magnets and arrange them to make an axial flux motor or generator.

While this approach will work, it is also inefficient since the force produced on the wires is always tangential to the flow of current and so a significant fraction of the force is not available to useful work.

Note that although all of the above examples are for core-less motors, the situation is unchanged for iron cored motors.

In my next post, the concept of flux linkage will be explored and used to understand why end turns still do not contribute to flux in conventional iron-core motors.

## 4. Conclusion

End turns in air-core motor designs are not considered in the calculation of torque as they are either outside of the region exposed to the permanent magnet's magnetic field or the force generated is in the wrong direction.

Equations were produced in this post with the help of arachnoid.com. If you have noticed any errors in the above article then please let me know.

## Saturday, May 11, 2019

### Understanding BLDC (PMSM) electric motors: Base speed, no load speed and torque vs speed

Search for the keywords "electric motor speed vs torque" and you will find many hundreds of images, with each looking different to the next.

 A selection of 'electric motor speed vs torque' charts
This has, understandably, lead to considerable confusion (myself included) when discussing 'hobbyist' grade BLDC (PMSM) motors and the topics of:
• base speed
• torque vs speed
where 'speed' is referring to how fast a motor is spinning (i.e. it's angular velocity). Therefore, this post will explore these three topics. For those wanting a quick summary see the brief overview below.

## 1.0 Brief overview

### 1.1 Voltage constraints

The amount of torque that the motor can produce continues to increase as you approach zero speed. At zero speed (i.e. shaft is stopped) the current flowing through the motor will be about equal to the supply voltage divided by the motor winding resistance. This is the 'stall current' of the motor and will typically be a very large value due to the low resistance of most brushless motors (i.e. < 0.1 Ohm) as is described below.

### 1.2 Thermal constraints

When in normal use, the torque output of a hobbyist brushless motor increases roughly in line with the current supplied to its windings. Therefore, if you double the current to the motor you also double the torque output. However, the resistive losses in the copper windings increase with the square of the current and so if you double the current and torque you are also producing four times as much waste heat.

Since most brushless motors have a very low winding resistance (i.e. 0.1 Ohm) it is possible to deliver extremely large currents to a motor even with a modest supply voltage. This is of course provided you have a big enough motor controller and power source to deliver the required current. However, the use of extremely large currents will quickly overheat and destroy a motor. Therefore, most hobbyist brushless motors are not limited by how much current you can deliver to the motor, but rather by how much current you can deliver before the motor overheats. That is to say, the torque output of most hobbyist brushless motors is thermally constrained.

Similarly, when a brushless motor is rotating it also produces waste heat due to core losses which are the losses generated in the iron and magnets. The faster you spin a motor, the larger these losses. Therefore your motor will overheat more quickly at higher speeds than at lower speeds for the same torque output.

This thermally constrained torque output of a motor is also often referred to as a motors rated torque. The maximum speed at which a motor can operate forever at its rated torque is known as a motors base speed.

It is possible to exceed a motors rated torque for short periods of time. The true maximum torque output of a motor is its peak torque. The peak torque output is usually defined as the maximum torque that can be sustained for some fixed period of time (i.e. 30, 60 or 180 seconds) without overheating (starting with a cold motor) or without damaging the magnets or mechanical components. The peak torque of a motor can be estimated from its peak current value. The peak current value is what is most commonly quoted by sellers of hobbyist brushless motors. Unfortunately, how the peak current is measured is usually not clear. As a rule of thumb, the rated current, and therefore torque output, of a hobbyist motor will be around one third that of the peak torque value, or around two thirds with aggressive forced air cooling.

In summary:
• The no-load speed of a motor is voltage constrained.
• Useful work can only be done below a motors no-load speed.
• The constant torque output of a motor is thermally constrained.
• The maximum speed which a constant torque can be produced is the base speed.
• A motor can operate above its constant torque value for short periods of time
For a more in-depth discussion read on.

## 2.0 In-depth analysis

As mentioned above, the no-load speed (\omega_{NL}) of a motor is the maximum rotational speed that a motor can achieve without any external load (i.e. no torque output) placed on it. The theoretical no-load speed is given by:

$$\omega_{NL} = K_{V} \times V_{supply}$$
where \omega_{NL} is the theoretical no-load speed in RPM, K_{V} is the motor's velocity constant in RPM/V and V_{\s\u\p\p\l\y} is the voltage supplied at the terminals of your power source such as a battery.

\omega_{NL} can be visualised as a single point on  a plot of speed vs torque

In practice, the theoretical no-load speed can never be reached. This is due to:
• Friction and windage: All motors experience frictional and windage losses. Frictional losses originate from the bearings while windage losses originate from the relative movement of the rotor to the stationary air surrounding it. These losses require torque to overcome and so the motor can never reach its true zero torque no-load speed. High rotational speed in-runner motors or out-runners with built-in cooling fans will experience larger frictional and windage losses.
• Core losses: The relative movement of the rotor magnets past the stator generates hysteresis and eddy current losses within the stator, even if there is no current being delivered to the motor windings. Similarly, differences in inductance seen by the rotor magnets as they move past a slotted stator also generated eddy currents within the magnets themselves. These losses must also be overcome by the generation of torque by the motor.
• Resistive voltage drop: As some small amount of torque is still required to overcome friction, windage and core losses at the no-load speed there will also be a corresponding current supplied to the motor. This current will result in a voltage drop across the motor and the motor controller, reducing the supply voltage available to overcome the motors back EMF.
• Motor controller dead time insertion: Most motor controllers have some level of dead time when switching their MOSFETs off and on so as to prevent a short between the high side and the low side of the DC bus. This dead time means that even if the motor controller is commanded to output a 100% PWM duty cycle the motor controller MOSFETs will still be off for some fraction of the time, preventing 100% utilisation of the supply voltage. It's also likely that there are other issues on the motor controller side at high electrical frequencies that may also reduce the top speed but which I don't understand at this time.
• Motor controller PWM frequency and control loop frequency. If you are using a motor controller that supports FOC (also called vector control) then it will be producing a sine wave via PWM. Theoretically, if the frequency of this sine wave (the motors electrical frequency) surpasses approximately 1/20th the PWM frequency (the exact value will depend on the capacitance and inductance of the system) then the quality of the sine wave will start to degrade. This will cause a fall-off in motor efficiency and if taken to the extreme, loss of motor control as commutation falls out of step with the rotor position. Also, if the electrical frequency of the motor approaches the control loop frequency then control will also be lost.
Note that those losses which are internal to the motor (frictional, windage and core losses) produce no net force outside of the motor and so are not something that can be measured directly with a torque sensor. This is because the stator is both applying the torque to turn the rotor and the braking torque (losses) to stop the rotor. In general, core losses are going to be far greater than those losses due to windage and bearing losses.

The no-load speed which is realistically achievable (\omega'_{NL}) can, therefore, be approximated by:

$$\omega'_{NL} = K_{V} \times (V_{supply} - V_{drop}) \times Modulation_{max} - \frac{P_{loss}}{\tau _{loss}}$$

where V_{drop} is the combined resistive voltage drop in the system, \M\o\d\u\l\a\t\i\o\n_{max} is the maximum PWM modulation achievable by the motor controller after accounting for dead time insertion and \frac{P_{loss}}{\tau _{loss}} is the speed reduction brought about by the power and subsequent torque required to be produced by a motor to overcome friction, windage and core losses.

The size of the difference between \omega_{NL} and \omega'_{NL} will depend greatly on the type of motor used, its associated losses and the supply voltage. In general, it would be reasonable to expect that a typical hobbyist motor and motor controller (ESC) can reach 85% of its theoretical no-load speed. This is provided that it is operated within its 'rated' voltage.

### 2.3 Thermal constraints and no-load speed

Beyond voltage constraints, there are also thermal and mechanical constraints. If you try to run a '6s' (~24V motor) off a 240V motor controller then don't expect to come anywhere near its theoretical base speed before you run into thermal or mechanical (i.e. exploding rotor) limitations.

Core losses generated within a motor at its base speed can be significant. For example, an unloaded 12n14p HobbyKing SK3 280 Kv motor produced over 100 W of losses at 10,000 RPM.
 Power draw as measured from the power supply for an SK3 280 Kv motor
The peak speed of 10,000 RPM only represents 75% of the maximum speed available with the 48V power source used. Despite this, the 100 W of losses within the motor at this speed is more than enough to heat the motor to the point of being damaged if operated for an extended period. Also note that the generation of core losses has nothing to do with the motor Kv and is instead determined by the pole and slot combination, with lower pole numbers, in general, producing lower losses. This is one of the reasons why a 4 pole in-runner can operate at much higher speeds than a 14 pole outrunner.

Due to the no-load losses within a brushless motor, the maximum continuous no-load speed \omega''_{NL} may be far lower than the maximum achievable no-load speed.

This could be an important consideration in some applications which require high speed but very little torque such as an engraving tool or mirror galvo optical assembly.

### 2.4 Base speed

Also mentioned above, the base speed (also called the rated or nominal speed) is the top speed at which a brushless motor can achieve at it's 'rated torque'. The 'rated torque' is the torque that a motor can output indefinitely without overheating. This speed is different to that of the achievable no-load speed because there will be a larger voltage drop across the motor and motor controller when the motor is producing its rated torque due to the requirement that a larger current be supplied to the windings.

In all other ways the base speed of the motor (\omega_{BS}) can be defined in the same way as the maximum achievable no-load speed:

$$\omega_{BS} = K_{V} \times (V_{supply} - V_{drop}) \times Modulation_{max} - \frac{P_{loss}}{\tau _{loss}}$$
with the exception that V_{drop} is now dependent upon the current required to produce the motors rated torque. Therefore, the base speed of a motor will fall off as the 'rated torque' is increased due to the additional voltage drop incurred at higher currents.

In the figure below \omega_{BS} represents the maximum speed at which a motor could produce a set level of torque. However, if less torque was required from the motor then a higher speed (following the light blue line) could be achieved.

The gradient of the light blue line will depend on the resistance of your motor and motor controller. Higher resistance will produce a shallower gradient.

### 2.5 Base speed and temperature rise

So far we have assumed a fixed motor and motor controller temperature, and therefore a fixed voltage drop per amp. However, as the temperature of a motor is increased so too does its winding resistance. For example, a 50C temperature increases the winding resistance by 20% and a 135C rise by 53% [1]. Therefore, a hot motor will have a larger voltage drop across the motor and a lower base and no-load speed.

Note that the above diagram is not to scale and the actual reduction in base speed for a 'hot' motor may be much smaller in reality.

While it's important to keep this point in mind, for simplicity the rest of this post will assume a constant motor and motor controller temperature.

### 2.6 Base speed and additional cooling

Adding additional cooling to a motor will allow you to increase the 'rated torque' that a motor can produce without overheating. Counterintuitively, adding additional cooling to a motor and raising its rated torque, and therefore current, will actually decrease its base speed since a larger voltage drop will produce over the motor and motor controller.

### 2.7 Momentary peak torque

Momentary peak torque is important for many applications. For example, rapidly reversing the direction of an axis on a CNC machine or jumping a quadruped robot into the air. The ability to produce a large momentary peak torque is one of the largest selling points for a brushless motor over stepper motors.

The peak torque output of a motor is generally limited by four factors:

1. Voltage constraints: As mentioned above, the maximum torque produced by a motor is reached at zero RPM (stall torque), where the current is equal to the supply voltage divided by the total resistance in the system. Therefore, a larger supply voltage will allow for a larger current to flow at zero RPM and greater peak torque.
2. Current constraints: If using a small power source then it may not be possible to deliver the stall current demanded by the motor without the supply voltage 'sagging' under the load.
3. Thermal constraints: Supplying a motor with a momentary large current will result in a small temperature rise. Supplying a 'momentary' large current repeatedly, or if the motor is already hot, and this temperature rise will accumulate to the point of damaging the motor. If you demand 100 A from a 10 A motor controller, then even a 'momentary' peak current and torque will likely cause it to overheat.
4. Material constraints: Ultimately, even if all other constraints are overcome, the peak torque output of a motor is limited by the materials used in its construction as described in more detail below.

• Material constraints: Stator saturation

The magnetically soft materials used to construct the stator, typically laminated Fe-Si steel, can only be magnetised (polarised) so far before it reaches magnetic saturation. As magnetic saturation is approached there is a loss of linearity between the current supplied to the motor and the torque it produces. i.e. a 1 A increase in current may produce a 0.1 N.m increase in torque when the motor is far from saturation, but only a 0.05 N.m increase in torque close to saturation.

The image below is taken from Ben Katz thesis and nicely shows a subtle ~12% fall off in the torque constant of a motor as its current is increased due to the stator approaching saturation.

 All credit to Ben Katz [2]. Click to enlarge

• Material constraints: Magnet demagnetisation
Supplying a motor with a very large current produces an equally large magnetic field. If the size of this magnetic field exposed to the rotor magnets starts to approach the coercivity of the magnets then there is a risk of demagnetisation. If the rotor magnets are significantly demagnetised then the motor is permanently ruined unless the magnets can be replaced. Thankfully, modern rare earth magnets have a very large coercivity. My own simulations suggest that supplying a current large enough to saturate a motor still does not produce a demagnetise risk for a typical grade N35 Nd-Fe-B magnets at room temperature.

However, the situation changes if the magnets become hot. Increase their temperature to 80C and their coercivity is about one half that at 20C. At this point, it becomes possible to demagnetise the magnets with a current approximately three times a typical brushless motor rated peak current.

 The magnet coercivity (x-axis crossing point) falls significantly with temperature.

• Material constraints: Mechanical limitations
Lastly, the shaft and fixings of a motor are all designed with its rated torque in mind. Exceed this value by too large a margin and you risk damaging the motor. i.e. slipping shaft where it connects to the rotor, twisting the frame etc.

### 2.8 Torque vs speed

From the discussion above it should now be clear why searching for a 'torque vs speed' plot yields so many results. The output torque of a brushless motor at a given speed depends on many different factors. However, these factors can be simplified by only considering a single situation which should apply to most hobbyist brushless motors and motor controllers. The assumptions are as follows:

• A powerful power source, such as a Lipo battery, is used which can supply current to the motor when peak torque is demanded without much in the way of voltage sag.
• Peak torque is infrequently and only momentary required from the motor so that the motor and motor controller is maintained at a reasonable temperature of around 60C at all times.
• The supply voltage is matched to a motor so that the motor will not overheat if run at its no-load speed indefinitely.
• The peak torque, and therefore peak current, supplied to the motor is still well below magnetic saturation of the stator and within the safe mechanical limits of the motor.
 A typical setup
In this scenario, the torque output of a motor below its base speed is thermally constrained while the maximum speed is voltage constrained. The peak torque that can be supplied to the motor, therefore, depends on the temperature of the motor and the length of time that torque is requied. If the motor starts at room temperature then a very large torque can be produced for a short period. Alternatively, more moderate peak torque can be supplied for a longer period.
However, if the motor was already at its safe upper limit of 60C then only the motors rated torque could be supplied without heating it to the point of damage. Therefore, the available speed vs torque curve of a brushless motor also depends on its thermal history. For this reason, industrial brushless motors will typically include a temperature sensor embedded in the windings so that the motor controller can actively limit the current if the motor becomes too hot. This is also why the brushless motors supplied by Odrive robotics come with an embedded thermistor.

Note that the power consumption of a motor producing its 'rated torque' at its base speed will be significantly higher than the same motor producing its 'rate torque' at the stall speed (zero RPM). Therefore, unless you have a powerful power supply or battery you may be unable to reach the base speed at rated torque of a motor. See this post for more details.

## 3.0 Conclusion

In summary, the torque that a motor can produce at a given speed depends on the motor used, how the motor is cooled, the thermal history of the motor and the power source used. For this reason, there is no single universal 'speed vs torque' plot for any given motor type, or even for two identical motors used under different conditions.

However, it is a safe assumption that any given brushless motor will have a flat constant torque output vs speed due to thermal constraints. Similarly, the peak torque output will, in general, be around 2 to 3 times the constant torque output for a short period (i.e. 30 s) when starting cold. Lastly, the base speed of a motor can be approximated as roughly 0.85 times the supply voltage multiplied by the motors Kv.

In the next few posts, I will be using a low-cost DIY motor dynamometer to characterise these torque vs speed curves, along with motor efficiency, for serval of hobbyist grade brushless motors.

Equations were produced in this post with the help of arachnoid.com. If you have noticed any errors in the above article then please let me know.

• [1] From  Design of Brushless Permanent Magnet Motors by Hendershot and J. R., Miller, T. J. E.
• [2] Katz, Benjamin G. A low-cost modular actuator for dynamic robots. Diss. Massachusetts Institute of Technology, 2018.

## Saturday, February 9, 2019

### How to model a BLDC (PMSM) motors Kv (velocity constant) and Kt (torque constant) in FEMM

In a previous post, the velocity constant (Kv) and torque constant (Kt) was estimated for a number of different 'hobbyist level' BLDC (PMSM) out-runner motors. In this post, I will be using FEMM to model one of these motors as accurately as possible. In the process, you will see how to use FEMM to model such a motor and how accurately its results compare to the real thing.

## Measuring key parameters prior to modelling

Before the motor can be simulated in FEMM we first need to determine its key parameters. The motor to be simulated is an Odrive Robotics N5065 motor which has previously been characterised in this post. Note that this motor is a pre-production unit that was sent to me for testing by Oskar of Odrive Robotics and differs slightly (5mm shorter stator, different winding number) from that offered in his shop today.

The key motor parameters are:
• The number of slots in the stator and the number of poles in the rotor
In order to make this, and the following steps, easier the motor was first disassembled. Disassembly was achieved by the removal of a single circlip from the shaft and pulling the rotor bell away from the stator. Some force is required due to the attraction of the magnets to the stator iron.

With the rotor bell removed it was simply a matter of counting the number of slots/teeth on the stator and the number of magnets on the rotor. This motor has 12 slots and 14 poles (12n14p). The likely reason that the manufacturer chose this number of slots and poles is explored in detail here.
• The stator lamination dimensions and tooth profile
The stator lamination dimensions and tooth profile can be difficult to measure directly. If you don't mind destroying your motor then you can remove all the windings and separate a single stator lamination (or just section the entire motor with a slow speed saw) and use a flatbed scanner to capture the dimensions of the lamination without any perspective errors. Including a ruler when scanning will allow you to set a correct scale. Then you could use a software package with edge detection (such as Inkscape, which is free) to produce a DXF file without the need for any other CAD software.

In my case, I didn't want to destroy or unwind the motor and so I instead opted to model it as best I could in F360 (free for hobby/student use) and then export the sketch as a DXF file for FEMM. By doing this in with parametric design software like F360 it also allows me to play around with motor parameters in future without the need to start from scratch.

The sketch was started with the inner and outer diameters and then a single tooth and its windings were added. This tooth was then copied using the 'circular pattern tool' to make a full stator.

 F360 sketch of the stator
While this method is far from ideal it does give a good approximation as we shall see.
• The rotor magnet dimensions and grade
The magnets are glued to the rotor bell with epoxy and are also held at the correct separation by segments machined into the end cap and also into an aluminium retaining ring at the open end of the rotor bell.

 The silver magnets are retained in the black rotor bell with the grey epoxy and the silver coloured aluminium ring in the open end
This approach gives good retention of the magnets but does make it difficult to estimate the actual thickness of these magnets. I made my best approximation by taking the outer bell diameter and the inner flat-to-flat measurement of the magnets and assuming the rotor bell thickness (the black walled region in the image above) is constant. Using this data, and the assumption that the magnets are flat and not curved the rotor was also modelled in F360.

 F360 sketch of the rotor
Since the magnets are assumed to be perfectly flat and are glued to the inside of a cylinder there is a small 'air space' (actually filled with epoxy) below the magnets. When I imported this sketch as a DXF into FEMM later these small gaps were removed by FEMM. If you wish to maintain them for accuracy then you can force FEMM to recognise them by including a small line segment between the rotor bell and the magnet that crosses the thin air gap region.

The magnet grade specified by the motor manufactures is N42. If you don't know the magnet grade of your motor then you can make a fairly safe guess by assuming N35 if its a cheaper model or N42 if its a higher end model. As there was in no material option in FEMM for N42, I went with N40 instead. In my case, FEMM suggests that using over N40 over N35 grade magnets leads to a ~7% increase in torque.

The rotor and stator sketches were exported as DXF files ready for importing into FEMM. You can find a copy of DXF files here. Note that FEMM doesn't know the difference between 'construction' geometry generated by F360 and normal sketch lines and so you may want to remove them before exporting as a DXF file.
• The diameter of the wire used, the turn number per tooth and the number of parallel strands
Since I didn't want to unwind the motor it was difficult to determine the number of strands and turns per tooth used. Isolating a single strand and measuring it with a calliper suggested the diameter to be around 0.35 mm and so the diameter of the wire without copper insulation was assumed to be 0.34 mm.

The number of parallel strands was estimated by isolating a segment of wire running from tooth to tooth and poking it with a pen until most of the individual strands were visible.

I count a total of 8 strands. Luckily for me, the number of turns per tooth (ten) was written right on the side of the stator stack.

However, it should be possible to estimate the number of turns with a reasonable margin of error from the winding resistance of the motor and the wire diameter if needed.
• The method of winding (i.e. star or delta)
In my previous post, I had assumed that this motor is star wound. This turns out not to be the case. As can be seen in the image above one of the leads coming from the motor has 16 parallel strands. This suggests that the motor is delta wound, with the ends of each phase winding (8 parallel strands) grouped together with the next phase for termination to give a total of 16 strands. This was further supported by the fact that that I can't see any region in the motor where a required join would be made for the natural point of a start wound motor. It is likely that all of the motors I previously look at are delta and not star wound.

Other motor parameters include the lamination thickness, the lamination hysteresis loss, skin and proximity effect of the windings, friction of the bearings etc. can be ignored here as we are only considering the motor under static conditions.

## Modelling a motor in FEMM

Getting started with FEMM is fairly straight forward if you can look past its outdated interface. Before starting the motor simulation below I would highly recommend following this quick tutorial so that you are familiar with the interface first.

When ready to model your motor, start a new magnetics problem. From the menu bar select 'problem' and check that you are using the correct units with respect to your DXF file and that have set up a planar simulation. You will also want to check that the 'depth' of your simulation matches the length of your rotor. In this case, my rotor was 30 mm long.

The next step is to import your DXF file(s) using the file -> import DXF option. This can take some time depending on the amount of geometry. If you have any 'construction' lines left over from your CAD model now is also time to remove those. Note that it is also important that the centre of your model be located at point (0,0) for torque estimation to be conducted later.

With your geometry imported its now time to assign material properties. In my case, the rotor bell was assumed to be M-22 Fe-Si steel which is inaccurate. I should have used hot rolled low carbon steel or similar. However, this is unlikely to have had a large impact on the results. The shaft was selected as hot rolled low-carbon steel. The laminations were chosen as M-22 Fe-Si steel which is a common grade used for electric motors.

Magnet placement was done using the Lua scripting language, which is natively supported by FEMM. You can find a copy of the script I used here. Simply enter the number of poles, the radius at which the magnets can be found and the angle (0 degrees being horizontal) that the first magnet can be found. The script will then place each material assignment at the correct orientation. If you wish to change the magnet grade just change the name of the block property in the script (i.e. "NdFeB 40 MGOe")  to that of the magnet you wish to use from the list in FEMM.

The windings were then assigned based on the number of turns and the number of strands. In my case, I created a new 'block property' using the property definition dialogue box. In this case, I selected plain stranded wire with a 0.34 mm diameter and 8 strands and then assigned the same properties used by the other copper conductors.

The actual order which the windings were placed was determined using this handy winding layout tool. By selecting 12 slots, 14 poles, two-layer winding (single winding per tooth as opposed to a single winding for every second tooth) the following is seen.

Using this as a template the winding phase and direction (either into or out of the page) was assigned. This was the most tedious process of creating the motor in FEMM. With all the materials properties assigned we see the following.
The regions filled by air were assigned (one for the air gap region inside the rotor bell and one for outside the motor) and the open boundary builder (circle within a circle button) was used to set the boundary conditions and default settings are fine here.

Finally, we need to set how much current is supplied to each winding. Since this motor is driven using field oriented control (FOC) we need only supply a current to one phase and a current opposite in sign, and twice as small, for the other two phases as per the figure below.

In addition, this motor is delta wound and so the current supplied to the phase will be equal to square root(3) of the line current as described here. In my previous post a phase current of 40 A produced a torque of 1.2 N.m. However, this value of 40 A is the line current (or phase current for a star wound motor) and so we need to divide that value by sqrt(3) to get the phase current that our delta wound motor saw. This gives a value of + 23.1 A for phase A, and - 11.55 A for phase B and C. These values are set in the circuit dialogue box with a series connection.

With this, the motor is now ready for simulating. The complete FEMM file used can also be found here. Note that in order to run this script the minimum angle setting (Problem -> Min. Angle.) needs to be set to 20.

## Torque constant estimation

The torque constant (Kt) of an electric motor can be estimated using FEMM with only a single simulation. However, you will want to make sure that your rotor is rotated so that it will produce its maximum torque per amp. In this case, peak motor torque is given when the rotors d-axis is 90 electrical degrees (360 electrical degrees is equal to 360 mechanical degrees divided by the number of pole pairs, which in this case is 7) ahead of the q-axis. A simpler 6n8p motor schematic below for illustration purposes.

Simulating the 12n14p N5065 motor with the rotor 90 elec. degrees ahead of the d-axis looks like the following.

To determine the torque produced the rotor needs to be selected using the area selection tool (green box inside four points) by clicking on each of the magnets and the back iron as shown below.

With rotor selected we can find the torque around point 0.0 using the integrate option and 'Torque via Weighted Stress Tensor' selection from the drop-down menu as recommended in the FEMM manual.

The result is a torque of 1.37 N.m.

For reference, the real motor measured 1.2 N.m. Therefore the FEMM value is an overestimation of about 15%. Considering the number of assumptions made with this model I was actually surprised that the FEMM estimation is as close as it is.

## Kv, back EMF estimation

The back EMF constant (Kv) can be estimated using only the torque constant and the motor current as described in this post. In this case, our line current of 40 A and FEMM modelled torque of 1.37 N.m gives

$$K_{V} \approx \frac{8.3 \times 40 }{\tau } \approx 242$$

The listed Kv of this motor was 270 Kv and so once again we are only off around 15 %.

In order to estimate the shape of the back EMF we need to do a series of simulations at different rotor angles and plot the torque output. This torque output can then be converted to Kv in the same way as shown in the equation above. In order to automate the rotor incrementing process, I used a Lua script which can be found here.  In short, it rotates group one by x degrees around point 0,0 and solves. At the completion of solving it prints the torque around point 0,0 to the console and takes a screenshot of the resulting output. Thefeore the use of this script means that you need to assign each line segment and materials property of the rotor to group one as is shown in this video. Don't forget to copy the contents of the console output into another file and the completion of your simulation or otherwise modify the script to output data to a file. To use this script, place it in the same folder as your motor's .femm script so that your file is listed in line 8 "open(mydir .. "Your file name here.FEM")". You can also set how many steps you wish to make and how many degrees of rotation you desire.

Note that solving for each position may take up to 30 seconds and so running hundreds of steps will take many hours to complete. Doing this for the D5065 motor model we see the sinusoidal torque output with rotor angle below.

This sinusoidal torque with rotor angle, and therefore back EMF as they are the same thing, agree well with the back EMF seen for this motor.

Using the excellent, and free, OBS studio in 'slide show' mode was also able to stitch all the separate solutions .bmp files that my script outputted and create this clip of the flux distribution for different rotor angles.

## Cogging torque estimation

The cogging torque (torque due to the attraction of the rotor magnets to the salient teeth of the stator) was estimated using the same method as the Kv estimation. The key difference being that the armature current was set to zero and a much finer step was used. The figure below shows the raw output data in black which is quite noisy and an FFT smoothed fit of the data in red.

As discussed in the last post, the number of cogging steps per rotation is given by the least common multiple (LCM) of the number of stator slots and rotor poles. The LCM of this 12n14p motor is therefore 84. A quick count of the peaks in the smoothed data above and extrapolated out to 1 rotation (360 mechanical degrees) and a cogging frequency of 84 is confirmed. As will be discussed in more detail in a future post, the cogging torque amplitude depends strongly on magnet placement, magnet shape, stator tooth shape etc. and so the magnitude and shape of your cogging torque estimated be FEMM will only be accurate if you have modelled your motor exactly.

## Modelling motors with symmetry

While it is not shown here, it is important to point out that it is often not necessary to model an entire motor to determine its properties. For motors with symmetry you may only need to model a small segment of the motor and from that, you can extrapolate the results of a full-size motor. There are also a few different tools available in FEMM that may be worth exploring.

## Alternatives to FEMM

As can be seen from the results above FEMM is capable of approximating some properties of an electric motor. However, there are many alternatives that are also worth considering. If you would like to use FEMM on platforms other than windows then take a look at the xFEMM project. While xFEMM does away with the user interface it is also considerably faster (up to 80x in some instances) and is still a free open source.

If you require simulations of 3D bodies or more complicated AC loss analysis etc. then you will need to move on to an expensive commercial FEA package such as the Maxwell from Ansys or Comsol Multiphysics

## Conclusion

Modelling an electric motor in FEMM is no easy undertaking. However, by carefully replicating all of the required components it is possible to get a reasonably accurate estimation of a motors Kt, Kv, back EMF shape and cogging torque parameters. Kt and Kv was shown to be within about 15% of that measured for a real motor and so FEMM can act as a nice starting point for verifying ideas when designing an electric motor.

## Sunday, January 27, 2019

### Selecting the best pole and slot combination for a BLDC (PMSM) motor with concentrated windings

The full spreadsheet used in this post is available for download. To modify the sheet just choose the 'Make a copy option' from the drop-down file menu.

## Introduction

Concentrated windings have a single coil per tooth and are commonly used on 'hobbyist' style BLDC (PMSM) out-runner motors.

A key advantage of concentrated windings is that they can be quickly and cheaply wound by machine as seen in the video below.

Even motors with very thick conductors can be wound by machine.

Other advantages include short end turns, which do not contribute torque to the motor and only increase the winding resistance, and space for effective air cooling. A major drawback of concentrated windings is that without careful consideration of the stator slot number and rotor pole number the performance of the motor will be poor.

This post will examine the advantages and disadvantages of different slot and pole combinations. The information is based off the paper 'Distribution, coil-span and winding factors for PM machines with concentrated windings' by S.E Skaar et al. and the book 'Design of Brushless Permanent-magnet Machines' by J. R. Hendershot and T. Miller. See the 'recommended reading' list above for more information on this book. The emotor.com glossary page is also useful for reference.

## The winding factor

An electric motors 'winding factor' (not to be confused with its copper fill factor) is a number between 0 and 1 which represents the fraction of the armature current which is used to produce torque.

From emotor.com's glossary page the winding factor is defined as the following:
The winding factor for a specific winding expresses the ratio of flux linked by that winding compared to flux that would have been linked by a single-layer full-pitch non-skewed integer-slot winding with the same number of turns and one single slot per pole per phase. The torque of an electric motor is proportional to the fundamental winding factor.
Despite being so fundamental, the calculation of the winding factor is not often discussed in textbooks. I found that the paper mentioned above has the easiest to follow description, although I had to rely on the emotor.com description of an unbalanced winding as the solution provided in the paper appears to not cover all scenarios.

I will not be going into detail about the calculation of the winding factor. However, you can find the spreadsheet used to calculate the values shown in this post here. This spreadsheet includes up to 108 slots and 70 poles and is easily extended further where needed.

The winding factor for a 3 phase machine with a non-skewed rotor can be seen listed in the table below where the top row is the number of poles and the left-hand column is the number slots. Since this is a three-phase machine the number of slots increases by three and since a magnet has two poles (a pole pair) the pole number increases by two.

It can be seen that some slot and pole combination have a winding factor of one while others are approaching zero, or even negative. However, not all of these slot and pole combinations can be used as is described below.

## 1. Exclude windings where q < 0.25

The slot to pole ratio with consideration for the number of phases is designated by the variable q. If q is less than 0.25 then the arc covered by a rotor pole is now less than half a stator tooth. This results in multiple north and south magnet poles interacting with each stator tooth and so the torque generated by the motor is reduced. Therefore, q values of less than 0.25 are generally not considered feasible and can be eliminated. The images below were created with this winding layout tool. Note that in reality, the gap between each stator tooth would be much smaller.

 Example winding layouts that give a q value that is not considered feasible.

 The calculated q-values for different slot and plot combinations.

 The winding factor of each slot and pole combination with those combinations that give q<0.25 removed.

## 2. Exclude windings where q > 0.5

Alternatively, if q is greater than 0.5 then it no longer makes sense to use a concentrated winding as a single rotor pole will span over multiple teeth. Instead, a distributed windings would be used. Therefore, these combinations can also be excluded for our purposes.

## 3. Exclude windings where Ns = Nm

If the number of slots (Ns) is equal to the number of poles (Nm) then the motor will produce large cogging torque and will no longer be self stating. This combination can therefore also be eliminated.

## 4. Removal of unbalanced windings and motors without symmetry

A motor with balanced windings will have the same number of coils for each phase per repeating segment of the motor. A motor must have balanced windings for operation. More information here

Also, it is ideal to have a motor that has a symmetry of at least 2 (i.e. the motor has two repeating sections) as this helps avoid unbalanced radial forces and noisy operation. In other words, a motor without any symmetry will produce its torque on only one side of the rotor. If a motor has a symmetry of two then its torque will be produced on opposite sides of the rotor, balancing the forces.

The slot and pole combinations that produced unbalanced windings and no symmetry can, therefore, be eliminated.

Based on the table above it can be seen that 12s10p and 12s14p are attractive combinations. This explains why these slot and pole combinations are so popular for 'hobbyist' out-runner electric motors.

## 5. Consideration of cogging torque

In addition to the winding factor, cogging torque is also an important consideration. Cogging torque creates vibration and noise during operation and acts to disturb the motor away from its desired position when used in servo applications. The cogging torque frequency is also closely correlated with the production of rotor losses both with and without current supplied to the armature.

The cogging torque frequency is given by the least common multiple (LCM) of the stator slots and the rotor poles. The figure below displays the possible slot and pole combinations with the winding factor replaced with the cogging torque frequency.

Selecting as high a cogging frequency as possible is desirable as it reduces the amplitude of the cogging torque. With this in mind, 24s22p would be an attractive option with its 264 cogging steps per rotation and 0.949 winding factor as opposed to only 84 cogging steps per rotation for 12s14p.

## 6. Higher order harmonics

So far we have only been considering the fundamental winding factor. However, different slot and pole combinations also affect the winding space harmonics. This topic will be covered in more detail in a future post once I have a better understanding of the topic.

From the last two tables above its easy to conclude that a 24s22p motor is 'better' than a 12s14p motor since it has both a higher winding factor and a higher cogging frequency. However, there are other factors that also need to be considered.

## 7. 1 Maximum electrical frequency

The electrical frequency (current sine-wave supplied to each phase) scales linearly with the number of pole pairs (Nm/2) of a motor. Therefore, doubling the pole count of a motor will increase its core losses by a factor of four since it scales with the square of the electrical frequency. Doubling the pole count will also double the back EMF produced and so twice the required voltage will be needed to drive the motor at the same RPM. However, rewinding of the motor with a lower number of conductors per tooth can be used to offset this increase in the back EMF and will not impact the efficiency of the motor. In addition, reducing the lamination thickness will help minimise eddy current losses, which scale with the square of the lamination thickness. See this spreadsheet for more details. At the extreme, the switching frequency of the motor controller used to drive the motor may also need to be increased in order to maintain a good approximation of a sine wave. This will incur additional dead time losses in the motor controller.

The electrical frequency required to operate the motor at 10,000 RPM is shown below.

## 7. 2 Mechanical winding considerations

Trying to squeeze more slots into a small motor is not always possible. Increasing the number of slots while keeping the motor diameter the same typically means that multiple smaller conductors must be used in order to bend around the tighter radii and make use of the available space. Smaller conductors have a larger fraction of their total cross section taken up by insulation and so the copper fill factor of the motor can be reduced even if the winding factor is increased. A lower copper fill factor means a higher current density in the winding and so higher I^2R losses.

## 7. 3 Cooling considerations

Effective cooling is critical for a high power density electric motor. Having a large number of slots and small conductors may congest the air flow path and reduce air cooling effectiveness. This, in combination with higher core losses due to an increase in the electrical frequency, can reduce the power density of a motor.

## 7. 4 Tooth-tip leakage

Squeezing more slots into a motor may also require that the gaps between the teeth be reduced. Reducing the size of this gap can allow some flux to jump (or zig-zag) from one tooth to the next, reducing the effective torque produced by the motor.

## 7. 5 Rotor skewing

If you plan on skewing your motor in an effort to reduce its cogging torque then it is important that the required skew does not lower the winding factor too much. The reduction in the winding factor is given by the skew factor discussed here. This topic will be covered in more detail in a future post.

## 7. 6 Reduced rotor and stator back-iron (yoke) mass

Increasing the pole and slot count of an electric motor has the advantage of reducing the flux density in the back-iron. This is because the higher number of stator slots and rotor poles means that flux does not need to travel as far around the stator or rotor to make a closed magnetic path. For every doubling of the pole or slot number, the flux density is halved and so the thickness, and therefore mass, of the back-iron can also be halved. This acts to increase the power density of the motor.

For example, Siemens appears to have increased the gravimetric torque density (torque per unit mass) by 50% for their 260 kW light aircraft electric motor by increasing its pole and slot count from an already high 36s30p to 72s60p. Note that this assumption is based solely on the images from their promotional material.

## Conclusion

The tables above provide guidance when selecting the number of slots and poles for a BLDC (PMSM) motor with concentrated windings. If the diameter of your motor is less than approximately 60 mm and if your maximum RPM is less than 10,000 then 10s14p or 12s14p may be a good choice. For larger diameter and/or motors that rotate more slowly a higher slot and pole count can be used, such as 24s26p with a reduced number of turns per tooth and thinner Fe-Si laminations. Higher pole count motors will have a higher electrical frequency, increasing losses, but will also have a lower mass since a reduced back-iron thickness can be used while also having a smaller cogging torque amplitude.