Estimating the Capacity of Urban Transportation Networks with an Improved Sensitivity Based Method

The throughput of a given transportation network is always of interest to the traffic administrative department, so as to evaluate the benefit of the transportation construction or expansion project before its implementation. The model of the transportation network capacity formulated as a mathematic programming with equilibrium constraint (MPEC) well defines this problem. For practical applications, amodified sensitivity analysis based (SAB)method is developed to estimate the solution of this bilevelmodel. The high-efficient origin-based (OB) algorithm is extended for the precise solution of the combined model which is integrated in the network capacity model. The sensitivity analysis approach is also modified to simplify the inversion of the Jacobian matrix in large-scale problems. The solution produced in every iteration of SAB is restrained to be feasible to guarantee the success of the heuristic search. From the numerical experiments, the accuracy of the derivatives for the linear approximation could significantly affect the converging of the SAB method. The results also show that the proposed method could obtain good suboptimal solutions from different starting points in the test examples.


Introduction
As the rapid growth of urbanization, the population and economy of most cities in the developing countries or regions are going through significant changes.The fast development of transportation infrastructures in these areas gives rise to a quick change on the travel behaviors.Along with the city's expansion, an inevitable thing is to design new transportation system which is capable of meeting the future development in land-using and population growth.However, under fastchanging of travel demand, the conventional forecasting methods (e.g., "four-step") could not provide a straightforward evaluation to the new designed transportation network such as on whether the new network is suitable for the travel demand in the target year or how many trips can be accommodated.Thus, to help the government to make decision on the expansion or construction project or preassessment of social benefits, the throughput, or "capacity, " of the given transportation network is of practical meaning to be estimated before planning implementation.
Capacity is a commonly used property to represent the maximum flows that can pass through the link or node in transportation system.In an attempt to address the question of what is the maximum attainable throughput of the given network, the concept of capacity is employed as an important measurement for transportation system evaluation.It is able to reflect how much traffic demand can be accommodated by a given transportation system.Thus, efficient policy for land use or traffic restraint and growth can be established in advance.According to the conventional network flow theory, the capacity problem is stated to find the maximum flows that can be sent from a specified source node to another specified sink node without exceeding the capacity of any link.This well-known problem is extended to the multicommodity and is widely used in freight transport.However, when modeling the capacity of urban transportation network, the problem becomes quite complex.Noted by Yang et al. [1], travelers in urban transportation network can choose their routes and their trip costs increase with increasing flow as a result of congestion.Besides, multiple origin and destination (O-D) pairs exist and the flows between distinct O-D pairs cannot be exchanged in passenger transportation system.These differences make the modeling of transportation network capacity complex, and the intriguing problem is also hard to solve.
The most popular formulation of the transportation network capacity is the bilevel programming model, which maximizes the traffic flows under the equilibrium constraint.Wong and Yang [2] first incorporated the reserve capacity concept into a traffic signal control network.This concept was widely extended in the study of signal controlled networks [3][4][5].The reserve capacity is defined as the largest multiplier applied to a given O-D demand matrix that can be allocated to a network, so the solution is significantly affected by the predetermined O-D matrix.However, it is unrealistic to assume that all O-D flows increase in the same rate, especially for the areas under rapid changing.If the predetermined distribution proportion is far from the future tendency, the solution will be of little use.Consequently, in order to reflect the differences of the future development of each urban subarea, Yang et al. [1] considered that the new increased O-D demand pattern should be variable both in level and in distribution, while the distribution of the current trips would be relatively fixed.Thereby they introduced the equilibrium trip distribution/assignment model with variable destination costs (ETDA-VDC) [6] to capture this characteristic for network capacity estimation.Based on this model, Kasikitwiwat and Chen [7] proposed the concepts and models of the ultimate and practical capacity.The former is used for the new network without any current flow, while the latter is the same as Yang's model.Then, Chen and Kasikitwiwat [8] used the practical network capacity model to describe the limited flexibility of transportation networks.According to the literatures, the concept of the practical capacity model is more fully functional and preferred, as it takes both the current demand pattern and the variability of future growth into consideration.
In order to solve the various bilevel capacity models, the SAB algorithm is generally employed.This SAB algorithm for bilevel programing was first presented in Friesz et al. [9].It is heuristic and depends absolutely on the derivative information produced by the sensitivity analysis of the lowerlevel problem [10].Benefiting from the rich achievements in the study of the sensitivity analysis for equilibrium models [11][12][13][14][15][16][17][18], SAB algorithm has been widely utilized in the optimization problems of equilibrium network flows, such as traffic signal control [3,5,19,20] and network design [21][22][23], as well as network capacity [1,[3][4][5]24].But due to the difficulty of the sensitivity analysis for ETDA-VDC model, Yang et al. [1] used an iterative estimation-assignment (IEA) algorithm [11] instead to solve the transportation network capacity problem.Later, Kasikitwiwat and Chen [7] and Chen and Kasikitwiwat [8] selected using a genetic algorithm to solve the problem in a very small network.However, since the complexity of the network capacity problem, the global optimization algorithms (e.g., genetic algorithm or simulated annealing) can hardly find the exact solution to the capacity problem in larger networks, and the computation time could be intolerably long.By contrast, SAB algorithm has the property of fast convergence which makes the computation terminate at a local optimum within a considerable time.Nevertheless, the calculation issue of the matrix inversion still limits the applications of the SAB method.To address this problem, we developed an effective method by simplifying the matrix inversion in the sensitivity analysis approach, which will take much less memory space, so the capacity of the real transportation networks could be estimated.
In this study, Yang's formulation [1] for the transportation network capacity model is employed to describe the practical capacity of the urban road system.In an attempt to estimate the capacity of the real road networks, a series of improvements are taken to the SAB method to make the heuristic search successfully converge to a relatively better suboptimal solution.Firstly, the OB algorithm [25] is modified for the solution of the lower-level ETDA-VDC model.Then, the restriction sensitivity analysis approach for the ETDA-VDC model [18] is employed in and improved on the expressions so as to deal with the large-scale problems.Besides, the solution update strategy is modified on the step-size adaption, which ensures the entire heuristic search to converge to a local optimum.Finally, numerical experiments are implemented to show the efficiency and capability of the proposed SAB method.

Road Network Capacity Model
It should be noted that the boldface type of the Notation section represents the corresponding column vectors in the remainder of this paper.

Model Formulation.
Conversional methods, like the reserve capacity model, evaluate the capacity of transportation networks by assuming that the travel demand increases with a determined distribution proportion, which is usually far from the regularity and underestimates the results.In order to evaluate zonal development potential and equilibrium network capacity more appropriately, Oppenheim's definitions on the behaviors of the existing travel demand and the additional demand are introduced as follows.
(i) The existing demand, denoted by   , has predetermined origins and destinations.The pattern of the existing demand is formed during the past long term, so its distribution is going to be relatively stable and can be regarded as fixed.The existing demand only changes routes to optimize the travel cost.
(ii) The additional demand,   , is variable.The new generated demands from residential area can decide their daily travels without the constraints of either destination or route choices.But the behavior of the additional demand still follows the rule that selects the destinations which maximize the "utility" of the trips.
The utility could include the destination attractiveness, the cost along traveling route, and other factors.The attractiveness of destination is determined by the congestion at destination and the expenses for the activity in that area.
In network capacity model, the utility from origin  to destination  is formulated as   = −(  +  ), in which   is the travel cost from  to , and   denotes the destination cost which could be a decrease function of the total additional trip attraction,   , at destination .Besides, the destination choices of the travelers at each origin  are assumed to have certain randomness.Thus, the conditional probability that an individual will choose destination  is derived by using the standard logit function, so the O-D travel demand is conducted by Thus, with the objective to maximize the additional demand under the above travel behavior regularity and certain physical constraints, the typical road network capacity model is formulated as the following bilevel programming problem.
Upper-level problem is as follows: where   (O) and V  (O) are obtained by solving the ETDA-VDC problem in lower-level problem.
Lower-level problem, ETDA-VDC model, is as follows: ≥ 0, ∀ ∈ ,  ∈   ,  ∈   , The upper-level problem defines a maximal trip production model.The objective is to maximize the summation of the additional trip production at origins.Equation (3) represents that the traffic flow on every link should not exceed its capacity.Constraints (4) and ( 5) are the limitation of the zonal trip production and attraction.They mean the number of trips generated and attracted at each traffic zone should be limited by some upper bounds, namely,  max  and  max  , respectively.
The lower-level problem is the ETDA-VDC model.The objective function (7) indicates the choice behavior of both the existing and additional travel demand.Constraint (9) shows that the amount of the existing flows is fixed for each O-D, while constraints ( 8) and (10) show that the additional flows are only restrained at the origin productions.The relationship between the link flow and route flow is represented in (11).All the variables must be nonnegative, that is, constraints ( 12)-( 14).The lower-level problem is a combined distribution and assignment model.
This bilevel model was first presented in work by Yang et al. [1].Because of the advantages on the formulation of the travel demand growth, it was continually used in later researches as a typical model for the road network capacity concept.The remaining part of this study focuses on the solution of this model in the real-sized road networks.

Sensitivity Based Heuristic Algorithm
This section presents an improved version of the SAB algorithm for the solution of the road network capacity model.To overcome the drawbacks of the conventional SAB algorithm [9] that cannot be applied to any real-sized network for capacity estimations, the following improvements are carried out.
(i) The lower-level ETDA-VDC model is fast solved by a modified OB algorithm to produce a high accurate solution.
(ii) The rectified sensitivity analysis method for the ETDA-VDC model is simplified on the calculation of the matrix inverse to be applicable for large-scale problems.
(iii) The solution update of the heuristic search is improved by step-size adaption in order to ensure that the SAB algorithm can converge to a local optimum.
Correspondingly, a series of techniques is proposed in this section, so our modified SAB algorithm will be capable to solve the bilevel road capacity model efficiently.

Origin-Based Algorithm for the ETDA-VDC Model.
In the standard SAB search, the lower-level ETDA-VDC model should be solved in every iteration to conduct an equilibrium traffic flow pattern, namely, the solution to the lower-level problem.According to recent researches on the traffic assignment problems, the OB algorithm is demonstrated to be one of the state-of-the-art algorithms [26].In addition, from the results of the OB algorithm, the set of all equilibrated routes can be easily extracted, which will be a precondition for the restriction sensitivity analysis approach in next step of SAB algorithm.
The OB algorithm uses the origin-based approach proportions,   = {  }, as its main solution variables, where   represents the proportion of flow that comes from origin  through link .The approach proportions are updated by shifting flows within the restricting subnetwork,   .The details of the OB algorithm can be referred to in Bar-Gera [26].Using the approach proportions, the OB algorithm is able to store the route flow information with significantly less memory than the route-based algorithms and achieve high-accuracy solutions compared to the link-based algorithms.
In the process of the OB algorithm, it starts with trees of minimum cost routes from origins as restricting networks.Then steps of updating restricting network and approach proportions are implemented for each origin separately.An "inner loop" is performed subsequently to accelerate convergence, in which the origin-based link flows are updated while keeping the restricting subnetworks fixed.The inner loop is useful, because the restricting subnetworks tend to stabilize fairly quickly but updating the restricting subnetworks requires more computational effort [26,27].
The OB algorithm for the combined trip distribution and assignment problems requires an additional step to update the O-D flows while keeping the route proportions fixed.Taking Evans' algorithm as a reference, we modified the original OB algorithm for the solution of the ETDA-VDC model, in which the step size of the O-D flow update is obtained by solving a one-dimension search problem.The proposed valgorithm is summarized as shown in Algorithm 1.
The performance of Algorithm 1 was testified on different test networks.The characteristics of these networks are obtained from http://www.bgu.ac.il/∼bargera/tntp/.The converging processes on two well-known networks are demonstrated in Figure 1, where the relative gap is defined as RG = |  −  −1 |/  (  is the value of objective function (7) at th main loop).It shows that the converging of the improved OB algorithm (IOBA) is faster and more stable than Evans' algorithm (EA) or the original OB algorithm (OBA) steps in Bar-Gera and Boyce [25].

Sensitivity Analysis of the ETDA-VDC Model.
On the basis of an equilibrium solution of the lower-level model, the derivatives of the lower-level decision variables with respect to the upper-level ones should be produced for linearly approximating the whole bilevel model in next step.The derivative expressions can be obtained by employing the restriction sensitivity analysis approach for the ETDA-VDC model [18].The restriction approach reduces the original network to a restricting one in which the nonuniqueness difficulty will be overcome, and thus the full derivatives of the link flows in the primal problem can be solved.The necessary results from the sensitivity analysis of the ETDA-VDC model are presented in this section without proof.

Initialization:
Set  := 0, determine the initial link cost t 0 and destination cost M 0 by setting x 0 = 0 and d 0 = 0, respectively.For each origin  ∈  do (i) Find the tree of minimum cost routes rooting from .Let   be the minimum route tree.Denote the minimum route cost to destination  ∈   by  0  , and choose a minimum cost route from p to q. (ii) Compute the initial variable O-D demands by , ∀ ∈   (iii) For  ∈   , assign the entire O-D demand   and  0  to the minimum cost route  from  to , and obtain initial link flow  0  .(iv) Update the link costs using the initial link flows.(v) Initialize the origin-based approach proportions  0  .

Main Loop:
Given the current variable O-D demand d  obtained in ( − 1)th-iteration:  In the restricting problem the derivatives of the model solutions with respect to the input parameters are derived from where J x is an invertible Jacobian matrix with respect to the solution variables x = (f  , h  , d  ,   ,   , u  )  , in which f  and h  are the basic variables of the additional and existing route flows, respectively.J O is the Jacobian matrix with respect to the additional zonal productions, O. Here, , , u are the Lagrange multipliers associated with the conservation constraints ( 8)- (10).O is referred to as an input parameter for the ETDA-VDC model.It is in terms of the perturbation .The expressions of J x and J O are as follows: where ) is a diagonal matrix with 1/  as its diagonal elements.The superscript "" represents the transposed matrix."" indicates the identity matrix.The link cost function, t(⋅), should be strongly monotonically defined in v, which guarantees the uniqueness of the equilibrium solution and is also a precondition for the sensitivity analysis approach.Utilizing the restriction approach, the derivative, ∇ O x(), is first conducted in the restricting problem and then ∇ O k, which is defined in the original one, can be obtained by Equation (18) indicates that the variations on the link flows can be just represented by the changes on the basic route flows, which is the rationale of the restriction sensitivity analysis approach.From ( 15), the inverse of the matrix J x is required to be worked out.However, the calculation of the inverse matrix is costly for both computation time and storage space in numerical computation for large-scale problems.The dimension of J x (from ( 16)) could be about five times of the number of O-D pairs.So the elements of J x could totally take approximately (25 4 ) space, where  is the amount of the traffic zones, and thus even a medium network with hundreds of origins will conduct a Jacobian matrix with over one hundred thousand dimensions.For example, in a small network with about 30 traffic zones (900 O-D pairs), nearly 400 MB RAM should be expended on the storage of the Jacobian matrix (if each element uses a double type).Consequently, the conventional sensitivity analysis approach for the combined model is very difficult to be utilized in the large-scale problems.Besides, the direct calculation for the inverse of a large matrix is pretty inefficient and inaccurate in practice.
To avoid the defects of deriving the inverse matrix directly, the usual way to solve the following linear equation as an alternative.Equation (15) which produces the sensitivity results can be rewritten as Then, ∇ O x is derived by solving the following series of linear equations: where b  is the th column vector in matrix of [−J O ] and z  is the th column vector in ∇ O x.Let e  be a unit vector with one in the th position and zeros elsewhere.The length of e  is equal to the column number of J O .Thus, where 0 in bold is zero vectors associate with the zero block in J O .
To solve (20), J x is premultiplied by the matrices K 1 and K 2 in sequence and postmultiplied by K  1 and K  2 simultaneously, which is equivalent to making the elementary transformation of matrix J x .Consider Thus, By using the above operations on (20), we obtain Obviously, Therefore, it is equivalent to solving the equation in terms of y and then deriving z  by Since it is only concerned with the first three subvectors of z  in the sensitivity result, we just need to find the value of y 1 , y 2 , y 3 , and y 6 .From (24), the following group of equations can be obtained: The above equation systems can be further simplified depending on whether the incidence matrix Λ h is square.
(1) If Λ h is square, consider the following.
Thus, Λ h will be constructed as an identity matrix; the value of y 2 and y 3 can be computed from ( 29) and (31), and thereby y 2 = y 3 = 0. From (32), Substituting ( 33) into ( 27), we obtain Combined with (30) we get where Because the matrix [ Λ f ] are leaner independent from the restriction approach.It can be easily conducted that matrix (36) is invertible.Let Solving (35), we can get To this extent, the only effort to produce the desired derivative results is to calculate y 1 , in which the inverse of matrices A and BA −1 B  should be derived.The dimension of matrix A is equal to the numbers of the basic routes used by the additional travel demand, which is a little more than the number of O-D pairs from our observation in computational experiments.The dimension of matrix BA −1 B  is the same as the number of origins and is not large.Consequently, by utilizing (33) and y 2 = y 3 = 0, the first three subvectors of z  can be obtained, which are separately corresponding to the derivative results, ∇ O f  , ∇ O h  , and ∇ O d, in ∇ O x.The expressions are Thus, (2) If Λ h is not square, consider the following.
In this case, y 2 = ∇ O h  will not be equal to zero.For further simplification, we rewrite Λ h as Let   denote the set of routes associated with Λ  h and the flows on these routes are represented by h  .Repeating the derivation from (33) to (35), we can get where the superscript "  " indicates the matrices and variables corresponding to the routes in   .Considering that the set of   is fairly small for network capacity problems, the equation system is slightly larger than (35), which is also easy to be solved.The derivative of the existing route flows is computed by y  2 .Consider And thus, From the above result, it can be noted that the primal problem to calculate the inverse ofthe entire matrix J x is replaced by solving a small system of linear equations.

