A practical model for the train-set utilization: The case of Beijing-Tianjin passenger dedicated line in China

As a sustainable transportation mode, high-speed railway (HSR) has become an efficient way to meet the huge travel demand. However, due to the high acquisition and maintenance cost, it is impossible to build enough infrastructure and purchase enough train-sets. Great efforts are required to improve the transport capability of HSR. The utilization efficiency of train-sets (carrying tools of HSR) is one of the most important factors of the transport capacity of HSR. In order to enhance the utilization efficiency of the train-sets, this paper proposed a train-set circulation optimization model to minimize the total connection time. An innovative two-stage approach which contains segments generation and segments combination was designed to solve this model. In order to verify the feasibility of the proposed approach, an experiment was carried out in the Beijing-Tianjin passenger dedicated line, to fulfill a 174 trips train diagram. The model results showed that compared with the traditional Ant Colony Algorithm (ACA), the utilization efficiency of train-sets can be increased from 43.4% (ACA) to 46.9% (Two-Stage), and 1 train-set can be saved up to fulfill the same transportation tasks. The approach proposed in the study is faster and more stable than the traditional ones, by using which, the HSR staff can draw up the train-sets circulation plan more quickly and the utilization efficiency of the HSR system is also improved.


Introduction
During the past few years, China has been in a period when high-speed railway (HSR) developed rapidly, with the operation mileage increasing from 404 kilometers in 2008 to over 19 thousand kilometers at the end of 2015. With the gradually formed HSR network, the mobility of society, including people and cargo, is increasingly dependent on the HSR. Even though the imbalance between transport supply and demand has been slightly mitigated by introducing HSR, great efforts are still needed to improve the transport capability of HSR, so as to satisfy the explosively increasing mobility demand. Generally, the transport capability of HSR can be improved from two aspects. Firstly, speeding up the infrastructure construction can raise the upper bound of capacity. Secondly, improving the operation and management of HSR is an effective way to make the most of the existing infrastructure. Considering the long construction period, it is impractical to increase the HSR operation mileage within a short time. Therefore, enhancing the utilization efficiency of HSR is an achievable and key solution to improve the transport capacity and mitigate the social contradiction of imbalance between transport supply and demand.
As the carrying tools of HSR, train-sets play a vital role in the process of completing society mobility task, the utilization efficiency of which is one of the most important factors influencing the transport capacity of HSR. Therefore, how to improve the train-set utilization efficiency has been a research hotspot in railway transport for several decades. The utilization efficiency of train-sets is directly dependent on the train-set utilization plan, which is made in a given train graph by a number of train-sets to complete the trip tasks. When a train-set utilization plan is drawn up, the amount of train-sets in demand as well as the working time and idle time of each train-set can be obtained. Obviously, by following a good train-set utilization plan, the proportion of train-sets' working time in each day can be dramatically increased. Thus, a number of train-sets can be saved to serve for other train schedules, which means more trips can be added to the train graph without considering the limitation of available train-sets.
During the past few decades, extensive research has been carried out worldwide to make an attempt to draw up a high-quality train-set utilization plan, leading to the development of various models and techniques. In Europe, the train-set utilization optimization problem mainly focused on the coupling and uncoupling arrangements of train-sets to carry out the trip tasks in the train graph. Schrijver [1] firstly paid attention to the train-set scheduling problem. He proposed a model based on the minimum cost flow theory and solved it by the software CPLEX. Abbink et al. [2] explored the train marshalling problem during the peak period, and the optimization goal was set to minimize the shortage number of train seats. Alfieri et al. [3] addressed the train-set utilization problem on a single train line in a single day. An integer programming model was described to obtain the circulation of rolling stock considering the order of the train units in the compositions. Fioole et al. [4] put forward a mixed integer programming model on basis of previous research achievements, and the improved branch and bound algorithm was applied to obtain the optimization solution. Barnhart et al. [5], Peeters and Kroon [6], and Holmberg & Yuan [7] used the Dantzig-Wolfe decomposition approach [8] to solve the train-set circulation problem, instead of the classical flow decomposition method used in other network flow problems. Maróti and Kroon [9] proposed an integer programming model with the consideration of maintenance constraint. Moreover, he designed a heuristic algorithm to solve the model and obtained a feasible solution. However, the operational mode of train-sets in China is quite different from that in European countries. In China, the train-sets are used as a whole, so there is no need to consider the problem of coupling and uncoupling. Zhao at el. [10] discussed the different utilization modes of train-sets comprehensively and they found that the mode of unrestrained section can improve the utilization efficiency. Nie et al. [11] proposed an integer programming model considering the deadhead trains condition, and applied the Hungary algorithm to solve the problem. Wang et al. [12] set up a train-set utilization model to complete the trip tasks in a cyclic train timetable, and a column generation algorithm was used to obtain the optimal solution. Then Wang et al. [13] improved this model, and designed a method of branch and price, and the findings have inspired latter research. Zhao [14] believed that the train-set utilization problem can be regarded as a typical NP hard problem, and it cannot be solved directly by ready-made software. Characterized by easy-to-understand and fast, a large number of heuristic algorithms have been applied to solve the train-set utilization problem. Yang et al. [15] built a train-set connection network model aiming to maximize the train-set utilization efficiency, and designed a genetic algorithm to solve the model. Chen et al. [16] proposed a multi-objective optimization model to maximize the equilibrium degree of train-set utilization. Considering the number of train-sets in use, a simulated annealing algorithm was designed to solve the model. Hong et al. [17] put forward an integer programming model considering the maintenance and train graph constraints, and a heuristic approach was applied to solve this problem.
Recent research on the train-set utilization problem mainly focuses on developing new methodologies and introducing rescheduling procedures. The train-set utilization problem has been proved to be a TSP (Traveling Salesman Problem) with some special constraints. Based on TSP, Zhang [18] built a train-set utilization model in a railway network by taking both the running time constraint and running distances constraint into consideration, and an ant colony algorithm was used to solve the model; Zhou [19] proposed an integer programming model to obtain the train-set utilization plan, considering the constraint of Level 1 maintenance standard, and an improved ant colony algorithm was designed to solve the model. Thorlacius [20] proposed an integrated rolling stock planning model that took into account all practical requirements simultaneously for rolling stock planning at DSB S-tog. By using a hill climbing heuristic algorithms to solve this model, the efficiency of new rolling stock plan was increased by 2%. Jørgen [21] designed two approaches, one of which is Mixed Integer Linear Program (MILP) solved by CPLEX and the other is column and row generation approach, and both approaches are sufficiently fast in real-time setting. Lai [22] put forward a hybrid heuristic process to solve this problem for Taiwan Railway Administration (TRA). The efficiency of rolling stock use was increased by about 5% and the solution time was reduced significantly from 3h to 11.2s. Tatsushi [23] proposed a column generation and Lagrangian relaxation heuristics, and good lower and upper bounds for 300 trips were found within reasonable computing time. Due to disturbances or disruptions from the real world, the train-set circulation plan need to be rescheduled in order to maintain a feasible rolling stock circulation, and generally, this procedure should be carried out in a short time. Some researchers studied the train-set utilization problem in the rescheduling procedures. Cadorso [24] proposed a new approach to obtain better and more robust circulations of the rolling stock train units, solving the rolling stock assignment while accounting for the train routing problem. By designing a heuristic based on the benders decomposition, more robust and efficient solutions were obtained. Wagenaar [25] proposed three different models for rescheduling, namely the Extra Unit Type model, the Shadow-Account model and the Job-Composition model. All models were tested on instances of Netherlands Railways. Results showed that the models were able to take maintenance appointments into account efficiently. Cacchiani [26] presented an overview of recovery models and algorithms for real-time railway disturbance and disruption management, and the application of these models and algorithms in real-life railway systems would be instrumental in improving the quality of provided railway services, leading to an increased utilization of the involved railway systems. Table 1 provides a systematic comparison of key model components and solution methods in the existing research.
Although there is a comprehensive body of literature on the train-set utilization problem, no results of previous research can be directly used into the practical production. There are two major reasons. First, limited studies take the maintenance constraints into consideration, which is actually a key factor influencing the compilation of the train-set utilization plan. Second, even if the maintenance is considered, the solutions of models are obtained by applying heuristic algorithms, which is unstable and inaccurate. Therefore, to realize the objective that the train-set utilization plan can be generated by computers automatically and used in practice, an integer programming model considering the maintenance constraints is built to maximize the train-set utilization efficiency. Specially, a two-stage method is proposed in the model solving process in this study to obtain a unique solution. With this method, the train-set utilization plan can be obtained by a computer-generated system.
The remainder of this paper is organized as follows. Section 2 describes the basic situation of the train-set utilization problem. In Section 3, the optimization model of the train-set utilization plan is developed. Section 4 introduces the solution method of the model. In Section 5, a case study based on the Beijing-Tianjin HSR line is performed to illustrate the model application. Finally, Section 6 provides the conclusions marks.

