This article proposes a mathematical model of an axial flux induction motor (AFIM) with one stator and one rotor. The model is based on the expression for the electromagnetic torque, which presents a function of two independent variables: the axial length of the stator core and the flux density in the air gap. This allows calculating the main dimensions of the motor with the highest possible torque density. Thus, developed model is suitable for designing the motor of specified volume with maximum torque, and solving the inverse problem of minimizing the machine volume with the specified torque. The detailed output of the model and the results of the calculations for the low-power engine powered by voltage of 7.35 V (RMS) are given. The results are validated using FEM in ANSYS software: with the outer motor diameter of 0.11 m, the flux density in it reaches 1.2 T.

expression of electromagnetic torque design optimization axial flux induction motor (AFIM) electromagnetic torque density analytical design electric vehicles rotating electrical machines
1. Introduction

The traditional system of designing electrical machines is based on the analysis of a previous experience of their creation, as shown in . With this approach, new samples inherit the advantages and disadvantages of the previous models. Electric machines of the traditional design have a big statistical base accumulated over the previous century of design, allowing for the design of efficient machines of any power. However, when the task is the design of an unusual electrical machine appears, this approach is not possible because there are no statistics available. There are various methods for designing such machines to produce the most effective results.

There are different indicators for evaluating the design quality of an electrical machine: power density , the volume torque density and weight torque density , the volume and weight of permanent magnets per Newton-meter (for permanent magnet synchronous motor, PMSM) , machine price per Newton-meter (for motors ), machine price per rated power (for generators ), and others [7,8]. When designing any electrical machine, one of the main tasks is obtaining a maximum or minimum one or several of these parameters. When solving problems of maximization of various parameters, different results can be obtained, as shown in .

There is still no single mathematical model of an electrical machine, which solves the problem of maximization for at least one of the quality indicators, as stated in . Various soft computing tools such as neural networks [3,9] and metaheuristics are used to solve such problems. The most effective methods for optimizing electrical machines are various matavistic methods. These are computational methods that provide a fairly good approximate solution to the problem in polynomial time.

Metaheuristics includes several methods, such as Ant Colony, evolutionary optimization, and genetic algorithms [5,10,11,12] as well as some combinations of them.

For example, L.S. Batista et al. in  use the Ant Colony metaheuristic method to solve the problems of optimization of the Interior Permanent Magnet Machine, namely, the task of minimizing the volume of permanent magnets per Newton-meter and the problem of improving the shape of the MMF curve in the air gap. The result of this study is two different topological forms of rotor design, providing solutions to two set tasks. Both structures are technologically feasible. However, a great disadvantage of the developed technique is that it is aimed only at optimizing the shape of the rotor and shape of stator slots but the stator core is not optimized.

A.C.F. Mamede and J.R. Camacho in  use an evolutionary optimization algorithm to maximize electromagnetic torque per total volume of a single-phase switched recovery machine and to decrease electrical losses in copper. The main parameter of motors of similar design is the angle of rotor and stator arc. The result of this study is the design of an engine having the maximum possible torque in a given volume. However, as with any other metaheuristic method, the result was achieved through repeated effort. In this case, the result was achieved in 50 iterations, but it should be noted that the work was carried out through the modernization of the already existing engine, i.e., the initial data for starting the evolutionary algorithm were known at the beginning. It would take many more iterations to calculate an engine without initial data.

P. Virtic and M. Vrazic in  apply genetic algorithms to the optimum design of the axial flux permanent magnet synchronous motor with two stators. Several optimization tasks were solved in this work: (1) maximized weight torque density, (2) maximized volume torque density, (3) minimized mass of PMs per Newton-meter, (4) minimized volume of PMs per Newton-meter, and (5) minimized machine price per Newton-meter. The result of the work was five completely different electric motor designs, solving different tasks. Each of the tasks required about 4000 iterations. With the use of computer equipment available to the authors, the calculation of each iteration took 145 s, and the complete calculation of the engine design took about 13.5 h.

T. Raminosoa and B. Blunier in  also use genetic algorithms to upgrade a high-speed switched reluctation machine in an air compressor.

In addition, there are combined optimization techniques that use neural networks and various metaheuristic algorithms for design. For instance, S. Meo and A. Zohoori in  use a combination of the metaheuristic method “particle swarm optimization” and four neural networks to optimize a permanent magnet flux switching generator for a low-power wind turbine, namely, to minimize the weight and price of the generator, as well as to improve the harmonic composition of the induced voltage. Neural networks are designed to model complex relationships between weight, volume, cost, and the harmonic distortion factor. The metaheuristic algorithm serves to find the best solution in the entire search space. FEM calculations were used to train neural networks. Such combined methods are very effective, however, training neural networks, as well as working out all iterations required a metaheuristic algorithm which takes a long time and computational power.

The main contribution of an article is that we propose a motor design procedure aimed for torque density optimization. Such a conception was not proposed before in any of the publications we know. Currently, the most common optimization approaches are based on metaheuristics or the application of neural networks. They provide an optimal solution, with a high degree of accuracy, but they require multiple iterations to be processed. This requires a lot of processing power and a lot of time. This article proposes a simplified mathematical model of an electric machine, allowing to design electric motors in one iteration with the maximum possible torque density. With this model, the Axial Flux Induction Machine (AFIM) with a short-circuited rotor will be designed in operation. However, this technique can be adapted to any type of a motor.

2. Materials and Methods

In the paper , M.N. Benallal gives the basic principles of the model of an induction motor proposed by M.A. Vaganov. The optimal value of the inner diameter of the stator and flux density in the air gap can be obtained as a result of using this model. Optimization is performed with a constant value of the outer diameter and axial length of the stator steel package. An optimum in the indicated variables ensures the maximum value of the electromagnetic torque. Thus, the model allows us to obtain the maximum value of the electromagnetic torque for a given volume of the machine. However, this model is proposed for the classic performance of an induction motor. It is necessary to formulate it for an axial flux induction motor.

2.1. Thought Experiment

The principle of this optimization can be illustrated by a simple thought experiment. Suppose that a machine occupies a fixed volume, i.e., its outer diameter Do and the total axial length L are unchanged. If, within the framework of a constant axial length L, the length of the stator core is changed by l1, then the length of the rotor core will be determined by the expression l2=Ll1.

Further, we assume that if we increase the length of the stator core l1 (as in Figure 1a) in the limit to the total axial length L, then the length of the steel stack of the rotor l2 tends to zero, and a motor of this design will develop an electromagnetic torque that also tends to zero. On the contrary, if the stator core length l1 (as in Figure 1c) is reduced to zero, then the length of the rotor core l2 will tend to the total length of the machine L, and the torque developed by such a machine will also tend to zero.

Based on Rolle’s theorem, the electromagnetic torque as a function of the length of the stator core must have at least one extremum in the considered interval. Thereby, the basis for the mathematical model is the expression of the electromagnetic torque of the induction motor:Tem=m1U12ωsR2/s(R1+c1R2/s)2+(X1+c1X2)2, where Tem is electromagnetic toque [Nm]; m1—a number of stator phases; ωs—electrical angular velocity; U1—the phase voltage of the stator [V, RMS]; R1—stator resistance [Ohm]; R2—the resistance of the rotor referred to the stator [Ohm]; X1—stator reactance [Ohm]; X2—the rotor reactance referred to the stator [Ohm]; and c1—coefficient taking into account voltage drop at full resistance of stator winding.

2.2. Model Expression Generation

Equation (1) needs to be converted. It is both multiplied and divided by the square root of the product of the phase resistance of the stator and the phase resistance of the rotor referred to the stator:Tem=m1U12ωsR1R2R1R2R2/s(R1+c1R2/s)2+(X1+c1X2)2.

The phase voltage U1 can be represented as the product of the elecromotive force (EMF) E10 [V, RMS] of the self-induction of the stator winding and the coefficient c1 taking into account the voltage drop on the winding resistance. U1=E10c1

