<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Energy Res.</journal-id>
<journal-title>Frontiers in Energy Research</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Energy Res.</abbrev-journal-title>
<issn pub-type="epub">2296-598X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1240473</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2023.1240473</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Energy Research</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Optimal design and control of permanent magnet assisted dual rotor motor</article-title>
<alt-title alt-title-type="left-running-head">Kong et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fenrg.2023.1240473">10.3389/fenrg.2023.1240473</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Kong</surname>
<given-names>Xiaoguang</given-names>
</name>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Yaowen</given-names>
</name>
<uri xlink:href="https://loop.frontiersin.org/people/2326142/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Xu</surname>
<given-names>Gepeng</given-names>
</name>
</contrib>
</contrib-group>
<aff>
<institution>School of Information Engineering</institution>, <institution>Shenyang University of Chemical Technology</institution>, <addr-line>Shenyang</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/333212/overview">Olusegun David Samuel</ext-link>, Federal University of Petroleum Resource Effurun, Nigeria</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1522298/overview">Andrea Toscani</ext-link>, University of Parma, Italy</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/615262/overview">Kenneth E. Okedu</ext-link>, Melbourne Institute of Technology, Australia</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Xiaoguang Kong, <email>kongxiaoguang@syuct.edu.cn</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>26</day>
<month>09</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>11</volume>
<elocation-id>1240473</elocation-id>
<history>
<date date-type="received">
<day>15</day>
<month>06</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>12</day>
<month>09</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Kong, Zhang and Xu.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Kong, Zhang and Xu</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>As high-performance motors, permanent magnet motors are widely used in a wide range of applications. It has become a consensus to mine reluctance torque in permanent magnet motors. The combination of permanent magnet motors and reluctance motors to generate higher output torque is one of the hotspots in motor research. A dual-rotor motor can be formed by connecting a coaxial connector or a concentric end disk, which can make the motor generate higher torque. However, although the motor torque has been improved, the cogging torque still affects the output torque of the motor. This paper describes a method to reduce the cogging torque of the permanent magnet rotor of the permanent magnet-assisted double rotor motor. By analyzing the motor power equation, it is concluded that the pole arc coefficient, the thickness of the magnetic steel, the length of the air gap, and the slot width of the stator have four influences on the teeth. For the parameters of the slot torque, the upper and lower limits of the parameter value are obtained according to the size of the motor. A certain parameter is taken as a fixed value, and the remaining parameters are uniformly valued. Use parametric scanning to determine the optimal value range of the parameter, and use Maxwell for parameterization. Simulation and analysis show that the cogging torque of the motor is reduced by 90% and the torque ripple is reduced by 50%. In order to simplify the motor control system, this paper designs a fuzzy controller based on granular functions, and the fuzzy rules of the fuzzy controller are to perform feature sampling and fit the response function, eliminating fuzzification and defuzzification, improving the response speed of fuzzy control, and simplifying the control system.</p>
</abstract>
<kwd-group>
<kwd>dual rotor</kwd>
<kwd>permanent magnet synchronous motor</kwd>
<kwd>cogging torque</kwd>
<kwd>parameter scan</kwd>
<kwd>granular function</kwd>
</kwd-group>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Process and Energy Systems Engineering</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>The advantages of permanent synchronous motors include their small size, light size, minimal loss, outstanding performance, and wide range of applications (<xref ref-type="bibr" rid="B18">Wang and Leng, 2018</xref>). They are widely used for various occasions and equipment. However, due to the high price and limited supply of rare earth magnets, the development of high-performance motors without rare earths has gradually attracted great attention. Reducing the use of permanent magnets in permanent magnet motors and using synchronous reluctance motors without permanent magnets are feasible alternatives. The synchronous reluctance motor has a reliable mechanical design, a low cost, and a multi-barrier construction that can efficiently enhance the output torque, but its applicability is also constrained by its drawbacks of poor torque density and low factor (<xref ref-type="bibr" rid="B8">Hofmann and Sanders, 2000</xref>).</p>
<p>The two ideal properties of a motor are high torque density and high efficiency (<xref ref-type="bibr" rid="B1">Agbo et al., 2021</xref>). Due to its high torque density and high space efficiency, the double-rotor motor has always been a research focus for scientists. In recent years, in order to further improve the motor efficiency and torque density, scientists from different countries have developed various types of double rotor motors, such as double rotor permanent magnet motors with different pole topologies (NS) (<xref ref-type="bibr" rid="B14">Qu and Lipo, 2003</xref>) and double rotor disc motors with modified stator core yoke (<xref ref-type="bibr" rid="B12">Liu and Li, 2022</xref>), steel sheets instead of part of the permanent magnets reduce the Manufacturing cost of the double rotor axial flux motor (<xref ref-type="bibr" rid="B6">Diao et al., 2023</xref>), and there is also a concentric switched double-rotor reluctance motor that integrates two rotors and two stators into a compact and robust structure (<xref ref-type="bibr" rid="B3">Aravind et al., 2011</xref>; <xref ref-type="bibr" rid="B2">Al-Ani, 2021</xref>), and a switched double-rotor reluctance motor with magnetic isolation ring (<xref ref-type="bibr" rid="B16">Shirzad, 2023</xref>), and yokeless double-rotor mutual coupling switched reluctance motor (<xref ref-type="bibr" rid="B7">Fu et al., 2023</xref>); magnetic field modulation brushless double-rotor machine (<xref ref-type="bibr" rid="B4">Bai et al., 2015</xref>), double-rotor in-wheel motor with multiple operating modes (<xref ref-type="bibr" rid="B10">Li et al., 2023</xref>) Permanent magnet reluctance dual rotor motor (<xref ref-type="bibr" rid="B11">Li et al., 2017</xref>), and so on. In general, it is mainly divided into permanent magnet double rotor motors and reluctance double rotor motors. In this paper, combining the advantages of a permanent magnet motor and a reluctance motor, a permanent magnet-assisted dual-rotor motor is proposed.</p>
<p>In order to achieve the goals of high torque density and wide constant-speed operating range, a permanent magnet assisted dual-rotor motor is proposed in this paper, which adopts a dual-rotor radial flux configuration, combining the synchronous reluctance rotor and the permanent magnet rotor with rare earth permanent magnets to form a single motor entity, the design can control reluctance torque and permanent magnet torque more flexibly, since the cogging torque determines the torque ripple, and the torque ripple affects the quality of the output torque, therefore the goal of the permanent magnet synchronous motorization of the outer rotor is to reduce the cogging torque, according to the definition of the cogging torque, determine the parameters that affect the cogging torque, such as the pole arc coefficient, considering that the parameters have little influence on each other, the parameter sweep method can be used. The range of parameter optimization is given according to the size constraints of the motor, and the cogging torque of the permanent magnet rotor is optimized. Through the simulation analysis of the optimized motor model, the effectiveness of this optimization scheme is verified. Finally, the permanent magnet-assisted dual-rotor motor designed in this paper has the characteristics of a magnetic isolation ring. The motor is split into internal and external parts for control. The internal and external motor control systems share a speed loop and are equipped with hysteresis controllers, respectively. The internal and external motors use torque ratios for current distribution and introduce granular functions to improve the constant fuzzy speed loop, which simplifies the control system and improves the rapidity and robustness of motor response.</p>
</sec>
<sec id="s2">
<title>2 Motor preliminary design</title>
<sec id="s2-1">
<title>2.1 Motor topology</title>
<p>The outer rotor of a permanent magnet assisted dual rotor motor is a surface-mounted permanent rotor, which shares several characteristics with permanent magnet dual rotor motors. The distinction is that the permanent magnet synchronous motor and the synchronous reluctance motor are integrated into a single motor, and the synchronous reluctance rotor is used to replace the internal permanent magnet rotor (<xref ref-type="bibr" rid="B22">Yunyun et al., 2012</xref>). The two rotors share a shaft and rotate at the same mechanical speed. The stator winding adopts the ring winding method, and the magnetomotive force direction of the stator winding of the outer motor is opposite that of the stator winding of the inner motor. In order to achieve magnetic insulation between the inner stator and the outer stator, a magnetic insulator used as a barrier ring is used in the middle of the stator yoke. The insulator is made of brass. <xref ref-type="fig" rid="F1">Figure 1</xref> is a three-dimensional diagram of the motor. The parameters of the motor are shown in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>3D model of motor.</p>
</caption>
<graphic xlink:href="fenrg-11-1240473-g001.tif"/>
</fig>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Main parameters of the motor.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="center">Parameter</th>
<th colspan="2" align="center">Data</th>
</tr>
<tr>
<th align="center">Internal motor</th>
<th align="center">External motor</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Rotor pole number</td>
<td align="center">6</td>
<td align="center">6</td>
</tr>
<tr>
<td align="center">Number of stator slots</td>
<td align="center">36</td>
<td align="center">36</td>
</tr>
<tr>
<td align="center">Stator radius/mm</td>
<td align="center">48.65</td>
<td align="center">63.95</td>
</tr>
<tr>
<td align="center">Outer radius of rotor/mm</td>
<td align="center">36</td>
<td align="center">74.65</td>
</tr>
<tr>
<td align="center">Inner radius of rotor/mm</td>
<td align="center">12</td>
<td align="center">66.65</td>
</tr>
<tr>
<td align="center">Air gap length/mm</td>
<td align="center">0.35</td>
<td align="center">0.5</td>
</tr>
<tr>
<td align="center">Length of iron core/mm</td>
<td align="center">70</td>
<td align="center">70</td>
</tr>
<tr>
<td align="center">Power rating/kw</td>
<td align="center">1</td>
<td align="center">2</td>
</tr>
<tr>
<td align="center">Speed/rpm</td>
<td align="center">1,600</td>
<td align="center">1,600</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>In order to ensure that the motor outputs ideal torque, the internal synchronous reluctance rotor should be properly designed. The typical number of rotor poles of a synchronous reluctance motor is four, but this requires the external surface-mounted permanent magnet rotor to also maintain a four-pole topology. This design method will cause the external permanent magnet rotor yoke to become thicker, and with added weight comes reduced performance. The external surface-mounted permanent magnet rotor is more suitable for the application of multi-rotor pole topology to reduce the thickness of the yoke and the weight of the motor, but this type of topology is not the best design for synchronous reluctance motors because, because of the multi-pole structure, the pole speed of the magnetization inductor of the synchronous reluctance rotor will decrease, and the performance of the synchronous reluctance rotor will also decrease under the condition of a high pole number. Therefore, the number of rotor poles selected for this type of motor is six, and the number of poles of the inner and outer rotors is the same, so as to balance the low number of magnetic poles of the synchronous reluctance rotor and the high number of magnetic poles of the surface-mounted permanent magnet rotor.</p>
<p>The design of the synchronous reluctance rotor flux barrier focuses on the number and position of each flux barrier layer. Generally, for the inner rotor synchronous reluctance motor, the torque ripple decreases with the increase of the number of flux barriers, but due to the limitations of the manufacturing process, the rotor flux barriers are generally set to 4 layers in the design of the synchronous reluctance motor (<xref ref-type="bibr" rid="B9">Jurca and Martis, 2017</xref>), so the number of rotor magnetic barriers in this paper is also set to 4 layers.</p>
<p>Different windings, cogging slots, permanent magnet arrangements, and reluctance rotor magnetic barrier shapes of permanent magnet-assisted dual-rotor motors have various topological structures. The permanent magnet rotor can not only be installed on the surface but also be buried in the ground, and the shape of the magnetic barrier of the reluctance rotor can also be LAL, ARC, segmented rectangle, etc. <xref ref-type="fig" rid="F2">Figure 2</xref> shows the permanent magnet rotor and a few examples of magnetic choke rotors.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Other types of permanent magnet assisted reluctance dual-rotor motor topology: <bold>(A)</bold> Hyperbolic magnetic barrier, <bold>(B)</bold> Plug-in permanent magnet, <bold>(C)</bold> Segmented magnetic barrier, <bold>(D)</bold> ARC barrier.</p>
</caption>
<graphic xlink:href="fenrg-11-1240473-g002.tif"/>
</fig>
</sec>
<sec id="s2-2">
<title>2.2 Motor power equation</title>
<sec id="s2-2-1">
<title>2.2.1 External permanent magnet motor</title>
<p>The rated output power of the motor is the product of the number of phases, voltage, and current of the motor. There must be conversion efficiency in power conversion, and the stator waveform of the motor is a rectangular wave rather than a sine wave in most cases. Therefore, it is necessary to multiply the motor efficiency coefficient <inline-formula id="inf1">
<mml:math id="m1">
<mml:mrow>
<mml:mi>&#x3b7;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> and the shape coefficient <inline-formula id="inf2">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, Therefore, after ignoring the loss caused by the external stator winding impedance and leakage inductance, the output power of the external permanent magnet motor <inline-formula id="inf3">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>p</mml:mi>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> can be expressed as:<disp-formula id="e1">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>p</mml:mi>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x3b7;</mml:mi>
<mml:mi>m</mml:mi>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
</p>
<p>Among them, <inline-formula id="inf4">
<mml:math id="m5">
<mml:mrow>
<mml:mi>&#x3b7;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the motor efficiency; <inline-formula id="inf5">
<mml:math id="m6">
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the number of phases; <inline-formula id="inf6">
<mml:math id="m7">
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the power form factor; <inline-formula id="inf7">
<mml:math id="m8">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the peak value of the back electromotive force of the single-phase winding; <inline-formula id="inf8">
<mml:math id="m9">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the peak value of the stator phase current. Since the stator current vector is aligned with the back EMF vector, the d-axis current is 0, and the value of <inline-formula id="inf9">
<mml:math id="m10">
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> can be set to 0.5.</p>
<p>The peak value of the back electromotive force of the external permanent magnet motor can be expressed as the product of the peak value of the air gap flux density, the effective length of the winding, and the speed at which the winding cuts the magnetic induction line, <inline-formula id="inf10">
<mml:math id="m11">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>B</mml:mi>
<mml:mi>L</mml:mi>
<mml:mi>V</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>. Among them, the peak value of the air gap magnetic field can be expressed by <inline-formula id="inf11">
<mml:math id="m12">
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, and the effective length of the winding can be expressed as the product of the number of turns <inline-formula id="inf12">
<mml:math id="m13">
<mml:mrow>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> of each phase winding, the winding coefficient <inline-formula id="inf13">
<mml:math id="m14">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and the axial length <inline-formula id="inf14">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, since the length of the air gap is very short, for an external rotor motor, the linear velocity of a point on the outer diameter of the stator can be used to represent the linear velocity of a certain point in the air gap, so the linear velocity of a certain point in the air gap can be calculated by the formula <inline-formula id="inf15">
<mml:math id="m16">
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> express, for an external rotor motor, the formula should be expressed as <inline-formula id="inf16">
<mml:math id="m17">
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, so the peak value of the back electromotive force of the external permanent magnet rotor <inline-formula id="inf17">
<mml:math id="m18">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> can be expressed as:<disp-formula id="e2">
<mml:math id="m19">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2217;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
<mml:mo>&#x2217;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x3c0;</mml:mi>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
</p>
<p>Among them, <inline-formula id="inf18">
<mml:math id="m20">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the winding coefficient; <inline-formula id="inf19">
<mml:math id="m21">
<mml:mrow>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the number of winding turns per phase; <inline-formula id="inf20">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the peak value of the surface-mounted permanent magnet rotor side air gap flux density; <inline-formula id="inf21">
<mml:math id="m23">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the power supply frequency; <inline-formula id="inf22">
<mml:math id="m24">
<mml:mrow>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the number of pole pairs; <inline-formula id="inf23">
<mml:math id="m25">
<mml:mrow>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the permanent The outer diameter of the stator of the magneto; <inline-formula id="inf24">
<mml:math id="m26">
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the axial length of the motor core, <inline-formula id="inf25">
<mml:math id="m27">
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the linear velocity of a certain point in the air gap, <inline-formula id="inf26">
<mml:math id="m28">
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the radius of the motor, and <inline-formula id="inf27">
<mml:math id="m29">
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the speed of the motor.</p>
<p>Use the electrical load value <inline-formula id="inf28">
<mml:math id="m30">
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> on the rotor side to represent the peak value of the stator phase current <inline-formula id="inf29">
<mml:math id="m31">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf30">
<mml:math id="m32">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> represents the current effective value of the electrical load, <inline-formula id="inf31">
<mml:math id="m33">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
<mml:mi>m</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula> represents the effective value of current that each phase needs to provide to the load, <inline-formula id="inf32">
<mml:math id="m34">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula> represents the effective value of current that each turn of wire needs to provide to the load, and the current peak value needs to be multiplied by the coefficient <inline-formula id="inf33">
<mml:math id="m35">
<mml:mrow>
<mml:msqrt>
<mml:mn>2</mml:mn>
</mml:msqrt>
</mml:mrow>
</mml:math>
</inline-formula>, so the peak value of the stator phase current can be expressed as <inline-formula id="inf34">
<mml:math id="m36">
<mml:mrow>
<mml:msqrt>
<mml:mn>2</mml:mn>
</mml:msqrt>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>, after adjusting the parameter position, peak current <inline-formula id="inf35">
<mml:math id="m37">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> can be obtained:<disp-formula id="e3">
<mml:math id="m38">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msqrt>
<mml:mn>2</mml:mn>
</mml:msqrt>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
<p>Substituting parameters in Eqs <xref ref-type="disp-formula" rid="e2">2</xref>, <xref ref-type="disp-formula" rid="e3">3</xref> into Eq. <xref ref-type="disp-formula" rid="e1">1</xref>, the output power of the motor can be expressed as:<disp-formula id="e4">
<mml:math id="m39">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>p</mml:mi>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msqrt>
<mml:mn>2</mml:mn>
</mml:msqrt>
<mml:msup>
<mml:mi>&#x3c0;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msubsup>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
</sec>
<sec id="s2-2-2">
<title>2.2.2 Internal reluctance motor</title>
<p>The stator winding must supply the d-axis and q-axis current components to generate the reluctance torque, so in the expression, the formula needs to be multiplied by an included angle coefficient <inline-formula id="inf36">
<mml:math id="m40">
<mml:mrow>
<mml:mi mathvariant="italic">cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b3;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>. The constraint equation of the synchronous reluctance rotor can be compared to the constraint equation of the permanent magnet synchronous motor. The difference is that the synchronous reluctance rotor does not have a permanent magnet that generates back electromotive force. As a result, the synchronous reluctance rotor&#x2019;s output power <inline-formula id="inf37">
<mml:math id="m41">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>y</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> can be stated as:<disp-formula id="e5">
<mml:math id="m42">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>y</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msqrt>
<mml:mn>2</mml:mn>
</mml:msqrt>
<mml:msup>
<mml:mi>&#x3c0;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi mathvariant="italic">cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b3;</mml:mi>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>.</mml:mo>
<mml:mi>R</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>R</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msubsup>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
</p>
<p>Among them, <inline-formula id="inf38">
<mml:math id="m43">
<mml:mrow>
<mml:mi mathvariant="italic">cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b3;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the power factor of the inner rotor motor, <inline-formula id="inf39">
<mml:math id="m44">
<mml:mrow>
<mml:mi>&#x3b3;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the angle between the q-axis component and the stator current; <inline-formula id="inf40">
<mml:math id="m45">
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>.</mml:mo>
<mml:mi>R</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the peak value of the air gap magnetic flux density of the synchronous reluctance rotor side; <inline-formula id="inf41">
<mml:math id="m46">
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>R</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the electrical load value on the rotor side of the synchronous reluctance; <inline-formula id="inf42">
<mml:math id="m47">
<mml:mrow>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the inner diameter of the stator of the synchronous reluctance motor.</p>
</sec>
<sec id="s2-2-3">
<title>2.2.3 Double rotor motor</title>
<p>By combining the constraint equations of the external permanent magnet motor and the internal reluctance motor, the complete constraint equation of the dual-rotor motor can be derived, the output power <inline-formula id="inf43">
<mml:math id="m48">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is expressed as:<disp-formula id="e6">
<mml:math id="m49">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msqrt>
<mml:mn>2</mml:mn>
</mml:msqrt>
<mml:msup>
<mml:mi>&#x3c0;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>.</mml:mo>
<mml:mi>R</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>R</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msubsup>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mi mathvariant="italic">cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b3;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>.</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msubsup>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
</p>
<p>Since the inner and outer motor stator windings of the motor are connected in series, the winding currents are equal, so it can be obtained:<disp-formula id="e7">
<mml:math id="m50">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>R</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>A simple transformation can be obtained<disp-formula id="e8">
<mml:math id="m51">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>R</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
</p>
<p>The relationship between the load on the inner rotor side and the outer rotor side can be expressed by the ratio of the inner diameter of the inner motor to the outer diameter of the outer motor:<disp-formula id="e9">
<mml:math id="m52">
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>R</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>R</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
</p>
<p>Putting <xref ref-type="disp-formula" rid="e9">Formula (9)</xref> into <xref ref-type="disp-formula" rid="e6">Formula (6)</xref>, after simplification, motor total power <inline-formula id="inf44">
<mml:math id="m53">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> can be expressed as:<disp-formula id="e10">
<mml:math id="m54">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msqrt>
<mml:mn>2</mml:mn>
</mml:msqrt>
<mml:msup>
<mml:mi>&#x3c0;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>.</mml:mo>
<mml:mi>R</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi mathvariant="italic">cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b3;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>.</mml:mo>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>R</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mrow>
<mml:msubsup>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mi>L</mml:mi>
</mml:mrow>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
</p>
</sec>
</sec>
</sec>
<sec id="s3">
<title>3 Cogging torque optimization</title>
<sec id="s3-1">
<title>3.1 Cogging torque and optimization method</title>
<p>A major contributor to torque ripple in permanent magnet motors is the cogging torque (<xref ref-type="bibr" rid="B17">Wang et al., 2011</xref>). In order to reduce the influence of the cogging torque on the surface-mount permanent magnet motor, <xref ref-type="bibr" rid="B13">Patel (2023)</xref> recommended slotting to reduce the gear toque by 48.78%. <xref ref-type="bibr" rid="B15">Rashid and Mohammed (2023)</xref> proposed a sawtooth moment mitigation method based on the radial slits of the pole pieces. <xref ref-type="bibr" rid="B20">Won et al. (2023)</xref> used neural network prediction to achieve the purpose of reducing cogging torque. <xref ref-type="bibr" rid="B19">Wang et al. (2023)</xref> adopted the method of combining in-phase unit (IPU) grouping and slotting angle offset so that the main harmonic component of the sawtooth torque is cancelled to achieve the purpose of reducing the cogging torque.</p>
<p>This paper chooses to determine the factors affecting cogging torque from the definition of cogging torque. Cogging torque <inline-formula id="inf45">
<mml:math id="m55">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>g</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is defined as the negative derivative of the magnetic energy <inline-formula id="inf46">
<mml:math id="m56">
<mml:mrow>
<mml:mi>w</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> with respect to the position angle <inline-formula id="inf47">
<mml:math id="m57">
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> when de-energized. The expression can be derived as follows:<disp-formula id="e11">
<mml:math id="m58">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>g</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3c0;</mml:mi>
<mml:mi>z</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi>R</mml:mi>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>&#x221e;</mml:mi>
</mml:msubsup>
<mml:mi>n</mml:mi>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi mathvariant="italic">sin</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mi>z</mml:mi>
<mml:mi>&#x3b1;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
</p>
<p>In the formula, <inline-formula id="inf48">
<mml:math id="m59">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is outside armature diameter of external motor, <inline-formula id="inf49">
<mml:math id="m60">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is stator inner diameter of external motor, <inline-formula id="inf50">
<mml:math id="m61">
<mml:mrow>
<mml:mi>z</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is number of stator slots, <inline-formula id="inf51">
<mml:math id="m62">
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the relative position angle between the stator and rotor, <inline-formula id="inf52">
<mml:math id="m63">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is vacuum permeability.</p>
<p>
<inline-formula id="inf53">
<mml:math id="m64">
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the Fourier coefficient of the residual flux density squared along the circumference of permanent magnets, and the formula is as follows:<disp-formula id="e12">
<mml:math id="m65">
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mn>2</mml:mn>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:msubsup>
<mml:mi>B</mml:mi>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>n</mml:mi>
<mml:msub>
<mml:mi>&#x3b1;</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
</p>
<p>In the formula, <inline-formula id="inf54">
<mml:math id="m66">
<mml:mrow>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the residual flux density in a permanent magnet, <inline-formula id="inf55">
<mml:math id="m67">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b1;</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is polar arc coefficient of permanent magnetic pole.</p>
<p>
<inline-formula id="inf56">
<mml:math id="m68">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the Fourier coefficient of the square relative permeability of the air-gap and the formula is as follows:<disp-formula id="e13">
<mml:math id="m69">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>&#x221e;</mml:mi>
</mml:msubsup>
<mml:mfrac>
<mml:mn>2</mml:mn>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
<mml:mrow>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>sin</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mi>&#x3c0;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mi>z</mml:mi>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
</p>
<p>In the formula, <inline-formula id="inf57">
<mml:math id="m70">
<mml:mrow>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is air gap length, <inline-formula id="inf58">
<mml:math id="m71">
<mml:mrow>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the thickness of the permanent magnet along the magnetization direction, <inline-formula id="inf59">
<mml:math id="m72">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the corresponding radian value of the width of the stator groove.</p>
</sec>
<sec id="s3-2">
<title>3.2 Determine optimization variables and parameter ranges</title>
<p>It can be seen from the above formula that the slot opening angle width, magnetic steel angle width, magnetic steel thickness, and air gap length all have a great influence on the cogging torque. Therefore, these four parameters are selected as the cogging torque optimization parameters in this paper to study their effect on cogging torque. <xref ref-type="table" rid="T2">Table 2</xref> lists the optimized variable parameter ranges.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Optimize the parameter range.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Parameter</th>
<th align="center">The lower limit</th>
<th align="center">The upper limit</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Polar arc coefficient</td>
<td align="center">0.6</td>
<td align="center">0.76</td>
</tr>
<tr>
<td align="center">Stator slot width/mm</td>
<td align="center">1.2</td>
<td align="center">3</td>
</tr>
<tr>
<td align="center">Permanent magnet thickness/mm</td>
<td align="center">1</td>
<td align="center">2.6</td>
</tr>
<tr>
<td align="center">Air gap length/mm</td>
<td align="center">0.2</td>
<td align="center">1.6</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3-3">
<title>3.3 Parameter optimization</title>
<p>In order to explore the law of cogging torque weakening caused by different design variables, each design variable is parameterized in the optimization process, and one parameter is selected as a variable. The upper and lower limits are given as shown in <xref ref-type="table" rid="T2">Table 2</xref>. Other parameters are constants; after obtaining the waveforms of the cogging torque in various situations, select the maximum value to draw the trend of the cogging torque changing with the parameters, and select the optimal parameters for simulation verification in Maxwell. If the optimal effect meets expectations, the optimization is complete. If the effect is not as expected, it is necessary to change the upper and lower limits of the parameters or to optimize the parameters. The specific optimization scheme flow chart is shown in <xref ref-type="fig" rid="F3">Figure 3</xref>. As we can see in <xref ref-type="fig" rid="F4">Figure 4</xref>, the optimal values of the four parameters can be easily found according to the graph.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Flow chart of cogging torque optimization scheme.</p>
</caption>
<graphic xlink:href="fenrg-11-1240473-g003.tif"/>
</fig>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>The maximum value of cogging torque when each parameter changes: <bold>(A)</bold> Maximum cogging torque at different stator slot widths, <bold>(B)</bold> Maximum cogging torque at different air-gap lengths, <bold>(C)</bold> Maximum cogging torque at different magnet thicknesses, <bold>(D)</bold> Maximum cogging torque at different polar arc coefficients.</p>
</caption>
<graphic xlink:href="fenrg-11-1240473-g004.tif"/>
</fig>
</sec>
<sec id="s3-4">
<title>3.4 Optimization results</title>
<p>During the cogging torque parameterization, different degrees of attenuation can be demonstrated by using different methods and changing the parameter range. It can be seen from <xref ref-type="fig" rid="F4">Figure 4</xref> that when the pole arc coefficient is between 0.72 and 0.74 and the thickness of the magnetic steel is between 1.6&#xa0;mm and 1.8&#xa0;mm, the calculated cogging torque is relatively lower. Considering that increasing the air gap will reduce the output torque of the motor, the length of the air gap should be kept between 0.6&#xa0;mm and 1&#xa0;mm. The cogging torque waveform before and after adjustment is shown in <xref ref-type="fig" rid="F5">Figure 5</xref>. The cogging torque is reduced from 1.5&#xa0;Nm to 0.15&#xa0;Nm, and the optimized output torque waveform is more stable, the torque ripple is lower, and the optimization effect is more obvious.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Waveform comparison before and after optimization: <bold>(A)</bold> Cogging torque waveform before optimization, <bold>(B)</bold> Cogging torque waveform after optimization, <bold>(C)</bold> Output torque waveform before and after optimization.</p>
</caption>
<graphic xlink:href="fenrg-11-1240473-g005.tif"/>
</fig>
</sec>
</sec>
<sec id="s4">
<title>4 Electromagnetic characteristic analysis</title>
<sec id="s4-1">
<title>4.1 Magnetic flux density</title>
<p>After the optimization is completed, use the finite element software Maxwell to establish a two-dimensional model of the motor, set the rated current on the inner and outer stator windings, and obtain the magnetic flux density cloud diagram of the permanent magnet assisted dual-rotor motor with load through finite element simulation. As shown in <xref ref-type="fig" rid="F6">Figure 6</xref>, the maximum value of the magnetic flux density appears at the position of the inner and outer stator yoke, which is about 1.8T, which does not exceed the design limit of 2T and meets the motor operation specification.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Motor magnetic field density map.</p>
</caption>
<graphic xlink:href="fenrg-11-1240473-g006.tif"/>
</fig>
</sec>
<sec id="s4-2">
<title>4.2 Air gap magnetic density</title>
<p>The air-gap magnetic field density is an important indicator of motor design. <xref ref-type="fig" rid="F7">Figure 7</xref> is the air-gap flux density curve of the internal and external motors under no-load and rated current. The amplitude of the air-gap flux density reflects the magnetic field strength and torque density of the motor, it can be seen that the waveform is closer to a sine wave. From the comparison of the two curves, we can see that the air gap flux density curve of the inner motor is more reasonable than the sinusoidal distribution of the outer motor, indicating that the inner motor can meet the requirements of torque output and power density. The peak values of the air gap magnetic density of the outer motor and the inner motor are 0.8T and 1&#xa0;T, respectively, which generally meet the design requirements.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Air gap magnetic: <bold>(A)</bold> Outer motor air gap flux density and <bold>(B)</bold> Inner motor air gap flux density.</p>
</caption>
<graphic xlink:href="fenrg-11-1240473-g007.tif"/>
</fig>
</sec>
<sec id="s4-3">
<title>4.3 Output torque</title>
<p>Since the motor has a double-air-gap structure, its total output torque is the sum of the output torques of the inner and outer motors, and the torque output waveforms of the inner and outer motors are shown in <xref ref-type="fig" rid="F8">Figure 8</xref>. According to the output torque results, the output torque of the inner rotor synchronous reluctance motor is about 4&#xa0;<inline-formula id="inf60">
<mml:math id="m73">
<mml:mrow>
<mml:mtext>Nm</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula>, and the torque ripple is about 28.8%. The output torque of the outer rotor permanent magnet synchronous motor is about 13&#xa0;<inline-formula id="inf61">
<mml:math id="m74">
<mml:mrow>
<mml:mtext>Nm</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula>, and the torque ripple is about 17.1%, which verifies the effectiveness of the optimized design, and the quality of the motor output torque has been significantly improved.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Output torque waveform of internal and external motors.</p>
</caption>
<graphic xlink:href="fenrg-11-1240473-g008.tif"/>
</fig>
</sec>
</sec>
<sec id="s5">
<title>5 Design of control system for dual rotor motor</title>
<sec id="s5-1">
<title>5.1 Motor mathematical model</title>
<p>The dual-rotor motor is composed of inner rotor, outer rotor, permanent magnet, magnetic isolation ring and other parts. It has the characteristics of strong coupling and nonlinearity. For the convenience of analysis, the assumptions are as follows: 1) Neglect eddy current and hysteresis loss; 2) Excitation current There is no response time; 3) The stator windings are star-connected, and the stator windings are connected in parallel. According to the above assumptions, it can be obtained that the dual-rotor motor can be split into a permanent magnet synchronous motor and a synchronous reluctance motor, and after transformation by <inline-formula id="inf62">
<mml:math id="m75">
<mml:mrow>
<mml:mtext>clark</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf63">
<mml:math id="m76">
<mml:mrow>
<mml:mtext>park</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> the voltage equation of the motor under the <inline-formula id="inf64">
<mml:math id="m77">
<mml:mrow>
<mml:mtext>dq</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> axis.<disp-formula id="e14">
<mml:math id="m78">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>q</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi mathvariant="italic">cos</mml:mi>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mi mathvariant="italic">cos</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mi mathvariant="italic">cos</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="italic">sin</mml:mi>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="italic">sin</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="italic">sin</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>B</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>
<disp-formula id="e15">
<mml:math id="m79">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>q</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi mathvariant="italic">cos</mml:mi>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mi mathvariant="italic">cos</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mi mathvariant="italic">cos</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="italic">sin</mml:mi>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="italic">sin</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="italic">sin</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>B</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(15)</label>
</disp-formula>
</p>
<p>Since the double-rotor motor is divided into two parts, the inner rotor and the outer rotor, the electromagnetic torque of the motor is expressed as the permanent magnet torque of the outer motor and the reluctance torque of the inner motor, and the sum of the inner and outer torques is the total electromagnetic torque of the motor.<disp-formula id="e16">
<mml:math id="m80">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>3</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mi>q</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>q</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>3</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>q</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>&#x3c8;</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(16)</label>
</disp-formula>
<disp-formula id="e17">
<mml:math id="m81">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>J</mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>&#x3c9;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>B</mml:mi>
<mml:mi>&#x3c9;</mml:mi>
</mml:mrow>
</mml:math>
<label>(17)</label>
</disp-formula>
</p>
<p>In these formulas, <inline-formula id="inf65">
<mml:math id="m82">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the logarithm of the motor poles, <inline-formula id="inf66">
<mml:math id="m83">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> <inline-formula id="inf67">
<mml:math id="m84">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>B</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> <inline-formula id="inf68">
<mml:math id="m85">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the motor stator three-phase current, <inline-formula id="inf69">
<mml:math id="m86">
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> <inline-formula id="inf70">
<mml:math id="m87">
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>q</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is current of <inline-formula id="inf71">
<mml:math id="m88">
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>q</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> axial, <inline-formula id="inf72">
<mml:math id="m89">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> <inline-formula id="inf73">
<mml:math id="m90">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>q</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is voltage of <inline-formula id="inf74">
<mml:math id="m91">
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>q</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> axial, <inline-formula id="inf75">
<mml:math id="m92">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is angle from <inline-formula id="inf76">
<mml:math id="m93">
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> axial to <inline-formula id="inf77">
<mml:math id="m94">
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>q</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> axial, <inline-formula id="inf78">
<mml:math id="m95">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the electromagnetic torque of a motor, <inline-formula id="inf79">
<mml:math id="m96">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is a load torque, <inline-formula id="inf80">
<mml:math id="m97">
<mml:mrow>
<mml:mi>B</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the damping coefficient, <inline-formula id="inf81">
<mml:math id="m98">
<mml:mrow>
<mml:mi>&#x3c9;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the mechanical speed of a motor, <inline-formula id="inf82">
<mml:math id="m99">
<mml:mrow>
<mml:mi>J</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the moment of inertia of the motor, <inline-formula id="inf83">
<mml:math id="m100">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c8;</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the permanent magnet linkage of a motor. <inline-formula id="inf84">
<mml:math id="m101">
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> equals to 1 is the reluctance part of the motor, <inline-formula id="inf85">
<mml:math id="m102">
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> equals to 2 is the permanent magnet portion of the motor.</p>
</sec>
<sec id="s5-2">
<title>5.2 Vector control</title>
<p>Vector control is based on the control idea of a DC motor, converts the three-phase current into excitation current and torque current, and realizes indirect control of the motor by changing the magnitude of the excitation current and torque current, that is, by controlling the axis current of the motor to achieve the purpose of controlling the motor. At present, the commonly used vector control methods of permanent magnet synchronous motors include <inline-formula id="inf86">
<mml:math id="m103">
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> control, <inline-formula id="inf87">
<mml:math id="m104">
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="normal">&#x3c6;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> control, maximum torque-current ratio control and flux weakening control; vector control methods of synchronous reluctance motor include maximum torque-to-current ratio control, maximum power factor control, maximum torque change rate control, and so on. The external motor of the double-rotor motor is a surface-mounted permanent magnet synchronous motor, and the <inline-formula id="inf88">
<mml:math id="m105">
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> control is equivalent to the maximum torque-to-current ratio control, therefore, the inner and outer motors of the dual-rotor motor in this paper are controlled by the maximum torque-current ratio, which simplifies the control system.</p>
</sec>
<sec id="s5-3">
<title>5.3 Torque distribution</title>
<p>Due to the existence of the magnetic isolation ring, the dual-rotor motor can be equivalent to a surface-mounted permanent magnet synchronous motor and a synchronous reluctance motor. Conventional control requires two sets of controllers that interact with each other. In order to solve this problem, torque distribution is introduced. Torque distribution is used to fix the ratio of the output torque of the internal and external motors, distribute the torque current of the internal and external motors, and simplify the control system.</p>
</sec>
<sec id="s5-4">
<title>5.4 Controller design</title>
<sec id="s5-4-1">
<title>5.4.1 Current controller design</title>
<p>The current controller in this design adopts the mode of hysteresis control, also called bang-bang control, which has the characteristics of real-time control, fast reaction, and strong robustness. The implementation method is to change the switching state of the inverter by comparing the given current value with the detected current value, and then output the actual current in the shape of a sawtooth wave. Although there are certain errors, this current control method is simple and easy to implement and does not rely on motor parameters.</p>
</sec>
<sec id="s5-4-2">
<title>5.4.2 Speed controller design</title>
<p>The speed control of this design adopts the granular function instead of the fuzzy controller. The granular function is based on the theory of granular computing. Granular computing (GrC) is usually defined in an informal way as a general computing theory that can be used effectively. To build an efficient computational model for complex applications with large amounts of data, information, and knowledge (<xref ref-type="bibr" rid="B5">Bargiela and Pedrycz, 2008</xref>), the granular function shifts from a machine-centered approach to a human-centered approach. It is one of the trends in granular computing research (<xref ref-type="bibr" rid="B21">Yao, 2010</xref>), and fuzzy set theory is an application of granular computing (<xref ref-type="bibr" rid="B23">Zarandi Baghini et al., 2023</xref>). The construction of granular functions in fuzzy control is divided into steps such as information sampling, information granulation, and fitting. Using granular functions instead of fuzzy controllers can save steps such as fuzzification and defuzzification. The construction process of a granular computing fuzzy controller is shown in <xref ref-type="fig" rid="F9">Figure 9</xref>.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Construction process of granular function optimization fuzzy controller.</p>
</caption>
<graphic xlink:href="fenrg-11-1240473-g009.tif"/>
</fig>
</sec>
</sec>
<sec id="s5-5">
<title>5.5 System simulation and analysis</title>
<sec id="s5-5-1">
<title>5.5.1 Control system block diagram</title>
<p>According to the description of the fuzzy control optimized by the above granular function, a dual-rotor motor control system frame can be built as shown in <xref ref-type="fig" rid="F10">Figure 10</xref>. The control system is composed of a traditional PID controller, a granular function optimization module, a torque distribution module, a current hysteresis controller, and a motor body.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Control system block diagram.</p>
</caption>
<graphic xlink:href="fenrg-11-1240473-g010.tif"/>
</fig>
</sec>
<sec id="s5-5-2">
<title>5.5.2 Simulation result analysis</title>
<p>As shown in <xref ref-type="fig" rid="F11">Figure 11</xref>, the highest speed for fuzzy control and granular function control is about 1,640 r/min, and after changing the load, fuzzy control and granular function control motors drop to 1,591 r/min. It can be concluded that the granular function is comparable to the fuzzy control in terms of control performance, but the granular function removes the complicated fuzzification and defuzzification processes, only needs one response function, and the simulation time is less than one-tenth of the fuzzy control.</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Motor output waveform: <bold>(A)</bold> The output waveforms of different control modes, <bold>(B)</bold> Partial enlarged view of the speed increase process. <bold>(C)</bold> Partial enlarged view of the speed drop process after.</p>
</caption>
<graphic xlink:href="fenrg-11-1240473-g011.tif"/>
</fig>
</sec>
</sec>
</sec>
<sec sec-type="conclusion" id="s6">
<title>6 Conclusion</title>
<p>This paper aims to address the problem that the cogging torque of a dual-rotor motor connected by concentric end disks affects the motor&#x2019;s output torque. By deriving the motor power equation constraints, the fixed-value scanning method is used for optimization analysis. After optimization, the motor cogging torque was reduced from 1.5&#xa0;Nm to 0.15&#xa0;Nm, a reduction of 90%; the external rotor torque ripple was reduced by 50%; the torque output waveform was improved; and the optimized electromagnetic performance was verified. Finally, in order to simplify the motor control system, this paper uses granular functions to improve the fuzzy controller, performs feature sampling on fuzzy rules, and uses response functions to replace the fuzzy controller, omitting the fuzzification and defuzzification links. Simulation shows that after granular function optimization, the response speed of the system is significantly improved.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s7">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/<xref ref-type="sec" rid="s12">Supplementary Materials</xref>, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s8">
<title>Author contributions</title>
<p>YZ and GX wrote and typed the manuscript, and XK provided writing guidance, checking and correction.</p>
</sec>
<sec id="s9">
<title>Funding</title>
<p>This work is supported in part by the National Natural Science Foundation of China under grant 51877139 and the Basic Scientific Research Project of Liaoning Provincial Department of Education of China under grant LJKMZ20220778.</p>
</sec>
<sec sec-type="COI-statement" id="s10">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s11">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s12">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fenrg.2023.1240473/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fenrg.2023.1240473/full&#x23;supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Table1.XLSX" id="SM1" mimetype="application/XLSX" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Agbo</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Samuel</surname>
<given-names>O. D.</given-names>
</name>
<name>
<surname>Amosun</surname>
<given-names>S. T.</given-names>
</name>
<name>
<surname>Oyejide</surname>
<given-names>O. J.</given-names>
</name>
<name>
<surname>Fayomi</surname>
<given-names>O. S. I.</given-names>
</name>
<name>
<surname>Bamisaye</surname>
<given-names>O. S.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Development of a solar energy-powered surface water pump</article-title>. <source>IOP Publ.</source> <volume>1107</volume> (<issue>1</issue>), <fpage>012002</fpage>. <pub-id pub-id-type="doi">10.1088/1757-899X/1107/1/012002</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Al-Ani</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Multi-physics design and analyses of dual rotor synchronous reluctance machine</article-title>. <source>eTransportation</source> <volume>8</volume>, <fpage>100113</fpage>. <pub-id pub-id-type="doi">10.1016/j.etran.2021.100113</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Aravind</surname>
<given-names>C. V.</given-names>
</name>
<name>
<surname>Norhisam</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Aris</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Marhaban</surname>
<given-names>M. H.</given-names>
</name>
<name>
<surname>Ahmad</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Nirei</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2011</year>). &#x201c;<article-title>Double-rotor switched reluctance machine (DRSRM): fundamentals and magnetic circuit analysis</article-title>,&#x201d; in <conf-name>2011 IEEE Student Conference on Research and Development</conf-name> (<publisher-name>IEEE</publisher-name>), <fpage>294</fpage>&#x2013;<lpage>299</lpage>.</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bai</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>A new magnetic-field-modulated brushless double-rotor machine</article-title>. <source>IEEE Trans. Magnetics</source> <volume>51</volume> (<issue>11</issue>), <fpage>1</fpage>&#x2013;<lpage>4</lpage>. <pub-id pub-id-type="doi">10.1109/tmag.2015.2445917</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bargiela</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Pedrycz</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Toward a theory of granular computing for human-centered information processing</article-title>. <source>IEEE Trans. Fuzzy Syst.</source> <volume>16</volume> (<issue>2</issue>), <fpage>320</fpage>&#x2013;<lpage>330</lpage>. <pub-id pub-id-type="doi">10.1109/tfuzz.2007.905912</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Diao</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Ding</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2023</year>). <source>Design and analysis of low-cost double rotor axial flux permanent magnet synchronous motor</source>.</citation>
</ref>
<ref id="B7">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Fu</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Lv</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Si</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2023</year>). <source>Analysis and comparison of novel yokeless double rotor mutually coupled switched reluctance motor</source>.</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hofmann</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Sanders</surname>
<given-names>S. R.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>High-speed synchronous reluctance machine with minimized rotor losses</article-title>. <source>IEEE Trans. Industry Appl.</source> <volume>36</volume> (<issue>2</issue>), <fpage>531</fpage>&#x2013;<lpage>539</lpage>. <pub-id pub-id-type="doi">10.1109/28.833771</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Jurca</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Martis</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2017</year>). &#x201c;<article-title>Analysis of outer rotor synchronous reluctance motor for low-speed applications[C]</article-title>,&#x201d; in <conf-name>2017 19th International Conference on Electrical Drives and Power Electronics (EDPE)</conf-name>, <conf-loc>Dubrovnik, Croatia</conf-loc>, <conf-date>04-06 October 2017</conf-date> (<publisher-name>IEEE</publisher-name>), <fpage>242</fpage>&#x2013;<lpage>247</lpage>.</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Ren</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Coordinated control strategy for drive mode switching of double rotor in-wheel motor based on MPC and control allocation</article-title>. <source>World Electr. Veh. J.</source> <volume>14</volume> (<issue>5</issue>), <fpage>132</fpage>. <pub-id pub-id-type="doi">10.3390/wevj14050132</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Bobba</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Sarlioglu</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Design and optimization of a novel dual-rotor hybrid PM machine for traction application</article-title>. <source>IEEE Trans. Industrial Electron.</source> <volume>65</volume> (<issue>2</issue>), <fpage>1762</fpage>&#x2013;<lpage>1771</lpage>. <pub-id pub-id-type="doi">10.1109/tie.2017.2739686</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Research on concentrated winding counter-rotating double rotor disk motor</article-title>. <source>Mech. Des. Manuf.</source> (<issue>03</issue>), <fpage>110</fpage>&#x2013;<lpage>113</lpage>. <pub-id pub-id-type="doi">10.19356/j.cnki.1001-3997.20211115.003</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Patel</surname>
<given-names>A. N.</given-names>
</name>
</person-group> (<year>2023</year>). <source>Slot opening displacement technique for cogging torque reduction of axial flux brushless DC motor for electric two-wheeler application</source>.</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qu</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Lipo</surname>
<given-names>T. A.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Dual-rotor, radial-flux, toroidally wound, permanent-magnet machines</article-title>. <source>IEEE Trans. industry Appl.</source> <volume>39</volume> (<issue>6</issue>), <fpage>1665</fpage>&#x2013;<lpage>1673</lpage>. <pub-id pub-id-type="doi">10.1109/tia.2003.818968</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rashid</surname>
<given-names>M. K.</given-names>
</name>
<name>
<surname>Mohammed</surname>
<given-names>A. M.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>A reduction method of cogging torque for magnetic gears</article-title>. <source>Iran. J. Electr. Electron. Eng.</source> <volume>19</volume> (<issue>2</issue>), <fpage>2752</fpage>. <pub-id pub-id-type="doi">10.22068/IJEEE.19.2.2752</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Shirzad</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2023</year>). <source>Calculation of flux density in air-gap for reluctance motor with two ports (Double-Stator, double-rotor) by fourier series</source>.</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Qiao</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Pei</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Jung</surname>
<given-names>S. Y.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Reducing cogging torque in surface-mounted permanent-magnet motors by nonuniformly distributed teeth method</article-title>. <source>IEEE Trans. Magnetics</source> <volume>47</volume> (<issue>9</issue>), <fpage>2231</fpage>&#x2013;<lpage>2239</lpage>. <pub-id pub-id-type="doi">10.1109/tmag.2011.2144612</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Leng</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2018</year>). &#x201c;<article-title>Summary on development of permanent magnet synchronous motor</article-title>,&#x201d; in <conf-name>2018 Chinese Control And Decision Conference (CCDC)</conf-name> (<publisher-name>IEEE</publisher-name>), <fpage>689</fpage>&#x2013;<lpage>693</lpage>.</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>An in-phase unit slot-opening shift method for cogging torque reduction in interior permanent magnet machine</article-title>. <source>Mathematics</source> <volume>11</volume> (<issue>7</issue>), <fpage>1735</fpage>. <pub-id pub-id-type="doi">10.3390/math11071735</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Won</surname>
<given-names>Y. J.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>J. H.</given-names>
</name>
<name>
<surname>Park</surname>
<given-names>S. H.</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>J. H.</given-names>
</name>
<name>
<surname>An</surname>
<given-names>S. M.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>D. Y.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Transfer learning-based design method for cogging torque reduction in PMSM with step-skew considering 3-D leakage flux</article-title>. <source>IEEE Trans. Magnetics</source>, <fpage>1</fpage>. <pub-id pub-id-type="doi">10.1109/tmag.2023.3294601</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yao</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Novel developments in granular computing: applications for advanced human reasoning and soft computation: applications for advanced human reasoning and soft computation</article-title>. <source>IGI Glob</source>.</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yunyun</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Xiaoyong</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Hua</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Electromagnetic performance analysis of double-rotor stator permanent magnet motor for hybrid electric vehicle</article-title>. <source>IEEE Trans. Magnetics</source> <volume>48</volume> (<issue>11</issue>), <fpage>4204</fpage>&#x2013;<lpage>4207</lpage>. <pub-id pub-id-type="doi">10.1109/tmag.2012.2206374</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zarandi Baghini</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Babaei</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Tabatabaei Mirhosseini</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Torkzadeh Tabrizi</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>New approach based on fuzzy hypergraphs in granular computing (an application to the urban vulnerability assessment)</article-title>. <source>Int. J. Nonlinear Analysis Appl.</source> <pub-id pub-id-type="doi">10.21203/rs.3.rs-1167172/v1</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>