Derive an expression for the drive torque of the stepper motor shown in Fig. 5.48. Assume that the rotor is a radially polarized multipole cylinder with a magnetization

Question

M = ± [latex] M_{s} \hat{r},[/latex] (5.280)

where ± takes into account the alternating polarity of adjacent poles. Further assume that each stator is operating below saturation in a linear region of its B-H curve with μ >> [latex] μ_{0}[/latex].

Accepted Answer

The torque is derived in two steps. First, we obtain a field solution in the interior of an energized stator with no rotor present. Second, we reduce the rotor to an equivalent current distribution and then compute the torque on it using a Lorentz force approach with the stator field as an external field.

Stator field: Let L and [latex]R_{s}[/latex] denote the length and inner radius of the stator, respectively. Each stator has [latex]N_{pole}[/latex] teeth, which equals the number of poles on the rotor. The stator is enclosed by an n-turn coil and a flux return ring. A cross section of the interior of a four-pole stator is shown in Fig. 5.49a. Let [latex]Φ_{t}[/latex] denote the angular span of each tooth and [latex]Φ_{g}[/latex] denote the angular expanse of the gap between each tooth. These angles are related by the following equation,

[latex]Φ_{g} = \frac{2π}{N_{pole}} -Φ_{t} [/latex]. (5.281)

We seek a field solution for the interior of the stator. To obtain an analytical solution, we assume that the stator is infinitely long. This assumption is not as radical as it might first seem because the stator tends to confine the flux and minimize leakage etc. In the interior of the stator,

∇ × H = 0.

Thus, we can represent H in terms of a scalar potential, H = - ∇ [latex]φ_{m }[/latex] (Section 3.4). The field problem reduces to a boundary value problem (BVP) in cylindrical coordinates (Section 3.6.2). Specifically, we solve

[latex]∇^{2} φ_{m }[/latex] = 0. (5.282)

The general solution to Eq. (5.282) that is well behaved at r = 0 is of the form

[latex]φ_{m }(r, Φ) = \sum\limits_{p=1}^{∞} a_{p} r^{p} sin(pΦ) + b_{p} r^{p} cos(pΦ).[/latex] (5.283)

We choose a coordinate system in which Φ = 0 coincides with the center of a stator tooth (Fig. 5.51a). In this coordinate system [latex]φ_{m }[/latex](r, Φ) must be an even function of Φ. Thus, we discard the sin(nΦ) terms in Eq. (5.283) and obtain

[latex]φ_{m }(r, Φ) = \sum\limits_{p=1}^{∞}b_{p} r^{p} cos(pΦ).[/latex] (5.284)

We determine the coefficients [latex]b_{p}[/latex] by matching the boundary conditions at the interior surface of the stator. We assume that the potential [latex]φ_{m}(R_{s}, Φ)[/latex] along this surface is constant across each tooth and alternates from [latex]φ_{max}[/latex] to -[latex]φ_{max}[/latex] on successive teeth. As the gap between neighboring teeth is small [latex](Φ_{g}<<Φ_{t})[/latex], we assume that the potential is linear between the teeth and this implies a profile for [latex]φ_{m }(R_{s}, Φ)[/latex] as shown in Fig. 5.51b. To solve the BVP we need a Fourier series representation for [latex]φ_{m }(R_{s}, Φ)[/latex]. First, we derive an expression for [latex]φ_{max}[/latex] in terms of the motor’s physical parameters. Let g denote the gap between adjacent stator teeth, and let [latex]H_{g}[/latex] denote the field in the gap region (Figs. 5.49b and 5.50). We assume that the potential is linear between adjacent teeth. Therefore,

[latex]H_{g} = - \frac{1}{R_{s}} \frac{Δφ_{m }}{ΔΦ} [/latex] (5.285)

= - [latex]\frac{2φ_{max}}{g},[/latex] (5.286)

where ΔΦ = [latex]Φ_{g}[/latex] and g = [latex]R_{s}Φ_{g}[/latex]. We apply Eq. (3.138) from Section 3.5 to the dotted path in Fig. 5.50 and obtain

Eq. (3.138): [latex] \oint_{c}H ⋅dl = I_{tot},[/latex]

[latex] H_{g} = \frac{ni}{g}.[/latex] (5.287)

Therefore,

[latex] φ_{max} =- \frac{ni}{2}.[/latex] (5.288)

Now that we know [latex]φ_{max}[/latex], we can represent the potential along the stator [latex]φ_{m }(R_{s}, Φ)[/latex] in terms of the Fourier series

[latex]φ_{m }(R_{s}, Φ) = \sum\limits_{k=1,3,5,...}^{∞}A_{k}cos(α_{k}Φ),[/latex] (5.289)

where

[latex]A_{k} = -\frac{8}{πk^{2} N_{pole}} \frac{φ_{max}}{Φ_{g}} [cos(kπ) - 1 ] cos(α_{k} \frac{Φ_{t}}{2}) [/latex]

= [latex]\frac{4}{πk^{2} N_{pole}} \frac{ni}{Φ_{g}} [cos(kπ) - 1 ] cos(α_{k} \frac{Φ_{t}}{2}), [/latex] (5.290)

and [latex]α_{k} = kN_{pole}/2.[/latex] We evaluate Eq. (5.284) at r = [latex]R_{s}[/latex], compare this with Eq. (5.289), and obtain

[latex]b_{p} = A_{k} R_{s}^{-p} \qquad (p = \frac{kN_{pole}}{2})[/latex].

Thus, the potential in the interior of the stator is

[latex]φ_{m }(r, Φ) = \sum\limits_{k=1,3,5,...}^{∞}A_{k} (\frac{r}{R_{s}})^{α_{k}}cos(α_{k}Φ).[/latex] (5.291)

The field components are obtained from B = - [latex]μ_{0}∇φ_{m}.[/latex] For example,

[latex] B_{r}(r, Φ) = -μ_{0}\sum\limits_{k=1,3,5,...}^{∞}\frac{α_{k}}{R_{s}^{α_{k}}} A_{k} r^{(α_{k} -1)} cos(α_{k}Φ)[/latex]. (5.292)

Drive torque: To compute the torque we represent the rotor as a distribution of equivalent currents. Recall that the rotor is radially polarized with

M = ± [latex]M_{s}\hat{r} [/latex].

Let [latex]R_{1}[/latex] and [latex]R_{2}[/latex] denote its inner and outer radius, respectively. Assume that it is initially positioned with one of its north poles centered about 0° (i.e., each pole is initially centered with respect to a stator tooth). The equivalent current densities are given by Eq. (3.95). We find that the volume current density is zero, [latex]J_{m}[/latex] = ∇ × M = 0, and that the surface current density is

Eq. (3.95):Equivalent Currents
[latex] J_{m} = ∇ × M (A/m^{2})[/latex] (volume current density)
[latex]j_{m} = M × \hat{n}[/latex] (A/m) (surface current density).

[latex]j_{m} = M × \hat{n}[/latex]

= ± [latex]M_{s}\hat{z} [/latex].

These surface currents are located at the edges of each sector at angular positions

[latex]Φ_{n} = (n + \frac{1}{2}) \frac{2π}{N_{pole}}[/latex] (n = 0, 1, 2, . . . , [latex]N_{pole}[/latex] - 1).

It is important to note that there are two such currents at these positions, one along the surface of each of the neighboring sectors. If the rotor is rotated by an angle θ, the positions of the surface currents change to

[latex]Φ_{n}(θ) = θ +(n+\frac{1}{2}) \frac{2π}{N_{pole}}[/latex] (n = 0, 1, 2, . . . , [latex]N_{pole}[/latex] - 1).

The surface current densities at these locations are

[latex]j_{m}(n) = \left\{\begin{matrix} M_{s}\hat{z} \qquad (n = 0, 2, 4, . . .) \ - M_{s}\hat{z} \qquad (n = 1, 3, 5, . . .).\end{matrix} ight.[/latex] (5.293)

The torque is given by Eq. (3.100) from Section 3.3, that is,

Eq. (3.100): T = [latex] \int_{v} r ×(J_{m} ×B_{ext}) dv + \oint_{s} r × (j_{m} × B_{ext}) ds [/latex],

T(θ) = [latex]2L \sum\limits_{n=0}^{N_{pole} -1} \int_{R_{1}}^{R_{2}} r× [j_{m}(n) × B_{stator} (r , Φ_{n}(Φ))] dr[/latex]. (5.294)

Factor 2 takes into account the contribution of each surface of the neighboring sectors. From Eq. (5.294) we find that the component of torque about the z-axis is

[latex]T_{z}(θ) = 2LM_{s} \sum\limits_{n=0}^{N_{pole} -1} (-1)^{n} \int_{R_{1}}^{R_{2}} B_{r}(r , Φ_{n}(Φ)) [/latex] r dr. (5.295)

Finally, substituting Eq. (5.292) into Eq. (5.295) gives

[latex]T_{z}(θ) = 2 μ_{0}LM_{s} \sum\limits_{n=0}^{N_{pole} -1} \sum\limits_{k =1,3,5,...}^{∞} (-1)^{n} \ × A_{k} ( \frac{α_{k}}{α_{k} +1} )\frac{(R_{2}^{(α_{k}+1)} - R_{1}^{(α_{k}+1)}}{R_{s}^{α_{k}}} cos[α_{k}Φ_{n}(θ)], [/latex] (5.296)

where L and [latex]R_{s}[/latex] are the length and interior radius of the stator and

[latex]A_{k} = \frac{4}{πk^{2} N_{pole}} \frac{ni}{Φ_{g}} [cos(kπ) - 1 ] cos(α_{k} \frac{Φ_{t}}{2}) [/latex].

Here,

[latex]α_{k} = \frac{kN_{pole}}{2}[/latex],

and

[latex] Φ_{n}(Φ)= θ + (n + \frac{1}{2}) \frac{2π}{N_{pole}}[/latex].

Calculations: We demonstrate Eq. (5.296) with some sample calculations. The following parameters are used in the analysis:

[latex]M_{s} = 4.3 × 10^{5}[/latex] A/m
[latex]N_{pole} [/latex] = 10
n = 100 turns
i = 0.25 A
[latex]R_{1} [/latex] = 6.0 mm
[latex]R_{2} [/latex] = 8.0 mm
[latex]R_{s} [/latex] = 8.1 mm
L = 2.0 mm.

The angular span of each stator tooth is taken to be [latex]Φ_{t}[/latex] = 27° (75% of the span of a rotor pole). The angular span of the gap between neighboring teeth is [latex]Φ_{g}[/latex] = 9° ( [latex]Φ_{g}[/latex] = 2π/ [latex]N_{pole} - Φ_{t}[/latex]). There are two phases to consider, which are labeled 1 and 2, respectively. Recall that the stator teeth of the two phases are offset from one another by an angular span equal to one half of the pole pitch, which in this case is 18°. The drive torque [latex]T_{z}[/latex] (θ) for each phase is shown in Fig. 5.52. This plot shows the torque on the rotor as it rotates the angular span of one pole (0 ≤ θ ≤ 36°). When θ = 0, the rotor is positioned as in Fig. 5.49c and d for phase 1 and phase 2, respectively. In this initial position, phase 1 exerts almost maximum torque, and phase two exerts zero torque. The peak torque for phase 1 is 1500 μN · m, which occurs at θ = 12°. When the rotor is rotated to θ = 18° (half a pole pitch) as shown in Fig. 5.49d, the torque of phase 1 is zero and the torque of phase 2 is near a maximum. For 18 ≤ θ ≤ 36°, the torque of phase 1 reverses direction as shown.

Question 5.14.1: Derive an expression for the drive torque of the stepper mot...

Related Answered Questions

Derive the equations of motion for the actuator shown in Fig. 5.53. ...

Verified Answer:

Develop a three-dimensional field solution for the gap region of an axial-field motor (Fig. 5.42). Assume that the magnet has a second quadrant demagnetization curve ...

Verified Answer:

Develop a two-dimensional model for the field in the gap region of an axial-field motor (Fig. 5.37). Assume the magnet has a linear second quadrant demagnetization curve of the form ...

Verified Answer:

Set up the equations of motion for the rotary bias actuator shown in Fig. 5.28. The magnet has a radius a and length h, with a << h. The coil has n turns and has a length lc. ...

Verified Answer:

Develop a model for the design of the axial-field resonant actuator shown in Fig. 5.23. ...

Verified Answer:

Determine the equations of motion for the axial-field actuator shown in Fig. 5.16. Assume that the coil has n turns. ...

Verified Answer:

Consider the cylindrical actuator shown in Fig. 5.14. Assume that the magnetic circuit is linear with μ >> μ0, and that there is no flux leakage. Set up the equations of motion for the actuator [2]. ...

Verified Answer:

Set up the equations of motion for the actuator shown in Fig. 5.15. Assume that all magnetic materials are linear with μ >> μ0. The nonmagnetic sleeves of thickness s provide low friction support for the moving plunger [4]. ...

Verified Answer:

Consider the magnetic actuator circuit shown in Fig. 5.13a. Assume that the magnetic circuit operates in a linear region of its B-H curve with μ >> μo, and that there is no flux leakage. Determine the force on the plunger [2]. ...

Verified Answer:

Set up the equations of motion for the actuator circuit of Fig. 5.11 [1]. ...

Verified Answer: