Spontaneous Speed Reversals in Stepper Motors

Similar documents
Stepping Motors. Chapter 11 L E L F L D

Step Motor Modeling. Step Motor Modeling K. Craig 1

Motor Info on the WWW Motorola Motors DC motor» /MOTORDCTUT.

Automatic Control Systems. -Lecture Note 15-

Overview of motors and motion control

FEEDBACK CONTROL SYSTEMS

Mathematical Modeling and Dynamic Simulation of a Class of Drive Systems with Permanent Magnet Synchronous Motors

MODELING AND HIGH-PERFORMANCE CONTROL OF ELECTRIC MACHINES

Positioning Servo Design Example

Dynamic Modeling of Surface Mounted Permanent Synchronous Motor for Servo motor application

Implementation of a Self Tuning Regulator and Optimal Design of Advanced Control Strategies with Full order Observer for Stepper Motors

UNIVERSITY OF BOLTON SCHOOL OF ENGINEERING BSC (HONS) MECHATRONICS TOP-UP SEMESTER 1 EXAMINATION 2017/2018 ADVANCED MECHATRONIC SYSTEMS

Robust Controller Design for Speed Control of an Indirect Field Oriented Induction Machine Drive

Mechatronics Engineering. Li Wen

CHAPTER 1 Basic Concepts of Control System. CHAPTER 6 Hydraulic Control System

Power System Stability and Control. Dr. B. Kalyan Kumar, Department of Electrical Engineering, Indian Institute of Technology Madras, Chennai, India

ECE 5670/6670 Lab 8. Torque Curves of Induction Motors. Objectives

Laboratory Exercise 1 DC servo

International Journal of Advance Engineering and Research Development SIMULATION OF FIELD ORIENTED CONTROL OF PERMANENT MAGNET SYNCHRONOUS MOTOR

E11 Lecture 13: Motors. Professor Lape Fall 2010

Dynamics of the synchronous machine

Chapter 4. Synchronous Generators. Basic Topology

GATE 2008 Electrical Engineering

Mathematical Modeling and Dynamic Simulation of DC Motors using MATLAB/Simulink Environment

CHAPTER 2 DYNAMIC STABILITY MODEL OF THE POWER SYSTEM

2.004 Dynamics and Control II Spring 2008

Synchronous Machines

Dynamic Performance Analysis of Permanent Magnet Hybrid Stepper Motor by Transfer Function Model for Different Design Topologies

Control of Wind Turbine Generators. James Cale Guest Lecturer EE 566, Fall Semester 2014 Colorado State University

Lecture 8: Sensorless Synchronous Motor Drives

R10 JNTUWORLD B 1 M 1 K 2 M 2. f(t) Figure 1

Analysis and Design of Control Systems in the Time Domain

Lesson 17: Synchronous Machines

Equal Pitch and Unequal Pitch:

Quanser NI-ELVIS Trainer (QNET) Series: QNET Experiment #02: DC Motor Position Control. DC Motor Control Trainer (DCMCT) Student Manual

Feedback Control Systems

Index. Index. More information. in this web service Cambridge University Press

Mechatronics Modeling and Analysis of Dynamic Systems Case-Study Exercise

INC 341 Feedback Control Systems: Lecture 3 Transfer Function of Dynamic Systems II

ECE 585 Power System Stability

Lecture 1: Induction Motor

Behaviour of synchronous machine during a short-circuit (a simple example of electromagnetic transients)

MODELING WITH CURRENT DYNAMICS AND VIBRATION CONTROL OF TWO PHASE HYBRID STEPPING MOTOR IN INTERMITTENT DRIVE

Electric Machines I Three Phase Induction Motor. Dr. Firas Obeidat

6) Motors and Encoders

3 d Calculate the product of the motor constant and the pole flux KΦ in this operating point. 2 e Calculate the torque.

Manufacturing Equipment Control

A Direct Torque Controlled Induction Motor with Variable Hysteresis Band

LIAPUNOV S STABILITY THEORY-BASED MODEL REFERENCE ADAPTIVE CONTROL FOR DC MOTOR

MAE143a: Signals & Systems (& Control) Final Exam (2011) solutions

C(s) R(s) 1 C(s) C(s) C(s) = s - T. Ts + 1 = 1 s - 1. s + (1 T) Taking the inverse Laplace transform of Equation (5 2), we obtain

Lecture 7: Synchronous Motor Drives

INDUCTION MOTOR MODEL AND PARAMETERS

Appendix A: Exercise Problems on Classical Feedback Control Theory (Chaps. 1 and 2)

Torques 1.0 Two torques We have written the swing equation where speed is in rad/sec as:

Sensorless DTC-SVM of Induction Motor by Applying Two Neural Controllers

Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science Electric Machines

Definition Application of electrical machines Electromagnetism: review Analogies between electric and magnetic circuits Faraday s Law Electromagnetic

ADMISSION TEST INDUSTRIAL AUTOMATION ENGINEERING

Frequency Response Techniques

Performance of Feedback Control Systems

R a) Compare open loop and closed loop control systems. b) Clearly bring out, from basics, Force-current and Force-Voltage analogies.

CPPM Mahine: A Synchronous Permanent Magnet Machine with Field Weakening

Robust Speed Controller Design for Permanent Magnet Synchronous Motor Drives Based on Sliding Mode Control

Reliability of Bulk Power Systems (cont d)

Equivalent Circuits with Multiple Damper Windings (e.g. Round rotor Machines)

PERFORMANCE ANALYSIS OF DIRECT TORQUE CONTROL OF 3-PHASE INDUCTION MOTOR

Texas A & M University Department of Mechanical Engineering MEEN 364 Dynamic Systems and Controls Dr. Alexander G. Parlos

SCHOOL OF ELECTRICAL, MECHANICAL AND MECHATRONIC SYSTEMS. Transient Stability LECTURE NOTES SPRING SEMESTER, 2008

Lab 3: Quanser Hardware and Proportional Control

Loss analysis of a 1 MW class HTS synchronous motor

JRE SCHOOL OF Engineering

Application Note #3413

Current Reduction in Stepping Motor Applications using an Adaptive PI controller based on Linearized Dynamics

Time Response Analysis (Part II)

Accurate Joule Loss Estimation for Rotating Machines: An Engineering Approach

Electrical Machines and Energy Systems: Operating Principles (Part 2) SYED A Rizvi

Comparative Analysis of Speed Control of Induction Motor by DTC over Scalar Control Technique

POWER SYSTEM STABILITY

Lecture 6: Control Problems and Solutions. CS 344R: Robotics Benjamin Kuipers

Analytical Calculation of Air Gap Magnetic Field Distribution in Vernier Motor

Introduction to Synchronous. Machines. Kevin Gaughan

Synchronous Machines

EE 410/510: Electromechanical Systems Chapter 4

The Effects of Machine Components on Bifurcation and Chaos as Applied to Multimachine System

Finite Element Analysis of Hybrid Excitation Axial Flux Machine for Electric Cars

NEPTUNE -code: KAUVG11ONC Prerequisites:... Knowledge description:

DEVELOPMENT OF DIRECT TORQUE CONTROL MODELWITH USING SVI FOR THREE PHASE INDUCTION MOTOR

(a) Find the transfer function of the amplifier. Ans.: G(s) =

DYNAMIC ANALYSIS OF DRIVE MECHANISM WITH FUNCTIONAL MODEL

A Novel Adaptive Estimation of Stator and Rotor Resistance for Induction Motor Drives

DETERMINING STATIC TORQUE CHARACTERISTICS OF PERMANENT MAGNET STEP MOTOR

H-INFINITY CONTROLLER DESIGN FOR A DC MOTOR MODEL WITH UNCERTAIN PARAMETERS

ENHANCEMENT MAXIMUM POWER POINT TRACKING OF PV SYSTEMS USING DIFFERENT ALGORITHMS

Chapter 14. Reluctance Drives: Stepper-Motor and Switched- Reluctance Drives

EE 742 Chapter 3: Power System in the Steady State. Y. Baghzouz

Chapter 3 AUTOMATIC VOLTAGE CONTROL

(b) A unity feedback system is characterized by the transfer function. Design a suitable compensator to meet the following specifications:

ROEVER COLLEGE OF ENGINEERING & TECHNOLOGY ELAMBALUR, PERAMBALUR DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING ELECTRICAL MACHINES I

Transcription:

Spontaneous Speed Reversals in Stepper Motors Marc Bodson University of Utah Electrical & Computer Engineering 50 S Central Campus Dr Rm 3280 Salt Lake City, UT 84112, U.S.A. Jeffrey S. Sato & Stephen R. Silver L-3 Communications 640N 2200W M/S F1J12 P.O. Box 16850 Salt Lake City, UT 84116, U.S.A. Abstract Experimental data shows that permanent magnet stepper motors can spontaneously reverse their direction of rotation when controlled in full step, open-loop mode. The paper shows that the reversal of speed can be explained, and that the reverse speed is equal to three times the forward speed. A linear approximation of the motor dynamics predicts that reverse running will arise when the undershoot of the single-step response exceeds 50% (or when the damping factor is less than 0.11). The observations of the paper suggest a revision of the conventional explanation that resonance problems occur because of high positive velocity at the stepping time. Instead, the data and the analysis of the paper suggest that resonance problems, including stalling and speed reversal, occur due to the growth of the undershoot in the response to repeated steps. Keywords: stepper motor, stepping motor, synchronous motor, resonance, stall, reverse motion, nonlinear dynamics. 1. Introduction The paper describes a phenomenon in which the speed of a permanent magnet (PM) stepper motor reverses unexpectedly. Fig. 1 shows data collected on a stepper motor at

Position (rad) 10 8 6 4 2 0-2 -4-6 0 1 2 3 4 5 6 Figure 1: Speed reversal of a stepper motor the University of Utah. The stepper motor is a two-phase permanent magnet motor with 7.5 step size (i.e., 12 pole pairs). The motor was commanded in open-loop, full stepping mode, with a voltage-source-amplifier and a sinusoidal reference profile. The figure shows several instances where the speed changes direction with respect to the command. A similar phenomenon was observed at L-3 Communications with a hybrid stepper motor (1.8 step size) under current command. Stepper motors controlled in open-loop are known to be susceptible to resonances at step rates equal to the frequency of resonance and to its subharmonics. A speed reversal phenomenon was reported in [5] for a variable-reluctance stepper motor. Fig. 1 indicates that this phenomenon can also be observed with PM stepper motors, although such motors are less susceptible to resonance problems. In this paper, we show how the phenomenon can be explained, and how its onset can predicted. The analysis suggests a revision to the common interpretation of low-speed resonance in stepper motors. 2. Explanation of speed reversal The possibility of speed reversal is counter-intuitive, since the stator currents create a 2

o o o 0 90 o 180 270 (a) Forward o 0 o -270 o -540 =-180 o -810 o =-90 o (b) Reverse Figure 2: Principle of speed reversal at three times the forward speed rotating magnetic field in a given direction. The existence of the phenomenon is tied to the stepping nature of the rotating field. Consider Fig. 2, where a permanent magnet motor with one pole pair is shown. Arrows show the directions of the rotor and of the stator magnetic fields. Assuming a two-phase motor, part (a) of the figure shows 4 steps of the forward motion of the rotor under full stepping commands. Part (b) shows how a reverse motion with the same equilibrium positions is theoretically possible. Note that the resulting reverse speed is 3 times the forward speed. The 3x reversal speed is consistent with data obtained at the University of Utah and with observations made at L-3 Communications on a separate system. The data reported in [5] (also reproduced in [4]) for a three-stack variable-reluctance motor indicates a 2x speed reversal. Note that our explanation of the speed reversal, appropriately adjusted, predicts a 2x reversal speed for a 3-phase motor. The same explanation also predicts the possibility of 5x, 9x,..., forward speeds, and 7x, 11x,... reverse speeds for a two-phase PM motor. The phenomenon is analogous to aliasing in sampled-data systems. Sustained operation at speeds other than the nominal forward 3

(a) (b) (c) (d) (e) (f) Figure 3: Intermediate positions of reverse motion speed, however, requires that the dynamics of the system permit transition from one equilibrium to the next along the pattern shown in Fig. 2. It is not immediately obvious how the rotor can move from the 0 to the 270 position, since the torque acts in the counterclockwise direction. Our investigation of the phenomenon has led us to conclude that, for reverse motion to be possible, the rotor must lag behind the stable equilibrium position by approximately 90 at the end of a step in the sequence. Fig. 3 explains how reverse motion can then be sustained. Stage (a) shows the position of the rotor at the end of a step, with the rotor lagging by slightly more than 90 behind the direction of the stator magnetic field. Torque is positive, but when the next step starts at stage (b), the angle becomes greater than 180, yielding a small negative torque. The motor will move in the negative direction. Stage(c)showsanintermediatepositionalongthatpath,wherethenegativetorquehas maximum magnitude. Note that it may take some time to reach position (c), since the torque is initially small. At (d), the torque is again zero, but the rotor has acquired a sufficient velocity that it overshoots the equilibrium. The torque becomes positive and the motor slows down while reaching position (e). At that time, the next step in the sequence is reached and the cycle repeats, with (f) equivalent to (b). In this manner, reverse motion with a negative average torque is possible. Note that the reverse motion is also possible if the initial angle is less than 90, provided that the motor has a small negative velocity, since the torque is small after switching. 4

It remains to be shown that the sequence of Fig. 3 can be initiated, i.e., that a 90 lag can be obtained. This is the subject of the next two sections. 3. Stepper motor model and characteristics Consider a two-phase permanent magnet (PM) stepper motor under current command. Denote the currents flowing in the stator windings i A and i B (in A). Let θ be the angle of the rotor (in rad), so that θ =0whenthePMflux in winding A is maximized. Assuming linear magnetics and a sinusoidal PM flux in the windings, the torque τ e (in N.m) produced by the motor is given by τ e = Ki A sin(n p θ)+ki B cos(n p θ) (1) where K is the torque constant (in N.m/A) and n p is the number of pole pairs (or rotor teeths) of the motor. In practice, the torque is not exactly sinusoidal, and reluctance effects (such as the detent torque) affect the overall torque. However, the approximation is adequate for our analysis. The torque can also be written as τ e = KI sin (θ e n p θ) (2) where I = q i 2 A + i2 B, θ e = ](i A,i B ) (3) I is the magnitude of the 2-dimensional current vector with components i A and i B,and θ e is the angle of the vector with respect to the x-axis. Open-loop control of position consists in applying currents such that the angle θ e is varied by discrete increments. The full step two-phases-on method corresponds to the sequence θ e i A i B 45 I 0 I 0 135 I 0 I 0 225 I 0 I 0 315 I 0 I 0 (4) 5

The sequence produces equilibrium positions at 45,135,..., separated by 90 electrical degrees. The dynamic response of a motor can be described by J d2 θ dθ = B dt2 dt + KI sin (θ e n p θ) (5) where J is the inertia of the motor, and B is a constant associated with viscous friction. Static friction and load effects (except for the load inertia that may be added to the inertia of the motor) have been neglected. (5) is a nonlinear differential equation. However, if the rotor stays close to the equilibrium position defined by the current vector (i.e., θ e n p θ is small), an approximate linear differential equation can be obtained The poles of the system are given by the roots of J d2 θ dθ = B dt2 dt + KI (θ e n p θ) (6) s 2 + B J s + KIn p J =0 (7) Because viscous friction is small in stepper motors, the roots are generally poorly damped. If θ e suddenly changes, the response of the motor is oscillatory. The resonant frequency (Hz) of the system can then be estimated from f r = 1 r KInp 2π J (8) Fig. 4 shows the two-phases-on step response of a stepper motor at the University of Utah. The response was obtained under voltage command, which generally provides more damping than current command. Yet, significant overshoot and multiple oscillations are observed. 4. Resonance and initiation of reverse motion 6

0.55 0.5 Position (rad) 0.45 0.4 0.35 0.3 0.38 0.4 0.42 0.44 0.46 0.48 0.5 Figure 4: Motor step response Consider a linear time-invariant system with poles at s = a ± jb and with unity DC gain H(s) = Y (s) X(s) = a 2 + b 2 s 2 +2as + a 2 + b 2 (9) The response of the system to a unit step input (x(t) =1)isgivenby y(t) =1 e at cos(bt) e at sin(bt) (10) For a/b small, the response is oscillatory. Peaking occurs for t 1 = π b, y(t 1)=1+e aπ b t 2 = 2π b, y(t 2)=1 e 2aπ b (11) and so on. The principle of resonance is shown in Fig. 5 for a =0.1, b = 1. If a second step is applied at time t 2 =2π/b, the response of the system for t>t 2 (shown as a solid curve) is the sum of the response due to the first step and the response due to the second step (both shown as dashed curves). Note that, due to the timing of the second step, the 7

3.5 3 2.5 Ouptut 2 1.5 1 0.5 0 0 2 4 6 8 10 12 14 Figure 5: Response of second-order system to double step at the frequency of resonance overshoot is considerably magnified after the second step. The positive peak becomes Asymptotically, the percent overshoot (PO) is given by y( 3π b )=2+e aπ b + e 3aπ b (12) PO(%) 100 = e aπ b + e 3aπ b +...= e aπ b 1 e 2aπ b (13) Fig. 5 shows that the undershoot is also magnified, and becomes asymptotically PU(%) 100 = e 2aπ b + e 4aπ b +...= e 2aπ b 1 e 2aπ b Note that the asymptotic undershoot reaches 100% if the first undershoot is 50%, that is, (14) if e 2aπ b =0.5, or a b =0.11 (15) The condition is numerically similar in terms of the damping factor ζ = a/b p 1+(a/b) 2 =0.11 (16) 8

In the context of full stepping of a stepper motor, a 100% undershoot means that the motor lags the equilibrium position by one full step. If the undershoot occurs when the next step is applied, the rotor will lag the next equilibrium position by two full steps (or 180 electrical degrees). Therefore, the condition for a 100% undershoot defines a limit above which loss of synchronism will occur. Reverse motion, as shown in Fig. 3 can begin. The nonlinear differential equation (5) predicts the growth of undershoot and the torque reversal. Fig. 6 shows the result of a simulation with multiple step inputs and with a step rate (in steps/s) equal to the resonant frequency (in Hz). One can clearly observe the build-up of the undershoot until the 180 position is reached and reverse motion starts. The progression is not exactly as predicted by the linear theory, because the oscillations take the system well into the nonlinear region (90 in the sin function). In practice, reverse motion may also be triggered, or delayed, by perturbations such as load or gear train oscillations. However, the linear approximation is useful to predict the likelihood of reverse motion. Once reverse motion is started, it may continue and show up as a backward-running motor. Often, however, the motor will switch between periods of forward and reverse motions. The resulting behavior will appear as stalling. Nevertheless, stalling and steady reverse motion are two resonance problems having the same origin. Our explanation of resonance in stepper motors is different from the common interpretation. Indeed, Fig. 7 shows a plot representative of diagrams found in [1], [3], [5]. The interpretation is that a step rate equal to the resonant frequency yields switching at a time where the rotor has high forward velocity. As the overshoot is amplified, a situation is reached where the torque changes sign and the rotor escapes in the forward direction. This interpretation does not match the data that we have collected or, for that matter, the data found in [5]. In fact, escape in the forward direction is unlikely to occur. In Fig. 7, the position at which the torque changes signatthetimeofstepcommand3isstep 9

Position (rad) 1.4 1.2 1 0.8 0.6 0.4 0.2 0 0 0.05 0.1 0.15 0.2 0.25 Figure 6: Simulation of stepper motor response to multiple steps position 5. Escape in the forward direction requires an overshoot of 200%. Condition (13) indicates such an overshoot is possible, although requiring a ratio a/b = 0.079 that is smaller than the condition for a 100% undershoot. Note that, for poorly damped systems, the first overshoot is close in magnitude to the first undershoot. Because escape in the reverse direction occurs for a 100% undershoot, instead of the 200% overshoot required in the forward direction, escape in the reverse direction will generally occur before it may occur in the forward direction. 5. Experimental results Data was collected on the motor whose step responsewasshownonfig. 4(Airpaxmotor #C82711-M1, 17VDC dual coil, 7.5 step). Note that the firstundershootisgreaterthan 50%, suggesting a high risk of speed reversal. The frequency of resonance was estimated to be f r =41.7Hz. For 12 pole pairs, a step rate equal to the frequency of resonance corresponds to a speed of rotation equal to ω RP M = 60f r 4np = 52rpm (17) 10

4 Step Positions 3 2 1 1 2 3 Step Commands Figure 7: Typical interpretation of resonance in stepper motors In experiments, resonances were found to occur in the range of 43 to 55rpm, approximately. Theresponsesinthemiddleoftherangewerecharacterizedbystalling. Fig. 8shows a typical example at 49rpm. One can see that stalling actually results from frequent switching between forward speed and reverse speed at 3 times the forward speed. Fig. 9 shows the response at 51rpm, which shows longer periods of forward and reverse motions. At the end of the plot, there is even a period where alternations of forward (1x) and reverse (3x) motions yield a reverse motion with an average speed close to the forward speed. As in [5], speed reversals were observed in the frequency bands bordering normal behavior, while stalling was observed in the middle of the resonance band. Also as in [5], reverse motion sometimes persisted for significant periods of time. In experiments at L-3 Communications, we have observed sustained reverse motion at 3x the forward speed for several revolutions of the motor. Fig. 10 shows speed reversals at 25rpm, which corresponds to the first subharmonic of the resonant frequency. The same theory shows that, for sufficiently low damping, build- 11

3 2.5 2 Position (rad) 1.5 1 0.5 0 0 0.5 1 1.5 2 2.5 Figure 8: Stepper motor response at 49rpm 5 4 3 Position (rad) 2 1 0 1 0 0.5 1 1.5 2 2.5 Figure 9: Stepper motor response at 51rpm 12

9 8 7 6 Position (rad) 5 4 3 2 1 0 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 Figure 10: Stepper motor response at 25rpm up of the second overshoot can result in loss of synchronism. The (approximate) condition is that the second overshoot must exceed 50%. Away from the resonant frequency and its first subharmonic, the motor was found to track the reference commands well. Speeds up to 250rpm could be reached in open-loop mode. A close look at the responses of the motor shows that reverse motion follows the stages of Fig. 3. Fig. 11 shows a detail of Fig. 1, with the smooth curve showing the response of the motor. The stable equilibrium positions corresponding to the commands are shown by the stepwise curves (separated by 2π/n p = π/6). The figure confirms that the reverse motion occurs at three times the forward speed, and that the motor briefly settles at the unstable equilibrium positions at the beginning of a new step. Interestingly, the onset of this reverse motion can be predicted by a relatively simple differential equation, although it strongly depends on the nonlinear nature of the differential equation. 6. Avoiding resonances and speed reversals Methods for avoiding low-speed resonances are well-known and include ([1], [2], [4]): 13

10 Position (rad) 9.5 9 8.5 8 4.05 4.1 4.15 4.2 4.25 4.3 4.35 Figure 11: Detail of reversal during multiple steps increasing mechanical damping, increasing damping through electrical means, microstepping, closed-loop control. Ofren, microstepping will be the easiest approach. Even a change from two-phaseson to one-phase-on can eliminate the problem. For the motor used in the experiments, the one-phase-on mode of stepping (with equilibrium positions at 0,90,..., insteadof 45,135,...) gives the step response of Fig. 12, which is considerably more damped. The undershoot is less than 50% and resonance problems have not been observed in that mode. The two-phases-on is often preferred because of its higher torque and because of the simplicity of the power electronic circuit. However, lower damping results not only fromthehighertorque,butalsofromthedetent torque that destabilizes the equilibrium directions of the two-phases-on mode. 14

0.6 0.55 Position (rad) 0.5 0.45 0.4 0.35 0.5 0.52 0.54 0.56 0.58 0.6 Figure 12: Step response in one-phase-on mode 7. Conclusions The paper showed that stepper motors controlled in open-loop may reverse direction unexpectedly. Although a similar phenomenon had been observed for variable-reluctance motors, the paper went further by providing an explanation for how this behavior could be initiated and sustained. Under a linear approximation, it was found that a condition for the initiation of a speed reversal is that the undershoot in the response exceeds 50% (or a damping factor less than 0.11). The resulting reversal of the torque is the source of resonance problems including both stall and speed reversal. While the usual interpretation suggests that switching at a time of high forward velocity is the underlying cause of the reversal, the analysis and data of this paper shows that the problem originates instead in the growing undershoot. 8. References [1] P. Acarnley, Stepping Motors: A Guide to Theory and Practice, 4th edition, Institution of Electrical Engineers, United Kingdom, 2002. 15

[2] M. Bodson, J. Chiasson, R. Novotnak, & R. Rekowski, High-Performance Nonlinear Feedback Control of a Permanent Magnet Stepper Motor, IEEE Trans. on Control Systems Technology, vol.1,no.1,1993,pp.5-14. [3] C.B. Gray, Electrical Machines and Drive Systems, Longman, United Kingdom, 1989. [4] T. Kenjo & A. Sugawara, Stepping Motors and Their Microprocessor Controls, 2nd edition, Oxford Science Publications, United Kingdom, 1994. [5] P.J. Lawrenson & I.E. Kingham, Resonance Effects in Stepping Motors, Proc. of the Institution of Electrical Engineers, vol. 124, no. 5, 1977, pp. 445-448. 16