Axial Flux Motor Design Tool

Size a coreless axial flux PM motor with radial-slice torque and EMF calculations. Compare winding layouts, copper losses, voltage demand and estimated efficiency.

#AFPM #axial-flux #concentrated-winding #coreless winding #coreless-stator #motor sizing #permanent magnets #quasi-3D

This axial flux motor design calculator estimates the electromagnetic and electrical performance of a three-phase permanent-magnet motor with a coreless stator between two magnet rotors. It is intended for engineers comparing layouts, students learning AFPM sizing, and builders assessing a direct-drive concept before detailed magnetic, thermal and mechanical analysis.

Start with the active outer diameter, inner-to-outer diameter ratio, pole count and operating speed. The active annulus excludes the housing and end-winding overhang. Choose overlapping coils or type-I concentrated coils, then enter turns per coil, bare strand diameter, strands in hand and parallel phase circuits. The tool derives the coil count and coil-side width for the selected balanced winding family. Calculated copper fill checks whether the proposed conductor can fit inside the coil pack; it is not an insulation or manufacturing specification.

For the most useful prediction, enter the peak fundamental axial flux density at the winding mid-plane from measurements or FEM. Inner and outer values can differ; the calculator interpolates between them and integrates across radial slices. These values already include both rotors. Alternatively, the ideal magnetic-circuit option estimates a field from magnet thickness, remanence, pole coverage and the nonmagnetic gap. This option omits leakage and fringing and needs independent verification.

Set either RMS phase current or a total copper-loss budget. Copper resistance uses the entered winding temperature, while an optional AC resistance factor accounts for additional winding resistance. Temperature is prescribed; the calculator does not solve a cooling system. Bearing, windage and additional rotational losses are user estimates at a reference speed. Their defaults illustrate a calculation and must be replaced when estimating a real machine's efficiency.

Read shaft output separately from developed electromagnetic torque. The required-current result compares the design with your requested shaft power. Voltage demand includes winding resistance and the supplied phase inductance, with a wye connection and ideal space-vector PWM bus limit. Warning messages flag insufficient voltage, insufficient shaft output, invalid parallel paths and excessive copper fill. The plots show current sensitivity and the contribution of each radial region; their curves do not impose a thermal rating.

The method combines radial-slice integration with published coreless winding geometry and sinusoidal PM energy conversion. Two cases compare torque and copper mass with Kamper, Wang and Rossouw's measured machines. A separate Parviainen case checks the 5 kW / 300 rpm shaft requirement only.

This AFPM sizing tool covers the stated coreless winding families. It does not predict cogging, harmonic torque, magnetic saturation, demagnetization, structural strength, bearing life, inverter losses or a continuous thermal rating. The drawing shows schematic proportions. Use the metric/US switch for supported mechanical units; electrical quantities retain their usual engineering units.

How it is calculated

The machine has one coreless three-phase wye stator and two PM rotors. This is a reduced analytical sizing model, not a reproduction of the full magnetic equivalent circuit, iron-loss or thermal model in Parviainen's dissertation.

The active annulus is divided into equal radial slices, evaluated at their midpoints, following the quasi-3-D principle described by Parviainen (2005), attached pp. 36–38. The prescribed fundamental field is interpolated linearly between inner and outer radii.

For total poles PP, overlapping windings use Q=3P/2Q=3P/2 coils and concentrated type-I windings use Q=3P/4Q=3P/4, one coil per phase group. Coil width and mean end-turn length follow Kamper et al. (2008), attached pp. 4–6: Eqs. (18)–(19), (24), (26), (30) and (32). The maximum side width can be reduced by the entered geometric width fraction. At each radius, θj=Pw/(2rj)\theta_j=Pw/(2r_j); the overlap factor is sin⁡(θj/2)/(θj/2)\sin(\theta_j/2)/(\theta_j/2). The concentrated factor also includes sin⁡[(θm−θj)/2]\sin[(\theta_m-\theta_j)/2]. This is an explicitly unskewed radial extension of their finite-width coil model; the concentrated active-side length correction is neglected.

The peak fundamental phase flux linkage is

λpk=4NsP∑jB1jkwjrjΔr,\lambda_{pk}=\frac{4N_s}{P}\sum_j B_{1j}k_{wj}r_j\Delta r,

where Ns=Nturn(Q/3)/aN_s=N_{turn}(Q/3)/a. Hence Eph=ωeλpk/2E_{ph}=\omega_e\lambda_{pk}/\sqrt2 and Tem=3(P/2)λpkIrms/2T_{em}=3(P/2)\lambda_{pk}I_{rms}/\sqrt2. The air-gap field is not multiplied by two: it already describes the field crossing the central winding. This sinusoidal power conversion is consistent with the general sizing framework of Huang et al. (1998), attached pp. 1–2. A 120-slice recalculation reports numerical convergence independently of the selected slice count.

Conductor area includes all strands in hand. With identical parallel paths,

R20=ρ20(Q/3)Nturnlturna2ACu,PCu=3Irms2RAC.R_{20}=\rho_{20}\frac{(Q/3)N_{turn}l_{turn}}{a^2A_{Cu}},\qquad P_{Cu}=3I_{rms}^2R_{AC}.

Copper temperature correction uses αCu=0.0038 K−1\alpha_{Cu}=0.0038\ \mathrm{K}^{-1} from Parviainen, attached p. 52, Eq. (2.31). Nominal annealed-copper resistivity and density are stated assumptions. The AC multiplier is supplied, not calculated from the skin-depth screening result.

Rotational loss allowances scale as speed (bearing), speed cubed (windage), and speed squared (additional eddy allowance). These are explicit user assumptions, not fitted paper correlations. Shaft power is converted power minus rotational losses, clipped to zero below drag torque. Efficiency divides that available output by converted power plus copper loss; it is defined as zero at zero electrical input. The raw shaft balance is retained for conservation checks.

Voltage uses Vph=(Eph+IRAC)2+(ωeLI)2V_{ph}=\sqrt{(E_{ph}+IR_{AC})^2+(\omega_eLI)^2} with current in phase with the PM EMF, and VLL=3VphV_{LL}=\sqrt3V_{ph}. The ideal linear SVPWM line-line RMS limit is Vdc/2V_{dc}/\sqrt2. Both the selected current and target shaft current are checked.

The optional ideal field is derived here from a one-dimensional linear magnetic circuit: B0=Br(2hm)/[2hm+μr(hw+2g)]B_0=B_r(2h_m)/[2h_m+\mu_r(h_w+2g)] and B1=(4/π)B0sin⁡(παp/2)B_1=(4/\pi)B_0\sin(\pi\alpha_p/2). It assumes rectangular alternating pole fields and is not a validated substitute for the field solution in either named paper.

Assumptions

  • Uniform round bare copper strands; all coils have equal turns and length. Copper resistivity at 20 °C is 1.724×10−8 Ω m1.724\times10^{-8}\ \Omega\,m and density is 8960 kg/m38960\ \mathrm{kg/m^3}. Leads, insulation and joints add mass and resistance in a real winding.
  • Three balanced wye phases; parallel paths are identical and connected with the correct polarity. The tool sizes two established coil-count families; it does not generate a winding connection diagram.
  • Sinusoidal fundamental field and current, aligned for torque production; no reluctance torque, field weakening or armature-reaction correction. Entered field is at operating magnet temperature and includes the effect of both rotor disks.
  • Each bare round strand must fit both the axial winding thickness and the coil-side width; the copper-area fill check also applies. These necessary geometric checks do not establish an achievable strand-packing arrangement or insulation clearance.
  • Radial active conductors are approximated as straight and unskewed. End-turn formulas are geometric approximations. Copper fill is copper area divided by gross coil-side cross-section.
  • Copper temperature, phase inductance, AC resistance factor and rotational losses are supplied design estimates. Default inductance (0.3 mH), bearing loss (8 W) and windage loss (2 W at 300 rpm) are illustrative assumptions, not measurements of the published prototypes.
  • With a known field, changing the air gap does not automatically change that field: update its inner and outer values. The gap still enters overlapping end-turn clearance. The ideal field option updates the field from gap and coil thickness.
  • No stator iron is present, so stator iron loss is zero. Rotor and winding eddy losses are not guaranteed to be zero; enter an additional allowance and/or an AC multiplier as appropriate. Inverter and supply losses are excluded.
  • The ideal field and loss scaling are derived/assumed here; they are not claimed to be Parviainen's complete quasi-3-D model. The static schematic needs no animation.

