Integrated Optimization of the Location–Inventory Problem of Maintenance Component Distribution for High-Speed Railway Operations

: Sustainable distribution network design for the maintenance components of electric multiple units (EMUs) is critical to reduce the problem of unreasonable resource allocation and capital occupation of high-speed railway (HSR) operations. Motivated by the above analysis, this study investigates the integrated optimization of the location and inventory of EMU maintenance component distributions. Aiming to improve the sustainable operation for high-speed railway, we proposed a corresponding nonlinear mixed-integer programming model to determine the location of the distribution center (DC) for EMU maintenance component delivery, inventory control strategy, and corresponding service level. The above optimization model is solved by an adaptive improved genetic algorithm. The proposed model and algorithm are applied to a real-world case study on China’s EMU maintenance components. The ﬁndings show that a higher service level is not better to achieve the lower total cost in the maintenance component distribution network. The ratios of transportation modes are signiﬁcant to balance the service level and total cost of the EMU distribution network. Furthermore, the unit out-of-stock cost and the service level both show great impacts on the total costs of the EMU distribution system. Finally, there exists an optimal ratio of di ﬀ erent transport modes, which ensures the least total cost of the EMU distribution system.


Introduction
As a new method of transportation, high-speed railways have attracted extensive attention in recent years. By the end of 2019, China had built a high-speed rail system of 35,000 km, which is more than 2/3 of the total length of all global high-speed rail systems [1,2]. Along with the construction of high-speed rail networks, the demand for electronic multiple unit (EMU) trains is increasingly strong. For example, as of early 2019, the China Railway Administration owned 2827 EMU trains, an increase of 11.3% year-on-year. The first high-speed line was the Beijing-Tianjin high-speed rail line, which was put into operation in 2008 [3]. According to the maintenance standards, the upcoming maintenance period will be a problem for a considerable number of EMU trains. The long-term plan for high-speed railway networks and the current macroeconomic background raise concerns for scholars and managers regarding the distribution network design of maintenance components for high-speed railway operation.
Inevitably, the construction of the distribution network for EMU maintenance components involves the distribution network design, the distribution transportation mode, and the corresponding inventory design. Moreover, our proposed problem aims to determine the optimal location selection of EMU maintenance component DCs and the combined transport modes, and considers also the stock shortage costs related to the corresponding service levels. Finally, our research findings show the relationship between the optimal service level and ratios of different transportation modes.
Finally, some features of related studies with the location-inventory problem are illustrated in Table 1. Thus, the main contributions of this paper to filling the research gap are as follows: (i) to construct a joint optimization model for the LIP considering multimodal transport and to minimize the cost of fixed operating, inventory, transportation, and stockout penalties. (ii) Because comprehensive location-inventory decisions are more complicated than individual decisions, we propose an improved adaptive genetic algorithm using matrix encoding. In addition, an efficient genetic operator is designed to address the model. (iii) According to the actual background of a high-speed railway, the model provides the foundation and strategy for solving the optimization problem of the EMU maintenance component distribution network.

Problem Description and Modeling Process
This paper focuses on a location-inventory model based on real situations. The China Railway Administration, as a company with a social welfare focus, must consider safe public travel. Thus, the distribution network needs to ensure the timely supply of EMU maintenance components. Otherwise, large-scale high-speed train outages and even social unrest may occur. Based on this characteristic, we added penalty costs, which are out-of-stock costs, into the location inventory formulation that distinguishes maintenance components of the EMU distribution network from regular components.

Problem Description
The considered three-layer distribution network is shown in Figure 1. It consists of multiple suppliers, potential distribution centers, and EMU depots. This problem aims to minimize the total cost of distribution network design. Location is a high initial investment and has long-term impact. Furthermore, inventory policy selection is imperative. Due to stochastic demand, determining inventory policy can solve problems such as how often to reorder, the economic order quantity (EOQ), and the level of safety stock to establish. It is significant to deal with the trade-off between the safety stock and the corresponding service level in supply chain management. To set a lower level of safety stocks will result in decreasing of the service levels and increasing of the corresponding penalty cost. On the contrary, the higher the safety stock, the higher the ratio of meeting the demand with uncertainty and the higher the service level. However, high safety stock will cause high capital occupancy, as well as the lower ratio of turnover of EMU components. Finally, it is important to choose transport modes. As we all know, different modes of transportation have different advantages. Air freight is efficient and fast but is greatly affected by weather. High-speed rail transport makes up for this shortcoming but is still in its infancy. The road transport cost is the smallest compared to that of the previous two modes, but the timeliness is poor. This study never considers the cost in transit, which distinguishes the advantages of different transportation modes by order cost and transportation cost. Therefore, their impact on the service level can be reflected in two aspects: the fixed order cost and the transportation cost. So, it is necessary to choose an optimal combined transport service scheme to reduce the total cost of the maintenance component distribution system.
Sustainability 2020, 12, x FOR PEER REVIEW 5 of 27 choose transport modes. As we all know, different modes of transportation have different advantages. Air freight is efficient and fast but is greatly affected by weather. High-speed rail transport makes up for this shortcoming but is still in its infancy. The road transport cost is the smallest compared to that of the previous two modes, but the timeliness is poor. This study never considers the cost in transit, which distinguishes the advantages of different transportation modes by order cost and transportation cost. Therefore, their impact on the service level can be reflected in two aspects: the fixed order cost and the transportation cost. So, it is necessary to choose an optimal combined transport service scheme to reduce the total cost of the maintenance component distribution system.

Supplier
Distribution Center EMU depot In summary, our model aims to minimize the total cost, which is the sum of the location cost, transportation cost, inventory cost, and penalty cost. The inventory cost in the model includes the fixed ordering, the holding cost of the cycle, and the safety stock in the distribution center. Through this mixed-integer nonlinear programming model of location inventory, we can finally determine (i) the facility location and number of distribution centers; (ii) the allocation plan from the suppliers to distribution centers and distribution centers to EMU depots, and the choice of transportation mode; and (iii) the distribution center economic order quantity, reorder point, and safety stock.

Assumptions and Notations
The study makes the following assumptions: (i) Every distribution center has a fixed construction cost; (ii) Normal distributions of known mean and variance are followed by the maintenance component demand of each EMU depot. There are independent demands for different EMU depots; (iii) There is only one supplier per distribution center and only one distribution center per EMU depot; (iv) Each distribution center adopts a continuous review (r, Q), which means that when the inventory level falls below the reorder point r, each distribution center orders a fixed quantity Q at once; (v) The unit transportation cost is proportional to distance and product shipment.
The model sets, parameters, decision variables, and intermediate variables used throughout the paper are listed as follows: Indices: In summary, our model aims to minimize the total cost, which is the sum of the location cost, transportation cost, inventory cost, and penalty cost. The inventory cost in the model includes the fixed ordering, the holding cost of the cycle, and the safety stock in the distribution center. Through this mixed-integer nonlinear programming model of location inventory, we can finally determine (i) the facility location and number of distribution centers; (ii) the allocation plan from the suppliers to distribution centers and distribution centers to EMU depots, and the choice of transportation mode; and (iii) the distribution center economic order quantity, reorder point, and safety stock.

Assumptions and Notations
The study makes the following assumptions: (i) Every distribution center has a fixed construction cost; (ii) Normal distributions of known mean and variance are followed by the maintenance component demand of each EMU depot. There are independent demands for different EMU depots; (iii) There is only one supplier per distribution center and only one distribution center per EMU depot; (iv) Each distribution center adopts a continuous review (r, Q), which means that when the inventory level falls below the reorder point r, each distribution center orders a fixed quantity Q at once; (v) The unit transportation cost is proportional to distance and product shipment.
The model sets, parameters, decision variables, and intermediate variables used throughout the paper are listed as follows: Indices: I: Index for EMU depots, i = 1, 2, . . . I J: Index for potential distribution center locations, j = 1, 2, . . . J K: Index for suppliers, k = 1, 2, . . . K T: Index for mode of transport, t = 1,2, . . . T Parameters: f j : Fixed operating cost of opening a distribution center j per unit time (¥) C t : Unit transportation cost by transport mode t (¥/unit·km) T k j : The distance from supplier k to distribution center j (km) S ji : The distance from distribution center j to EMU depot i (km) µ i : Mean daily demand at EMU depot i (unit) σ 2 i : Variance of daily demand at EMU depot i (unit 2 ) χ: Days worked in a year (day) LT j : Lead time of the distribution center j to order from supplier k (day) OC t j : Fixed ordering cost at distribution center j by transport mode t (¥/order) HC j : Holding cost at distribution center j (¥/unit) SC j : Shortage cost at distribution center j (¥/unit) α j : Service level for distribution center j Z α : Safety factor related to service level and standard normal distribution, P(Z ≤ Z α ) ≤ α Intermediate Variables: D j : Mean of daily demand allocated to potential distribution center j (unit) V 2 j : Variance of daily demand allocated to potential distribution center j (unit 2 ) Q j : Optimal order quantity for distribution center j (unit) r j : Reorder point for distribution center j (unit) Decision Variables: X j ∈ {0, 1}, Takes the value 1 when located at distribution center j; otherwise, the value 0; Y t k j ∈ {0, 1}, Takes the value 1 when distribution center j is allocated to supplier k by transport mode t; otherwise, the value 0; Z t ji ∈ {0, 1}, Takes the value 1 when EMU depot i is allocated to distribution center j by transport mode t; otherwise, the value 0.

Location-Inventory Model
The location-inventory problem of the maintenance components was formulated as a mixed-integer nonlinear programming model. The objective function of the model was to minimize the total cost, including the fixed operating cost, transportation cost, inventory cost, and penalty cost as follows: (1) The fixed operating cost of opening distribution centers, given as j f j X j ; (2) The transportation cost from the supplier to the opened distribution center, given by k j t C t T k j χµ i Z t ji Y t k j , and from the opened distribution center to the EMU depot, given by We assume that the transportation cost is related to maintenance component shipment, transportation distance, and transport mode; (3) The inventory cost includes the ordering cost, the holding cost of working inventory and safety stock. The mean of daily demand at EMU depot i follows a normal distribution with N µ i , σ i 2 .
Because of unrelated demand, the mean of daily demand at distribution center j also obeys the normal To facilitate the model solution, Shen et al. [34] proposed approximating the demand of a normal distribution by a Poisson distribution. It is common to assume that σ 2 i /µ i = γ, ∀i, where γ is the ratio of variance and mean. Based on this assumption, safety stock can be obtained as ss j = Z α γD j LT j , ∀ j, and the reorder point is r j = D j LT j + Z α γD j LT j , ∀j; Sustainability 2020, 12, 5447 7 of 24 (4) For each distribution center j, we set the penalty cost caused by stock shortages. Miranda and Garrido [31] assumed a normal demand distribution for stochastic demand during lead time, SD j LT j , and the density function of SD j LT j is computed by (1): (1) Expected out-of-stock demand is given by (2).
and +∞ r j The penalty cost for a distribution center j in a single ordering period can be obtained by (5).
In short, the inventory and penalty cost of the distribution center j can be expressed by (6), where χD j /Q j is the expected order frequency in a year.
The first term of (6) is the fixed ordering cost, the second term represents the holding cost with working inventory and safety stock, and the third term expresses the out-of-stock penalty cost.
Consistent with [31,33], this paper assumes that the EOQ model determines the optimal order quantity of each distribution center j, Q j = 2OC t j χD j /HC j . Thus, the inventory and penalty cost for each distribution center j can be rearranged by (7).
Extracting the constant term in (7), the inventory cost and penalty cost can be rewritten by (8) and (9), respectively. where The service level determines the safety stock on the one hand, and is related to the stock shortage on the other hand. Penalty costs are incurred due to out-of-stock inventory, which shows great impact on the total cost of the entire distribution system because the China Railway Administration should Sustainability 2020, 12, 5447 8 of 24 meet the requirement of EMU maintenance component supply in a timely manner, so as to provide a strong guarantee for on-time operation of high-speed train operators.
As described above, this paper analyzes the inventory cost of the distribution center j; thus, we propose a location-inventory model for the optimization of the distribution network of the EMU maintenance components. The model is formulated as follows (10)- (17): The objective function (10) minimizes the total cost of the distribution network for EMU maintenance components and warrants further explanation. The first term is the fixed cost of opening a DC to operate at site j. The second term indicates the inventory cost at the opening DCs, which includes fixed ordering cost, holding cost of the cycle and the safety stock for each DC. The third term is the penalty cost of stockout. The last two terms represent the transportation cost from supplier k to opening DC j and from opening DC j to EMU depot i.
Constraints (11) state that a DC is only serviced by a supplier; constraints (12) show that each EMU depot is only serviced by a DC. Constraints (13) stipulate that the DC-serviced EMU depots are open. Constraints (14) and (15)

Adaptive Genetic Algorithm
Heuristic algorithms are the most appropriate approaches that are used to solve the proposed mixed-integer nonlinear programming (MINLP). There are some classical heuristic methods, such as simulated annealing (SA), tabu search (TS), and genetic algorithm (GA). Each algorithm has its own strengths and weaknesses. GA, based on survival of the fittest and population genetics, is the most widely used in optimal NP-hard problems. Since the LIP is essentially an NP-hard problem, this paper adopted a GA to solve it. Starting from a set of feasible solutions, a GA finds better solutions by simulating the mechanism of selections and genetics in the natural world, such as crossover and mutation processes, and evolves according to the regulation of nature-survival of the fittest. The GA encodes the problem parameters into genes in the chromosome for optimization, and the resolution process starts with a population rather than an individual. Thus, one of the best features of a GA is that with implicit parallel search, the possibility of falling into a local optimal solution is reduced. Therefore, in this paper, we chose the GA to solve this INLP because of its performance in parallel search and avoiding location-optimal solutions. However, undesirable effects and low efficiency are always found when using a GA to solve the LIP because of the initial solution, imperfect coding of chromosomes, and inappropriate crossover and mutation rates. In our adaptive genetic algorithm, we randomly generated the initial solution and adopted multi-matrix coding of chromosomes to improve the speed and enhance the legitimacy of the chromosomes to reduce the probability of illegal solution.

Chromosome Representation
In GAs, the conversion method that converts the feasible solution of a problem from a model variable to the search space that the GA can handle is called coding. The encoding and decoding of chromosomes determine the rate of solution. In this paper, the decision variables of the proposed model consist of three binary variables, and two vectors and four matrices make up the code for each chromosome. To facilitate the description of the encoding and decoding process, Figure 2 is presented with the scale of the three suppliers, five distribution centers, and five EMU depots. The two vectors represent the binary allocation decisions of the supplier-distribution center and distribution center-EMU depot. The four matrices represent the economic order quantity between suppliers and distribution centers (S-DC Quantity), the transport mode between suppliers and distribution centers (S-DC Transport mode), the economic order quantity between distribution center and EMU depot (DC-E Quantity), and the transport mode (DC-E Transport mode). A visual description of these matrices is shown in Figure 2; the former number represents quantity and the latter number represents transport mode. The rows and columns of each matrix are determined by suppliers, distribution centers, or EMU depots. The volume of EMU maintenance components is a real number, and the mode of transportation is considered only as an integer.
Sustainability 2020, 12, x FOR PEER REVIEW 9 of 27 and adopted multi-matrix coding of chromosomes to improve the speed and enhance the legitimacy of the chromosomes to reduce the probability of illegal solution.

Chromosome Representation
In GAs, the conversion method that converts the feasible solution of a problem from a model variable to the search space that the GA can handle is called coding. The encoding and decoding of chromosomes determine the rate of solution. In this paper, the decision variables of the proposed model consist of three binary variables, and two vectors and four matrices make up the code for each chromosome. To facilitate the description of the encoding and decoding process, Figure 2 is presented with the scale of the three suppliers, five distribution centers, and five EMU depots. The two vectors represent the binary allocation decisions of the supplier-distribution center and distribution center-EMU depot. The four matrices represent the economic order quantity between suppliers and distribution centers (S-DC Quantity), the transport mode between suppliers and distribution centers (S-DC Transport mode), the economic order quantity between distribution center and EMU depot (DC-E Quantity), and the transport mode (DC-E Transport mode). A visual description of these matrices is shown in Figure 2; the former number represents quantity and the latter number represents transport mode. The rows and columns of each matrix are determined by suppliers, distribution centers, or EMU depots. The volume of EMU maintenance components is a real number, and the mode of transportation is considered only as an integer.  First, under the condition that an EMU depot is only serviced by one distribution center, the distribution center randomly serves the EMU depots in the DC-E echelons and repeats this allocation process until the demand of each EMU depot is assigned and the relevant transport modes are selected.
Subsequently, in the S-DC echelons, the similar allocation procedure between suppliers and distribution centers is repeated, and the volume of EMU maintenance components is the total demand of each EMU depot serviced by the open distribution center. In this echelon, the supplierserviced open distribution center and the corresponding transport flows and transportation modes are determined.
This form of representation facilitates the manipulation, search, and expression of chromosomes in the evolutionary process.

Genetic Operator
Genetic operators are an important part of GAs, including selection, crossover, and mutation operators used to produce the offspring. To diversify the population, the three major operators were employed to select individuals with good fitness, crossover, or mutate in a current population with a certain probability.

Selection
In this paper, we used a combination of roulette and optimal individual preservation to select excellent individuals from the population and eliminate individuals with low fitness. This combination selection strategy can inherit the contemporary optimal individual into the next generation population. In this way, we could avoid the influence of crossover or mutation on its First, under the condition that an EMU depot is only serviced by one distribution center, the distribution center randomly serves the EMU depots in the DC-E echelons and repeats this allocation process until the demand of each EMU depot is assigned and the relevant transport modes are selected.
Subsequently, in the S-DC echelons, the similar allocation procedure between suppliers and distribution centers is repeated, and the volume of EMU maintenance components is the total demand of each EMU depot serviced by the open distribution center. In this echelon, the supplier-serviced open distribution center and the corresponding transport flows and transportation modes are determined.
This form of representation facilitates the manipulation, search, and expression of chromosomes in the evolutionary process.

Genetic Operator
Genetic operators are an important part of GAs, including selection, crossover, and mutation operators used to produce the offspring. To diversify the population, the three major operators were employed to select individuals with good fitness, crossover, or mutate in a current population with a certain probability.

Selection
In this paper, we used a combination of roulette and optimal individual preservation to select excellent individuals from the population and eliminate individuals with low fitness. This combination selection strategy can inherit the contemporary optimal individual into the next generation population. In this way, we could avoid the influence of crossover or mutation on its chromosome structure, while the remaining individuals who survive as the fittest are determined by roulette to enter the next generation. This method can not only make the individual generation of the largest fitness in the population better than the first generation but also ensure the individual diversity in the population and choose whether to inherit the next generation according to the ratio of the fitness function value.
The roulette operator is also called the proportional selection operator, assuming that u f u is the sum of the fitness function values of a certain population, and f u is the fitness function value of the u individual; then, P u is the probability that the individual is selected and their genes are passed to the next generation for (18): That is, the probability of selection and inheritance in the next generation is directly proportional to the individual fitness value. By simulating the roulette, the number of times each chromosome is present in the population is obtained, and the genetic purpose is achieved. Combined with the best individual preservation strategy, individuals with the largest fitness function value in the population can be directly copied to the next generation, and to ensure their superiority, they do not participate in crossover and mutation operations.

Crossover
In this paper, a partially matched crossover (PMX) was used to randomly select two intersections in a chromosome, thereby obtaining a partial matching region and then exchanging the region between two intersections in the parent to obtain two offspring chromosomes.
Step 1: Determine whether an individual has to crossover; Step 2: Determine the crossover section; Step 3: Determine the crossover position, namely, which columns to exchange; Step 4: Modify the relationship between individual fragments. If there are sections that do not meet the condition, then reconstruct upstream, and the process is similar to the initial solution.
The crossover procedure is shown in detail in Figure 4.

Mutation
In this section, we used the method of interchange mutation to exchange new genes to obtain new chromosomes.
Step 1: Determine whether an individual has to mutate; Step 2: Determine the mutated section; Step 3: Determine the mutated position, namely, determine which lines to exchange; Step 4: Modify the relationship between individual fragments. If there are sections that do not meet the condition, then reconstruct upstream, and the process is similar to the initial solution.
The mutation procedure is shown in detail in Figure 3.

Adaptive Probabilities of Crossover and Mutation
In the genetic algorithm, the choice of the crossover probability c P and the mutation probability m P will affect the performance and effectiveness of the algorithm. Unreasonable parameter selection may result in the algorithm not converging to the global optimum and falling into the local optimum. The traditional genetic algorithm presets the crossover and mutation probability to a specific value, which is calculated by the probability in the whole genetic process. When the iteration reaches a certain level, it is obviously unreasonable to cross or mutate with a certain probability. To affect the performance of the algorithm, Srinivas and Patnaik [50] proposed an adaptive genetic algorithm through a large number of experimental studies. In recent years, many scholars have improved the genetic probability and proposed different adaptive genetic algorithms. In this paper, the improved formulas for calculating the adaptive crossover ratio and mutation ratio are shown in (19) and (20):

Process of the Genetic Algorithm
Based on the above description, we can conclude the main process of the adaptive genetic algorithm:   Figure 3. Crossover process.

Mutation
In this section, we used the method of interchange mutation to exchange new genes to obtain new chromosomes.
Step 1: Determine whether an individual has to mutate; Step 2: Determine the mutated section; Step 3: Determine the mutated position, namely, determine which lines to exchange; Step 4: Modify the relationship between individual fragments. If there are sections that do not meet the condition, then reconstruct upstream, and the process is similar to the initial solution.
The mutation procedure is shown in detail in Figure 4.

Adaptive Probabilities of Crossover and Mutation
In the genetic algorithm, the choice of the crossover probability P c and the mutation probability P m will affect the performance and effectiveness of the algorithm. Unreasonable parameter selection may result in the algorithm not converging to the global optimum and falling into the local optimum.
The traditional genetic algorithm presets the crossover and mutation probability to a specific value, which is calculated by the probability in the whole genetic process. When the iteration reaches a certain level, it is obviously unreasonable to cross or mutate with a certain probability. To affect the performance of the algorithm, Srinivas and Patnaik [50] proposed an adaptive genetic algorithm through a large number of experimental studies. In recent years, many scholars have improved the genetic probability and proposed different adaptive genetic algorithms. In this paper, the improved formulas for calculating the adaptive crossover ratio and mutation ratio are shown in (19) and (20): where k 1 , k 2 , k 3 , k 4 , k 5 and k 6 all lie in the interval (0, 1) and k 1 > k 2 > k 3 , k 4 > k 5 > k 6 ; F represents the individual fitness value; F min , F max and F avg represent the minimum, maximum, and average values of all the individual fitness functions in the contemporary population, respectively; and F represents the larger individual fitness function value at the intersection of individual operations.

Process of the Genetic Algorithm
Based on the above description, we can conclude the main process of the adaptive genetic algorithm: Step 1: Initialize parameters. Set the population size, encoding and decoding methods, and the maximum number of iterations of the genetic algorithm.
Step 2: Generate the initial population. Two vectors and four matrices are used to randomly generate the initial population and obtain the chromosomes in line with the population size.
Step 3: Obtain the fitness function value. The fitness function is used to calculate the fitness value fu, u = 1,2, . . . popsize.
Step 4: Selection operation. The combination of elite strategy and roulette strategy is adopted to determine whether the individual performs the selection operator. If the elite strategy is satisfied, the individual will be directly inherited to the next generation without crossover or mutation; otherwise, roulette will decide whether to enter the next step or not.
Step 5: Crossover operation. Use the improved adaptive crossover probability formula to obtain the current crossover probability P c and compare it with the number generated randomly between [0, 1] to decide whether to perform the partial matching crossover operation.
Step 6: Mutation operation. Adopt the improved adaptive variation probability formula to obtain the current variation probability P m and compare it with the number generated randomly between [0, 1] to decide whether to perform the swap mutation operation.
Step 7: Termination. Determine whether the maximum number of population evolutions is satisfied. If so, the final solution is obtained; if not, return to Step 4.
The algorithm flow chart is shown in Figure 5.

Computational Implementations
In this section, we considered 52 EMU depots of the China Railway Administration as customers. Meanwhile, because the EMU maintenance component distribution center is usually built close to the EMU depot, 52 EMU depots were also regarded as distribution centers. Four suppliers provided the EMU maintenance components. The above optimal model of n in Section 3.3 and the corresponding algorithm mentioned in Section 4 were applied to solve this two-echelon EMU components distribution network. The corresponding parameters are shown in Table 2.

Computational Implementations
In this section, we considered 52 EMU depots of the China Railway Administration as customers. Meanwhile, because the EMU maintenance component distribution center is usually built close to the EMU depot, 52 EMU depots were also regarded as distribution centers. Four suppliers provided the EMU maintenance components. The above optimal model of n in Section 3.3 and the corresponding algorithm mentioned in Section 4 were applied to solve this two-echelon EMU components distribution network. The corresponding parameters are shown in Table 2. To verify the effectiveness of the algorithm, we set six experiments with different sizes by standard genetic algorithms (SGAs) to compare with proposed adaptive genetic algorithms (AGAs). The basic parameter settings of the algorithm are shown in Table 3. The computational tests were conducted on the MATLAB R2012a platform, the Intel (R) Core (TM) i5-8250U CPU @ 1.60 GHz CPU, 8.00 GB RAM, 64-bit operating system laptop.
Experiments were named as follows: S-DC-E, S, DC and E, representing the number of suppliers, distribution centers, and EMU depots, respectively. Because the EMU depots discussed in this paper were both distribution centers and themselves, the values of DC and E were the same. Comparing the SGA with the AGA, the following results were obtained. Table 4 and Figure 6 show that when the scale was small, the results of the SGA and the AGA were relatively similar, but the larger the scale, the better the AGA was over the SGA. On the scale of 4-60-60, the difference between the AGA and the SGA was 969.32 × ¥10,000. Although the running time of the AGA was longer than that of the SGA, the difference was acceptable within 30 s.   The geographical locations of four suppliers and 52 EMU depots are shown in Figure 7. The figure also shows the allocation and transportation modes from supplier to distribution center and distribution center to EMU depots. Assume that transportation modes 1 to 3 are air, high-speed rail, and road transportation, respectively.
In the actual study, there were four suppliers and 52 EMU depots.  The geographical locations of four suppliers and 52 EMU depots are shown in Figure 7. The figure also shows the allocation and transportation modes from supplier to distribution center and distribution center to EMU depots. Assume that transportation modes 1 to 3 are air, high-speed rail, and road transportation, respectively.    In the actual study, there were four suppliers and 52 EMU depots. Four suppliers: Tianyi Shangjia High-tech Materials Co., Ltd.; Knorr Vehicle Equipment Co., Ltd.; CRRC Qishuyan Locomotive Co., Ltd.; Changsha Boke Metal Products Co., Ltd. were respectively represented by suppliers 1-4. Fifty-two EMU depots: Hongqiao, Nanjing South, Nanxiang, Shanghai South, Hangzhou East, Hefei South, Nanjing, Guangzhou South, Shenzhen North, Changsha South, Huizhou, Guangzhou East, Foshan Xike, Foshan West, Changzhutan, Guangzhu, Sanya, Beijing South, Tianjin, Shijiazhuang, Beijing West, Beijing, Xiamen North, Fuzhou South, Nanchang West, Fuzhou, Nanchang, Longyan, Shenyang North, Dalian North, Shenyang South, Changchun West, Chengdu East, Guiyang North, Chongqing North, Wuhan, Hankou, Jinan, Qingdao North, Qingdao, Nanning, Guilin, Harbin West, Xi'an North, Zhengzhou East, Zhengzhou, Taiyuan South, Kunming South, Urumqi, Lanzhou West, Xining, Hohhot were respectively represented by DC 1-52 and Depot 1-52. See Tables 5 and 6 for detailed allocation. Table 5. Situations of suppliers serving distribution centers.

Suppliers Served Distribution Centers and Transportation Mode
Supplier 1 DC19 by air, DC20 by high-speed rail, DC21 by high-speed rail, DC44 by high-speed rail Supplier 2 DC1 by air, DC2 by highway, DC5 by highway, DC35 by highway, DC51 by high-speed rail Supplier 3 DC7 by air, DC13 by high-speed rail, DC40 by high-speed rail, DC43 by highway Supplier 4 DC10 by high-speed rail, DC15 by high-speed rail, DC18 by highway, DC25 by highway, DC26 by highway, DC39 by high-speed rail, DC41 by highway, DC47 by high-speed rail According to the results of Tables 5 and 6, in the first stage from a supplier to a distribution center, because shipment sizes were large and the distances were long, the road transportation mode was selected 10 times, followed by the high-speed rail transportation mode eight times, and the air transportation mode was selected the least, three times. In the second stage, from a distribution center to an EMU depot, there were fewer shipment sizes. Even if the distance was long, the road was selected more often than the high-speed railway.
From the perspective of the total cost, Figure 8 shows that the proportion of the fixed operating cost (FC) account was the largest, which was 53.44%, followed by the penalty cost (PC), which was 22.37%, and the transportation cost (TC) and inventory cost (IC), which were 21.11% and 3.08%, respectively. It can be seen that with the exception of the fixed operating costs, the penalty costs caused by stock shortages were the largest. Next, we will discuss the parameters related to the stock shortage cost.
From the perspective of the total cost, Figure 8 shows that the proportion of the fixed operating cost (FC) account was the largest, which was 53.44%, followed by the penalty cost (PC), which was 22.37%, and the transportation cost (TC) and inventory cost (IC), which were 21.11% and 3.08%, respectively. It can be seen that with the exception of the fixed operating costs, the penalty costs caused by stock shortages were the largest. Next, we will discuss the parameters related to the stock shortage cost.

Impact of Service Level for Total Cost and Transportation Mode
To explore the impact of penalty costs on total costs and modes of transportation, this section used service levels as variables to derive management decisions about the impact of changing service levels on decisions. Then, we focused on how transportation mode decisions were subject to different service levels for each DC and EMU deport. Figure 9 shows the optimal solution for the model at different service levels. As the service level increased, the total cost first decreased and then increased. The service level at the lowest total cost of the model was 95%.

Impact of Service Level for Total Cost and Transportation Mode
To explore the impact of penalty costs on total costs and modes of transportation, this section used service levels as variables to derive management decisions about the impact of changing service levels on decisions. Then, we focused on how transportation mode decisions were subject to different service levels for each DC and EMU deport. Figure 9 shows the optimal solution for the model at different service levels. As the service level increased, the total cost first decreased and then increased. The service level at the lowest total cost of the model was 95%. Combining Table 7 with Figure 10, in the first stage from a supplier to a distribution center, the proportions of choosing air and high-speed rail transportation modes increased slightly, while the proportion of road transportation gradually decreased. Table 7 and Figure 11 show that in the second stage from a distribution center to an EMU depot, because of the medium distances, the high-speed rail transportation mode was faster than the road, the cost was lower than the air mode, the advantages appeared, the proportions of the air and road transportation modes gradually decreased, and the proportion of the high-speed rail transportation mode gradually increased. However, the inexpensive road transportation mode still accounted for approximately half.  Combining Table 7 with Figure 10, in the first stage from a supplier to a distribution center, the proportions of choosing air and high-speed rail transportation modes increased slightly, while the proportion of road transportation gradually decreased. Table 7 and Figure 11 show that in the second stage from a distribution center to an EMU depot, because of the medium distances, the high-speed rail transportation mode was faster than the road, the cost was lower than the air mode, the advantages appeared, the proportions of the air and road transportation modes gradually decreased, and the proportion of the high-speed rail transportation mode gradually increased. However, the inexpensive road transportation mode still accounted for approximately half.  The findings show that (i) adjusting the service level can achieve the purpose of reducing the total cost of the distribution network, but the concept of safety first does not apply to the inventory, and the higher service level is not better. Although the higher the service level is, the smaller the stockout cost, the transportation cost is higher. (ii) In addition to changing the shipment size and allocation situation, it is feasible for the China Railway Administration to appropriately adjust the transportation mode in the maintenance component distribution network. When choosing the ratio of air, high-speed rail, and road transportation as approximately 2:3:5, the total cost is the smallest.

Impact of Unit Stockout Cost on Total Cost and Transportation Mode
To explore the impact of unit stockout cost on the total cost and the choice of transportation method, we performed a simulation analysis to study the sensitivity of the cost to changes in different service levels and unfulfilled demand costs. Assume that the stockout costs are 0.4, 0.7, and 1.0 × ¥10,000, respectively. The results are shown in Figure 12. The findings show that (i) adjusting the service level can achieve the purpose of reducing the total cost of the distribution network, but the concept of safety first does not apply to the inventory, and the higher service level is not better. Although the higher the service level is, the smaller the stockout cost, the transportation cost is higher. (ii) In addition to changing the shipment size and allocation situation, it is feasible for the China Railway Administration to appropriately adjust the transportation mode in the maintenance component distribution network. When choosing the ratio of air, high-speed rail, and road transportation as approximately 2:3:5, the total cost is the smallest.

Impact of Unit Stockout Cost on Total Cost and Transportation Mode
To explore the impact of unit stockout cost on the total cost and the choice of transportation method, we performed a simulation analysis to study the sensitivity of the cost to changes in different service levels and unfulfilled demand costs. Assume that the stockout costs are 0.4, 0.7, and 1.0 × ¥10,000, respectively. The results are shown in Figure 12.  As we expected, Figure 12 shows that the optimal service level increased with increasing unit stockout cost. When the stockout cost was 0.4, the optimal service level was 0.93. When the unit stockout costs were 0.7 and 1.0, the best service level was 0.95.
Except for changes in total cost, we analyzed the proportions of different transportation modes at the optimal service level for different stockout costs. Table 8 and Figure 13 show that with the increase in the unit stockout cost, in the first stage, the proportion of choosing aircraft transportation increased greatly. Although the proportions of both high-speed rail and road transportation modes declined, the decreased proportion of the latter was greater than that of the former. This implies that the operator of China High-Speed Railway would rather choose air transportation with expensive transport costs and ordering costs than high-speed rail or road transportation if there existed a higher cost of out-of-stock maintenance components because the penalty cost on the out-of-stock cost was much more than the total of the above two costs (i.e., transport cost and ordering cost). Table 8 and Figure 14 show that, in the second stage, the proportions of transportation modes did not change significantly. Consistent with the conclusions reached in Section 5.2.1, the air, highspeed rail, and road transportation methods remained at a ratio of 2:3:5, respectively, with minimal total cost. However, when the maintenance components reached the DC, the distances from DC to the EMU depot served were not far and the demand batch was not large, the impact of the unit shortage cost was less than that of the transportation cost and the order cost. Therefore, when the As we expected, Figure 12 shows that the optimal service level increased with increasing unit stockout cost. When the stockout cost was 0.4, the optimal service level was 0.93. When the unit stockout costs were 0.7 and 1.0, the best service level was 0.95.
Except for changes in total cost, we analyzed the proportions of different transportation modes at the optimal service level for different stockout costs. Table 8 and Figure 13 show that with the increase in the unit stockout cost, in the first stage, the proportion of choosing aircraft transportation increased greatly. Although the proportions of both high-speed rail and road transportation modes declined, the decreased proportion of the latter was greater than that of the former. This implies that the operator of China High-Speed Railway would rather choose air transportation with expensive transport costs and ordering costs than high-speed rail or road transportation if there existed a higher cost of out-of-stock maintenance components because the penalty cost on the out-of-stock cost was much more than the total of the above two costs (i.e., transport cost and ordering cost). Table 8 and Figure 14 show that, in the second stage, the proportions of transportation modes did not change significantly. Consistent with the conclusions reached in Section 5.2.1, the air, high-speed rail, and road transportation methods remained at a ratio of 2:3:5, respectively, with minimal total cost. However, when the maintenance components reached the DC, the distances from DC to the EMU depot served were not far and the demand batch was not large, the impact of the unit shortage cost was less than that of the transportation cost and the order cost. Therefore, when the unit shortage cost changed, the proportions of the different transportation modes had fewer changes.  The findings show that (i) there is a correlation between unit stockout cost and optimal service level. The higher the stockout cost is, the higher the optimal service level, the greater the proportion of air transportation options with high transportation costs and high order costs. (ii) Even at the optimal service level, the higher the unit stockout cost is, the larger the proportion of choosing expensive and fast transportation methods. The China Railway Administration can improve the  The findings show that (i) there is a correlation between unit stockout cost and optimal service level. The higher the stockout cost is, the higher the optimal service level, the greater the proportion of air transportation options with high transportation costs and high order costs. (ii) Even at the optimal service level, the higher the unit stockout cost is, the larger the proportion of choosing expensive and fast transportation methods. The China Railway Administration can improve the supply chain. If a low stockout cost is set, the total cost can be reduced by reducing the optimal service level.

Conclusions and Future Work
This paper addressed an integrated optimization model of the location-inventory problem for EMU component distribution that incorporates the location cost of a DC, the inventory cost, the linear transportation cost from the distribution center to the EMU depot, and the penalty cost. The proposed problem needed to solve the following three key questions: (i) the number and locations of distribution centers; (ii) the allocations from suppliers to the distribution centers and from the distribution centers to the EMU depots, and the choice of transportation modes; and (iii) the distribution center optimal order quantity, reorder point, and safety stock.
To solve the above optimal model, an adaptive genetic algorithm was given. The performance of the proposed adaptive genetic algorithm was evaluated by comparing it with a standard genetic algorithm. Moreover, a real-world case study based on the EMU maintenance component distribution network design of the China High-Speed Railway was employed to verify the effectiveness of the above optimization model and its algorithm.
The main management insights of this study are shown as follows: (i) The higher service level is not necessarily the better, although the higher the service level is, the smaller the stockout cost, but the transportation cost is higher. When higher service levels are required, faster transportation methods are needed, so the cost of the entire supply chain will rise. (ii) In addition to the optimal service level, the China Railway Administration must determine the unit stockout cost; the higher the unit stockout cost is, the larger the proportion of choosing expensive and fast transportation methods. (iii) It is reasonable to choose the mode of transportation and set the proportions of transportation modes. Air, high-speed rail, and road transportation modes are optimal at a ratio of 2:3:5.
Although we used our best efforts to reflect the actual situation of the EMU maintenance component distribution network design, in order to facilitate the theoretical solution we simplified it appropriately.
There are many extensions that should be studied in the future. First, we hope to develop the cases with multiple items. The China Railway high-speed train has nine types of body models composed of more than 140 systems and approximately 40,000 parts. Second, we would like to consider the problem of location inventory transportation and make strategic, tactical, and operational decisions simultaneously. This needs to take into account the schedule of various modes of transportation or vehicle routing problems. Moreover, we need to combine the actual situation, considering the sharing of maintenance components between railway bureaus. Finally we would like to improve the performance of the AGA by embedding with other heuristic algorithms (e.g., tabu search, local search, simulated annealing) and to investigate more comparison or benchmarking instances in the future studies.