Mathematical modeling of stepping motor and vibration torque mechanism research on its different operations

Hao Zhao1 , Hao Feng2

1Jiaxing University, Jiaxing Zhejiang 314001, China

2Hangzhou Dianzi University, Hangzhou Zhejiang 310018, China

2Corresponding author

Journal of Vibroengineering, Vol. 15, Issue 3, 2013, p. 1102-1111.
Received 5 June 2013; accepted 3 September 2013; published 30 September 2013

Copyright © 2013 Vibroengineering This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Abstract.

Vibration torque are existence obviously during operation of stepping motor, it is a periodic structure vibration problem. In this paper, the motion model of motor is deduced and is compared with the motion mathematic model of pendulum; then actual vibration torque testing of stepper motor is executed through a self made vibration torque sensor, including single-step operation, low frequency continuous operation, low frequency lost step oscillation and continuous operation, the experimental results show that the stepping motor vibration torque of single three shot operation is more bigger than the six shot, and obvious reply oscillation existence in two operations; sharp vibration torque is generated on low frequency lost step oscillation; the vibration torque of single step operation is greater than the continuous operation, and vibration torque is decreasing with the frequency increasing; last, the generation of these experimental phenomena are analyzed, and the relationship between the stepper motor vibration torque peak value and frequency of the normal continuous operation are found, it can provide a certain reference for controller or control algorithm designation of stepper motor.

Keywords: stepper motor, vibration torque, turbulence detection, mechanism analysis.

1. Introduction

Stepper motor is a kind of micro motor which convert electric pulse into angular displacement or line displacement, its angular displacement is proportional to the number of pulses and has no relation to voltage or out condition in the load capacity scope, and it is used mainly as the implementation components in the open-loop control system, such as printers, NC machine tools, robots and other places where needs to be located accurately [1-4]. But the stepper motor exist obvious oscillation in operation process, it is manifested as the motor speed fluctuation (resonance or jitter), even cause displacement and vibration on extremely cases, and it is the main factors for step motor speed stability, because the torque angle characteristic of stepper motor is nearly sine function, and the magnetic field of the stator turning as jumping.

In recent years, the torque fluctuation or speed fluctuation of stepper motor is studied by some experts and scholars. According to the starting out-of-step and mechanical impact of stroke termination effectively on linear stepper motor, a control system based on DSP2407 was proposed by HongMin Zheng [5]. A new type of stepper motor controller is designed by Bo Qu [6], the performance of the microprocessor STM32F103RBT6 and the driving principle of stepper motor driver chip L6208 are analysed. To solve the problems existed in traditional stepper speed control algorithm, a novel algorithm based on space control instead of time control was introduced by Baoshan You [7]. The measurement of the stepper motor rotation stability by optical encoder and its evaluation criterion were presented by Pan Jin Yu [8], it is based on the analysis of effect for rotation performance of stepping motor in low-speed scan. Based on repetitive control of time-varying periodic signals, a new method of restraining toque ripple for 5-phase hybrid stepping motor system is proposed [9], and the principle of torque ripple is analyzed. By analyzing two-phase hybrid type stepping motor mechanical structure, a variety of factors impacting on the starting process of stepping motor were deduced theoretically, and a set of measurement system of stepping motor performance based on the photoelectric sensor was constructed [10].

In this paper, a self-made vibration torque sensor is used for measuring the stepper motor torque fluctuation [11]. The sensor is able to detect instantaneous torque fluctuation relative to the grating and photoelectric encoder method in the reference. Last, the experimental results are analyzed and studied, it can provide a certain reference for controller or algorithm designation of stepping motor.

2. Stepper motor motion model

The reluctance torque of reactive step motor produced by different reluctance of teeth, and the static torque is refer to a electromagnetic torque, which generated on conditions that a windings going through DC current, the expression of stepping motor static electromagnetic torque T can be induced through the principle of energy balance relationship, which exist in energy conversion between mechanical and electrical. When A phase windings of stepper motor going through DC current ia:

(1)
U a = r a i a + d ( L a i a ) d t .

In Eq. (1), Ua is control voltage applied in A phase windings; ra, La are resistance and inductance of A phase windings respectively; the input energy of windings in dt is:

(2)
U a i a d t = i a 2 r a d t + i a d L a i a .