Restriction of the Equilibrium Network.
When the restriction sensitivity analysis approach was proposed [10], the restricting problem was conducted by finding an extreme point of the feasible region of the equilibrium route flow.An equivalent linear programming was provided for solution.However, this method has been verified usually failed, since the linear programming may find a degenerated solution which causes the fact that the reduced restricting problem cannot reserve the essential information and conducts a wrong result [16].A rectified method is to find a maximum set of the linear independent columns in the coefficient matrix of the following equation system [18]: The above equation system defines the feasible region of the equilbirum route flow under unique euqilibrium solution of v * and d * .Note that the notations are distinguishing from the ones in (16).f 0 and h 0 correspond to the routes that carry positive flows, which are the equilibrated (or minimumcost) routes according to Assumption 2 in Du et al. [18], and can be obtained from the results of the OB algorithm.The superscript "0" is associated with the variables in a reduced problem, which only consists of the equilibrated routes.
From the coefficient matrix of (47), the maximum set of the linear independent columns corresponds to the set of the equilibrated and linear independent routes (ELI, denoted by   ) and can be found using the row echelon form of the coefficient matrix.We can utilize the blocks Λ f and Λ h which are incidence matrices with one in each column (because one route can only serve one O-D).Since there may be more than one equilibrated route between an O-D pair, Λ f and Λ h are generally not column full rank.Thus, the coefficient matrix can be rewritten as where Λ f and Λ h are square matrices, which consist of exactly one equilibrated route for each O-D, so they are equal to identity matrices.Thus, the question is reduced to find and eliminate the linear dependent columns in a submatrix.Consider The columns of the above matrix are only related to the alternative routes between every O-D.Therefore, the set   will include the routes corresponding to Λ f , Λ h , and the maximum set of the independent columns in (49).