Problem description
The train-set utilization plan is drawn up to identify the work arrangements of train-sets according to the given train graph as well as the rules and regulations of maintenance. For train-sets, except for some necessary operations in the railway stations, the main work includes two parts: one is to complete all the trip tasks in the given train graph; and the other is to carry out maintenances. Therefore, the train-set utilization plan should figure out the following aspects: the number of needed train-sets; the specific train-set undertaking its corresponding trip; the sequence of the trips undertaken by the same train-set; and the time and location each train-set is maintained. Different train-set utilization plans can be made to complete the trip tasks in a given train graph. Therefore, in this study, the optimization problem of the train-set utilization plan is to complete all the trip tasks in a given train graph by utilizing the least train-sets on the premise of satisfying the requirements of the maintenance. Fig 1 illustrates the framework of the train-set utilization optimization problem. Firstly, some basic information should be prepared as the input parameters, including: the train graph, the rules and regulations of maintenance and station working. Secondly, in order to model the train-set utilization plan, a train-set utilization network is represented by a directed graph. Thirdly, an integer programming model is established to minimize the total connecting  time, considering the spatial constraints, the time constraints, the maintenance constraints, and the unicity constrains. Finally, a two-stage approach is designed to solve the problem, and an optimal train-set utilization plan is obtained as the output. Fig 2 illustrates the relationship between the train graph and the train-set utilization plan. Fig 2(a) is a simple sample of a given train graph, including three stations and eight trips. After adding the train-sets' connection relationships and maintenance information, Fig 2(b) represents the train-set utilization plan. In this sample, to complete the eight transportation tasks, it needs at least three train-sets. Particularly, Trip G82 and G89 are undertaken by one train-set; Trip G84, G87, and G88 are undertaken by the second train-set; and Trip G83, G86, and G91 are undertaken by the last train-set. The train-set undertaking Trip G83 should be maintained at the inspection and repair depot near Station B after completing the task of Trip G86.

Train-set utilization plan modeling Train-set utilization network representation
In order to model the train-set utilization plan, a train-set utilization network is represented by a directed graph G(V, E). The node set V is a union of subsets V T and V M , which represent trips and maintenances, respectively. Each node i 2 V T is defined as a tuple Each node k 2 V M is defined as a tuple ðφ k ; t a k Þ, where φ k indicates the inspection and repair depot, t a k indicates the time train-set departing from the inspection and repair depot to arrival station after maintenance. As shown in Fig 3, there are three kinds of arcs in the train-set utilization network, including connection arcs (i.e., trip!trip), maintenance arcs (i.e., trip!maintenance) and arcs after maintenance (i.e., maintenance!trip). Therefore, the arc set E consists of three disjoint sets E C , E Q and E H , where E C represents directional arcs between two trip nodes, i.e., connection arcs, E Q represents directional arcs from trip node to maintenance node, i.e., maintenance arcs and E H represents directional arcs from maintenance node to trip node, i.e., arcs after maintenance. The weight of connection arcs o C ij means the connection time between Trip i and Trip j. Generally, when the arrival station of Trip i is different from the departure station of Trip j, the connection time should be defined as infinite. It means that that the two trips cannot be connected due to the huge waste of train-set empty running. As shown in Fig 4, when two trips can be connected, the values of o C ij is determined by the sequential order of the two trips' departure time and arrival time. If the departure time of Trip j is after the arrival time of Trip i, the o C ij is equal to t d j À t a i , otherwise, the o C ij is equal to 1440 þ t d j À t a i , i.e., the train-set will undertake Trip j the next day. Based on above analyses, the values of o C ij can be calculated as Eq (1).
The weight of maintenance arcs o Q ik means the duration from the moment train-set arriving at its arrival station, t a i , to the moment train-set returning back to its arrival station after maintenance, t a k . As shown in Fig 5, when t 0, it means that the train-set can only undertake Trip j the next day, and the value of o Q ik is equal to 1440 þ t a k À t a i . Therefore, the values of o Q ik can be calculated as Eq (2).
( The weight of arcs after maintenance o H ki means the duration from the moment train-set returning back to its arrival station after maintenance, t a k , to the moment train-set departing from the station, t d i . As shown in Fig 6, when t d i À t a k ! 0, the value of o H ki is equal to t d i À t a k . When t d i À t a k 0, it means that the train-set can only undertake Trip i the next day, and the value of o H ki is equal to 1440 þ t d i À t a k . Therefore, the values of o H ki can be calculated as Eq (3).

Analyses of constraints
To undertake a trip task, the train-set should satisfy four specific requirements, including the spatial constraint, the time constraint, the maintenance constraint and the unicity constraint. The detail analyses are as follows: 1. Spatial constraints. As mentioned above, an empty running train-set needs to be dispatched to connect two trips. When the departure station of the latter trip is different from the arrival station of its former trip, it will cause huge waste. Therefore, as shown in Eq (4), when a train-set undertake two adjacent trips, the departure station of the latter trip, s d j , is restricted to be the same with the arrival station of the former trip, s a i , where i j.
2. The time constraints. The time constraints means that the duration of connection arcs, maintenance arcs and arcs after maintenance should respectively satisfy the requirement of the minimal necessary working procedure duration between two tasks, as shown in Fig 7. When a train-set arrives at the station after completing the trip task, some preparation work need to be carried out before departing for the next trip, such as cleaning, pollution discharge, etc. The duration needed for such preparation work is defined as operation time, T o . Therefore, to ensure that two trips can be undertaken by one train-set, the connection time o C ij of the two trips should be longer than the operation time T o , as shown in Eq (5).
When a train-set completes a trip task and reaches to the maintenance standard, it should be maintained in the inspection and repair depot. As shown in Fig 8, four essential procedures will be followed by the train-set to complete the maintenance task. Table 2 lists the detailed information about the four procedures along with their duration definitions.
To ensure a trip node can be connected to a maintenance node, there should be enough time for the train-set to go through these procedures. Therefore, Eq (6) restrains that the weight of maintenance arc o Q ik should be larger than the duration of the process Once completing the maintenance task, train-sets can be utilized to undertake the trip tasks. However, as mentioned above, a series of preparation work should be carried out before departing to ensure the train-set working normally, as shown in Fig 9. Therefore, in order to connect a maintenance node with a trip node, Eq (7) restrains that the weight value of arc after maintenance o H ki should be larger than the duration of preparation work o H ki . o H ki ! T c ð7Þ 3. The maintenance constraints. For safe operation concern, the train-sets must be strictly complied with the maintenance rules and standards. In China, the train-sets' maintenance standards are divided into five levels according to the accumulated operation time and mileages. In this study, it is assumed that the train-set utilization plan is drawn up by treating one day as a cycle. Generally, the operation time standard of Standard Level 2 to Standard Level 5 is more than a month. Therefore, only the Standard Level 1 needs to be considered in the study. The Standard Level 1 requires that when the accumulated operation time reaches up to 48 hours or the accumulated operation mileages reaches up to 4,000 km, the train-set must be maintained. Additionally, in practical applications, a 10% fluctuation compared with standard is allowable (e.g., Train-sets should be maintained when the accumulated operation mileages are within a range from 3600 km to 4400 km).
When undertaking several trips, the accumulated operation time and mileages of the trainset are  The train-set will return back from the inspection and repair depot to the station.

Running time T b
Total Duration of the process T z = T a + 2T b + T m https://doi.org/10.1371/journal.pone.0175698.t002 4. The uniqueness constraints. The uniqueness involves two aspects. One is that after completing a trip task, each train-set should either undertake another trip or go for maintenance. The other one is that for each trip in a given train graph, only one train-set can be allocated to complete the task.
The binary variable x C ij is defined to describe whether the trip node i and trip node j are connected by the connection arc. The binary variable x Q ik is defined to describe whether the trip node i and maintenance node k is connected by the maintenance arc. The binary variable x H ki is defined to describe whether the maintenance node k and trip node i are connected by the arc after maintenance. If the two nodes are connected by arcs, the value of the binary variable is 1, otherwise 0, as shown in Eqs (10)-(12). Eq (13) restrains when a train-set completing Trip i, it should either undertake Trip j or go for Maintenance k. Eq (14) restrains that Trip j should be undertaken by a train-set either just completing Trip i or Maintenance k.
Optimization goal Railway passenger transport corporations aim to realize profit maximization on the premise of ensuring passengers' life and property security. Therefore, the cost is a key factor influencing the operation managers to make decisions. Due to the high acquisition cost and maintenance cost of train-sets, it is an effective measure to improve the utilization efficiency of train-sets. Thus, to complete the same number of trip tasks in a given train graph, less train-sets are needed. In a train-set utilization plan, the required amount of train-sets is equal to the value that 1440 minutes, i.e., a cycle, divides the sum of the total connecting time and trips' running time. The running time of all trips is a fixed value given in the train graph, so the required amount of train-sets is only determined by the total connecting time. Therefore, the optimization goal of train-set utilization plan is to minimize the total connecting time, which includes the following three parts: the connecting times of all connection arcs (

Model solution: A two-stage approach
The train-set utilization planning model aims to identify whether arcs should be connected between trips and maintenances. Among the three kinds of arcs, the connection arcs between two trips is relatively easy to obtain due to the static characteristics that all trips are given in the train graph. Inversely, when and where the maintenances should be carried out are indeterminate according to the accumulated operation time and mileages. Therefore, the difficulty in solving the train-set utilization planning problem lies in the solving process of maintenance arcs and arcs after maintenance. To solve the above difficulty, a two stage approach is proposed. As shown in Fig 10, in the first stage, all the trips in the given train graph are allocated into a set of segments and no arcs are connected between trip nodes and maintenance nodes.
In the second stage, the maintenance arcs and arcs behind maintenance will be connected between different segments. By designing a solving algorithm in each stage, the optimal and unique solution of the train-set utilization planning model can be obtained, and so the compilation of train-set utilization plan can be realized in a computer-generated system.

The first stage: Segment generation
With the introduction of the segment concept, a directed graphĜðV ;ÊÞ should be defined to represent the train-set utilization network. The node setV is a union of ordered subsetsV SE andV M , which represent segments and maintenances, respectively. Each node i 2V SE is defined as a tuple (p i ; The first stage aims to obtain the connection relationship between trip nodes and prepare for the next stage simultaneously. Therefore, when allocating trips into segments, the following requirements should be met.
1. For each segment, to ensure no maintenance arcs and arcs after maintenance are included, the running time and distance of all trips should be within the maintenance standards. Based on the above analyses, the segment generation algorithm is designed as follows.

The second stage: Segment connection
There are three decision variables in the train-set circulation planning model, including x C ij , x Q ik and x H ki . In the first stage, the decision variable x C ij has been solved since the connection relationships between trip nodes are obtained. Therefore, the objective of the second stage is to calculate the decision variables x Q ik and x H ki , i.e., the connection relationship between trips and break; end maintenance. With the concept of segment, the arcs between segments is equivalent to the combination of the maintenance arc and arc after maintenance. Therefore, the second stage aims to identify the segment connection relationship.
When node i 2V SE is connected to node j 2V SE , it means that the train-set should be maintained after completing the last trip in Segment i, and then undertake the first trip in Segment j. Therefore, to ensure that node i 2V SE can be connected to node j 2V SE , the following requirements should be satisfied: 1. The arrival station of the last trip, l s i in Segment i should be the same with the departure station of the first trip, f s i in Segment j.
2. The time interval, o S ij , between the last trip in Segment i and the first trip in Segment j should be longer than the duration of a series of working procedures, which is equal to the sum of T Z and T C . Particularly, the value of o S ij can be obtained by Eq (16).

Case studies
In this section, the case of the train graph of Beijing-Tianjin passenger dedicated line was adopted to test the performance of the proposed train-set utilization model and solution approach. Several evaluation indexes were applied to evaluate the quality of the computation results. Furthermore, the comparison between the ant colony algorithm and the two-stage approach was conducted to evaluate the advancements. The proposed approach is coded in C# and solved by Microsoft Visual Studio, and the computational experiments were performed on a 3.4 GHz Core i7 PC with 16 GB of RAM. The proposed model can be solved within a few seconds.
Main features of the case study
Step 2. Classify the segment setV SE according to the last trip arrival station, l s i ¼ st, and the first trip departure station, f s l s j ; l t j ; time j ; distance j Þg, sort the two sets by time ascending.
Step 3. Calculate o S ij according to the above rule, get the matrix {o S ij }.
Step 4. Gain a set of optimal solutions by the Hungarian Algorithm, get the matrix {o 0 S ij }, add the zero corresponding arc into set E.
Step 5. If st = m, end the algorithm; else let st = st + 1 and go to step 1.
• Duration of preparation work before train-set departing in station T C : T BJSS • Distance range of maintenance standard S one : 3, 600 kilometers S one 4, 400 kilometers.
• Time range of maintenance standard T one : 1.8 day T one 2.2 day. Computational results of the train-set utilization plan  Table 3 lists the results of some important evaluation indexes.
The results show that to complete the above 174 trip tasks, the minimal amount of needed train-set N is 12. The total maintenance times M of all the 12 train-sets is 6 each day. The average running distances of each train-set per day D average = D all / N is 1740 kilometers. Each day, the rail line overhaul is carried out during the period 0:00-6:00. Therefore, the train-sets cannot work for six hours per day, and the maximal running time of each train-set can be calculated as follows: 1440-6 Ã 60 = 1080 minutes. The average utilization efficiency of each trainset U = T all / (1080 Ã N) reaches 46.9%.
Theoretically, the minimum connection time of this case is obtained by applying Hungarian algorithm, which is equal to 25,542 minutes when maintenance constraints are not considered. Moreover, the total maintenance times of all train-sets is 4.7 each day, which can be calculated as D all / (1 + 10%)S one . However, the theoretical results cannot be reached in practical production, because train-sets must be maintained according to the rules and regulations. Also, not all train-sets need to be maintained until their running distances reaching the upper limit of the standard. Therefore, the connection time 31,055 minutes and maintenance times 6 obtained by the proposed model and approach are very close to the theoretical values. These results prove that the proposed model and approach are feasible and the train-sets are utilized with high efficiency.

Comparison with the ant colony algorithm
Many researchers have applied the ant colony algorithm (ACA) to solve train-set utilization problem. In order to highlight the performance of this study in improving the train-set utilization efficiency, the ant colony algorithm is applied to carry out the same experiments with the same computer, and the detailed comparison results are as listed in Table 4. The results show that obviously the two-stage approach has the advantages over the traditional ACA in improving the solution quality and shortening the computation time.
On one hand, the utilization efficiency of train-sets is much higher by applying the twostage approach. First, the amount of train-sets in use is reduced, which will save the acquisition cost significantly. Second, the total maintenance times is decreased from eight to six, saving much invalid running time from the railway stations to the inspection and repair depot. Third, the average running distances of each train-set per day is 134 kilometers longer than that of the traditional ACA.
On the other hand, the computation time of the two-stage approach is 2302 milliseconds, which is almost one-third of the traditional ACA. What should be noted is that, as the convergence speed of the traditional ACA is indeterminate, an average of 20 computation times is used. Furthermore, the computation of the two-stage approach is much more stable. Only one optimal train-set utilization plan can be obtained with two-stage approach while different results are achieved by using the traditional ACA. There are certain scope and limitations in the study. Firstly, the proposed approach is applicable for a short-term train-set circulation plan formulation. In China, there are five levels of maintenances in total, and the maintenance of level one is considered only in the study. As the maintenance of Level 2 to Level 5 take more than one month usually, it is beyond the scope of consideration in the study. Secondly, only one type of train-sets is studied by using the proposed approach, so the approach is appropriate for the new-built HSR with single type trainset, instead of the multi-type train-sets situation. As China has developed standard train-sets and tried to use the standard train-set as the main carrier in the future, the present approach will show important referential meaning in the utilization of standard train-set in the future.

Conclusion
The train-set utilization plan is drawn up to identify the work arrangements of train-sets, which is profoundly affected by the given train graph as well as the rules and regulations of maintenance. To optimize the train-sets' utilization efficiency, in this study, an integer programming model is proposed considering the spatial constraints, the time constraints, the maintenance constraints, and the unicity constraints. In the process of model solution, a twostage approach is designed and an optimal train-set utilization plan is obtained as the output. Using the model, the case of the real train graph of Beijing-Tianjin passenger dedicated line was carried out. The results show that the model and approach proposed in this study is practical and reasonable. By comparing the results obtained by the two-stage approach and the ant colony algorithm, it could be found that to complete the same trip tasks, the former one needs less train-sets, less maintenance times, and longer average running distances for each train-set per day. In short, the utilization efficiency of train-sets is much higher by applying the twostage approach. Moreover, the fast computation speed and good stability can contribute to the realization of a train-set utilization plan computer-generated system.