Research Article  Open Access
Xiaoling Luo, Wenbo Fan, Yangsheng Jiang, Jun Zhang, "Optimal Design of Bus Stop Locations Integrating Continuum Approximation and Discrete Models", Journal of Advanced Transportation, vol. 2020, Article ID 8872748, 10 pages, 2020. https://doi.org/10.1155/2020/8872748
Optimal Design of Bus Stop Locations Integrating Continuum Approximation and Discrete Models
Abstract
Although transit stop location problem has been extensively studied, the two main categories of modeling methodologies, i.e., discrete models and continuum approximation (CA) ones, seem have little intersection. Both have strengths and weaknesses, respectively. This study intends to integrate them by taking the advantage of CA models’ parsimonious property and discrete models’ fine consideration of practical conditions. In doing so, we first employ the stateoftheart CA models to yield the optimal design, which serves as the input to the next discrete model. Then, the stop location problem is formulated into a multivariable nonlinear minimization problem with a given number of stop location variables and location constraint. The interiorpoint algorithm is presented to find the optimal design that is ready for implementation. In numerical studies, the proposed model is applied to a variety of scenarios with respect to demand levels, spatial heterogeneity, and route length. The results demonstrate the consistent advantage of the proposed model in all scenarios as against its counterparts, i.e., two existing recipes that convert CA modelbased solution into real design of stop locations. Lastly, a case study is presented using real data and practical constraints for the adjustment of a bus route in Chengdu (China). System cost saving of 15.79% is observed by beforeandafter comparison.
1. Introduction
Transit route design problem can be divided into two categories: transit network design and single transit route design [1–5]. Welldesigned transit routes constitute as the basic bricks to the big transit network in many cities for defending the wide spread of roadway traffic congestion. The design of a single transit route mainly concerns the locations of stops/stations and the service headways/frequencies during the operation periods. Being physically inflexible (at least for a short term), transit stop locations affect the service accessibility to potential patrons as well as their experienced level of service in terms of, e.g., commercial speed. On the supply side, the design of stop locations also influences transit agency’s operation efficiency in terms of vehicle fuel cost and vehicle fleet size, for instance.
The transit stop location problem has been extensively studied in the literature. Methodologically, two categories can be identified: discrete models and continuum approximation (CA) models. A majority of studies belong to the discretemethod category. For instance, Vuchic and Newell and Vuchic [6, 7] may be the two pioneering works. They sought to find the optimal interstation spacings of a rapid transit corridor to minimize passenger travel time and maximize number of passengers, respectively. Later on, Gleason [8] developed a set covering approach for locating bus stops. This work was extended by Murray with a hybrid set covering model, which determined the stop locations of an existing route segment as well as the locations of new stops for route extension in unserved areas [9]. Furth and Rahbee [10] optimized bus stop locations from a set of prespecified candidate stops in a bus route of Boston. Similarly, Chien and Qin [11] identified a set of demand points as candidate locations and proposed to minimize the total system cost through finding the optimal number and locations of bus stops. Recently, Ceder et al. [12] integrated the impacts of uneven topography into a bus stop location model to more precisely account for users’ walkingaccess speed and vehicles’ acceleration performance.
In the second category, CA models had been developed as an alternative option in locating transit stops. Instead of based on dozens of location variables, these models were built upon a single stop density/spacing variable or function. This parsimonious property endows the CA model with the highefficient finding of the global optimum solution, or sometimes the closedform solution. The first endeavor in this vein was made by Newell [13, 14]. Later, Wirasinghe and Ghoneim [15] proposed a more general CAbased model for determining bus stop spacing (expressed as a function of location). Hurdle and Wirasinghe and Wirasinghe and Seneviratne [16, 17] analyzed the influence of stop spacing and line length with the objective function of system cost minimization. Medina et al. [18] applied a similar CAbased model to locate bus stops considering multiperiod demand in Santiago, Chile. Mostly recently, Su et al. [19] incorporated environmental factors into CA models for an ebus stop location problem.
The CA models, however, have been criticized being too idealized with unrealistic assumptions, such as a continuous space for locating stops at anywhere along the route. Thus, it is recognized that the designs offered by CA models are not ready for implementation. Endeavors had been made to enhance the applicability of CA models. In Wirasinghe and Ghoneim and Medina et al. [15, 18], the continuous stop density/spacing function was discretized into specific locations via the integral method. Yet their models still lack the consideration of realistic street layout and practical location restrictions, e.g., intersections, bridges, and natural obstacles, where no bus stops should be placed.
This paper intends to fill the gap. We propose an optimization framework that integrates CA models with discrete ones for locating bus stops with respect to location constraint. The idealized design of the CA model serves as input to the discrete model, which accordingly defines a given number of stop location variables and formulates the location constraint. The corresponding problem is a nonlinear multivariate optimization problem. A heuristic solution algorithm is presented to find the optimal solution. To the best of our knowledge, this is the first work connecting CA and discrete models so as to furnish implementationready transit route designs.
The remainder of the paper is organized as follows. The next section introduces the existing CA and discretization models. After that, a novel optimization model is proposed for locating bus stops. In Section 3, the solution method is developed to solve the bus stop location problem model. Section 4 presents numerical studies of various experiments in a hypothetical route and a case study in Chengdu (China). Conclusions are drawn in the final section.
2. Models
Section 2.1 presents the stateoftheart CA model of bus route design, followed by the existing recipes that discretize the solutions of the CA model into real designs. Sections 2 and 3 propose our discretization recipe that offers the improved designs and admit practical constraints on stop locations. Table 1 summarizes the notation used in the paper.

2.1. Continuum Approximation Model
Consider a linear bus route with length km. The daily operation time can be divided into periods, e.g., indicating peak and offpeak hours. For each period , the duration time is denoted by hours. The CA model of bus route design can be expressed as the following minimization problem with the decision variables/functions being headways and stop density (as a function of location , or equivalently stop spacing function ) (Medina et al. [18]):which is subject towhere is the total generalized system cost, which is the sum of bus users’ costs and the agency’s costs. The integrands , , and are patrons’ access/egress time cost, waiting time cost, and invehicle travel time cost at location during period , respectively. The , , and are agency’s distancebased cost (irrelevant to location ), timebased cost, and amortized infrastructure cost during period . The is the maximum vehicle load of the bus route, which is restricted from exceeding vehicle capacity, patrons/vehicle, by constraint (1b). Constraint (1c) is the stop capacity constraint that guarantees the maximum amount of waiting patrons does not exceed the stop’s capacity patrons/stop. Constraint (1d) dictates decision variables/functions being nonnegative. The , , , , , , , and can be approximated by and , of which detailed expressions are referred to our previous work [19] and omitted here for the sake of brevity.
From the firstorder conditions of (1a)–(1d), the following relationship can be derived for the optimal and [18, 19]:where and are interdependent, i.e., and ; , , and are demand functions of boarding density, alighting density, and onboard flow at during period , respectively. The is the maximum headway obtained from the vehicle capacity constraint (1b). The operator in (2a) guarantees no larger than , and thus, vehicle load never exceeds the maximum capacity. The is the minimum bus stop density obtained from the bus stop capacity constraint (1b). The max operator in (2b) guarantees no less than , and thus, the number of waiting patrons never exceeds the bus stop capacity. The detailed expressions of and can be found in Su et al. and Medina et al. [18, 19].
Based on the above analytical results, the efficient algorithm can be readily developed using the iteration method to find the optimal solution (see again in Su et al. and Medina et al. [18, 19]). The solution to (1a)–(1d) is, however, still not real design. The is a continuous function in space, as illustrated in Figure 1, and needs to be discretized into specific stop locations (see the next section for the discretization method).
2.2. Discretization Recipes in Literature
In the literature of CA transit route design models, we found two discretization recipes for translating into real designs, namely, the “midpoint” and “endpoint” approaches, as demonstrated in Figure 2 (Medina at al. and Wirasinghe and Ghoneim, [15, 18]). The underlining logic is straightforward: when the integral of the stop density function yields an integer, one stop should be located in the integral interval, e.g., and in Figure 2. Specifically, the midpoint approach locates the stop in the middle of the integral interval, while the endpoint approach locates at the end, as shown in Figure 2.
The two methods are formulated as follows. First, define as the boundaries that let yields integer , where is the total number of stops along the route. Then, bus stop locations, , are determined by
The endpoint method further includes a default stop at , i.e., . Given above stop locations, we obtain the demand coverage of stop by the left and right boundaries and :
Consequently, discrete system metrics can be computed: e.g., boarding and alighting volumes at each stop by , respectively, and patrons costs and agency costs in the next section.
It is worth noting that although (3) and (4) produce real stop locations, the two discretization recipes have flaws. For instance, they cannot guarantee that the discrete stops are optimally located. This is because both midpoint and endpoint methods neglect the locally nonuniform demand distribution, which apparently impacts the specific locations of stops. In addition, the existing recipes are blinded by ignoring practical location restrictions. The consequence may be improper stop locations that cannot be directly implemented in practice.
2.3. Proposed Discretization Recipe
Other than arbitrarily determining stop locations, we propose a multivariate optimization model to do so and admit constraint of stop locations. Given the knowledge of the total number of stops obtained from the CA model, we accordingly define variables of stop locations, . We also specify the restricted locations to be avoided from being stop locations, . Thus, we can formulate the following optimization problem of minimizing the system cost with respect to :which is subject towhere , , , , and are the corresponding cost items derived based on . Constraint (6b) defines the feasible space of . Constraint (6c) restricts stops from being located in the domain of any restricted locations, i.e., . Constraints (6d) and (6e) are the capacity constraints of bus vehicles and stops.
Computations of , , , , , and and are straightforward, and their expressions are given below:where is the onboard flow passing stop , obtained by , and = max is bus dwelling delay at stop due to passenger boarding and alighting:
Note in (6a)–(6e) that patrons’ waiting cost is discarded from the total system cost because it is irrelevant to stop locations. Also note that in (7a)(7b)(8a)–(8e), the headways take the optimal obtained by the CA model (1a)–(1d).
3. Solution Method
Problem (6a)–(6e) is a nonlinear optimization problem with respect to decision variables, . The interiorpoint of the barrier method can be used to solve this problem. For any inequality constraint in problem (6a)–(6e), we can use a barrier function in objective function to replace the inequality constraint . As problem (6a)–(6e) is a minimization problem, the used barrier function should produce 0, when the constraint is satisfied; otherwise, the barrier function produces . Therefore, the barrier function can be expressed approximately bywhere is a parameter in the approximated barrier function . The larger is the value of , the better is the approximated function. Therefore, an iteration process can be used to update the variable of by increasing the parameter [20]. In this paper, we directly employ the interiorpoint algorithm of ‘fmincon’ function in Matlab 2018a. The “fmincon” function is a builtin program in Matlab to solve the nonlinear problems.
Admittedly, the above solution method does not guarantee a globallyoptimal solution due to the nonconvex nature of (6a)–(6e). Thus, we repeated the solutionsearching procedure 10 times for each instance examined in the following numerical studies. Each time, the optimization started with an initial solution generated from the ideal solution to (1a)–(1d) by randomly adjusting stop locations that validate the space constraints (6b) to the neighborhood area. We found that each repetition of the solution procedure always produced the same final solution and thus reckoned that the global optima were attained. Similar treatment can also be found in Wu et al. and Fan et al. [21, 22].
4. Numerical Studies
To demonstrate the effectiveness of the proposed model, Section 4.1 compares two existing discretization recipes via a variety of experiments in a hypothetical bus corridor. Section 4.2 illustrates an application of the proposed model in a case study of a bus route in Chengdu city (China).
4.1. Experimental Comparisons
Following Vaughan and Cousins [23], we consider an arbitrary demand density function as follows:where is the total demand of the corridor; the distributions of trip origins and destinations, and , are assumed to follow a truncated normal distribution, denoted by and with means of 0 and km, respectively, and variance of being the same. The symmetric setting is purposely made to isolate the findings with regards to the spatial variance of demand. Larger indicates lower spatial variation, and vice versa. Other parameter values are all retrieved from Li [24]. They are summarized in Table 2.

Experiments are conducted under a variety of demand scenarios with respect to km, and passenger/h. Table 3 summarizes the cost savings of our model as compared to that of midpoint and endpoint methods. It is observed that the proposed method always leads to positive system cost saving in a range of 0.12% to 2.95%. Comparatively, slightly more savings are found as against the endpoint approach than against the midpoint approach. Closer observation shows that for the midpoint approach, the comparative savings in terms of user cost and agency cost may be negative, but the average savings remain mostly positive. For the endpoint approach, the savings in user cost are always negative, which is opposite for the agency cost. This result can be partly explained by the additional default stop in endpoint design (see Section 2.2), which leads to less user cost but higher agency cost.
 
and are set to be 10 km and 300 × , respectively; and are set to be 10 km and 300 × , respectively; and and are set to be 10 km and 10 km, respectively. 
Note that the values in Table 3 may look small. This is because the result of the CA model is quite flat at the optimal solution, Estrada et al. [25]. Such benefit will accumulate in the daytoday operation and become substantial for both bus patrons and agency.
4.2. Case Study in Chengdu
We apply the proposed model to bus route no. 3 in Chengdu (China), as depicted in Figure 3. According to our survey, the bus route is of length 18.85 km and the operation time is between 6 : 15 am and 11 : 00 pm on weekdays. The peak period is 6 hours on each weekday. Based on historical records, bus vehicles’ cruising speed during peak and offpeak periods are 20 km/h and 30 km/h, respectively. The headways during the peak period and offpeak period are 3 min and 6 min, respectively. The generation of candidate bus stops follows the policy that the bus stop should keep a minimum distance from restricted locations (e.g., intersections) [26]. Other parameters take the same values in Table 2 as above. Along the bus route contains 43 intersections, whose locations were measured in Google Map and given in Table 4.

In preparation for bus route adjustment, boarding and alighting demand at 35 stops were surveyed on April 10th, 2018. Correspondingly, the density functions of boarding and alighting demand are fitted using spline interpolation, as shown in Figure 4.
Based on the CA model and proposed discretization recipe, we redesign the current transit service. Figure 5(a) presents the optimized bus stop density along the corridor as well as the specific bus stop locations with and without the consideration of location constraint. The stop density ranges from 1.25 stops/km to 4.42 stops/km. After stop density being discretized, 44 bus stops are determined in the corridor. Figure 5(b) exemplifies two stops, i.e., the 10th and 23th stops. The location constraint is violated by the idealized design while the optimal design with location constraint can guarantee the bus stop location satisfies the location requirement.
(a)
(b)
Table 5summarizes the current transit service, optimized transit service with and without considering the location constraint. The results show that the optimized result without considering the location constraint is the best but cannot be used into practice due to the constraint violation. Besides, the cost difference between the results with considering and without considering location constraint is very small. The optimized result addressing the location constraint can save 35.81% agency cost at the expense of 4% increasing passenger cost, forming 15.79% system cost saving when compared with the current transit service.
 
The difference is the calculated by , where are the costs of passenger, agency, and system with the consideration of location constraint, respectively. are the costs of passenger, agency, and system without the consideration of location constraint, respectively. 
5. Conclusion
This paper proposes a modeling framework that connects continuum approximation methods and discrete ones in optimizing bus stop locations. To our best knowledge, this is the first work in the transit route design literature. Our model is no longer limited by the given set of candidate stop locations as the conventional discrete models. Meanwhile, our design outreaches the idealized design of CA models and explicitly addresses practical stop locating restrictions. The proposed hybrid model not only bears the solution efficiency of CA models due to the parsimonious property but also produces implementationready designs as do by discrete models. Numerical studies of various scenarios demonstrate the effectiveness of the proposed model. A case study in Chengdu (China) illustrates how the model is applied to bus line redesign/adjustment in reality.
Of note, the present study still has several limitations. For instance, more realistic concerns (e.g., socioeconomic and political ones) are involved in locating bus stops, which may require further fine tuning. The local conditions (e.g., design and safety) of streets may also influence the decision of bus stop location [27]. To account for these constraints, it is expected to develop a decisionsupport platform based on the proposed modeling framework and integrate other computer aided tools to facilitate designers’ operation.
Data Availability
The boarding and alighting data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This study was funded by the National Nature Science Foundation of China (NSFC 51608455), Sichuan Provincial Science & Technology Innovation Cooperation Funds (2020YFH0038), Doctoral Innovation Fund Program of Southwest Jiaotong University (DCX201826), Chongqing Municipal Transportation Engineering Key Laboratory Open Project (2018TE04), and National Key R&D Program of China (2018YFB1601100).
References
 M.c. Tan, C. O. Tong, S. C. Wong, and J.m. Xu, “An algorithm for finding reasonable paths in transit networks,” Journal of Advanced Transportation, vol. 41, no. 3, pp. 285–305, 2007. View at: Publisher Site  Google Scholar
 G. Laporte, A. Marín, J. A. Mesa, and F. Perea, “Designing robust rapid transit networks with alternative routes,” Journal of Advanced Transportation, vol. 45, no. 1, pp. 54–65, 2011. View at: Publisher Site  Google Scholar
 Q. K. Wan and H. K. Lo, “A mixed integer formulation for multipleroute transit network design,” Journal of Mathematical Modelling and Algorithms, vol. 2, no. 4, pp. 299–308, 2003. View at: Publisher Site  Google Scholar
 C. Ma and D. Yang, “Public transit network planning in small cites considering safety and convenience,” Advances in Mechanical Engineering, vol. 12, no. 1, Article ID 1687814020902351, 2020. View at: Publisher Site  Google Scholar
 C. Ma and X. D. Xu, “Providing spatialtemporal priority control strategy for BRT lanes: a simulation approach,” Journal of Transportation Engineering, Part A: Systems, vol. 146, no. 7, Article ID 04020060, 2020. View at: Publisher Site  Google Scholar
 V. R. Vuchic, “Rapid transit interstation spacings for maximum number of passengers,” Transportation Science, vol. 3, no. 3, pp. 214–232, 1969. View at: Publisher Site  Google Scholar
 V. R. Vuchic and G. F. Newell, “Rapid transit interstation spacings for minimum travel time,” Transportation Science, vol. 2, no. 4, pp. 303–339, 1968. View at: Publisher Site  Google Scholar
 J. M. Gleason, “A set covering approach to bus stop location,” Omega, vol. 3, no. 5, pp. 605–608, 1975. View at: Publisher Site  Google Scholar
 A. T. Murray, “A coverage model for improving public transit system accessibility and expanding access,” Annals of Operations Research, vol. 123, no. 1–4, pp. 143–156, 2003. View at: Publisher Site  Google Scholar
 P. G. Furth and A. B. Rahbee, “Optimal bus stop spacing through dynamic programming and geographic modeling,” Transportation Research Record: Journal of the Transportation Research Board, vol. 1731, no. 1, pp. 15–22, 2000. View at: Publisher Site  Google Scholar
 S. I. Chien and Z. Qin, “Optimization of bus stop locations for improving transit accessibility,” Transportation Planning and Technology, vol. 27, no. 3, pp. 211–227, 2004. View at: Publisher Site  Google Scholar
 A. Ceder, M. Butcher, and L. Wang, “Optimization of bus stop placement for routes on uneven topography,” Transportation Research Part B: Methodological, vol. 74, pp. 40–61, 2015. View at: Publisher Site  Google Scholar
 G. F. Newell, “Dispatching policies for a transportation route,” Transportation Science, vol. 5, no. 1, pp. 91–105, 1971. View at: Publisher Site  Google Scholar
 G. F. Newell, “Scheduling, location, transportation, and continuum mechanics: some simple approximations to optimization problems,” SIAM Journal on Applied Mathematics, vol. 25, no. 3, pp. 346–360, 1973. View at: Publisher Site  Google Scholar
 S. C. Wirasinghe and N. S. Ghoneim, “Spacing of busstops for many to many travel demand,” Transportation Science, vol. 15, no. 3, pp. 210–221, 1981. View at: Publisher Site  Google Scholar
 V. F. Hurdle and S. C. Wirasinghe, “Location of rail stations for many to one travel demand and several feeder modes,” Journal of Advanced Transportation, vol. 14, no. 1, pp. 29–45, 1980. View at: Google Scholar
 S. C. Wirasinghe and P. N. Seneviratne, “Rail line length in an urban transportation corridor,” Transportation Science, vol. 20, no. 4, pp. 237–245, 1986. View at: Publisher Site  Google Scholar
 M. Medina, R. Giesen, and J. C. Muñoz, “Model for the optimal location of bus stops and its application to a public transport corridor in Santiago, Chile,” Transportation Research Record: Journal of the Transportation Research Board, vol. 2352, no. 1, pp. 84–93, 2013. View at: Publisher Site  Google Scholar
 Y. Su, X. Liu, G. Lu, and W. Fan, “Optimal design of bus routes for different vehicle types considering various driving regimes and environmental factors,” Transportation Research Record: Journal of the Transportation Research Board, vol. 2673, no. 6, pp. 299–316, 2019. View at: Publisher Site  Google Scholar
 S. Boyd, S. P. Boyd, and L. Vandenberghe, Convex Optimization, Cambridge University Press, Cambridge, UK, 2004.
 W. Fan, Y. Mei, and W. Gu, “Optimal design of intersecting bimodal transit networks in a grid city,” Transportation Research Part B: Methodological, vol. 111, pp. 203–226, 2018. View at: Publisher Site  Google Scholar
 L. Wu, W. Gu, W. Fan, and M. J. Cassidy, “Optimal design of transit networks fed by shared bikes,” Transportation Research Part B: Methodological, vol. 131, pp. 63–83, 2020. View at: Publisher Site  Google Scholar
 R. J. Vaughan and E. A. Cousins, “Optimum location of stops on a bus route,” in Proceedings of the 7th International Symposium on Transportation and Traffic Theory, Kyoto, Japan, August 1977. View at: Google Scholar
 H. Li, “Optimization and case study of transit lines based on continuum approximation,” Southwest Jiaotong university, Chengdu, China, 2018, Master’s Thesis. View at: Google Scholar
 M. Estrada, M. RocaRiu, H. Badia, F. Robusté, and C. F. Daganzo, “Design and implementation of efficient transit networks: procedure, case study and validity test,” Procedia—Social and Behavioral Sciences, vol. 17, pp. 113–135, 2011. View at: Publisher Site  Google Scholar
 K. Li, Code for Design of Urban Road Public Transportation stop, terminus and Depot Engineering, China Architecture & Building press, Beijing, China, 2011.
 Z. Feng, M. Yang, W. Zhang, Y. Du, and H. Bai, “Effect of longitudinal slope of urban underpass tunnels on drivers' heart rate and speed: a study based on a real vehicle experiment,” Tunnelling and Underground Space Technology, vol. 81, pp. 525–533, 2018. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Xiaoling Luo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.