Solution of the Maximum Trip Production Problem.
When the sensitivity results are derived, the SAB algorithm will use this information to represent the implicit relationships between v, d, and O. Thus, the bilevel problem is firstorder approximated at the given point, Let k(O * ) and d(O * ), respectively, denote the solutions to the lower-level model at O * .The relationship can be represented by using the Taylor expansion: where the derivatives ∇ O k and ∇ O d are obtained from the sensitivity analysis of the ETDA-VDC model.Therefore, the bilevel problem can be reformulated as The solution of this linear programming can be easily derived using the simplex method.However, because the linear programming problem is just locally approximated, the new solution might be infeasible to the original problem.At an infeasible point of O * the lower-level trip distribution and assignment results may not satisfy the upper-level constraints.In extreme cases, for example, some link flow V *  may be much greater than its capacity   , which could cause the capacity constraints to fail for any nonnegative O.In consequence, the new linear approximation will have no solution.
In consideration of this flaw, we modified the conventional SAB algorithm by restraining the solution to be always located within the feasible region.A trip distribution and assignment step is implemented at the solution of the above linear programming, O * , and then the results are checked with the capacity conditions.If any capacity constraint is violated, the new solution should be updated by a convex combination of the solution from last iteration, O () , and O * () .Therefore, let Ô be the solution to the approximate linear problem above; the maximum step size  = 2 − ( = 0, 1, 2, . ..) is chosen which ensures the capacity constraint is satisfied at any link or traffic zone.The new solution is updated by (52)

A Heuristic Solution
Process.In summary, the modified SAB algorithm involves an iterative process between the

Optimal solution
Step 5 Convergence criterion

Capacity check
Solving the lower-level EDTA-VDC model

Sensitivity analysis Restriction approach
Local linear approximation Simplex method Lower-level distribution and assignment and check Step-length setup

Updating solution Yes No
Step 1 Step 2 Step 3 Step 4 No Yes Initialization Initial feasible solution • Zonal production O (0)   Step 0 upper-level and lower-level problems and can be summarized as follows.
Step 1 (solving lower-level problem).Solve the ETDA-VDC model in lower-level for the given O () , by which the equilibrium link flows v () and O-D demand d () are obtained.The set of all equilibrated routes can be derived at the equilibrium point.
Step 2 (sensitivity analysis).Calculate the partial derivatives, ∇ O k and ∇ O d, using the sensitivity method for the ETDA-VDC model.
Step 3 (local linear approximation).Formulate local linear approximations of the upper-level capacity constraints using the derivative information and solve the approximate linear programming problem to produce an auxiliary trip production Ô() .
The above process is a modification of the conventional SAB approach.The flowchart is shown in Figure 2. It should be noted that without the step-size adaption in Step 4 the SAB method only works in very small networks but could fail in larger examples.The set of all equilibrated routes is used to apply the restriction approach for sensitivity analysis.In this study, the modified OB algorithm for the ETDA-VDC problem is used for this purpose.In addition, the derivative results of the sensitivity analysis can give a precise local linear approximation to the upper-level capacity constraints, which is very important to make the heuristic search converge to a considerable good solution.
In addition, given the nonconvexity of the bilevel problem, the SAB algorithm will converge to a local optimal point [2].However, as shown in Yang et al. [28], for MPEC models, if the upper-level objective function is a linear function of its decision variables, the heuristic algorithm can at least find a noninferior optimal solution.So the SAB algorithm is supposed to be able to find a satisfying solution for the road network capacity problem in practice.

Numerical Experiments
Numerical experiments are conducted in this section.In the numerical experiments the link cost function employs the Bureau of Public Roads (BPR) function: The destination cost function is defined, referring to Chen and Kasikitwiwat [8], as follows: where   is a scaling factor between demand and service cost;   is a dimensionless parameter related to the severity of congestion; and   represents a fixed attraction at destination .
The experiments were implemented on Intel Core i5 CPU 3.20 GHz, 4 GB RAM, using the Microsoft Windows 7 operating system.All of the coding was carried out in Visual C#.The solution accuracy is measured in each iteration  by the relative error of the trip productions O  (RE) as follows:  1.The stopping tolerance is set to be 10 −7 .
From the tests, although SAB algorithm shows a slower speed to reach a local optimum, it always produces much better solutions compared to IEA algorithm.In the tests on the Anaheim network, IEA has not converged and tended to zigzag after a few cycles of computation.Since IEA can quickly produce the approximate derivatives at current point, the CPU time taken in each iteration is much less than SAB algorithm.However, IEA could easily stop at a nonoptimal solution or not converge.This indicates that the accuracy of the derivative results has a significant influence on the convergence of the heuristic search for the network capacity model.The converging processes of the two algorithms on the two networks are illustrated in Figure 3.Note that The final solutions from the five computations are shown in Figure 5.The experiment result in Figure 4 indicates that the solutions by our SAB algorithm can reach a sufficient high precision.All the solutions in Figure 5 are very close, which reflects that the SAB algorithm is comparatively stable for the network capacity problem.Consequently, the proposed SAB algorithm is testified to probably converge to a suboptimal solution from any appropriate feasible start point for the network capacity model.Since the bilevel problem is not convex, the suboptimal solution can be regarded as a satisfying approximation.This solution might be quite close to the global optimum according to the aforementioned characteristics of the network capacity model.Since it only takes a few (less than 30) times of iterations to reach the local optimum, the whole computation of the SAB algorithm can be completed in a considerable CPU time, which is quite superior to the global optimization, for example, genetic algorithm.Thus, the modified SAB method could be a valuable tool for the estimation of the road network capacity in practice.

