Modeling and design optimization of a shaft-coupled motor and magnetic gear

: This paper presents the modeling and design of an actuator consisting of an electrical motor and a magnetic gear. To minimize the overall actuator dimensions, both of the electromagnetic devices need to be optimally designed and matched. An issue in performing a simultaneous design as such arises from a high number of design variables that signiﬁcantly increases the complexity of the optimization problem. A method to reduce the design variables is discussed in this paper, which is the application of response surface methodology (RSM) to represent the optimized torques of the electrical motor and magnetic gear as polynomial functions of their respective dimensions. Prior to the application of RSM, optimization problem statements are deﬁned for the electrical motor and magnetic gear, for which the optimization objective and constraint functions are derived from analytical electromagnetic models of the considered electromagnetic devices.


Introduction
A common type of actuator used in robotic applications is a mechanically-geared electrical motor, which has a high torque within a relatively small volume. For specific applications such as soft robotics [1], however, it is not desirable to use mechanical gears due to the presence of static friction that degrades the force/torque control performance [2]-not to mention that mechanical gears are vulnerable in the event of torque overload.
A promising solution to address the previous shortcomings of mechanical gears is to replace them with magnetic gears [3], whose topology is depicted in Figure 1. In a magnetic gear, the torque between the rotating members (high-speed and low-speed rotors) are transmitted without mechanical contact, resulting in a major reduction of the static friction. Furthermore, the rotating members would slip when the load torque exceeds a certain level, giving an inherent overload protection.
There are two possibilities to integrate an electrical motor and a magnetic gear. The first one is a conventional shaft-coupled motor and gear as shown in Figure 2a. The second one is the so-called "pseudo" direct-drive motor [4] illustrated in Figure 2b, which integrates all the electromagnetic components belonging to an electrical motor and a magnetic gear. In this paper, the shaft-coupled motor and gear is the chosen one to be investigated as it has a smaller outer diameter, which is a desirable aspect in the considered application.
In robotic/mechatronic applications, the optimal selection of motor and (mechanical) gear transmission ratio is necessary for the overall system performance [5]. For this particular optimization problem, generally, a number of off-the-shelf candidate motors and gears are considered, from which a motor-gear pair is selected based on criteria such as dimension, mass and efficiency [5][6][7].
However, when stringent design requirements are considered, it would be preferable to design and optimize both motor and gear specifically for the intended application, as it would provide flexibility in selecting design parameters such as outer dimensions.
This paper presents a rapid method to simultaneously optimize an electrical motor and a magnetic gear for an application. For this purpose, high-accuracy analytical electromagnetic models of an electrical motor and a magnetic gear are developed based on the harmonic modeling method [8,9]. The models are used to derive the objective and constraint functions from optimization problem statements for the considered electromagnetic devices. Since the optimization of the motor and magnetic gear are treated simultaneously, there is a high number of design variables that are involved, resulting in a slow and computationally-intensive numerical optimization task. For that reason, a solution is proposed in this paper to reduce the number of design variables by representing the optimized torques of the motor and magnetic gear as polynomial functions of their respective dimensions, i.e., outer diameter and axial length. The polynomial functions are approximated by the application of response surface methodology (RSM) [10] on the optimization routines that are defined for the motor and magnetic gear. Based on these polynomial models and the reduced design variables, the simultaneous design optimization of the motor and magnetic gear can be performed rapidly.

Harmonic Modeling Method
The goal of applying the harmonic modeling method [9] is to obtain an accurate estimation of the static magnetic field behavior within different material regions of an electromagnetic device, from which quantities such as magnetic flux density and torque, can be calculated. Note that the following assumptions are adhered in the modeling of the rotary electromagnetic devices discussed in this paper:

1.
The electromagnetic problem can be described in a 2D polar coordinate system (r, ϕ) 2.
For a given region, the material has linear and homogenous magnetic properties in the r-direction.

