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.

#axial segmentation #circumferential segmentation #eddy current loss #magnetic diffusion #permanent magnets #PWM harmonics #rotor heating #skin effect

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

∇2B^=jωμσB^,B^∣∂Ω=Bpk.\nabla^2\hat B=j\omega\mu\sigma\hat B,\qquad \hat B|_{\partial\Omega}=B_{\rm pk}.

Here Ω\Omega 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 d=min⁡(a,b)d=\min(a,b), r=max⁡(a,b)/dr=\max(a,b)/d, Q=ωμσd2Q=\omega\mu\sigma d^2, An=(nπ)2A_n=(n\pi)^2, Zn=rAn+jQ/2Z_n=r\sqrt{A_n+jQ}/2, with odd nn:

D=∑n=1,3,…8AnjQAn+jQ(1−tanh⁡ZnZn).D=\sum_{n=1,3,\ldots}\frac{8}{A_n}\frac{jQ}{A_n+jQ}\left(1-\frac{\tanh Z_n}{Z_n}\right).

The volume-average normalized field is 1−D1-D. Integration by parts gives

Ph=VtotalωhBpk,h22μIm⁡Dh.P_h=\frac{V_{\rm total}\omega_h B_{{\rm pk},h}^2}{2\mu}\operatorname{Im}D_h.

This equals the volume integral of ∣J^∣2/(2σ)|\hat J|^2/(2\sigma). 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 min⁡(Pbudget,(Tlimit−Tsink)/Rθ)\min(P_{\rm budget},(T_{\rm limit}-T_{\rm sink})/R_\theta). 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 ΔT=PRθ\Delta T=P R_\theta.

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

InputSymbolUnitDefault
Circumferential magnet widthWmm30
Axial magnet lengthLmm100
Radial magnet thicknesshmm5
Parent magnets in rotor
Count the full blocks before either segmentation cut.
N_mag8
Circumferential segmentsN_c3
Axial segmentsN_a4
Resistivity at operating temperature
Use the grade's measured value at the expected magnet temperature.
rho_pmuohm*m1.5
Relative recoil permeabilitymu_r1.05
Active harmonics
Enter distinct rotor-frame frequencies and peak radial fields at the magnet surface.
n_h3
Harmonic 1 frequency
Slot passing frequency is slots × mechanical revolutions per second.
f1Hz1000
Harmonic 1 peak radial field
Use the alternating harmonic amplitude, excluding the rotor's synchronous DC field.
B1mT20
Harmonic 2 frequencyf2Hz10000
Harmonic 2 peak radial fieldB2mT5
Harmonic 3 frequencyf3Hz20000
Harmonic 3 peak radial fieldB3mT2
Harmonic 4 frequencyf4Hz30000
Harmonic 4 peak radial fieldB4mT0
Harmonic 5 frequencyf5Hz40000
Harmonic 5 peak radial fieldB5mT0
Harmonic 6 frequencyf6Hz50000
Harmonic 6 peak radial fieldB6mT0
Rotor magnet loss budgetP_budgetW20
Rotor heat sink temperatureT_sink°C60
All-magnet thermal resistance
Use the thermal resistance from all rotor magnets together to the heat sink.
R_thetaK/W1
Magnet temperature limit
Use the grade's limit at the actual magnetic load line.
T_limit°C120
Search circumferential count up toN_c_max12
Search axial count up toN_a_max12
Minimum segment dimension
The search excludes smaller widths or axial lengths.
d_minmm2
Diffusion series terms
Double this to check convergence when the skin depth is small.
N_terms256

Results

With the default inputs above, the calc gives:

ResultValue
Selected rotor magnet loss32.806 W
Unsegmented rotor magnet loss52.534 W
Loss reduction (%)37.55
Mean magnet temperature rise32.81 K
Mean magnet temperature92.81 °C
Selected designSelected above target
Candidate circumferential segments5
Candidate axial segments1
Candidate pieces in rotor40
Candidate rotor magnet loss19.114 W
Candidate mean temperature79.11 °C
Search resultFewest pieces meeting target
Lowest loss in search2.996 W
Lowest-loss circumferential count12
Lowest-loss axial count12

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.

ResultSourceThis calc
Selected rotor magnet loss0.09599 W0.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.

ResultSourceThis calc
Selected rotor magnet loss15.42 W15.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.

ResultSourceThis calc
Selected rotor magnet loss33.544 W33.544 W
Unsegmented rotor magnet loss15.42 W15.421 W

Related Resources

Responses

Deine E-Mail-Adresse wird nicht veröffentlicht. Erforderliche Felder sind mit * markiert

L ading...