Hindawi Publishing Corporation Modelling and Simulation in Engineering Volume 2012, Article ID 189757, 5 pages doi:10.

1155/2012/189757

Research Article DC Motor Parameter Identication Using Speed Step Responses

Wei Wu
Flight Control and Navigation Group, Rockwell Collins, Warrenton, VA 20187, USA Correspondence should be addressed to Wei Wu, wu [email protected] Received 8 April 2012; Revised 31 August 2012; Accepted 3 September 2012 Academic Editor: F. Gao Copyright 2012 Wei Wu. 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. Based on the DC motor speed response measurement under a step voltage input, important motor parameters such as the electrical time constant, the mechanical time constant, and the friction can be estimated. A power series expansion of the motor speed response is presented, whose coecients are related to the motor parameters. These coecients can be easily computed using existing curve tting methods. Experimental results are presented to demonstrate the application of this approach. In these experiments, the approach was readily implemented and gave more accurate estimates than conventional methods.

1. Introduction
DC motors have wide applications in industrial control systems because they are easy to control and model. For analytical control system design and optimization, sometimes a precise model of the DC motor used in a control system may be needed. In this case, the values for reference of the motor parameters given in the motor specications, usually provided by the motor manufacturer, may not be considered adequate, especially for cheaper DC motors which tend to have relatively large tolerances in their electrical and mechanical parameters. General system identication methods [14] can be applied to DC motor model identication. In particular, various methods have been applied to DC motor parameter identication; that is, [5, 6] used the algebraic identication method, [7] used the recursive least square method, [8] applied the inverse theory, [9] used the least square method, and [10] applied the moments method. Identied DC motor models are often subsequently used for controller design and/or optimization, for example, [6, 9, 11]. Without expensive testing apparatus and a long testing cycle, a quick and eective system identication approach based on the motor input and output is desirable and valuable, especially for the eld applications and quick controller prototyping. In this paper, a DC motor parameter identication approach based on the Taylor series expansion of the motor speed response under a constant voltage input is presented. The relationships between the motor parameters

and the coecients of the Taylor series are established. In the implementation, the motor speed response under a constant voltage is sampled, then t the samples to obtain the coecients of power terms in the Taylor series. Then, the DC motor mechanical and electrical time constants, back-EMF, and the friction can be computed using these coecients. With the knowledge of these parameters, a precise motor model is obtained for the subsequent controller design. For application point of view, this approach requires only a speed/position sensor, such as an optical encoder, and a voltage power supply, no current measurement is needed and the motor is run in open loop; thus it is practical and cost eective. The curve tting can be performed using many existing methods, such as the least square method, and these optimization methods are widely available in commercial computing packages such as Matlab and LabVIEW.

2. Main Results
Consider the following DC motor governing equations: L di + iR + kb = V , dt

(1) d = kt i + Td , J dt where is the motor speed, V is the motor terminal voltage, i is the winding current, kb is the back-EMF constant of the motor, kt is the torque constant, R is the terminal resistance,

2 L is the terminal inductance, J is the motor and load inertia, and Td is the disturbance torque. Td is a combination of the cogging torque, Tcog , the kinetic friction, T f , and the viscous friction (viscous damping force): Td = Tcog + T f + c, (2)

Modelling and Simulation in Engineering Combining the three terms together, we have the total speed response: (t) = V0 1 2 1 3 1 0 t + 1 t + 2 t 4 + , kb 2 6 24 (11)

where c is the damping coecient. According to (1), the velocity response in the Laplace domain is (s) = (1/J)tm (te s + 1) 1/kb V (s) + Td (s), (3) 2+t s+1 tm te s tm te s2 + tm s + 1 m

where 0 = ab, 1 = ab(a + b), and 2 = ab(a2 + ab + b2 ). According to (6), ab = Thus, we have tm = 1 2, 0 te = 0 . 1 (13) 1 , tm te a+b = 1 . te (12)