Conclusions
This paper has presented an effective SAB method for the solution of the transportation network capacity problem.An OB algorithm for ETDA-VDC problem is presented to solve the lower-level problem rapidly and efficiently.The primal problem of sensitivity analysis has been simplified, so the inverse calculation of a full-size matrix has been converted to solving corresponding linear systems of equations.To ensure the heuristic search can proceed properly, the solution from the local linear approximation is rectified to be restrained in its feasible region.Based on these improvements, the SAB method is able to be applied to estimate the capacity of the large-scale real networks.The performance of the proposed SAB method has been demonstrated on two experimental networks.
Based on the method in this paper, the network capacity model could be used to provide an ideal travel demand pattern for a given transportation system, by which the traffic resources can be utilized in maximum.From the model results, the traffic engineers can evaluate the design scheme of a transportation network more accurately before its construction.The results can also be used as a reference for the layout of land use in the new developing areas.In addition, the traffic flow pattern at full saturation may reflect the potential bottleneck of the given network.Further applications of the network capacity model could be detected based on the proposed estimation method.

Figure 2 :
Figure 2: Flowchart of sensitivity analysis based algorithm for network capacity model.

D flows, retain origin-based approach proportions
1 to number of main iterations ( Main ) for each  in  do Update restricting network   Update origin-based approach proportions   end for :Given the origin-based approach proportions    , link flows v  and link costs t  obtained in the steps above: (i) For each  ∈ , compute the destination cost    by O-D demand    and   .(ii) Find the set of auxiliary trip demands d  by solving the following logit distribution model:

Table 1 :
Computational result test networks.Since the characteristic of the heuristic search, little can be said theoretically about the convergence of the SAB algorithm.In this study, we implement SAB algorithm from different start points to observe the converging process.On both the Sioux Falls and Anaheim networks, SAB algorithm can find a sufficient good solution given any arbitrary initial solution.We select the converging processes from five different start points on Anaheim network and plot the results in Figure4.All the computations were implemented in 20 times of iterations.
Set of nodes in the network   : Restricting subnetwork for origin    : Capacity of link    : Existing trip attraction at destination    : Additional trip attraction at destination   max  : Upper limit of trip attraction at destination    : Cost at destination    : Existing trip production at origin    : Additional trip production at origin   max  : Upper limit of trip production at origin    : Set of routes from  to  : Set of zones in the network   : Set of all destination for origin    : Set of all origins for destination  : Link index,  ∈    : Additional demand from origin  to destination    : Existing demand from origin  to destination     : Flow on route  from  to  associated with elastic demand ℎ   : Flow on route  from  to  associated Route index,  ∈     : Travel cost on link  and the function of all link flows x  free  : Free flow travel time on link  V  : Flowonlink V  : Traffic flow on link  from origin    : Approach proportion of link  from origin  : Impedance parameter for trip distribution   : The minimum route cost from  to    , : Link/route incidence indicator, 1 if link  is on route  from origin  to destination ; 0, otherwise Δ: Link/route incidence matrix; Δ f and Δ h correspond to route flows f and h, respectively Λ: O-D/route incidence matrix; Λ f and Λ h correspond to route flows f and h, respectively Φ: Origin/O-D incidence matrix Ψ: Destination/O-D incidence matrix. :