3.
The ferromagnetic material is infinitely permeable. Consequently, no analytical expression of the magnetic flux density can be obtained within the ferromagnetic material.
As with the modeling of any electromagnetic devices, Maxwell's equations govern the field behavior. From magnetostatic Maxwell equations, the following Poisson equation [8] is obtained where A = A z is the magnetic vector potential, µ 0 is the vacuum magnetic permeability, µ is the material permeability, M 0 = [M r , M ϕ ] is the residual magnetization vector and J = J z is the current density vector. Considering a 2D polar coordinate system, the Poisson Equation (1) can be expressed as Within the investigated rotary electromagnetic devices, (circular) spatial periodicities are exhibited by the source terms in Equation (2), allowing them to be expressed by Fourier series where N is the number of harmonics to be modeled, M r0 , M ϕ0 and J z0 are DC components of the corresponding source terms (M r , M ϕ , J z ); M rs , M rc , M ϕs , M ϕc , J zs and J zc are Fourier coefficients of the source terms, τ k is half of the tangential width (ϕ-direction) of a material region k (within the electromagnetic device) and ϕ k is the tangential position within the region k. Due to Equations (3)-(5), the solution of the Poisson Equation (2), A z , appears as a Fourier series. Given that B = ∇ × A (Gauss' law of magnetic field), the magnetic flux density in the r and ϕ-directions for a given material region can also be expressed as Fourier series where B rs , B rc , B ϕs , B ϕc and B ϕ0 are obtained from the solution of the Poisson Equation (2), given by: B ϕc = a n r where a n , b n , c n , d n and A 0 in Equations (8)- (12) are obtained by considering the boundary conditions between the different material regions of the modeled electromagnetic device. The solutions for G rs , G rc , G ϕs and G ϕc are available in [8].

Harmonic Modeling of an Electrical Motor
The selected motor topology in this paper is the fractional-slot concentrated-windings PM motor [11] shown in Figure 3, which has favorable characteristics such as high torque density and low cogging [12]. For modeling purpose, the motor structure is represented by Figure 4. The structure is characterized by the material regions (I, II, etc.) described in Table 1 and parameters r 0 , ...r 7 . Note that the slot winding regions IV 1 ,...,IV Q are simplified from its actual shape in Figure 3, as illustrated in Figure 5, in which w t denotes the tooth width. The equivalent current density of the winding J eq in Figure 5 is calculated as where S actual and S simple are the cross-section areas of the actual and simplified slot winding, respectively, w t shown in Figure 5 is the tooth width, and J 1 and J 2 denote any adjacent current densities within a slot.    To obtain the expression for the magnetic flux density distribution throughout the investigated motor, Equations (6) and (7) are defined for each of the regions (Table 1) and the corresponding unknown coefficients in Equations (8)- (12) are solved through a set of linear equations arising from the following boundary conditions: •

Neumann boundary condition
This specifies the tangential magnetic field strength to be H ϕ = 0 [8] at the boundary between two regions of which one has µ = ∞, giving 1.
Between region I and rotor ferromagnetic core (r = r 1 ) 2.

Continuous boundary condition
By applying Maxwell equations at the interface between different regions [13], the boundary conditions B r,k = B r,k+1 and H ϕ,k = H ϕ,k+1 are obtained, where k and k + 1 indicate adjacent regions having a finite permeability (µ = ∞). Considering regions with the same tangential width, the following continuous boundary conditions are identified: 1.
Between regions I and II (r = r 2 ) B r,I = B r,I I (16) 2.
Between regions III q and IV q and stator ferromagnetic core (r = r 4 ), for q = 1, ..., Q B r,I I I q = B r,IV q (18) •

