Analytical Design of Sculpted Rotor Interior Permanent Magnet Machines

: A computationally efﬁcient design of interior permanent magnet (IPM) motor rotor features is investigated utilizing analytical methods. Over the broad operating range of IPM machines, interactions of MMF sources, permeances, and currents result in torque harmonics. The placement of traditional rotor features along with sculpt features are utilized to minimize torque ripple and maximize average torque. We extend the winding function theory to include the IPM rotor’s primary and secondary reluctance paths and the non-homogeneous airgap of the rotor sculpt features. A new analytical winding function model of the single-V IPM machine is introduced, which considers the sculpted rotor and how this model can be used in the design approach of machines. Results are validated with ﬁnite elements. Rotor feature trends are established and utilized to increase design intuition and reduce dependency upon the lengthy design of experiment optimization processes.


Introduction
The IPM motor is increasingly being utilized throughout industry as a primary source of propulsion due to its good efficiency, torque and power density. Examples include the development of battery electric vehicle traction motors [1,2], plugin hybrid electric vehicles [3], and hybrid electric vehicles [4]. Ideally, the traction machine provides an average torque produced from a sinusoidal distribution of the airgap flux density. In reality, embedding the magnet within the salient structure of the rotor lamination and distributing windings in discrete locations result in airgap flux density harmonics. These harmonics result in increased torque ripple, radial forces, losses, and other unwanted phenomena.
In this paper, an approach to minimize torque ripple with rotor features is presented, based upon analytically modeling the machine features. The analytical modeling approach enables efficient use of computational resources, without the sacrifice of harmonic content, prior to the use of more expensive finite element methods. The calculation of the IPM machines' spatially-dependent torque harmonics is performed through the extension of the winding function method. New to the winding function framework is a method to model the equipotential nature of the rotor's salient features and rotor surface modifications. The non-homogeneous airgap of rotor surface modifications is included in the model through an additional MMF term. Unique to this analytical method, both the constituents and aggregates of the torque harmonics are found. A detailed investigation into the rotor geometry design space to minimize torque harmonics while managing average torque design trade-off is presented.
Inherent to the design of IPM machines, torque ripple is a persistent problem. Design choices to increase torque density or decrease manufacturing cost are often at odds with minimizing the torque ripple [5]. Rotor features, including surface modifications or sculpt The third and final feature type is sculpting the rotor surface at the airgap to redirect flux. The single flat magnet IPM machine cogging torque was reduced in [29] and experimentally verified. A grid on/off optimization of the rotor surface was conducted on the single flat magnet IPM using finite elements in [30,31], resulting in an asymmetric rotor surface with reduction in torque ripple and maintaining average torque. The double V-magnet IPM torque ripple was minimized with both rotor core and surface sculpted features in [27]. The delta-magnet IPM machine torque ripple was minimized with rotor surface sculpt features in [32]. Then, a general analytical expression for torque harmonics was developed and utilized to optimize the solution with finite elements. This paper presents a detailed analytical model of the sculpted rotor IPM machine. The model allows for a break down of flux and torque into magnet, primary reluctance, secondary reluctance, and sculpt features. Multiple sculpt features configurations are demonstrated to achieve similar torque harmonic reductions. Results are validated with finite elements and utilized to improve the torque harmonic characteristics of an existing industrialized machine. This is because finite elements accurately predict experimental results across a broad range of machines [33][34][35][36] and have been utilized to evaluate and compare machine types [37][38][39] and validate analytical solutions [40][41][42][43][44][45][46][47][48][49][50]. Section 2 introduces the topics relevant to the design of IPM motor construction and control. Section 3 provides details on how to model an IPM motor magneto motive force (MMF), permeance, and linear current density in order to model the machine geometry, flux, and torque harmonics. The model developed in Section 3 is applied to that of a well-known industrial IPM machine in Section 4. Design features are explored in Section 5. Contributions of this paper include a novel analytical winding function-based IPMSM model, the analytical description of rotor sculpt features, and modeling of magnet and reluctance torque component alignment due to asymmetric sculpt features. In addition, this paper demonstrates the torque effects of reluctance path pole arc, sculpt feature type (symmetrical/asymmetrical), sculpt feature location, and sculpt feature depth and sculpt feature width.

Flux Distribution and Control of IPM Machines
Performance, harmonics, and control are all dependent upon the distribution of flux density within the machine. Permanent magnets provide a constant source of flux density, which enable efficient torque production but can limit high speed operation. Reluctance features provide a source of torque dependent upon armature current at high current angles, useful for extending operation at high speeds. The ratio of magnet and reluctance torque is balanced to enable the machine to stay within its operation constraints while efficiently using the voltage and current available at the terminals. This section provides a brief overview of the machine's construction, design features, and flux paths. In addition, the necessary framework for control is introduced. Figure 2 shows a two-dimensional illustration of a four pole IPM machine. The IPM motor is fundamentally constructed of a stator and a rotor. The stator is the mechanically grounded part of the machine. It is constructed of slots, teeth, a yoke, and the three phase windings. The stator teeth and yoke are constructed of a magnetically permeable iron alloy. The teeth and yoke allow for easy flow of magnetic flux to and from the airgap of the machine. The slots allow space for the copper windings. The windings are distributed within the slots to produce a current dependant magneto motive force (MMF), which in turn creates the radial magnetic flux density. The placement of the windings also creates a current density along the bore of the stator, resulting in a tangential component of flux density. When arranged and controlled properly, the currents in the windings produce a rotating set of fields to produce torque. The rotor, the mechanically rotating part, is constructed of iron ribs, air-filled barriers, iron bridges, and permanent magnets. The iron core is constructed of ribs and bridges. The ribs control the distribution of flux density while the bridges mechanically couple all parts together. Barriers provide air pockets, which assist the ribs in directing flux, and also contain embedded magnets. The permanent magnets are embedded within the rotor and produce an MMF, which is independent of current.

Flux Distribution
From the perspective of the rotor, a direct axis (d-) and a quadrature axis (q-) of the machine are electromagnetically aligned to the rotor characteristics. The daxis is the primary axis of which the permanent magnet flux density flows. The magnet flux density flows through the magnet into the central rib of the magnet pole, into the airgap, through the stator teeth and yoke and returns into the adjacent opposite magnet pole. This permanent magnet flux density path is shown in Figure 2 as a red ellipse. The qaxis is the axis in which the armature-induced flux flows through the rotor. This flux is produced from the armature MMF and the reaction of the rotor/stator permeance. Two armature-induced paths result; one through the primary reluctance path, and a second through the secondary reluctance path. The primary reluctance path flux is shown as the solid blue ellipse, and the second reluctance path flux is shown as the dashed blue ellipse in Figure 2.

Control of IPM Machine
Control must be considered in the design of the IPM motor. The steady state torque and voltage equations are shown in Equations (1)-(3). These are the fundamental starting points to develop the necessary analysis for the control of electric machines. The equations are based upon dand qaxis voltages, v d and v q , currents, i d and i q , inductance λ d and λ q , magnet flux linkage, λ m , and phase resistance R s . Magnet offset, δ, as shown in Figure 3, is included to account for magnet alignment relative to the reluctance path, which may be caused by rotor sculpting features [51,52]. Traditional IPM alignment would feature δ = 0 • , with the d-axis aligned to the maximum of magnet flux linkage. For purposes of this paper, the q-axis remains aligned to the minimum reluctance of the rotors first reluctance path. By inspection, the resistance or loss terms do not have an effect on the torque and only affect the voltage and electrical power. It may be useful to assume the phase resistance is negligible.
Performance assessments require including an analysis of the phase current constraint, I max , and phase voltage V max . For a wye connected machine, I max is equal to the phase current, I ss , and is limited by the power devices of the inverter and the electric machines thermal capability. The voltage limit is the maximum phase voltage that the inverter can apply, limited by the specific pulse width modulation (PWM) technique used. The voltage limit for space vector PWM is . The maximum modulation index is set to MI = 0.95 to account for cable and device voltage drops.

Analytical Model: MMF, Permeance, Flux, and Torque
This section develops the necessary analytical winding function model for the IPM machine idealized to focus on the effects of the rotor geometry. Assumptions include closed stator slots and no saturation leading to infinite permeability, leaving the permeability of the airgap assumed to be that of free space µ • = 4π10 −7 H m . Stator conductors are modeled by discrete current sheets along the stator bore inner diameter and phases are assumed to be wye connected.
Focused on the MMF interaction with the second reluctance path, the analytical model describes MMF and permeance functions. These winding function-based MMFs, F x (θ, φ), and permeance functions, Λ(θ, φ), express the harmonic content of stator and rotor features as Fourier series. Relationships between the rotor spatial coordinate θ and rotor position φ and current angle β are included. Flux densities are computed using Equation (4), and contributions of the stator and rotor harmonic interactions to the torque ripple are determined.

Permeance Functions
The general form of the permeance functions Λ x , shown in Figure 4, can be written as Equation (5), where the amplitude Y and phase γ define the location of features relative to the d-axis of the machine.
Permeance functions, along with rotor and stator construction, are shown in Figure 4. The salient features of the rotor begin with the definition of the primary reluctance path, which assumes a small airgap l g , aligned with the minimum reluctance of the q-axis, and a large airgap considering barrier and magnet dimensions, l m . Above the magnet, a secondary reluctance path exists, which reacts as equipotential salient iron to the armatures MMF. The permanent magnet permeance path describes the total amount of air the magnet must push its flux through, including its thickness in the same region of the second reluctance path. Coefficients of Equation (5) are derived from the local definition of the permeance Equation (6).
The defining airgaps of the permeance functions are listed in Table 1.

Permeance Term Minimum Airgap Maximum Airgap
First Reluctance Path l g l g + l m Second Reluctance Path l g ∞ Permanent Magnet Path l g + l m ∞ Figure 4. Conductor locations and permeability functions, (red) phase A conductor, (blue) phase B conductor locations, (green) phase C conductor location, (gold) first reluctance path permeability function, (lavender) second reluctance path permeability functions, (purple) magnet path permeability function, (gray) equivalent permanent magnet conductor locations.

Magneto Motive Forces (MMF)
Equation (7) gives the general form of the winding function composed of the turns function n(θ, φ) and its mean < n(φ, θ) >, where the turns function is a result of the closed path integral of the conductor [53]. MMF functions, F(θ, φ), are comprised of the turns function multiplied by a current.
Expanded into a Fourier series, the winding takes the form as follows (9): where Y w f x is the coefficient, γ w f x (n) is the phase, mechanical order n, rotor position φ, and initial rotor position φ i .

Stator
Stator MMF is formed from the interaction of the phase currents (I a , I b , I c ) and the phase winding functions (N a , N b , N c ). The summation of the three phases creates a rotating MMF, F abc , which directly acts upon the primary reluctance path.

Magnet
The permanent magnet is represented as an equivalent current I PM , which is related to the remnant flux B r , permeability of free space u • , relative permeability µ r , and magnet thickness. Magnet MMF for the IPM machine interacts with the magnet path permeance, where τ m is the magnet's salient iron pole pitch, and w m is the combined width of the magnets for a single pole.

Second Reluctance Path Modification
The equipotential nature of the rotor's second reluctance path reacts only to the regional harmonics of the stator MMF. In this case, all looping flux that enters the second reluctance path pole arc through the airgap must exit through the same airgap. Modification to the armature MMF by removing its mean satisfies this condition and is made possible through Equation (13). The symbols F <abc> and < F abc (θ, φ)N PM (θ) > represent the modified MMF, which interacts with the second reluctance path and the mean of the MMF across this same boundary. As a matter of convenience, the permanent magnet winding function N PM is also used to consider the stator MMF in the region of the second reluctance path, invert it, remove the mean, and revert to the original polarity.

