Analytical Modeling of Magnetic Field Distribution at No Load for Surface Mounted Permanent Magnet Machines

: This paper presents an analytical study of the air-gap magnetic ﬁeld of a Surface Permanent Magnet (SPM) linear machine under no load. By means of the method of images, the complex expression of the magnetic ﬁeld of a conductor, inside the air gap between two smooth iron surfaces, is retrieved. Then, integrating the conductor expression, the formulation of the magnetic ﬁeld of a current sheet and thus the one of a SPM, using two vertical current sheets, is obtained. At last, the no-load magnetic ﬁeld expression, for a generic time instant, of a slotless machine is retrieved. The novelty of the proposed approach is the availability, due to a different calculation approach, of a unique closed-form formulation for the slotless machine air-gap ﬁeld, a quantity that, in literature, is usually present in Fourier series formulation. Additionally, as a means to calculate integral quantities and show the goodness of the method a complex slotting function is introduced to account for the slotted geometry. Finally, starting from Lorenz’s force formulation, the expression of the Maxwell tensor in complex form is retrieved and the contribution of forces, integral of the complex stress tensor quantity, will be calculated and compared with FEM simulations, showing a good agreement also with the analytical slotted model.


Introduction
The study of analytical techniques to predict the magnetic field distribution inside the air gap of the electrical machines has been, and still is a topic of interest, also after the improvement experienced by numerical methods.It is so because analytical methods not only help us to better understand the core parameters which influence the physical phenomena; but also because they are generally faster than their counterpart.Thus, even if under simplifying hypotheses, they are helpful, especially during the early part of a design project to predict motor performances or to carry out optimizations.Through the study of the air-gap magnetic field, it is possible to calculate quantities like flux linkages, back electromotive forces, and magnetic forces, hence, investigating the cause of noise and vibration inside electrical machines.
Among the analytical methods concerning the flux density distributions inside the air-gap of the electrical machines present in literature, it is found the method based on the solution of the analytical 2D Laplacian quasi-Poissonian field equations assuming infinitely permeable iron [1,2].Those types of models were introduced for slotless SPM motors in [3] and later extended to slotted geometries, throughout the use of the relative air-gap permeance [4], for the normal flux density component only, which, expressed in term of Fourier series and multiplied to the slotless flux density distribution gives the flux density profile of the slotted machine.Nevertheless, the lack of tangential flux density component of the slotted geometry, developed in [4] does not allow to retrieve quantities like the cogging torque, by integrating the tangential component of the Maxwell stress tensor.To this purpose, [5] introduced the complex permeance function allowing us to also retrieve the tangential component of the flux density air-gap field distribution and analytically calculate the torque, using the Maxwell stress tensor [6]; this process, unfortunately, requires the solution of a nonlinear complex equation, leading to higher computational time.The methods just mentioned together with the work done by Carter [7] and the Schwarz-Christoffel mapping method [8][9][10] fall all into the conformal mapping category.
Another category of models concerns the subdomain models, which solve the flux density field inside the slotted machine, by dividing the entire domain into subdomains, enclosing air, magnets, slots, and slot openings; solving Maxwell's equations into it using the separation of variables [11][12][13].Subdomains models can accurately predict flux density into magnetic air gap, magnets, and slot opening with very good accuracy, and they appear to be better also for the cogging torque calculation compared with permeance models [14].
Methods that rely on the equivalent magnetic circuit methods (EMC) can be also found and consist in drawing a magnetic circuit model of the machine with lumped or distributed reluctances parameters to model the machine's magnetic flux tubes and retrieve quantities like flux linkages and average flux density values [15,16].A variety of hybrids methods of the previously mentioned ones are also present [17][18][19].
Lastly, a less used and known category of analytical methods based on Lord's Kelvin method of images as [20,21] can be found and the presented paper belongs to it.The method here introduced models the slotless field of an SPM machine, not solving the 2D Laplacian quasi-Poissonian field equations like Zhu did in [3], but using a different approach based on the method of images and duality between the electric and magnetic fields.
In Sections 2.1 and 2.2 the flux density air-gap expression produced by the powered conductor stated by [22] and then formulated in complex form [23], has been extended to an indefinitely long current sheet energized by a linear current density: then, two vertical parallel current sheets have been used, in order to model a surface permanent magnet.
In Section 2.3 the totality of permanent magnets present on the machine periphery has been modeled introducing also the time variable.
Then, in Section 3.1, to account for the slotting effect, similarly to [5], a more compact formula of the relative air-gap permeance function has been introduced, starting from [24].
In Section 3.2, the flux density profile of the slotted machine, product of the complex slotting function, and the slotless flux density vector is retrieved and compared with FEM simulations, showing an outstanding agreement.
In Section 4.1, starting from the Lorentz's force formula, the Maxwell stress tensor in complex form has been obtained and in Section 4.2 it is used to retrieve both the normal and tangential component of the total force of the slotted machine.
Finally in Section 4.3 the computational time required for the proposed method is shown in comparison with FEM simulations.
What distinguishes the proposed model from the other ones is that: the slotless formulation is not expressed by means of a solution in the Fourier series but by means of a closed-form analytical formulation.Furthermore, the high level of accuracy of the flux density normal and tangential air-gap components calculation is in the same range as the sub-domain models.
Among the advantages of the above-mentioned method, there is the simplicity of implementation and fast response, bringing accurate results also for the Maxwell stress tensor calculation, which is not an easy task as it is well known in the literature.On the other side, the method shows the weak points of assuming unity relative recoil permeability and neglect saturation.However, modern SPM machines are usually equipped with PM materials with pu recoil permeability typically around 1.09 [25][26][27], so quite close to vacuum permeability.Moreover, it is possible to show that the actual air-gap flux density distribution taking into account saturation can be obtained from the unsaturated one, multiplied by suitably calculated saturation correcting factors.

Magnetic Field between Smooth Ferromagnetic Surfaces 2.1. Linear Conductor with Infinitesimal Section
The magnetic field distribution produced by an indefinite current-carrying conductor, cutting the paper plane perpendicularly and placed into the air gap between two smooth ideal ferromagnetic surfaces, can be modeled analytically in Cartesian coordinates, using the method of images described in [22].The aforementioned method makes use of the following assumptions and limitations: i Stator and moving surfaces are supposed to be smooth and separated by an air gap of uniform width.ii Magnetic field is just 2D and lies in the paper plane.iii The machine is supposed to be so long into the paper plane direction, so end effects are negligible.iv Relative recoil permeability is considered to be equal to the vacuum permeability Iron permeability is assumed to be infinite so the method of superimposition can be used.vi The iron parts are perfectly laminated thus no eddy current can be induced in the iron by the varying magnetic field.vii The conductors placed inside the air gap extend indefinitely perpendicular to the paper plane.viii Skin effect due to alternating currents flowing in the conductors is neglected.
By representing the magnetic field vector quantities as complex numbers, made by a real and imaginary part, representative respectively of the tangential and normal vector components, let's focus on Figure 1, where a conductor powered by a current I, flowing perpendicularly outwards the paper plane xy, is positioned into the air-gap complex plane at a generic point p of coordinates z p = x p + jy p , between two smooth infinitely permeable iron surfaces.To be able to calculate the field strength vector at a certain point z between two smooth iron surfaces, the aforementioned configuration can be exchanged with an equivalent one in which the iron parts are removed and substituted with two infinite series of images.The conductors from the first and second series of images are placed at positions z p + j2ng and zp + j2ng, with n ∈ Z ranging from minus and plus infinity.By using the theorem of superimposition the field strength inside the air-gap can be expressed as the sum of two mathematical series representative of the two image series, a formulation a little complicated which has been simplified in [23] giving rise to the formulation: where I is the current flowing inside the conductor, g the air-gap length, z p and zp the position of the conductor in complex coordinates and its conjugate.

Permanent Magnet Current Sheet Model
As it is well known already in the literature, the magnetic field produced by a permanent magnet can be modeled by the use of current sheets.In the following subsection, the approach which models the magnetic field produced by a conductor set between two smooth ideal iron surfaces previously introduced will be extended to current sheets in order to model the magnetic field produced by surface permanent magnets.
In Figure 2 current sheet indefinitely extended into the paper plane and vertically oriented, inside the air gap of two smooth iron surfaces are considered, distributed along the points p, with coordinates z p = x p + jy p .This distribution has extreme points z p1 = x p + jy p1 and z p2 = x p + jy p2 .Let's assume that the current sheet is powered by a constant current density of ∆ = I/h cs , flowing outwards the paper plane, where I is the total current and h cs the current sheet finite length.By taking into consideration an infinitesimal current sheet element ∆ • dy p , it is possible to express the infinitesimal flux density air-gap complex vector as: Then, let's integrate (2) between y p1 and y p2 , with y p1 set to zero (the current sheet starts from the lower iron surface), the air-gap flux density vector, produced by a current sheet of length h cs , with z p1 = x p + j0 and z p2 = x p + jy p2 its ends points and powered by a current density ∆ flowing outwards the paper plane is retrieved: Making use of (3), the field of an SPM, leaning on the lower iron surface and with the north pole oriented positively towards y, can be modeled by the use of two current sheets vertically oriented, whose current density equals the PM coercivity The first current sheet will be positioned on the right, inside the first quadrant, powered by a current density ∆ = −H o , flowing inwards the xy plane, whose ends are situated at z p1 = b m /2 + j0 and z p2 = b m /2 + jh m , being b m and h m the width and height of the PM.The second one will instead be positioned symmetrically to the y axis on the left, inside the second quadrant, powered by a current density ∆ = H o which flows outwards the plane and whose ends point equals the negative conjugates of points z p1 and z p2 , giving: wherein the real part B PM x (z) = Re{B PM (z)} and imaginary part B PM y (z) = Im{B PM (z)} represent the tangential and normal components of the vector flux density function of a generic complex point z, produced by an SPM leaning on the lower iron surface, with y as an axis of symmetry and north oriented positively towards it.
The theoretical case of a SPM, between two smooth iron surfaces, indefinitely extended along x, that lies on top of the lower iron surface, is now investigated making use of the FEM analysis and the analytical vertical current sheets model previously introduced.The magnet, positioned in the center of the coordinate system, has peripheral size b m = 102.05mm and height h m = 10 mm, and its north points toward y positive.In Figure 3a,b the analytical and FEM normal and tangential flux density distribution components on the middle of the mechanical air gap at y M = g − δ/2 over the permanent magnet peripheral extension is plotted.FEM simulations have been carried out with infinitely permeable iron and relative recoil permeability set to one.Both the normal and tangential distributions are expressed in p.u. with respect to B i = µ o • U/g, which is the ideal flux density flat profile reference value, with U the air-gap magnetic voltage drop and g the distance between the ferromagnetic surfaces.is important to mention that, for the considered ideal case, where the magnet is between two smooth iron surfaces extended indefinitely along x, the magnetic voltage drop in the part of air-gap where the flux comes back, equals zero.Being the crossing area of the returning flux of infinite extension, as it is clearly shown in Figure 3a, everywhere the magnetic voltage drop is positive, thus

Magnetic Field due to a Sequence of Permanent Magnets
The SPM current sheet model introduced in the previous subsection can be extended to a series of permanent magnets, to model the magnetic field produced by an SPM machine during no load working conditions.For the sake of simplicity, the method will be now applied to a linear SPM machine, but it can also be generalized for cylindrical geometries using conformal mapping, as it is done for a series of concentrated conductors in [23].It is worth mentioning that the following method works also for more complex magnet shapes for example the Halbach ones.The machine type examined for the purpose is shown in Figure 5a whose main data are listed in Table 1; the machine sizes listed below are retrieved using a pre-design procedure applied to an axial flux PM machine, linearized along the average radius [28].Let's consider first the slotless machine pictured in Figure 5a, where A and B are both filled with iron.Considering as moving the PM slider at a constant speed equal to 2τ m /T, (4) can be also expressed as a function of time and generalized for a k th magnet belonging to the machine periphery as: where the term ∆ kx represents the shift in space, function of time, along the negative coordinate x of a generic k th SPM belonging to the machine periphery.Being in magnetic linearity, from (5), the contribution of the flux density at the point z into the air-gap, produced by the SPMs totality situated on the moving structure of the considered configuration can be finally stated as: in which the term (−1) k justifies the presence of the alternating north and south adjacent polarities.Figure 6a,b shows the flux density distributions components of the total SPM model stated in (6) on the middle air-gap line in comparison with a magneto-static FEM simulation with infinitely permeable iron for a PM slider position corresponding to a time instant t = T/15 (T = 16.67 ms).Note that, the inequality between adjacent peaks for both the FEM and analytical patterns in Figure 6b isn't real but just a plotting effect due to the spatial discretization of the machine periphery.

Slot Opening Model
In the subsequent sections, the presence of the slot opening and its influence on the magnetic field distribution inside the air gap will be analyzed.Using the approach based on the solution in Fourier series of the Laplace equation [24], a unique and more compact expression of the complex relative air-gap permeance function is obtained, effective for any segment positioned inside the air-gap of the machine.
Figure 7 exhibits the flux lines map, produced by a slot with two rectangular energized active sides with a total slot m.m.f of constant value M s , inside a linear ferromagnetic slotted structure, for a slot opening size b so = 6 mm.The focus will be placed according to the classical theory on the area regarding the unenergized slot, the slot where the m.m.f contribute is null, being the energized slot area object of another work.
Considering a generic point z = x + jy inside the air-gap domain, it is possible to define a complex function β(z) = β x (z) + jβ y (z) whose components are: ) with B i = µ o M s /2g.Recalling the flux density vector B(z) between smooth iron surfaces the slotted flux density vector can be stated as: where λ(z) is the relative complex permanence function [5].Due to the fact that the slotless flux density vector has only the vertical component (B(z) = jB i ) from ( 9): Hence it is possible to define a complex slotting function β(z) which transforms the slotless flux density vector B(z) to a slotted vector according to: and finally, starting from the approach explained in [24] and through some manipulations, the complex slotting function can be formulated as: This slotting complex function is periodic along the coordinate x over a machine tooth pitch τ t and in respect to [24] it has the advantage to be just one compact complex formula generalized for a generic complex point z = x + jy within the air-gap domain.
Now bringing the attention to the schematic shown below in Figure 8, where a generic exploration segment γ e is horizontally positioned at y e in the range 0 ≤ y e ≤ g, the slotted flux density vector along segment γ e can be obtained as follows:

Magnetic Field due to a Sequence of Permanent Magnets for Slotted Surfaces
It will be now taken into consideration the slotted machine in Figure 5a where A is filled with iron and B is filled with air.In Section 3.1 behavior of the flux density vector inside the mechanical air-gap, in proximity of an unpowered slot has been investigated, defining a complex slotting function analogous to the complex permeance function.As soon as the value of the exploration coordinate is chosen, the slotted flux density vector on the horizontal line at y = y e as a function of time can be retrieved: with ( 12) also taken along a horizontal line at y = y e shifted by τ t /2 to align with the coordinate system of Figure 5a.It can be easily noticed that, being the slot opening of the linear machine positioned on the stator, the complex slotting function is fixed and does not depend on time.
Figures 10a and 11b show the normal and tangential flux density distributions on the middle air-gap line z = x + jy M for a time instant t = T/15, where t = 0 corresponds to the time in which the north magnet (white color PM in Figure 5a is aligned with the central tooth n th = 5.The analytical pattern is perfectly superimposed almost everywhere to the FEM one, except in some areas where the two patterns slightly differ from each other, as it can be seen in Figures 10b and 11b.The reason is explained by the fact that, in the area of interest, a magnet edge lies in front of the slot opening, slightly modifying, because of fringing, the flux lines patterns from the ideal ones represented by the slot opening complex function.Notice that the method here presented models the no-load air-gap flux density pattern of a machine with semi-closed slots but it can be applied to model the air-gap flux density of machines with open slots as well.

Complex Vectors and Mechanical Actions
In the antecedent sections, the analytical expression of the magnetic field inside the air gap of a smooth infinitely permeable iron machine has been introduced and subsequently, by the use of ( 12), the complex slotting function was obtained to account for the presence of the slots.Now being aware of the analytical field distribution, into the air-gap domain, it is possible to determine integral quantities like magnetic fluxes and mechanical forces.Here attention will be devoted only to the mechanical forces generated inside the machine.
Let's consider the infinitesimal portion of a conductor, with infinitesimal cross-section dA and length d , immersed in a magnetic field and powered by a current density J shown in Figure 12.Writing the well-known Lorenz's expression, the infinitesimal contribution of force acting on the conductor is: Moving the vector information from the vector infinitesimal element d to the current density J and dividing by the infinitesimal volume dv = dAd allows to write: which is the expression of the force per unit volume, that, integrated over the conductor volume, gives the entire mechanical force acting on it.Consequently ( 16) can be integrated in the volume of interest: Because the field is planar, lying entirely on the xy plane, (17) can be written as a surface integral over the area of interest A multiplied by the axial length : At this stage, before to proceed with further manipulations of the aforementioned equations, a small mathematical statement must be postulated: Considering a generic complex vector W(z(x, y)) = W x + jW y , lying on the xy plane and z = x − jy, being the conjugate of a generic complex point z, being ∂ z ∂x = 1 and ∂ z ∂y = −j it is possible to write: and from (19a): Introducing another generic complex vector Ω(z(x, y)) = Ω x + jΩ y , function of the complex variable z it can be demonstrated that: thus looking at (19b) and (19c), setting Ω = ∇, the curl and divergence of the generic complex vector W can be expressed as follows: Then, according to what is stated above in (19d) and (19e), the divergence and curl of the vector H can be expressed as the real and imaginary part of twice the derivative of H in respect to z: where k is the versor of the axis perpendicular to the complex plane.However, due to the solenoidal property of the field and being inside a constant permeability medium (air), (20) equals zero and the derivative of H in respect to z is purely imaginary: consequently ( 16) can be stated as: Finally by the use of ( 23), remembering the property explained in (19a)-(19e) the total force acting on the conductor expressed into (18) can be stated as a force per unit length F = F/ as follows: If the generic closed curve C oriented counterclockwise shown in Figure 13 is taken, by applying the divergence and curl theorem: wherein the first term on the first line of ( 25) is the flux of B 2 /(2µ o ) across C and the second term is the circulation of B 2 /(2µ o ) along C respectively.Then the final expression of the complex force vector per unit length is obtained.The quantity B 2 /(2µ o ) is the conjugate of the flux density vector squared, divided by twice the permeability of vacuum.Finally deriving (25) in respect to dp the Maxwell stress complex tensor expression can be retrieved: It is worth mentioning that the equations described in this subsection are of general validity, also with currents flowing in the stator windings.

Maxwell Stress Complex Tensor and Forces
The Maxwell stress tensor is used in classical electromagnetism to calculate forces; the volume integral of J × B, throughout some mathematical manipulations is reduced to a surface integral of the stress tensor valid for any enclosing surface in the air of the considered volume, simplifying the calculations.If, into (26) the slotted flux density vector stated in ( 14) is inserted, it is possible to write: obtaining the Maxwell tensor complex expression, the function of time for the slotted configuration.Figure 14a wherein the first term on the first line of ( 25) is the flux of B 2 /(2µ o ) across C and the second term is the circulation of B 2 /(2µ o ) along C respectively.Then the final expression of the complex force vector per unit length is obtained.The quantity B 2 /(2µ o ) is the conjugate of the flux density vector squared, divided by twice the permeability of vacuum.Finally deriving (25) in respect to dp the Maxwell stress complex tensor expression can be retrieved: It is worth mentioning that the equations described in this subsection are of general validity, also with currents flowing in the stator windings.

Maxwell Stress Complex Tensor and Forces
The Maxwell stress tensor is used in classical electromagnetism to calculate forces; the volume integral of J × B, throughout some mathematical manipulations is reduced to a surface integral of the stress tensor valid for any enclosing surface in the air of the considered volume, simplifying the calculations.If, into (26) the slotted flux density vector stated in ( 14) is inserted, it is possible to write: obtaining the Maxwell tensor complex expression, the function of time for the slotted configuration.Figure 14a  As concerns the calculation of the forces per unit length, there are two different targets: • evaluation of the forces applied to each tooth, typically considered in NVH calculations; • evaluation of the forces applied to the global system.
With reference to the linear system of Figure 5a, the force on the tooth n th can be expressed as: Equation ( 28) is integrated on the position y = y M , on the middle mechanical air-gap line, over the tooth pitch of the n th tooth, centered on the tooth axis x n th value, with n th = 1, 2, ..., N th .The real part and imaginary part of (28) represent the tangential and normal components of the force per unit length acting on the n th tooth tip.The force applied to the global stator can be written as: but it can also be expressed in the following manner: Equations ( 27)-( 30) have general validity, also in case of flowing currents in the stator windings.If (28) is applied to the central tooth n th = 5 under no load conditions, the normal and tangential forces per unit length are the ones shown in Figure 15.
Energies 2023, 1, 0 14 of 19 As concerns the calculation of the forces per unit length, there are two different targets: • evaluation of the forces applied to each tooth, typically considered in NVH calculations; • evaluation of the forces applied to the global system.
With reference to the linear system of Figure 5a, the force on the tooth n th can be expressed as: Equation ( 28) is integrated on the position y = y M , on the middle mechanical air-gap line, over the tooth pitch of the n th tooth, centered on the tooth axis x n th value, with n th = 1, 2, ..., N th .The real part and imaginary part of (28) represent the tangential and normal components of the force per unit length acting on the n th tooth tip.The force applied to the global stator can be written as: but it can also be expressed in the following manner: Equations ( 27)-( 30) have general validity, also in case of flowing currents in the stator windings.If (28) is applied to the central tooth n th = 5 under no load conditions, the normal and tangential forces per unit length are the ones shown in Figure 15.The following remarks can be highlighted: • FEM and analytical waveforms are in fair agreement; • the normal component amplitude is predominant with respect to the tangential one; • the force waveform acting on a generic tooth tip has a frequency that is double the fundamental frequency (1/T).
It can be realized that the no load force waveforms are equal for each tooth tip but shifted in time by a quantity multiple of T • τ t /(2τ m ).It is possible then to calculate the total force acting on the stator structure from (28) referred to a single machine tooth tip; for convenience the tooth whose axis is aligned with the y axis of the reference frame n th = 5: The following remarks can be highlighted: • FEM and analytical waveforms are in fair agreement; • the normal component amplitude is predominant with respect to the tangential one; • the force waveform acting on a generic tooth tip has a frequency that is double the fundamental frequency (1/T).
It can be realized that the no load force waveforms are equal for each tooth tip but shifted in time by a quantity multiple of T • τ t /(2τ m ).It is possible then to calculate the total force acting on the stator structure from (28) referred to a single machine tooth tip; for convenience the tooth whose axis is aligned with the y axis of the reference frame n th = 5:  The following considerations can be done: • the diagrams shown in Figure 16a,b in which both analytical and FEM calculations are based on Maxwell tensor, can be also obtained using (29); • (29) implies a greater computation time in respect to (31) due to a wider spacial integration interval; • on this peripheral extension, corresponding to the space period, FEM calculation can be also performed by using the Virtual Work Principle; • the analytically calculated normal component shows a better agreement with respect to the FEM result; • a wider difference appears for the tangential component: the reason for that could be a massive difference in the order of magnitude of these two components, making the numerical noise error predominant in the analytical calculation of the minor quantity; moreover, the tangent component calculation is affected by local inaccuracies due to fringing effect as shown in Figure 11b.
Figure 17 shows the FEM model at t = 0 which has been used in order to validate the introduced analytical model together with a picture of the mesh inside the equivalent air-gap region.The two dotted lines B.C. indicate the boundary conditions that allow us to obtain periodicity in the field patterns.With the aim of carrying out an equal comparison of the introduced model with the numerical one, the relative recoil permeability of the PMs of the FEM model was set to one (µ r = 1) and both the rotor and iron relative permeabilities were set to a very high numerical value (µ f e = 10 6 ) to recreate, numerically speaking, the theoretical condition of infinite iron permeability and neglect the saturation effect.Magnetic saturation is not included in the model described in this paper, but already developed procedures demonstrated that the field in the air gap under saturated conditions can be obtained starting from the unsaturated distribution using suited correcting coefficients.The following considerations can be done: • the diagrams shown in Figure 16a,b in which both analytical and FEM calculations are based on Maxwell tensor, can be also obtained using (29); • (29) implies a greater computation time in respect to (31) due to a wider spacial integration interval; • on this peripheral extension, corresponding to the space period, FEM calculation can be also performed by using the Virtual Work Principle; • the analytically calculated normal component shows a better agreement with respect to the FEM result; • a wider difference appears for the tangential component: the reason for that could be a massive difference in the order of magnitude of these two components, making the numerical noise error predominant in the analytical calculation of the minor quantity; moreover, the tangent component calculation is affected by local inaccuracies due to fringing effect as shown in Figure 11b.
Figure 17 shows the FEM model at t = 0 which has been used in order to validate the introduced analytical model together with a picture of the mesh inside the equivalent air-gap region.The two dotted lines B.C. indicate the boundary conditions that allow us to obtain periodicity in the field patterns.With the aim of carrying out an equal comparison of the introduced model with the numerical one, the relative recoil permeability of the PMs of the FEM model was set to one (µ r = 1) and both the rotor and iron relative permeabilities were set to a very high numerical value (µ f e = 10 6 ) to recreate, numerically speaking, the theoretical condition of infinite iron permeability and neglect the saturation effect.Magnetic saturation is not included in the model described in this paper, but already developed procedures demonstrated that the field in the air gap under saturated conditions can be obtained starting from the unsaturated distribution using suited correcting coefficients.

Computational Efficiency
The time required by the analytical and FEM for the calculation of forces is now compared.Table 2 shows the time necessary to compute (29) and (31) with respect to transient FEM simulations using the Virtual Work Principle (VWP) and the Maxwell tensor (MT).The calculator is an i7-10850H @ 2.70GHz 32GB RAM, the analytical model has been implemented using MATLAB R2019b while the FEM simulations have been carried out with Ansys Electronics Desktop 2021 R2.In the FEM environment, the entire portion of the machine periphery has to be modeled, due to the fact that τ t and τ m are not equal no smaller equivalent models exploiting spatial periodicity can be used.The FEM MT requires a much higher computational time with respect to the FEM VWP because, for each time step of the simulation, the field map must be saved.Both the analytical methods are faster than FEM and comparing slowest analytical (ANA1) with the fastest FEM (FEM VWP), analytical is more than 100 times faster.This makes the presented method a good choice for an easily implemented and fast analytical calculation of the forces acting on the machine under no load conditions, opening also the possibility to develop future models for the study of the armature reaction and the on load condition.

Conclusions
The presented paper introduces a new method to analyze the flux density generated into the air gap of an SPM machine under no-load conditions.Starting from the model of a single conductor between two smooth iron surfaces, the analytical formulation of two vertical parallel current sheets have been retrieved and used to model an SPM, showing a very good agreement with magneto-static FEM simulations, for the slotless machine but also, after the application of the slotting function, for the slotted machine.
Subsequently, the parameter time has been introduced into the model; and the Maxwell tensor in the complex formulation has been retrieved starting from Lorentz's force equation, adding the possibility to calculate the motor total forces, by the integration of the complex tensor along the middle air-gap machine periphery using two different approaches.
The proposed method is applied to a classical SPM configuration although it works very well for more complex PM configurations too, like the Halbach one, as it will be shown in future work.Future research will focus on the slotting function in front of the energized slots and will aim to deal with saturation and leakage which is one of the main limitations of the proposed model, especially under loaded conditions.By using conformal mapping, the method will be adapted to cylindrical machines.Moreover, saturated conditions will be modeled by suited correction factors that, applied to unsaturated air-gap flux density distributions, give the saturated ones.z p , zp generic point p of complex coordinates x p + jy p where the conductor is positioned and its conjugate.z, z generic point of coordinate x + jy inside the air-gap complex plane and its conjugate.β(z), λ(z) slotting, permeance complex function.
γ e , γ M generic exploration segment at y = y e , middle mechanical air-gap segment at y = y M inside the air-gap.ξ, ρ so p.u. ratios x/τ t and b so /τ t .τ m , τ t PM pitch, tooth pitch.

Figure 1 .
Figure 1.Conductor powered by a current I, flowing outwards the paper plane, between two smooth infinitely permeable iron surfaces, in the complex plane.

Figure 2 .
Figure 2. Vertically oriented current sheet of length h cs , powered by a current density ∆ flowing outwards the paper plane, whose extremes are situated at z p1 and z p2 between two smooth iron surfaces (g = 11.5 mm and δ = 1.5 mm).

Figure 3 .
Figure 3. Analytical and FEM flux density distributions, on the middle mechanical air-gap line, caused by an SPM between two smooth iron surfaces (−b m ≤ x ≤ b m N x = 1000).(a) p.u. normal; (b) tangential component.

Figure 4 .
Figure 4. Analytical and FEM flux density distributions on the middle PM line, caused by an SPM between two smooth iron surfaces (−b m ≤ x ≤ b m N x = 1000).(a) p.u. normal; (b) tangential component.

Figure 5 .
Figure 5. (a) Machine case of study: (1) slotless: zones A and B both filled with iron; (2) slotted: zone A filled with iron and B with air.(b) Zoom with main sizes.

Figure 6 .
Figure 6.Analytical and FEM no-load flux density distributions on the middle mechanical air-gap line of the slotless machine of Figure 5a, with PM slider position corresponding to a time instant t = T/15 (−4τ m ≤ x ≤ 4τ m N x = 1000).(a) normal; (b) tangential component.

Figure 7 .
Figure 7. Flux lines inside energized and not energized slot (slot opening b so = 6 mm).Even if not shown in the figure, the magnetic structure is symmetrical with respect to the energized slot axis.

Figure
Figure 9a,b show the p.u. flux density distributions retrieved with (12) compared with the FEM on an horizontal segment at y = y M in front of a slot opening b so = 6 mm.

Figure 10 .
Figure 10.Analytical and FEM no-load flux density normal component on the middle mechanical air-gap line of the slotted machine of Figure 5a, with PM slider position corresponding to a time instant t = T/15 (−4τ m ≤ x 4τ m N x = 1000, N h = 150); (b,c) are two different selective zooms of (a).

Figure 11 .
Figure 11.Analytical and FEM no-load flux density tangential component on the middle mechanical air-gap line of the slotted machine of Figure 5a, with PM slider position corresponding to a time instant t = T/15 (−4τ m ≤ x ≤ 4τ m N x = 1000, N h = 150); (b,c) are two different selective zooms of (a).

Figure 12 .
Figure 12.Infinitesimal conductor portion powered by a current density J immersed in a magnetic field of flux density value B.

Figure 13 .
Figure 13.Generic closed curve C inside the complex plane embracing an area A.

Figure 13 .
Figure 13.Generic closed curve C inside the complex plane embracing an area A.

Figure 14 .
Figure 14.Analytical and FEM no load Maxwell tensor distributions on the middle mechanical air-gap line of the slotted machine of Figure 5a, with PM slider position corresponding to a time instant t = T/15 (−4τ m ≤ x ≤ 4τ m N x = 1000, N h = 150).(a) normal; (b) tangential component.

Figure 14 .
Figure 14.Analytical and FEM no load Maxwell tensor distributions on the middle mechanical air-gap line of the slotted machine of Figure 5a, with PM slider position corresponding to a time instant t = T/15 (−4τ m ≤ x ≤ 4τ m N x = 1000, N h = 150).(a) normal; (b) tangential component.

Figure 16 .
Figure 16.Analytical and FEM no load force per unit length of the slotted machine of Figure 5a (0 ≤ t ≤ T/18 N t = 120).(a) normal component; (b) tangential component.The analytical calculation was based on (31) and FEM curve was calculated using the Maxwell tensor.

Figure 16 .
Figure 16.Analytical and FEM no load force per unit length of the slotted machine of Figure 5a (0 ≤ t ≤ T/18 N t = 120).(a) normal component; (b) tangential component.The analytical calculation was based on (31) and FEM curve was calculated using the Maxwell tensor.

Figure 17 .
Figure 17.FEM model of the machine case of study and selective zoom of the mesh inside the equivalent air-gap region.

Table 1 .
Machine Main Data.
31)wherein the real part F o x (t) = Re{F o (t)} and imaginary part F o y (t) = Im{F o (t)} represent the tangential and normal components of the force per unit length acting on the stator structure, under no load conditions, at a generic time instant t.Computing (31) for N t samples of time with 0 ≤ t ≤ T, allows us to retrieve a transient analytical calculation of the machine's instantaneous force acting on the stator.Figure16a,bshow the normal and tangential component of (31) over a period of the cogging force, which corresponds to 1/18 of the electrical period T.wherein the real part F o x (t) = Re{F o (t)} and imaginary part F o y (t) = Im{F o (t)} represent the tangential and normal components of the force per unit length acting on the stator structure, under no load conditions, at a generic time instant t.Computing (31) for N t samples of time with 0 ≤ t ≤ T, allows us to retrieve a transient analytical calculation of the machine's instantaneous force acting on the stator.Figure16a,bshow the normal and tangential component of (31) over a period of the cogging force, which corresponds to 1/18 of the electrical period T.