Combination of Neumann and continuous boundary conditions
Each slot air region III 1, ..., Q has different tangential widths with respect to the radially adjacent airgap region II. This gives rise to a combination of both Neumann and continuous boundary conditions which hold at certain intervals between regions II and III q , for q = 1, . . . , Q: B r,I I I q = B r,I I , for 0 ≤ ϕ jq ≤ α s (20) The slots introduce tangential Neumann boundary conditions on region III q , resulting in H r = 0 at the tangential boundaries ϕ jq = 0 and ϕ jq = α s . Consequently, the cosine component of B r in in region III q vanishes, resulting in B rc = 0 and B ϕs = 0. Therefore, the boundary condition Equation (20) can be expressed as where M denotes the number of modeled harmonics in region I I I q . To simplify Equation (22), a correlation technique is applied [8], resulting in The boundary condition Equation (21) can be expressed as which can be simplified through the correlation technique as follows: where s , c , κ c , κ 0 , ζ c and ζ 0 are the correlation terms (see [8] for the solutions).

Model Verification
The analytical model is verified by comparing the calculated spatial magnetic flux density distribution with that obtained from a 2D nonlinear FEM (Finite Element Method) model. Table 2 presents the motor parameters used for the model verification. As evident from Figure 6, there is an excellent agreement between the calculated airgap (region II) magnetic flux density from the analytical and FEM models. Furthermore, the simulation duration of the analytical model is significantly shorter than that of linear and nonlinear FEM models by factor 20 and 120, respectively, as apparent from Table 3. Table 2. Motor parameters for the model verification.

Parameter Value
No. of PM pole pairs, p 7 No. of slots, Q 12 Inner shaft radius, r 0 2.7 mm Inner ferromagnetic core outer radius, r 1 3.7 mm Inner PM outer radius, r 2 5.7 mm Inner airgap outer radius, r 3 6 mm Slot air outer radius, r 4 6.5 mm Slot winding outer radius, r 5 10.5 mm Outer stator radius, r 6 12.5 mm Axial length, L 18 mm Slot opening, α s 5 • Tooth width, w t 2 mm PM pole-arc to pole-pitch ratio, τ m 1 PM remanence, B r 1.39 PM relative permeability, µ r 1.

Harmonic Modeling of a Magnetic Gear
The considered magnetic gear topology in this paper is depicted in Figure 1, which for modeling purpose is represented by the structure in Figure 7. The modeled material regions (I, II, etc.) of the magnetic gear are described in Table 4.   As in the case of the electrical motor modeling, the magnetic flux density in the magnetic gear is obtained by considering the following boundary conditions:

•
Neumann boundary condition at r = r 1 and r = r 6 • Continuous boundary condition at r = r 2 and r = r 5 • Combination of Neumann and continuous boundary conditions at r = r 3 (between regions II and III 1, ..., Q ) and r = r 4 (between regions IV and III 1, ..., Q ) • Conservation of the magnetic flux around the pole pieces This boundary condition concerns Gauss' law for magnetic field given by By applying the the above to the pole-piece depicted in Figure 8, the following is obtained for which q = 1, . . . , Q − 1, with Q − 1 being the number of independent equations. An extra equation is required to solve all the unknowns, obtained from the application of Ampere's law at the pole-pieces radial boundary highlighted in Figure 8, as follows: where h c is the infinitesimal height of the boundary between region IV and III q . Since B ϕ0,IV = 0, the following is obtained:  Table 5 presents the magnetic gear parameters used for the model verification. As evident from Figure 9, there is an excellent agreement between the calculated airgap magnetic flux density distributions from the analytical and FEM models.

Parameter Value
No. of inner PM pole pairs, p in 2 No. of outer PM pole pairs, p out 5 No. of pole-pieces, Q = p in + p out 7 Transmission ratio, N mg = Q P in 3.5 Inner shaft radius, r 0 2.5 mm Inner ferromagnetic core outer radius, r 1 4.5 mm Inner PM outer radius, r 2 5.5 mm Inner airgap outer radius, r 3 6 mm Pole-piece outer radius, r 4 8.5 mm Outer airgap outer radius, r 5 9 mm Outer PM outer radius, r 6 10 mm Stator outer radius, r 7 12.5 mm Inner PM pole-arc to pole-pitch ratio, τ m,in 1 Outer PM pole-arc to pole-pitch ratio, τ m,out 0.9 Pole-piece arc to pitch ratio, τ Q 0.5 PM remanence, B r 1.39 PM relative permeability, µ r 1.05 Ferromagnetic core material (for FEM) Steel 1010

