16 Linear DC PM Brushless Motors

16.1 Introduction

The linear dc homopolar-PM motor has been developed from the voice-coil speaker-microphone linear electromagnetic device; however, in most cases, it uses a flat (rather than tubular) topology and is applied for progressive (rather than oscillatory) linear motion; also, the travel may reach 100 mm or more.

As it is based on Lorenz force, to change direction of motion the polarity of the current is reversed [1–3]. The main asset is the simplicity of the control with full force starting from any position when supplied from a 4 quadrant dc–dc converter, for a rather long (100 mm or even more), but limited, travel. Short coil mover with long PM or long coil stator with short PM-mover combinations may be used, depending on the application.

The copper losses per N of thrust are not small as the coil is in the air, to avoid any cogging (zero current) force and thus allow for precision control.

Recording pens, optical disk drives, focusing systems in digital video cameras, xy chart recorders, etc., are typical applications of linear dc PM brushless motors [4–9].

Linear dc PM brush motors have also been proposed [1], but today these are of mere historical interest; consequently, only brushless topologies are considered.

Finally, as tubular configurations will be investigated for linear oscillatory machines in a subsequent chapter, only flat configurations will be treated in this chapter.

Design issues of iron-core-backed or air-core-backed PM topologies are studied in some detail; then iron losses and geometrical optimization design are illustrated; finally, sliding mode precision control is investigated, as applied in focusing a digital video camera.

16.2 Topology Aspects

Two basic topologies are shown in Figure 16.1a and b [3,8].

The single-coil mover in Figure 16.1a—intended for travels above 10 mm—allows for a larger force density than the dual-coil mover with PM in air (Figure 16.1b), but the latter has a larger force/weight, so crucial in some applications, when the mover has a small weight (1–2 g). For larger force per volume, single-coil dual-core configurations may be used (Figure 16.2).

As the thin magnetic fixture may be replaced even with a nonmagnetic rugged one, the inductance of the coil is rather small and thus the force variation with mover position is very small, because the core magnetic saturation under load is not large, provided the PM flux does produce a reasonable flux density in the core. Also, the iron losses are to be smaller as they are produced mainly by the local variation of the core magnetic flux density due to armature current.

Out of the four coil sides, this time two are active (only one side is active in Figure 16.1); this justifies the claims of larger force density and lower copper losses per N of force. Lateral nonmagnetic plates (Figure 16.2b) may provide better PM stiffness over PM length up to 0.1–0.15 m travel or even more.

To double the thrust for the single coil, the PM thickness has to be doubled (roughly), but the core weight is only slightly larger (20% or so) while the coil weight is also only slightly increased (20% or so). The configuration in Figure 16.2 shows further potential for improvements in linear dc PM brushless topologies.

Images

FIGURE 16.1 Basic topologies: (a) with one-coil mover and iron cores, (b) with PM in air and dual-coil mover.

16.3 Principle and Analytical Modeling

Let us consider here the analytical model of the configuration in Figure 16.1a, illustrated in more detail in Figure 16.3.

With N turns per coil, wc, active coil length and BPM, PM flux density, considered uniform here along x (direction of motion), the Lorenz force on the coil Fx is

Fx=Kfi;Kf=BPMNwc(16.1)

Images

FIGURE 16.2 Single-coil dual-core configuration.

Images

FIGURE 16.3 Geometry and flux paths of linear dc PM mover. (After Wakiwaka, H. et al., IEEE Trans., MAG-32(5), 5073, 1996.)

At center point, the PM flux and the armature (coil mmf) flux lines and flux density distribution are given in Figure 16.4.

As visible in Figure 16.4, the armature reaction increases the core flux on the right side of mover and decreases it on the left side, much like in a standard rotary dc brush motor (transversal armature reaction).

Now, with the mover at extreme position, the coil (armature) flux stays the same (within the core) while the PM fringing (leakage) at its ends leads to an even smaller force. So Kf is expected to vary with mover position due to magnetic saturation (Figure 16.4c) because the field of the coil mmf closes mainly within the magnetic core.

This reality would become evident with 3D FEM (even with 2D FEM) analysis, but an approximate analysis, valid for preliminary design, would be in order here.

So let us consider the center position when the flux produced in the core by the coil is

ϕ=Bawctc;Ba=Niμ(Bc)2(Ls+Lc+ls+tp+g+tc)(16.2)

Images

FIGURE 16.4 Flux lines and flux densities at (a) center position, (b) extreme position, and (c) force versus position at given current. (After Boldea, I. and Nasar, S.A., Linear Motion Electromagnetic Devices, Taylor & Francis, New York, 2001.)

On the other hand, the PM flux in the core right corner ϕc is

ϕc=Bcmaxtcwc=Bawctc+BPM(Ls+Lc)2wc(16.3)

Bcmax in the core is limited by the core material B(H) curve and may be assessed by combining (16.2) and (16.3):

ϕcmax=wctcBcmax=[μ(Bc)Ni2(Ls+Lc+ls+lp+g+tc)+BPM(Ls+Lc)2tc]wctc(16.4)

Equation 16.4 may be solved iteratively, if first

BPMBrtptp+g11+Kfringe(16.5)

The fringe factor Kfringe depends on mover position to some extent (for start Kfringe = 0.1–0.2).

Equation 16.4 suggests that, for given Bcmax, for more force, more core area (wctc) and more Ls + Lc (total length) are required. Conversely, if the length of travel Ls increases, by specifications, then, for more force, the core area (wctc) should also increase.

The force per watt of copper losses is also to be considered. The ideal coil thickness (g = hc + 2gw) has to be smaller or equal to PM thickness to secure the highest force density:

hcg=NinLcjconKfill(16.6)

Kfill is the copper fill factor

jcon is the design current density

The coil length lcoil is

lcoil2(wc+g)(16.7)

and thus the coil resistance Rc is

Rc=1σcolcoilinjconN(16.8)

σco is the copper electrical conductivity.

So the design copper losses pcon can be written as

pcon=2σcon(wc+NinLcjconKfill)jconNin(16.9)

Finally, the design force/watt of copper losses is

fx=Fxnpcon=BPMσco2(1+NinwcLcjconKfill)jcon(16.10)

Now the PM flux in the coil ΨPM becomes

ψPM=NBPMwcx11+Kfringe(x)(16.11)

With Kfringe ≈ 0.1–0.2 ≈ const, x = 0 in the center position.

So the emf is

E=dψPMdxdxdtNBPMwc1+Kfringe(x)u;uisthespeed(16.12)

Finally, the coil inductance is

Lc(1+Kl)μ(i)N2wctc2(Ls+Lc+ls+tp+g+tc)(16.13)

The core permeability μ(i) is an average one; in an analytical model, it may be a weighted average of three values in three key points of the core, while the actual magnetization curve B(H) of the core has to be considered. The factor Kl refers to the leakage inductance, which, in any case, is less than 0.2 even for, say, a soft magnetic composite core; not so for an air-core device.

Note on core losses: An analytical model of core losses may be obtained by refining the previously mentioned analytical model via the magnetic equivalent circuit model, where quite a few zones of uniform flux density variations have to be identified.

Then, by knowing the flux density variation in time, for the given linear motion profiles, the core loss may be calculated, by standard expressions, in each region; then they are added together.

Alternatively, 3D FEM may be used in magnetostatic mode, and then with the flux density/time variations in each finite element, the total core losses are calculated. For a tubular configuration, 2D FEM may be used directly, considering eddy current losses and even hysteresis losses [5].

16.4 Geometrical Optimization Design By FEM

For the configuration in Figure 16.1a, detailed, for the geometrical variables, in Figure 16.5, the optimization variables are obtained as L1, L2, L3, L4, and L5. Two inequality constraints in mm [7] are visible in Figure 16.5:

23688L3L720.5†††2L3+L536(16.14)

A few equality constraints (specifications) are also given [7]:

  • Current density jcon = 24.5 A/mm2

  • Coil volume = 280 mm3

  • Mechanical gap G1 = 0.4 mm (on one side of the coil)

  • G2 = 0.5 mm—gap at the end of travel

With a given initial variable vector [L1, L2, L3, L4, L5]T, a 2D FEM (at least), and an optimization method, to change Li in order to get maximum thrust, the optimization design is completed.

Images

FIGURE 16.5 Definition of geometrical design variables. (After Boldea, I. and Nasar, S.A., Linear Motion Electromagnetic Devices, Taylor & Francis, New York, 2001.)

We may augment the objective function by adding, say, temperature constraints, based on a simplified thermal model, etc.

In our case, with maximum force as the objective function, but also with given coil volume and current density, the copper losses are fixed. Therefore, in fact we are seeking here maximum thrust/watt of copper losses.

In Ref. [7], the Rosenbrock optimization method was used for the scope. The initial and final variable vectors are shown in Table 16.1 and the results are given in Table 16.2 [7].

The nonuniform distribution of force versus position is proven by FEM (Figure 16.6) but is not severe (Table 16.1).

During the 39 iterations (for a 0.01 mm convergence criterion), the variable vector was changed by a given step of 10% of initial values and the thrust was calculated by 2D FEM. If the thrust increased with respect to the previous step, the variables were changed by

ΔLj(k+1)=αLj(k);α>1.5otherwise,ΔLj(k+1)=βLj(k);β<0.75(16.15)

TABLE 16.1
Dimensions of Initial and Final Shapes

Images

TABLE 16.2
Turns of Coil

Thrust (N)

Position of Coil (mm)

Initial Shape

Optimal Shape

Ratio Increase (%)

−8.5

0.82

0.99

20

0

0.88

1.13

29

8.5

0.86

1.05

22

Images

FIGURE 16.6 Initial force versus position. (After Takahashi, S., Optimal design of linear dc motor using FEM, Record of LDIA-1995, Nagasaki, Japan, pp. 331–334, 1995.)

A 20% improvement in the thrust was obtained to reach 0.23 N/W [7], which is still modest, but the fixed current density was 24.5 A/mm2 and the coil volume was 280 mm3.

16.5 Air-Core Configuration Design Aspects

For an optical disk drive application, the PM-in-air configuration (Figure 16.1b) was adopted [8]. A virtual path magnetic circuit model of the machine proved to compare acceptably with 3D FEM results (Figure 16.7), except for small airgap values. Even for this air-core machine, the PM flux density in the coil area is about 0.3 T.

The experimental results (Figure 16.8) for the mini linear dc PM motor of Table 16.3 [8] show a wide frequency response band, suitable for the application.

Images

FIGURE 16.7 Force versus airgap length between the two PMs. (After Park, J.H. et al., IEEE Trans., MAG-39(5), 3337, 2003.)

Images

FIGURE 16.8 Frequency response of the prototype. (After Park, J.H. et al., IEEE Trans., MAG-39(5), 3337, 2003.)

TABLE 16.3
Discretization Data

Number of elements

3.456

Number of nodes

1.741

Number of iterations of Rosenbrock’s method

39

CPU time (s)

201

16.6 Design for Given Dynamics Specifications

For the iron-core linear dc brushless PM motor, in a small power application (a digital video camera focuser), a more complicated optimal design should include not only the energy conversion ratio and the battery consumption but also the rising time (travel time for given travel length).

If the space is restricted to volume (lx × ly × lz), Figure 16.9, we may end up with a few main variables [6]:

γ=lmlw=PMthicknesscoil†thickness;N=round(lwϕ)round(lpϕ)(16.16)

lp is the coil length

ϕ is the coil conductor diameter

Images

FIGURE 16.9 Linear dc PM focuser for digital video camera. (After Yu, H.C. and Liu, T.S., IEEE Trans., MAG-43(11), 4048, 2007.)

The three performance indexes are then

  • Rising time tr:

    dmax=0trudt(16.17)

    u is the speed

    dmax is the travel length

  • Battery consumption (in J) E0:

    E0(γ,ϕ)=0tri(t)v(t)dt(16.18)

  • Energy efficiency:

    ηe(γ,ϕ)=E00tri2(t)R(γ,ϕ)dtE0(16.19)

All variables in (16.17) through (16.19) are functions of γ and ϕ; to use them in the optimization design, we have to solve for current i and speed u from voltage and motion equations:

didt=v(t)i(t)R(γ,ϕ)Ku(γ,ϕ)u(t)L(γ,ϕ)(16.20)

dudt=Ku(γ,ϕ)i(t)FloadBu(t)M(16.21)

where

L is the coil inductance

M is the mover mass

Ku is the emf (and force) constant (see Equation 16.11):

Ku(γ,ϕ)=dΨPMdx=Fxi=Kf(16.22)

If the force Fx is calculated by 2D FEM, then Ku(γ, ϕ) may be determined from (16.22); also, the inductance L may be calculated from FEM, while the resistance Rc is

Rc=4σcoNlcoil(γ)πϕ2(16.23)

With the given voltage, R, L, Kf, Ku, Bm, dmax, and M and Fload, the three performance indexes may be calculated after solving for U and i in (16.20 and 16.21).

For M = 2 g, v = 3 V, imax = 0.03 A, dmax = 5.2 mm, Bm = 0.005 N s/m, and Fload = 0.05 gw (gram weight), to compromise the three performance indexes finally [6], γ = 2.5, ϕ = 0.07 mm, for ηe = 4.9% (energy efficiency!), E0 = 1.3 mJ, and rising time tr = 44 ms [6]. This rising time is considered much better (6 times smaller) than that obtained with commercial rotary-stepper-plus transmission actuator. To fully characterize the design, Rc = 32.8 Ω, L = 1.2 mH, and Ku = Kf = 42.3 gw/A (measured) [6].

Though the energy efficiency ηe seems low, we have to consider it in comparison with existing solutions for such a small device.

16.7 Close-Loop Position Control for a Digital Video Camera Focuser

We are extending here the digital video camera focuser (Figure 16.9) analysis by dealing with the close-loop position control (open loop dynamics were investigated in the previous paragraph).

At this time, the voltage and motion equations parameters are as given previously.

The movable part is composed of the lens holder, the coil, and a linear magnetic strip mounted on a side tube. The stator includes the PMs, the yoke, a steel plate, and two motion-guiding pads. The position encoder includes the magnetoresistive (MR) sensor and a linear magnetic strip with 0.88 mm pole pitch; the final position estimation precision is around 3 μm. In addition, speed estimation is calculated from position feedback.

Two position control methods are used for comparison:

  • PID control

  • Sliding mode control

The PID control contains two cascaded loops, one for position and one for speed (Figure 16.10).

For sliding mode control, a sliding mode functional S(e), e = d* − d (position error) should be defined [9]:

S(e)=e¨+c2e˙+c3e(16.24)

If the sliding mode poles are p1,2 = −α ± iβ, then

c2=2α,c3=(α+β)2(16.25)

Images

FIGURE 16.10 PID position control.

To be stable, the system has to stay close to the sliding mode surface S(e) = 0, but first to reach the sliding mode condition (surface)

SS<σ(S);σ>0(16.26)

To verify the stability and define a control law, the Lyapunov criterion is used, after a Lyapunov (energy) function V is chosen [9]:

V=S2(16.27)

But

V˙=dVdt=2SS˙<2σ(s)=2σV<0(16.28)

As V > 0, V < 0 and thus the control will be asymmetrically stable if (16.26) is satisfied.

The control law may be defined, based on (16.24), as

U=Qsat(s)L31c2e¨L31c3e˙(16.29)

With Q > 0, the maximum boundary of system states and disturbance is

Sat(S)=sign(S);if|S|>εS/ε;if|S|ε(16.30)

Images

FIGURE 16.11 Sliding mode position control.

Images

FIGURE 16.12 Step position response: (a) in PID control and (b) by SM control. (After Yu, H.C. and Liu, T.S., IEEE Trans., MAG-43(11), 4048, 2007.)

The Sat function was used, instead of Sign function, to provide for low chattering response around the target. The close-loop control scheme is illustrated in Figure 16.11. The parameter Q may be adjusted based on minimization of error energy function En(e)=e2dt.

It has to be noticed that

  • There is no cascaded control (as for PID control).

  • Quite a few position error time derivatives are required, but, as position is measured, the two derivatives may be obtained by adequate observers; also if an inertia accelerometer were made available, it would smooth notably the ė and ë observers (as done in MAGLEVs).

Sample results for the case in point, for PID and SM control, are given in Figure 16.12.

Both focusing time (from 120 to 60 ms) and steady-state precision (from 21 to 7 μm) are reduced for SM in comparison with PID control, as expected (SM control is known to provide fast and robust response with small overshooting).

16.8 Summary

  • The linear dc homopolar brushless PM motor has been developed from the voice-coil speaker microphone, but it develops control for progressive linear motion rather than oscillatory motion.

  • It is based on the Lorenz force and is essentially a single-phase linear synchronous brush-less PM motor with limited travel.

  • The practical travel length goes from a few mm to 100 mm or more.

  • Flat configurations with air core and iron core are typical for the linear dc brushless PM motor, while tubular ones are used in most speaker-microphone applications.

  • The Lorenz force sign changes with current polarity.

  • For limited travel, no brushes are required and the single (dual)-coil mover is fed from a flexible power cable.

  • The PM-mover configuration is feasible, but then the copper Joule loss in the stator coils is larger.

  • Iron-core configurations show higher force/volume while air-core configurations show slightly better force/weight but, also, higher frequency response band (due to lower electric time constant and inductance).

  • For iron-core flat configurations, analytical field model for circuit model parameters, and performance including magnetic saturation, produce satisfactory results for preliminary design purposes.

  • Secondary effects such as force reduction at travel ends may be ascertained properly with FEM analysis.

  • Force density and force/copper losses are two conflicting performance indexes that have to be reconciled in a proper design.

  • Core losses may be calculated by analytical or numerical field modes and are important when the duty cycle of back and forth motion is large.

  • For air-core topologies, virtual magnetic circuit models have been proven acceptable (validated by 3D FEM), except for a small distance between PMs (when two or three of them are provided); notable computer time is saved this way for optimization design.

  • Direct FEM geometrical optimization design has been successfully proven on linear dc brushless PM motors with 0.8 N force for 17 mm travel length at 0.23 N/W of copper losses as the objective function.

  • Direct FEM geometrical optimization may be exercised with some sophisticated performance indexes such as rising time, travel time per travel length, total energy consumption per rising time, and energy efficiency for the same. For a 2 g coil mover, a rising time of 44 ms, at total energy consumption of 1.3 mJ and energy efficiency of 4.9%, was obtained in view of focusing a digital video camera.

  • For the same application with given circuit parameters R = 32.8 O, L = 1.2 mH, Ku = Kf = 42.3 gw/A, and M = 1.8 g, both cascaded PID position and speed control and sliding mode (SM) position control have been investigated for a magnetostriction position sensor (5 µm in precision); SM control reduced the settling time for a 5 mm travel at 60 ms for a steady-state positioning precision of 7 µm, adequate for design video camera focusing; the existing rotary microstepper needs 400 ms for the same job.

  • Due to ruggedness, simplicity, and its rather long travel length, the linear dc brushless PM motor will find more and more applications when high-frequency band motion (in the hundreds of Hz range) is needed; the hot spots in R&D will perhaps be geometrical FEM, optimization design, and robust high-precision fast positioning control (For more, see Refs. [10,11]).

References

1. S.A. Nasar and I. Boldea, Linear Motion Electric Machines, John Wiley, New York, 1976.

2. A. Basak, Permanent Magnet Linear Dc Motors, Oxford University Press, 1996.

3. I. Boldea and S.A. Nasar, Linear Motion Electromagnetic Devices, Taylor & Francis, New York, 2001, Chapter 7.

4. H. Wakiwaka, H. Yajima, S. Senoh, and H. Yamada, Simplified thrust limit equation of linear Dc motor, IEEE Trans, MAG-32(5), 1996, 5073–5075.

5. M. Hippner, H. Yamada, and T. Mizuno, Iron losses in linear dc motor, IEEE Trans, MAG-38(5), 1999, 3715–3717.

6. H.-C. Yu, T.-Y. Lee, S.-J. Wang, M.-L. Lai, J.-J. Ju, D.-R. Huang, and S.-K. Lin, Design of a voice coil motor used in the focusing system of a digital video camera, IEEE Trans, MAG-41(10), 2005, 3979–3981.

7. S. Takahashi, Optimal design of linear dc motor using FEM, Record of LDIA-1995, Nagasaki, Japan, pp. 331–334.

8. J.H. Park, Y.S. Baek, and Y.P. Park, Design and analysis of mini-linear actuator for optical disk drive, IEEE Trans, MAG-39(5), 2003, 3337–3339.

9. H.C. Yu and T.S. Liu, Output feedback sliding mode control for a linear focusing actuator in digital video cameras, IEEE Trans, MAG-43(11), 2007, 4048–4050.

10. P. Fang, F. Ding, Q. Li, and Y. Li, High response electromagnetic actuator with twisting axis for gravure systems, IEEE Trans, MAG-45(1), 2009, 172–175.

11. C.-S. Liu, P.-D. Lin, S.-S. Ke, Y.-H. Chang, and J.-B. Horng, Design and characterization of miniature auto focusing voice coil motor actuator for cell-phone camera applications, IEEE Trans, MAG-45(1), 2009, 155–159.

..................Content has been hidden....................

You can't read the all page of ebook, please click here login for view all page.
Reset
18.118.227.69