where te = L/R is the electrical time constant, tm = RJ/kt kb is the mechanical time constant, and s is the Laplace variable. Based on these equations, we would like to know tm , te , Td , J, and so forth, by measuring the velocity response under a known, controlled voltage input. In this paper, we consider two application situations: the rst situation is that the disturbance torque is negligible, while in the second one, the disturbance needs to be considered. 2.1. Estimation without the Disturbance Torque. When the voltage speed response dominates; for example, the input voltage is large, we can ignore the disturbance torque in the speed response see (3). In this case, we can consider the following DC motor model: 1/kb (s) = . V (s) tm te s2 + tm s + 1 The transfer function can be factorized into (s) 1/kb , = V (s) tm te (s + a)(s + b) where a, b = 1 1 4te /tm . 2te (6) (5) (4)

The above equation allows us to calculate the mechanical and electrical time constants tm and te using the coecients of the power series in (11). These coecients can be obtained by curve tting the motor speed step response data using power functions. 2.2. Estimation with the Disturbance Torque. Consider that the disturbance torque in the DC motor is not negligible. The disturbance transfer function is (1/J)tm (te s + 1) (s) = . Td (s) tm te s2 + tm s + 1 (14)

Assumption. It is assumed here that there are two distinct real poles; that is, tm > 4te . For a constant voltage input V (s) = V0 /s, the speed response is (s) = where 1 = V0 , kb 2 = V0 b , kb a b 3 = V0 a . kb b a (8) V0 /kb 1 = + 2 + 3 , tm te s(s + a)(s + b) s s+a s+b (7)

Disturbance torque generally consists of the cogging torque and the friction torque. The cogging torque is quite complicated and is not addressed here. Both the kinetic and viscous frictions are considered and are assumed to be constant on average under a constant motor speed. Given a constant motor terminal voltage V (s) = V0 /s and the constant disturbance (ignore the cogging torque or consider the average cogging torque eect on speed over one revolution is zero) Td (s) = T0 /s, the speed response is (s) = V0 (1/J)tm (te s + 1) T0 1/kb + . tm te s2 + tm s + 1 s tm te s2 + tm s + 1 s (15)

As in the previous section, applying the partial fraction expansion of the step response in the Laplace domain, then expanding the exponential terms in the time domain using the Taylor series, we obtain the total step response in the time domain: (t) = 0 t + 1 t 2 2 t 3 + 3 t 4 . Based on these coecients, we have ab =
2 182 243 1 2 , 30 2 + 21

(16)

Consider the three terms in the step response one at a time. 1 /s is a step function in the time domain; both 2 /(s + a) and 1 /(s + b) are exponential functions in the time domain and can be expanded using the Taylor series. Expanding the term 2 /(s + a), we get 1 V0 b 1 1 at + a2 t 2 a3 t 3 + . kb a b 2 6 Expanding the term 3 /(s + b), we get V0 a 1 1 1 bt + b2 t 2 b3 t 3 + . kb b a 2 6 (10) (9)

62 0 ab , a+b = 21 and another equation for a + b: a+b = 123 + 1 ab . 22

(17)

(18)

Modelling and Simulation in Engineering

Table 1: RK370CA parameter values. Parameter Terminal resistance Terminal inductance Torque constant Mass moment of inertia Counter-electromotive force Value 17 15% N/A 18.3 18% 9.0 0.0233 Unit Henry mNm/A gcm2 volt/(rad/sec) Parameter kt tm te T0 /J
1.3 1.2 1.1 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 0

3
Table 2: RK370CA test results. w/o dist. 20 v 0.0238 0.0407 0.00554 N/A w/dist. Spec. (meas.) Unit 2 v/10 v 0.0207/0.0169 0.0183 18% Nm/A 0.0359 sec 0.0211/0.0203 sec 0.00122/0.00134 (0.00122) N/A Nm/kgm2 10.551/115.758

Then, we can express the motor parameters as tm = te = a+b , ab 1 , a+b

T0 = 0 , J kb = ab V0 . 21

(19)

Velocity (ips)

100 200 300 400 500 600 700 800 900 1000 Time (ms)

In practice, t the measured motor speed step response using power functions according to (16); then calculate the motor parameters using (19). Remark 1. Another relationship useful for checking the algorithm is based on the steady-state response of (15), expressed by the following equation: V0 + 0 tm = ss , kb (20)

Figure 1: Approach w/consideration of disturbance under 2 volt input: black, measurement: red, power series tted.

where 0 = T0 /J and ss is the motor steady-state angular speed.

3. Implementation and Results

The proposed approaches were rst applied to a Mabuchi RK370CA motor, then a Mabuchi FC130 motor. To implement the algorithms, a LabVIEW program was created to interface a pulse width modulated (PWM) motor drive and an optical encoder with quadrature digital outputs mounted on the motor shaft. The determinism of the sample time was assured by the LabVIEW real-time module. And, a national instrument (NI) LabVIEW FPGA (eld programmable gate array) card was utilized to process the digital quadrature encoder signals to obtain the motor speed and to control the motor PWM drive. Values of the motor parameters given in the motor specications for reference are presented in Table 1. Note that the Back-EMF and torque constant are not equal (although it should be theoretically). Inductance value is not given and was measured as 20.25 Henry. The resistance was measured as 16.4 . Thus we calculated te = L/R = 0.00122 sec. First, apply the algorithm for no-disturbance torque. To apply this algorithm, the speed response part due to the voltage input is assumed to dominate. To meet this condition,

for example, the speed variation at the steady-state is small compared to the steady-state speed, we send a large voltage to the motor drive, V = 20 volt. Next, we apply the approach for disturbance torque. The disturbances, that is, friction, eects on the speed response are signicant when the input voltage is small. To demonstrate the eectiveness of the algorithm, we sent two voltages, 2 volt and 10 volt, to the motor. Driving the motor at two dierent voltage levels can demonstrate that the viscous friction varies with the speed, also can allow us to calculate the viscous damping coecient. Usually te is very small compared to tm a good estimate of both te and tm at the same time is dicult. Because tm is usually much larger than te , tm and te were estimated separately using dierent data collected with dierent sample rates and dierent time durations. For estimating tm , the motor speed in both the transient phase and the steady-state was sampled at 1 kHz for one second; for estimating te , the motor speed in the transient phase was sampled at 8 kHz for 200 msec. In each test, the motor was driven multiple times and parameter estimates were averaged. Results are summarized in Table 2. Column two gives the values estimated using the algorithm for no disturbance, and column three gives the values obtained using the algorithm considering disturbance; values in the fourth column are computed using values from Table 1. Note that R = 17 , J = 9 gcm2 , kb = 0.0233 volt/(rad/sec), and kt = 0.0183 Nm/A are used to calculate tm in the fourth column of the table. According to Table 2, the estimates of kt , tm , and te are in good agreement with those given by the motor specications. Time responses sampled at 1 kHz for 1 s are given in Figures 1, 2 and 3. In these gures, red curves represent the power series, n=1 xi t i , resulting from curve tting the motor i

4
Table 3: FC130 test results. Parameter kt tm te w/dist. 0.0137 0.0208 0.251 Spec. 0.0127 10% 0.024 0.214 Unit Nm/A sec msec
1.3 1.2 1.1 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 0

Modelling and Simulation in Engineering

T(s) =

9078 . s2 + 334.6s + 18860 te = 3 msec,

(21)

Using the model coecients, we get tm = 0.0177 sec, (22)

kt = 0.031 volts/(rad/sec).

Velocity (ips)

speed responses. Comparing these gures, it is obvious that the approach with disturbance consideration approximates the measurements much better, because of the existence of the linear term, 0 t, in the power series due to the presence of the constant disturbance in the motor. To further demonstrate the eectiveness of the proposed algorithms, we compared them to conventional identication approaches. First, we drove the motor using random voltage input (10 volts maximum) and measured the motor speed at a sampling rate of 10 kHz. Then, the motor/drive frequency response function was calculated through spectral analysis. Based on the calculated frequency response data, we used Matlab system identication toolbox to identify a second-order model. Various methods, that is, subspace approach in the system identication toolbox, were tried and compared. The best model found was

Velocity (ips)

100 200 300 400 500 600 700 800 900 1000 Time (ms)

Figure 2: Approach w/o consideration of disturbance under 2 volt input: black, measurement: red, power series tted.
14 13 12 11 10 9 8 7 6 5 4 3 2 1 0 0 100 200 300 400 500 600 Time (ms) 700 800 900 1000

These estimates are bad, especially the electrical time constant te due to the very small time scale as alluded to earlier. Remark 2. T0 /J may be used to calculate the friction (both kinetic and viscous) if J is known. First, calculate the viscous friction coecient c = (T1 T0 )/(1 0 ). Then, calculate the dynamic friction, T f = T0 c0 . For example, 0 = 1.21 ips under 2 volt, 1 = 6.274 ips under 10 volt, J = 9.0 gcm2 , and it renders c = 0.0187 mNm/ips. Remark 3. The number of terms in the power series included for tting the data was determined through trial and error. When disturbance was not considered, twenty-ve terms were included; when disturbance was considered, including forty terms gave the best results. Since the coecients were calculated using the polynomial curve tting function from the math library provided inside LabVIEW, it was not difcult and time consuming to try dierent number of terms. Including more terms does not necessarily improve the parameter estimation accuracy. A Mabuchi FC130 motor was tested as well. It is a smaller motor compared to RK370. Good results were obtained again this time; see Table 3. Note the very small te in this small motor. Algorithm considering disturbance torque was applied. In the testing, 10 volts was used as the motor drive input. For tm estimation, the speed response was sampled at 1000 Hz for 500 samples, while for te estimation, it was sampled at 6000 Hz for 850 samples.

Figure 3: Approach w/o consideration of disturbance under 20 volt input:Black, measurement: Red, power series tted.

4. Conclusions
A convenient, eective system identication approach is proposed to estimate the DC motor torque constant, mechanical time constant, electrical time constant, and friction. This approach was implemented on two Mabuchi motors, and the great test results were presented. This open-loop method requires little hardware, only a speed/position sensor and a voltage supply. The estimated motor parameters can be used to verify the DC motor performance or be used to build a model of the motor for the subsequent controller design or system optimization. This approach is especially suited to quick eld applications.

Appendix
Coecients for no disturbance case are as follows: 0 = ab, 1 = ab(a + b), 2 = ab a2 + ab + b2 . (A.1)

Modelling and Simulation in Engineering Coecients for disturbance case are as follows: 0 = 1 = 2 = T0 , J 1 V0 ab, 2 kb 1 V0 T ab(a + b) 0 a2 + ab + b2 6 kb J + 3 = T0 tm ab(a + b) , J

1 V0 T ab a2 + ab + b2 0 a3 + a2 b + ab2 + b3 24 kb J + T0 tm ab a2 + ab + b2 J . (A.2)

References
[1] L. Ljiung, System Identication: Theory for the User, Prentice Hall, 2nd edition, 1999. [2] H. Unbehauen and G. P. Rao, A review of identication in continuous-time systems, Annual Reviews in Control, vol. 22, pp. 145171, 1998. [3] G. F. Franklin, J. D. Powell, and M. L. Workman, Digital Control of Dynamic Systems, Addison Wesley, 2nd edition, 1990. [4] J. C. Basilio and M. V. Moreira, State-space parameter identication in a second control laboratory, IEEE Transactions on Education, vol. 47, no. 2, pp. 204210, 2004. [5] G. Mamani, J. Becedas, H. Sira-Ramirez, and V. Feliu Batlle, Open-loop algebraic identication method for DC motors, in Proceedings of the European Control Conference, Kos, Greece, 2007. [6] G. Mamani, J. Becedas, and V. Feliu-Batlle, On-line fast algebraic parameter and state estimation for a DC motor applied to adaptive 16 control, in Proceedings of the World Congress on Engineering, London, UK, 2008. [7] R. Krneta, S. Antic, and D. Stojanovic, Recursive least square method in parameters identication of DC motors models, Facta Universitatis, vol. 18, no. 3, pp. 467478, 2005. [8] M. Hadef and M. R. Mekideche, Parameter identication of a separately excited DC motor via inverse problem methodology, in Proceedings of the Ecologic Vehicles andRenewable Energies, Monaco, France, 2009. [9] M. Ruderman, J. Krettek, F. Homan, and T. Betran, Optimal state space control of DC motor, in Proceedings of the 17th World Congress IFAC, pp. 57965801, Seoul, Korea, 2008. [10] M. Hadef, A. Bourouina, and M. R. Mekideche, Parameter identication of a DC motor via moments method, International Journal of Electrical and Power Engineering, vol. 1, no. 2, pp. 210214, 2008. [11] A. Rubaai and R. Kotaru, Online identication and control of a dc motor using learning adaptation of neural networks, IEEE Transactions on Industry Applications, vol. 36, no. 3, pp. 935942, 2000.