Magnet Eddy Current Loss and Rotor Magnet Segmentation
Solve 2D magnetic diffusion for slot and PWM harmonics. Compare circumferential and axial magnet segmentation, find the fewest pieces meeting a loss target, and estimate magnet temperature.
This magnet eddy current loss calculator helps engineers and students choose circumferential and axial segmentation for rectangular rotor magnets. Enter the magnet dimensions, material properties and alternating field spectrum in the rotor frame. The tool solves two dimensional magnetic diffusion inside each insulated segment, then searches the integer combinations that satisfy your loss and temperature targets.
Start with the dimensions of one parent magnet before cutting it. Circumferential segmentation divides its width; axial segmentation divides its length. The parent magnet count scales the answer to all magnets in the rotor. The selected counts let you investigate a proposed layout independently of the search limits. Set the minimum segment dimension to exclude pieces that cannot be manufactured or assembled.
For each active harmonic, enter its frequency as seen by the rotating magnet and its peak radial flux density at the segment perimeter. Slot passing frequency is the slot count multiplied by mechanical revolutions per second. PWM components must also be expressed in the rotor frame. Obtain field amplitudes from a field calculation or measurement; neither motor current nor DC magnet remanence can be entered as an alternating harmonic amplitude. Combine contributions at the same frequency as phasors before entering them. Up to six distinct frequencies can be summed.
The calculation retains the reaction field of the eddy currents. At low frequency, resistance limits the circulating current and smaller insulated pieces usually reduce PM rotor loss. When skin depth becomes small, the interior is screened from the alternating field. A few cuts can expose more material to that field and increase loss. The integer sweep bars and frequency response show this segmentation anomaly instead of assuming an inverse-square reduction.
The candidate is the arrangement with the fewest total pieces that meets both targets within the search range. Equal piece counts are resolved by lower loss. If none qualifies, the candidate is the lowest-loss arrangement in the range and the verdict says so. The loss map contains discrete design options; its axes count segments, not cuts. Increase a search limit if the useful region lies beyond the plotted range. The separate lowest-loss result is bounded by those limits and is not an economic optimum.
Temperature rise is total rotor magnet loss multiplied by the thermal resistance of all magnets together to the stated heat sink. Use resistivity at the expected operating temperature and update it if the predicted temperature changes substantially.
This tool does not calculate the machine's harmonic field, full three dimensional end fields, sleeve or rotor steel losses, electrical contact across cuts, or demagnetization margin. It holds the prescribed surface spectrum fixed as segmentation changes. Confirm a final high-speed machine design with its coupled electromagnetic and thermal model, especially when voltage-driven PWM, external magnetic reluctance or local hot spots affect the result.
How it is calculated
For each distinct rotor-frame harmonic, the radial phasor field satisfies
Here is the circumferential–axial rectangle of a single insulated segment. A constant field along its boundary implies zero normal eddy current there. The field and in-plane return currents are solved together; the radial extrusion has constant properties and field.
The implementation derives a separated series for this boundary-value problem. Put , , , , , with odd :
The volume-average normalized field is . Integration by parts gives
This equals the volume integral of . The factor 1/2 is required for peak phasors. The infinite-slab portion is summed analytically; only the end correction is truncated. This exact rearrangement improves convergence for strong skin effect. The complex arithmetic is evaluated through equivalent real and imaginary components; decaying exponentials avoid hyperbolic overflow. A small-argument expansion avoids subtracting nearly equal numbers.
Each pair of segment counts is solved for all active harmonics. Their time-average losses add. The target is . The search minimizes piece count among feasible designs, breaking ties by loss; it also finds the minimum loss in the bounded search. Mean rise is .
Ruoho et al. (2009) motivate finite current-return length; Huang et al. (2010) and Königs & Löhlein (2023), §§4–7, discuss reaction fields and nonmonotonic segmentation. This is a direct diffusion solution derived here, not their fitted circuit model. The axial segmentation paper by Ede et al. is the 2007 journal article (DOI below), rather than the 2002 control paper.
Assumptions
- Rectangular, identical parent magnets; equal insulated segments with negligible kerf and unchanged total magnet volume.
- Linear, isotropic operating resistivity and recoil permeability; no radial variation. The material defaults are illustrative NdFeB values, consistent with Pacific PAC's typical 150 µΩ·cm resistivity and relative recoil permeability 1.05.
- Each harmonic imposes the same peak radial field on all four segment edges. The supplied spectrum is maintained after cutting; different frequencies are orthogonal in the time average.
- The synchronous permanent magnet field is DC in the rotor frame and contributes no eddy loss.
- One steady mean temperature and one user-supplied thermal resistance for all rotor magnets. Resistivity is held at its entered operating value; temperature feedback and local hot spots are not solved.
Where it applies
The result is a solution of the stated 2D diffusion problem, not a universal machine loss prediction. Actual surface fields can vary with position, segment gaps and external magnetic reluctance. In particular, voltage-fed PWM and the surrounding core can change harmonic amplitudes when segmentation changes. Chen, Wang & Nair (2016), pp. 1–2 of the accepted manuscript, discuss the limitations of radial-only and reduced-dimensional assumptions. Use a coupled 3D model when radial currents, tangential harmonics, axial end fields, curvature, sleeves, proximity effects or conductive coatings matter.
Six distinct field harmonics are supported. Omitted spectral components and spatial harmonics are not synthesized from the switching frequency. Electrical bridges between segments invalidate the insulated-boundary assumption. A mean temperature below the entered limit does not establish local demagnetization margin; the limit depends on magnet grade and magnetic load line.
The default sum uses 256 odd modes per harmonic, with a doubled-term convergence comparison for the selected and uncut magnets. The integer search solves every pair, without smoothing the design map. Increase the series terms if its convergence warning appears; extreme skin depths may require more terms and longer runtime. The frequency response uses 101 solved samples at a fixed 10 mT peak field; it is a transfer comparison, not the entered spectrum. Schematic loops illustrate closed current paths rather than their local magnitudes.
Independent five-point finite-difference benchmarks verify the diffusion solution in resistance-limited and skin-limited regimes; they are numerical verification, not measured motor validation. The search assumes unchanged material volume, harmonic forcing and thermal resistance. The fewest-piece recommendation is only within the entered count and size limits; no unique finite loss-only optimum exists if arbitrarily fine, cost-free segmentation is allowed.
References
- Ruoho, Santa-Nokki, Kolehmainen & Arkkio (2009), Modeling Magnet Length in 2-D Finite-Element Analysis of Electric Machines, IEEE Trans. Magn. 45, 3114–3120. https://doi.org/10.1109/TMAG.2009.2018621
- Ede, Atallah, Jewell, Wang & Howe (2007), Effect of Axial Segmentation of Permanent Magnets on Rotor Loss in Modular Permanent-Magnet Brushless Machines, IEEE Trans. Ind. Appl. 43, 1207–1213. https://doi.org/10.1109/TIA.2007.904397
- Huang, Bettayeb, Kaczmarek & Vannier (2010), Optimization of Magnet Segmentation for Reduction of Eddy-Current Losses in Permanent Magnet Synchronous Machine, IEEE Trans. Energy Convers. 25, 381–387. https://doi.org/10.1109/TEC.2009.2036250
- Königs & Löhlein (2023), Lumped model for the calculation of harmonic eddy current losses in permanent magnets for homogeneous flux distributions considering eddy current reaction flux. https://doi.org/10.1007/s00502-023-01147-z
- Chen, Wang & Nair (2016), An analytical method for predicting 3D eddy current loss in permanent magnet machines based on generalized image theory. https://doi.org/10.1109/TMAG.2015.2500878
- Pacific PAC Technologies, NdFeB typical mechanical properties: resistivity and recoil permeability. https://www.ppactech.com/resources/ndfeb/
- Independent finite-difference benchmarks: Python 3 / SciPy 1.15.3, squared-gradient loss integration on refined 2D grids; attached verification source s5.
Inputs
| Input | Symbol | Unit | Default |
|---|---|---|---|
| Circumferential magnet width | W | mm | 30 |
| Axial magnet length | L | mm | 100 |
| Radial magnet thickness | h | mm | 5 |
| Parent magnets in rotor Count the full blocks before either segmentation cut. | N_mag | 8 | |
| Circumferential segments | N_c | 3 | |
| Axial segments | N_a | 4 | |
| Resistivity at operating temperature Use the grade's measured value at the expected magnet temperature. | rho_pm | uohm*m | 1.5 |
| Relative recoil permeability | mu_r | 1.05 | |
| Active harmonics Enter distinct rotor-frame frequencies and peak radial fields at the magnet surface. | n_h | 3 | |
| Harmonic 1 frequency Slot passing frequency is slots × mechanical revolutions per second. | f1 | Hz | 1000 |
| Harmonic 1 peak radial field Use the alternating harmonic amplitude, excluding the rotor's synchronous DC field. | B1 | mT | 20 |
| Harmonic 2 frequency | f2 | Hz | 10000 |
| Harmonic 2 peak radial field | B2 | mT | 5 |
| Harmonic 3 frequency | f3 | Hz | 20000 |
| Harmonic 3 peak radial field | B3 | mT | 2 |
| Harmonic 4 frequency | f4 | Hz | 30000 |
| Harmonic 4 peak radial field | B4 | mT | 0 |
| Harmonic 5 frequency | f5 | Hz | 40000 |
| Harmonic 5 peak radial field | B5 | mT | 0 |
| Harmonic 6 frequency | f6 | Hz | 50000 |
| Harmonic 6 peak radial field | B6 | mT | 0 |
| Rotor magnet loss budget | P_budget | W | 20 |
| Rotor heat sink temperature | T_sink | °C | 60 |
| All-magnet thermal resistance Use the thermal resistance from all rotor magnets together to the heat sink. | R_theta | K/W | 1 |
| Magnet temperature limit Use the grade's limit at the actual magnetic load line. | T_limit | °C | 120 |
| Search circumferential count up to | N_c_max | 12 | |
| Search axial count up to | N_a_max | 12 | |
| Minimum segment dimension The search excludes smaller widths or axial lengths. | d_min | mm | 2 |
| Diffusion series terms Double this to check convergence when the skin depth is small. | N_terms | 256 |
Results
With the default inputs above, the calc gives:
| Result | Value |
|---|---|
| Selected rotor magnet loss | 32.806 W |
| Unsegmented rotor magnet loss | 52.534 W |
| Loss reduction (%) | 37.55 |
| Mean magnet temperature rise | 32.81 K |
| Mean magnet temperature | 92.81 °C |
| Selected design | Selected above target |
| Candidate circumferential segments | 5 |
| Candidate axial segments | 1 |
| Candidate pieces in rotor | 40 |
| Candidate rotor magnet loss | 19.114 W |
| Candidate mean temperature | 79.11 °C |
| Search result | Fewest pieces meeting target |
| Lowest loss in search | 2.996 W |
| Lowest-loss circumferential count | 12 |
| Lowest-loss axial count | 12 |
Checked against sources
Each case below was calculated by this calc and agrees with the value from its source (within the tolerance where one is shown).
Resistance-limited finite-difference benchmark
Source: Independent Python/SciPy 1.15.3 finite-difference PDE solution; 400 × 1333 cells. Loss integrated from squared edge gradients, not the calculator's separated series. See attached numerical reference; verifies this BVP, not a measured motor.
Inputs: Active harmonics = 1; Harmonic 1 frequency = 50 Hz; Circumferential segments = 1; Axial segments = 1; the others as above.
| Result | Source | This calc |
|---|---|---|
| Selected rotor magnet loss | 0.09599 W | 0.096 W |
Skin-limited uncut finite-difference benchmark
Source: Independent refined finite-difference reference, 400 × 1333 cells; peak 4 mT at 20 kHz. Squared-gradient integration captures screened interior and skin current. Mesh doubling changes power by 0.0254%.
Inputs: Active harmonics = 1; Harmonic 1 frequency = 20,000 Hz; Harmonic 1 peak radial field = 4 mT; Circumferential segments = 1; Axial segments = 1; the others as above.
| Result | Source | This calc |
|---|---|---|
| Selected rotor magnet loss | 15.42 W | 15.421 W |
Segmentation anomaly finite-difference benchmark
Source: Independent refined finite-difference reference, 400 × 800 cells in each 10 × 20 mm segment. The 3 × 5 split increases the total loss at 20 kHz despite reducing the current-loop size. Mesh doubling changes segmented power by 0.0052%.
Inputs: Active harmonics = 1; Harmonic 1 frequency = 20,000 Hz; Harmonic 1 peak radial field = 4 mT; Axial segments = 5; the others as above.
| Result | Source | This calc |
|---|---|---|
| Selected rotor magnet loss | 33.544 W | 33.544 W |
| Unsegmented rotor magnet loss | 15.42 W | 15.421 W |
Related Resources
Related articles
-
Computational Modeling of Multi-Scale Plasma Flow and Magnetic Flux Dynamics in Tokamak System Using Matlab -
Modeling and Visualization of Dipole Magnetic Fields Using Streamlines and Runge–Kutta Method in Matlab -
Parametric Sensitivity Analysis of Permanent Magnet Synchronous Motors (PMSM) for Minimizing Electromagnetic Vibration in MATLAB
Related courses
Induction motor simulation using RMXPRT tool of ANSYS MAXWELL
IEC Electrical Schematics and Industrial Control Panels

Responses