In turn, E10=π2f1kwW1Φm, where f1 is the frequency of the voltage of the stator winding; W1—the number of turns connected in series of the stator winding phase; Φm—magnetic flux [Wb]; and kw—winding factor.

If the active length of the stator conductors is taken to be equal to l1=(DoDi)/2 and the pole pitch of the stator is equal to τp=π(Do+Di)/2, then the amplitude of the rotating magnetic flux of mutual flux density in the air gap is Φm=2πτpl1Bδ=Bδ(Do2Di2)4p, where Di is inner diameter [m]; Do—outer diameter [m]; τp—pole pitch [m]; Bδ—magnetic flux density in air gap [T]; and p—the number of pole pairs.

Of greatest interest in Equation (2) is the first coefficient. For the further formation of the expression of the mathematical model, it is necessary to obtain expressions for the resistances of the rotor and the stator windings, as well as expressions for the cross-sectional area of the slots of the stator and rotor.

2.2.1. Stator Resistance

The phase resistance can be defined as the ratio of the doubled number of turns of the winding W1 multiplied by the average length of the half-turn of the coil l1ht and the temperature factor kϑ1 to the specific electrical conductivity γ1 [S/m] of the winding material and the cross-sectional area of the winding conductor qa1 [m2]. R1=2W1l1htkϑ1γ1qa1

In the process of generating an expression for the average half-turn of the stator winding, it must be taken into account that the stator has two coil spans of different lengths along the outer diameter of the stator πDokp/4p and along the inner diameter of the stator πDikp/4p:l1ht=kew2(DoDi)+π4p(Do+Di)kp, where the radial length factor of the end winding kew shows the radial length of the end winding in units of the length of the active part and kp is pitch factor.

We introduce the notation kD=Do+DiDoDi and transform Equation (6):l1ht=k1*DoDi2k1*=kew+πkpkD2p, where coefficient k1* describes the length of the end winding.

Next, we transform Equation (5) for the resistance of the stator winding, taking into account Equation (8):R1=W1(DoDi)k1*kϑγ1qa1.

2.2.2. Rotor Cage Resistance

According to the work in , the resistance of the short-circuited rotor winding related to the stator winding of a conventional cylindrical induction motor is determined by the expression:R2=4m1(kwW1)2z2Rbar+Rring2sin2πpz2, where z2 is the number of the rotor slots.

The length of the bar is defined as lbar=(DoDi)/2. Therefore, the resistance of the rod is Rbar=kϑ2(DoDi)2γ2qbar, where kϑ2 is the thermal factor of the rotor; γ2—conductivity of rotor conductors [S/m]; and qbar—area of the rotor rods [m2].

In the short-circuited winding of the rotor of the end induction motor there must also be two short-circuited rings: the outer ring with diameter Do and the inner with diameter Di. Without significant error, the diameters of the rings can be equal to the outer and inner diameters of the rotor core, respectively, and then Rring=πkϑ2z2γ2qring(Do+Di), where qring is area of short-circuited rings [m2].

Equation (10) can be transformed as R2=4m1(kwW1)2z2kϑ2γ2qbarDoDi2+π(Do+Di)z2qbarqring2sin2πpz2.

The cross-sectional areas of the bar qbar and the short-circuit ring qring can be represented by the ratio of the currents in them Ibar and Iring to the corresponding current densities Jbar and Jring.

The decrease in current density in the short-circuited rotor rings relative to the current density in the rods (to improve the cooling conditions of the rods) can be taken into account using the coefficient kJ=Jring/Jbar. Thus, the ratios qbar and qring can be compiled taking into account the fact that Ibar/Iring=2sin(πp/z2) according to the work in . qbarqring=IbarIringJringJbar=2kJsin(πp/z2)

Substitute this in Equation (11). R2=2m1(kwW1)2kϑ2γ2z2qbar(DoDi)+kJπ(Do+Di)z2sin(πp/z2)

For the final form of expression for R2, (7) should be taken into account. R2=2m1(kwW1)2kϑ2k2*γ2z2qbar(DoDi),k2*=1+πkJkDz2sin(πp/z2)

2.2.3. Slot Cross-Sectional Area

The expression for the stator slot area Qs1 [m2] is the product of the total number of stator winding conductors Nc1=2W1m1 and the cross-sectional area of the conductor without insulation qa1, taking into account the space factor for copper kCu. Qs1=Nc1qa1kCu=2W1m1qa1kCu

A similar expression can be written for the rotor, taking into account the fact that with a squirrel cage rotor, only one bar with a cross section qbar is located in the slot. Then, the total area Qs2 [m2] is required to place z2 rods in it with a space factor for aluminum kAl. Qs2=z2qbarkAl

2.2.4. General Model Expression

Now Equation (2) can be transformed by substituting Equations (3), (4), (9), and (12), taking into account (13) and (14):Tem=πf1kwc12γ1γ2kCukAl32kϑ1kϑ2(Do+Di)(Do2Di2)Bδ2Qs1Qs2k1*k2*1pF.

The last two coefficients of Equation (2) can be denoted as a “parametric factor” F. F=R1R2R2/s(R1+c1R2/s)2+(X1+c1X2)2

The first factor of (15) is a generalized coefficient kF, which includes several quantities that are constant or almost constant, i.e., their numerical values are known or set in accordance with the task for engine design. kF=πf1kwc12γ1γ2kCukAl32kϑ1kϑ2

The remaining part of (15) can be called a “functional factor” FM describing the longitudinal-transverse geometry of an AFIM. This expression will be further transformed taking into account the specific shape of the slot of the stator and rotor. FM=(Do+Di)(Do2Di2)Bδ2Qs1Qs2k1*k2*1p

Moreover, the equation of the model can be written as Tem=kFFMF.

2.2.5. Area of Rectangular Slots

For simplicity, rectangular open slots on the stator are considered. The width of the slot bs1 can be determined from the expression for the area of the slot in the plane of the air gap Q˜s1. If we make the assumption that the arcs of the circles that limit the slot along the outer and inner surfaces of the stator core are equal to the corresponding chords, then bs1=Q˜s1/(DoDi)/2.

In its turn, the stator slot area can be found as the difference between the area of the tooth pitch and the tooth Qs1=QτQt1. The area corresponding to the tooth pitch of the stator (hatched area 1 in Figure 2, see also isometric view in Figure 3) is defined as Qτ=πDo242π/z12ππDi242π/z12π=π4z1Do2Di2, where z1 is the number of the stator slots.

If flux density in the air gap Bδ is evenly distributed within the tooth pitch, then the magnetic flux corresponding to the tooth pitch Φt1=QτBδ=QtBtkCu.

The dependence of the area of tooth pitch and tooth can be found as Qt1=QτBδBt1kFe, where kFe is space factor for iron; Bt1—magnetic flux density in stator teeth [T].

Then, the tooth pitch and stator slot area will be written as Qτ=π4z1Do2Di2BδBt1kFe, Q˜s1=QτQt1=π4z1Do2Di2(1BδBt1kFe).

Now it is possible to determine the width of the stator tooth. bs1=Q˜s1(DoDi)/2=π2z1Do+Di(1BδBt1kFe)

To find the height of the stator core hc1 and tooth hs1, it is necessary to convert the expression for the amplitude of the magnetic flux. It is desirable that only half of the magnetic flux passes through the core. In such a case, the expression for the magnetic flux (Equation (4)) is written as Φm=Bδ(Do2Di2)4p=Bc1(DoDi)2hc1kFe, where Bc1 is magnetic flux density in stator core [T].

After the magnetic flux conversion, it is possible to find the stator core thickness as well as the height of the stator slots. hc1=(Do+Di)4pBδBc1kFe

The height of the stator tooth can be determined as the difference between the thickness of the stator steel stack and the thickness of the stator core. hs1=l1hc1=l1(Do+Di)4pBδBc1kFe

In the simplest case, it can be assumed that the rotor slots are also rectangular, and the outer and inner diameters of the rotor and stator are equal (Do1=Do2=Do;Di1=Di2=Di). Then, the geometric dimensions of the slots and core of the rotor can be found in the same way as the slots and core of the stator, i.e., Equations (20)–(22) simply change indexes from 1 to 2.

2.2.6. Relative Dimensions

In order to obtain a general solution to this problem, we will introduce a system of relative dimensions, taking as the basic total axial length of the engine L. Under this condition, the main relative dimensions of the engine can be written as Do*=Do/L, Di*=Di/L, l1*=l1/L.

In such a case, the geometric dimensions of the stator slots can be converted as bs1=L(π2z1Do*+Di*(1BδBt1kFe)), hs1=L(l1*(Do*+Di*)4pBδBc1kFe).

The total area of all stator slots can be written as:Qs1=Z1bs1hs1=L2ks1, ks1=π2(Do*+Di*)(1BδBt1kFe)(l1*(Do*+Di*)4pBδBc1kFe).

The axial length of the rotor core can be written as l2=Ll1δ, where δ is thickness of air gap. Provided that the rotor slots are also rectangular, their total area is written similarly. Qs2=Z2bs2hs2=L2ks2 ks2=π2(Do*+Di*)(1BδBt2kFe)(l1*(Do*+Di*)4pBδBc2kFe)

Given (25) and (26), the Equation (18) is written as FM=(Do*+Di*)(Do*2Di*2)L5Bδ2Ks1Ks2k1*k2*1p.

This conversion allows us to rewrite the expression for the electromagnetic torque in relative units:Tem=kFFM*FL5, FM*=(Do*+Di*)(Do*2Di*2)Bδ2Ks1Ks2k1*k2*1p.

3. Results 3.1. Calculation Example 3.1.1. Design Setup

As mentioned in Section 1, the developed mathematical model allows us to design optimal AFIM by setting some initial data. In order to illustrate the developed method, we solve the problem of AFIM design with maximum possible electromagnetic torque density in the specified volume, i.e., with hard-limited axial length and outer diameter of stator and rotor steel packages. However, note that by using the same patterns and equations, it is possible to solve the inverse problem, i.e., the task of building an engine with specified torque in minimum volume. In any case, the result of this technique is an electric motor with the maximum possible torque density.

The results below will be for solving the problem of designing a low-voltage AFIM with four pole pairs. Such motors can be widely used in small electric transport [17,18,19,20], for example, in electric bicycles or electric kick scooters, i.e., in devices with severely limited consumption. Initial data for the design are given in Table 1.

3.1.2. Applications of the Analytical Model

The most significant variables defining all the characteristics of the machine are the axial length of the stator, which determines the relationship between the stator and the rotor (see Section 2.1) and flux density in the air gap. These variables mainly determine the value of the electromagnetic torque (19). As these variables are independent of each other, one can find the extremes of this expression from them.

It should be noted that the values of the parametric factor (16) and generalized coefficient (17) do not depend on these variables, i.e., extremes of electromagnetic torque are at the same points as extremes of the expression of functional factor (18).

After substituting the initial data in the model expression, the following dependencies were obtained (Figure 4).

From Figure 4, it can be seen that the electromagnetic torque takes the maximum value at the following values of variables: l1=13.9mm;Bδ=1.13T. These values provide the maximum torque density in the designed engine.

Rotor steel pack thickness can be found as the difference of total engine thickness, stator steel pack thickness, and air gap values. l2=Ll1δ=15.1mm

Substituting in Equations (22) and (20) initial data and obtained values, dimensions of stator and rotor slots can be found geometrically. The following results were obtained hs1=5.4mm;hs2=6.3mm;bs1=3.7mm;bs2=2.8mm. After substitution of all obtained values in all model expressions, the maximum possible value of electromagnetic torque was obtained Tem=0.84Nm.

3.2. Validation of Calculated Results

The calculated engine was validated by using FEM, the analytical method based on Maxwell’s equations, and the measured torque. The external view of the designed engine is shown in Figure 5.

Electromagnetic torque is presented in Figure 6, in the steady state electromagnetic torque pulsations are observed, however this is normal for AFIM . The pattern of magnetic flux density in the air gap is presented in Figure 7.

Figure 7 clearly shows four pairs poles of motor; in addition, it is clear that the amplitude of magnetic flux density in the air gap reaches 1.1–1.2 T, which corresponds to the optimal value of the amplitude of magnetic flux density in the air gap (Figure 4). For AFIM with one stator and one rotor, the amplitude of magnetic flux density in air gap usually does not exceed 1.2 T, as shown in [22,23].

4. Discussion

This article shows the possibility of optimizing an electric machine in a closed volume. Classic electric machine models, such as those based on the output equation , optimize the machine for the inner diameter of the stator without controlling its outer diameter. A distinctive feature of our technique is taking into account the external dimensions of the machine and optimizing the dimensions of the stator and rotor already inside this given volume.

In our opinion, another important scientific result of this work is obtaining a mathematical model of the motor as a design object. This makes it possible to conduct research on the influence of various design parameters on the characteristics of an electric machine. Note that for a more reliable description of the electric machine, the model should be improved. The model must be supplemented with characteristics of steel magnetization, heat transfer and heat removal processes, the ability to set an efficiency class, and other additional functions.

5. Conclusions

The article proposes an analytical model of AFIM, with which it is possible to design an electric motor with the maximum possible electromagnetic torque density. By analyzing the expression of the electromagnetic torque, the model allows us to obtain an optimal value of magnetic flux density in the air gap, as well as geometric parameters of the motor, namely, axial lengths of the stator and rotor steel packs and the height and width of the stator and rotor slots.

In the articles considered in the introduction, similar results are achieved through various soft computing techniques, namely, neural networks and metaheuristic algorithms. A common disadvantage of these methods is that they require working out multiple iterations to achieve the goal. It takes a lot of time and computing power for processing these iterations. The model presented in this article allows us to solve the problem in a single iteration.

To illustrate the operability of the model, we calculated a small low-voltage AFIM with four pairs of poles. In addition, the designed engine was validated using the FEM.

This model not only provides an algorithm for calculating the engine for a single iteration, but is also a valuable research tool. The influence of various parameters on the characteristics of an electric machine can be obtained in the context of its analysis. Therefore, in the following publications, the analysis of the influence of the number of poles and other parameters on the characteristics of the machine will be checked.

Author Contributions

Conceptualization, G.B. and T.K.; data curation, V.O.; formal analysis, G.B. and A.V.; investigation, A.Z. and A.V.; methodology, G.B.; project administration, G.B. and T.K.; resources, A.Z., T.K., and A.V.; software, A.Z. and T.K.; supervision, G.B.; validation, V.O. and A.V.; visualization, A.Z. and V.O.; writing—original draft, G.B.; writing—review and editing, V.O. and T.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript: AFIMAxial flux induction motorFEMFinite element methodpthe number of pole pairsTemelectromagnetic toque [Nm]m1number of stator phasesU1the phase voltage of the stator [V], RMSωselectrical angular velocityR1stator resistance [Ohm]R2the resistance of the rotor referred to the stator [Ohm]c1Coefficient taking into account voltage drop at full resistance of stator windingX1stator reactance [Ohm]X2the rotor reactance referred to the stator [Ohm]E10electromotive force (emf) [V], [RMS]W1the number of series-connected turns of the stator windingΦmmagnetic flux [Wb]Diinner diameter [m]Doouter diameter [m]τppole pitch [m]l1the stator thickness [m]l2the rotor thickness [m]δthickness of air gap [m]Lthe total thickness of the machine [m]Bδmagnetic flux density in air gap [T]Bt1magnetic flux density in stator teeth [T]Bt2magnetic flux density in rotor teeth [T]Bc1magnetic flux density in stator core [T]Bc2magnetic flux density in rotor core [T]kϑ1the thermal factor of the statorkϑ2the thermal factor of the rotorγ1conductivity of stator conductors [S/m]γ2conductivity of rotor conductors [S/m]qa1area of the stator conductor [m2]qbararea of the rotor rods [m2]qringarea of short-circuited rings [m2]Jbarcurrent density of rotor rods [A/m2]Jringcurrent density of short-circuited rings [A/m2]Qτthe area of the coil pith [m2]Q˜s1the area of the stator slot [m2]Qs1the total area of the all stator slots [m2]Qt1the area of the stator tooth [m2]Q˜s2the area of the rotor slot [m2]Qs2the total area of the all rotor slots [m2]Qt2the area of the rotor tooth [m2]kppitch factorkewradial length factor of the end windingkwwinding factorkFespace factor for ironkCuspace factor for copperkAlspace factor for aluminumz1the number of the stator slotsz2the number of the rotor slotsbs1the width of the stator slot [m]bs2the width of the rotor slot [m]hs1the height of the stator slot [m]hc1the height of the stator core [m]hs2the height of the rotor slot [m]

References Dineva A. Mosavi A. Faizollahzadeh Ardabili S. Vajda I. Shamshirband S. Rabczuk T. Chau K.W. Review of soft computing models in design and control of rotating electrical machines Energies 2019 12 1049 10.3390/en12061049 Virtic P. Pisek P. Marcic T. Hadziselimovic M. Stumberger B. Analytical analysis of magnetic field and back electromotive force calculation of an axial-flux permanent magnet synchronous generator with coreless stator IEEE Trans. Magn. 2008 44 4333 4336 10.1109/TMAG.2008.2001528 Krasopoulos C.T. Beniakar M.E. Kladas A.G. Multicriteria PM motor design based on ANFIS evaluation of EV driving cycle efficiency IEEE Trans. Transp. Electrif. 2018 4 525 535 10.1109/TTE.2018.2810707 Virtic P. Pisek P. Hadziselimovic M. Marcic T. Stumberger B. Torque analysis of an axial flux permanent magnet synchronous machine by using analytical magnetic field calculation IEEE Trans. Magn. 2009 45 1036 1039 10.1109/TMAG.2009.2012566 Virtič P. Vražić M. Papa G. Design of an axial flux permanent magnet synchronous machine using analytical method and evolutionary optimization IEEE Trans. Energy Convers. 2015 31 150 158 10.1109/TEC.2015.2477319 Muselli M. Notton G. Louche A. Design of hybrid-photovoltaic power generator, with optimization of energy management Sol. Energy 1999 65 143 157 10.1016/S0038-092X(98)00139-X Bramerdorfer G. Zăvoianu A.C. Surrogate-based multi-objective optimization of electrical machine designs facilitating tolerance analysis IEEE Trans. Magn. 2017 53 1 11 10.1109/TMAG.2017.2694802 Bramerdorfer G. Tapia J.A. Pyrhönen J.J. Cavagnino A. Modern electrical machine design optimization: Techniques, trends, and best practices IEEE Trans. Ind. Electron. 2018 65 7672 7684 10.1109/TIE.2018.2801805 Meo S. Zohoori A. Vahedi A. Optimal design of permanent magnet flux switching generator for wind applications via artificial neural network and multi-objective particle swarm optimization hybrid approach Energy Convers. Manag. 2016 110 230 239 10.1016/j.enconman.2015.11.062 Hejra M. Mansouri A. Trabeisi H. Optimal design of a permanent magnet synchronous motor: Application of in-wheel motor Proceedings of the 2014 5th International Renewable Energy Congress (IREC) Hammamet, Tunisia 25–27 March 2014 1 5 Raminosoa T. Blunier B. Fodorean D. Miraoui A. Design and optimization of a switched reluctance motor driving a compressor for a PEM fuel-cell system for automotive applications IEEE Trans. Ind. Electron. 2010 57 2988 2997 10.1109/TIE.2010.2041133 Knypiński Ł. Pawełoszek K. Le Menach Y. Optimization of Low-Power Line-Start PM Motor Using Gray Wolf Metaheuristic Algorithm Energies 2020 13 1186 Frederico L.S.B.F.C. Min G.G.J.A.R. Ant colony optimization for the topological design of interior permanent magnet (IPM) machines Electron. Eng. 2014 26 1324 1345 Mamede A.C.F. Camacho J.R. Evolutionary algorithms for optimization of 4/4 single phase switched reluctance machine IEEE Lat. Am. Trans. 2018 16 1684 1691 10.1109/TLA.2018.8444387 Benallal M. Vaganov M. Pantouhov D. Ailam E. Hamouda K. Optimal value of air gap induction in an induction motor Proceedings of the XIX International Conference on Electrical Machines—ICEM 2010 Rome, Italy 6–8 September 2010 1 4 Pyrhonen J. Jokinen T. Hrabovcova V. Design of Rotating Electrical Machines John Wiley & Sons Hoboken, NJ, USA 2013 Yang Y.P. Shih G.Y. Optimal design of an axial-flux permanent-magnet motor for an electric vehicle based on driving scenarios Energies 2016 9 285 10.3390/en9040285 Du-Bar C. Design of an Axial Flux Machine for an In-Wheel Motor Application Master’s Thesis Chalmers Reproservice Göteborg, Sweden 2011 Gao P. Gu Y. Wang X. The design of a permanent magnet in-wheel motor with dual-stator and dual-field-excitation used in electric vehicles Energies 2018 11 424 10.3390/en11020424 Kreim A. Schäfer U. An approach to an optimal design of permanent magnet synchronous machines for battery electric vehicles World Electr. Veh. J. 2013 6 673 683 10.3390/wevj6030673 Nobahari A. Darabi A. Hassannia A. Axial flux induction motor, design and evaluation of steady state modeling using equivalent circuit Proceedings of the 2017 8th Power Electronics, Drive Systems & Technologies Conference (PEDSTC) Mashhad, Iran 14–16 February 2017 353 358 Valtonen M. Parviainen A. Pyrhonen J. Influence of the air-gap length to the performance of an axial-flux induction motor Proceedings of the 2008 18th International Conference on Electrical Machines Vilamoura, Portugal 6–9 September 2008 1 5 Nasiri-Gheidari Z. Lesani H. Design optimization of a single-phase axial flux induction motor with low torque ripple Przegląd Elektrotechniczny 2012 88 168 172 Figures and Table

