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

B = μH + $μ_{0} M_{s},$ (5.261)

where

μ = $\frac{B_{r}}{H_{c}}$ . (5.262)

Question

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

B = μH + [latex] μ_{0} M_{s},[/latex] (5.261)

where

μ = [latex]\frac{B_{r}}{H_{c}} [/latex]. (5.262)

Accepted Answer

We reduce the three-dimensional motor geometry to an approximate two-dimensional geometry [21]. This is done by introducing a cylindrical cutting plane at the mean radius of the magnets. This cylinder is then imagined to be unrolled into a two-dimensional surface of infinite extent into and out of the page, as illustrated in Fig. 5.40a. This two-dimensional geometry can be reduced further by exploiting the symmetry that results from the repeating magnetic structure. Specifically, the field is symmetric about the vertical center lines of each pole. The field solution for the reduced geometry is similar to that of the periodic geometry of Example 3.6.2. Therefore, we follow closely the solution procedure presented there. As depicted in Fig. 5.40b, there are three regions to consider. In this figure, l is the arc length from the center of one pole to the next at the mean radius, and x is the distance along the ‘‘straightened’’ arc. Therefore, if [latex]R_{1}[/latex] and [latex]R_{2}[/latex] are the inner and outer radii of the magnet, respectively, and [latex]N_{pole}[/latex] = the number of poles, then

l = [latex] \frac{π(R_{1} + R_{2})}{N_{pole}}[/latex].

The variable x is related to the angular measure Φ along the arc as follows:

Φ = [latex]\frac{x}{l} \frac{2π}{N_{pole}} [/latex]. (5.263)

The magnetization for the three regions is

[latex]M_{s} = \left\{\begin{matrix} M_{s}\hat{y} \qquad (Region 1)\ -M_{s}\hat{y} \qquad (Region 2) \qquad (5.264)\ 0 \qquad (Region 3). \end{matrix} ight. [/latex]

Because [latex]M_{s}[/latex] is constant in each region, the equation for the magnetic scalar

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

for all three regions. The three regions under consideration are collectively bounded by four lines that form a rectangular area. The vertical sides of the rectangle, which are defined by the y-axis and the line x = l, represent lines of symmetry for the field, along which the normal component [latex]H_{n}[/latex] is zero. Moreover, the boundaries that form the top and bottom of the rectangle are bordered by soft-magnetic materials (flux return plates). For the sake of simplicity, we assume that the flux plates have infinite permeability (μ = ∞). Consequently, the tangential component of the field vanishes along these boundaries. When we impose these boundary conditions we find that the general solutions for the magnet and gap regions are potential reduces to

[latex]φ_{mag}(x,y) = \sum\limits_{n=1}^{∞} C_{mag,n} sinh (\frac{nπ(y-h)}{l}) cos(\frac{nπx}{l}) [/latex] (5.266)

and

[latex]φ_{gap}(x,y) = \sum\limits_{n=1}^{∞} C_{gap,n} sinh (\frac{nπy}{l}) cos(\frac{nπx}{l}) [/latex], (5.267)

respectively, where [latex]C_{mag,n}[/latex] and [latex]C_{gap,n}[/latex] are constants that need to be determined. At the interface between the magnet and the gap, the tangential component of H must be continuous:

[latex]H_{mag,t}(x,g) = H_{gap,t}(x,g)[/latex] (0 < x < l). (5.268)

This condition, along with the fact that the functions sin(nπx/l) are orthogonal on (0, l), implies that

[latex]C_{mag,n} = C_{gap,n} \frac{sinh(nπg/l)}{sinh(nπ(g-h)/l)}.[/latex] (5.269)

In addition, the normal component of B must also be continuous

[latex]B_{mag,n}(x,g) = B_{gap,n}(x,g)[/latex] (0 < x < l). (5.270)

As

[latex]B_{mag,n} = \left\{\begin{matrix}μ H_{mag,n} + μ_{0} M_{s} \qquad (0 < x < l/2)\ μ H_{mag,n} - μ_{0} M_{s} \qquad (l/2 < x < l), \end{matrix} ight.[/latex] (5.271)

then Eq. (5.270) reduces to

[latex] \sum\limits_{n=1}^{∞} \frac{nπ}{l} C_{gap,n} K (n, h, g, l, μ) cos (\frac{nπx}{l}) = \left\{\begin{matrix} \frac{μ_{0}}{μ} M_{s} \qquad (0 < x < l/2)\ - \frac{μ_{0}}{μ} M_{s} \qquad (l/2 < x < l)\end{matrix} ight.[/latex]

where

K (n, h, g, l, μ) = [latex]-\frac{μ_{0}}{μ} cosh(\frac{nπg}{l}) + sinh(\frac{nπg}{l}) coth(\frac{nπ}{l} (g - h))[/latex]. (5.272)

Again, by exploiting the orthogonality of the functions cos(nπx/l) on (0, l), it follows that

[latex]C_{gap,n} = \frac{4lμ_{0}M_{s}}{μ} \frac{(-1)^{(n-1)/2}}{(nπ)^{2} K (n, h, g, l, μ)} [/latex] (n = 1, 3, 5, . . .). (5.273)

Substituting Eq. (5.273) into Eq. (5.267) gives the complete solution for the gap region

[latex]φ_{gap}(x,y) = \frac{4lμ_{0}M_{s}}{μ}\sum\limits_{n = 1, 3, 5, . . .}^{∞} \frac{(-1)^{(n-1)/2}}{(nπ)^{2} K (n, h, g, l, μ)} sinh(\frac{nπy}{l}) cos (\frac{nπx}{l}). [/latex]

Therefore, the vertical or axial component of the field in this region is

[latex]B_{gap,y}(x,y) = \frac{4μ_{0}^{2}M_{s}}{μ}\sum\limits_{n = 1, 3, 5, . . .}^{∞} \frac{(-1)^{(n-1)/2}}{(nπ) K (n, h, g, l, μ)}[/latex]

× [latex] cosh(\frac{nπy}{l}) cos (\frac{nπx}{l})[/latex], (5.274)

where x equals the horizontal distance along a circumferential arc at the mean radius of the magnet, y equals the vertical distance above the lower flux plate, g is the gap height, [latex]t_{m}[/latex] is the thickness of the magnet, h = g + [latex]t_{m}[/latex], l = π [latex](R_{1} + R_{2})/N_{pole}[/latex], and [latex]R_{1}[/latex] and [latex]R_{2}[/latex] are the inner and outer radii of the magnet.

Calculations: We apply Eq. (5.274) to a motor with the following parameters:

[latex]M_{s} = 4.0 × 10^{5} [/latex] A/m
g = 6.0 mm
[latex]t_{m}[/latex] = 4.0 mm
[latex]R_{1} = \frac{40}{π} [/latex] mm
[latex]R_{2} = \frac{120}{π} [/latex] mm
l = 20 mm
[latex]N_{pole}[/latex] = 8. (5.275)

Field values are computed with μ = 1.0 [latex]μ_{0}[/latex], 1.5 [latex]μ_{0}[/latex], and 2.0 [latex]μ_{0}[/latex], and with the vertical position in the gap set to y = 2 mm. Ten values are computed for each value of μ. Specifically, field values are computed at Φ = 0, 5, 10, . . . , 45° (i.e., from the center of one pole to that of its neighbor). These data are shown in Fig. 5.41. Notice that the field decreases with increasing μ.

Question 5.13.1: Develop a two-dimensional model for the field in the gap reg...

Related Answered Questions

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

Verified Answer:

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 ...

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:

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: