Purpose

This study investigates the impact of flagship trains on high-speed railway capacity utilization and develops a brand value-oriented optimization framework that balances service quality enhancement with operational efficiency.

Design/methodology/approach

A mathematical optimization model based on integer programming is developed, incorporating flagship train constraints into capacity optimization. Case studies compare scenarios with and without flagship train considerations using the Beijing–Shanghai High-Speed Railway data across 20 experimental groups.

Findings

Operating flagship trains with hourly departure constraints results in an average decrease of 0.9 trains and an 8.4% reduction in capacity utilization rate. When scheduling 2 flagship trains within a 2-h timeframe, capacity utilization decreases from 86.43% to 83.73%, quantifying the trade-off between brand positioning and operational capacity.

Originality/value

This research provides the first quantitative framework for brand value-oriented railway capacity optimization, establishing clear definitions for flagship trains and mathematical foundations for evaluating service quality versus efficiency trade-offs. The findings offer practical decision support for railway operators balancing competitive positioning with capacity maximization.

The rapid development of high-speed railway systems worldwide has fundamentally transformed intercity transportation, offering unprecedented travel speeds and service quality. As competition intensifies among different transportation modes and passenger expectations continue to rise, railway operators face the dual challenge of maximizing infrastructure capacity utilization while enhancing service quality and brand value. This challenge is particularly pronounced in densely operated high-speed railway corridors, where the strategic scheduling of premium services can significantly impact both operational efficiency and market competitiveness. The concept of brand value in railway operations has emerged as a critical factor, influencing not only passenger choice behavior but also the fundamental approaches to capacity planning and timetable design. This research addresses the intricate relationship between brand-oriented service design and railway capacity utilization, proposing a novel framework that balances operational efficiency with service quality enhancement in high-speed railway systems.

Railway transportation optimization is a comprehensive optimization problem encompassing a series of optimization processes, including passenger flow demand, network planning, train operation scheme optimization, timetable optimization, rolling stock utilization optimization, and crew scheduling optimization. In railway operations, high-speed trains with high travel speeds, fewer stops, and serving intercity routes with large passenger volumes are commonly referred to as “flagship trains” by field practitioners. During the timetable compilation phase, China Railway Corporation currently adopts a phased approach, first scheduling same-line and cross-line flagship trains, followed by other trains. Flagship trains have gradually become premium services in high-speed railway transportation organization, possessing significant brand value. However, there is currently no clear definition for existing flagship trains.

To further quantify the impact of flagship trains on timetable compilation and capacity utilization, this study proposes the concept of “flagship trains” and investigates the issue of carrying capacity utilization from a brand value-oriented perspective. The operation of flagship trains is regarded as key to enhancing brand value. This study characterizes the number of train paths consumed by flagship trains under different transportation organization schemes, and further explores the impact of operating clock-face flagship trains—a service enhancement measure—on the actual capacity “deduction.” This research establishes non-value-oriented carrying capacity utilization as the foundation, and further develops a brand value-oriented high-speed railway carrying capacity utilization model. By comparing capacity optimization results under different train stopping pattern alternatives, case studies are conducted on the Beijing South-Jinan West section and the entire Beijing-Shanghai High-Speed Railway line. The results reveal the average changes in capacity indicators such as train numbers and capacity utilization rates under the brand value-oriented constraint of fixed clock-face departures for flagship trains.

Due to their characteristics of fewer stops, higher average speed, superior service quality, operation on high-speed grade lines, and departure times at regular hourly intervals that are convenient for passengers to remember and recognize, flagship trains have become the preferred travel option for some passengers. With their higher occupancy rates, they should possess significant operational value. However, their scheduling requires sacrificing the number of operating train paths, which also has a notable impact on capacity utilization.

Lu (2019) analyzed the number of non-stop trains deducted (deduction coefficient) on high-speed railway parallel operation diagrams from the perspective of flagship train overtaking characteristics, and based on the analysis results, proposed relevant suggestions for the Beijing-Shanghai high-speed railway timetable planning. Zhang, Fan, Wang, and Guo (2022) considered the periodic train operation structure and utilized the close relationship between section passing capacity and train stops to optimize the number of flagship trains and their travel speeds. Xia, Li, Zhang, and Zhang (2024) constructed an integer linear programming model based on Event Activity Network (EAN) with the objectives of OD service balance and capacity utilization. The model considered basic timetabling constraints such as flagship train section running time and connection time, combined with special constraints including flagship train departure time constraints and grouping constraints, providing reference for high-speed railway flagship train timetable compilation and optimization.

Qu, Shi, and Xu (2022) explored the increase in the number of high-speed railway flagship trains, adopting an hourly balance-centralized scheduling approach to explore improvements in flagship train service quality. When analyzing the factors affecting the passing capacity of corridor-type high-speed railways in existing train timetables, they identified the pattern that flagship trains generally depart at regular hourly intervals to improve service quality. Fan, Shen, Liu, Li, and Liang (2024) analyzed from a systematic perspective, simplifying train path efficiency as the ratio of technical indicators to invested transportation resources to objectively and effectively evaluate the efficiency of various types of train paths. Through DEA projection analysis, they provided targeted suggestions for improving train path efficiency. The results showed that the comprehensive efficiency of trains adopting high-speed grade operating standards was higher than other categories of trains in the entire timetable.

Mei, Lai, Zhang, Pan, and Cui (2024) examined flagship trains' operational characteristics through comprehensive evaluation metrics. Their analysis using the improved TOPSIS method revealed that flagship trains such as G4, G3, and G8 ranked among the top performers on the Beijing-Shanghai HSR, characterized by high speed, high grade, high mileage rates, and fewer stops. Sun, Jing, Li, Tian, and Tian (2024) proposed a method for predicting passenger flow distribution in high-speed railway train riding sections based on an interpretable machine learning framework to explain the impact of passenger service product characteristics on passenger flow distribution prediction. They focused on analyzing the relationships between features such as benchmark trains, ticket prices, travel time, date, day of the week, train number, departure time, and passenger flow. The results showed that date and benchmark trains exhibited a nonlinear positive correlation with passenger flow distribution prediction, while travel time showed a nonlinear negative correlation. Benchmark trains with low travel time, high ticket prices, and departure times at regular hourly intervals had a positive impact on passenger flow distribution prediction, with benchmark trains having a significant influence on both passenger flow and revenue. Wu, Shan, and Zhao (2025) identified that flagship trains exhibit relatively low stop ratios in their optimization model for the Beijing-Shanghai high-speed railway, as these trains primarily serve longer travel distances and high-level stations. Their analysis showed that a small number of flagship trains maintain stop ratios below 0.6, distinguishing them from other train types that make more frequent stops at county-level stations to serve local passenger demand.

In other studies regarding timetables, capacity, and service elements, Jensen, Schmidt, and Nielsen (2020) addressed railway infrastructure capacity determination for strategic long-term planning without requiring fixed timetables. They formulated capacity determination as an optimization problem to find the maximal train set that minimizes average waiting time between trains of the same type, testing feasibility through train sequence sampling rather than fixed timetable generation. Zhong, Zang, Yang, Sheng, and Xu (2024) investigated the impact of timetable elements on railway line capacity, establishing the relationship between consumed capacity and timetable elements as the signed sum along a “critical path.” They demonstrated that changing timetable elements impacts capacity with neighboring trains in a combinatorial manner, analyzing stop plans and extending findings to capacity factors including average speed and heterogeneity. Pouryousef and Lautala (2015) investigated applying European timetable management methodologies to the less structured U.S. rail system through a hybrid simulation approach combining RTC and RailSys tools. They demonstrated that optimizing dwell time constraints can significantly improve level of service parameters while managing the trade-off with capacity utilization, showing the potential benefits of structured timetable management in traditionally unstructured rail networks. Gu, Wei, and Xia (2025) transformed the train timetable scheduling problem into a node routing problem where nodes represent train paths and directed arcs depict sequential connections between paths. They constructed an optimization model maximizing the number of trains under the constraint that total arc weight does not exceed available time-space resources, revealing how train path sequences impact section carrying capacity and proposing optimization strategies for improving capacity utilization of high-speed railways. Xu et al. (2023) addressed railway capacity optimization through quantum computing by developing a space-time network-based binary integer programming model for train timetabling. They demonstrated that quantum simulators, particularly the coherent Ising machine (CIM), can effectively solve capacity-constrained timetabling problems, offering a novel computational approach for maximizing line capacity utilization in complex high-speed rail networks where traditional methods face computational limitations. Gao, Chen, Xu, Liu, and Chen (2025) developed a bi-objective optimization model for high-speed train rescheduling under disruptions that enhances capacity recovery and service quality by minimizing both timetable deviation and impacts on transferring passengers. Their computational intelligence framework achieved significant service improvements by preserving line capacity and service levels during disruptions. D'Ariano, Pacciarelli, and Pranzo (2008) introduced flexible timetables as an alternative to traditional time reserves for improving railway punctuality without reducing line capacity. Their approach plans less rigidity in baseline schedules while empowering real-time management to resolve inter-train conflicts dynamically. Jensen, Landex, Nielsen, Kroon, and Schmidt (2017) developed a strategic planning framework for calculating railway infrastructure occupation and capacity consumption independent of specific timetables. Their model generates and assesses different train sequences to achieve timetable independence, employing stochastic delay simulation to determine capacity consumption.

Existing research has not clearly defined the quantitative relationships between service levels for different types of trains, including flagship trains, cross-line trains, and other mainline trains. Research on flagship train timetable compilation methods is relatively limited overall. Existing studies indicate that operating flagship trains is an effective way to improve service levels, with their establishment being driven by both the necessity and public welfare nature of macro policies on one hand, and the convenience and comfort of passenger travel on the other. However, in existing research, there is limited analysis of the impact of flagship trains on capacity, making it difficult to provide effective analysis for capacity utilization.

Before discussing flagship train classification, this study first explores the characteristics and connotations of flagship trains. In railway operations, high-speed trains with high travel speeds, fewer stops, and serving high-traffic intercity routes are commonly referred to as flagship trains. Train timetabling involving flagship trains represents a special category with distinct objectives and constraints compared to conventional timetabling methods.

Flagship trains feature fast travel speeds, fewer stops, operation at maximum speed standards, stops mainly at key provincial capitals, and typically whole-hour departures. Taking the Beijing-Shanghai high-speed line in Q4 2022 as an example, five of the 28 trains (G1-G28) depart at non-whole-hour times. Through analyzing current flagship train data, existing research has found it difficult to identify common characteristics. This paper identifies the following key issues:

  1. Consideration of whole-hour departures

While most whole-hour trains on the Beijing-Shanghai line are flagship trains, it remains unclear whether non-whole-hour departure trains can be flagship trains and whether whole-hour departure necessarily defines a flagship train.

  1. Consideration of stops

Flagship trains on the Beijing-Shanghai line stop at 2–5 stations, but the specific threshold for “few stops” among 23 stations remains undefined. Whether trains stopping at 6–7 stations qualify as flagship trains is debatable. While flagship trains primarily serve provincial capitals, some also stop at high-traffic stations like Xuzhou East, making stopping standards difficult to define.

  1. Consideration of travel speed

Whether “flagship trains” should exist among lower-grade trains operating below maximum speed, and whether lower-grade trains with whole-hour departures and few stops qualify as flagship trains requires further discussion.

Given these conceptual uncertainties, this study proposes the concept of flagship trains and examines their impact on capacity utilization.

This study defines flagship trains as a special category of flagship trains that: depart at whole hours, operate at maximum speed grade, stop at no more than one-quarter of total line stations, and serve major prefecture-level cities with high passenger flow. Key characteristics include short travel times and balanced whole-hour departures, which are crucial for optimizing capacity utilization in high-speed railway systems.

  1. Short travel time

All flagship trains operate at 350 km/h with limited stops (≤25% of line stations), enabling higher speeds and reduced stopping time, thereby improving passenger travel efficiency and satisfaction.

  1. Balanced whole-hour departures

Whole-hour departure times provide regularity and predictability, facilitating passenger planning while enhancing reliability and optimizing transportation organization efficiency.

Flagship trains are primarily characterized as flagship trains—high-quality trains with strict timetabling requirements that serve as reference standards for timetable optimization.

Based on operating characteristics, flagship trains can be classified as:

  1. Different stopping patterns

Flagship trains can be divided into direct and limited-stop services. Different stopping patterns help evaluate operational performance by comparing travel time, station spacing, and passenger demand satisfaction to select optimal stopping schemes.

  1. Different time periods

Peak and off-peak flagship trains may have different timetables and stopping schemes to accommodate varying passenger demand and traffic flow.

This classification helps planners evaluate different train schemes to improve railway efficiency and service quality.

As reference standards in high-speed railway systems, flagship trains possess significant value:

  1. Brand value

Flagship trains provide high-quality transportation services that meet rapid travel demands, promote economic exchange, and drive regional integration. They attract passenger flow from other modes and serve as the flagship brand of high-speed rail products.

  1. Economic value

Short travel times reduce time costs and promote business/tourism development. Whole-hour departures improve reliability and capacity utilization while attracting more passengers, increasing revenue. Flagship trains enhance economic benefits through improved efficiency and optimized organization.

  1. Structural value

In timetables, flagship trains exhibit hourly pseudo-periodic patterns that provide structural value for overall timetable layout and compilation, significantly impacting capacity optimization.

The value of flagship trains encompasses brand significance, economic benefits, and structural contributions, providing efficient travel options while serving as references for system planning and management to improve timetable efficiency and capacity utilization.

This chapter addresses the throughput capacity utilization problem for high-speed railways considering the operation of flagship trains. Flagship trains, characterized by their fixed departure times at round hours and premium service levels, represent a critical component of brand value in railway operations. The model optimizes the selection and scheduling of trains from a candidate set while ensuring flagship trains operate at their designated hourly departure times.

The following assumptions are made for the capacity utilization model:

  1. Known candidate train set: The candidate train set is predetermined and serves as input to the optimization problem. The set includes trains of different speed grades (300km/h and 350km/h) with various stopping patterns, where flagship trains represent a premium subset with mandatory whole-hour departures and limited stops. Each candidate train is defined by its origin-destination pair, stopping scheme, and service characteristics.

  2. Sufficient rolling stock and crew resources: The model assumes adequate availability of train units and crew, without explicitly optimizing their allocation.

  3. Train-kilometer conversion coefficient as objective: The optimization uses weighted train-kilometer conversion coefficients instead of simple train counts, allowing proper consideration of short-haul services converted to long-haul equivalent capacity.

As shown in Figure 1, train b has the longest running distance. Without considering other factors, the train-kilometer coefficient of train b should be 1, while the train-kilometer coefficients of train a and train c are converted to 0.4 and 0.6 respectively.

Figure 1
A time-space diagram for train scheduling shows the movement of three trains between three stations.The time-space diagram shows the movement of three trains across three stations. The vertical axis is labeled “Stations,” showing “Station i” at the bottom, “Station i plus 1” in the middle, and “Station i plus 2” at the top. The distance between “Station i” and “Station i plus 1” is marked “6 kilometers.” The distance between “Station i plus 1” and “Station i plus 2” is marked “4 kilometers.” The horizontal axis is labeled “Time.” A green rectangular grid forms the plot background, with vertical lines and a horizontal green line dividing the area. Three trains labeled “Train a,” “Train b,” and “Train c” are shown as diagonal lines moving upward and to the right. “Train a” travels from “Station i plus 1” to “Station i plus 2,” covering 4 kilometers. “Train b” passes between “Station i” and “Station i plus 2,” covering 10 kilometers. “Train c” moves from “Station i” to “Station i plus 1,” covering 6 kilometers.

Illustration of the train-kilometer conversion coefficient. Source: Author’s own work

Figure 1
A time-space diagram for train scheduling shows the movement of three trains between three stations.The time-space diagram shows the movement of three trains across three stations. The vertical axis is labeled “Stations,” showing “Station i” at the bottom, “Station i plus 1” in the middle, and “Station i plus 2” at the top. The distance between “Station i” and “Station i plus 1” is marked “6 kilometers.” The distance between “Station i plus 1” and “Station i plus 2” is marked “4 kilometers.” The horizontal axis is labeled “Time.” A green rectangular grid forms the plot background, with vertical lines and a horizontal green line dividing the area. Three trains labeled “Train a,” “Train b,” and “Train c” are shown as diagonal lines moving upward and to the right. “Train a” travels from “Station i plus 1” to “Station i plus 2,” covering 4 kilometers. “Train b” passes between “Station i” and “Station i plus 2,” covering 10 kilometers. “Train c” moves from “Station i” to “Station i plus 1,” covering 6 kilometers.

Illustration of the train-kilometer conversion coefficient. Source: Author’s own work

Close modal

The sets, indices, parameters, and decision variables used in the model are shown in Table 1 and Table 2.

Table 1

Set, indexes, parameters of the model

SymbolsDefinitions
NSet of stations on the line
KSet of trains
KnSet of low-speed grade trains (300km/h) KnK
KbSet of high-speed grade trains (350km/h) KbK
KjSet of flagship trains, KjKb
i,jStation index, i,jN
mTotal number of trains running on the line
k,kTrain index, k,kK
NkkSet of stations commonly passed by trains k and k, k,kK
okOrigin station of train k
dkTerminal station of train k
STkiStopping status of train k at station i, 0 means passing without stopping, 1 means stopping
FiMinimum number of stops required at station i within the specified time
Tri,i+1kPure running time of train k in section (i,i+1)
tikDeparture time of train k, generally taken as whole hours, i.e. 420min, 480min, etc., tikTI
SkiStarting additional time of train k at station i
BkiStopping additional time of train k at station i
TsiminMinimum dwell time of trains at stations
TsimaxMaximum dwell time of trains at stations
TiminUpper limit of operation time window at station i
TimaxLower limit of operation time window at station i
IddMinimum train departure interval
IdpMaximum train passing interval
IpdMinimum passing-to-departure interval
IaaMinimum arrival-to-arrival interval
IapMinimum arrival-to-passing interval
IpaMinimum passing-to-arrival interval
IppMinimum passing-to-passing interval
MA large positive number
LsvkTrain-kilometer conversion coefficient for train k
TISet of whole hour times
Source(s): Author’s own work
Table 2

Variables of the model

SymbolsDefinitions
zkBinary variable, indicating whether train k is operated. If operated, zk=1; otherwisezk=0
aikInteger variable, the arrival time of train k at station i
dikInteger variable, the departure time of train k at station i
yik,kBinary variable, indicating the running order of trains k and k in section (i1,i). If yik,k=1, train k is the preceding train of train k; if yik,k=0, train k is the following train of train k
Source(s): Author’s own work

The capacity utilization optimization model considering flagship trains proposed in this section is an integer programming model based on transportation resource request conflicts. The model is constructed using the Big M method, which has certain similarities with the capacity calculation model based on train diagram optimization graphical methods. The optimization objective is adjusted from the number of trains to maximizing the weighted train-kilometer conversion coefficient, that is, establishing the optimization objective function by maximizing the sum of weighted train-kilometer conversion coefficients of selected trains to be operated.

(1)

The model includes the following constraints:

  1. Station Service Frequency Constraint

To ensure passengers' basic travel needs are met, it is necessary to limit the total number of trains stops at each station. The minimum service frequency at each station should be set based on passenger demand, using either passenger boarding/alighting numbers or passenger travel demand as approximation methods.

(2)
  1. Train Dwell Time Constraint

To ensure smooth passenger boarding/alighting and necessary technical operations, minimum dwell times must be established for trains at each station. To minimize negative impacts on passenger experience, train dwell times should also not be excessive. Therefore, both minimum and maximum time limits are necessary. Since train k's operation is uncertain, both sides of the equation are multiplied by zk.

(3)
(4)
(5)

Due to the binary decision variable for train operation, the constraints are linearized using the Big-M method:

(6)
(7)
(8)
(9)
  1. Inter-station Running Time Constraint

The inter-station running time includes the travel time between adjacent stations. For train k in section (i,i+1), the running time equals the pure running time plus the starting additional time at station i and the stopping additional time at station (i+1).

(10)

Similarly, linearized using the Big-M method:

(11)
(12)
  1. Train Departure Interval Constraint

To ensure safe operating distances between trains, adjacent trains must maintain specified minimum departure intervals. These intervals are affected by the stopping strategies of adjacent trains at stations. Four types of intervals must be considered: minimum stop-to-stop departure interval Idd (preceding train stops, following train stops); minimum stop-to-pass departure interval Iap (preceding train stops, following train passes); minimum pass-to-stop departure interval Ipa (preceding train passes, following train stops); minimum pass-to-pass departure interval Ipp (both trains pass). The specific interval type depends on both trains' stopping decisions at the station and their sequence order in the section.

When both trains k and k operate and stop at station i: if yi+1k,k=1 (train k precedes k in section (i,i+1)), equation (13) applies; if yi+1k,k=0 (train k follows k), equation (14) applies.

(13)
(14)

When train k stops at station i while k passes:

(15)
(16)

When both trains pass station i:

(17)
(18)
  1. Train Arrival Interval Constraint

Similar to departure intervals, adjacent trains must maintain specified arrival intervals at stations. Four types of arrival intervals are considered: Iaa, Iap, Ipa, and Ipp.

(19)
(20)
(21)
(22)
(23)
(24)
  1. Time Window Constraint

To ensure safe train operations on high-speed railways, time window constraints such as maintenance windows must be established. This study uses parallelogram (diamond) shaped time windows, where each train's arrival and departure times at each station must fall within the time window.

(25)
(26)
  1. Section Running Order Constraint

This constraint defines the running order of any two trains in a specific section. On the same track section, one train must be the preceding train while the other follows.

(27)
  1. Flagship Train Hourly Departure Constraint

Flagship trains provide new approaches for brand value-oriented capacity optimization. All flagship trains belong to high-grade train categories and must depart at the top of the hour.

(28)
(29)
  1. Decision Variable Constraints

(30)
(31)

The Beijing-Shanghai High-Speed Railway is one of the busiest high-speed rail lines in China, comprising 23 stations in total. This case study focuses on the southbound Beijing-Jinan section (from Beijing South Station to Jinan West Station), which includes 6 stations. The time range selected is a two-hour daytime period from 7:30 to 9:30.

The minimum dwell time for trains at intermediate stations is set at 2 minutes, with a maximum of 10 minutes. The departure headway between trains at each station is 5 minutes, the arrival headway is 4 minutes, the passing headway is 3 minutes, the passing-departure headway is 3 minutes, the departure-passing headway is 5 minutes, the arrival-passing headway is 4 minutes, and the passing-arrival headway is 3 minutes.

Two speed categories of trains operate on this line: high-speed trains (350 km/h) and low-speed trains (300 km/h). The running times of trains with different speed levels in each section are shown in Table 3.

Table 3

Jingji section line and parameter information

StationMileage (km)Section pure running time (min)Additional time for train starting/topping (min)Min service frequency
350km/h300km/hStartingStopping
BJN02310
LF591516231
TJN1311113235
CZX2191618233
DZD3271820235
JNX41917192310
Source(s): Author’s own work

The models proposed in this paper are implemented using Python 3.11 with the Branch-and-Cut algorithm. The computing environment consists of a personal computer with an Intel Core i7-11700 CPU @ 2.50 GHz and 32 GB of RAM. The construction of candidate train sets follows the approach of Zhang et al. (2022), whose research was based on a quarterly timetable from 2022. After providing a basic set of candidate trains, the candidate stopping patterns were continuously supplemented, using a 2-h time range as the research example. This case study establishes 20 different candidate train sets, each maintaining a fixed number of 22 candidate trains without further increases. Among these, 3 are high-speed train candidates and 19 are low-speed train candidates. The main difference between each candidate train set lies in the low-speed train candidate sets, which are configured based on different schemes derived from actual operational timetables. The specific generation rules are shown in Figure 2, with low-speed train stopping patterns detailed in Appendix A.

Figure 2
A visual of train stop patterns restructured into experimental schemes using selected low-grade trains.The train stopping scheme diagram visually represents the stopping patterns of trains at various stations along a route, divided into two main sections. The left section is labeled “Low-grade actual train stopping scheme.” It lists “G 101,” “G 103,” “G 105,” “G 107,” and “G 901” with horizontal lines representing their routes. Each line includes several circular stops. The middle of each line has a small open circle, and each route extends horizontally from left to right. The right section is labeled “Form 20 sets of alternative experimental schemes.” It shows two sections. Each section contains train routes labeled with stations “B J N,” “L F,” “T J N,” “C Z X,” “D Z D,” and “J N X” along the top. Horizontal lines with circular nodes represent different train stopping patterns. Inside both sections, these schemes are divided into two: “Low grade” and “High grade.” High-grade points are placed inside the red box, and low-grade points are placed inside the green box. The top section is labeled “Train stopping alternative set Scheme 1.” High grade at top, arrows pointing from left to right: 3, B J N, J N X. Low grade, starting from top to bottom, and arrows pointing from left to right: 1, B J N, L F, T J N, D Z D, J N X. 1, B J N, L F, T J N, J N X. 1, B J N, L F, C Z X, J N X. 1, B J N, T J N, D Z D, J N X. 3, B J N, C Z X, D Z D, J N X. 3, B J N, C Z X, J N X. 2, B J N, D Z D, J N X. 1, B J N, L F, T J N, C Z X, D Z D, J N X. 2, B J N, D Z D. 2, B J N, T J N, D Z D. 1, B J N, C Z X, D Z D. 1, T J N, C Z X, D Z D, J N X. The bottom section is labeled “Train stopping alternative set Scheme 20.” High grade at top, arrows pointing from left to right: 3, B J N, J N X. Low grade, starting from top to bottom, and arrows pointing from left to right: 2, B J N, L F, C Z X, J N X. 2, B J N, T J N, C Z X, J N X. 2, B J N, T J N, D Z D, J N X. 2, B J N, C Z X, D Z D, J N X. 2, B J N, C Z X, J N X. 2, B J N, D Z D, J N X. 2, B J N, J N X. 1, B J N, C Z X, D Z D. 1, B J N, D Z D. 1, T J N, C Z X, D Z D, J N X. A blue arrow is labeled “Select 19 trains from low-grade actual train,” and an orange arrow is labeled “Select another 19 trains from low-grade actual train.” Four blue arrows are pointing from the left section to the right section's top, labeled “Train stopping alternative set Scheme 1.” Four orange arrows are pointing from the left section to the right section's bottom, labeled “Train stopping alternative set Scheme 20.”

Candidate train set generation rules. Source: Author’s own work

Figure 2
A visual of train stop patterns restructured into experimental schemes using selected low-grade trains.The train stopping scheme diagram visually represents the stopping patterns of trains at various stations along a route, divided into two main sections. The left section is labeled “Low-grade actual train stopping scheme.” It lists “G 101,” “G 103,” “G 105,” “G 107,” and “G 901” with horizontal lines representing their routes. Each line includes several circular stops. The middle of each line has a small open circle, and each route extends horizontally from left to right. The right section is labeled “Form 20 sets of alternative experimental schemes.” It shows two sections. Each section contains train routes labeled with stations “B J N,” “L F,” “T J N,” “C Z X,” “D Z D,” and “J N X” along the top. Horizontal lines with circular nodes represent different train stopping patterns. Inside both sections, these schemes are divided into two: “Low grade” and “High grade.” High-grade points are placed inside the red box, and low-grade points are placed inside the green box. The top section is labeled “Train stopping alternative set Scheme 1.” High grade at top, arrows pointing from left to right: 3, B J N, J N X. Low grade, starting from top to bottom, and arrows pointing from left to right: 1, B J N, L F, T J N, D Z D, J N X. 1, B J N, L F, T J N, J N X. 1, B J N, L F, C Z X, J N X. 1, B J N, T J N, D Z D, J N X. 3, B J N, C Z X, D Z D, J N X. 3, B J N, C Z X, J N X. 2, B J N, D Z D, J N X. 1, B J N, L F, T J N, C Z X, D Z D, J N X. 2, B J N, D Z D. 2, B J N, T J N, D Z D. 1, B J N, C Z X, D Z D. 1, T J N, C Z X, D Z D, J N X. The bottom section is labeled “Train stopping alternative set Scheme 20.” High grade at top, arrows pointing from left to right: 3, B J N, J N X. Low grade, starting from top to bottom, and arrows pointing from left to right: 2, B J N, L F, C Z X, J N X. 2, B J N, T J N, C Z X, J N X. 2, B J N, T J N, D Z D, J N X. 2, B J N, C Z X, D Z D, J N X. 2, B J N, C Z X, J N X. 2, B J N, D Z D, J N X. 2, B J N, J N X. 1, B J N, C Z X, D Z D. 1, B J N, D Z D. 1, T J N, C Z X, D Z D, J N X. A blue arrow is labeled “Select 19 trains from low-grade actual train,” and an orange arrow is labeled “Select another 19 trains from low-grade actual train.” Four blue arrows are pointing from the left section to the right section's top, labeled “Train stopping alternative set Scheme 1.” Four orange arrows are pointing from the left section to the right section's bottom, labeled “Train stopping alternative set Scheme 20.”

Candidate train set generation rules. Source: Author’s own work

Close modal

This study examines 3 high-speed trains, with 2 designated as flagship trains scheduled to depart on the hour. A capacity optimization case considering flagship trains (FO) is designed, with 20 experimental groups named FO1, FO2, ..., FO20, to further analyze the capacity results.

Under the condition of scheduling 2 flagship trains, the optimization results are shown in Figure 3. The figure illustrates that two high-speed trains operate as flagship trains departing at 8:00 and 9:00 sharp. In subfigures 3(1), (3), and (6), the third high-speed train is not selected for operation. The specific solution results for case FO are presented in Table 4.

Figure 3
A panel of six graphs shows multiple comparative train speed data trends.Six graphs are arranged in a 3 by 2 grid. The graphs labeled from top left “(1) Result of case F O 1,” “(2) Result of case F O 2,” “(3) Result of case F O 3,” “(4) Result of case F O 4,” “(5) Result of case F O 5,” “(6) Result of case F O 6,” Each row contains two graphs. In all graphs, the vertical axis is labeled with categories from bottom to top as follows: “J N X,” “D Z D,” “C Z N,” “L F N,” “L F,” and “B J N,” and the horizontal axis ranges from “7:00” to “12:00,” with increments of 1:00 hour. In all cases, multiple pink and blue lines slope diagonally downward from the upper left to the lower center and lower right sides. The blue lines are inside and overlapping the pink lines. In graphs (1), (3), and (6), 2 blue lines overlap the pink lines. In graphs (2), (4), and (5), 3 blue lines overlap the pink lines. Each graph begins with lines at the top in the range of “7:30” to “9:30” and ends with lines toward the bottom in the range of “9:10” to “11:10” across all six plots. Note: All numerical data values are approximated.

Timetable example of the Case FO. Source: Author’s own work

Figure 3
A panel of six graphs shows multiple comparative train speed data trends.Six graphs are arranged in a 3 by 2 grid. The graphs labeled from top left “(1) Result of case F O 1,” “(2) Result of case F O 2,” “(3) Result of case F O 3,” “(4) Result of case F O 4,” “(5) Result of case F O 5,” “(6) Result of case F O 6,” Each row contains two graphs. In all graphs, the vertical axis is labeled with categories from bottom to top as follows: “J N X,” “D Z D,” “C Z N,” “L F N,” “L F,” and “B J N,” and the horizontal axis ranges from “7:00” to “12:00,” with increments of 1:00 hour. In all cases, multiple pink and blue lines slope diagonally downward from the upper left to the lower center and lower right sides. The blue lines are inside and overlapping the pink lines. In graphs (1), (3), and (6), 2 blue lines overlap the pink lines. In graphs (2), (4), and (5), 3 blue lines overlap the pink lines. Each graph begins with lines at the top in the range of “7:30” to “9:30” and ends with lines toward the bottom in the range of “9:10” to “11:10” across all six plots. Note: All numerical data values are approximated.

Timetable example of the Case FO. Source: Author’s own work

Close modal
Table 4

Analysis of the results of the case FO

CaseNumber of trainsNumber of high-priority trainsNumber of low-priority trainsNumber of constraintsNumber of variablesObjective functionSolution time (s)
FO120/222/318/1910,5502,58218.71212.24
FO219/223/316/199,6412,57217.94246.67
FO320/222/318/1910,7342,58419.17218.19
FO418/223/315/1910,7632,58417.77193.42
FO520/223/317/1910,3672,58018.78186.38
FO619/222/317/1911,5292,59218.7238.88
FO720/223/317/199,6412,57218.71208.99
FO819/223/316/199,9972,57618.24204.35
FO919/223/316/1910,7342,58418.47262.39
FO1020/222/318/1910,5462,58218.94250.49
FO1119/223/316/1910,5892,58218.24226.73
FO1219/223/316/1911,1262,58818.47250.7
FO1319/223/316/1910,1762,57818.4256.71
FO1419/222/317/1911,1372,58818.4215.31
FO1518/222/316/1911,3482,59018252.84
FO1619/223/316/1910,5462,58218.47248.43
FO1718/223/315/1910,5892,58217.47189.02
FO1819/223/316/1910,7382,58418.47242.74
FO1919/223/316/1910,5612,58218.7279.33
FO2019/223/316/1910,1762,57818.4254.33
Source(s): Author’s own work

To compare the optimization results with those considering flagship trains, a capacity utilization optimization model without flagship train considerations (NFO) is designed using the same candidate train sets. Twenty experimental groups are conducted, named NFO1, NFO2, ..., NFO20, to further analyze the capacity results. The optimization outcomes are shown in Figure 4.

Figure 4
Six comparative graphs plot red and blue time events across seven labeled rows.Six graphs are arranged in a 3 by 2 grid. The graphs are labeled, from the top left, “(1) Result of case N F O 1,” “(2) Result of case N F O 1,” “(3) Result of case N F O 1,” “(4) Result of case N F O 4,” “(5) Result of case N F O 5,” “(6) Result of case N F O 6” two graphs in each row. In all graphs, the vertical axis is labeled with categories from bottom to top as follows: “J N X,” “D Z D,” “C Z N,” “L F N,” “L F,” and “B I N,” and the horizontal axis plots time from “7:00” to “12:00,” with increments of 1:00 hour. Across all graphs, multiple pink and blue diagonal lines slope downward from left to right. The blue lines are positioned overlapping the pink lines. In graphs (1), (5), and (6), 2 blue lines overlap the pink lines. In graphs (2) and (3), 3 blue lines overlap the pink lines. In graph (4), 1 blue line overlaps the pink lines. Each graph begins with lines at the top in the range of “7:30” to “9:30” and ends with lines toward the bottom in the range of “9:10” to “11:10” across all six plots. Note: All numerical data values are approximated.

Timetable example of the Case NFO. Source: Author’s own work

Figure 4
Six comparative graphs plot red and blue time events across seven labeled rows.Six graphs are arranged in a 3 by 2 grid. The graphs are labeled, from the top left, “(1) Result of case N F O 1,” “(2) Result of case N F O 1,” “(3) Result of case N F O 1,” “(4) Result of case N F O 4,” “(5) Result of case N F O 5,” “(6) Result of case N F O 6” two graphs in each row. In all graphs, the vertical axis is labeled with categories from bottom to top as follows: “J N X,” “D Z D,” “C Z N,” “L F N,” “L F,” and “B I N,” and the horizontal axis plots time from “7:00” to “12:00,” with increments of 1:00 hour. Across all graphs, multiple pink and blue diagonal lines slope downward from left to right. The blue lines are positioned overlapping the pink lines. In graphs (1), (5), and (6), 2 blue lines overlap the pink lines. In graphs (2) and (3), 3 blue lines overlap the pink lines. In graph (4), 1 blue line overlaps the pink lines. Each graph begins with lines at the top in the range of “7:30” to “9:30” and ends with lines toward the bottom in the range of “9:10” to “11:10” across all six plots. Note: All numerical data values are approximated.

Timetable example of the Case NFO. Source: Author’s own work

Close modal

In the capacity utilization optimization case NFO without flagship train considerations, trains are scheduled within the parallelogram time window (2-h pink area) with the objective of maximizing the weighted train-kilometer conversion coefficient. Analysis reveals that in Figure 4(b) and (c), high-speed trains appear in “paired” tracking patterns, which is a common approach demonstrating reasonable capacity utilization. In Figure 4(a), (d), (e), and (f), high-speed and low-speed trains make reasonable use of spatial-temporal resources through stopping patterns and overtaking arrangements. These cases typically feature more low-speed trains with high-speed trains appearing in a dispersed manner. Detailed calculation results are shown in Table 5.

Table 5

Analysis of the results of the case NFO

CaseNumber of trainsNumber of high-priority trainsNumber of low-priority trainsNumber of constraintsNumber of variablesObjective functionSolution time (s)
NFO120/222/318/1910,5462,58218.94232.43
NFO220/222/318/199,6372,57218.41342.7
NFO320/222/318/1910,7302,58419.17225.25
NFO420/222/318/1910,7592,58419.01201.2
NFO520/223/317/1910,3632,58019.01266.67
NFO619/221/318/1911,5252,59218.7183.19
NFO722/223/319/199,6372,57219.8889.52
NFO820/222/318/199,9932,57618.71216.92
NFO920/222/318/1910,7302,58419.17211.38
NFO1020/222/318/1910,5422,58219.01223.35
NFO1120/223/317/1910,5852,58219.01240.16
NFO1219/222/317/1911,1222,58818.7187.4
NFO1320/222/318/1910,1722,57818.64252.18
NFO1419/220/319/1911,1332,58818.4276.3
NFO1520/221/319/1911,3442,59019.31224.44
NFO1620/223/317/1910,5422,58219.17245.33
NFO1720/221/319/1910,5852,58218.55257.22
NFO1820/222/318/1910,7342,58419.17228.65
NFO1920/223/317/1910,5572,58219.4261.31
NFO2021/223/318/1910,1722,57819.94219.08
Source(s): Author’s own work

In the headers of Tables 4 and 5, the numbers to the left of the slash in the fields for train count, high-speed train count, and low-speed train count represent the number of trains selected by the model optimization, while the numbers to the right represent the provided candidate train count. Comparative analysis between case FO (considering flagship trains) and case NFO (without flagship train considerations) reveals that cases with flagship train operating constraints select fewer trains on average. The following sections provide a comparative analysis of the two cases based on train count and capacity utilization rate indicators.

  1. Changes in Train Numbers

As shown in Figure 5, the average number of trains in case FO is lower than that in case NFO. The average number of trains in case NFO is 20, while case FO has an average of 19.1 trains, representing a decrease of 0.9 trains on average. In experiments 1, 3, 5, 6, 10, 12, and 14, the number of trains remains equal, indicating that under identical stopping pattern alternatives, both models select the same number of trains. Research analysis reveals that under a given set of alternative train options, operating flagship trains may reduce the total number of trains. In the 20 experimental groups, under the given stopping pattern alternatives and other case parameters, when two flagship trains with hourly departures are operated within a 2-h period, the number of trains decreases by an average of 0.9.

Figure 5
A line graph shows train number versus cases, comparing two curves for Case F O and Case N F O.The vertical axis is labeled “Train Number,” ranging from 15 to 23 in increments of 2. The horizontal axis is labeled “Cases,” ranging from 1 to 20 in increments of 1. The “Case F O” curve, a solid line with circular markers, starts at (1, 20), shows a steady pattern until (5, 20), falls slightly to (6, 19), rises sharply to a peak at (7, 22), then drops again to (8, 20). It remains nearly stable with minor fluctuations between (11, 20) and (15, 20), maintains a steady trend until (19, 20), and ends with a small rise near (20, 21). The “Case N F O” curve, a solid line in a different color with circular markers, begins at (1, 20), falls gradually to (2, 19), rises again around (3, 20), and then shows several minor rises and falls between (4, 18) and (18, 19) before ending near (20, 19). Note: All numerical data values are approximated.

Train Number Change of the Case FO and NFO. Source: Author’s own work

Figure 5
A line graph shows train number versus cases, comparing two curves for Case F O and Case N F O.The vertical axis is labeled “Train Number,” ranging from 15 to 23 in increments of 2. The horizontal axis is labeled “Cases,” ranging from 1 to 20 in increments of 1. The “Case F O” curve, a solid line with circular markers, starts at (1, 20), shows a steady pattern until (5, 20), falls slightly to (6, 19), rises sharply to a peak at (7, 22), then drops again to (8, 20). It remains nearly stable with minor fluctuations between (11, 20) and (15, 20), maintains a steady trend until (19, 20), and ends with a small rise near (20, 21). The “Case N F O” curve, a solid line in a different color with circular markers, begins at (1, 20), falls gradually to (2, 19), rises again around (3, 20), and then shows several minor rises and falls between (4, 18) and (18, 19) before ending near (20, 19). Note: All numerical data values are approximated.

Train Number Change of the Case FO and NFO. Source: Author’s own work

Close modal
  1. Changes in Capacity Utilization Rate

The capacity utilization rate here is calculated as the current weighted train-kilometer conversion coefficient divided by the total number of alternative stopping trains, which equals the model's optimization objective value divided by the upper limit of alternative stopping trains. As shown in Figure 6, the capacity utilization rate of the optimization model considering flagship trains is lower on average than that of the capacity utilization model without value distinction. The average capacity utilization rate for case NFO is 86.43%, while case FO achieves 83.73%, showing an average decrease of 8.4% in capacity utilization rate when flagship trains are considered. In experimental cases FO3, FO6, and FO14, the capacity utilization rates remain identical, indicating that under these experimental parameters, both models produce consistent train operation optimization schemes. The capacity utilization rate calculated using train-kilometer conversion coefficients can reflect capacity changes in long and short routes that cannot be distinguished by train numbers alone. Research analysis indicates that under a given set of alternative trains, operating flagship trains as a brand value behavior may reduce capacity utilization rate by 8.4%.

Figure 6
A line graph shows train number versus cases, comparing Case F O and Case N F O curves.The vertical axis is labeled “Train Number,” ranging from 78 percent to 92 percent in increments of 2 percent. The horizontal axis is labeled “Cases,” ranging from 1 to 20 in increments of 1. The “Case F O” curve, a solid line with circular markers, starts at (1, 86), falls slightly to (2, 83.941), rises to (3, 87.113), and dips marginally to (6, 85.013), then rises sharply to a peak at (7, 90.356). After a short fall to (8, 85.241), it fluctuates between (9, 87.102) and (15, 88.572), dips again at (17, 84.304), and climbs gradually to end near (20, 91.032). The “Case N F O” curve, also a solid line in a different color with circular markers, begins at (1, 85.012), declines to (2, 81.876), and rises again near (3, 87.274). It falls steadily to (4, 80.745), recovers to (5, 85.112), and fluctuates slightly through (6, 84.926) and (7, 84.702). The curve dips to (8, 83.112), rises again to (9, 86.045), and continues alternating between minor peaks and troughs until it reaches (17, 79.043), the lowest point. It then rises modestly and ends at (20, 83.646). Note: All numerical data values are approximated.

Capacity utilization of the case FO and NFO. Source: Author’s own work

Figure 6
A line graph shows train number versus cases, comparing Case F O and Case N F O curves.The vertical axis is labeled “Train Number,” ranging from 78 percent to 92 percent in increments of 2 percent. The horizontal axis is labeled “Cases,” ranging from 1 to 20 in increments of 1. The “Case F O” curve, a solid line with circular markers, starts at (1, 86), falls slightly to (2, 83.941), rises to (3, 87.113), and dips marginally to (6, 85.013), then rises sharply to a peak at (7, 90.356). After a short fall to (8, 85.241), it fluctuates between (9, 87.102) and (15, 88.572), dips again at (17, 84.304), and climbs gradually to end near (20, 91.032). The “Case N F O” curve, also a solid line in a different color with circular markers, begins at (1, 85.012), declines to (2, 81.876), and rises again near (3, 87.274). It falls steadily to (4, 80.745), recovers to (5, 85.112), and fluctuates slightly through (6, 84.926) and (7, 84.702). The curve dips to (8, 83.112), rises again to (9, 86.045), and continues alternating between minor peaks and troughs until it reaches (17, 79.043), the lowest point. It then rises modestly and ends at (20, 83.646). Note: All numerical data values are approximated.

Capacity utilization of the case FO and NFO. Source: Author’s own work

Close modal

This paper focuses on solving the capacity utilization optimization problem under the condition of flagship train hourly departure scheduling, supporting research on capacity utilization with brand value considerations. The study identifies brand value with the operation of “flagship trains” as high-value popular services, and based on the operational characteristics of flagship trains, treats the scheduling of fast hourly trains as a means of brand value enhancement. A brand value-oriented capacity utilization optimization model for flagship train operation is constructed and compared with high-speed rail capacity utilization models that do not consider flagship trains as the benchmark for value-oriented capacity utilization optimization results. Analysis of cases NFO and FO demonstrates that under given alternative train sets, operating 2 flagship trains within a 2-h timeframe results in an average capacity utilization rate decrease of 8.4%. This finding provides decision support for brand value-oriented high-speed rail capacity utilization research.

The supplementary material for this article can be found online.

D’Ariano
,
A.
,
Pacciarelli
,
D.
, &
Pranzo
,
M.
(
2008
).
Assessment of flexible timetables in real-time traffic management of a railway bottleneck
.
Transportation Research Part C: Emerging Technologies
,
16
(
2
),
232
245
. doi: .
Fan
,
J.
,
Shen
,
H.
,
Liu
,
L.
,
Li
,
B.
, &
Liang
,
Y.
(
2024
).
Evaluation of electric multiple unit train path efficiency based on DEA model
.
Railway Transport and Economy
,
46
(
7
),
23
.
Gao
,
T.
,
Chen
,
J.
,
Xu
,
H.
,
Liu
,
Y.
, &
Chen
,
Z.
(
2025
).
High-speed train rescheduling under disruptions: A computational intelligence framework for passenger transfer optimization
.
Computers & Industrial Engineering
,
209
, 111459. doi: .
Gu
,
Y.
,
Wei
,
Y.
, &
Xia
,
Y.
(
2025
).
Optimization method for maximizing capacity utilization of high-speed railway train timetable based on node routing model
.
Journal of Rail Transport Planning & Management
,
36
, 100548. doi: .
Jensen
,
L. W.
,
Landex
,
A.
,
Nielsen
,
O. A.
,
Kroon
,
L. G.
, &
Schmidt
,
M.
(
2017
).
Strategic assessment of capacity consumption in railway networks: Framework and model
.
Transportation Research Part C: Emerging Technologies
,
74
,
126
149
. doi: .
Jensen
,
L. W.
,
Schmidt
,
M.
, &
Nielsen
,
O. A.
(
2020
).
Determination of infrastructure capacity in railway networks without the need for a fixed timetable
.
Transportation Research Part C: Emerging Technologies
,
119
, 102751. doi: .
Lu
,
Y.
(
2019
).
Study on the capacity of beijing-shanghai high speed railway under overtaking condition
.
Railway Economics Research
,
05
,
34
.
Mei
,
T.
,
Lai
,
J.
,
Zhang
,
Z.
,
Pan
,
J.
, &
Cui
,
Z.
(
2024
).
Research on the benefit and efficiency evaluation of the high-speed railway operation plan oriented to the railway network
. In
International Conference on Smart Vehicular Technology, Transportation, Communication and Applications
(pp. 
193
206
)
Singapore
:
Springer Nature Singapore
.
Pouryousef
,
H.
, &
Lautala
,
P.
(
2015
).
Hybrid simulation approach for improving railway capacity and train schedules
.
Journal of Rail Transport Planning & Management
,
5
(
4
),
211
224
. doi: .
Qu
,
S.
,
Shi
,
J.
, &
Xu
,
R.
(
2022
).
Assessment and analysis of carrying capacity of 350 km/h EMUs operating on Beijing-shanghai HSR
.
China Railway
,
02
,
36
.
Sun
,
G.
,
Jing
,
Y.
,
Li
,
H.
,
Tian
,
Z.
, &
Tian
,
X.
(
2024
).
An interpretable machine learning framework-based approach for predicting passenger flow distribution in train riding sections
.
Journal of Transportation Systems Engineering and Information Technology
,
24
(
02
),
249
.
Wu
,
J.
,
Shan
,
X.
, &
Zhao
,
S.
(
2025
).
Optimization of passenger train line planning adjustments based on minimizing systematic costs
.
Inventions
,
10
(
4
),
64
. doi: .
Xia
,
M.
,
Li
,
B.
,
Zhang
,
B.
, &
Zhang
,
X.
(
2024
).
Research on benchmark train working diagram based on equilibrium and capacity utilization
.
Journal of Railway Science and Engineering
,
21
(
06
),
2162
.
Xu
,
H. Z.
,
Chen
,
J. H.
,
Zhang
,
X. C.
,
Lu
,
T. E.
,
Gao
,
T. Z.
,
Wen
,
K.
, &
Ma
,
Y.
(
2023
).
High-speed train timetable optimization based on space–time network model and quantum simulator
.
Quantum Information Processing
,
22
(
11
),
418
. doi: .
Zhang
,
X.
,
Fan
,
J.
,
Wang
,
B.
, &
Guo
,
Y.
(
2022
).
Collaborative optimization of train stop plan and timetable of Beijing-shanghai high speed railway based on capacity utilization
.
Railway Transport and Economy
,
44
(
10
),
1
.
Zhong
,
Q.
,
Zang
,
L.
,
Yang
,
R.
,
Sheng
,
Z.
, &
Xu
,
R.
(
2024
).
Analyzing the impact of timetable elements on railway line capacity
.
Journal of Rail Transport Planning & Management
,
30
, 100450. doi: .
Published in Railway Sciences. Published by Emerald Publishing Limited. This article is published under the Creative Commons Attribution (CC BY 4.0) licence. Anyone may reproduce, distribute, translate and create derivative works of this article (for both commercial and non-commercial purposes), subject to full attribution to the original publication and authors. The full terms of this licence may be seen at Link to the terms of the CC BY 4.0 licence.

Supplementary data

or Create an Account

Close Modal
Close Modal