A Review of Magnetostatic Field Derivation Techniques in Reluctance Motors and Possible Extensions to Segmented Design

1. Introduction

Reluctance motors are a promising family of electric machines due to their inherent properties such as a lack of rotor windings, fault tolerant driver design, and the zero need for magnets in their fundamental structure. Research into reluctance motors includes driver optimization [1], structural modification, and design optimization to counteract some of its drawbacks for general use [2,3]. Machine learning methods for design optimization are being applied as well [4,5]. A major design modification is the segmentation of the rotor and/or stator, which introduces a different set of advantages and issues to solve [6,7]. All these would benefit from reliable motor models to evaluate their performance before a laboratory test is performed. Computational methods such as boundary element analysis (BEA) and especially finite element analysis (FEA) can produce outputs close to laboratory results even when end effects are excluded as in 2D simulations. This accuracy extends to nonconventional machine designs as well [7,8]. However, they can be time intensive without access to high end simulation servers, whose results will only be valid to that particular motor design. Analysis methods, such as multi-objective optimization, may need a number of input sample data from FEA simulated outputs which would further incur significant computational burden. Sensitivity analysis may be performed beforehand to lessen this load [9]. FEA is opaque in that it may not provide obvious numerical insights between quantities such as input current and the change in flux flow per rotor position, unlike laid out mathematical expressions. Optimization that instead utilizes a surrogate mathematical model leads to less computational load, in exchange for possibly less accuracy, as exemplified in [10].

There is then an interest in investigating analytical methods derived from first principles or use iterative processes [11]. For reluctance motors, magnetic flux and inductance are necessary to obtain the machine’s phase current to torque output relationship. Therefore, finding the magnetostatic field distribution in the motor is a primary step in obtaining an analytical solution. While the differential equation solutions obtained may not necessarily be of closed form, they usually only require the motor dimensions and input parameters. Solutions in series form may be implemented quickly in software once the base equations have been derived. Compared to the need to re-simulate a new motor model for a new FEA run which would take new computing power and time, analytical and open form equations could provide faster solution times. This could be an acceptable compromise given that analytical solutions may only be suited for specific conditions. Analytical models or techniques can also provide mathematical insights into which machine geometry or design methodology has a measurable impact on the motor output under observation.

Given this, most works on reluctance motor analytical modeling techniques often lay out the mathematical techniques to be used then focus on the results of such methods when applied to the reluctance machine design under test. However, there appears to be limited clarification or discussion on how exactly those techniques originate from first principles such as Maxwell’s equations and what assumptions are necessary to make the solution more tractable and realistic for motor designs. It is not often clear whether the initial equations should be taken for granted or why certain assumptions were made that resulted in such a method of derivation of the magnetostatic field. A more nuanced discussion of magnetostatic field derivation techniques can open possible avenues of attack in developing magnetic field analysis techniques for reluctance machines.

Thus, the purpose of this work is to examine major analytical techniques and assumptions used to obtain the magnetic field distribution in conventional reluctance machines and to see how they could be extended for segmented motor designs. This work defines an analytic solution or derivation as one that does not require an iterative process or discretization of the machine (e.g. as in FEA or in some magnetic equivalent circuit methods) and provides a set of equations that relate the motor dimensions to output quantities such as phase inductance or static torque. They do not rely on equations that are determined ex post facto by using datapoints from prior computer simulations, as in interpolation [12,13] or design-of-experiment based methods [9,14]. Infinite series can still be involved as they usually arise in differential equations and their Fourier solutions, but the overall equations that solve the magnetic field of the motor can be immediately established.

Three magnetostatic analytical methods in solving magnetic field B around reluctance motors with conventional structures are examined. In this work, they are referred to as the air gap permeance method, magnetic potentials method, and the conformal mapping technique and will be discussed in Sections 2 to 4, respectively. The implementation of these techniques in solving the magnetic field in reluctance motors in published literature will be discussed along with their underlying assumptions. Their common traits will be examined, as it could be possible to use insights from one technique to support another. A summary of these techniques is presented in Section 6. Other related techniques are presented in Section 7. Finally, Section 8 will discuss possible ways to apply those analytical methods for segmented reluctance motor designs. This segmentation introduces a set of boundary conditions that are not necessarily the same as in machines with conventionally continuous designs and so presents a further area of analytical study for reluctance motors.

2. Air Gap Permeance

2.1. Theoretical Background

Air gap permeance uses the geometric properties of the air gap around the rotor to approximate the path to be taken by the magnetic flux in the air gap region of a machine. In basic magnetic field or circuit calculations, magnetic fields across air gaps are assumed to be uniform and thus no fringing effects are considered. However, for reluctance motors this is not necessarily always a valid assumption especially when the rotor salient poles are too far from the active stator pole. Therefore, the problem of characterizing how flux jumps through an air gap considering fringing is a central issue when deriving analytical models for reluctance motors.

Various incarnations of Ampere’s circuital law are then invoked to obtain the air gap magnetic field as a function of position around the motor. Often, very high or infinite iron permeance is also assumed so that all energy is stored in the air gaps of the machine. As will later be presented, some popular definitions used in published works of the quantity referred to as permeance are inconsistent with the textbook definition as relates to magnetic circuit theory.

In a previous work [15], it was asserted that the shape and slotting of the air gap due to the salient poles around synchronous machines is the source of torque, force, and noise production In particular, it is claimed to provide a novel Fourier series solution to the variation in air gap length around a motor in a doubly slotted machine which essentially what a conventional switched reluctance motor (SRM) is. This slotting notation is shown in Figure 1. This series representation is common when solving Maxwell’s equations in SRMs as they usually naturally appear as solutions to differential equations, with specific assumed boundary conditions. No results for the magnetic field or torque production were provided in [15].

Assuming that the iron is infinitely permeable, it is presented that the magnetic field B(γ) (where γ is the called the “peripheral position around the air gap” in radians) is provided to be of the form:

$$B\left(\gamma \right)={\mu }_0{\displaystyle \frac{MMF\left(\gamma \right)}{g\left(\gamma \right)}}={\mu }_{0}MMF\left(\gamma \right)p\left(\gamma \right)$$

(1)

This equation and its variations appear to be the core equation of works that use air gap permeance to solve for the magnetic field. In fact, p(γ) is the quantity called permeance in some definitions. This equation is related to the magnetic field of an idealized solenoid with radius R and of length L with L≫R in vacuum: $B={\mu }_0{\displaystyle \frac{NI}{L}}$. From another perspective, the equation can be visualized as infinitesimally small solenoids (whose magnetic fields are along their long axis), distributed radially around the motor whose lengths are equal to the air gap distances due to the slotting and rotor position. However, unlike that idealized solenoid visualization, the magnetic field at one section of the machine might be affected by the magnetic field in a different section of the machine. One example cause of this dependence is the position of the rotor machine, which would shorten the path while elongating another depending on the rotor pole’s location. Thus, the relationship with respect to γ must take this dependence into consideration.

At its root, the solenoid equation comes from Ampere’s law which has the following form in magnetostatic conditions:

$$\oint \mathit{B}·dl={\mu }_{0}NI$$

(2)

Here, B is constant along a chosen path and the reluctance from the air gaps with permeability ${\mu }_0$ dominate; the dl term separates into g(γ) as in (1). In effect, most of the MMF is assumed to drop in the air gap regions. Also, the work done by the magnetic field from which torque may be further claimed and shown to be of the form:

$$W_{mag}={\displaystyle \frac{{\mu }_0{R}_{m}l}{2}}{\int }_0^{2\pi }{\left[MMF\left(\gamma \right)\right]}^{2}p\left(\gamma \right)d\gamma$$

(3)

where $R_m$ is the mean radius to air gap and l is the (longitudinal) length of the air gap. Using either input-output power considerations or co-energy analysis, the basic equation for torque production in an SRM can be shown to be dependent on the rate of change of the phase inductances per rotor position θ and the square of the current i, so that excitation polarity does not matter in terms of the sign of the torque T [16]:

$$T={\displaystyle \frac{1}{2}}{i}^2\left(t\right){\displaystyle \frac{dL\left(\theta ,t\right)}{d\theta }}$$

(4)

In (1), the MMF is the distribution of the current excitation from the coils and g(γ) is referred to as the “effective air gap length for a doubly-slotted revolving air gap” by the original work. Of interest then is the Fourier series representation of p(γ), the reciprocal of the air gap length distribution. The solution to this is by first finding g(γ), which is done by finding the depth of the stator and rotor slots from an imaginary smooth cylindrical surface with radius $R_m$.

The resulting Fourier series expansions depths can then be added to produce the doubly slotted air gap series from which p(γ) is derived. The method of solving for p(γ) lays a foundation upon which most other analyses regarding magnetic flux path length in the air gap are based upon and thus been referenced in other material.

2.2. Using Air Gap Permeance for Reluctance Machine Analysis

In a similar approach [17], the flux linkage and back emf is solved from the magnetic field expression obtained for a 10 stator-8 rotor pole field excited flux-switching machine at no-load. A difference in this particular application from [15] is the assumption of infinitely deep stator slots and circular paths for magnetic fields from the sides of the slots, so that the net length of the magnetic field is longer than what is considered previously. Figure 2 shows their assumed paths. The equation used for the magnetic field in the radial direction conceptually the same as (1):

$$B_r(\upsilon ,\theta )={\mu }_{0}MM{F}_{exc}(\upsilon ,\theta )\Lambda (\upsilon ,\theta )$$

(5)

The magnetomotive force is then further provided as:

$$MM{F}_{exc}\left(\nu ,\theta \right)~{\kappa }_c\left(\nu ,\theta \right)+\sum _{n=\mathrm{1,3},5\dots }^{\infty }{b}_n\mathrm{s}\mathrm{i}\mathrm{n}\left(K\nu \right)$$

(6)

where $b_n$ and K are constants for the Fourier series of the MMF distribution function shown in Figure 3. The quantity κ_c is claimed to be the rotor scalar magnetic potential of the form:

$${\kappa }_c\left(\nu ,\theta \right)=-{\displaystyle \frac{{\int }_0^{2\pi }MM{F}_{exc}\left(\nu ,\theta \right)\mathsf{\Lambda }\left(\nu ,\theta \right)d\nu }{{\int }_0^{2\pi }\mathsf{\Lambda }\left(\nu ,\theta \right)d\nu }}$$

(7)

The concept of magnetic scalar potential will be discussed later, but for now it refers to the magnetic field distribution in a volume without current carrying conductors (NI = 0). The purpose of κ_c as claimed in the work is to act as a compensation so that ∇∙B=0. Referring to Figure 3, it serves to effectively cancel the offset of the MMF function so that the average value for the entire motor surface is zero.

In (5) and (7), Λ(υ,θ) is the inverse of the air gap distance distribution between the rotor and stator with the aforementioned circular path assumption included. This is thus effectively p(γ) in (1). υ and θ are the angular position around the air gap and the rotational position of the rotor, respectively. It has the form:

$$\Lambda \left(\upsilon ,\theta \right)={\displaystyle \raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$g\left(\upsilon ,\theta \right))$}\right.}$$

(8)

where g(υ,θ) considers the stator and rotor slotting along with the nominal air gap length e:

$$g\left(\upsilon ,\theta \right)=e+{\delta }_s\left(\upsilon \right)+{\delta }_r(\upsilon ,\theta )$$

(9)

The choice for the assumed flux path (in Figure 2) in the motor slots may be vindicated especially when a rotor pole is near aligned with the stator poles, given that most of the flux especially in the aligned position is near the tips of the poles as shown in Figure 4. Λ(υ,θ) is referred to as the air-gap permeance in [3]. Given that this air gap permeance has units of inverse length $lengt{h}^{-1}$, it deviates from the typical definition of permeance in electromagnetic theory as the inverse of reluctance $mas{s}^{-1}lengt{h}^{-2}tim{e}^{2}curren{t}^2$ or ${\displaystyle \frac{length}{{\mu }_0{\mu }_{r}area}}$) [18,19].

In Figure 2, the flux paths from the two slot sides along angular position $\nu$ are assumed to follow ${\Gamma }_1$ and ${\Gamma }_2$which are quarter circular paths with lengths:

$$\left\{\begin{array}{c}{\Gamma }_1=\pi R_{si}sin\left({\displaystyle \frac{\nu }{2}}\right)\\ {\Gamma }_2=\pi R_{si}sin({\displaystyle \frac{{\beta }_s{\theta }_s}{2}}-{\displaystyle \frac{\nu }{2}})\end{array}\right.$$

(10)

It is then taken that flux flows through both these paths, such that they are in parallel ${\Gamma }_1\left|\right|{\Gamma }_2$. Hence the effective stator path length ${\delta }_s\left(\upsilon \right)$ is computed to be:

$${\delta }_s\left(\nu \right)={\displaystyle \frac{\pi R_{si}}{2}}{\displaystyle \frac{\mathrm{s}\mathrm{i}\mathrm{n}\left(\frac{\nu }{2}\right)\mathrm{s}\mathrm{i}\mathrm{n}(\frac{{\beta }_s{\theta }_s}{2}-\frac{\nu }{2})}{\mathrm{s}\mathrm{i}\mathrm{n}\left(\frac{{\beta }_s{\theta }_s}{4}\right)cos(\frac{{\beta }_s{\theta }_s}{4}-\frac{\nu }{4})}}$$

(11)

where rotor path length ${\delta }_r\left(\upsilon \right)$ is derived in the same manner, completing equation (9). There appears to be an inherent assumption with this parallel path approach in this stator slot air gap calculation that the flux is coming or going to both walls of the infinitely deep stator slot in Figure 2. This may be true when the rotor pole is unaligned with a stator pole, resulting in the flux going through both ${\Gamma }_1$ or ${\Gamma }_2$, but this may not be applicable when the machine is either in the unaligned position or when only one stator pole currently has an energized coil. For example, in Figure 4c where the rotor of a 6-4 SRM is at 30° out of aligned position, the paths already appear to be elliptical. Granted, this may be an extreme case due to the low number of stator-rotor poles resulting in huge air gaps when the rotor is close to the unaligned position. The circular path assumption could be acceptable for machines such as the 10-8 reluctance motor that this analysis was originally applied upon [17]. The parallel equation does allow the shorter path to dominate so this can be enough to counteract this slight inaccuracy. Overall, this approach of finding the effective path lends itself to other possible formulations or even simplifications for different slot shapes. As an example, are the paths followed by the flux further away from the air gap still circular? Or would they elongate or form elliptical arcs of various lengths? By combining with other techniques such as insights from FEA results or finding the magnetic potentials in these regions, better approximations for path shapes may be obtained.

In an approach whose authors call the air gap magnetomotive force permeance model [21], the MMF F(θ,t) is composed of winding functions $N_X\left(\theta \right)$ (sometimes also called the turns function) from a certain phase, expressed as Fourier sums [22]. Equation (12) shows this relationship. The winding currents from coils A, B, and C are also considered in the total excitation of the motor. Thus, the effect of multiphase excitation is more naturally included in the rotating rotor analyses.

$$F\left(\theta ,t\right)=N_A\left(\theta \right)I_A\left(t\right)+N_B\left(\theta \right)I_B\left(t\right)+N_C\left(\theta \right)I_C\left(t\right)$$

(12)

This work has a slightly different definition of permeance was also obtained from [23]:

$${\mathsf{\Lambda }}_{total}(\theta ,t)={\displaystyle \frac{g_0·{\mathsf{\Lambda }}_s\left(\theta \right)·{\mathsf{\Lambda }}_r(\theta ,t)}{{\mu }_0}}$$

(13)

where t = time (relates rotor speed and thus angular position) and θ = spatial position around the air gap. $g_0$ in (12) is the nominal air gap distance of the machine. ${\Lambda }_s$and ${\Lambda }_r$ are the permeances modeled using the approach similar to both [15] and [17]. Given this similar definition of Λ, it is seen that the quantity used here is a scaling factor to $g_0$instead of being additive quantity as in (9). Other works such as [24] use this multiplicative definition as well. Observe that the computed the effective path length for ${\Lambda }_s$and ${\Lambda }_r$ are the same as in (10) and (11). The magnetic field is then:

$$B_g\left(\theta ,t\right)=F\left(\theta ,t\right)·{\mathsf{\Lambda }}_{total}\left(\theta ,t\right).$$

(14)

Finding the total flux is a relatively simple integration of the B along the chosen area around the air gap. However, the presented torque equation is provided is not as straightforward:

$$T_e\left(t\right)=-l_{stack}{R}_{si}{\int }_0^{2\pi }{B}_g\left(\theta ,t\right)d{F}_s(\theta ,t)$$

(15)

$$F_s\left(\theta ,t\right)=F\left(\theta ,t\right)·{\mathsf{\Lambda }}_s\left(\theta \right)·{\displaystyle \frac{g_0}{{\mu }_0}}.$$

(16)

The integration is in terms of the stator MMF through $d{F}_s(\theta ,t)$, instead of the spatial position around the rotor, for example. Detailed derivation is found on [25], which starts from the usage of a stator electrical loading quantity $K_s\left(\theta ,t\right)$ and then applying Lorentz force law. Specific electrical loading is usually defined as the density of the current carrying conductors around the periphery of the electric machine whose SI units are in Ampere-turns/m. [16]. $K_s\left(\theta ,t\right)$ is then provided as having a form:

$$K_s\left(\theta ,t\right)={\displaystyle \frac{1}{R_{si}}}{\displaystyle \frac{dF\left(\theta ,t\right)}{d\theta }}.$$

(17)

Thus, it can be shown that (18) follows, which in turn simplifies to (15) when the dθ terms cancel.

$$T_e\left(t\right)=-R_{si}{L}_{stk}{\int }_0^{2\pi }{B}_g\left(\theta ,t\right)K_s\left(\theta ,t\right)R_{si}d\theta .$$

(18)

Figure 5 shows the predicted torque of [21] which applies equations (12) to (18). There appears to be a maximum error of around 3 Nm at certain rotor positions, while the averages differ by 1.14 Nm, representing a +4.17% error from the simulated value.

In the 6-4 SRM example, it can be supposed that semicircular paths for the flux may be a safe assumption for small deviations from the aligned rotor position, depending on the stator-rotor pole relationship. An investigation into determining at which angle the errors start to become significant, and what a significant error means in the first place, may be of interest. For Figure 4c, it appears that at even relatively small angles the paths have already become elliptical. Permeance computations using elliptical paths show that the eccentricity of the path is dependent on the air gap distance, as verified in FEA calculations [26]. Referring to Figure 6, it was specifically found that in an SRM with an air gap distance of $l_g$, the semimajor axis is a multiple kt of the minor axis of length t. Moreover, k is a nonlinear function only of the ratio $t/l_g$.

The root equation (1) implicitly presumes that the path length L modulates the field around the motor. These techniques are reminiscent of the derivation of the Carter coefficients [27,28], although Carter initially used conformal mapping techniques that produce more justifiable flux paths instead of rather arbitrary but easier to handle choices discussed here. Carter’s coefficients are used to scale the base magnetic field value at the appropriate locations around the motor by considering the effect of the slotting of the machine in the flux flow, much like air gap permeances solved above. While it is not technically a flux path-based parameter, finding air gap lengths are certainly important considerations for the coefficients’ derivation [29].

In general, air gap permeance techniques treat permeance as a quantity that considers the effect of the slotting of the machine, which also affects the magnitude of the magnetic field. The initial use of Ampere’s law itself requires strong assumptions regarding the chosen closed loop around the current carrying conductor, but authors do claim good agreement with simulated results Due to the already defined concept of magnetic permeance, it might be more appropriate to call the quantities p(γ) or ${\Lambda }_{total}$ an air gap modulator or coefficient, since they adjust how the computed MMF flows around the assumed paths in the motor [21,30]. Despite nomenclature inconsistencies and path simplifications, the works claim small errors with the predicted outputs such as torque production compared to FEA results.

3. Magnetic Potentials

3.1. Theoretical Background of Magnetic Potentials

Magnetic potential is analogous to electric potential and arises from the scalar and vector fields embedded in Maxwell’s equations. The magnetic vector potential’s (represented by A) existence in a region with non-zero current flow is guaranteed by Gauss’s law, Helmholtz’s decomposition theorem and the fact that the divergence of a curl of a vector field is zero:

$$\nabla ·\mathit{B}=0=\nabla ·\left(\nabla \times \mathit{A}\right).$$

(19)

The magnetic scalar potential, on the other hand, exists if there is no current source within the volume discussed [31]. While not as commonly utilized compared to the magnetic vector potential [32], it is useful for determining the magnetic field in zero-current volumes such as around a permanent magnet or, in this discussion, in the air gap regions of reluctance motors. Aside from stator-rotor air gap regions, the magnetic scalar potential has been used to design flux barriers in synchronous reluctance machines as well [33]. For a volume in air without excitation current J:

$$\mathit{B}={\mu }_0\mathit{H}$$

(20)

$$\nabla \times \mathit{H}=\mathit{J}=0={\mu }_0\left(\nabla \times \mathit{B}\right).$$

(21)

These equations are only valid for magnetostatic fields. From the fact that the curl of a gradient of a conservative field is zero in a simply connected region or volume, some corresponding magnetic scalar potential φ in this zero-current region must exist [32,34] (or, in other words if the curl of a vector field is zero, then it is a gradient of a conservative potential field):

$$\mathit{H}=-\nabla \mathsf{\phi }.$$

(22)

Simple connectedness in a region D implies that there are no holes in D; that is, every closed curve entirely within D only encloses points that are also included in D. The rotor slot, stator slot, and air gap regions as shown in Figure 7 are indeed simply connected and the magnetic scalar potential equations described above should apply. It is a sufficient and necessary condition to guarantee the existence of such a scalar potential (for completeness’s sake, it can be further pointed out that the statement “a vector field with zero divergence can be written as a curl of a potential field” is more strictly true only if the so-called second de Rahm cohomology of the vector field vanishes, but that is beyond the scope of this work and does not apply to real life motor geometries. Simple connectedness should suffice for our purposes). If the entire air gap of the machine is used instead of piecewise as in [20], then there will be a rotor shaped hole in the region assuming that the iron is infinitely permeable (no drop in the MMF inside the iron sections). To make the model as accurate as possible, the entire air gap region, including the slots, must be accounted for, so the borders of the chosen subregions must be well considered. The magnetic potentials as used relate to the MMF around the motor. The units of φ can be shown to be in ampere-turns [32]. Thus, the use of the MMF from the surface used as boundary conditions for the differential equations is justified.

In cylindrical coordinates, the magnetic scalar is then found to be of the form in (23) assuming uniformity in the z-direction, as safe assumption for motors in general especially if end effects are excluded:

$$\mathit{H}=-\nabla \mathsf{\phi }=-\left({\displaystyle \frac{\partial \mathsf{\phi }}{\partial r}}{\mathit{a}}_{\mathit{r}}+{\displaystyle \frac{1}{r}}{\displaystyle \frac{\partial \mathsf{\phi }}{\partial }}{\mathit{a}}_{\mathit{\theta }}\right).$$

(23)

This new definition can be applied to (20) and (21) and get the Laplace equation needed to be solved to obtain the magnetic flux density B from the magnetic scalar potential around the motor:

$${\displaystyle \frac{1}{{\mu }_0}}\nabla ·\left({\mu }_0\nabla \mathit{\phi }\right)={\nabla }^2\mathsf{\phi }={\displaystyle \frac{\partial \mathsf{\phi }}{\partial {\mathrm{r}}^2}}+{\displaystyle \frac{1}{r}}{\displaystyle \frac{\partial \mathsf{\phi }}{\partial \mathrm{r}}}+{\displaystyle \frac{1}{{\mathrm{r}}^2}}{\displaystyle \frac{{\partial }^2\mathsf{\phi }}{\partial {\mathsf{\theta }}^2}}.$$

(24)

Similar logic is applied to obtain B from the magnetic vector potential, the difference being that the result is a Poisson Equation due to the non-zero current density ${\mathit{J}}_{\mathit{z}}$ enclosed by the region.

$$\mathit{B}=\nabla \times \mathit{A}={\displaystyle \frac{1}{\mathrm{r}}}{\displaystyle \frac{\partial {\mathrm{A}}_{\mathrm{z}}}{\partial \mathsf{\theta }}}{\mathit{a}}_{\mathit{r}}-{\displaystyle \frac{\partial {\mathrm{A}}_{\mathrm{z}}}{\partial \mathrm{r}}}{\mathit{a}}_{\mathit{\theta }}$$

(25)

$${\mu }_0\nabla \times \left({\displaystyle \frac{1}{{\mu }_0}}\nabla \times \mathit{A}\right)={\nabla }^2{A}_z={\displaystyle \frac{{\partial }^2{\mathrm{A}}_{\mathrm{z}}}{\partial {\mathrm{r}}^2}}+{\displaystyle \frac{1}{r}}{\displaystyle \frac{\partial {\mathrm{A}}_{\mathrm{z}}}{\partial \mathrm{r}}}+{\displaystyle \frac{1}{r^2}}{\displaystyle \frac{{\partial }^2{\mathrm{A}}_{\mathrm{z}}}{\partial {\mathsf{\theta }}^2}}=-\mu {\mathit{J}}_{\mathit{z}}.$$

(26)

The only assumptions applied so far are those regarding the 2D magnetostatic scenario and infinitely permeable iron. The equations are not yet simplified to ease the derivation of the final expression. The realization here is that magnetic potential can be related to the magnetomotive force provided by the windings which are more readily known or designed. At the same time, it is based solely on Maxwell’s equations; the only assumptions here would come from how the boundary conditions of the MMF are imposed, e.g. magnetic potential is zero at the surface and inside iron sections. There is no need to assume certain flux path shapes, as this will naturally arise from the solution to the differential equations and the boundary conditions chosen. Neglect of end effects and assumption of uniform z-axis B distribution is taken for granted. Thus, this approach enables one to get a more theoretically grounded solution for B given only the MMF distribution around the machine, with reasonable assumptions for the boundary conditions.

The flux linkage of a phase coil with N number of turns is found by simply applying the classic definition of linkage (number of turns multiplied by the flux in a loop of a coil), as seen in (27). Care must be taken to properly include either $H_r$ or the $H_{\theta }$component in computing the total flux around the air gap; which one depends on the specific method. (27) explicitly excludes tangential components and so can affect the accuracy. Inductance at a given a current I is then simply the ratio ${\displaystyle \frac{\lambda }{I}}$. There are also other works [26] that have a permeance P based definition of the form $L=NP$, where N is the number of turns of the coil.

$$\lambda =2N{\mu }_0{R}_{rout}{\int }_{{\theta }_1}^{{\theta }_2}H\left(r,\theta \right)d\theta .$$

(27)

The solutions to the homogeneous Laplace’s equations are derived using separation of variables and is of the form:

$$\phi \left(r,\theta \right)=\sum _{n=1}^{\infty }{(A}_n{r}^n+B_n{r}^{-n}\left)\right(C_n\cos m_n\theta +D_n\sin m_n\theta )$$

(28)

where the constants $A_n$ to $D_n$ are to be determined from the assumed boundary conditions. The summation implicitly takes care of Fourier series usage to express the boundary conditions.

3.2. Applying Magnetic Potentials to Reluctance Motors

The three regions assumed for the 4-2 switched reluctance motor is shown in Figure 7. This general approach can be thought of as an expansion of the air gap analysis in Section 2 in that it also includes the stator and rotor slot areas in the analysis. Instead of assuming paths of different shapes and lengths within the stator and rotor slots, the magnetic field distribution naturally arises from the equations and chosen boundary conditions.

In assuming near zero reluctance of the iron core sections, it is implied that the iron surfaces have zero magnetic scalar potential $\phi =0$ if no magnetic flux comes out from that surface. Handling surfaces with magnetic flux flows on the other hand requires analysis. It is also put forward that the magnetic scalar potential can be derived from the MMF distribution around the motor. Whether these assumptions are supported or validated by theory can be a subject of future studies or at least be justified with FEA results.

A major simplification is due to the nature of the chosen general solution to Laplace’s and Poisson’s equations. The separation of variables used to gain (28) eases derivations if the boundary conditions are also separated in terms of r and θ. Thus, calculations will be relatively simpler if the walls of the stator and rotor slots and the current conducting coils are assumed to be radial from the center, because this will remove the θ-dependence of the boundaries. Observing Figure 7, the stator slot walls are not directed radially towards the center of the motor and instead are parallel to each other. Thus, if cylindrical coordinates are used, the equation for that wall is not solely in terms of r and θ because they will vary simultaneously in both coordinates. This of course is a common motor geometry as it eases the insertion of windings. This is also mostly true for the rotor slots of conventional SRMs although in this case the rotor slot walls are indeed directed radially. Thus, for this approach, some distortion of the stator slot region is required to obtain a tenable solution.

A similar adjustment must be applied to the current carrying coils or windings while maintaining the total current flow. Coils are commonly modeled in 2D as rectangular boundaries with area $a_{coil}$ with some corresponding current area density as in Figure 7. However, given the (r,θ) dependence of the solutions, solving equations (25) and (26) is simplified by thus distorting the coils while preserving the total current. It can reshaped to be a sector of a circle with the same height as the previous coil model height; this way, the region with nonzero current is bounded simply between some constant angles ${\theta }_1\le \alpha \le {\theta }_2$ and constant radii $r_1\le r\le r_2$. With the original undistorted current density scaled such that it is $J_m$ when the motor radius is 1meter, the following current distribution is then imposed:

$$J\left(r,\theta \right)=\left\{\begin{array}{c}{\displaystyle \frac{J_m}{r}},{\theta }_1<\theta <{\theta }_2\\ 0,otherwise.\end{array}\right.$$

(29)

This is used so that the current density linearly varies as with the radius. Using the angle as a parameter may result in the involvement of trigonometric functions (see (27)) in the downstream solutions and so is not preferred. The inductance of an active phase with current I is then obtained from its classic definition, using the flux linkages found for each air region:

$$L={\displaystyle \frac{{\lambda }_{air}+{\lambda }_{rotorslot}+{\lambda }_{statorslot}}{I}}.$$

(30)

All these assumptions regarding the solutions and motor geometry are claimed to result in a fairly accurate reproduction of FEA results in [35] as shown in Figure 8. It is expected that the least accurate result would be when the rotor pole face does not fully line up with the stator pole face, resulting in more flux flowing through the slots instead of the more predictable short air gap between the pole faces. The assumed modifications to the motor geometry are merely mathematical devices to untangle the differential equation solutions for cylindrical coordinates and so are not strictly necessary. Without them, the simultaneous variation in the r and θ coordinates of the motor slots and conductors would be more challenging to obtain. It is then worth pointing out that despite those changes, FEA results are claimed to be reproduced quite well.

The results of magnetic potential analysis and air gap permeances are based on the geometry of the motor. Despite the significant assumptions or spatial distortions of the regions involved in the solution, results seem to be faithful to simulated outputs. As the previous derivation shown, analyses that simplify or distort motor geometry still produce solutions that approximate numerical results. This poses a potential problem of determining which reluctance motor dimension introduces the most significant contribution to the output. From the plot in Figure 8, it seems that the method that can reproduce the most accurate answers at the extrema positions might be the most useful, since inductances in the intermediate rotor position are already reproduced well. The solutions produced by the magnetic potentials approach rely heavily on the chosen boundary conditions. Magnetic potentials has fewer base assumptions necessary compared to the air gap length approach, though it does need some geometry distortion to ease out results. Also, given that this approach includes slot regions, it is of interest to how this could be applied to other reluctance motor topologies such as segmented stator and/or rotor machines. Given that segmentation, the bottom of the slots will only have air that can be seen as extending to infinity. What would the boundary condition for that then be?

4. Conformal Mapping

4.1. Theoretical Background of Conformal Mapping

Inherent to discussions in Section 2 and Section 3 are some assumptions regarding the paths or shape of magnetic flux and/or potentials. These are adopted due to the complicated shape of salient poles of electric motors combined with the variation with rotor position, especially when the rotor and pole faces are unaligned. While these appear to produce analytical results with acceptably small errors compared to FEA, the effect of such deviations from the theoretical path might be minimized by applying conformal mapping techniques.

Conformal mapping is the transformation between functions in the complex plane that preserves the angles between two vectors while not necessarily preserving length. As an example, in Section 3, complicated boundaries introduced by the stator and rotor slots result in difficult-to-solve potential fields, which prompted some design simplifications to ease solvability. Through conformal mapping, these boundaries are transformable to a different geometry with a much easier to solve potential field (e.g. parallel plates separated by a distance) using the appropriate transform equations [36,37], and vice versa.

The basis of the SC transformation is that the upper half of the complex plane can be mapped into the interior Γ (different variable from Section 2) of polygon P with interior angles $a_1\pi$ to $a_n\pi$ using the conformal map f(z):

$$f\left(z\right)=A+C{\int }^z\prod _{k=1}^{n-1}{\left(\zeta -z_k\right)}^{a_k-1}d\zeta .$$

(31)

P has vertices $w_1$ to $w_n$. The result will be such that $w_k=f\left(z_k\right)$ for k = 1,..n-1 [37]. f(z) takes the locus of points in P as input and gives out a different set of points, hence a transformation. Obtaining the values of A and C may not be trivial and may depend on the shape and boundary conditions of the problem. (31) partially stems from the fact that a multiplication of a complex position vector by another complex number in the complex plane can be imagined to result in a rotation and/or scaling of that vector; thus a locus of such points can be transformed to a different curve with the proper complex number. This is why (31) is a product $\Pi$ of several terms.

The derivative of f(z) in (31) is equal to:

$${\displaystyle \frac{df\left(z\right)}{d\zeta }}=C{\left(\zeta -z_k\right)}^{a_k-1}.$$

(32)

Given that integration is usually more difficult than derivatives, this form could be easier to manipulate especially when the derivative in the left side is easy to obtain from the chosen transformation; in fact, this is done in [20]. For example, in Figure 9, $B_s{C}_s$ and $A_s{D}_s$ in the S-plane are parallel and circular arcs with different radii. By using the function $\mathit{ln}\left(z\right)$ as the transformation from S-plane to Z-plane, those arcs are transformed to straight lines $B_z{C}_z$ and $A_z{D}_z$since they are a constant distance from the center in the S-plane[20,40]. In the S-plane, therefore the result of using this transformation equation is:

$${\displaystyle \frac{dz}{dS}}={\displaystyle \frac{1}{S}}.$$

(33)

This can then be substituted back in (32) along with the properly determined vertices $z_s$.

4.2. Conformal Mapping in Reluctance Motor Analysis

Conformal mapping is not a machine design tool per se but is instead a mathematical device arising from the properties of the complex plane. This tool is suited in solving electric motor fields because the slotting and air gap geometry of most reluctance motors can be somewhat simple. An early application of conformal mapping in electric motors is in the original derivation of Carter coefficients [27]. The Schwarz-Christoffel (SC) transformation is commonly used and has been applied to electric motors as well [38,39,40]. It can be used to quantify the saliency [41] or slotting of machines in parallel using Carter coefficients [42]. An air gap use-case not utilized for slotting considerations is in the design of flux barriers for synchronous reluctance machines [43]. Of note in [43] is the usage of circular geometry for the barriers to take advantage of the relatively simple mapping of circular sections as described in (30) and (31). Flux barriers in synchronous reluctance motors can be more complex than simple slotting in the switched versions but are analytically solvable with this method. Another example of its application is performed during an extension of the work in [35] to include SC transformations to better compute for the magnetic potentials in the chosen regions.

The usage of magnetic potentials with and without conformal mapping has been tried in [35] and [20], respectively, using the same 6-4 motor design. The conformal transformation process is shown in Figure 9. The inductance results of the transformed motor in [20] are shown in Figure 10 which can be compared to the result from [35] that wa s shown Figure 8. It can be observed that their conformal mapping implementation reasonably approximates the static inductance of the machine as in the previous one without any conformal mapping. It is further claimed that in terms of percentage errors, the conformally mapped solution produced an error within approximately 1 to 8% of FEA benchmark values compared to the previous work’s range of 7-25%.

The utility of the SC transformation in magnetic potential analysis is that it can be proved that once the magnetic field B is solved in the transformed polygon, there is no need to reverse-transform the results back to the original plane [20,40]. When the correct coordinates and transformed boundaries are considered, the results obtained are numerically equivalent, and the magnetic potential values for the original geometry are found as well. The exact forms for $f_k$ are determined from the chosen distribution of the vertices of the original polygon to the half plane, and so careful analysis of how to distribute the vertices of the old polygon to the new region must be performed. As can be seen in the SC equation above, the resulting integrals or derivatives may not be trivial to solve for, and numerical methods may be necessary. With careful consideration to the geometry as in [20], the transformation f(z) can have a form and solution that is still analytic.

The choices of transformation to be made are rather arbitrary, but the choice should ease the derivation for (29) and (30). This choice of half plane vertex locations is important so as not to distort shape, so to speak. For example, if the side of a polygon is a straight line, the choice of $f_k$ for that side must be such that the transformed shape is also straight in the new coordinate system. This can be seen again in Figure 9. Given the possible difficulty of obtaining the integral in (29), this method may cease being analytic if numerical techniques are involved in solving the integral.

Conformal transformation, especially the SC transform, presents a promising approach that combines mathematical rigor with minimized simplifications regarding the geometry and flux path of the machine. Applying finite magnetic permeability values could be an additional step in its application as well. There has also been research that combines this with magnetic equivalent circuit analysis for synchronous reluctance motors [44]. In fact, given its nature it could be applicable to any technique that utilizes magnetic potentials as a step in its derivation.

6. Summary of Major Techniques

Table 1 shows a summary of the major magnetic field solution techniques presented in this work. These approaches are not seen as mutually exclusive; insights from one approach can be used to reinforce assumptions in another. They can then also be further extended, such as with the inclusion of saturation by using MEC analysis. As was discussed, they ultimately are based upon Maxwell’s equations and in their attendant equations. They differ in how the quantities required are computed from the motor. End effects are not included in their analyses, and the cores are assumed to have very low reluctance. Once B is derived, quantities such as inductance, force, and torque can be relatively straightforward to compute with Lorentz force law, though the exact formulations vary.

Obtaining precise solutions to Maxwell’s equations for actual motor structures is difficult and can result in overtly complicated solutions. Simplifications due to geometric symmetry, coil/phase current distribution, and boundary conditions for the reluctance motors eases the solutions but ultimately affects the accuracy of the results to some extent. FEA results are an accepted reference point for determining the accuracy of proposed techniques’ solutions in published literature. Magnetostatic studies claim good agreement, as in the works previously mentioned in this paper.

7. Other Techniques Used in Reluctance Motor Research

There are other techniques used in mathematically analyzing the magnetic behavior of reluctance machines. The following techniques do not directly solve for B and are not strictly for analytically deriving it. However, they provide additional means of extending the work previously described to areas not allowed by their assumptions. These areas include realistic BH curves, magnetic saturation, and force/vibration production in the entire motor body. Thus, a complete review of magnetostatic derivations must consider the following methods.

7.1. Magnetic Equivalent Circuits

The air gap permeance method can be seen as a subset of the more general magnetic equivalent circuit (MEC) approach. Usually assuming a specific value of the relative permeability ${\mu }_r$, the magnetic flux, reluctance, and MMF are used analogously to the electric current, resistance, and voltage, respectively. By assigning those values to the different sections of the motor, a magnetic equivalent circuit can then be solved. An example of such a circuit for an SRM is shown in Figure 11.

It is commonly used in reluctance motor analysis and can also be used in conjunction with the previously described methods, especially when saturation and realistic ${\mu }_r$ is to be included in the analyses, as in Figure 12 [45]. Other possibilities are using fluxes solved from air gap permeance techniques [21] or conformal analysis [20] at various excitation currents and rotor positions and updating the values of the magnetic circuit. An example in Figure 13 shows magnetic potential usage with MEC [20]. It can even be applied to segmented motor designs [45,46] since the reluctances can be assigned freely. This freedom also enables choices regarding which sections to ignore or how to include leakage fluxes due to fringing [47]. It may also be used when nonlinearity or saturation is considered [48].

The circuit may be solved analytically but it is often used in an iterative manner until some chosen error criteria is passed. It generally does not provide a specific picture of how the magnetic field is distributed around the motor – in fact, the reverse approach is taken. The flux is assumed to flow in certain iron and air gap region first, upon which a reluctance element is applied to that path. This is unlike analytical methods where this distribution develops from the assumed boundary conditions and geometry of the machine. MEC does relatively easily accommodate changing permeability conditions around the machine.

7.2. Reluctance Mesh and Flux Tubes

The reluctance mesh method can be viewed as middle ground between MEC analysis and finite difference techniques. The iron structure and air gaps of the motor are discretized into mesh elements as exemplified in [49,50] and shown in Figure 14.

The potentials assigned at the nodes follow from the discretization of magnetic scalar potential in (22). A related element in this analysis is the flux tube [49]. This is an element in the discretized electric motor with arbitrary shapes as in Figure 15, from whose side no leakage flux escapes. It has also been employed in assuming flux flow in the air gaps at the aligned and unaligned positions [51], as in Figure 16. This is effectively a more flexible element whose shape can be assigned based on assumed shapes of the flux paths. This includes circular sectors and trapezoidal volumes [11].

The discretization and solution steps of this method can conceivably be performed manually, but it best lends itself to programmatic approach as in numerical techniques for differential equations. While it does not provide an expression or equation for the resulting magnetic field distribution in the motor, it is useful in modelling nuances in the magnetic field distribution similar to MEC but with more degrees of freedom due to the possible shapes of flux tubes.

7.3. Maxwell Stress Tensor

The Maxwell stress tensor (MST) is an expression for the strain tensor used primarily for computation of forces due to the electromagnetic field in a certain volume. There are a variety of equivalent forms of expressing the tensor $\mathit{\sigma }$, and one is shown in (34) in cartesian coordinates containing only a magnetostatic field [52].

$$\mathit{\sigma }=\left[\begin{array}{cc}{B}_x^2-{\displaystyle \frac{1}{2}}{\left|B\right|}^2& B_x{B}_y\\ B_y{B}_x& B_y^2-{\displaystyle \frac{1}{2}}{\left|B\right|}^2\end{array}\right].$$

(34)

This can be integrated along the surface S of volume V to find the torque T exerted by the magnetic field [53]:

$$\mathit{T}={\oint }_S\mathit{r}\times \mathit{\sigma }·\mathit{d}\mathit{S}.$$

(35)

MST can be quite difficult to solve by hand for physical systems and was not often employed before widespread use of computers [54]. The choice of the path or surface in which the tensor integral is evaluated is one area of interest in the context of reluctance motors, especially considering the slotting of reluctance machines [55,56].

As the name implies, it is most useful when computing for forces and torques present in reluctance motors [11,57]. It is not generally used to solve for the magnetic field itself; it is instead assumed that the Bfield in the region to be solved is already known. It is also mostly used in regions with ${\mu }_r=1$, and co-energy techniques might be better suited for regions with varying permeability. An advantage is that it is directly obtained from Maxwell’s Equations and deals with how the electromagnetic field itself exerts force. Thus, it can be used to solve forces, torques, or vibrations on a sounder theoretical basis once the magnetic field has been solved through other means.

8. Potential Applications to Segmented Reluctance Motor Design

The techniques discussed so far have been applied to reluctance motors with conventional structures; that is, they have the classic cylindrical stator and rotor with a continuous body made of essentially one piece of iron. Flux flows throughout the machine’s body and may enter the rotor through multiple poles as in Figure 17a [58].

However, given the fundamental nature of how switched reluctance motors operate, this conventional structure is not strictly necessary if the phases are properly commutated, much like a circularly arranged set of electromagnets consecutively attracting an iron rotor. For SRMs, it is possible to segment either or both the stator and rotor, with a nonmagnetic material securing them in place. Its advantages include smaller amount of electrical steel used [60], efficient flux path distribution, torque improvements [61], and even increase fault tolerance through the employment of modular motor designs [62]. With a segmented stator, less MMF might be needed because the stator back iron will have small reluctance (if finite ${\mu }_r$is assumed). Granted, segmented structures also introduce new mechanical issues to be surmounted such as ensuring the concentricity of the parts [63] and possible vibration due to imbalance if the rotor is segmented, but the advantages are enough to motivate further work. As an aside, segmentation can also be applied to motor drivers simultaneously because each SRM phase can be operated individually with no need of a common node as in wye or delta connected drivers. There is also the fact that torque direction is independent of current direction. Asymmetric bridge converters (ABCs) can then be used to drive them [16].

Research into analytic methods in solving the magnetic fields in such machines using the discussed analytical methods here are quite few. Magnetic equivalent circuits have been used in analysis of conventional stator, segmented rotor (CSSR) SRM designs [61]. Most existing work focuses on using the FEA-obtained rotor position-current-flux relationship of the machine together with the classic voltage equation ${\displaystyle \frac{d\lambda }{dt}}=V_t-Ri\left(t\right)$, from which the motor performance under different mechanical loading or motor drive scenarios can be modelled [64,65,66]. The usage of analytical techniques in segmented motors will have the same advantages as before such as less computer power reliance and better mathematical insight into the input-output relationships of the machine. However, the discontinuous nature of the machine must impose a different set of necessary boundary conditions before the previously discussed methods can be applied. This different set of boundary conditions, combined with the usual reliance on simulation techniques for advantageous segmented machines, constitutes a problem suitable for mathematical analysis magnetostatic field of segmented motor designs. The resulting techniques can result in further advancement of this motor topology.

The flux flow is constrained to the nearest stator poles even when only the rotor is segmented. This can be seen in Figure 17b, which shows a CSSR design [58]. The concentration of flux is even more pronounced in the aligned position for a segmented stator, segmented rotor (SSSR) structure shown in Figure 18. In both cases, the rotor segments act as shunts for the flux leaving the nearest stator poles so that the position of highest inductance for the phase involved is achieved.

This near-total isolation between the magnetic sections of the SSSR motor can both simplify and introduce complications in applying the air gap permeance techniques. The flux produced by a phase will only be restricted to the pole tips of the rotor or stator segments involved, especially for SSSR motors. It is observed that the face of the rotor pole is almost always within nominal air gap distance until the rotor is at the unaligned rotor position, as depicted in Figure 19. Thus, the resulting permeance function $\Lambda (\theta ,t)$ could be almost perfectly uniform so there are few changes to the air gap distance along the air gap periphery though this depends on the precise shape of the stator and/or rotor segments. The MMF distribution function may have a slightly different expression for its distribution depending on how the segmented stator is wound (especially if the phases are electrically isolated from each other), but it should mostly emerge out of the stator pole face close to rotor segments.

Fundamentally it does not appear necessary in air gap permeance analyses to have a continuous motor structure as long as the MMF distribution and air gap paths are properly considered. It does appear that the flux flows through all the segments when the rotor is not aligned in SSSR SRM as in Figure 18b. Given the uniformity of $\Lambda (\theta ,t)$ it may then appear that the torque variation comes from the variation in the MMF along the air gap, or that the effect of this slotting due to the separated segments may be enough to induce an observable torque. This hypothesis merits investigation in future studies.

Applying magnetic potential techniques may be tricky, as the first question would be how to assign the boundary conditions. Theoretically, there should be no difference with the assignment as in a conventional machine. As discussed in Section 3, the potential is zero at the walls of the slots where the flux does not emanate and $\pm NI$ at the pole faces. The problem of how to assign boundaries in the regions between the stator or rotor segments must be addressed; in a conventional machine the rotor slots are surrounded on at least three sides by iron and an air gap that could also have an MMF. In a segmented motor as in Figure 19, a rotor slot is bounded on two sides by iron, another by the air gap with the stator, while the last can be seen as an air gap or slot that extends through the rotor. This last air gap may be viewed as infinitely deep from the point of view a rotor slot magnetic scalar potential. It is possible to set the machine center to be at zero potential. The opposite is the case for the intervening regions between the stator segments, whose slots (which can be seen as infinitely deep due to the absence of a bottom iron) face outwards from the motor. To set realistic boundary conditions for these situations, insights might be gained from preliminary FEA magnetic potential results. The stator segments do have similarly shaped salient poles as in conventional structures (Figure 7 vs. Figure 18a) and so similar distortions as discussed in Section IIIB may be made to simplify computations. The rotor segments, on the other hand, usually have trapezoidal shapes [59,61,63], which complicate the resulting differential equation solutions. Solving the slot Laplace equations (22) and (23) for spaces between the segmented rotors then requires further examination. The supporting tabs of the segments could be ignored, since very little flux flows in them especially in the aligned rotor position, as in Figure 18a and 19.

Conformal transformation (particularly the SC transform), however, may be especially suited for this segmentation boundary condition problem because of the underlying linearity of magnetic potentials and Maxwell’s equations. Each stator segment, which is a source of MMF, can be individually analyzed. For example, the SC transformation could be applied to a stator segment with the proper f(z) into a rectangular domain similar to the process schematically shown in Figure 9. This may allow the derivation of the magnetic scalar potential field produced by that single segment, upon which the effect of the other segments that are actively excited can then be added due to linearity. A caveat is that Figure 9 shows the process of transformation for the motor slots, not the stator or rotor segment or poles themselves as in the situation just described. This slot-based transformation may therefore also be performed; in fact, the region between the stator segments can be imagined to be akin to a capacitor in that the walls of each segment (capacitor plates) are separated by a dielectric (air or aluminum), with the one segment actively with MMF (positive voltage) and other commutated off (ground). The potential field between two plates is an example problem commonly solved in conformal transformations, and could be applicable in a segmented stator context, thus avoiding complicated assumptions in assigning boundary conditions when basic magnetic potential analysis is applied. On the other hand, that scenario does not yet consider the effect of the nearby rotor segments which will provide a flux path with less reluctance. In fact, it can already be seen in Figure 18b how the flux flows through the neighboring stator and rotor segments.

Overall, the application of the aforementioned techniques in the mathematical analysis of the magnetostatic field in segmented motor design is relatively unexplored, given that segmented motor analysis in literature often involve FEA, the usage of reluctance motor properties as exemplified by the classic flux linkage-voltage relation ${\displaystyle \frac{d\lambda }{dt}}=V_t-Ri\left(t\right)$ and (4), or interpolation techniques. While approximate techniques for this machine topologies already exist [66], the hope is to provide a systematic, if not general, expression of inductance and torque production for segmented machines that are based on the fundamental distribution of the magnetic field.

9. Conclusions

Analytical methods for solving the magnetic field in a reluctance motor are derived from Maxwell’s equations, and there are various methods of making the problem tractable for the geometry of such a machine. Previous works focus on the application of these methods to various motor models instead of the analytical methods themselves. Mathematical models may be more desirable given the computational needs of numerical or simulation-based modelling.

In this review, the major techniques that rely on obtaining the magnetostatic field distribution using fundamental principles are discussed and analyzed. Other methods that also solve for the B field but are iterative are also shown. There are a multitude of possible ways to heuristically approach the motor magnetostatic field distribution employing different simplifications. However, the major techniques can be boiled down to using air gap permeance, magnetic potentials, and conformal transformations and these appear to be able to support each other in reluctance motor research as well. Their mathematical details and assumptions are indeed examined here, and it is shown that it is not impossible to link them together, making them more accurate by strengthening their underlying assumptions using ideas from other techniques. They are all able to reasonably approximate simulated results, however it appears that for saturated iron situations, magnetic equivalent circuit analysis is still the most relevant process.

As segmented stator and rotor designs have also emerged, possible applications of those analysis techniques in segmented motor designs are investigated. Avenues of approaching the analysis of such machines using the discussed derivations are thus conjectured and suggested as areas of future research. Air gap permeance must be set with different considerations with segmented rotor designs because of the near-constant air gap lengths between the stator-rotor pole faces. It is believed that conformal mapping is well suited for those with segmented stator designs, while magnetic potentials may allow for the solution of the magnetic field between the stator and rotor segments which are only bound by iron on two sides. This machine topology has properties that are not present in conventional structures and whose systematic study must then pose an interesting area of mathematical analysis.