Next Article in Journal
Study of Surface Roughness Effect on a Bluff Body—The Formation of Asymmetric Separation Bubbles
Previous Article in Journal
Development of Steelmaking Slag Based Solid Media Heat Storage for Solar Power Tower Using Air as Heat Transfer Fluid: The Results of the Project REslag
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Energy-Efficient Speed Profile Optimization and Sliding Mode Speed Tracking for Metros

School of Electrical Engineering, Southwest Jiaotong University, Chengdu 610031, China
*
Author to whom correspondence should be addressed.
Energies 2020, 13(22), 6093; https://doi.org/10.3390/en13226093
Submission received: 26 October 2020 / Revised: 18 November 2020 / Accepted: 18 November 2020 / Published: 20 November 2020
(This article belongs to the Section A1: Smart Grids and Microgrids)

Abstract

:
Nowadays, most metro vehicles are equipped with an automatic train operation (ATO) system, and the speed control method, combining cruise speed planning and proportional-integral-derivative (PID) control, is widely used. The automation is achieved, and the energy-efficient can be improved. This paper presents an improved artificial bee colony algorithm for speed profile optimization with coast mode and an adaptive terminal sliding mode method for speed tracking. Specifically, a multi-objective optimization model is established, which considers energy consumption, comfortableness, and punctuality. Then, a novel artificial bee colony algorithm named regional reinforcement artificial bee colony (RR-ABC) is designed, to search the optimal speed profile with coast mode, in which some improvements are made to speed up convergence and to avoid local optimal solutions. For speed-tracking control, the adaptive terminal sliding mode controller (ATSMC) is used to improve the speed error, robustness, and energy saving. In addition, a disturbance observer (DOB) is designed to improve the anti-interference ability of the system and further improve the robustness and anti-disturbance, which are also conducive to speed error and energy saving. Finally, the line and train data of the Qingdao Metro Line 6 are used for simulation, which proves the effectiveness of the study. Specific to the energy saving rate, and compared with normal algorithms, RR-ABC with coast mode is approximately 9.55%, and ATSMC+DOB is 7.58%.

1. Introduction

In the modern metro system, the quality of operation, including energy saving, punctuality, and comfortableness, is the core factor for operation [1]. With the popularization of the Automatic Train Operation (ATO) system, the automation of train operation has been achieved. However, the performance of the algorithm in ATO system is not good enough, and it can be further improved [2]. Inaccuracy and high energy consumption seriously affect the quality of train operation.
The frequent switching of operation mode is needed for the metro system because the start–stop is frequent. The ATO system needs to refer to the planned speed profile as a reference to issue control commands, such as acceleration, braking, and coasting [3]. In practice, the PID controller is widely used to realize speed control, and the fluctuation of control force occurs when the coasting mode is needed [4]. This problem eventually leads to larger speed errors and poor energy saving effects. To improve the tracking quality, a large number of papers are published to try to find a perfect substitute of the PID controller. Predictive controller is presented to track the speed profile. In Reference [5], an integral sliding mode control (ISMC) is used to design distributed adaptive control strategies. In Reference [6], a tracking controller is formulated to resolve errors caused by unmolded effects. In References [7,8], a robust adaptive controller is proposed to adapt unknown system parameters and nonlinear traction/braking notches. Wang et al. proposed an approach to design a distributed and fault-tolerant controller which is applied to develop braking control schemes for high-speed trains subject to traction and braking failures [9]. In Reference [10], to deal with input saturation and unknown disturbance in high-speed trains, a neural adaptive fault-tolerant controller is designed. In all the abovementioned studies, either the convergence cannot be proved by Lyapunov theory or the convergence speed is not good enough. In addition, large speed errors in those controllers also result in large energy consumption. Therefore, they are not suitable for metro trains that require fast convergence and high-precision tracking.
Recently, sliding mode control is becoming a research hotspot. Because it is easy to design and did not require accurate parameter identification, it has been used in actual physical systems, such as power systems and industrial robot manipulators [11,12,13]. In addition, the controller with a nonlinear sliding surface has been proven to converge in a limited time [14]. Especially for the adaptive terminal sliding mode controller (ATSMC), in addition to the advantages of fast convergence, it can also be adaptive to the internal parameter disturbance and external interference of the system [15]. Inspired by the above articles, we introduce the additional resistance caused by slope and the basic running resistance of the train into the sliding mode controller, in the form of disturbance, and an adaptive terminal sliding mode controller for train-speed-tracking control (ATSMC) is designed in this paper. In order to further suppress the disturbance and enhance the control accuracy, which are important to energy saving, comfortableness, and punctuality, this paper introduces a disturbance observer (DOB) on the basis of ATSMC to track the speed profile; finally, ATSMC+DOB is designed, and the convergence is proven by Lyapunov theory.
As the premise of speed-tracking control, the reference speed profile needs to be generated by the speed optimization algorithm. Similar to the speed-tracking controller, the speed optimization algorithm used in engineering is not good enough either. As a widely used cruise optimization algorithm (cruising mode), it is easy to calculate, in theory. However, it is not only poor in energy-saving, but also difficult to achieve precise control in industry. Therefore, the coast optimization algorithm is adopted by a large number of studies, to reduce energy consumption. It provides a reference speed curve that is convenient for the controller to track. The energy saving of coast optimization has been proven with Pontryagin maximum principle by some studies [16,17,18]. In References [19,20], based on coast optimization rules, the genetic algorithm is used to search the optimal coast point. Chuang et al. used an optimization method based on neural networks to calculate the optimal solution [21]. Khmelnitsky et al. used the maximum–minimum ant system (MMAS) in the ant colony optimization algorithm (ACO) to optimize the train’s running speed trajectory [22]. In Reference [23], a dual heuristic programming method is used to optimize the coast time and dwell time at the same time. In Reference [24], a kind of dual speed profile optimization is designed, to save energy. Cheng et al. proposed an improved ant colony optimization method to study train operation [25]. Rocha et al. used the coast mode to design a n approach for real time train optimization [26]. In References [27,28], multi-objective optimization is considered to study speed profile. Shang et al. proposed an online energy-saving driving strategy for metro [29]. In Reference [30], a switched nonlinear model predictive controller (NMPC) is designed for collaborative ecodrive control of railway vehicles. Furthermore, to use more renewable energy, a dissension-based adaptive law (DAL) is proposed, to adjust the parameters of the NMPC cost. In Reference [31], a kind of online distributed cooperative model predictive controller is proposed, to generate a speed profile for multiple high-speed train movement. In Reference [32], a kind of energy-efficient driving strategy based on transmission losses is proposed.
In the above research, most optimization models focus on energy saving. The lack of comfortableness and punctuality is fatal to the metro. Therefore, a multi-objective optimization model is established in this paper. In this model, in addition to energy consumption, comfortableness and punctuality are taken into consideration as penalty factors to comprehensively optimize operation quality.
In addition, for the intelligent evolutionary algorithm mentioned above, there are always two problems to be solved. Firstly, maintaining good population diversity to avoid a local optimal solution while evolving rapidly is always a problem [33]. Secondly, the initial value of the population has a significant impact on the optimization result, so how to obtain a better initial value is worth studying. To solve these problems, the artificial bee colony (ABC) algorithm, which shows good global search ability, may be used in speed profile optimization [34]. However, though good at global search, the ABC algorithm is poor in local search. When it is close to the optimal solution, the search efficiency drops sharply in a local region. To make up for this deficiency, a regional reinforcement artificial bee colony (RR-ABC) algorithm is proposed, to optimize the speed profile. By using chaotic mirror initialization, local search strategy, and the elite group guidance strategy, the RR-ABC algorithm maintains better population diversity and enhances search efficiency.
The specific structure is as follows. Section 2 shows the dynamic model and some special restrictions of metro. Section 3 establishes a speed profile optimization calculation model which considers energy consumption, comfortableness, and punctuality. In Section 4, the RR-ABC algorithm is proposed, to optimize the speed profile. In Section 5, to track the speed, we propose the ATSMC+DOB controller. In Section 6, the validity of the method is verified based on the practical data of Qingdao Metro Line 6. Finally, we conclude the study in Section 7.

2. Dynamic Characteristics of Metro Trains

There are two common dynamic models: multi-mass point model and single-particle model. For the multi-particle model, each vehicle of the train is treated as an independent mass point. As a result, the force deviation caused by train length and the power transmission between different vehicles are considered in the model. However, this model has a complicated modeling process and is also difficult to calculate [35]. Due to the lighter weight and shorter length of the metros, this accuracy is not necessary. Compared with multi-particle model, the single-particle model regards the entire train as a mass point which has a simple modeling process. A large number of related studies have adopted the single-particle model [36,37].
Train single-particle dynamics model in Reference [36] is used in this paper. The schematic diagram of the train force is shown in Figure 1. In this model, m is the train mass, which is considered as a constant value throughout the journey. About the decomposition of the forces, denotes the traction force, F b denotes the braking force, R c denotes the additional resistance of the line that mainly comes from the slope, and R b denotes the basic running resistance of the train, mainly from the track and air. The dynamics model can be written as follows:
x ˙ = v m 1 + γ v ˙ = F t v F b ( v ) R b v R c v
where v denotes the speed, and x is the position; F t v is the force of traction, F b v is the force of braking, R c v is the additional resistance, and R b v is the basic resistance; and γ is the rotary mass coefficient which is usually a given constant.
About the rotary mass coefficient, λ , it is a coefficient to compensate the moment of inertia. The train kinetic energy has two parts, one is the horizontal movement energy of the train, and the other is the moment of inertia of the wheels. Therefore, when only the horizontal movement is concerned, the moment of inertia needs to be compensated. This value can be obtained by dividing the converted weight of the moment of inertia by the total weight. Since this value is difficult to determine, empirical values are generally used in simulation. Therefore, m 1 + γ is the total weight of the train after correction.
About the traction, the maximum that the train can provide is as follows:
F t max = 203 0 v 51.5 0.002 v 3 + 0.49 v 2 42.12 v + 1343 51.5 < v 90
F t v = μ t F t max v , 0 < μ t < 1
where F t ( v ) is the real force, and μ t control factor.
About the braking, similar to traction, the calculation formula can be written as follows:
F b max = 166 0 v 77 0.134 v 2 25.07 v + 1300 77 < v 90
F b v = μ b F b max v , 0 < μ b < 1
About the basic running resistance, it can be written as follows:
R b v = r 1 + r 2 v + r 3 v 2
where r 1 , r 2 , and r 3 are empirical parameters determined by experiment.
The additional resistance can be written as follows:
R c v = w i + w r + w s
where w i is the additional resistance of the slope, w r is the additional resistance in the curve, and w s is the additional resistance in the tunnel. Generally, the resistance in the curve and tunnel are not considered in the study and the slope resistance is given as:
w i = m g sin θ
where m is the train mass, g is the acceleration of gravity, and θ is the slope value.
The external characteristic curves (max force) of train traction and braking force are shown in Figure 2, and Equations (2), (4), and (6) are examined by the actual experiment from Qingdao Metro Line 6.

3. Fitness Calculation Function about Speed Profile Optimization

This part establishes the optimization model of the profile, including the energy consumption calculation model, penalty function about punctuality, and comfortableness. These factors will be discussed separately below.

3.1. Energy Consumption Calculation Function

As shown in Figure 3, because of the speed limit, the operation of the metro on the entire line is divided into N intervals, where v i 0 i = 1 , 2 , , N represents the starting speed of the train in the interval, i , and it also represents the ending speed in the interval, i 1 ; and v 1 0 and v N + 1 0 represent the starting speed and the arrival speed, respectively; obviously, v 1 0 = v N + 1 0 = 0 .
As shown in Figure 4, inside the interval, i , we split the single one into n subintervals, which have same running time, Δ t . In subinterval n , the description of the symbols is shown in Table 1.
Therefore, the length, L i , and the time, T i , of the interval, i , can be descripted as follows:
L i = j = 1 n X j T i = n Δ t
The total length, S t o t a l , and the total time, T t o t a l , of all intervals are as follows:
S t o t a l = i = 1 N L i T t o t a l = i = 1 N T i
The traction energy consumption, E j , of subinterval can be calculated as follows:
E j = F j X j
In Equation (11), the length of subinterval X j can be calculated as follows:
X j = v i j 1 Δ t + 0.5 a j Δ t 2
Thus, the traction energy consumption of the interval, i , can be calculated as follows:
E = i = 1 N E i E i = j = 1 n E j

3.2. Punctuality Penalty Function

Punctuality is a key factor in evaluating the quality of metro operation. The current trains run on time according to the timetable issued by the ATS system. The running time between stations is already specified, and the error between the actual running time and the specified running time needs to be limited within a required error. According to the actual operation requirements of the metro, the maximum allowable time error is generally 6 s.
Define the punctuality judgment function of the train as follows:
T e r r o r = T D T R
where T e r r o r is the arrival time error, T D is the planned running time, and T R is the actual running time.
As shown in Figure 5, drawing on the Gaussian model, the punctuality penalty function is established as follows:
R t = e T e r r o r T E ρ 2
where R t is the time penalty factor, T E is the maximum allowable time error (set 6 s here), and ρ is the sensitivity control coefficient, which can be used to control penalty level.

3.3. Comfortableness Penalty Function

Comfortableness is another important indicator of subway operation. Generally, the rate of change of acceleration is used to evaluate comfortableness.
Define the acceleration change rate calculation model as follows:
J e r k = d a d t
where J e r k is the acceleration rate, a is acceleration, and t is running time.
As shown in Figure 6, similar to punctuality penalty function, the comfortableness penalty function is established as follows:
R c o m = e J e r k J d ρ 2
where R c o m is comfortableness penalty factor, J d is maximum allowable acceleration rate of change, and ρ is the sensitivity control coefficient.
As shown in Table 2, according to international standard IS02631 [38], set the value of J d to 1. When the acceleration change rate is larger, the penalty for comfort is greater, so that the non-conforming values in the result can be eliminated.

3.4. Fitness Calculation Function Based on Double Penalty Mechanism

As considered in Section 3.1, Section 3.2 and Section 3.3, the fitness calculation function based on double penalty mechanism is summarized as follows:
M i n   f i t n e s s = R t + R c o m + E

4. Speed Profile Optimization with RR-ABC

4.1. Emergency Braking Intervention (EBI) Curve and Warning Speed Limit

Resulted from some safety reasons, the metro line should consider various speed-limit protections.
Firstly, the metro line has a fixed speed limit from the line. In addition, the metro also has a security protection mechanism to ensure that emergency braking is feasible in any position that the train cannot cross the fixed speed limit; this mechanism is called emergency braking intervention (EBI). To ensure that the metro does not exceed the fixed speed limit during emergency braking under any circumstances, the EBI protection considers such an extreme case:
  • At first, the train is in the process of maximum acceleration. After the train issued an emergency braking command, the train still maintained the maximum acceleration during the transmission delay and traction cutoff delay, because the traction cannot be cut off instantly and transmission delay exists objectively at any moment.
  • Since the braking system needs some time to receive instructions and gradually generate braking force, it usually takes some time to reach the maximum force from 0. During this period, the train is considered to keep coasting.
Therefore, as shown in Figure 7, the emergency braking process is divided into three stages:
Stage1: In this stage, the train accelerates with the maximum traction acceleration, and the acceleration time is the sum of the traction cutoff delay and the braking system action delay.
Stage2: After the train traction is cut off, the emergency braking force needs a certain time to reach the nominal value, and the train is considered to coast during this period.
Stage3: The train follows the emergency braking curve to brake.
Finally, as a result of the single-mass model, we ignore the size and length of the train. However, the ignorance will bring security risks to trains, so we design the warning speed limit based on the EBI curve.
As shown in Figure 8, when the train moves from the low-speed-limit zone to the high-speed-limit zone, the train needs to reserve a buffer equal to the length of the train. When the train moves from the high-speed-limit zone to the low-speed-limit zone, the same buffer is required.
According to the above principle, taking Qingdao Metro Line 6 grabbing Zhuomashan station to Heluobu station as an example, these kinds of speed limit curves are obtained as shown in Figure 9.

4.2. Energy-Saving Strategies and Coast Choices

According to actual operating conditions and theory, coast travel is an effective strategy for train energy saving. The coast strategy we designed is shown in the Figure 10, in which two different coast intervals are set.

4.2.1. Coast Interval x 1

  • Features:
    (1)
    Driving from high speed limit zone to low speed limit zone
    (2)
    Maintain constant speed at the next low speed limit
  • Strategy selection:
From the coast point, let the train maintain the coast condition until point S 1 , then we get the speed of S 1 point v S 1 .
(1)
If v S 1 > v 3 , then brake to v 3 ;
(2)
If v S 1 < v 3 , then when v = v 3 , let the train drive at a constant speed v 3 .

4.2.2. Coast Interval x 2

  • Features:
    (1)
    Driving from high speed limit zone to low speed limit zone;
    (2)
    Braking in the next low speed limit.
  • Strategy selection:
From the coast point, let the train maintain the coast condition until point S 2 , then we get the speed of S 2 point v S 2 .
(1)
If v S 2 > v 4 , then brake to v 4 ;
(2)
If v S 2 v 4 , then keep coasting to maintain the coast condition until it contacts the braking curve.

4.3. Regional Reinforcement Artificial Bee Colony (RR-ABC) Algorithm

The ABC algorithm simulates the process of the bee colony searching for the honey source with the highest honey content. It has few control parameters, and it has good global convergence. However, when the local search is approached to converge to the optimal solution, the search efficiency drops sharply [39,40]. With the deepening of evolution, population diversity is seriously lacking, and the search efficiency is significantly reduced at the end of evolution. To balance the global and local search capabilities and to improve efficiency of the ABC, a more efficient ABC algorithm is proposed called RR-ABC algorithm. The principles and steps of the RR-ABC algorithm are as follows:
(1) 
Model initialization: Chaos mirror initialization
Because of the randomness of the chaotic map and the sensitivity to the initial conditions, we extracted the decision space information to grow the diversity of the group. In this part, the sine formula is used to iteratively generate chaotic variables. The chaotic initialization Equation is as follows:
c h 1 = r a n d 0 , 1 c h k + 1 = sin π c h k x i j = x i j min + c h k j x i j max x i j min
where c h 1 is random number between 0 and 1; k is number of iterations; k = 1 , 2 , , D 1 , and j = 1 , 2 , , D ; x i j is the j-dimensional parameter of the i-th honey source; and x i j max , x i j min are the top and bottom limits of the source.
Use Equation (19) to perform the chaos initialization operation to obtain the i th honey source x i = ( x i 1 , x i 2 , , x i D ) , then mirror it to get the mirror honey source m i r i = x i D , , x i 2 , x i 1 , then compare the fitness value of the initial solution one by one, and select the better individual as the initial population.
(2) 
Employed bee stage: Evolutionary dimension adjustment strategy
In the stage of employed bee and onlooker bee, they need to update the honey source through Equation (20); the original honey source has only one dimension for each change. Obviously, this method is extremely inefficient. To solve the problem, we use a new parameter W , to control the dimensions that need to be changed in each evolution. The new search Equation is as follows:
x i w = x i w + R i w × x i w x k w
where w = 1 , 2 , , W , parameter W controls the number of dimensions of each individual change, and W = 1 , 2 , , D .
Every time an employed bee searches for a new honey source, it will estimate the potential of the source (calculate the fitness value of the objective function):
f i t = 1 1 + f i
where f i is the objective function value of the i-th honey source; f i t i is the fitness value of i-th honey source. The employed bee chooses the better honey source by comparing the fitness value of X i and X i (Greedy choice). The larger the fitness value, the smaller the corresponding objective function value, indicating that the honey source is also better, and with a higher probability of getting better honey sources around it.
X i = X i         f i t ( X i ) > f i t ( X i ) X i         f i t ( X i ) f i t ( X i )
(3) 
Onlooker bees stage: Search radius adjustment strategy
As mentioned in “Model of employed bees”, the employed bee will search within a certain radius of the initial solution. In the early stage of evolution, because the honey source x i is far from the real solution, the use of a larger search radius can accelerate the convergence rate; when the honey source evolves to the final stage, it is closer to the real solution, and reducing the search radius can improve the search accuracy. Therefore, using a variable search radius and adjusting the size of the search radius can significantly speed up the convergence rate.
In this paper, the strategy of changing the search radius with the individual threshold, that is, as the threshold, L i m i t i (the digit of times that honey source x i was not selected), increases, the search radius continues to decrease. The Equation is as follows:
r i w = ϕ + cos π L i m i t i w 2 L i m i t
where r i w is search radius honey source x i w , ϕ is the basic value of search radius, L i m i t i w is the threshold of the honey source x i w , and L i m i t is the maximum allowed threshold. The new search Equation is as follows:
x i w = x i w + r i w × x i w x k w
Remark 1.
For the traditional ABC algorithm, with the deepening of evolution, more nectar sources are continuously explored on a global scale. However, due to the unknown quality of the preliminary search results, each explored nectar needs to be searched more deeply, which seriously affects the search efficiency in the later stages of evolution, even if it can avoid the generation of local optimal solutions. In addition, for classic biomimetic algorithms, the quality of the initial value usually greatly affects the evolution process, while a simple random process can easily produce low-quality initial values. Through the flowchart about the comparison of classic method and reinforcement method (Figure 11), we can see that dynamically adjusted evolutionary dimensions and search radius are adopted. In the later stages of the search, fewer high-quality nectar sources will be paid attention, and a better efficiency can be got.

4.4. Speed Profile Optimization Model

Based on the above discussion, the speed curve optimization model is summarized as follows:
M i n           f i t n e s s = E + R t + R c o m
s.t.
x ˙ = v m 1 + γ v ˙ = F t v F b v R b v R c v 0 v v s a f e t y

5. Adaptive Terminal Sliding Mode Controller

This section needs to design a speed tracking controller to track the speed profile optimized before. As mentioned before, compared with PID controller, which is widely used in modern industrial control, the sliding mode controller is quick to respond and completely robust to system parameter changes and external interference. However, sliding mode control will produce chattering when dealing with uncertain factors. This kind of high-frequency chattering will not only affect the accuracy of the control, but also cause the system to oscillate or become unstable.
Adaptive terminal sliding mode control not only has strong robustness, but also lets the system state converge to the ideal trajectory within a limited time. In addition, it also has a strong parameter adaptive processing function, to ensure that the system will not have unnecessary continuous switching when the parameters are uncertain. In order to design a speed-tracking controller with excellent comprehensive performance, this paper designs an adaptive terminal sliding mode controller with disturbance observer.

5.1. Dynamics Model of Speed Tracking

This part models and transforms the known fixed time delay. Here, the time delay is compensated by the sliding mode control coefficient, so that the system can be used as a non-delay system to design the sliding mode control rate in the subsequent chapters.
As mentioned in Section 2, the ideal dynamics model is Equation (1); for convenience of expression, the model with speed v as its independent variable is rewritten as the model with time, t , as its independent variable. Furthermore, the traction force factor, F t v , and braking force factor, F b v , are replaced by the actual generated traction/braking force, F u t , and the model is summarized as follows:
x ˙ = v m 1 + γ v ˙ = F u t R b t R c t
In the actual control system, due to response delay, the model needs to be written as a first-order lag model:
a ˙ c t = 1 τ a c t + 1 τ a t t T d
where a c is the control acceleration, which is produced by the controller; a t is the target acceleration; τ is the response time constant; and T d is controlling transmission delay.
For subsequent kinetic description, we define the following Equation:
a c = F u t m 1 + γ
Therefore, the acceleration formed by the basic resistance of the train and the additional resistance can be expressed as follows:
d t = W t a b v c v 2 W t = R a t m 1 + γ a = r 1 1 + γ b = r 2 1 + γ c = r 3 1 + γ
Furthermore, according to Equation (30), we can get the following:
v ˙ t = a c t + d t
where d t is the resistance acceleration.
The existing metro controllers are mostly stepless control modes, and the control system model is shown in Figure 12.
In order to facilitate the design of subsequent controllers, pade approximation is used to rewrite the delay parameters T d , in Equation (28):
e T d s s + λ s + λ
where λ is a constant related to delay, T d .
Laplace transform of Equation (28):
a c = τ λ a ¨ c λ τ + 1 λ a ˙ c 1 λ a ˙ t + a t
In order to facilitate expression, Equation (33) is rewritten as follows:
a c = q 1 a ¨ c q 2 a ˙ c q 3 a ˙ t + a t q 1 = τ λ q 2 = λ τ + 1 λ q 3 = 1 λ
As we can see, if the response time constant τ and the delay T d are known and kept unchanged, then q 1 , q 2 , and q 3 are all constants. In the following content q 1 , q 2 , and q 3 are used to design the sliding mode control rate, so the delay will be compensated. The time delay is compensated to the sliding mode control rate coefficient through the abovementioned method, and the effect of delay is eliminated in this way.

5.2. Design of Sliding Mode Terminal Controller

We define the state error as follows:
e 1 = x x r e 2 = v v r
where e 1 is the train position error, e 2 is the train speed error, x r is the reference position of optimized profile, and v r is the reference speed of optimized profile.
Differentiate Equation (35) and substitute it into Equation (31), and the resulting error state space Equation is as follows:
e ˙ 1 = e 2 e ˙ 2 = a c + d v ˙ r
To improve the tracking accuracy, the designed sliding mode surface needs to introduce the train position error e 1 , and speed error e 2 , to ensure fast error convergence. The designed terminal sliding mode function is as follows:
s m t = β e 1 + e p / q 2
where β > 0 ; and p and q are positive odds, 1 < p q < 2 .
Design the sliding mode controller below.
Differentiate the sliding mode function of (37) as follows:
s ˙ m t = β e 2 + p q a c + d v ˙ r e 2 p / q 1
According to Equations (30) and (31), we can obtain the following:
a c = v ˙ + a + b v + c v 2 + W t a ˙ c = v ¨ + b v ˙ + 2 c v v ˙ + W ˙ t a ¨ c = v + b v ¨ + 2 c v v ¨ + 2 c v ˙ 2 + W ¨ t
Substitute Equation (39) into Equation (34), since q 1 b q 2 , ignore q 1 b v ¨ , and get the following:
a c = q 1 v + 2 c v ˙ 2 + 2 c v v ¨ + W t q 2 v ¨ + b v ˙ + 2 c v v ˙ + W ˙ t q 3 a ˙ t + a t
By combining Equations (39) and (40), and bringing Equation (40) into the first derivative (38) of the sliding mode function, we can obtain the following:
s ˙ m t = β e 2 + p q e 2 p / q 1 a b v c v 2 W t q 1 W ¨ t q 2 W ˙ t v ˙ r                           + p q e 2 p / q 1 q 1 v 2 q 1 c v ˙ 2 2 q 1 c v v ¨ q 2 v ¨ q 2 b v ˙ 2 q 2 c v v ˙ q 3 a ˙ t + a t
Based on the principle of sliding mode controller, the sliding mode control input has the following form:
u t = u e q + u n
where u e q is the equivalent control items of the system, and u n is the nonlinear switching control items of the system.
Make Equation (41) approach zero, to get the following:
u e q = q p β e 2 2 p / q + W t + q 1 W ¨ t + q 2 W ˙ t + a + b v + c v 2                   + v ˙ r + q 1 v + 2 q 1 c v ˙ 2 + 2 q 1 c v v ¨ + q 2 v ¨ + q 2 b v ˙ + 2 q 2 c v v ˙ + q 3 a ˙ t
Make D t = W t q 1 W ¨ t q 2 W ˙ t a b v c v 2 , since W t is a smooth resistance and has a boundary; then its first and second derivatives are bounded, and the basic running resistance of the train f b t is bounded, so there is D D max , and D max is the boundary. According to the design principle of sliding mode control, the nonlinear switching term is designed as u n = k sgn s m , and then the terminal sliding mode control input is as follows:
u t = q p β e 2 2 p / q + q 1 v + q 2 v ¨ + q 3 a ˙ t + 2 q 1 c v ˙ 2 + 2 q 1 c v v ¨                         + q 2 b v ˙ + 2 q 2 c v v ˙ D + v ˙ r k sgn s m
where k > 0 is the control gain.
In addition, since the resistance items of the train are unknown in the actual operation of the train, that is D is unknown, but D is bounded, the actual controller should be changed to the following:
u t = q p β e 2 2 p / q + q 1 v + q 2 v ¨ + q 3 a ˙ t + q 4 v ˙ 2                         + 2 q 4 v v ¨ + q 5 v ˙ + 2 q 6 v v ˙ + v ˙ r k sgn s m
where q 4 = q 1 c ; q 5 = q 2 b ; q 6 = q 2 c ; and k > D max .
Select the Lyapunov function of the following form:
V t , x = 1 2 s m 2 x
Let k 0 = k D max , derivate Equation (46), and bring the controller formula, Equation (45), into it:
V ˙ = s m s ˙ m = s m k 0 p q sgn ( s m ) p q D max sgn ( s m ) D e 2 p / q 1                               k 0 p q e 2 p / q 1 s m < 0  
The system is stable, but in order to reduce the control burden of switching items and reduce the possibility of chattering, a parameter adaptation mechanism is introduced.
First modify the controller to the following form:
u t = q p β e 2 2 p / q + q ^ 1 v + q ^ 2 v ¨ + q ^ 3 a ˙ t + 2 q ^ 4 v ˙ 2                         + 2 q ^ 4 v v ¨ + q ^ 5 v ˙ + 2 q ^ 6 v v ˙ + a ^ + b ^ v + c ^ v 2 + v ˙ r k sgn s m
where q ^ 1 , q ^ 2 , q ^ 3 , q ^ 4 , q ^ 5 , q ^ 6 , a ^ , b ^ , c ^ are the estimated values of parameters q 1 , q 2 , q 3 , q 4 , q 5 , q 6 a , b , c , respectively.
Construct the following Lyapunov function:
V ˙ = s m s ˙ m + i = 1 6 1 λ i Δ q i Δ q ˙ i + 1 λ 7 Δ a Δ a ˙ + 1 λ 8 Δ b Δ b ˙ + 1 λ 9 Δ c Δ c ˙
where
Δ a = a ^ a , Δ a ˙ = a ^ ˙ Δ b = b ^ b , Δ b ˙ = b ^ ˙ Δ c = c ^ c , Δ c ˙ = c ^ ˙ Δ q i = q ^ i q i , Δ q ˙ i = q ^ ˙ i         i = 1 , 2 , , 6
Bring Equation (48) into Equation (49), and we can get the following:
V ˙ = p q e 2 p / q 1 q ^ 1 q 1 + v ¨ q ^ 2 q 2 + a ˙ t q ^ 3 q 3 s m + p q e 2 p / q 1 2 v ˙ 2 q ^ 4 q 4 + 2 v v ¨ q ^ 4 q 4 + v ˙ q ^ 5 q 5 s m + p q e 2 p / q 1 2 v v ˙ 2 q ^ 6 q 6 + a ^ a + v b ^ b + s m c ^ c p q e 2 p / q 1 s m k + G max s m + i = 1 6 1 λ i Δ q i Δ q ˙ i + 1 λ 7 Δ a Δ a ˙ + 1 λ 8 Δ b Δ b ˙ + 1 λ 9 Δ c Δ c ˙
To meet the stability rules of the Lyapunov formula, organize Equation (50), and the available parameter adaptation law is as follows:
q ^ ˙ 1 = λ 1 p q e 2 p / q 1 s m , q ^ ˙ 2 = λ 2 p q e 2 p / q 1 s m v ¨ q ^ ˙ 3 = λ 3 p q e 2 p / q 1 s m a ˙ 1 , q ^ ˙ 4 = λ 4 p q e 2 p / q 1 s m v ˙ 2 + v v ¨ q ^ ˙ 5 = λ 5 p q e 2 p / q 1 s m v ˙ , q ^ ˙ 6 = 2 λ 6 p q e 2 p / q 1 s m v v ˙ a ^ ˙ = λ 7 p q e 2 p / q 1 s m b ^ ˙ = λ 8 p q e 2 p / q 1 s m v c ^ ˙ = λ 9 p q e 2 p / q 1 s m v 2
As showed in Equations (50) and (51), discontinuous switching function k sgn ( s m ) , will result in unconnected control signals to keep the train running on the sliding surface. However, due to resistance deviation, measurement deviation, and other factors, the train cannot stay on the sliding surface, which will eventually cause chattering. To avoid chattering, replace signum function sgn ( s m ) , in Equation (47) with saturation function, s a t ( s m ) and ϕ is the saturation width.
Generally speaking, the distinguishing feature of typical sliding mode control is the use of a signum function. After the introduction of a saturation function, the controller will no longer be a classic sliding mode controller. According to References [41,42], the sat function implies a pseudo sliding mode in fact.
s a t ( s m ) = 1               s m > ϕ s m         s m ϕ 1         s m < ϕ

5.3. Disturbance Observer

The introduction of the saturation function will reduce the robustness of the system, but the line interference may be greater during train operation, especially on ramps. Disturbance observer (DOB) has been one of the most widely used robust control tools which can effectively improve robustness and anti-disturbance [43,44]. Therefore, a disturbance observer is introduced in this part.
The Equation (48) is introduced into the disturbance observation value and transformed into the following form:
u t = q p β e 2 2 p / q + q ^ 1 v + q ^ 2 v ¨ + q ^ 3 a ˙ t d ^ q ^ 1 d ^ + v ˙ r k sgn s m
As shown in Figure 13, the actual model and the nominal model are as follows:
G = e T d s τ s + 1 G n = e T n s τ n S + 1
According to the principle of disturbance observer, by considering the input term as the transfer function of u r and d , we can get the following:
G d d ^ s = d ^ d = Q s 1 + Q s G s G n s Q s G u r d ^ s = d ^ u r = G s Q s G n s Q s 1 + G s Q s G n s Q s
Furthermore, the expression of the disturbance observation can be obtained as follows:
d ^ = G u r d ^ s u r + G d d ^ s d
According to Equation (56), if G n s = G s , the required disturbance value can be accurately observed, but in the actual system, there may be a deviation between G n s and G s . From Equation (56), the deviation of the estimated value mainly comes from the error of the model:
Δ d ^ = Δ G s d u r 1 Δ G s
where Δ G s = G n s G s , Q s is the low pass filter, and Q s = 1 .
According to Equation (57), if there is an error in the model, the obtained observations will also have a certain deviation. The magnitude of the model error determines the magnitude of the deviation of the observations. After introducing the disturbance observer, the gain of the switching function will be weakened to a certain extent. In general, the error of the train model is small, so the disturbance observer can improve the anti-interference ability of the system and reduce the switching gain, to avoid chattering.

6. Experimental Result

To verify the validity, this section uses the interval between Zhuomashan station and Heluobu station in Qingdao Metro Line 6 as the simulation object.

6.1. Basic Test of RR-ABC

The vehicle data and line data (slope and speed limit) used are shown in Table 3, Table 4, and Table 5, respectively.
To verify the effectiveness of the AA-ABC algorithm, we compare it with the standard ABC algorithm and the GA, respectively. We set the given interval running time to 147 s and use different evolutionary algorithms for testing. For each algorithm, we conducted three repeated experiments. The results are shown in Figure 13.
Figure 14a is the comparison of the RR-ABC algorithm and the standard ABC algorithm; the line in the figure represents the changing trend of fitness. We can observe the following two results:
(1)
The initial value of the RR-ABC algorithm is better than that of the standard ABC algorithm;
(2)
Compared with the standard ABC algorithm, RR-ABC maintained good evolutionary performance, both in the early and late stages; especially in the late stage, the evolution efficiency of RR-ABC is significantly higher than that of ABC. This proves that the RR-ABC algorithm has good local search capabilities.
Figure 14b is the comparison of the RR-ABC algorithm and GA; we can observe that GA gets a local optimal solution twice. Relatively speaking, ABC has good global search capabilities
By combining Figure 14a,b, we can find that RR-ABC algorithm has better global search capabilities compared with ordinary evolutionary algorithms. Furthermore, it also improves the local search capabilities, which are poor in standard ABC algorithm. RR-ABC algorithm can quickly obtain the optimal solution.
Figure 15a is the optimal speed profile, and Figure 15b is the energy-consumption curve. Figure 15a shows the fixed speed limit curve, the economic speed limit curve, and the optimal speed curve. The curve in Figure 15b shows the cumulative upward trend of energy consumed by the train as the vehicle moves. We can observe the following results:
(1)
Running with the safety speed limit curve, the train can safely travel under the fixed speed limit curve;
(2)
The speed profile contains two coast sections;
(3)
When the train is in traction section and cruising section, the energy consumption increases. Furthermore, the rate during traction is higher than that in cruise. In contrast, the energy consumption remains a level-out during coast section and braking section.

6.2. Advanced Test about Different Interval Time and Different Mode

To verify the performance of the algorithm under different running times, set the given interval time to 141 and 153 s, respectively. The results are shown in Figure 16 and Figure 17.
In Figure 16, because of the reduction of running time, compared with the speed profile shown in Figure 15, we can observe that coast sections are shortened and the energy consumption is higher. In Figure 17, resulted from the growth of running time, we can observe that coast sections are extended and the energy consumption is lower.
The detailed data of the comparison are shown in Table 6. We can see that, with the increase of the given time, the energy consumption of the train gradually decreases, and the amount of decrease continues to fall. The reason is that the train can use more coast sections with increasing given time. However, the energy-saving effect of coast section is not linear. As the coast section becomes longer, the energy-saving effect of coast section will get worse. In addition, we can also see that the indexes about time and acceleration are within the allowable range, which means that the train maintains good punctuality and comfortableness while performing energy-saving driving.

6.3. Advanced Test about Different Driving Mode

In order to further test the performance of the algorithm, we compared the profile and consumption of optimal mode with that of cruising mode, which is widely used in engineering practice.
As shown in Figure 18, if we set the given interval time to 147 s, both for optimal mode and cruising mode, the profile of cruising mode maintains a level-out after a period of acceleration, which is definitely easy to plan. In comparison, following with the acceleration, the optimal profile conducts a coast condition. In addition, the energy consumption of the optimal profile is less than the cruising mode.
Detailed data are listed in Table 7; we can observe the following two results:
(1)
The energy consumption of optimal mode is 9.55% lower than that of cruising mode, so optimal mode is at energy saving;
(2)
The average acceleration rate of cruising mode is much lower than that of optimal mode, and although the index of optimal mode is qualified, the cruising mode is much better at comfortableness.

6.4. Basic Tracking Test about ATSMC+DOB Controller under Different Fixed Delays

To verify the effect of delay compensation, four groups of time delay experiments with different values are set up. The delay parameters of each group are shown in Table 8. Usually, T is less than 0.9 s and τ is less than 0.3 s in modern metro control system.
As shown in Figure 19a,c and Figure 20a,c, after a period of adjustment, the tracking curve can always converge to the reference curve. As shown in Figure 19b,d and Figure 20b,d, when the time delay becomes larger, the system applies more acceleration to keep it tracked. As a result, in all the given cases, the system has achieved accurate tracking. Since the time delay of the actual system generally does not exceed 0.9 s, the sliding mode controller is considered effective for the compensation of the known fixed time delay.

6.5. Tracking Test about ATSMC+DOB Controller

To verify the effect of the ATSMC+DOB, it is compared with the PID controller and the STSMC controller, respectively. Select the 153 s speed curve as the reference curve, and the comparison results are shown in Figure 21 and Figure 22.
Compared with the ATSMC algorithm, as shown in the Figure 21, the disturbance observer significantly reduces the control error, especially when the train is running on a steep slope. Furthermore, compared with the PID algorithm, as shown in Figure 22, the ATSMC+DOB algorithm significantly reduces the speed error, which will greatly improve the comfort, punctuality, and parking accuracy of the train.
In addition, as the focus of our attention, as shown in Figure 23, due to the precise tracking of the ATSMC algorithm, although its energy consumption is slightly higher than that of the reference curve, it still much better than the PID algorithm, which has bad energy saving due to a large tracking error. In conclusion, the ATSMC+DOB algorithm performs good tracking accuracy.
The detail data are listed in Table 9. The ATSMC+DOB algorithm performs better than the PID algorithm and the ATSMC algorithm in all indicators; PID algorithm is the worst, far worse than ATSMC+DOB algorithm and ATSMC algorithm. In terms of energy consumption, ATSMC+DOB is 7.58% lower than PID and 0.7% lower than ATSMC. Generally, the ATSMC+DOB algorithm perfectly meets the requirements of the metros for energy saving, punctuality, and comfortableness.
In addition, as we can see in Figure 21b, ATSMC+DOB has less speed errors than ATSMC, and ATSMC+DOB can converge faster. The detailed data in the table also prove this. The observer brings better anti-disturbance.

7. Conclusions

This paper carried out a series of researches on metro speed profile optimization and speed tracking control, and further designed corresponding algorithms and conducted simulation verification based on real data. The relevant conclusions are summarized as follows:
(1)
A multi-objective optimization model that considers energy consumption calculation and takes punctuality and comfortableness as penalty factors is established to optimize the train speed profile. This model takes the comfortableness of the metro and the punctuality into consideration when optimizing energy saving.
(2)
An optimization strategy that considers the metro EBI speed limit and the actual volume of the train (warning speed limit) is proposed. An improved ABC algorithm named RR-ABC algorithm is proposed for speed profile optimization. Compared with ordinary algorithms, the RR-ABC algorithm not only has good global search ability to avoid the local optimal solutions, but it also has excellent local search ability to improve the evolutionary efficiency.
(3)
A terminal sliding mode controller with disturbance observer (ATSMC+DOB) is designed by introducing parameter adaptation mechanism and disturbance observer. The controller has better robustness and anti-disturbance which brings minor speed tracking error and good energy saving.
(4)
The real data from Qingdao Metro Line 6 were used for the verification of the research. The simulation test results prove that the research on speed-profile optimization and speed-tracking control is effective.
The optimizations based on coast optimization and various speed limits proposed in the paper are targeted and only effective for a specific single interval, but for most metro lines, the conditions of each interval are different. Assuming that we design a different strategy for different intervals, it must not be a good way. Therefore, in future research, a more general optimization strategy is worthy of being studied.

Author Contributions

X.W. and Z.X. were responsible for modeling. P.S. and M.C. were responsible for data analysis and sorting. Q.W. revised the article. X.F. provided the data source. X.W. proposed the algorithm and wrote the article. All authors have read and agreed to the final version of the manuscript.

Funding

This work was supported by the National Key Research and Development Program of China, under Grant 2016YFB1200502.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lu, S.; Hillmansen, S.; Ho, T.K.; Roberts, C. Single-train trajectory optimization. IEEE Trans. Intell. Transp. Syst. 2013, 14, 743–750. [Google Scholar] [CrossRef]
  2. Li, Z.; Tang, T.; Gao, C. Long short-term memory neural network applied to train dynamic model and speed prediction. Algorithms 2019, 12, 173. [Google Scholar] [CrossRef] [Green Version]
  3. Su, S.; Li, X.; Tang, T.; Gao, Z.Y. A Subway Train Timetable Optimization Approach Based on Energy-Efficient Operation Strategy. IEEE Trans. Intell. Transp. Syst. 2013, 14, 883–893. [Google Scholar] [CrossRef]
  4. Wang, B.; Yang, J.; Jiao, H.N.; Zhu, K.; Chen, Y.Q. Design of auto disturbance rejection controller for train traction control system based on artificial bee colony algorithm. Measurement 2020, 160, 9. [Google Scholar] [CrossRef]
  5. Guo, X.; Wang, J.L.; Liao, F.; Teo, R. Distributed adaptive sliding mode control strategy for vehicle-following systems with nonlinear acceleration uncertainties. IEEE Trans. Veh. Technol. 2017, 66, 981–991. [Google Scholar] [CrossRef]
  6. Tan, C.K.; Wang, J.; Paw, Y.C.; Ng, T.Y. Tracking of a moving ground target by a quadrotor using a backstepping approach based on a full state cascaded dynamic. Appl. Soft Comput. 2016, 47, 47–62. [Google Scholar] [CrossRef]
  7. Bidikli, B. An observer-based adaptive control design for the maglev system. Trans. Inst. Meas. Control 2020, 42, 2771–2786. [Google Scholar] [CrossRef]
  8. Song, Y.D.; Song, Q.; Cai, W.C. Fault-tolerant adaptive control of high-speed trains under traction/braking failures: A virtual parameter-based approach. IEEE Trans. Intell. Transp. Syst. 2014, 15, 737–748. [Google Scholar] [CrossRef]
  9. Wang, Y.; Song, Y.; Gao, H.; Lewis, F.L. Distributed fault-tolerant control of virtually and physically interconnected systems with application to high-speed trains under traction/braking failures. IEEE Trans. Intell. Transp. Syst. 2016, 17, 535–545. [Google Scholar] [CrossRef]
  10. Lin, X.; Dong, H.; Yao, X.; Bai, W. Neural adaptive fault-tolerant control for high-speed trains with input saturation and unknown disturbance. Neurocomputing 2017, 260, 32–42. [Google Scholar] [CrossRef]
  11. Chen, Z.; Huang, F.H.; Chen, W.J.; Zhang, J.H.; Sun, W.C.; Chen, J.W.; Gu, J.; Zhu, S.Q. RBFNN-based adaptive sliding mode control design for delayed nonlinear multilateral telerobotic system with cooperative manipulation. IEEE Trans. Ind. Inform. 2020, 16, 1236–1247. [Google Scholar] [CrossRef]
  12. Qiao, L.; Zhang, W.D. Trajectory tracking control of AUVs via adaptive fast nonsingular integral terminal sliding mode control. IEEE Trans. Ind. Inform. 2020, 16, 1248–1258. [Google Scholar] [CrossRef]
  13. Liu, X.D.; Yu, H.S.; Yu, J.P.; Zhao, L. Combined speed and current terminal sliding mode control with nonlinear disturbance observer for PMSM Drive. IEEE Access 2018, 6, 29594–29601. [Google Scholar] [CrossRef]
  14. Cao, W.J.; Xu, J.X. Nonlinear integral-type sliding surface for both matched and unmatched uncertain systems. IEEE Trans. Automat. Contr. 2004, 49, 1355–1360. [Google Scholar] [CrossRef]
  15. Van, M.; Mavrovouniotis, M.; Ge, S.S. An adaptive backstepping nonsingular fast terminal sliding mode control for robust fault tolerant control of robot manipulators. IEEE Trans. Syst. Man Cybern. Syst. 2019, 49, 1448–1458. [Google Scholar] [CrossRef] [Green Version]
  16. Yang, X.; Chen, A.; Ning, B.; Tang, T. A stochastic model for the integrated optimization on metro timetable and speed profile with uncertain train mass. Transp. Res. Part B Methodol. 2016, 91, 424–445. [Google Scholar] [CrossRef]
  17. Howlett, P. The optimal control of a train. Ann. Oper. Res. 2000, 98, 65–87. [Google Scholar] [CrossRef]
  18. Howlett, P. An optimal strategy for the control of a train. Anziam J. 1990, 31, 454–471. [Google Scholar] [CrossRef] [Green Version]
  19. Koopialipoor, M.; Armaghani, D.J.; Haghighi, M.; Ghaleini, E.N. A neuro-genetic predictive model to approximate overbreak induced by drilling and blasting operation in tunnels. Bull. Eng. Geol. Environ. 2019, 78, 981–990. [Google Scholar] [CrossRef]
  20. Su, S.; Tang, T.; Li, X. Driving strategy optimization for trains in subway systems. Proc. Inst. Mech. Eng. Part F 2018, 232, 369–383. [Google Scholar] [CrossRef]
  21. Chuang, H.J.; Chen, C.S.; Lin, C.H.; Hsieh, C.H.; Ho, C.Y. Design of optimal coasting speed for MRT systems using ANN models. IEEE Trans. Ind. Appl. 2009, 45, 2090–2097. [Google Scholar] [CrossRef]
  22. Khmelnitsky, E. On an optimal control problem of train operation. IEEE Trans. Autom. Control 2000, 45, 1257–1266. [Google Scholar] [CrossRef]
  23. Sheu, J.W.; Lin, W.S. Energy-saving automatic train regulation using dual heuristic programming. IEEE Trans. Veh. Technol. 2012, 61, 1503–1514. [Google Scholar] [CrossRef]
  24. Song, Y.; Song, W. A novel dual speed-curve optimization-based approach for energy-saving operation of high-speed trains. IEEE Trans. Intell. Transp. Syst. 2016, 17, 1–12. [Google Scholar] [CrossRef]
  25. Cheng, R.J.; Song, Y.D.; Chen, D.W.; Ma, X.P. Intelligent positioning approach for high speed trains based on ant colony optimization and machine learning algorithms. IEEE Trans. Intell. Transp. Syst. 2019, 20, 3737–3746. [Google Scholar] [CrossRef]
  26. Karaboga, D.; Basturk, B. A powerful and efficient algorithm for numerical function optimization: Artificial Bee Colony (ABC) algorithm. J. Glob. Optim. 2007, 39, 459–471. [Google Scholar] [CrossRef]
  27. Youneng, H.; Chen, Y.; Shaofeng, G. Energy optimization for train operation based on an improved ant colony optimization methodology. Energies 2016, 9, 626. [Google Scholar]
  28. Rocha, A.; Araújo, A.; Carvalho, A.; Sepulveda, J. A new approach for real time train energy efficiency optimization. Energies 2018, 11, 2660. [Google Scholar] [CrossRef] [Green Version]
  29. Tan, Z.; Lu, S.; Bao, K.; Zhang, S.; Wu, C.; Yang, J.; Xue, F. Adaptive partial train speed trajectory optimization. Energies 2018, 11, 3302. [Google Scholar] [CrossRef] [Green Version]
  30. Farooqi, H.; Incremona, G.P.; Colaneri, P. Railway collaborative eco drive via dissension-based switching nonlinear model predictive control. Eur. J. Control 2019, 50, 153–160. [Google Scholar] [CrossRef] [Green Version]
  31. Yan, X.H.; Cai, B.G.; Ning, B.; Wei, S.G. Online distributed cooperative model predictive control of energy-saving trajectory planning for multiple high-speed train movements. Transp. Res. Part C 2016, 69, 60–78. [Google Scholar] [CrossRef]
  32. Chen, M.; Xiao, Z.; Sun, P.; Wang, Q.; Jin, B.; Feng, X. Energy efficient driving strategies for multi-train by optimization and update speed profiles considering transmission losses of regenerative energy. Energies 2019, 12, 3573. [Google Scholar] [CrossRef] [Green Version]
  33. Cao, Y.; Wang, Z.C.; Liu, F.; Li, P.; Xie, G. Bio-inspired speed curve optimization and sliding mode tracking control for subway trains. IEEE Trans. Veh. Technol. 2019, 68, 6331–6342. [Google Scholar] [CrossRef]
  34. Shang, F.; Zhan, J.; Chen, Y. An online energy-saving driving strategy for metro train operation based on the model predictive control of switched-mode dynamical systems. Energies 2020, 13, 4933. [Google Scholar] [CrossRef]
  35. Huang, Y.N.; Bai, S.; Meng, X.H.; Yu, H.Z.; Wang, M.Z. Research on the driving strategy of heavy-haul train based on improved genetic algorithm. Adv. Mech. Eng. 2018, 10, 16. [Google Scholar] [CrossRef] [Green Version]
  36. Liu, D.; Zhu, S.Q.; Xu, Y.X.; Liu, K. Train operation optimization with adaptive differential evolution algorithm based on decomposition. IEEJ Trans. Electr. Electron. Eng. 2019, 14, 1772–1779. [Google Scholar] [CrossRef]
  37. Li, X.; Lo, H.K. An energy-efficient scheduling and speed control approach for metro rail operations. Transp. Res. Part B-Methodol. 2014, 64, 73–89. [Google Scholar] [CrossRef]
  38. ISO. ISO 2631-5: 2018. Mechanical Vibration and Shock—Evaluation of Human Exposure to Whole-Body Vibration of International Standard; ISO: Geneva, Switzerland, 2018. [Google Scholar]
  39. Banharnsakun, A.; Achalakul, T.; Sirinaovakul, B. The best-so-far selection in artificial bee colony algorithm. Appl. Soft. Comput. 2011, 11, 2888–2901. [Google Scholar] [CrossRef]
  40. Kiran, M.S.; Hakli, H.; Gunduz, M.; Uguz, H. Artificial bee colony algorithm with variable search strategy for continuous optimization. Inf. Sci. 2015, 300, 140–157. [Google Scholar] [CrossRef]
  41. Utkin, V.; Guldner, J.; Shi, J. Sliding Mode Control in Electro-Mechanical Systems, 2nd ed.; CRC Press: New York, NY, USA, 2009; pp. 159–203. [Google Scholar]
  42. Ferrara, A.; Incremona, G.P.; Cucuzzella, M. Advanced and Optimization Based Sliding Mode Control: Theory and Applications; Society for Industrial and Applied Mathematics (SIAM): Philadelphia, PA, USA, 2019. [Google Scholar]
  43. Sariyildiz, E.; Oboe, R.; Ohnishi, K. Disturbance observer-based robust control and its applications: 35th anniversary overview. IEEE Trans. Ind. Electron. 2020, 67, 2042–2053. [Google Scholar] [CrossRef] [Green Version]
  44. Huang, J.; Zhang, M.S.; Ri, S.; Xiong, C.H.; Li, Z.J.; Kang, Y. High-order disturbance-observer-based sliding mode control for mobile wheeled inverted pendulum systems. IEEE Trans. Ind. Electron. 2020, 67, 2030–2041. [Google Scholar] [CrossRef]
Figure 1. Force analysis of single particle model.
Figure 1. Force analysis of single particle model.
Energies 13 06093 g001
Figure 2. External characteristic curve of traction force and braking force.
Figure 2. External characteristic curve of traction force and braking force.
Energies 13 06093 g002
Figure 3. Train multi-interval operation planning.
Figure 3. Train multi-interval operation planning.
Energies 13 06093 g003
Figure 4. Subinterval division.
Figure 4. Subinterval division.
Energies 13 06093 g004
Figure 5. Punctuality penalty function.
Figure 5. Punctuality penalty function.
Energies 13 06093 g005
Figure 6. Comfortableness penalty function.
Figure 6. Comfortableness penalty function.
Energies 13 06093 g006
Figure 7. Three stages in emergency braking intervention (EBI) protection.
Figure 7. Three stages in emergency braking intervention (EBI) protection.
Energies 13 06093 g007
Figure 8. Safety speed limit.
Figure 8. Safety speed limit.
Energies 13 06093 g008
Figure 9. Speed-limit curve.
Figure 9. Speed-limit curve.
Energies 13 06093 g009
Figure 10. Coast strategy.
Figure 10. Coast strategy.
Energies 13 06093 g010
Figure 11. Comparison of classic method and reinforcement method.
Figure 11. Comparison of classic method and reinforcement method.
Energies 13 06093 g011
Figure 12. Controller model.
Figure 12. Controller model.
Energies 13 06093 g012
Figure 13. Disturbance observer structure.
Figure 13. Disturbance observer structure.
Energies 13 06093 g013
Figure 14. Comparison of evolutionary performance: (a) Comparison about RR-ABC and ABC; (b) Comparison about RR-ABC and GA.
Figure 14. Comparison of evolutionary performance: (a) Comparison about RR-ABC and ABC; (b) Comparison about RR-ABC and GA.
Energies 13 06093 g014
Figure 15. Optimal speed profile: (a) The optimal speed profile made by optimization; (b) Cumulative energy consumption growth curve.
Figure 15. Optimal speed profile: (a) The optimal speed profile made by optimization; (b) Cumulative energy consumption growth curve.
Energies 13 06093 g015
Figure 16. Optimal speed profile. (a) The optimal speed profile made by optimization; (b) Cumulative energy consumption growth curve.
Figure 16. Optimal speed profile. (a) The optimal speed profile made by optimization; (b) Cumulative energy consumption growth curve.
Energies 13 06093 g016
Figure 17. Optimal speed profile. (a) The optimal speed profile made by optimization; (b) Cumulative energy consumption growth curve.
Figure 17. Optimal speed profile. (a) The optimal speed profile made by optimization; (b) Cumulative energy consumption growth curve.
Energies 13 06093 g017
Figure 18. Comparison of different driving modes: (a) Comparison about coast mode and cruise mode in speed profile; (b) Comparison about coast mode and cruise mode in cumulative energy consumption curve.
Figure 18. Comparison of different driving modes: (a) Comparison about coast mode and cruise mode in speed profile; (b) Comparison about coast mode and cruise mode in cumulative energy consumption curve.
Energies 13 06093 g018
Figure 19. Experiment 1 and Experiment 2.
Figure 19. Experiment 1 and Experiment 2.
Energies 13 06093 g022
Figure 20. Experiment 3 and Experiment 4.
Figure 20. Experiment 3 and Experiment 4.
Energies 13 06093 g023
Figure 21. Comparison of adaptive terminal sliding mode controller (ATSMC) and ATSMC+disturbance observer (DOB): (a) Comparison in speed curve; (b) Comparison in speed error.
Figure 21. Comparison of adaptive terminal sliding mode controller (ATSMC) and ATSMC+disturbance observer (DOB): (a) Comparison in speed curve; (b) Comparison in speed error.
Energies 13 06093 g019
Figure 22. Comparison of PID and ATSMC+DOB: (a) Comparison in tracking curve; (b) Comparison in speed error.
Figure 22. Comparison of PID and ATSMC+DOB: (a) Comparison in tracking curve; (b) Comparison in speed error.
Energies 13 06093 g020
Figure 23. Comparison of energy consumption of different controllers.
Figure 23. Comparison of energy consumption of different controllers.
Energies 13 06093 g021
Table 1. Symbol description.
Table 1. Symbol description.
SymbolDescription
F n traction force
E n energy consumption
a n acceleration
X n length of the subinterval
Table 2. International comfort standards.
Table 2. International comfort standards.
Comfort LevelAcceleration RateEvaluation
Level 1<0.315Very comfortable
Level 20.315–0.63Comfortable
Level 30.63–1.0Relatively comfortable
Level 41.0–1.6Uncomfortable
Level 51.6–2.5Very uncomfortable
Level 6>2.5Extremely uncomfortable
Table 3. Vehicle data.
Table 3. Vehicle data.
Vehicle ParametersEvaluation
Formation6 trains
Mass-AW2 (ton)339.6
Length (m)120
Basic resistance parameters ( r 1 )9.067
Basic resistance parameters ( r 2 )0
Basic resistance parameters ( r 3 )0.001334
Table 4. Ramp data.
Table 4. Ramp data.
Location (m)Slope (‰)
0–970
97–3772.357
377–108710
1087–192617.682
1926–20450
Table 5. Speed-limit data.
Table 5. Speed-limit data.
Location (m)Speed Limit (km/h)
0–17560
175–90080
900–183170
1831–203065
Table 6. Detailed data for comparison.
Table 6. Detailed data for comparison.
Given Time (s)Real Time (s)Consumption (kJ)Average Acceleration Rate (m/s3)
141143.6118,1470.0456
147148.1106,9230.0331
153155.3101,2400.0237
Table 7. Detailed data about the comparison.
Table 7. Detailed data about the comparison.
ModeGiven Time (s)Real Time (s)Consumption (kJ)Average Acceleration Rate (m/s3)
Optimal147144.6118,1470.0456
Cruising147146.3130,6260.0171
Table 8. Values about delay experiments.
Table 8. Values about delay experiments.
Experiment Number T (seconds) τ (seconds)
10.50.1
20.80.2
310.3
41.20.4
Table 9. Detailed data about the comparison.
Table 9. Detailed data about the comparison.
Performance IndicesATSMC+DOBATSMCPID
Arrival time error (s)1.21.46.7
Total consumption (kJ)102,528103,231110,207
Average accelerated change rate (m/s3)0.03710.03990.55
Average speed error (m/s)0.052760.19330.8731
Stop error (m)0.2680.9412.459
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, X.; Xiao, Z.; Chen, M.; Sun, P.; Wang, Q.; Feng, X. Energy-Efficient Speed Profile Optimization and Sliding Mode Speed Tracking for Metros. Energies 2020, 13, 6093. https://doi.org/10.3390/en13226093

AMA Style

Wang X, Xiao Z, Chen M, Sun P, Wang Q, Feng X. Energy-Efficient Speed Profile Optimization and Sliding Mode Speed Tracking for Metros. Energies. 2020; 13(22):6093. https://doi.org/10.3390/en13226093

Chicago/Turabian Style

Wang, Xiaowen, Zhuang Xiao, Mo Chen, Pengfei Sun, Qingyuan Wang, and Xiaoyun Feng. 2020. "Energy-Efficient Speed Profile Optimization and Sliding Mode Speed Tracking for Metros" Energies 13, no. 22: 6093. https://doi.org/10.3390/en13226093

APA Style

Wang, X., Xiao, Z., Chen, M., Sun, P., Wang, Q., & Feng, X. (2020). Energy-Efficient Speed Profile Optimization and Sliding Mode Speed Tracking for Metros. Energies, 13(22), 6093. https://doi.org/10.3390/en13226093

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop