Abstract

To reduce the traction energy consumption of urban rail trains by regenerating energy, a train traction and braking model was designed based on algebraic graph theory and train dynamics theory and the following consistency model of energy-saving operation of urban rail trains was constructed based on the conditions of the coordination coefficient and the operating condition conversion of trains. Not only to constantly update the consistency controller and save communication resources but also to optimize the energy-saving effect, the consistency algorithm of event triggering was used to search for the optimal operating conditions of trains and the energy-saving operation scheme for multiple trains was established. Taking the train diagram of a subway line in Jinan as an example, the energy-saving control scheme of four trains was solved by MATLAB simulation. The simulation results show that the model can not only ensure parking accuracy and punctuality but also energy savings effectively; that is, the proportion of the total regenerative energy used by the follower train in the actual energy consumption is increased from 3.32% to 10.76%, and the actual total energy consumption of the train is reduced by 9.23%.

1. Introduction

With the continuous expansion of the construction scale of urban rail transit, the energy consumption of urban rail trains shows a rapid growth trend. Among all kinds of energy consumption of urban rail transit systems, the proportion of traction energy consumption of trains in the total energy consumption is the most prominent. How to reduce traction energy consumption by optimizing the train operation is a common concern of operation companies and scholars [1]. Due to the short distance between urban rail transit stations, frequent braking during train operation will produce considerable regenerative energy, which can be fed back to the catenary for immediate use by other traction trains. Therefore, efficient and full utilization of renewable energy is an important way to reduce the energy consumption of urban rail transit systems [2].

At present, to improve the utilization of renewable energy, many scholars at home and abroad mainly study the optimization of train schedules and train coordinated control. The optimization of train schedules can improve the utilization rate of renewable energy and reduce actual energy consumption from the perspective of the entire system [3], mainly including optimizing the arrival and departure times, departure intervals, and interstation operation redundancy time allocation [46]. From the perspective of train cooperative control, some scholars have studied how to increase the heavy fall time of train traction and braking conditions through the cooperative operation process of multiple trains, so as to improve the utilization rate of renewable energy [710]. Although some methods mentioned in the literature [310] have improved the utilization rate of renewable energy, they have had some impact on the operation schedule of urban rail trains. The train cooperative control is extended to the general agent cooperative control, and its control basis is the consistency problem, which refers to that multiple control systems update their own information by the information about themselves and leaders, so as to make a certain negotiation variable tend to be consistent [11]. The control systems communicate by the network in practical application, and because of the bandwidth and other constraints of the network itself, it brings higher challenges to multisystem cooperative control [12]. Therefore, Tabuada and other scholars proposed a new method of event-triggered communication [13]. At present, the research on event-triggered multisystem consistency has made many achievements [1417]. Chen et al. [18] introduced an event trigger mechanism to avoid excessive consumption of communication resources. Li et al. [19] proposed an event trigger strategy to reduce the communication burden of the two channels from the device and model to the controller. Zhou and Chang [20] applied an elastic trigger mechanism to reduce the communication channel load and save network resources.

As an important consistency problem, following consistency has also attracted the attention of more and more scholars. In this kind of problem, the leader is not affected by other followers, and under the action of consistency agreement, a negotiation variable changes with the leader [2123]. The works [1123] have studied the application of event-triggering mechanism and consistency protocol in other fields, which provides a reference for this study.

Under the condition of not changing the train schedule [24], we studied the application of the following consistency event in the energy-saving operation of urban rail trains, introduced the event trigger mechanism into the following consistency collaborative control, took the weighted acceleration of the train as the collaborative variable, changed the secondary traction time of the follower train through the event trigger mechanism, and increased the renewable energy utilization rate of the leader train, and it not only reduced the traction energy consumption but also ensured the punctuality of train operation at each station.

2. Problem Description

Consider that the leading urban rail vehicle is the leader and the following vehicle is the follower, a multiagent system is formed. The dynamic models of the leader and follower are shown inwhere and are the states of the leader and follower, respectively, are the system matrices, is the control input of the follower, is the traction control coefficient, and represents a symmetric input saturation function, as shown in

According to the actual situation of urban rail transit, under the control input saturation constraint of the follower, a consistency control strategy is designed to reduce the update frequency of the controller, so as to reduce network consumption and avoid network congestion and delay the network’s life. Therefore, it is assumed that the leader and follower trains are installed with an event monitor to determine the time when communication is required and the control input is updated.

The research of this paper is based on the following assumptions:(1)The running time and dwell time of each station interval of the subway train are constant and cannot be adjusted.(2)The operation curve of the leader train is determined by the optimal energy-saving operation scheme of the traditional single train operation, which is not variable, that is, the maximum traction force is adopted for traction under traction conditions, the maximum braking force is adopted for braking under braking condition, and the coasting condition is adopted in the middle [25]; The operating condition of the follower train is to draw to the set speed with the maximum traction force and then coasting. During this coasting, consistent coordination is accepted. According to the coordination or coasting to the last moment that can meet the operation time and distance, retraction (traction equivalent to the braking force)-coasting braking (maximum braking force) occurs.(3)Assuming that the train is in the same power supply arm, the regenerative braking electric energy of the leader train can be used by other following trains immediately and the remaining energy can be consumed by trackside resistance. The utilization process is shown in Figure 1. In Figure 1, train 1 is in the traction condition stage and train 0 is in the braking condition stage. The regenerative energy generated by train 0 braking can be immediately utilized by train 1, which is reflected in the feasibility of path 1. In addition, when the electric energy required by train 1 during traction is less than the regenerative energy generated by train 0, the excess regenerative energy is consumed by trackside resistance, that is, path 2 is feasible; in addition, when the electric energy required by train 1 is greater than the regenerative energy generated by train 0, the regenerative energy generated by train 0 is feasible. Path 4 will not exist. At this time, the traction substation will transmit enough electric energy for train 1 to ensure the normal operation of train 1, which is reflected in the feasibility of path 3.(4)To ensure the accuracy of train operating time, it is assumed that the train will not be affected by coordinated control under the outbound traction condition and inbound braking condition, that is, the maximum traction force and maximum braking force will remain unchanged.(5)Due to the influence of the urban rail transit operating environment and operation requirements, the maximum traction acceleration of the train is set as , the maximum braking acceleration is set as , and the maximum speed limit is set as ,(6)The consistency of multitrain energy-saving cooperative operations is studied. The train is regarded as a particle, and the operation conditions such as ramps and tunnels of the line are not considered for the time being.

3. Train Operation Model

3.1. Force Analysis

The train is mainly affected by traction, braking force, and resistance during operation, as shown inwhere is the resultant force of the train, is the train traction force, is the train braking force, and is the train resistance; when the train is in traction condition, the value of is zero and negative, the values of and are zero, is negative in the coasting condition, the value of is zero, and at the same time, and are negative in braking conditions.

Train resistance refers to the force opposite to the running direction of the train due to the influence of friction, impact, vibration, the plane and section of the line, and other external conditions during train operation. According to the causes, train resistance is divided into basic resistance, additional resistance, and starting resistance [26]. The following train operation processes under traction, coasting, and braking conditions are shown inwhere , , and are the acceleration, speed, and driving distance of the train in the traction condition, respectively, and and are the initial speed and time when the train starts. , and are the acceleration, speed, and driving distance of the train during coasting, and are the initial speed and time of the train during coasting, , and are the acceleration, speed, and driving distance of the train under braking condition, and and are the initial speed and time of the train during braking.

3.2. Model Solving

According to the assumptions in the previous section of this paper, the leader train adopts the maximum traction force for traction, the maximum braking force for braking, and the coasting condition in the middle. The follower train adopts the maximum traction force for traction under the traction condition, enters the first coasting condition when reaching the set speed, and according to the coordination condition or coasting to the last moment that can meet the running time and running distance, retraction (traction equivalent to the braking force)-coasting braking (maximum braking force, parameter setting and solution steps are shown in Figure 2.

In Figure 2, and tt2, respectively, represent the initial speed and time of secondary traction of urban rail trains; zstan and n represent the total number of stations completed of urban rail trains and the total number of target operation stations, respectively. According to the operation steps of the leader train and the follower train, the premise of solving the train objective function is the operation sequence, that is, the solving equations. Once the operation sequence is changed, the train operation objective function needs to be recalculated and adjusted.

4. Consistency Agreement and Analysis

4.1. Consistency Agreement

Definition 1. If the state of the agent finally meets under any initial conditions, it is said that the system obtains master-slave consistency. The purpose of energy savings is achieved by the immediate utilization of braking regenerative energy. The total energy consumption of train i in a certain station is Ei, which is shown as follows:where is the total traction energy consumption of the train, is the braking regenerative energy of the leader train, is the braking regenerative energy of the leader train absorbed by the train, is the traction instantaneous power of trains, is the braking instantaneous power of the leader train, is the traction force of the train at the moment, and is the braking force of the leader train at the moment; they can be expressed as follows:where is the mass of trains, is the acceleration of trains, is the mass of the leader train, is the acceleration of the leader train (negative value under braking condition), and and are the basic running resistances of the follower train and the leader train, respectively; they can be expressed as follows: is the unit of the basic running resistance of the train, is the gravitational acceleration, and , , and are the empirical parameters determined by the experiment, respectively. Multitrain energy-saving cooperative following consistency can be described as follows:where is the acceleration of the leader train, is the acceleration of train i, and and are both energy-saving synergy coefficients; generally, is any real number, . is defined as the difference between the weighted acceleration value of the following train i and the weighted acceleration value of the leader train. is called the collaborative deviation degree, which is recorded as . During the actual operation of the train, we do not need the acceleration of each train to be completely consistent with the acceleration of the leader. When the acceleration state meets our requirements, we can also consider that a collaborative agreement is reached, that is, ; the multitrain system reaches coordination, where is the synergetic error.
Based on the consistent design objective, the system regulation coefficient is introduced and the event driving function is defined, as shown in the following equation:If equation (17) does not hold, the system will update the controller; otherwise, the system will not be driven. is the error function of the system at time t. The system regulation parameter shall meet .
A control protocol with input saturation can be designed by the driving function, as shown in the following equation:Since the negotiation variable is train acceleration, (18) can be converted into the following equation:where , and are the event-driven event sequences of the train and the leader train, respectively, and are positive integers, is the control gain matrix, and represents constants greater than 0.
Since the leader train is not affected by the consistency protocol, only the follower train is considered in the research process and the control structure of the follower train i is shown in Figure 3.

4.2. Consistency Analysis

Lemma 1 (see [27]). It is supposed that for a set , where and are vector elements in , there is a dead zone nonlinear function , so that the inequality holds, where is any positive definite diagonal matrix. The condition provides a feasible basis for transforming the system stability problem into an LMI solution problem.

Theorem 1. Considering equations (1) and (2), it is supposed that the topological graph of the system is a strongly connected graph and contains a directed spanning tree, if there is a positive definite matrix , any diagonal positive definite matrix , matrix , matrix , a positive number , and two positive scalars and , ; these satisfy the following equations:where , , and . Then, the error trajectory of the follower and the leader will be included in this ellipsoid; that is, the follower can achieve the desired consistency.
Prove: in vector , represents the state error root of the ith follower and leader; then, the error system can be shown in the following equation:where and the dead zone nonlinear function will be introduced to replace the saturated nonlinear function; if , then (22) can be converted into the following equation:From event-driven function (17), the following can be obtained:Therefore, it can be obtained that (26) is valid.Lyapunov function is selected for stability analysis;we assume that (27) holds.where P is a positive definite symmetric matrix because it is stable, and the equation has a unique positive definite solution. To prove the stability of (22), it is necessary to prove that function (27) is satisfied; that is, it is necessary to prove that for any and e are established under equation (17), satisfying the driving conditions.
The S-procedure method can transform the problems that do not belong to convex constraints into linear matrix inequality constraints, so it is convenient to solve them by using linear matrix inequality [28].Order: is a real valued objective function defined on the linear vector space V, with the following two conditions. S1: for any , and the constraint condition is ; S2: if , then and . SupposeS1: S2: According to the ss-procedure method, the premise of proving S1 is to prove S2. Therefore, S2 can be simplified to an equation.where ; suppose , ; for any can satisfy the conditions in Lemma 1, that is ; the establishment of LMI can ensure the ellipsoid [27]. To prove that is true, only (29) needs to be proved to be true.Namely,Suppose ; we just prove that . Among them,where ; therefore, according to Lyapunov’s asymptotic stability theorem, the error system is asymptotically stable; that is, the follower and the leader achieve the desired consistency.
To show that there is no Zeno phenomenon in the whole event-triggering control process, reference [16] gives the following definitions.

Definition 2. If the event driving interval meets the conditions , the system can eliminate the Zeno phenomenon;that is, the system will not be infinitely triggered within the limited driving time.

Theorem 2. Consider master-slave systems (1) and (2), in which communication topology is strongly connected. Based on the event-driven control protocol (19) and driving condition (17), each train can obtain the initial state of other trains at the initial time and each train can transmit the driving time state information to other trains; then, all trains applying control protocol (19) and driving conditions (17) will not generate the Zeno phenomenon.
Prove: because , we can get ; from , we will getSince there is (33) and , (32) can be changed intoFormula (35) is valid.When the system is driven, (36) can be established.Because exists, the system will have when driven, that is,

To sum up, if the event driving interval is greater than zero, the system can eliminate the Zeno phenomenon and Theorem 2 is proved.

5. Simulation and Analysis

A MATLAB simulation is used to verify the correctness of the theoretical scheme. The number of trains in the system is set to 4, including 1 leader and 3 followers. The train operation diagram refers to a part of a certain line of a subway operation company. The distance between stations of train operations is shown in Table 1, and the train operation parameters are shown in Table 2.

It can be seen from Tables 1 and 2 that, under the premise of not changing the system operation time, the research does not need to consider train safety. Set train 0 as the leader train (the leading train) and train 1, train 2, and train 3 as the following train (the following train). The leader train 0 operates in the optimal operation sequence of a single train, and the follower trains 1, 2, and 3 operate in the cooperative energy-saving optimal operation sequence. The set parking accuracy is not greater than 0.3 m, and the punctuality is less than 1 s. The acceleration is regarded as a collaborative variable, and a consistency protocol is introduced, along with its energy-saving synergy coefficient and . The value of ui is shown in Table 3.

It can be seen from Table 3 that only under first inert working condition can the follower train carry out the following consistency coordinated control.

Then, the operation curve of 4 trains can be obtained, as shown in Figure 4.

As can be seen from Figure 4, the following train 1 runs between station 8 and station 9, train 2 runs between station 1 and station 2, between station 2 and station 3, between station 6 and station 7, and between station 8 and station 9, and train 3 runs between station 1 and station 2, between station 2 and station 3, and between station 6 and station 7. The follower train is affected by the acceleration change of the leader train between the corresponding station areas, and its weighted acceleration follows consistently. With the change of its working condition transition point, the secondary traction is carried out at the event trigger time (the leader train braking time), and the follower train starts to absorb and utilize the regenerative braking energy of the leader train. Energy-saving optimization is carried out on the premise of ensuring operation time and train parking accuracy. The variation trend of consistency coordination deviation is shown in Figure 5.

As shown in Figure 5, the follower train realizes coordination and the trigger time of the follower train is shown in Figure 6, which greatly reduces the system communication frequency and avoids the waste of network resources.

According to equations (8)–(11), the energy consumed and regenerative energy of the follower train before and after consistency control to complete the above simulation events are given, as shown in Table 4.

Figure 7 shows the energy consumption of the follower train.

It can be seen from Figure 7(a) that the newly added follower train 1 absorbs and utilizes the regenerative energy of the leader train when it operates between stations 8 and 9. It can be seen from Figure 7(b) that the newly added follower train 2 absorbs and utilizes the regenerative energy of the leader train when it operates between stations 1 to 2 and 2 to 3. The newly added follower train absorbs and utilizes the regenerative energy of the leader train when it operates between stations 6 to 7 and 8 to 9, as shown in Figure 7(c). It can be seen that when follower train 3 operates between stations 6 and 7, its absorbed regenerative energy has been improved in detail, and the newly added regenerative energy is absorbed and utilized by the leader train when it operates between stations 1-2 and 2-3. This is consistent with the operation interval of the follower train and the leader train in Figure 4. It can be seen from Table 4 and Figure 7(d) that, on the premise of ensuring the parking accuracy and punctuality of the original train diagram, after the consistency optimization, the regenerative energy absorbed and utilized by each follower train has been improved to varying degrees and the proportion of the total utilized regenerative energy in the actual energy consumption has increased from 3.32% before optimization to 10.76% after optimization.

6. Conclusions

Taking the urban rail train as the research object, this paper establishes the energy-saving operation model of the urban rail train and studies the following consistency of urban rail train cooperative control. It is assuming that each urban rail vehicle is an agent, the synergy coefficient and synergy deviation degree are introduced, and the train acceleration is taken as the negotiation variable to study the feasibility of the follower train to achieve synergy. To save network resources and the amount of calculations of train coordination consistency, when designing the follower train consistency protocol, make each follower vehicle update the cooperative controller only at the time of event triggering. Through the consistency analysis, it is proved that the multitrain system can achieve the leader following consistency under the consistency protocol, and the Zeno phenomenon is eliminated. Finally, it is verified by referring to the train diagram of a certain line of a subway operation company. Through simulation analysis, based on the energy-saving method proposed in this paper, on the premise of ensuring the original parking accuracy and punctuality, the proportion of the total renewable energy used by the follower train in the actual energy consumption is increased from 3.32% to 10.76%, reducing 9.23% of the actual total energy consumption of the train. The proposed method improves the utilization rate of regenerative energy of the leader train. Considering the leader following consistency under switching topology is the next key research content, which can expand the utilization range of renewable energy and further save traction energy consumption.

Data Availability

The data used to support the findings of this study were supplied under license and so cannot be made freely available. Requests for access to these data should be made to the corresponding author.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (NSFC) under Grants 72061021 and 62063013, Science and Technology Plan Project of Gansu (20JR10RA251 and 21JR7RA284), Youth Fund Project of Lanzhou Jiaotong University (2021018).