Where it applies

Use for preliminary comparison of dual-rotor coreless AFPM motors with the two supported balanced winding families. Concentrated windings require pole counts divisible by four; overlapping windings require even pole counts. Parallel circuits must divide the coils per phase. For overlapping windings the inner radius must exceed the total magnet-to-magnet nonmagnetic gap used by the end-turn approximation.

The numerical input ranges are editing limits, not validated operating envelopes. Published machine validation is at 400 mm active diameter, inner/outer ratio 0.7, 24 poles, peak fundamental field 0.527 T, and 120 W copper loss. Prediction versus measurement is 45.49 versus 43.9 N·m for overlap and 44.87 versus 44.6 N·m for concentrated coils. Copper-mass tolerances are wider because simplified end turns omit leads and detailed construction. These are limited validation points, not a universal accuracy specification.

The third reference verifies only the rated shaft requirement in Parviainen Table 3.1 (5 kW, 300 rpm, 159 N·m). The iron-cored dissertation prototype has a different topology; its measured 211 V EMF and efficiency are not claimed as predictions of this calculator. The Kamper checks use torque and copper mass, not the Table VI rated-voltage entries: its voltage labeling is not used to establish independent back-EMF accuracy.

Verify field, end turns and effective resistance for a final design. Low pole counts, severe radial field variation, thick coreless windings, large gaps, skew and unusual magnet shapes can exceed the simplified model. Refine radial slices until the reported integration difference is small; this only checks numerical integration. The ideal one-dimensional field neglects saturation, flux leakage, radial end effects and fringing.

Current sweeps are theoretical at fixed field and temperature; they do not enforce a continuous thermal rating or inverter current rating. Below drag torque, available motoring output and efficiency are zero while the raw power balance is negative. No stall solution, rotor stress, magnet retention, bearing force, thermal equilibrium, demagnetization or fault behavior is calculated.

References

  • Parviainen, A. (2005). Design of axial-flux permanent-magnet low-speed machines and performance comparison between radial-flux and axial-flux machines. LUT dissertation, Acta Universitatis Lappeenrantaensis 208. Attached pp. 36–38 (radial planes), 51–53 (resistance/losses), 73 (rated prototype). https://lutpub.lut.fi/handle/10024/31185
  • Kamper, M. J., Wang, R.-J., Rossouw, F. G. (2008). Analysis and Performance of Axial Flux Permanent-Magnet Machine With Air-Cored Nonoverlapping Concentrated Stator Windings. IEEE Transactions on Industry Applications 44(5), 1495–1504. DOI: 10.1109/TIA.2008.2002183. Author repository PDF: https://research.ee.sun.ac.za/emr/staff/rwang/papers/kamper_08.pdf . Attached pp. 3–7 (equations), 9–10 (Tables VI–VII).
  • Huang, S., Luo, J., Leonardi, F., Lipo, T. A. (1998). A General Approach to Sizing and Power Density Equations for Comparison of Electrical Machines. IEEE Transactions on Industry Applications 34(1), 92–97. DOI: 10.1109/28.658727. Open institutional copy: https://minds.wisconsin.edu/handle/1793/9930 . This is the general sizing framework; the later AFPM comparison paper is 1999, DOI: 10.1109/60.766982.

Inputs

