DC Motor Parameter Identification Using Speed Step Responses

Abstract

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 coefficients are related to the motor parameters. These coefficients can be easily computed using existing curve fitting 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.

(ProQuest: ... denotes non-US-ASCII text omitted.)

Recommended by F. Gao

Flight Control and Navigation Group, Rockwell Collins, Warrenton, VA 20187, USA

Received 8 April 2012; Revised 31 August 2012; Accepted 3 September 2012

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.

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 specifications, 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 identification methods [ 1- 4] can be applied to DC motor model identification. In particular, various methods have been applied to DC motor parameter identification; that is, [ 5, 6] used the algebraic identification 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. Identified 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 effective system identification approach based on the motor input and output is desirable and valuable, especially for the field applications and quick controller prototyping. In this paper, a DC motor parameter identification 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 coefficients of the Taylor series are established. In the implementation, the motor speed response under a constant voltage is sampled, then fit the samples to obtain the coefficients 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 coefficients. 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 effective. The curve fitting 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: [figure omitted; refer to PDF] where ω is the motor speed, V is the motor terminal voltage, i is the winding current, _{k b} is the back-EMF constant of the motor, _{k t} is the torque constant, R is the terminal resistance, L is the terminal inductance, J is the motor and load inertia, and _{T d} is the disturbance torque. _{T d} is a combination of the cogging torque, _{T cog} , the kinetic friction, _{T f} , and the viscous friction (viscous damping force): [figure omitted; refer to PDF] where c is the damping coefficient. According to ( 1), the velocity response in the Laplace domain is [figure omitted; refer to PDF] where _{t e} = L / R is the electrical time constant, _{t m} = RJ / _{k t k b} is the mechanical time constant, and s is the Laplace variable.

Based on these equations, we would like to know _{t m} , _{t e} , _{T d} , J , and so forth, by measuring the velocity response under a known, controlled voltage input. In this paper, we consider two application situations: the first 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: [figure omitted; refer to PDF]

The transfer function can be factorized into [figure omitted; refer to PDF] where [figure omitted; refer to PDF]

Assumption 2.1.

It is assumed here that there are two distinct real poles; that is, _{t m} >4 _{t e} .

For a constant voltage input V (s ) = _{V 0} /s , the speed response is [figure omitted; refer to PDF] where [figure omitted; refer to PDF]

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 [figure omitted; refer to PDF] Expanding the term _{α 3} / (s +b ) , we get [figure omitted; refer to PDF]

Combining the three terms together, we have the total speed response: [figure omitted; refer to PDF] where _{β 0} =ab , _{β 1} = -ab (a +b ) , and _{β 2} =ab ( ^{a 2} +ab + ^{b 2} ) .

According to ( 6), [figure omitted; refer to PDF] Thus, we have [figure omitted; refer to PDF]

The above equation allows us to calculate the mechanical and electrical time constants _{t m} and _{t e} using the coefficients of the power series in ( 11). These coefficients can be obtained by curve fitting 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 [figure omitted; refer to PDF] 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 ) = _{V 0} /s and the constant disturbance (ignore the cogging torque or consider the average cogging torque effect on speed over one revolution is zero) _{T d} (s ) = _{T 0} /s , the speed response is [figure omitted; refer to PDF]

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: [figure omitted; refer to PDF]

Based on these coefficients, we have [figure omitted; refer to PDF] and another equation for a +b : [figure omitted; refer to PDF] Then, we can express the motor parameters as [figure omitted; refer to PDF]

In practice, fit 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: [figure omitted; refer to PDF] where _{β 0} = _{T 0} / J and _{ω ss} is the motor steady-state angular speed.

3. Implementation and Results

The proposed approaches were first 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 (field 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 specifications for reference are presented in Table 1.

Table 1: RK370CA parameter values.

Parameter	Value	Unit
Terminal resistance	17 ± 15%	Ω
Terminal inductance	N/A	Henry
Torque constant	18.3 ± 18%	mNm/A
Mass moment of inertia	9.0	gc m 2
Counter-electromotive force	0.0233	volt/(rad/sec)

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 _{t e} = 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, effects on the speed response are significant when the input voltage is small. To demonstrate the effectiveness of the algorithm, we sent two voltages, 2 volt and 10 volt, to the motor. Driving the motor at two different voltage levels can demonstrate that the viscous friction varies with the speed, also can allow us to calculate the viscous damping coefficient.

Usually _{t e} is very small compared to _{t m} a good estimate of both _{t e} and _{t m} at the same time is difficult. Because _{t m} is usually much larger than _{t e} , _{t m} and _{t e} were estimated separately using different data collected with different sample rates and different time durations. For estimating _{t m} , the motor speed in both the transient phase and the steady-state was sampled at 1 kHz for one second; for estimating _{t e} , 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 gcm² , _{k b} =0.0233 volt/(rad/sec), and _{k t} =0.0183 Nm/A are used to calculate _{t m} in the fourth column of the table. According to Table 2, the estimates of _{k t} , _{t m} , and _{t e} are in good agreement with those given by the motor specifications.

Table 2: RK370CA test results.

Parameter	w / o dist. 20 v	w /dist. 2 v/10 v	Spec. (meas.)	Unit
k t	0.0238	0.0207/0.0169	0.0183 ± 18%	Nm / A
t m	0.0407	0.0211/0.0203	0.0359	sec
t e	0.00554	0.00122/0.00134	(0.00122)	sec
T 0 /J	N/A	10.551/115.758	N/A	Nm / kg m 2

Time responses sampled at 1 kHz for 1 s are given in Figures 1, 2and 3. In these figures, red curves represent the power series, ^{∑ i =1 n} _{x i} ^{t i} , resulting from curve fitting the motor speed responses. Comparing these figures, 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.

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

[figure omitted; refer to PDF]

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

[figure omitted; refer to PDF]

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

[figure omitted; refer to PDF]

To further demonstrate the effectiveness of the proposed algorithms, we compared them to conventional identification 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 identification toolbox to identify a second-order model. Various methods, that is, subspace approach in the system identification toolbox, were tried and compared. The best model found was [figure omitted; refer to PDF] Using the model coefficients, we get [figure omitted; refer to PDF] These estimates are bad, especially the electrical time constant _{t e} due to the very small time scale as alluded to earlier.

Remark 2.

_{T 0} / J may be used to calculate the friction (both kinetic and viscous) if J is known. First, calculate the viscous friction coefficient c = ( _{T 1} - _{T 0} ) / ( _{ω 1} - _{ω 0} ) . Then, calculate the dynamic friction, _{T f} = _{T 0} -c _{ω 0} . For example, _{ω 0} =1.21 ips under 2 volt, _{ω 1} =6.274 ips under 10 volt, J =9.0 gcm² , and it renders c =0.0187 mNm/ips.

Remark 3.

The number of terms in the power series included for fitting the data was determined through trial and error. When disturbance was not considered, twenty-five terms were included; when disturbance was considered, including forty terms gave the best results. Since the coefficients were calculated using the polynomial curve fitting function from the math library provided inside LabVIEW, it was not difficult and time consuming to try different 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 _{t e} in this small motor. Algorithm considering disturbance torque was applied. In the testing, 10 volts was used as the motor drive input. For _{t m} estimation, the speed response was sampled at 1000 Hz for 500 samples, while for _{t e} estimation, it was sampled at 6000 Hz for 850 samples.

Table 3: FC130 test results.

Parameter	w /dist.	Spec.	Unit
k t	0.0137	0.0127 ± 10%	Nm / A
t m	0.0208	0.024	sec
t e	0.251	0.214	msec

4. Conclusions

A convenient, effective system identification 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 field applications.

References

[1] L. Ljiung System Identification: Theory for the User , Prentice Hall, 1999., 2nd.

[2] H. Unbehauen, G. P. Rao, "A review of identification in continuous-time systems," Annual Reviews in Control , vol. 22, pp. 145-171, 1998.

[3] G. F. Franklin, J. D. Powell, M. L. Workman Digital Control of Dynamic Systems , Addison Wesley, 1990., 2nd.

[4] J. C. Basilio, M. V. Moreira, "State-space parameter identification in a second control laboratory," IEEE Transactions on Education , vol. 47, no. 2, pp. 204-210, 2004.

[5] G. Mamani, J. Becedas, H. Sira-Ramirez, V. Feliu Batlle, "Open-loop algebraic identification method for DC motors," in Proceedings of the European Control Conference, Kos, Greece, 2007.

[6] G. Mamani, J. Becedas, 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, D. Stojanovic, "Recursive least square method in parameters identification of DC motors models," Facta Universitatis , vol. 18, no. 3, pp. 467-478, 2005.

[8] M. Hadef, M. R. Mekideche, "Parameter identification 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. Hoffman, T. Betran, "Optimal state space control of DC motor," in Proceedings of the 17th World Congress IFAC, pp. 5796-5801, Seoul, Korea, 2008.

[10] M. Hadef, A. Bourouina, M. R. Mekideche, "Parameter identification of a DC motor via moments method," International Journal of Electrical and Power Engineering , vol. 1, no. 2, pp. 210-214, 2008.

[11] A. Rubaai, R. Kotaru, "Online identification and control of a dc motor using learning adaptation of neural networks," IEEE Transactions on Industry Applications , vol. 36, no. 3, pp. 935-942, 2000.

[]

Appendix

Appendix

Coefficients for no disturbance case are as follows: [figure omitted; refer to PDF]

Coefficients for disturbance case are as follows: [figure omitted; refer to PDF]