Definition of Optimization Problem Statements
To optimize the shaft-coupled motor and magnetic gear, optimization problem statements are defined separately for the two considered electromagnetic devices. A general optimization problem statement [14] is given as follows: Minimize the objective function f ( X) (32) subject to inequality constraints g a ( X) ≤ 0 for a = 1, ..., x equality constraints h b ( X) = 0 for b = 1, ..., y and bound constraints X l c ≤ X c ≤ X u c for c = 1, ..., z where X = [X 1 X 2 ... X o ] is the design variable vector, X l c and X u c are the lower and upper bounds of a design variable, respectively. Based on the design requirements, the functions Equations (32)-(34) are derived from the model of the investigated electromagnetic device. In this paper, the chosen optimization algorithm is the interior-point method [15] that can deal with nonlinear constraints and known for its speed and robustness [16].

•
Objective function The considered application requires that the motor and magnetic gear are compact and lightweight. An objective function that handles these requirements is the inverse of the mass torque density where m mot is the motor mass andT mot is the motor peak torque, defined as the highest level of torque that can be sustained for a period of time while the winding temperature rise does not exceed its limit. A thermal model is therefore developed, which is used to estimate the motor winding temperature and to derive the respective inequality constraint function (further described in the next part of this paper subsection).
The motor torque T mot can be calculated using Maxwell stress tensor as follows: where L is the axial length, r in = r 2 +r 3 2 is the radius of the middle of airgap region II in Figure 4, B r,I I and B ϕ,I I are the radial and tangential magnetic flux densities in region II.

Winding temperature
Based on Equation (36), the following constraint on the peak torque is defined where Θ w,10s is the winding temperature that can be sustained for 10 s (given that its initial temperature at time = 0 s is the same as the ambient temperature), while the peak winding copper lossP w is maintained. Note that core (iron) losses also contribute to the temperature increase, although this is neglected in the considered application since the motor maximum speed is 3000 rpm, which is relatively low for the investigated topology and this occurs only in a transient task. Ultimately, constraint function Equation (38) imposes a limit on the winding current density when the motor peak torque is applied. A simple transient thermal equivalent circuit (TEC) in Figure 10 is developed to estimate the temperatures in the winding (Θ w ) and stator (Θ s ) as functions of time, which are obtained by solving where R tw , R y,1 , R y,2 are the radial conductive thermal resistances of the winding, stator yoke (inner and outer part), respectively, and R h is the convective thermal resistance (natural cooling is assumed). C w and C s are the thermal capacitances of the winding and stator, respectively. The calculation of thermal resistances and capacitances are based on the method described in [17]. The winding copper loss P w is calculated as where J is the current density, l coil is the single-turn length of a winding coil loop, S coil is the coil cross-section area and k f is the slot fill-factor. Figure 11 illustrated the temperatures of the motor winding and stator dynamically respond to the changes in copper loss, at ambient temperature of 20 • C. A good agreement between TEC and FEM is apparent.  -

Magnetic flux density in the ferromagnetic cores
The used ferromagnetic core material has a typical saturation point at B = 1.5 T in its B − H characteristic. Thus, to maintain a linear current-torque relation in the motor, the following constraint functions are defined where B t and B s are the magnetic flux density in the stator tooth and stator yoke, respectively, given by where w t is the tooth width, φ ti and φ sj are the flux that goes through the i-th tooth and j-th yoke part (see Figure 12), respectively, estimated as where mod(j + i − 1, Q) in (46) is a modulo operation, giving the remainder after division between j + i − 1 and Q, while the integration limits ϕ i,1 and ϕ i,2 are depicted in Figure 12. -

Torque ripple
A smooth torque characteristic is required in the considered application. For that reason, the ripple in the motor torque shown Figure 13 is constrained by the following function whereT mot is the motor peak torque andT mot,avg is the average value of motor peak torque. •

Equality constraint functions
A series of optimization tasks will be performed on the motor. For a given optimization task, fixed values of motor outer dimensions are assigned. Therefore, the following equality constraint on outer diameter D mot is introduced Meanwhile, the motor axial length L mot is assigned as a design parameter. •

Design variables and bound constraints
The design variable vector X mot consists of motor geometric parameters (see Figure 4) and current density for which each variable adheres to the following lower and upper bound vectors ( X l mot and X u mot , respectively),

Optimization Problem Statement for the Magnetic Gear
• Objective function Similar to the electrical motor optimization, the defined objective function of the magnetic gear is the inverse of mass torque density where m mg is the magnetic gear mass andT mg,out is the peak value of the outer rotor torque, which has a sinusoidal torque-position characteristic. This is shown in Figure 14 together with the torque characteristics of the inner rotor and stator. The outer rotor torque T mg,out is calculated as follows:  Figure 15 shows the variations of torque and discrepancy between analytical and FEM models for different values of the constraints B sat,PP and B sat,SC , belonging to the pole-pieces and stator core, respectively. The constraints B sat,PP = 3 T and B sat,PP = 3 T as calculated by the analytical model are selected based on the previous consideration on torque and model accuracy; thus, the following constraint functions are defined: where B SC is the magnetic flux density in the stator core, estimated through where φ SC is flux that enters the stator core over a range, from ϕ = 0 • to ϕ = α SC as illustrated in Figure 16, for which the angle α SC is determined from the zero crossing of the vector potential spatial distribution A z,V in region V. Meanwhile, B PP is the magnetic flux density in the pole-piece where A z,PP is the vector potential distribution throughout regions III 1 ,...,III Q and the pole pieces, obtained through the the following linear interpolation where A z,I I and A z,IV are the inner and outer airgap vector potential distributions, respectively, while r PP = r 3 +r 4 2 (see Figure 7 for r 3 and r 4 ). As can be seen from Figure 17, there is a good agreement between the analytical and FEM estimated values of B SC and B PP .  •

Equality constraint functions
Following Equation (48), the following constraint on the magnetic gear outer diameter D mg is introduced Additionally, the magnetic gear axial length D mg is fixed as a design parameter.

Modeling
The torque T mmg and speed ω mmg of the shaft-coupled electrical motor and magnetic gear are expressed in terms of motor torque T mot Equation (37), speed ω mot and magnetic gear transmission ratio N mg as follows: where where p in and p out are the numbers of inner and outer PM pole pairs of the magnetic gear, respectively. Meanwhile the combined length of the integrated actuator, L mmg , is given by where L mot and L mg are the axial lengths of the motor and magnetic gear, respectively, and D mmg is the diameter of the integrated actuator, which is a factor to account for the extended length due to the shaft coupling and housing.

Design Requirements
The shaft-coupled electrical motor and magnetic gear will be used as an actuator in a robotic application that has to be compact and lightweight. These qualitative requirements are taken care by optimization objective functions Equations (36) and (52). An important design aspect to note related to the magnetic gear is the choices of inner PM pole pairs p in and transmission ratio N mg , which significantly influence the optimized volume torque density as can be seen in Figure 18. As the numbers of PM pole pairs p in and p out (see Equation (65)) increase, the ferromagnetic cores become less prone to saturation and therefore they can be made thinner to improve the torque within a given magnetic gear volume. However, this improvement has its limitation since increasing the number of magnets implies smaller magnets that deliver less flux, and therefore the volume torque density is reduced again for higher N mg and p in . Based on this and the characteristics in Figure 18, the values p in = 3 and 3.33 ≤ N mg ≤ 9.67 are chosen for the design.
As for the quantitative requirements, the actuator has to deliver a peak torque of 1 Nm and has a maximum speed of 300 rpm. Note that from Equation (64), the required maximum speed can be achieved regardless of the magnetic gear transmission ratio, since the motor maximum speed is 3000 rpm (see Section 3.1). Additionally, there are following requirements related to dimensions

•
The diameter of motor is equal to that of magnetic gear.

•
The maximum axial length of the actuator is twice its diameter.

Response Surface Methodology
To reduce the computational burden associated with the simultaneous electrical motor and magnetic gear optimization, a major reduction of design variables, from 21 Equations ((49) and (60)) to five variables, is performed by expressing the electrical motor and magnetic gear peak torques as functions of their respective outer dimensions and transmission ratio, i.e.,T mot (D mot , L mot ) and T mg,out (D mg , L mg , N mg ). The functions are approximated through the application of response surface methodology from a series of motor and magnetic gear optimization tasks.
Response surface methodology (RSM) is a technique to approximate a functional relationship between a response y and a number of input variables x 1 , x 2 , ..., x k [10] of a given system. The approximated function is in the form of a first-order or the following second-order polynomial model: where β 0 , β i , β ij , β ii are regressions coefficients and is a random error. Equation (67) can be expressed as where X is a n × p matrix containing p variables from Equation (67) obtained from n experiments, β is a vector of regression coefficients and is a vector of random errors from n experiments. The estimated value of the regression coefficients vector,β, is obtained from the least-square method as follows [10]: In the context of this paper, an experiment refers to the set-up of the optimization design variable and parameter values that are chosen based on the Design of Experiment (DOE) method. Several DOE methods for first-and second-order polynomial models are described in [10]. The chosen DOE method in this paper is the popular one for approximating second-order models, i.e., Central Composite Design (CCD) [18,19].
To apply the CCD, a range of input variable values are initially selected. Based on pre-calculations, the following input variable range of values are expected to fulfil the design requirements • 20 mm ≤ D mot ≤ 40 mm, 5 mm ≤ L mot ≤ 15 mm • 20 mm ≤ D mg ≤ 40 mm, 10 mm ≤ L mg ≤ 20 mm, 3.33 ≤ N mg ≤ 9.67.
Next, CCD design matrices are defined for the motor and magnetic gear, as denoted by D CCD,mot and D CCD,mg in Equation (70) , respectively. Each matrix contains the normalized values of the previous input variable values, which are used as design parameters in the optimization problem statements in Sections 3.1 and 3.2. Note that for a given optimization task, the electrical motor assumes the same values of p, Q, B r , µ r as in Table 2 while for the magnetic gear, the same values of B r , µ r in Table 5.
By applying the previously described RSM on a series of optimization tasks arising from the CCD design matrices Equation (70), the following polynomial models are obtained whereT mot andT mg are the optimized peak torques of the motor and magnetic gear, respectively, and the regression coefficients vectorsβ mot andβ mg are given bŷ The motor and magnetic gear peak torque values from the series of optimization tasks are shown in Figure 19. There is a good agreement between the torque values calculated from the RSM and analytical model. Based on this finding, the simple polynomial functions obtained from RSM, i.e., T mot (D mot , L mot ) andT mg,out (D mg , L mg , N mg ), are used for the optimization of the shaft-coupled motor and magnetic gear.

Definition of Optimization Problem Statement
Multiple objectives are considered in the optimization of the shaft-coupled motor and magnetic gear, namely for the minimization of outer diameter D mmg and transmission ratio N mg . A single objective function can be formulated by the scalarization [20] of the two objective functions as follows: where w 1 and w 2 are linear weights (with 0 ≤w 1 ≤ 1 and 0 ≤w 2 ≤ 1) that will be varied in a series of optimization tasks to obtain a Pareto front. Meanwhile, the following constraints are defined based on the design requirements discussed in Section 4.2: Additionally, the following inequality constraint is imposed on the magnetic gear peak torque where the additional 20% in Equation (77) is a margin provided above the magnetic gear peak torque to prevent rotor slip. The design variable vector X mmg is given by for which the lower and upper bounds of the variables are based on the RSM input variable range of values in Section 4.3.

Results
Figure 20 depicts a Pareto front constructed from the objective function Equation (73), which is useful to evaluate design trade-offs between the magnetic gear transmission ratio N mmg and the outer diameter of the shaft-coupled motor and magnetic gear D mmg . As the ratio between the weights in Equation (73), w 1 w 2 increases, the objective function minimization is emphasized on the diameter D mmg , resulting in the decrease of motor and magnetic gear outer dimensions (see Equations (75) and (76)), while the magnetic gear transmission ratio N mmg increases. This illustrates how the magnetic gear ultimately contributes in reducing the motor size. On the other hand, it is also evident from Figure 20 that there is little improvement in motor size reduction for higher values of magnetic gear transmission ratios. This is caused the decrease of magnetic gear torque density as the transmission ratio increases after a certain value, as apparent from Figure 18. Additionally, from a practical point of view, a high magnetic gear transmission ratio implies more difficult manufacturing due to the increasing number of the required permanent magnets (see Equation (65)). Therefore, a specific optimum solution has to be chosen from Figure 20 to account for the trade-off between the transmission ratio and overall size of the shaft-coupled motor and magnetic gear.
The chosen optimum solution is marked in Figure 20, corresponding to the variable values N mg = 5.3, D mot = D mg = 24 mm, L mot = 9.8 mm and L mg = 14.2 mm. The resulting magnetic gear transmission ratio has to match the transmission ratio Equation (65), with p i = 3. Furthermore, a specific value of the number of outer PM pole pairs p o has to be selected such that the resulting cogging torque is low [3,21]. Based on these considerations, N mg = 5.33 (which coincidentally is very close to the selected solution in Figure 20) is chosen as the transmission ratio value. A verification of the optimum solution is performed by separately re-optimizing the motor and magnetic gear using the analytical model as discussed in Section 3, based on the previous selected outer diameters, axial lengths and transmission ratio. The optimization results are shown in Figure 21, the average value of the motor peak torque is 0.19 Nm (resulting inT mmg = 1.01 Nm), while the magnetic gear outer rotor peak torque is 1.29 Nm, which fulfil the design requirements. By comparing the outer dimensions and volume of the shaft-coupled motor and magnetic gear, with those of an optimized electrical motor (having the same design requirements) in Table 6, there is an apparent benefit in terms of size reduction if the former actuator type is used.  Optimized peak motor torqueT mot and magnetic gear outer rotor torque T mg,out characteristics.

Conclusions
An actuator based on a shaft-coupled electrical motor and magnetic gear is described in this paper. High-accuracy analytical electromagnetic models of the electrical motor and magnetic gear are developed and verified. Based on these models, design optimization objective and constraint functions are formulated in detail to account for electromagnetic and performance aspects such as magnetic saturation of the ferromagnetic cores, torque ripple/cogging and torque density.
For the simultaneous optimization of the electrical motor and magnetic gear, response surface methodology (RSM) is employed to reduce the number of design variables originating from the two electromagnetic devices, by representing their optimized torques as polynomial functions of their respective outer dimensions and transmission ratio (for the magnetic gear). A factor four of design variables reduction is achieved, and, as the polynomial functions are used to model the actuator, optimizations can be performed rapidly without sacrificing the accuracy of their results. From the constructed Pareto front, it is evident that the actuator size can be reduced by increasing the magnetic gear transmission ratio up to a certain level, after which the magnetic gear torque density limitation becomes more apparent and no further actuator size reduction can be achieved.