InputSymbolUnitDefault
Winding layout
Three-phase, wye-connected coreless stator between two magnet rotors. Concentrated coils use Kamper type I, one coil per phase group.
windingConcentrated, 3 coils per 4 poles
Active outer diameter
Diameter of the active magnet / coil annulus, excluding end turns and housing.
D_outmm400
Inner / outer diameter ratio
Defines the active inner diameter. End turns extend beyond the active annulus.
radius_ratio0.7
Total magnet poles
Use an even number; concentrated layout requires a multiple of four. This is poles, not pole pairs.
poles24
Rotor speed
Positive motoring speed. Stall and field weakening are outside this model.
n_rotrpm300
Required shaft output
Used for the required torque and current comparison; it does not force the calculated output to match.
P_targetW1300
Set the operating current by
Current is in phase with the fundamental back-EMF (id = 0); no reluctance torque.
drive_modeRMS phase current
RMS phase current
Total terminal current per phase, equal to line current for wye connection.
I_setA16
Total copper loss budget
All three phases, using AC resistance at the entered copper temperature.
P_cu_budgetW120
Available DC bus voltage
Ideal linear space-vector PWM limit: line-line fundamental RMS ≤ Vdc / √2. Allow practical control margin.
V_dcV120
Phase synchronous inductance
Supply an estimate or measured value. The default is illustrative, not a prediction from geometry.
L_phasemH0.3
Turns per coil
All coils have the same turn count.
N_turn25
Parallel winding circuits
Must divide the coils per phase. Assumes identical balanced paths; strand count is a separate input.
a_paths1
Bare strand diameter
Round copper strand; enamel and inter-coil insulation are excluded. Each strand must fit both winding thickness and coil-side width; allow additional insulation clearance.
d_wiremm0.85
Strands in hand per turn
Parallel strands forming one conductor in each coil; not parallel phase circuits.
n_strand6
Axial winding thickness
Gross thickness of the coreless coil pack; used for fill and the optional field estimate.
h_coilmm8.4
Fraction of maximum coil-side width
Geometric coil-side width fraction, not copper fill factor. Reducing it increases calculated copper packing.
width_fill1
Copper temperature
Specified operating temperature. No thermal equilibrium is solved.
T_cudegC80
AC / DC resistance factor
Use 1 for negligible skin and proximity effects, or supply a measured / calculated multiplier.
k_ac1
Fundamental field source
Prefer measured or FEM fundamental field at the winding mid-plane. Ideal estimate neglects leakage, fringing and saturation.
field_modeKnown field at inner and outer radii
Peak fundamental field at inner radius
Fundamental peak axial field, not magnet remanence or an RMS value. Interpolated linearly across the annulus.
B_innerT0.527
Peak fundamental field at outer radius
Includes effects of both rotors; the field is not doubled in the torque calculation.
B_outerT0.527
Mechanical air gap on each side
Used in the ideal field estimate and overlapping end-turn geometry. With known field, update that field separately when changing the gap.
g_airmm1.8
Magnet thickness on each rotor
Two identical surface magnet layers in the ideal circuit.
h_magmm6
Magnet remanence at operating temperature
Enter remanence at magnet temperature; no automatic temperature or demagnetization model.
B_remT1.2
Magnet recoil relative permeability
Linear recoil line used only by the ideal circuit.
mu_rec1.05
Magnet arc / pole pitch
Ideal rectangular pole coverage. The fundamental uses (4/π) sin(πα/2).
alpha_p0.8
Speed of the supplied loss estimates
Reference speed for the three user-supplied rotational loss components.
n_lossrpm300
Bearing loss at reference speed
Illustrative default: constant drag torque, so loss scales with speed.
P_bearW8
Windage loss at reference speed
Illustrative default: loss scales with speed cubed. Not a disk-drag correlation.
P_windW2
Additional rotational loss at reference speed
Optional rotor / conductor eddy-loss allowance, scaling with speed squared. Default zero does not prove absence of eddy loss.
P_eddyW0
Number of radial slices
Midpoint integration. Compare 20, 40 and 80 slices when radial field or coil pitch varies strongly.
N_rad40

Results

With the default inputs above, the calc gives:

ResultValue
Available shaft output1.412 kW
Available shaft torque44.95 N·m
Developed electromagnetic torque45.27 N·m
Estimated motor efficiency89.8 %
RMS current for target14.74 A
Selected RMS phase current16 A
Copper loss at operating point150 W
Open-circuit line-line EMF51.32 V
Required line-line RMS voltage56.81 V
Line-line RMS voltage for target56.38 V
Ideal line-voltage margin28.04 V
Winding copper mass2.874 kg
Coil-side copper fill41.5 %
Electrical frequency60 Hz

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).

Kamper overlapping winding: measured torque and copper mass

Source: Kamper, Wang & Rossouw (2008), repository PDF pp. 9–10, Tables VI–VII: 400 mm diameter, radius ratio 0.7, 24 poles, 36 coils, 11 turns, 7 × 0.8 mm strands, Bp = 0.527 T, 120 W copper loss. Table VII measured torque 43.9 N m, measured copper mass 3.484 kg. Entries read visually because table cells are raster images. Tolerance 2.0 N m (4.6%) and 0.30 kg (8.6%) cover simplified end-turn length and omitted leads; no fitted coefficients. Checks developed torque, not generator shaft torque or iron-cored designs.

Inputs: Rotor speed = 200 rpm; Copper temperature = 20 °C; Set the operating current by = Copper loss budget; Bearing loss at reference speed = 0 W; Windage loss at reference speed = 0 W; Required shaft output = 900 W; Winding layout = Overlapping, 3 coils per pole pair; Turns per coil = 11; Bare strand diameter = 0.8 mm; Strands in hand per turn = 7; the others as above.

ResultSourceThis calc
Developed electromagnetic torque43.9 ± 2 N·m45.486 N·m
Winding copper mass3.484 ± 0.3 kg3.209 kg
Electrical frequency40 Hz40 Hz

Kamper concentrated winding: measured torque and copper mass

Source: Kamper, Wang & Rossouw (2008), repository PDF pp. 9–10, Tables VI–VII: same rotors as overlap case, 18 type-I concentrated coils, 25 turns, 6 × 0.85 mm strands, 120 W copper loss. Table VII measured torque 44.6 N m and copper mass 2.977 kg; 2-D FEA torque 45.13 N m. Raster table entries inspected visually. Tolerances 2.0 N m (4.5%) and 0.15 kg (5.0%) account for approximate end turns and ignored active-side skew. Uses standard copper resistivity at 20 °C without calibration.

Inputs: Rotor speed = 200 rpm; Copper temperature = 20 °C; Set the operating current by = Copper loss budget; Bearing loss at reference speed = 0 W; Windage loss at reference speed = 0 W; Required shaft output = 900 W; the others as above.

ResultSourceThis calc
Developed electromagnetic torque44.6 ± 2 N·m44.872 N·m
Winding copper mass2.977 ± 0.15 kg2.874 kg
Electrical frequency40 Hz40 Hz

Parviainen 5 kW / 300 rpm rated shaft requirement

Source: Parviainen (2005), Table 3.1, attached PDF p. 73: rated output 5 kW, speed 300 min−1, torque 159 N m. Tests only power-to-shaft-torque conversion. This tool does not reproduce the dissertation's one-rotor/two-iron-stator prototype or claim its measured 211 V back-EMF and efficiency.

Inputs: Required shaft output = 5,000 W; the others as above.

ResultSourceThis calc
Required shaft torque159 ± 0.2 N·m159.2 N·m
Mechanical angular speed31.416 rad/s31.416 rad/s

Kamper overlapping coil pitch: 180 electrical degrees

Source: Kamper, Wang & Rossouw (2008), Section II-A, attached PDF p. 3: the normal overlapping winding assumes a full electrical coil pitch θm = π rad (180°). The reference machine inputs use their overlapping prototype from Tables VI–VII; this case checks the stated winding geometry.

Inputs: Rotor speed = 200 rpm; Copper temperature = 20 °C; Set the operating current by = Copper loss budget; Bearing loss at reference speed = 0 W; Windage loss at reference speed = 0 W; Required shaft output = 900 W; Winding layout = Overlapping, 3 coils per pole pair; Turns per coil = 11; Bare strand diameter = 0.8 mm; Strands in hand per turn = 7; the others as above.

ResultSourceThis calc
Electrical coil pitch3.142 rad3.142 rad

Related Resources

Responses

Your email address will not be published. Required fields are marked *

L ading...