Research Article  Open Access
Jianhua Liu, Rui Nie, Zheshu Li, "Multimotors Coordination Control Method for Consistent Total Traction Torque of Electric Locomotives", Journal of Advanced Transportation, vol. 2020, Article ID 2092463, 11 pages, 2020. https://doi.org/10.1155/2020/2092463
Multimotors Coordination Control Method for Consistent Total Traction Torque of Electric Locomotives
Abstract
Multimotors traction is a common driving mode for electric locomotives. The complexity and variability of the operating environment often cause difficulties for one or some traction motors to output the distributed traction torque, which requires the redistribution of the torque output of multimotors to maintain the total traction torque of all motors at a given constant. To this end, a multimotors coordination control method was proposed for consistent total traction torque of electric locomotives. Through this approach, a traction torquecoordinated control structure was constructed on the basis of energy minimisation with an optimal dynamic adjustment mode. The sliding mode tracking control was developed for each independent motor to improve its dynamic response performance. A simulation verification of motors fault types in two different cases was conducted. The results showed that the proposed method could track the total traction torque in a short time, and the system fluctuations in independent motor coordination process were very small.
1. Introduction
Most of the existing electric locomotives use multiple motors together to provide traction torque [1]. For example, the traction power of the SS_{8}type electric locomotive is provided by four motors together. When the locomotive is under traction condition, the microcomputer control system calculates the total traction torque required according to the driver’s handling instruction and then distributes it to multimotors mounted on the bogie according to the present coordination strategy [2–5].
The locomotive can be safely and smoothly operated only when the sum of all the traction torque provided by multimotors is consistent with the total traction torque required [6]. However, the complexity and variability of the operating environment often cause one or some traction motors to lose traction performance. For example, when the locomotive suddenly enters the rail covered with ice and snow from the normal rail surface, some power wheels often lose traction due to slippage [7–9], resulting in the decrease of the total traction torque. Especially for locomotive on uphill slopes, the backward failure may occur if the required total traction torque cannot be maintained [10]. This phenomenon needs to solve two problems: one is to coordinate the output traction torque among multimotors ensuring that the vector sum of all the traction torque provided by multimotors is consistent with the total traction torque required, and another is to achieve the optimal dynamic performance for independent motors.
In the field of multimotors coordination control methods, scholars at home and abroad have conducted considerable research work regarding decentralised coordination, adaptive robust, and circular crosscoupling control methods [11–13]. In [14], the multimotors coordination control system was considered as an integrated system with multiple dynamic interval subsystems, and a decentralised coordination control strategy was established to solve the system uncertainty. In [15], a coordinated selfadaptive robust control method was proposed to weaken the influence of load variation and system uncertainty. In [16], a consistency algorithm was introduced into the multiaxis synchronous control system to achieve the highperformance synchronous control of the multimotors system. In [17], a multimotors synchronous process coupling control system was designed to improve the synchronisation performance by changing the synchronisation and tracking coefficients in real time in speed feedback and synchronisation compensation. However, according to the latest information mastered by the author, the existing multimotors coordination control method is mostly used to ensure the coordination of a certain state variable, such as position or/and speed. For the problem of consistent total traction torque, although the author and his colleges have made an innovative attempt, many significant scientific problems must be solved because a detailed allocation algorithm of the torque output of the traction system is unavailable in [18].
The dynamic response performance of an independent motor is directly related to the overall performance of the multimotors coordination process. In [19], a suitable quasisliding mode observer with low sampling frequency was designed to estimate the rotor position and stator flux; then the tracking control of the maximum power point of wind turbine was realised. In [20], a fault reconstruction method based on sliding mode variable structure was proposed. In [21], the direct torque and flux control without the beat was introduced into the motor control, which realised the minimisation control of rapid torque and dynamic loss of the motor. In [22], a feedback linearization direct torque control method based on space vector modulation was proposed to reduce the influence of electromagnetic torque and stator flux pulsation.
For electric locomotives, the traction torque of all motors is in the same direction; when one or some motors lose traction performance, there are an infinite number of adjustment methods to coordinate the torque output of the remaining motors, such that the total traction torque is constant. But the optimal dynamic adjustment procedure is minimizing the energy function composed of the traction torques [23–25]. For this purpose, the SS_{8}type electric locomotive with four motors that provide traction torque was taken as the research object in this work; a multimotors coordination control framework with consistent total traction torque was constructed. A multimotors coordination control model based on energy minimisation was proposed to ensure the minimum energy of the multimotors dynamic adjustment process through the optimisation of the weight matrix. Meanwhile, the traction torque control based on sliding mode variable structure was constructed to improve the dynamic response performance of each independent motor.
The structure of this paper is as follows. In Section 2, the SS_{8}type electric locomotive is taken as the research object, and a multimotor coordination control strategy for consistent total traction torque is proposed. Section 3 elucidates the traction torque output coordination method based on energy minimisation. Section 4 designs the sliding mode variable structure controller for each independent motor. Section 5 includes the experimental and data analysis and Section 6 is the conclusion.
2. Multimotors Coordination Control Strategy
The SS_{8}type electric locomotive is an AC/DC transmission phasecontrolled electric locomotive [26]. The locomotive adopts B0–B0 axis structure, which belongs to the fouraxis passenger locomotive, and uses AC–DC traction mode, with each motor equipped with four traction motors and transmits traction torque through a flat pull rod [27]. When the locomotive is running, the locomotive traction control initially calculates the required total traction power after receiving the driver’s handle command and then assigns it to the four motors mounted on the bogie according to a predetermined strategy. In fact, the effect of output traction torque on drive shaft is to provide longitudinal force; this means that the output traction torque of each motor acts in the same direction, so the basic condition for the safe and stable operation of electric locomotives indicates that the sum of all the traction torque provided by all motors should be consistent with the total traction torque required for the locomotive to operate; that is,where is the output traction torque of motor j and is the required total traction torque according to the handle command.
So, a multimotors coordination control strategy is constructed to ensure that the total traction torque is consistent when one or some motors lose traction torque output; meanwhile, the dynamic response performance of each independent motor is improved (see Figure 1).
Figure 1 shows that the coordination controller integrates the traction torque output of each motor to compare it with the given total traction torque required. Based on the principle of minimized energy function, the realtime traction torque output state of each motor is redistributed to ensure that the sum of the traction torque of each motor is consistent with the total traction torque . At the same time, a sliding mode controller (SMC) is proposed to ensure that each directcurrent motor (DCM) has preferable dynamic response performance.
The mathematical model of permanent magnet brushless DCM can be expressed as [28]where and represent the resistance and inductance of the armature circuit, respectively; is the motor current; is the output angular velocity of motor gearbox j; is the voltage of the input port of the armature circuit; is the back electromotive force constant; is the gear ratio of the gearbox; and are the moments of inertia of the motor and gearbox; and are the viscous friction coefficients of the motor and the gearbox; is the motor torque constant; is the load torque; and is the output torque of the motor.
Set , , and ; the state equation of (2) can be expressed aswhere is the equivalent moment of inertia; is the equivalent viscous damping constant; and , , and ; besides, are the timevarying disturbances.
The load torque of the motor cannot be infinitely large; that is, . The axis end of the motor is mechanically connected to the wheelset through the gearbox. Moreover, the load torque can be considered slow, timevarying, or substantially constant; that is, . Therefore, the timevarying disturbance is bounded; that is, , where are all positive constants.
3. Multimotors Coordination Based on Energy Minimization
The multimotors coordination procedure is the dynamic change procedure of multimotors output traction torque, in which the energy function is composed of each traction torque reaching the minimum, which indicates the optimal dynamic adjustment process [29]. The energy function of traction torque can be expressed aswhere is the energy function that is composed of each traction torque; is the traction torque matrix; and is the weight matrix, and . The main diagonal element corresponds to the weight coefficient of each torque, and .
For the electric locomotive traction system that multiple motors provide power together, due to the different axle weights of the driving wheels, the traction torque provided by corresponding motor may have different effects on traction performance. The purpose of using the weight coefficient a_{j} is to reduce the traction fluctuation amplitude of motor that will cause greater degree of influence and, at the same time, to increase the traction adjustment amplitude of motor with a smaller degree of influence. So, the multimotors traction torque coordination method based on energy minimisation can be expressed as a constrained optimisation problem; that is,where is a positive constant.
Formula (5) can be solved by using the Karush–Kuhn–Tucker (KKT) conditions. The equality constraint of (5) can be expressed as ; the inequality constraints of it can be expressed as and , respectively. Moreover, the inequality constraint of (5) can be transferred to equality constraint by introducing slack variables,where and are the slack variables, respectively.
So the Lagrange equation is [30]where is the Lagrange multiplier and and are the KKT multipliers.
After the derivation of partial derivatives in (7),
The KKT condition indicates that the energy function must be satisfied when it reaches the minimum, and the derivations of partial derivatives in (8) are 0. Setting , that is,where the detailed explanation can be found in [31]; thus,
Also, with the percentage traction torque loss of the fourth motor, each traction torque matrix in (4) becomes , and is the changed traction torque matrix. At the same time, by changing the initial weight matrix , let it becomes the changed weight matrix . Finally, (10) can still be established.
Therefore, when the traction torque changes, the total traction is maintained at a constant under the constraint condition and the energy function value can be minimized by changing the weight coefficient matrix .
4. Sliding Mode Tracking Controller for Independent Motors
In a variable structure control system, the motion of the system can be divided into two stages. The first stage is the stage of reaching motion, that is, the approaching process in the sliding mode control in which the reaching condition ensures that the system state trajectory reaches the switching surface from any initial state within a finite time. The second stage is that the error converges to zero under the influence of control law. This control method makes the system state slide along the sliding surface by switching the control quantity, thereby making the system invariant when subjected to parameter perturbation and external disturbance; thus, it has the advantages of rapid dynamic response and strong robustness in multimotors system coordination control, aircraft motion control, and other fields [32]. For this purpose, a variable structure control was introduced for improving the dynamic response performance of each independent motor.
The tracking error between the output torque of motor j and the setting command is to be ; that is,
After the derivation of (11),
Formula (3) is substituted into (12); that is,
Integral sliding surface is defined aswhere is the positive constant to be designed and is a sign function.
The integral sliding mode control law is designed based on the nonlinear system shown in (3) and the integral sliding surface shown in (14), which can be expressed as follows:where and is the positive constant to be designed.
The controller contains unknown interference . Thus, the influence of can be offset by selecting a sufficiently large switching gain to converge the error to 0 in a finite time.
Proof. Lyapunov function is defined asAfter the derivation of ,Then,Set ; if is met, thenAccordingly,Meanwhile, the error gradually converges to 0 from the initial state, such that the output torque of each independent motor can be tracked to the reference state in a limited time.
This completes the proof.
5. Simulation Analysis
The traction system that consists of four motors was taken as the object, and the parameters are shown in Table 1.

The initial weight coefficient in (4) is selected as . The parameter to be designed on the switching function in (15) is selected as , and the parameter to be designed in (15) is selected as .
In order to reflect three processes including the acceleration, uniform speed, and braking, the corresponding reference command is set as
Four types of noise signals with high frequency, mutation, slowvarying, and uniformity are superimposed on the four traction motors, respectively, at 12 s to verify the dynamic response performance and noise interference suppression performance of the proposed method (see Figure 2).
(a)
(b)
(c)
(d)
The simulation is conducted on the system, that is, at 15 s, 100% loss and 30% loss of traction torque for a motor.
5.1. 100% Traction Torque Loss of a Motor
When the fourth motor loses 100% traction torque, the changing curve of the total traction torque and each motor traction torque are shown in Figure 3. The solid black line indicates the total given traction torque , the red dotted line indicates the sum of the traction torque of the four motors , and the lower coloured short line indicates each traction torque , respectively.
Figure 3 shows that, in the system acceleration phase of 0–10 s, the proportional coefficient assigned to the four motors is unchanged, and the sum of the output traction torques of the four motors is consistent with the total traction torque required. At 12 s, the total traction torque fluctuates with the maximum error of 0.12%. At 15 s, the traction torque output of the fourth motor is 0 due to the loss of traction, and the required traction torque is distributed to the other three motors. The maximum tracking error during the tracking of total quantity is 6%, and the tracking time is 0.3 s. Thus, the system can achieve total traction torque when a motor loses 100% traction torque.
To reflect the tracking performance of each independent motor more clearly, the traction torque change curve of the four motors in Figure 3 is shown in Figure 4. The solid black line indicates the assigned command signal, and the red dotted line indicates the output torque of each independent motor.
(a)
(b)
(c)
(d)
Figure 4 shows that, under the given noise signal interference at 12 s, each motor has different degrees of fluctuation, and the fluctuation quantities are 0.03%, 0.05%, 0.01%, and 0.1% with tracking time of 0.15, 0.3, 2.8, and 0.2 s, respectively. At 15 s, the traction torques of the three remaining motors show different degrees of overshoot during the transition to the torque required for redistribution due to the complete loss of traction of the fourth motor, and the overshoot is 10.7%, 11.6%, and 10.6% with stability time of 0.04 s, 0.15 s, and 0.15 s, respectively. Notably, after the fourth motor loses traction torque, the torque does not immediately become 0 and the fluctuation existed due to the influence of sliding mode. Thus, when noise exists, the motor torque output has high tracking performance using the proposed method.
5.2. 30% Traction Torque Loss of a Motor
When the fourth motor loses 30% of traction torque, the changing curve of the total traction torque and the each motor traction torque is shown in Figure 5. The solid black line indicates the given total traction , the red dotted line indicates the sum of the traction torque of the four motors , and the lower coloured short line indicates each traction torque .
Figure 5 shows that, at 15 s, the maximum tracking error during the tracking process of total traction torque is 1.3% due to the 30% traction loss of the fourth motor, and the tracking time is 0.1 s. Thus, the system can achieve total traction torque when a motor loses 30% traction torque.
The torque change curve of the four motors in Figure 5 is shown in Figure 6. The solid black line indicates the assigned command signal, and the red dotted line indicates the output torque of each independent motor.
(a)
(b)
(c)
(d)
Figure 6 shows that when the fourth motor loses 30% traction torque at 15 s, the output torque of the fourth motor is set to 70% of the original output torque, and the required traction torque is mainly compensated by the remaining motors. In the process of output torque coordination, the output torque fluctuations of the four motors are 3.21%, 3.24%, 3.29%, and 9.43% with tracking time of 0.04, 0.15, 0.2, and 0.2 s, respectively. Thus, when the motor loses part of the traction, the system still has preferable dynamic response performance.
6. Conclusions
The consistent total traction torque in the operation of electric locomotives is the premise of safe and stable operation of locomotives. A multimotors coordination control method for consistent total traction torques is proposed to ensure that the sum of multimotors output traction torque of electric locomotives is consistent with the total traction torque required. The main innovations of this study are as follows: firstly, a multimotors coordination control framework is constructed for consistent total traction torque. Secondly, a multimotors allocation control strategy, which provides an idea for optimising the multimotors dynamic adjustment process, is proposed based on minimisation of energy function. Thirdly, a traction torque control strategy based on sliding mode variable structure is constructed for each independent motor. Finally, the simulation environments under two different conditions are constructed to verify the proposed method. The results show that the proposed method shows good performance in terms of consistency of total traction torque, as well as dynamic response performance of an independent motor.
Data Availability
The survey data used to support the findings of this study are available from the corresponding author upon request.
Disclosure
The views expressed here are solely the responsibility of the authors.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported in part by the Natural Science Foundation of China (Grant no. 61773159), Excellent Youth Research Project of Education Department of Hunan Province (18B303), and Hunan Provincial Natural Science Foundation of China (Grant no. 2018JJ2100).
References
 Y. Liu, B. Zou, A. Ni, L. Gao, and C. Zhang, “Wheel/rail adhesion state identification of heavyhaul locomotive based on particle swarm optimization and kernel extreme learning machine,” Journal of Advanced Transportation, vol. 2020, Article ID 8136939, 6 pages, 2020. View at: Publisher Site  Google Scholar
 A. A. Zarifian and P. G. Kolpahchyan, “Computer modeling of electric locomotive as controlled electromechanical system,” Multibody System Dynamics, vol. 22, no. 4, pp. 425–436, 2009. View at: Publisher Site  Google Scholar
 J. U. Duncombe, “Infrared navigation—part I: an assessment of feasibility,” IEEE Transactions on Electron Devices, vol. 11, no. 1, pp. 34–39, 1959. View at: Google Scholar
 A. Hoffrichter, P. Fisher, J. Tutcher, S. Hillmansen, and C. Roberts, “Performance evaluation of the hydrogenpowered prototype locomotive “Hydrogen Pioneer” ,” Journal of Power Sources, vol. 250, pp. 120–127, 2014. View at: Publisher Site  Google Scholar
 K. Song, G. Konstantinou, W. Mingli, and P. Acuna, “Windowed SHEPWM of interleaved fourquadrant converters for resonance suppression in traction power supply systems,” IEEE Transactions on Power Electronics, vol. 32, no. 10, pp. 7870–7881, 2017. View at: Publisher Site  Google Scholar
 C. Zhang, Z. Lin, J. Liu, and J. He, “Consensusbased totalamount cooperative tracking control for multimotor locomotive traction system,” Journal of the Franklin Institute, vol. 356, no. 2, pp. 819–834, 2019. View at: Publisher Site  Google Scholar
 W.C. Cai, D.Y. Li, and Y.D. Song, “A novel approach for active adhesion control of highspeed trains under antiskid constraints,” IEEE Transactions on Intelligent Transportation Systems, vol. 16, no. 6, pp. 3213–3222, 2015. View at: Publisher Site  Google Scholar
 W. Zhang, J. Chen, X. Wu, and X. Jin, “Wheel/rail adhesion and analysis by using full scale roller rig,” Wear, vol. 253, no. 12, pp. 82–88, 2002. View at: Publisher Site  Google Scholar
 R. Galas, M. Omasta, I. Krupka, and M. Hartl, “Laboratory investigation of ability of oilbased friction modifiers to control adhesion at wheelrail interface,” Wear, vol. 368369, pp. 230–238, 2016. View at: Publisher Site  Google Scholar
 J. H. Feng, “Research on the permanent magnet synchronous motor drive system for highspeed EMUs,” Journal of Electric Drive for Locomotives, no. 4, pp. 1–5, 2016, in Chinese. View at: Google Scholar
 N. R. Abjadi, J. Askari, J. Soltani, and G. R. Arab Markadeh, “Nonlinear slidingmode control of a multimotor webwinding system without tension sensor,” IET Control Theory & Applications, vol. 3, no. 4, pp. 419–427, 2009. View at: Publisher Site  Google Scholar
 C. Zhang, M. Niu, J. He, and K. Zhao, “Robust synchronous control of multimotor integrated with artificial potential field and sliding mode variable structure,” IEEE Access, vol. 5, pp. 197–207, 2017. View at: Publisher Site  Google Scholar
 M.T. Yan, M.H. Lee, and P.L. Yen, “Theory and application of a combined selftuning adaptive control and crosscoupling control in a retrofit milling machine,” Mechatronics, vol. 15, no. 2, pp. 193–211, 2005. View at: Publisher Site  Google Scholar
 H. Hou, X. Nian, H. Xiong, Z. Wang, and Z. Peng, “Robust decentralized coordinated control of a multimotor webwinding system,” IEEE Transactions on Control Systems Technology, vol. 24, no. 4, pp. 1495–1503, 2016. View at: Publisher Site  Google Scholar
 J. Liao, Z. Chen, and B. Yao, “Performanceoriented coordinated adaptive robust control for fourwheel independently driven skid steer mobile robot,” IEEE Access, vol. 5, pp. 19048–19057, 2017. View at: Publisher Site  Google Scholar
 C. Zhang, H. Wu, J. He, and C. Xu, “Consensus tracking for multimotor system via observer based variable structure approach,” Journal of the Franklin Institute, vol. 352, no. 8, pp. 3366–3377, 2015. View at: Publisher Site  Google Scholar
 T. Shi, H. Liu, Q. Geng, and C. Xia, “Improved relative coupling control structure for multimotor speed synchronous driving system,” IET Electric Power Applications, vol. 10, no. 6, pp. 451–457, 2016. View at: Publisher Site  Google Scholar
 C. Zhang, Z. Lin, S. X. Yang, and J. He, “TotalAmount synchronous control based on terminal slidingmode control,” IEEE Access, vol. 5, pp. 5436–5444, 2017. View at: Publisher Site  Google Scholar
 Z. Zhang, Y. Zhao, W. Qiao, and L. Qu, “A spacevectormodulated sensorless directtorque control for directdrive PMSG wind turbines,” IEEE Transactions on Industry Applications, vol. 50, no. 4, pp. 2231–2341, 2014. View at: Publisher Site  Google Scholar
 C. Zhang, H. Liao, X. Li, J. Sun, and J. He, “Fault reconstruction based on sliding mode observer for current sensors of PMSM,” Journal of Sensors, vol. 2016, Article ID 9307560, 9 pages, 2016. View at: Publisher Site  Google Scholar
 W. Xu and R. D. Lorenz, “Dynamic loss minimization using improved deadbeatdirect torque and flux control for interior permanentmagnet synchronous machines,” IEEE Transactions on Industry Applications, vol. 50, no. 2, pp. 1053–1065, 2014. View at: Publisher Site  Google Scholar
 Y. S. Choi, H. C. Han, and J. W. Jung, “Feedback linearization direct torque control with reduced torque and flux ripples for IPMSM drives,” IEEE Transactions on Power Electronics, vol. 31, no. 5, pp. 3728–3737, 2015. View at: Publisher Site  Google Scholar
 F. A. Rahimi, M. G. Lauby, J. N. Wrubel, and K. L. Lee, “Evaluation of the transient energy function method for online dynamic security analysis,” IEEE Transactions on Power Systems, vol. 8, no. 2, pp. 497–507, 2002. View at: Publisher Site  Google Scholar
 M. Levitt, M. Hirshberg, R. Sharon, and V. Daggett, “Potential energy function and parameters for simulations of the molecular dynamics of proteins and nucleic acids in solution,” Computer Physics Communications, vol. 91, no. 1–3, pp. 215–231, 1995. View at: Publisher Site  Google Scholar
 R. Yousefian, A. Sahami, and S. Kamalasadan, “Hybrid transient energy functionbased realtime optimal widearea damping controller,” IEEE Transactions on Industry Applications, vol. 53, no. 2, pp. 1506–1516, 2017. View at: Publisher Site  Google Scholar
 H. Shao, H. Jiang, H. Zhang, and T. Liang, “Electric locomotive bearing fault diagnosis using a novel convolutional deep belief network,” IEEE Transactions on Industrial Electronics, vol. 65, no. 3, pp. 2727–2736, 2017. View at: Google Scholar
 C. Mayet, J. Pouget, A. Bouscayrol, and W. Lhomme, “Influence of an energy storage system on the energy consumption of a dieselelectric locomotive,” IEEE Transactions on Vehicular Technology, vol. 63, no. 3, pp. 1032–1040, 2014. View at: Publisher Site  Google Scholar
 F. BeltranCarbajal, A. ValderrabanoGonzalez, J. C. RosasCaro, and A. FavelaContreras, “An asymptotic differentiation approach of signals in velocity tracking control of DC motors,” Electric Power Systems Research, vol. 122, pp. 218–223, 2015. View at: Google Scholar
 Y. D. Song, J. N. Anderson, and H. Y. Lai, “Robust control of multirobot systems: the generalized energy accumulation approach,” International Journal of Robust & Nonlinear Control, vol. 4, no. 5, pp. 673–696, 2010. View at: Publisher Site  Google Scholar
 C. Mahapatra, Z. Sheng, P. Kamalinejad, V. C. M. Leung, and S. Mirabbasi, “Optimal power control in green wireless sensor networks with wireless energy harvesting, wakeup radio and transmission control,” IEEE Access, vol. 5, pp. 501–518, 2016. View at: Publisher Site  Google Scholar
 M. Li, “Generalized Lagrange multiplier method and KKT conditions with an application to distributed optimization,” IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 66, no. 2, pp. 252–256, 2019. View at: Publisher Site  Google Scholar
 S. Laghrouche, F. Plestan, and A. Glumineau, “Higher order sliding mode control based on integral sliding mode,” Automatica, vol. 43, no. 3, pp. 531–537, 2007. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Jianhua Liu et al. 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.