Inductance is magnetic chain number produced by unit current, then La=ψia , so we can receive:

(3)
U a i a d t = i a 2 r a d t + i a d ψ .

It can indicated that the energy of A phase windings one part is converted into heat energy ia2radt, it is consumption of the windings resistance ra, the other part is converted into electromagnetic energy stored in the magnetic field of A phase windings, it is the source of torque to do foreign power. When the response electromagnetic torque T drive the rotor turned the deflection angle dθ, the stepping motor output mechanical energy Tdθ which is converted from the storage magnetic field energy, so dWm=Tdθ according to the energy balance relationship, then:

(4)
T = d W m d θ .

According to the magnetic energy expression:

(5)
W m = 1 2 L I 2 = 1 2 Φ F .

By the definition of inductance and the following relation:

(6)
ψ = N φ = N 2 N I R m = 2 N 2 I Λ .

In Eq. (6), ψ is magnetic chain of A phase windings; 2NI is magnetomotive force of magnetic circuit; Rm is the magnetic resistance of A phase windings magnetic circuit; Λ is the magnetic conductivity of A phase windings magnetic circuit, it can be considered as the air gap magnetic conductivity when the magnetic circuit unsaturated; so the variation of storage magnetic energy is:

(7)
d W m = d 1 2 L I 2 = N I 2 d Λ .

The rotor deflection angle dθ is converted into electric angle, then dθe=Zrdθ, namely:

(8)
d θ = d θ e Z r .

Substituted Eq. (7) and (8) into Eq. (4):

(9)
T = ( N I ) 2 Z r d Λ d θ e = F δ 2 Z r d Λ d θ e .

In the equation, Fδ is magnetomotive force of each phase windings, and Fδ=NI; Zr is the teeth number of rotor; dΛdθe is the change rate of air gap magnetic conductivity for rotor angle deflection.

Stepping motor air gap Λ can be showed as ratio magnetic conductance λ, it is defined as one blade space magnetic conductance between stator and rotor on unit core length of motor, so air gap magnetic conductivity can be written as:

(10)
Λ = Z S l λ .

In this equation, Zs are numbers of small teeth of each stator pole; l is the length of core. Due to the λ is a periodic function of rotor position θe, and the period are teeth number θte. The λ can be shown in the form of Fourier series:

(11)
λ = λ 0 + n = 1 λ n cos n θ s e .

In the type, λ0 is the average of air gap ratio magnetic conductance; λn is n harmonic amplitude in air gap ratio magnetic conductance; the Eq. (11) can be expressed as Eq. (12) if the high harmonics are ignored:

(12)
λ = λ 0 + λ 1 c o s θ s e .

Substituted Eq. (12) and (10) into Eq. (9):

(13)
T = - F δ 2 Z S Z r l λ 1 s i n θ e = - T m a x s i n θ e .

Visibly, the static torque is produced when single phase windings of motor going through DC current, and it is sine function of disorder angle θe if ignoring the influence of high harmonics. Its action are always make the rotor disorder angle tending to zero, and the direction is contrast to θe increasing. The static torque T is proportional to the secondary of current I under the condition of stable structure and unsaturated magnetic circuit.

The stepping motor is effected by air friction drag torque when in actual no load running, according to the torque balance equation:

(14)
T - γ d θ d t = J d 2 θ d t 2 = T V i b r a t i o n T o r q u e d 2 θ d t 2 .

In this equation, γ is coefficient of air resistance moment; J is rotary inertia of system. Substituted Eq. (8) and (13) into Eq. (14):

(15)
T m a x s i n θ e + γ Z r d θ e d t + J Z r 2 d 2 θ e d t 2 = 0 .

3. Torsional pendulum motion model

Input an initial force to the torsion pendulum, the spiral spring turned angle θl, and generated the restoring moment Ts=-kθl, and k is the stiffness coefficient of spiral spring. The initial force disappeared, and the whole system begin to reverse swing under the action of restoring moment, and in this process, the system is also effected by air resistance moment T2s=-γldθldt, γl is the air resistance coefficient. Set the inertia of torsion pendulum system is Jl, according to torque balance equation:

(16)
J l d 2 θ l d t 2 + γ l d θ l d t + k θ l = 0 .

Set 2β=γlJl and ωn2=kJl, β is system damping coefficient, ωn is system natural oscillation frequency, so:

(17)
d 2 θ l d t 2 + 2 β d θ l d t + ω n 2 θ l = 0 .

Then we know θl(t)=θ01-ξ2e-βtcos(ωdt-ϕ), θ0 is the initial amplitude of torsional pendulum angular displacement; ωd=ωn1-ξ2 is system damped oscillation frequency; damping ratio is ξ=γl2kJl=βωn; and ϕ=arcsinξ.

Visible, the regularity of torsional pendulum angle is sine, and the amplitude decreased with time, thus the recovery torque of torsion pendulum system is:

(18)
T s = - k θ l = - T m a x l c o s ω d t - ϕ .

Thus, the mathematical expressions form of torsion pendulum system restoring torque is similar to stepper motor static torque, and its properties are also drive system reaches equilibrium state, namely the drive torque is zero on this position.

The differences are Tmaxl is proportional to initial angle θ0, and torsion angle always decreased with time, so the restoring torque is reducing in the whole process. Stepper motor static torque amplitude Tmax is proportional to square of motor windings average current I, and the rotation angle is also decreasing with time, but the static torque is not always reduced, such as the initial angle is big than π2 , the static torque changed into big first and then into small.

4. Experimental results and mechanism analysis

The experimental platform is made by Fangyuan Industrial Technology Co. Ltd. of Zhejiang University, its type is NMCL-II, stepping motor model is M10, the DC resistances of each phase windings are 45 Ω, and the supporting drive power type is NMEL-10.

4.1. Experiment 1

Stepping motor and homemade vibration torque sensor are connected coaxially, three phase three clap and six clap of single step operation without load, the vibration torque of motor are shown in Figures 1 and 2 respectively.

Fig. 1. Vibration torque on three phases single three claps

 Vibration torque on three phases single three claps

Fig. 2. Vibration torque on three phases six claps

 Vibration torque on three phases six claps

From the experimental results, the vibration torque amplitude of stepper motor operation in three-phase three beat is greater than six beat, the reasons are shown in Fig. 3.

Fig. 3. Torque angle characteristic of stepping motor

Torque angle characteristic of stepping motor

In Fig. 3, sine wave A is the torque angle characteristic of stepping motor when A phase windings going through DC current; sine wave B is the torque angle characteristic of stepping motor when B phase windings going through DC current; sine wave AB is the torque angle characteristic of stepping motor when A phase windings and B phase windings going through DC current at the same time; then:

(19)
T A = - T m a x s i n θ e , T B = - T m a x s i n θ e - 120 , T A B = - T m a x s i n θ e - 60 .

Energizing sequence for stepping motor working on single three clap is AB; and six clap is AAB. The staring position is a, which is the balance position of angle characteristic A, and this time θe=0; if B phase windings going through DC current, the angle characteristic change into B instantaneously, and the static torque is:

(20)
T B = - T m a x × s i n - 120 = T m a x 3 2 .

If A and B phase windings going through DC current simultaneously, the angle characteristic change into AB instantaneously, and the static torque is:

(21)
T A B = - T m a x × s i n - 60 = T m a x 3 2 .

Therefore, if stepping motor operation on single step, the torque fluctuation amplitude of three clap and six clap are same in theoretically when take changing electric pulses, but due to the rotor must be turning and then the sensor output signals, it is namely the stepping motor rotor has turned certain angles. As shown in Fig. 3, it is known that the static torque of three clap is bigger than six clap, then leading to the large vibration torque.

In addition, when the stepper motor operating on no-load, as shown in Fig. 3, the acceleration interval of three clap is 120 electrical degrees, and the acceleration is increased firstly and then decreased; the acceleration interval of six clap is 120 electrical degrees, and the acceleration is always decreased.

The rotor is effected not only by the static torque, but also by the air friction resistance torque and bearing friction resistance torque during actual operation, so the acceleration interval is less than the theoretical value, and we still know that acceleration interval of three clap is greater than six clap. So, the rotor speed at balance position b of three clap operation is greater than at the balance position ab of six clap operation. Due to the inertia, the rotor is not stop immediately when reaches the equilibrium position, and its rotation direction will away from the equilibrium position, which is similar to the pendulum motion process. The function of static torque is drive the rotor to return the balance position, so stepping motor rotor will produce reciprocating motion, and the six clap operation will generate obvious reciprocating oscillation than three clap operation.

4.2. Experiment 2

Stepping motor and homemade vibration torque sensor are connected coaxially, three claps of single step operation without load, and the frequency of pluses are about 2 Hz, the vibration torques are shown in Fig. 4.

According to Fig. 4, the vibration torque of rotor do not influence each other on low frequency continuous operation, because the next step is always coming when the forth step reaches steady state.

4.3. Experiment 3

Stepping motor and homemade vibration torque sensor are connected coaxially, three claps of single step operation without load, and the frequency of pluses are about 170 Hz, the vibration torque is shown in Fig. 5.

The reason of stepper motor lost steps is shown in Fig. 6, the stepper motor rotor is located in the equilibrium position a when the A phase windings going through DC current, then the A phase windings power cut off and B phase windings power turn on, the rotor is turning to balance position b under the affection of static torque, and it is an accelerated state in the process; when the rotor reaches balance position b, a reply oscillation will be generated on the rotor due to the inertia; if the rotor is at the position d, then B phase windings power cut off and C phase windings power turn on, the rotor is subjected to static torque which direction is negative as shown in Fig. 10, and the speed at the same direction, so the rotor is to acceleration at negative direction, which causes the stepping motor generates lost step; and according to the graph 6, there is a static torque jump at position d,which changes static torque positive into negative, and it is always negative in the location process de; if the A phase windings going through DC current at this process, there will generate another static torque jump, which is driving the static torque from negative to positive; and the stepper motor is generated intense torque fluctuation if this processes is repeated.

Fig. 4. Vibration torque on low frequency continuous operation

 Vibration torque on low frequency continuous operation

Fig. 5. Vibration torque on low frequency oscillation of lost step

 Vibration torque on low frequency oscillation of lost step

4.4. Experiment 4

Stepping motor and homemade vibration torque sensor are connected coaxially, three claps of single step operation without load, and the frequency of pluses are about 300 Hz, 375 Hz and 460 Hz, the vibration torques are shown in Fig. 7.

Fig. 6. Principle analysis of low frequency lost step oscillation

 Principle analysis of low frequency lost step oscillation

Fig. 7. Vibration torque on differences pulse frequency

 Vibration torque on differences pulse frequency

Fig. 8. Stepping motor driver circuit of each phase windings

 Stepping motor driver circuit of each phase windings

According to experimental results, the stepper motor torque fluctuation is reduced gradually with increasing the frequency, and compared Fig. 1 with Fig. 7, the torque ripple amplitude on continuous operation is smaller than on the single-step operation. The reason for this is shown in Fig. 8.

Fig. 8 is the driving circuit of each phase windings of the stepping motor, uKe is the driving voltage pulses; V1 is the triode, in order to amplify the windings current, Rf1 is the limiting resistor; the windings discharge circuit is formed by V2 and Rf2; U is the DC voltage of windings; when a pulse through the triode V1, the current i is generated in phase windings.

Due to the phase windings existence inductance, the phase windings average current will be reduced when the pulse frequency is increased, as the Eq. (13) shows, the maximum static torque is proportional to the square average current of windings, so the amplitude of motor torque angle characteristic is decreased with the pulses frequency increased. Besides, the stepper motor speed is up with improving the pulse frequency, then the friction resistance moment becomes larger, and according to the torque balance Eq. (14), the torque fluctuation amplitude will be reduced.

Table 1. Experiments data

Frequency of stepper motor / Hz
Vibration torque sensor peak voltage / V
Average voltage / V
200
0.198
0.192
0.194
0.195
225
0.136
0.138
0.14
0.138
250
0.094
0.092
0.094
0.093
275
0.08
0.078
0.08
0.079
300
0.062
0.064
0.062
0.063
325
0.05
0.052
0.05
0.051
350
0.04
0.042
0.042
0.041
375
0.032
0.036
0.03
0.034
400
0.026
0.024
0.024
0.025
425
0.018
0.016
0.018
0.017
450
0.012
0.01
0.008
0.01

Fig. 9. Stepping motor vibration torque peak value-frequency curve

 Stepping motor vibration torque peak value-frequency curve

In order to acquire the relationship between the stepper motor vibration torque peak value and frequency of the normal continuous operation, experiment 4 is operated several times, and the results data are shown in Table 1, the experiments data of Table 1 were fitted by a polynomial function as shown in Fig. 9, and the corresponding mathematical expressions are obtained as Eq. (22):

(22)
f x = i = 1 5 p i x 5 - i .

In this equation p1=1.499e-010,p2=-2.196e-007,p3=0.00012,p4=-0.02937,p5=2.786 and the SSE is 0.0000542; R-square is 0.9983; adjusted R-square is 0.9972; RMSE is 0.0003005.

5. Conclusions

The static torque of stepping motor is induced through the law of conservation of energy in this article, so the stepper motor mathematical model is gained, and it is studied by comparing with the torsion pendulum motion mathematical model. Stepper motor torque fluctuations were tested, including single-step operation, low frequency operation, low frequency lost step operation and high frequency operation, the test results are analyzed, and the corresponding conclusion are received. In order to reduce the torque fluctuation of motor on the practical operation, the following work and research is to design a new step motor controller, which controls the winding current input according to the torque fluctuation, it can improve the performance and extend the application range of stepping motor.

Acknowledgements

The work was sponsored by in part by the Jiaxing Science and Technology Research Project of China under Grant (2012AY1021), in part by the Zhejiang Provincial Department of Education Scientific Research Project of China under Grant (Y201226082), in part by the Zhejiang Province Major Science and Technology Project of China under Grant (2008C01011-4) and in part by the Jiaxing University Nanhu College Key Scientific Research Project.

References

  1. P. Krishnamurthy, F. Khorrami TriM: an ultra-accurate high-speed six degree of freedom manipulator using planar stepper motors. Journal of Intelligent and Robotic Systems, Vol. 51, Issue 2, 2008, p. 137-157. [Search CrossRef]
  2. Maki K. Rashid, Zahi A. Khalil Configuration design and intelligent stepping of a spherical motor in robotic joint. Journal of Intelligent and Robotic Systems, Vol. 40, Issue 2, 2004, p. 165-181. [Search CrossRef]
  3. P. Melin, O. Castillo Intelligent control of a stepping motor drive using a hybrid neuro-fuzzy approach. Soft Computing, Vol. 8, Issue 8, 2004, p. 546-555. [Search CrossRef]
  4. Tomoaki Yano Simulation results of a hexahedron octahedron based spherical stepping motor. Journal of Mechanical Science and Technology, Vol. 24, Issue 1, 2010, p. 33-36. [Search CrossRef]
  5. Hong Min Zheng, Wen Feng Cui Design of controlling system for linear stepper motor based on DSP2407. Advances in Electronic Commerce, Web Application and Communication Advances in Intelligent and Soft Computing, Vol. 149, 2012, p. 335-341. [Search CrossRef]
  6. Bo Qu, Hong Lin Design of the stepper motor controller. Electronics and Signal Processing, Lecture Notes in Electrical Engineering, Vol. 97, 2011, p. 525-531. [Search CrossRef]
  7. Baoshan You, Liying Pei A novel stepper speed control algorithm based on FPGA. Advances in Computer Science and Education Applications, Communications in Computer and Information Science, 2011, p. 202-209. [Search CrossRef]
  8. Pan Jin-Yu Analysis and determination of stepping motor rotation stability. Aerospace Shanghai, Vol. 22, Issue 5, 2005, p. 54-57. [Search CrossRef]
  9. Sun Yao-Jie, Zuo He, Kang Long-Yun Time-varying repetitive control of restraining torque ripple for hybrid stepping motor system. Proceedings of ME CSEE, Vol. 24, Issue 11, 2004, p. 18-17. [Search CrossRef]
  10. Zheng Xiao-Dong, Zou Kun, Sun Yi-Ze Rotating speed fluctuation and response speed of two-phase hybrid type stepping motor. Journal of Donghua University (Natural Science), Vol. 35, Issue 4, 2009, p. 463-487. [Search CrossRef]
  11. Zhao Hao, Feng Hao Research of a new permanent magnet angular acceleration sensor based on electromagnetic induction. Chinese Journal of Sensors and Actuators, Vol. 25, Issue 9, 2012, p. 1257-1261. [Search CrossRef]