Integrated Optimization of Stop Location and Route Design for Community Shuttle Service

: The community shuttle system plays an important role in serving communities with a heavy travel demand for the metro service. Stop location and route design are the two main decisions of planning a community shuttle service. Those two decisions are interrelated and interact, and are strongly related to the user cost and operating cost. The optimal stop location and route can help to reduce the walking distance of passengers and the route length. To make a trade-off between the walking distance of passengers and route length, we propose a discrete optimization problem. A single integrated formulation is established to optimize stop location and route design. Planners can decide the stop location and route design of the community shuttle system simultaneously based on this formulation. Then, we present a non-dominated sorting genetic (NSGA-II) based algorithm to obtain the non-dominated solutions of the discrete optimization formulation. The numerical experiments and a case study based on real-world data are used to demonstrate that the proposed solution method can yield a set of plans of stop location and route in a reasonable time. We also ﬁnd that when the maximum tolerable walking distance is set to 418 m, the trade-off between the total walking distance of passengers and route length can be obtained.


Introduction
A heavy demand for the metro service has arisen along with the continuous expansion of communities in China. In the urban public transport system, the community shuttle service plays an important role in transporting these ever-increasing passengers from their origins to metro stations or from metro stations to their destinations. Obviously, systematic and proper planning of the community shuttle system is necessary before this public transport system is in operation. The planning process of the public transport service includes four components: route and network design, timetabling, vehicle scheduling, and crew scheduling [1]. Among these components, route and network design is a basic one which can influence the overall operational efficiency of the public transport system. There have been many studies concentrating on this component [2][3][4][5][6]. However, most of them are carried out for conventional bus or feeder bus network design, and few scholars focus on community shuttle route optimization.
A distinctive feature of the community shuttle system is that it mainly provides feeder service for the travel demand for metro station access. The travel demand presents a many-to-one or one-to-many travel characteristic. This kind of community shuttle service has been commonly implemented in practice. For example, Figure 1 shows typical community shuttle route access to a metro station (Zhuxinzhuang subway station) in Beijing, China. Community shuttle routes have the following   features: (1) there are a number of demand points (e.g., residential areas, schools, hospitals, commercial areas) to be served by community shuttle routes; (2) there are many candidate stops from which some stops can be selected; (3) among the candidate stops, there is a transfer stop located close to a metro station and passengers can transfer between the metro and community shuttle at this stop; and (4) in general, only one community shuttle route is available for metro station access in each community. These features are visually presented in Figure 2. Here, the stops and route of a community shuttle line are represented by circles and arrowed solid lines, respectively. Usually, a community shuttle route is circular. A transfer stop, denoted as a square, is usually located near a metro station. Hexagons represent demand points which are served by stops selected from candidate stops. Dash lines show the situation of matching demand points to stops.   From Figure 2, it is clear that stop location and route design are the two main decisions in the first component of the community shuttle service planning process. Route design determines route length, which influences fleet size and vehicle travel cost. Moreover, route length is directly related to the in-vehicle time of passengers. Operators and passengers usually hope to reduce route length. However, the community shuttle system aims to improve the accessibility of the public transport system. Therefore, the walking distance of passengers, related to stop location, should be taken into account when the community shuttle system is designed, more than just the route length. The walking distance of passengers is expected to be as short as possible and should be within an acceptable range. Generally, reducing the walking distance of passengers often means more detours or stops, whereas reducing the route length generally increases the walking distance of passengers, as stops may not be allowed to be located near passengers' origins or destinations. Hence, a balance or a trade-off between the two objectives is critical. To analyze this trade-off, a properly integrated method for stop location and route design is necessary for planners.
Considering the features of the community shuttle system, integrated optimization of stop location and route design for the community shuttle service is a multi-objective location-routing problem, which is different from the conventional bus network design problem. When planning conventional bus network, frequent transfers and complicated route plans require planners to consider passengers' route choice behaviors [7].
Our goal is to define an approach to tackle the integration of community shuttle stop location and route design. A multi-objective integer programing formulation is proposed to minimize the total walking distance of passengers, as well as the route length. This formulation takes into account the route length constraints, the stop spacing constraints, and the maximum tolerable walking distance constraint. An algorithm based on the non-dominated sorting genetic algorithm II (NSGA-II) is designed to solve the formulation. The correctness and validity of the formulation and algorithm are verified by numerical examples and a real-world problem in Shanghai, China.
The remainder of this paper is organized as follows. Section 2 reviews the literature on route design, stop location, and the integrated problem of route design and stop location. In Section 3, the problem studied is described and a mathematical formulation for the integrated optimization of stop location and route design for the community shuttle service is proposed. Section 4 introduces the algorithm to solve our formulation. Section 5 presents numerical experiments. A case study is given in Section 6. Conclusions are presented in Section 7.

Literature Review
This section presents a review on three major problems of community shuttle service research: route design, stop location, and the integrated optimization of route design and stop location.

Route Design
Route design is one of the key steps of feeder bus network planning [8] and is a well-studied subject. Many researchers consider route length or travel time associated with the feeder bus route in their works. Shrivastava and Dhingra [9], Verma and Dhingra [10], Shrivastava and O'Mahony [11,12], and Song and Liu [13] generate feeder bus routes using a heuristic algorithm and the route length limit is considered. Their approaches are based on an assumption, i.e., stop location is an input and bus stops are demand points in their works, which is reasonable if planners hope to design feeder bus routes based on existing bus stops or adjust existing feeder bus routes. When planners hope to make new feeder bus route plans, their approaches are not applicable [14]. A similar assumption is used by Szeto and Wu [15], and the travel time of buses is a part of the objective of a mathematical formulation in their work. Teng et al. [16] take into account the travel time constraint, but ignore stop locations based on the assumption that if residents live in the traffic zones passed by the feeder bus, they can be served. This assumption indicates that if a traffic zone is not passed by feeder bus routes, passengers there cannot be served, even if they are near by the bus routes. This may not be reasonable.
The access cost of passengers is not considered by all the above studies. Transit operators and passengers prefer short community shuttle or feeder bus routes to reduce travel time and in-vehicle time. Passengers also hope to easily access the community shuttle or feeder bus service. In addition, operators hope to decrease the access cost of passengers to attract more passengers. Therefore, it is necessary to consider the access cost of passengers. Access cost is considered by some scholars. Chien and Yang [17] and Jerby and Ceder [18] assumed that the travel demand is uniformly distributed along route links or within traffic zones. Hence, travel demand is related to bus routes and locations of bus stops are neglected. This assumption may be easily violated, because buses cannot enter gated communities in China. Instead, passengers will gather at some pick-up/drop-off points [19]. Xiong et al. [20] and Zhu et al. [21] regard access distance as the distance between the center point of each traffic zone and the corresponding road link. However, in fact, buses cannot stop anywhere on bus routes. In addition, Lin and Wong [22] present a multi-objective programming formulation and the maximum service distance of a feeder bus route is considered. They assume that stop spacing is commonly within a reasonable walking distance, thus stop location is not considered. Nevertheless, this assumption is not suitable for communities in the suburbs or large communities.

Stop Location
Few studies have concentrated on the stop location of a community shuttle or feeder bus service. Liu et al. [23] propose a stop planning method for an airport shuttle bus. Their objective is to improve the accessibility of airport shuttle bus stops. However, stop location can also influence travel time due to the extra detours, more than just the accessibility of feeder bus service. Realizing the above critical problem, Zhao and Chien [24] aim at minimizing user cost (access time, wait time, and in-vehicle travel time) and operator cost. However, they consider that travel demand is uniformly distributed along the bus route, which may be unreasonable in China for the same reason described in Section 2.1.
Some studies have been conducted on the stop location of a conventional bus system. Alonso et al. [25] and Moura et al. [26] use a bi-level optimization model to obtain optimal bus stop locations on a macroscopic scale. Their model considers user cost and operator cost; however, the maximum walking distance of passengers is ignored. We believe that it is significant for a public transport service, especially for a community shuttle service, to ensure that no passenger is connected to a bus stop further away than the maximum walking distance. Chen et al. [27] select the optimal stops from candidate stops along bus routes and the bus routes are given as inputs. This may not be valid if planners hope to make bus stop plans from scratch or radically different plans from the existing one. For stop location, the most similar to our model is Jahani et al. [28], who propose a multi-objective bus stop location model to determine the optimal stop locations and the match situation of demand points to stops. Their numerical cases show that the model can generate sensible stop results.

Integrated Optimization of Stop Location and Route Design
Some scholars concentrate on the optimization of stop location and route design in a successive way. Xiong et al. [29] develop a solution method for the optimization of a community shuttle route. A heuristic algorithm is proposed for locating stops on a community shuttle route. This community shuttle route is an input of the heuristic algorithm. Li and Chen [30] and Leksakul et al. [31] determine stop locations by clustering algorithms first, and then bus route are constructed based on bus stops. A successive approach optimizes stop location or route design first, while the results are the input of the other problem. According to its definition, a successive approach may result in sub-optimal solutions.
Other scholars solve the two problems in a completely integrated way. They propose a formulation or a solution methodology to determine the stop location and route design simultaneously. For example, Perugia et al. [14] propose a cluster routing approach to model both bus stop location and route design for a home-to-work bus service in a metropolitan area and the model is solved by a tabu search algorithm. The main difference between the method of Perugia and ours is that we simultaneously determine the number of bus stops to be served, which is considered as an input in their research.
When decision makers hope to make a new plan from scratch or a radically different plan with respect to the existing one, considering the number of bus stops as an input may not be valid. Chen et al. [7] consider how to design a suburban bus route for airport access. Their approach can generate the bus route and stop location simultaneously. However, the number of stops to be served is also given as an input. In order to consider the access cost of passengers, we take into account the locations and travel demand of demand points, which is ignored by Chen et al. [7].
In summary, the access cost of passengers should be considered when optimizing the community shuttle service. In addition, to our knowledge, most of the existing community shuttle or feeder bus service studies ignore the relationship between stop location and route design, and few studies are known about the simultaneous optimization of stop location and route design. However, as a matter of fact, stop location is directly related to route design, thus the integrated optimization of the two problems will lead to a better understanding of the planning of the community shuttle service.
In this paper, we propose a multi-objective integrated formulation to simultaneously decide the stop location and route design of the community shuttle, which is characterized by two aspects that make our problem different from other similar problems in the literature. First, considering that gated communities are common in China, we believe that travel demand is related to demand points. Our approach can decide the community shuttle route, stop location, and the match situation of demand points to stops simultaneously. Second, we can simultaneously decide the number of stops to be served and stop locations. To the best of our knowledge, these features cannot be found in other works (e.g., Perugia et al. [14], Chen et al. [7]) in which the number of stops to be served is an input and the situation of matching demand points to stops is ignored.

Problem Description and Assumptions
The problem considered in this paper can be described as follows. A metro station is located near a community and connects this community with other areas in the same city. There are a set of demand points in the community. The travel demand for the community shuttle service is an input. There is also a set of candidate stops to be considered. Some candidate stops will be chosen as stops and passengers can only be picked up and dropped off at stops which are located within the maximum tolerable walking distance of their demand points. The problem is deciding where stops are located and which demand points can be served by each stop. Furthermore, the shortest loop route based on these stops is determined simultaneously. The objectives are to minimize the total walking distance of passengers and the community shuttle route length.
In order to simplify the problem and make the formulation reflect the practical situation, the following assumptions are made.

1.
The locations of all demand points and candidate stops are known.

2.
The travel demand for the community shuttle service of each demand point is given. 3.
The distances between different candidate stops and the distances between demand points and candidate stops are known.

Mathematical Formulation
To describe the formulation, the following notations are introduced in Table 1.
The integer programing formulation is given below. The objectives (1) and (2) minimize the total walking distance of passengers and the route length, respectively. Constraint (3) defines the degree of the transfer stop vertex. Constraint (4) ensures flow conservation at the candidate stop vertices. Constraints (5) and (6) mean that a candidate stop can serve demand points under the condition that the candidate stop is visited by the community shuttle route. Constraint (7), in conjunction with constraint (8), indicates that if a candidate stop is selected as a stop, the candidate stop must serve at least one demand point. Constraint (9) ensures that each demand point is served by a stop. Constraint (10) ensures that the route length is within a range. Constraint (11) ensures that stop spacing lies within the maximum and minimum values. The maximum tolerable walking distance of passengers is enforced by constraint (12). Constraint (13) indicates that there is only one loop route. Constraints (14) and (15) restrict the feasible domains of decision variables.

Solution Method
The formulation presented in Section 3 is a bi-objective programing formulation which aims to minimize the total walking distance of passengers depending on stop location and, at the same time, to minimize the route length based on route design. Assigning weight to each objective function is a classical solution method of the multi-objective optimization problem. Multiple objectives are converted into a single one with this method. However, we believe that it is of great practical significance to provide decision makers with trade-off curves between values of different objective functions (i.e., the Pareto optimal front). Furthermore, it is difficult for decision makers to determine the weight of each objective reasonably. Therefore, we plan to search for the Pareto optimal front and each Pareto optimal solution is a plan of community shuttle stop location and route design. To this end, an algorithm based on the non-dominated sorting genetic algorithm II (NSGA-II) is designed to find the Pareto optimal front. As one of the most effective multi-objective evolutionary algorithms (MOEAs), NSGA-II has been used in urban public transport planning and operation problems successfully, including the public transport network design problem [32], timetabling problem [33], and vehicle scheduling problem [34]. Deb et al. [35] introduce details about NSGA-II and indicate that solutions obtained by NSGA-II are closer to the true Pareto optimal front compared to two other MOEAs, i.e., PAES and SPEA.
A brief description of NSGA-II is as follows. Firstly, an initial population P 0 of size Z is created. Then, all solutions in the initial population P 0 are sorted into different non-dominated levels by the fast non-dominated sorting procedure. Each non-dominated front is assigned a rank which is equal to its non-dominated level and represents the fitness of solutions in this non-dominated front. Thereafter, the selection operator, crossover operator, and mutation operator are applied to obtain a child population Q 0 of the same size Z. After the child population Q 0 is generated, P 0 and Q 0 are associated to create a new population R 0 of size 2Z. The best Z solutions are chosen from R 0 to fill a new population P 1 . Then, the procedure continues, iterating until the number of iterations I is reached. Finally, the Pareto optimal front is obtained.
Based on the general framework of NSGA-II, we modify chromosome coding and an infeasible solution processing method to better suit our problem. The algorithm based on NSGA-II in this paper is described in detail below.

Chromosome Coding
An integer coding method is used for the chromosome. The length of the chromosome is |D|, which is the number of demand points. Each gene presents which candidate stop serves the corresponding demand point. For example, a gene "5" means that the corresponding demand point is served by the fifth candidate stop. Logically, this candidate stop is chosen as a stop.

Solution Evaluation
The fitness of each solution is calculated based on its values of the two objective functions (i.e., the total walking distance of passengers and the route length). Each solution shows a plan of stop location and the situation of matching demand points to stops. Based on the stop location, the length of the optimal community shuttle route can be solved by formulation M1, shown as follows.
(M1) min C 2 subject to (3) and (4), (10) and (11), (13)- (15) and y kl ∈ {0, 1} ∀k, l ∈ H M1 is a linear programing formulation and can be solved by ILOG CPLEX. The total walking distance of passengers can be solved based on the distance between demand points and candidate stops and the situation of matching demand points to stops presented by chromosome coding.
Then, the fitness of a solution is represented by the non-domination rank and crowding distance. If a solution's non-domination rank is less than other solutions, or if its crowding distance is larger than other solutions with the same non-dominated rank, this solution is fitter than other solutions. For the calculating method of crowding distance, readers can refer to Deb et al. [23].

Infeasible Solution Processing Method
Crossover operator and mutation operator may result in a mass of infeasible solutions, which affects the efficiency and quality of the algorithm. Therefore, the following measures are adopted to solve this problem. Firstly, infeasible solutions are allowed to exist in crossover and mutation operators. If the passengers' walking distance shown by a gene of an infeasible solution exceeds the maximum tolerable walking distance, the passengers' total walking distance of the solution is set to a large constant. If model M1 has no solution, the route length of the infeasible solution is set to a large constant. Infeasible solutions of these two cases will be treated as feasible solutions in the algorithm. Then, the initial population is created under the premise that all solutions in the initial population are feasible. Algorithm 1 can be used to generate the initial population.

Algorithm 1:
Step 1 For i ∈ P 0 , set i = 0. P 0 is the initial population whose size is Z.
Step 2 For solution i, set j = 0 and j ∈ D.
Step 3 Define a set L j , and put each candidate stop whose distance to demand point j does not exceed the maximum tolerable walking distance into L j .

Step4
Choose a candidate stop from L j randomly to serve demand point j. Set j = j + 1, go to Step 3, until j = |D|, then solution i can be created.
Step 5 For solution i, solve model M1, if the optimal solution of M1 exists, set objective value of this optimal solution as the route length of solution i. Otherwise, go to Step 4 to create solution i again.
Step 6 Calculate the total walking distance of passengers for solution i. Set i = i + 1, go to Step 2, until i = Z. Output the initial population.

NSGA-II Based Algorithm
The NSGA-II-based algorithm designed in this paper works as Algorithm 2: Algorithm 2: Step 1 Initialization. Set i = 0, and determine population size Z, crossover probability p c , mutation probability p m , and the number of iterations I. Create the initial population P 0 using Algorithm 1.
Step 2 Apply the fast non-dominated sorting procedure to P i to obtain a series of non-dominated fronts.
Step 3 Perform selection, crossover, and mutation operators for P i , and then the child population Q i with size Z can be obtained. Calculate each solution's route length and total walking distance of passengers in Q i .
Step 4 Combine P i and Q i to generate a new population R i = P i UQ i with size 2Z. Implement the fast non-dominated sorting procedure for R i . Then, a series of non-dominated fronts can be obtained and solutions whose non-dominated rank is r are in Pareto front F r .
Step 5 Elitism. Set P i+1 = ∅, r = 1. Execute Step 5.1-5.3: Step 5.1 If the total number of solutions in F r and P i+1 does not exceed Z, put all solutions of F r in P i+1 .
Step 5.2 If the total number of solutions in F r and P i+1 exceeds Z, calculate crowding distances of solutions in F r . Then, these solutions are put in P i+1 according to the crowding distances from large to small, until the number of solutions in P i+1 is Z.

Step 5.3
If the number of solutions in P i+1 does not exceed Z, r = r + 1, go to Step 5.1. Otherwise, go to Step 6.
Step 6 Apply the fast non-dominated sorting procedure to P i + 1 to obtain a set of non-dominated fronts.
Step 7 Update the Pareto optimal front. Solutions in the first non-dominated front are the latest Pareto optimal front.
Step 8 Judge whether stop. If i = I, output the Pareto optimal front and stop the algorithm. Otherwise, set i = i + 1, go to Step 3.

Computational Experiments
In this section, numerical results solved by the NSGA-II-based algorithm are presented. A randomly generated network example is used to illustrate the efficiency of the algorithm. Moreover, the effects of the maximum tolerable walking distance on the total walking distance of passengers and route length are discussed.

Network Configuration
In the example network, a metro station is located near a community and connects the community with other areas in the same city. There are 20 demand points and 25 candidate stops (the first one is the transfer stop and near the metro station) in the community. The randomly generated distances between demand points and candidate stops and the randomly generated distances between candidate stops are presented as Tables A1 and A2 in Appendix A. The travel demand for the community shuttle service of the demand points (as shown in Table 2) is also known. The minimum route length l min is 3 km, and the maximum route length l max is 12 km. The minimum and maximum stop spacing (i.e., s min and s max ) values are 300 m and 840 m, respectively. The maximum tolerable walking distance d max is 400 m.

Parameters Setting
The key parameters used in the algorithm are given as follows. Crossover probability p c and mutation probability p m are 0.9 and 0.5, respectively. The number of iterations I is set to 500. The NSGA-II-based algorithm is implemented in C# language and ILOG CPLEX 12.4. All tests are executed on an Intel Core (TM) i5 processor at 2.27 GHz under Windows 7 using 4 GB of RAM.

Computational Results
For the example network shown in Section 5.1, the non-dominated solutions obtained by the NSGA-II-based algorithm are presented in Table 3. The computing time is 55.09 min. Considering that we are in the planning phase and this computational experiment is on an actual size, the amount of computing time is reasonable. As can be seen from Table 3, we can obtain more than one non-dominated solutions by running the NSGA-II algorithm once. Decision makers can choose a solution as the final plan of stop location and route design. The demand points served by each candidate stop are shown in Table 4. Figure 3 shows the non-dominated solutions for this sample network visually.
It is important to select optimization objectives reasonably for multi-objective optimization problems. If there is a positive correlation relationship between two minimized objectives, the two objectives are not conflicting, which means that only one objective is effective. It can be seen from Figure 3 that for the optimization objective values of these non-dominated solutions, the shorter the total walking distance of passengers, the longer the route length for the non-dominated solutions, which means that the two objectives are conflicting. Hence, the two objectives are reasonable in this paper.

Effect of the Maximum Tolerable Walking Distance
The main aim of the community shuttle service is to reduce the walking distance of passengers. Therefore, the maximum tolerable walking distance is an important parameter reflecting the service level of the community shuttle system. The total walking distance of passengers and route length are normalized according to Equations (17) and (18), respectively, so that the values of the two objectives are between 0 and 1. Therefore, the two objectives can be described uniformly. We adjust the maximum tolerable walking distance between 300 m and 500 m for the example network represented in Section 5.1. The total walking distance of passengers and route length after normalization are presented in Figure 4.
As shown in Figure 4, when the maximum tolerable walking distance increases from 300 m to 500 m, the total walking distance of passengers increases and the route length decreases. When the maximum tolerable walking distance increases from 300 m to 400 m, the total walking distance grows slowly. When the maximum tolerable walking distance grows from 400 m to 500 m, the total walking distance of passengers increases significantly. When the maximum tolerable walking distance increases from 300 m to 500 m, the route length decreases continuously. If the maximum tolerable walking distance is set to 418 m, the trade-off between the total walking distance of passengers and the route length can be achieved. Therefore, the optimal maximum tolerable walking distance is 418 m. Decision makers can set this parameter to 400 m to make it convenient to be applied in practice. Decision makers can also measure the total walking distance of passengers and route length to determine the maximum tolerable walking distance.

Effect of the Maximum Tolerable Walking Distance
The main aim of the community shuttle service is to reduce the walking distance of passengers. Therefore, the maximum tolerable walking distance is an important parameter reflecting the service level of the community shuttle system. The total walking distance of passengers and route length are normalized according to Equations (17) and (18), respectively, so that the values of the two objectives are between 0 and 1. Therefore, the two objectives can be described uniformly.
where indexes i = 1, 2, 3, 4, and 5 correspond to cases that the maximum tolerable walking distances are 300 m, 350 m, 400 m, 450 m and 500 m, respectively; n WD i is the total walking distance of passengers of the ith case; and f WD i is the value of the total walking distance of passengers after normalization of the ith case.
where indexes i share the same meaning as that in Equation (17); n LL i is the route length of the ith case; and f LL i is the value of route length after normalization of the ith case. We adjust the maximum tolerable walking distance between 300 m and 500 m for the example network represented in Section 5.1. The total walking distance of passengers and route length after normalization are presented in Figure 4.
As shown in Figure 4, when the maximum tolerable walking distance increases from 300 m to 500 m, the total walking distance of passengers increases and the route length decreases. When the maximum tolerable walking distance increases from 300 m to 400 m, the total walking distance grows slowly. When the maximum tolerable walking distance grows from 400 m to 500 m, the total walking distance of passengers increases significantly. When the maximum tolerable walking distance increases from 300 m to 500 m, the route length decreases continuously. If the maximum tolerable walking distance is set to 418 m, the trade-off between the total walking distance of passengers and the route length can be achieved. Therefore, the optimal maximum tolerable walking distance is 418 m. Decision makers can set this parameter to 400 m to make it convenient to be applied in practice. Decision makers can also measure the total walking distance of passengers and route length to determine the maximum tolerable walking distance.

Case Study
This section demonstrates the effectiveness of the formulation and algorithm by applying them to a community in Shanghai, China. As shown in Figure 5, there are nine demand points, one candidate stop (transfer stop, H1) near a metro station, and 12 other candidate stops in the community. The connection between the demand points and the candidate stops and the connection between the candidate stops are shown by the topological diagram in Figure 6. The travel demand for the community shuttle service of each demand point is shown in Table 5 (refer to Wu [36]). Actual data, including the distances between candidate stops and demand points and the distances between candidate stops, are shown as Tables A3 and A4 in Appendix A. The minimum route length is 2 km and the maximum route length is 8 km. The minimum and maximum stop spacing values are 150 m and 1000 m, respectively. The maximum tolerable walking distance of passengers is 500 m. The key parameters used in the algorithm are the same as those in Section 5.2.

Case Study
This section demonstrates the effectiveness of the formulation and algorithm by applying them to a community in Shanghai, China. As shown in Figure 5, there are nine demand points, one candidate stop (transfer stop, H 1 ) near a metro station, and 12 other candidate stops in the community. The connection between the demand points and the candidate stops and the connection between the candidate stops are shown by the topological diagram in Figure 6. The travel demand for the community shuttle service of each demand point is shown in Table 5 (refer to Wu [36]). Actual data, including the distances between candidate stops and demand points and the distances between candidate stops, are shown as Tables A3 and A4 in Appendix A. The minimum route length is 2 km and the maximum route length is 8 km. The minimum and maximum stop spacing values are 150 m and 1000 m, respectively. The maximum tolerable walking distance of passengers is 500 m. The key parameters used in the algorithm are the same as those in Section 5.2.
Symmetry 2018, 10, x FOR PEER REVIEW 12 of 21 Figure 4. Relation between the maximum tolerable walking distance, total walking distance of passengers, and route length.

Case Study
This section demonstrates the effectiveness of the formulation and algorithm by applying them to a community in Shanghai, China. As shown in Figure 5, there are nine demand points, one candidate stop (transfer stop, H1) near a metro station, and 12 other candidate stops in the community. The connection between the demand points and the candidate stops and the connection between the candidate stops are shown by the topological diagram in Figure 6. The travel demand for the community shuttle service of each demand point is shown in Table 5 (refer to Wu [36]). Actual data, including the distances between candidate stops and demand points and the distances between candidate stops, are shown as Tables A3 and A4 in Appendix A. The minimum route length is 2 km and the maximum route length is 8 km. The minimum and maximum stop spacing values are 150 m and 1000 m, respectively. The maximum tolerable walking distance of passengers is 500 m. The key parameters used in the algorithm are the same as those in Section 5.2.    As shown in Table 6, there are two non-dominated solutions obtained by the NSGA-II-based algorithm. The demand points served by each candidate stop are shown in Table 7. The optimal route and stops of the community shuttle system in the topological graph and realistic network are presented in Figures 7-10, respectively.   Figure 6. Topological structure of demand points and candidate stops. Table 5. Travel demand of demand points for the community shuttle service (trip/day). No.
Travel demand  673  553  572  582  542  627  605  614  590 As shown in Table 6, there are two non-dominated solutions obtained by the NSGA-II-based algorithm. The demand points served by each candidate stop are shown in Table 7. The optimal route and stops of the community shuttle system in the topological graph and realistic network are presented in Figures 7-10, respectively.      There is a community shuttle line (i.e., Line 1604) through this community now. Residents in this community often take the section between Gucun Park and Jutai Road-Baotai Road of Line 1604.   There is a community shuttle line (i.e., Line 1604) through this community now. Residents in this community often take the section between Gucun Park and Jutai Road-Baotai Road of Line 1604. There is a community shuttle line (i.e., Line 1604) through this community now. Residents in this community often take the section between Gucun Park and Jutai Road-Baotai Road of Line 1604. The route and stops of this section are shown in Figure 11. The length of this section is 4.8 km. If passengers choose the nearest stops, the total walking distance is 2597.75 km. As can be seen from Table 6, the proposed method generates a longer route than Line 1604. However, for the total walking distance of passengers, the solutions obtained by the proposed method are shorter than Line 1604. Therefore, compared with Line 1604, the solutions obtained by the proposed method are more convenient for passengers. Operators can make decisions by considering the trade-off between route length and the total walking distance of passengers. The route and stops of this section are shown in Figure 11. The length of this section is 4.8 km. If passengers choose the nearest stops, the total walking distance is 2597.75 km. As can be seen from Table 6, the proposed method generates a longer route than Line 1604. However, for the total walking distance of passengers, the solutions obtained by the proposed method are shorter than Line 1604. Therefore, compared with Line 1604, the solutions obtained by the proposed method are more convenient for passengers. Operators can make decisions by considering the trade-off between route length and the total walking distance of passengers.

Conclusions
In a good community shuttle system, a proper compromise exists between the passengers' walking distance and the route length. The passengers' walking distance and the route length are closely related to stop location and route design, respectively. Nevertheless, stop location and route design, although directly related, have been rarely considered simultaneously in literature, which makes the relationship between them difficult to analyze. In this paper, a formulation aiming at minimizing the total walking distance of passengers and the route length is developed to determine the plans of stop location, route design, and the situation of matching demand points to stops simultaneously. The route length constraint, stop spacing constraint, maximum tolerable walking distance constraint, and situation of matching demand points to stops are taken into account. Based on the formulation, a NSGA-II-based algorithm is devised to obtain non-dominated solutions. Results of the computational experiment and the case study show that decision makers can obtain more than one solution through the formulation and algorithm within a reasonable time, which verifies the correctness and effectiveness of the formulation and algorithm. In addition, decision makers can find a compromise between the total walking distance of passengers and the route length and know how long the passengers' walking distance will be reduced if buses travel a longer distance. Therefore, the proposed integrated approach is a significant method for obtaining plans of stop location and route design of a community shuttle service.
In conclusion, contributions of this paper are threefold. First, a multi-objective mathematical formulation is proposed to simultaneously determine the stop location, route design plans of a community shuttle service, and the situation of matching demand points to stops, which makes it possible to study the trade-off between the total walking distance of passengers and the route length.

Conclusions
In a good community shuttle system, a proper compromise exists between the passengers' walking distance and the route length. The passengers' walking distance and the route length are closely related to stop location and route design, respectively. Nevertheless, stop location and route design, although directly related, have been rarely considered simultaneously in literature, which makes the relationship between them difficult to analyze. In this paper, a formulation aiming at minimizing the total walking distance of passengers and the route length is developed to determine the plans of stop location, route design, and the situation of matching demand points to stops simultaneously. The route length constraint, stop spacing constraint, maximum tolerable walking distance constraint, and situation of matching demand points to stops are taken into account. Based on the formulation, a NSGA-II-based algorithm is devised to obtain non-dominated solutions. Results of the computational experiment and the case study show that decision makers can obtain more than one solution through the formulation and algorithm within a reasonable time, which verifies the correctness and effectiveness of the formulation and algorithm. In addition, decision makers can find a compromise between the total walking distance of passengers and the route length and know how long the passengers' walking distance will be reduced if buses travel a longer distance. Therefore, the proposed integrated approach is a significant method for obtaining plans of stop location and route design of a community shuttle service.
In conclusion, contributions of this paper are threefold. First, a multi-objective mathematical formulation is proposed to simultaneously determine the stop location, route design plans of a community shuttle service, and the situation of matching demand points to stops, which makes it possible to study the trade-off between the total walking distance of passengers and the route length. Second, we can know how many bus stops should be set up based on our approach. Third, an efficient NSGA-II-based algorithm is provided, which is suitable for our problem. This algorithm is capable of handling numerical examples of a realistic size.
One of the future research topics is to consider scheduling decisions, such as headway, while planning stop location and route design, for the reason that headway is important to the service level and operating cost. Another topic of future work is to design the stop location and route plan of a community shuttle system, considering more than one metro stations near a community.