Sculpt Feature Description: Equivalent Magnetizing Dipole Current
MMF-permeance methods allow for flux density for the IPM smooth rotor homogeneous airgap to be calculated. Rotor sculpting affects both magnet flux distribution and reluctance flux distribution. This smooth rotor MMF-permeance theory does not adequately describe stator slots and rotor sculpt features as it has been developed with a constant airgap dimension [53]. In this section, an extension of winding function theory is developed, which can be used in the description of both slots and sculpts based upon equivalent magnetic currents (EMC) [54] and the equivalent magnetic dipole [55]. These non-homogeneous airgap features are represented by additional MMF terms utilizing the description of a magnetic dipole and its equivalent magnetic currents. The redistribution of flux density and MMF is possible with the use of the equivalent dipole concept.
The magnetic dipole in free space is formed by a loop of radius b and current of I. The solution at far fields, when R >> b, solved in spherical coordinates, using the magnetic vector potential, A, is shown in Equation (15), where the magnetic dipole moment m is written as m = a z Iπb 2 [55].
This dipole in free space can be used to explain the magnetism at the atomistic level, where small circulating currents are formed by the process of magnetization. This magnetization aligns the individual atomic dipoles and modifies the orbital spin of the electrons for each atom.
The macroscopic volume density of magnetization, M, with units of A/m, is computed through a sum of the individual microscopic dipoles. Shown in [55], the magnetization vector M is equivalent to both a volume current density, J m with units of A m 2 , and a surface current density J ms with units of A m .
J ms = M × a n (17) Given M, the flux density B can be found by computing both J m and J ms . These values are used to determine the magnetic vector potential A. Uniform M within a magnetic material will result in no volume current density and only a surface current density J ms on its borders. If space variations of M exist within a material, a net volume current density will exist. Hence, a magnetic dipole inside a material with constant magnetization M can be represented by a current loop in the air, formed at the exterior boundary material. Figure 5 illustrates the process of analyzing the rotor sculpt feature effects. Rather than account for the changing flux density over the sculpt feature, the assumption of homogeneous flux density holds when breaking the geometry into smaller discrete dipoles (i) of fixed width. For purposes of this analysis, it is assumed the sculpt features have a constant depth l ms . A sufficient number of points (i) must be defined in order to hold the assumption of homogeneous flux density. For each point (i) contained within the sculpt feature, the first current I ms1 (i) is applied based upon prior analysis of the magnetic dipole.
A second dipole counter current in the adjacent point I ms2 (i + 1) = −I ms1 (i). The net effective dipole current for each point, I ms3(i) is formed through summation dipole currents I ms3(i) = I ms1(i) + I ms2(i) . A third dipole current I ms3 , or summed current, becomes the current-turns function for the sculpted feature, in which an equivalent MMF for the rotor surface features can be determined through the use of a winding function. It should be noted that the sculpt features analyzed with this process do not create flux but only distribute flux away from the sculpt feature.

Flux Density
Total airgap flux density, B tot , is a result of the flux creating sources and flux distributing features. Primary reluctance path, secondary reluctance path, and the permanent magnets create and distribute flux. The sculpt and slotting features serve to redistribute flux, with the assumption of small features. Flux densities can be computed for each individual component or in summation.

Torque
The Maxwell stress tensor, Equation (20), allows for torque computation given the radial flux density and the conductor's linear current density. Torque due to individual components of flux density can be separated or combined for the net effects.

Rotation and Convolution
In this section, a computationally efficient method to evaluate the phenomenon of rotation and the ensuing interaction of the time and space domain harmonics is reviewed. To simplify analysis of the magnetic fields, the rotor is left frozen while the stator rotates in the counter direction. The Fourier spatial coefficients of the permeance and winding functions are determined at the initial position of the rotor through use of the Fast Fourier Transform (FFT). For each rotor position, φ(i), a complex rotation, Equation (21), is used to form a rotation vector, Equation (22). The complex rotation vector is used to transform the FFT coefficients of permeance and winding functions at each position.
In place of the convolution of the Fourier coefficients to determine flux density, the permeance and winding function is reconstructed with the inverse fast Fourier transform.

Comparative Analysis to Recent Analytical Methods
This section will compare the analytical methods developed in this paper to three recent methods [45,49,50].
The analytical MEC model is used to design a reduced magnet cost single V consequent pole (CP) machine with the same average torque as single V IPM in [45]. The models are developed based on zones and regions, allowing for an assumption of the open circuit flux density distribution. Two flux sources (magnets) and six reluctance network paths are used to create the open circuit single V IPM model. Open circuit flux density is assumed to take a trapezoidal waveform where the flux density is determined spatially from the fluxes in the MEC reluctance network. Similarly, the single V CP IPM network consists of one magnet flux source and five reluctance paths. In both cases, the rotors reluctance path reaction to the armature loading is not computed. The open circuit flux densities are used to quickly determine an equivalent single V CP IPM fundamental to that of the single V IPM fundamental. In order to guide design, the method is used to find an equivalent consequent pole open circuit flux density fundamental to the traditional IPM flux density. Finite elements are relied on to complete the study of the torque performance.
Multi-barrier synchronous reluctance and Permanent Magnet Assisted Synchronous Reluctance Machines (PMSynRM) are modeled using conformal mapping and magnetic equivalent circuits in [50]. The slotless U-type IPM machine open circuit flux density is analytically modeled with a subdomain method solving Laplace's and Poisson's equations in [49]. Analytical equations are derived and presented for each subdomain. Results are validated against finite elements. The model is divided into four regions, which consist of the airgap and magnets. The governing system of partial differential equations is developed, along with simplifications, interface and boundary conditions. A separation of variables is used to develop the general solutions of the PDEs, and they are written as a Fourier series. The system of equations is solved and compared to FE. Strong agreement of the radial and tangential flux density is shown between the FE and the subdomain methods.
The analytical models discussed were developed for multiple purposes. The MEC method is used in [45] to quickly estimate the open circuit flux density fundamental of the single-V IPM and single-V CP IPM machines. The armature reaction of the reluctance features is not considered by the model, and finite elements are used to finish the designs. Conformal mapping is used in [50] to determine the reluctances of a multi-barrier PMSynRM and further evaluated using a MEC network. Both open circuit and loaded conditions are evaluated for airgap flux densities and torque performance and compared to finite elements. The analysis is not extended to the design. The subdomain methodology is employed in [49] for the analysis of the U-shape IPM machine open circuit conditions. Both tangential and radial flux densities are shown to match finite element results. The model requires further extension to consider the torque performance due to a loaded armature. The single-V sculpted rotor IPM winding function model developed in this paper considers both open circuit and loaded conditions. The model is constructed such that computational efficiency is possible without the sacrifice of spatially dependent harmonic content to drive the design. The choice of which analytical model to develop is dependent upon the application, computational resources, and intended purpose. In this case, the new winding factor IPM model was developed.

Application of Analytical Model to Example Machine
The method developed in Section 3 is validated with a well-known industrialized IPM machine. Details of the 2004 Toyota Prius traction motor are included in Table 2, and the geometry is modeled analytically within Matlab and finite elements within Ansys Maxwell. Focusing on the effects of rotor, the stator geometry has been idealized with no slots. Both sinusoidally distributed stator windings and the production configuration of discretely placed windings are modeled. Only the stator winding harmonics interactions with rotor geometry harmonics are considered. With the assumption of infinite permeability, the bridge features are omitted. Airgap fringing in the region of the magnet barrier is not considered.
Design parameters are studied within this section using the analytical winding function model previously validated. Rotor sculpt features are included along with their additional MMF term developed in Section 3.

Sculpting Geometry
The sculpted rotor IPM machine geometry design space to be explored is shown on a single pole of the example machine in Figure 6. Rotor primary and secondary reluctance paths are shown in green with no bridge features. The stator, shown in gray, continues to have omitted its slot features, and distributed windings, orange, are placed within the airgap. The primary design parameters are centered on the rotor effects, which include the ratio of primary and secondary reluctance path and sculpt features. The magnet pole arc width is varied. Up to two sets of symmetrically placed sculpt features are placed on the second reluctance path. The symmetrical feature span locations τ 1 and τ 2 define the symmetrical location of the feature in terms of its percentage of the magnet pole arc span. Single asymmetrical features are described with a similar parameter but with only one sculpt feature on the pole. In this single asymmetrical case, the feature location, τ, is set to be positive for right hand side placement and negative for the left hand side placement. The depths D 1 and D 2 are measured from the outer surface of the rotor to the root of the sculpt feature. The widths W 1 and W 2 are measured in terms of a single feature percentage of the pole span.
A Fourier analysis of design feature effects on torque harmonic analysis is presented. The phase of the torque harmonics is set with respect to the negative zero crossing of phase A back EMF, i.e., the north pole of the machine. Torque components are separated in terms of total, primary reluctance, secondary reluctance and magnet torque with both the torque amplitude and its corresponding phase.

Model Implementation
Implemented in Matlab, the analytical model in Section 3 has been utilized to explore the design space. Rotor permeances of primary reluctance, secondary reluctance, and magnet and the permanent magnet winding function are modeled spatially as the rectangular waves. Stator phase winding functions and linear turn densities are spatially modeled, and the spatial harmonic coefficients are determined with the FFT. For each rotor position, the complex rotation vector is created and applied to the stator winding and linear turn functions. With the current applied, the modified secondary reluctance path is determined from Equation (13). Radial flux density is now determined for the primary reluctance, second reluctance, and magnet paths. A broad range of phase current, control angle, and rotor sculpting features are studied. Maxwell's stress tensor in conjunction with Equation (19) provides the ability to separate components of torque.
The analytical model is implemented utilizing both the FFT and the inverse FFT algorithms to efficiently move between frequency and spatial domains. Rotation is best performed within the frequency domain, Equation (22), and convolution is performed within the spatial domain. This provides the most efficient use of resources by a factor of five times. Results are computed over an entire electrical cycle with a sufficient number of points to provide a smooth torque waveform. The analytical model is executed within Matlab in 6.7 s, and the finite element results are executed in 20 min. Ansys Maxwell was also used to perform the finite element analysis.

Model Validation: Radial Flux Density
This section compares the radial flux density results of the analytical model and finite elements while varying: (1) winding type, (2) current, (3) control angle, and (4) rotor sculpt features. Figures 7-10 plot the flux density along the rotors spatial coordinate, θ, over a single pole pair. Both sinusoidal and distributed windings are compared. The q-axis, which is aligned to the rotor minimum reluctance, occurs at θ elec = 90 • and θ elec = 270 • . The d-axis is aligned, which is aligned to the smooth rotors permanent magnet maximum flux linkage, occurs at θ elec = 0 • , θ elec = 180 • , and θ elec = 360 • . In all results, the finite element and the analytical model result in comparable flux densities.
Flux densities shown for sinusoidal windings, Figure 7, illustrate the changing airgap flux density harmonics with current and control angle. When the phase current is set to zero, only the permanent magnet field is present. As current and current angle increases, the flux density becomes more jagged, with the case of a fully negative d-axis current displaying the most harmonic content. It is clear that as the negative d-axis current becomes dominant, so do the reluctance path harmonics. Harmonic effects of the discretely distributed windings are shown in Figure 8. As current increases, so do the airgap reluctance harmonics. In all cases, the analytical model and finite element results agree with reasonable accuracy.  The reluctance flux density single sculpt feature is plotted in Figure 10. Similar to the symmetrical sculpt features, the flux density is reduced in the region of the sculpt feature. In all cases, the analytical model and finite element results agree with reasonable accuracy.

Model Validation: Torque Ripple
Torque ripple of the smooth rotor IPM, Table 2, is compared between finite element and the analytical model in Figure 11. Good agreement between the finite element and analytical models is observed. The torque ripple effects of two symmetrical rotor sculpt features are demonstrated in Figure 12, directly calculated by the analytical model, whereas the finite element model requires two runs, once with and once without sculpt features, to determine the sculpt feature effects.
Good correlation between the model and finite elements is demonstrated and shown.

Torque Ripple Components
In this section, the torque ripple results of the analytical model are studied while varying current and control angle. Figure 13 plots the torque ripple for a complete electrical cycle of the example machine. The torque components for the first reluctance path, second reluctance path, total reluctance torque, magnet torque, and total machine torque are plotted. Magnet torque and its harmonics are dominant at lower currents, but the reluctance paths cannot be ignored. As current is increased, the reluctance torque increases relative to the magnet torque. The stronger field weakening currents cause the contribution of the reluctance features torque to increase. The dominant torque harmonic orders are the 6th and 12th electrical orders. In the design, both the torque harmonic amplitudes and phases of each of the components need to be considered as the sculpt feature design will provide the counter torque at the counter phase.

Investigation of Design Features
In this section, the effects of design features are demonstrated to influence both torque harmonic amplitudes and phases. Carefully applied, these effects are used to design counter torque harmonics. Second reluctance path pole arc, sculpt feature type (symmetrical/asymmetrical), sculpt feature location, sculpt feature depth, and sculpt feature width can all be used to design an appropriate counter torque to reduce the machine's torque harmonics. While mildly affecting average torque, the second reluctance path pole arc, τ p , strongly affects the phase of the 12th electrical order torque harmonic. A single pair of symmetrical sculpt features placed upon the second reluctance path pole arc reduce the average torque. Feature position provides 12th electrical order torque harmonic phasing, and the feature width and depth directly affect the torque harmonics amplitude. The single asymmetrical feature is shown to increase average torque when placed on a specific side of the second reluctance path pole arc. The asymmetrical feature placement can also be used to modify the phase of both 6th and 12th electrical order torque harmonics. Finally, feature phasor summation is shown to be effective in combining the effects of multiple design features, further providing the ability to design both the amplitude and phase of these minimizing torque harmonics. These relationships provide the necessary intuition to reduce computationally intensive design steps.

Magnet Pole Arc
Magnet pole arc span, τ m , effects upon the torque harmonics, without rotor sculpt features, are explored. Figure 14 shows the average, 6th, and 12th harmonics of torque as a function of τ m . For this case, magnet torque is dominant. Although not always the case, it is just as important to follow the trends of individual torque components. Total, first reluctance, and magnet average machine torque are reduced as the pole arc is increased, and only the second reluctance torque increases the average torque. The 12th order torque harmonic is dominant, with primary contributions from the magnet and the second reluctance path, whereas the 6th order torque harmonic is mostly contributed to by the primary reluctance and magnet. Rotor geometry has a strong influence on the phase of the 12th order torque harmonics, whereas the 6th harmonic is less affected by rotor geometry.

Single Pair Symmetrical Rotor Sculpt Feature
A single symmetrical rotor sculpt feature torque is studied in Figures 15-17. Only the effects of the sculpt feature torque are plotted. In Figures 15 and 16, a single sculpt feature position is varied, with fixed width, W 1 , and fixed depth, D 1 , along the magnet pole arc. Rotor sculpt features have a negative effect upon average torque, as the phase is 180 • out of phase with the smooth rotor average torque. Figure 15 compares the analytical model to finite elements and shows precise agreement with the phase and matching trends for torque amplitude. Using a wider feature width, D 1 , Figure 16 translates amplitude and phase plots to a phasor representation. The 12th harmonic is clearly the dominant torque in both amplitudes and choice of phase.   Figure 17 includes the sculpt features width and depth effects. Increasing sculpt feature width, W 1 , and/or the sculpt features depth, D 1 , increases the amplitude of the sculpt features torque harmonic. The primary influence on the torque harmonic amplitude is the width of the sculpt feature. Sculpt feature width and depth have no effect on the torque harmonics phase.

Single Asymmetrical Rotor Sculpt Feature
A single sculpt feature resulting in asymmetrical placement upon the rotor surface is studied in this section. Similar parameters D 1 , W 1 , and τ 1 are used to describe the features width, depth, and location. In the asymmetrical case, the location, τ 1 , is described with the same location parameter, where in this case, a positive τ 1 results on the right side of Figure 6 and a negative τ 1 results in sculpt feature placement on the left hand side. Figure 18 compares the model to finite element and shows precise agreement with phase and matching trends for torque amplitude. In the asymmetric sculpt feature case, a torque improvement is possible due to the aligned axis effect from τ 1 > 0. The placement of the sculpt feature allows for placement of the torque harmonic phase angle across a broad range of phases. Negative values of τ 1 result in the largest amplitudes of the 12th electrical torque harmonic.

Two Symmetrical Rotor Sculpt Features
More than one symmetrical rotor sculpt feature can be used. In this section, it is shown that the components of a first symmetrical feature can be combined with that of the second symmetrical rotor sculpt feature. The MMF-permeance model is validated by comparing to finite element results in Figure 12.
To illustrate this concept, the parameters of the two sculpt feature sets of Figure 6 are shown in Table 3. Through vector summation, the two vectors were used to create a 12th order counter torque with a phase of −116 • .  Figure 19 illustrates the first (red) and second 12th order electrical torque (green) phasors. The two phasors combine to create the effective total phasor (blue). This phasor summation is plotted along with the torque complex mapping of the previous single feature design sweep. These symmetrical rotor sculpt features are designed to mitigate the 12th order electrical torque harmonics to near zero. A single feature or multiple features can be designed to minimize the torque ripple. The sculpt features are not without consequence, as the average torque is negatively affected.

Conclusions
This paper has presented an analytical modeling and design approach to reduce torque ripple with rotor sculpt features. By carefully placing rotor sculpt features and rotor barrier features, average torque can be maintained while minimizing torque harmonics. Contributions of this paper include: • A new analytical winding factor modeling approach for the single V IPM machine relating the rotor's first reluctance feature, second reluctance feature, permanent magnet features, sculpt features, and stator windings to the resulting torque harmonics; • An analytical modeling approach accounting for both symmetrical and torque aligning asymmetrical rotor sculpt features; • Results from the analytical model providing valuable insights for identifying rotor feature design improvements; • Design approach for placement of rotor sculpt features to minimize torque ripple while maintaining average torque; • Demonstration of close agreement of radial flux density and torque harmonics results between the analytical model and that of finite element results.
These results enable better design insight and an efficient design process through use of an analytical model.