Illustration of a thought experiment. (a) Illustration of the Axial Flux Induction Machine (AFIM) in which the stator occupies 80% of the axial length, (b) the stator occupies 50% of the axial length, and (c) the stator occupies 30% of the axial length.

Stator tooth zone in the plane of the air gap.

Stator teeth in isometric.

Constraints determining electromagnetic torque.

Engine model in isometry.

Electromagnetic torque.

Magnetic flux density in air gap.

wevj-12-00024-t001_Table 1

Initial data of design process.

Designation Name Value Unit
U Line voltage, RMS 7.35 V
U m a x Maximum line voltage 36 V
f Rated frequency 50 Hz
f m a x Maximum frequency 100 Hz
D o Outer diameter 0.11 m
D i Inner diameter 0.05 m
L Total thickness 0.03 m
m 1 The number of phases 3 -
p The number of pole pairs 4 -
q The number of slots per pole and phase 1 -
y 1 Step of span in slot pitches 3 -
z 2 The number of rotor slots 32 -
k C u Space factor for copper 0.5 -
k F e Space factor for iron 0.98 -
k A l Space factor for aluminum 0.97 -
k ϑ 1 Thermal factor 1.32 -
δ Thickness of air gap 1 mm
k m Overload capability 2 -
B c 1 Inductance of stator core 1.4 T
B c 2 Inductance of rotor core 1.3 T
B t 1 Inductance in stator teeth 1.8 T
q a 1 Cros-section of 1 conductor 2.8 × 10 7 m 2

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.