Research Article  Open Access
Mostafa Ameli, JeanPatrick Lebacque, Ludovic Leclercq, "Computational Methods for Calculating Multimodal Multiclass Traffic Network Equilibrium: Simulation Benchmark on a LargeScale Test Case", Journal of Advanced Transportation, vol. 2021, Article ID 8815653, 17 pages, 2021. https://doi.org/10.1155/2021/8815653
Computational Methods for Calculating Multimodal Multiclass Traffic Network Equilibrium: Simulation Benchmark on a LargeScale Test Case
Abstract
This study reviews existing computational methods to calculate simulationbased dynamic network equilibrium. We consider a tripbased multimodal approach for the dynamic network loading. Mode and path choices are carried out at the same level; therefore, travel times depend on the travel path and the mode attributes of travelers. This study develops a multiclass model with several parameters per class. Two different categories of algorithms (heuristic and metaheuristic) are considered in order to solve the discrete dynamic traffic assignment (DTA) problem. Finally, we analyze the equilibrium in a largescale multimodal DTA test case (Lyon 6th + Villeurbanne) in order to investigate the performance of different optimization approaches to solve tripbased DTA. The results show that, in a multimodal and heterogeneous setting, the metaheuristic methods provide better solutions than the heuristic methods in terms of optimality and computation time. These improvements are even more significant than in a homogeneous setting.
1. Introduction
Network user equilibrium is the situation where all users’ travel costs in all routes used are equal and lower than those related to the not chosen route [1]. Dynamic traffic assignment (DTA) is a process to assign a given set of trips to a possible path (set of routes and modes) following the network user equilibrium by considering the dynamics of the transportation system. In a largescale network, when we have numerous paths per origindestination (OD) pairs and also large numbers of origins and destinations, the problem becomes almost intractable analytically because of multiple flow exchanges at nodes. The simulationbased approach, introduced by [2], is an effective way to address the question of DTA in largescale networks [3]. Simulationbased DTA models rely on a network performance module called dynamic network loading (DNL). In DNL, the dynamic traffic simulator does the network loading to calculate network states and travel times; then the traffic assignment component determines the route and path flow distribution [4]. This combination, for largescale networks, allows us to consider sophisticated traffic models in order to incorporate disaggregate information into the estimates of travel costs [5, 6].
Traffic simulators can be divided into two classes: flowbased models, which consider the flow of each path, and tripbased models, which define how many travelers take each path. Macroscopic traffic flow models fall into the first category, while microscopic models belong to the second. The macroscopic approach and flowbased models are usually fast in equilibrium calculation as the path flow discipline is more flexible (flows are not necessarily equivalent to traveler units). However, they are less realistic for OD pairs without adding integrality constraints with low demand as vehicles are split into small particles in practice. In this study, we decide to focus on the tripbased approach in which each traveler is reproduced individually. The tripbased approach can be considered as a multiagent system, wherein each trip impacts the system by its route (including path and mode(s)) choice. Microscopic traffic simulators are now widely used for operational studies, and we have chosen to focus on DTA performance for this kind of model. Tripbased DNL attempts to assign particlediscretized timedependent origin/destination flows in a dynamic network equilibrium framework [7]. The goal of this study is to design and solve a largescale tripbased traffic network equilibrium model.
The user equilibrium situation, theoretically, is the solution of a fixedpoint problem [8]. From a practical point of view, it means that each user perceives no delay or extra cost compared to other users with the same origin and destination [9]. This refers to minimizing the total gap between users’ travel cost of each OD pair and the corresponding shortest path (the path(s) with minimum travel cost) travel cost [10]. In the context of DTA, the Method of Successive Average (MSA) algorithm remains by far the most widely used solution method in simulationbased DTA [11]. The simplicity of the MSA implementation and the nonrequirement of derivative information are the main reasons for its widespread use [12]. However, Yang and Jayakrishnan [13] showed that the convergence is not guaranteed to be the user equilibrium solution with the traditional MSA algorithm. Based on this indication, several heuristic approaches are proposed to improve the determination of the step size (e.g., [14–16]). Ameli et al. [17] did a benchmark on most of the heuristic methods in the literature. Here, in addition to heuristic algorithms, we also considered metaheuristic approaches (proposed by [18]) in order to compare different methods to solve the DTA problem.
The travel cost is the criterion for each user to select the travel route and mode(s). It is important to consider generalized travel costs in simulationbased DTA in order to include different features of travelers’ behavior. Generally, travelers consider different criteria in selecting their optimal properties of travel, including travel route and travel mode(s). Thus, the total demand of the system can be clustered into different user classes based on the definition of their travel cost function. Many attempts have been made on modeling travelers’ behavior in the traffic network. Most multiclass flow models classify travelers by identifying traveler attributes. In the literature, this classification process identifies eight categories of information about travelers: (1) choice model, (2) cost function, (3) economic attributes, (4) generalized cost function, (5) knowledge level of the network, (6) risk taking attitude, (7) social class, and (8) vehicle class. Moreover, some research uses hybrid classification (see, e.g., [19–21]). In this study, we attempt to define a generic cost function for users in order to represent their behavior accurately. For the convenience of the reader, we discuss all categories and identify the different approaches in Appendix A (Table 1) in order to clarify our contribution compared to the previously mentioned papers.

Multiclass DTA models are widely investigated in order to take into account the heterogeneity of the users [22]. While the model becomes more accurate for representing the users’ travel cost function, finding the equilibrium becomes more difficult computationally. There are many studies about developing algorithms to find user equilibrium [9]. This study aims to find out which category of solution algorithms is more efficient for the multiclass DTA problem. There are two main categories of solution methods to solve simulationbased DTA with computational methods: heuristic methods based on fixedpoint theory and metaheuristic methods that have been recently introduced [18]. In the literature, these methods are applied individually to networks with different sizes and different levels of details about users and network dynamics [23]. In this study, we apply solution methods with the multiclass model to a real test case (Lyon 6th + Villeurbanne: about ) in order to benchmark all the solution methods. The goal is to answer which solution method is the most efficient in generalized and advanced numerical simulation settings.
In the next section, the multiclass model is presented. The Wardrop condition for the network equilibrium model is discussed for tripbased DTA in Section 2. The benchmarked algorithms that find the network equilibrium are explained in Section 3. In Section 4, the numerical experiment and the simulator are presented. Section 5 presents the numerical results and discusses them based on several indicators. Finally, we conclude this study in Section 6.
2. Dynamic Traffic Assignment Model
2.1. Multiattribute Generalized Travel Cost
In order to describe the traveler’s behavior, a generalized cost function that takes into account four types of information mentioned in Section 1: (2), (3), (4), and (7), is used. Typically, the generalized cost function will integrate travel time and travel cost . Also, will take into account the classdependent value of time (VOT). Zhan and Ukkusuri [24] denoted for class and a vector of economic and social class parameters. Therefore, we need to define, for class and path (: set of paths), the path generalized cost .
The vector of parameters includes the traveler mode(s) and social attributes. For instance, one class in the network can be students who pay less to use the transit network and who have the specific value of time. The demand of OD pair (: set of OD pairs) is fixed andwhere is the user class index , is a set of paths for , and is the number of users of class on path . After assigning the user to the path between OD pair , (the number of users on link ) iswhere is the link index (: set of links in the network), is the number of users of class on link , and if path contains link and 0 otherwise. The travel time and cost functions depend on the congestion of the link. Therefore, both functions are calculated by the following formulas:
and, consequently, the link generalized cost function is calculated as follows:
Equation (6) calculates the generalized cost of each link for each class of users. Finally, the cost of path is calculated by the following equation:
Now, we can build a DTA model based on the generalized cost defined by equation (7).
2.2. Equilibrium Conditions
Consider a network with a finite set of nodes and a finite set of directed links . The demand is timedependent. The period of interest (planning horizon) of duration is discretized into a set of small time intervals indexed by ( and ). is the duration of the time intervals. In an interval , travel time and traffic conditions do not change. The important notations to introduce the dynamic equilibrium model and optimization algorithms are as follows.
For each class , for all OD pairs and for all paths , the dynamic traffic network equilibrium condition with fixed travel demands for the aforementioned traffic network equilibrium problem is
Lu et al. [10] reformulated the problem as a nonlinear problem in order to minimize the gap function. The gap function is defined as the gap between average path travel time and the shortest path travel time. Therefore, the solution of this problem is equivalent to finding the solution to the following variational inequality:where denotes the flow constraints based on . Note that , a feasible assignment, is the optimal assignment, if Equation (9) holds for all feasible . In Equation (9), both and are decision variables and hence the gap function is a function of both variables. Existence of solution(s) to Equation (9) requires the continuity of the generalized cost functions of paths with respect to the numbers of users on paths. In this case, the unicity of the solution requires the strict monotonicity of the OD paths generalized cost functions with respect to OD path flows (i.e., the number of cars) [25]. In general urban transportation networks, these conditions do not hold [22].
3. Solution Algorithm to Solve TripBased Dynamic Traffic Assignment
The development of efficient algorithms for the numerical computation of equilibria is a topic as important as the qualitative analysis of equilibria. The complexity of equilibrium problems, coupled with their increasing scale, is precluding their resolution via closedform analytics. Note that we will use a tripbased simulationbased algorithm to look for largescale user equilibrium. As mentioned before, the solution space is then discrete rather than continuous. The optimizer (algorithm) is embedded in the optimization framework of the DTA problem. First, we present the optimization framework, which comes from the literature. Then we focus on the optimizer part, wherein the path assignment is updated by heuristic and metaheuristic methods.
Finding an equilibrium through simulation (when no closedform analytical solution is available) typically involves a solution scheme that relies on an iterative procedure. In largescale application, this iterative procedure contains three costly steps: traffic simulation (network loading), timedependent shortest path calculation, and path flow optimization. In order to focus on the path flow optimization, we choose column generation approach.
Figure 1 presents the optimization framework based on the column generation approach. The framework separates the three costly steps in order to minimize the number of shortest path calculations which is very costly in largescale traffic networks. In other words, the idea of the column generation approach is based on the generation of the set of paths to which users may be assigned. Classical approach starts from an empty set and augments it in every iteration as required. The column generation approach contains two loops: outer loop and inner loop. The outer loop is responsible for path discovery, while the inner loops implement the path flow optimization for a given path set. Note that the classic approach executes both steps in a single loop. In largescale network problems, it is extremely costly to keep the data of all the possible paths between each OD pair, so, the column generation approach definitely reduces memory requirements [26]. Lu et al. [10] implement and examine the column generation approach in DTA context and show that it not only reduces memory requirements but also outperforms other algorithms in convergence with a designed swapping algorithm embedded in the inner loop.
This approach provides us with a framework to compare different path flow optimization algorithms to address which of them is more efficient for the multiclass multimodal DTA problem. Therefore, we are going to focus on the green modules in Figure 1. Before presenting the framework in detail, we need to define some solution quality indicators to measure the distance between each solution and optimal solution (UE path flow distribution).
The total gap as a quality indicator in time interval is
The quality of the solution is evaluated by two indicators. The first indicator that calculates the travel cost gap per user in seconds is
The total gap (Equation (10)) and the first quality indicator (Equation (11)) are calculated for each solution with the updated shortest path. The perfect user equilibrium state has zero ; thus, the best optimization algorithm obtains minimum . In other words, measures the distance between the algorithms’ solution and UE. The second indicator is the assignment violation. The violation indicator is defined in order to count users assigned on nonoptimal path(s). The following steps calculate the violation indicator:(1)Calculate the user violation: It is defined by considering the gap of each user. The following function defines the user violation: where is the experienced travel time of the user who travels for OD pair and is the shortest path of OD pair . In other words, if the gap between the user perceived travel time and shortest path travel time is bigger than 10% of the shortest path travel time, the user is in violation.(2)Compute the OD violation: The OD pair is in violation when there are more than 10% of the users on in violation. Particularly, function defines the OD violation: where is the set of users who travel for OD pair .(3)The violation indicator of network is the share of ODs which are in violation. Therefore, the second indicator for the quality of the solution is defined as follows:
As the first indicator evaluates the overall quality of the solution, the assignment violation evaluates the solution’s quality at the OD level. Based on and , we can measure the quality of the solution at each iteration. The indicator is used for checking the convergence. The steps of the column generation for simulationbased DTA are detailed in Algorithm 1.

The framework has two convergence tests to stop the calculation process in order to save computation time. In both steps, first, the iteration index is checked with the maximum number of loops ( and ) that we fix for each loop. Suppose that it is not reached ( or ); then we check the solution’s quality at the current loop, but there are three differences between the outerloop convergence test and innerloop convergence test. First, in the outer loop, we calculate the quality indicator (Equation (11)) based on the new shortest path(s), but the path set is fixed during the innerloop operations, and only the shortest path is updated in each iteration. Second, in the inner loop, we calculate the relative total gap (i.e., the difference of the total gap between two iterations) to check the improvement of the process and compare it with . is a threshold for checking the exploration rate of the optimization algorithm in the inner loop. But, in the outer loop, the goal is to monitor the quality of the solution. Therefore, is calculated and compared with the fixed value . If we have a solution with equal to or less than threshold , we terminate the algorithm. Third, in the outer loop, we have one more condition other than checking the solution’s quality. The condition is the shortest path checking. If we have no new shortest path per OD and good quality of the solution, then the process is converged. It means we will continue the process when the quality of the solution is not good until we no longer have a new shortest path.
The optimization algorithm is embedded in the inner loop of the process (Figure 1) which updates the assignment plan in each iteration. The heuristic and metaheuristic algorithms are defined in this step.
3.1. Heuristic Methods
Generally, the heuristic methods are pathswapping descent direction methods that decompose the problem into many (origindestination, departure time interval, and user class) subproblems. They solve each of the subproblems by adjusting timevarying OD flows between all nonshortest paths and the shortest path(s). We did a crosscomparison of heuristic algorithms in [17] with homogeneous users. Here we present the best algorithm for a large scale, which may have a good performance in multiclass and multimodal contexts.
3.1.1. Method of Successive Average with Ranking Method
The first method is an extension of the Method of Successive Average (MSA). The MSA algorithm is the classical method to solve the traffic assignment problem [27] and it is widely used in theory and application for DTA problems [28]. Equation (15) presents the swapping formula of MSA algorithm [29]:where is extracted from the auxiliary path assignments obtained by the allornothing procedure; that is, everyone is placed on the shortest path. Note that the step size in this method is . It means that, during iteration for this method, a fraction of users on each nonshortest path swaps to the one of the shortest path in the current iteration. MSA ranking adds a ranking method to choose the users from nonshortest paths to the shortest path(s). In other words, it ranks all users by their generalized cost and swaps the users with higher to the shortest path(s). The step size is the same as the MSA method equal to of all users on nonshortest paths.
3.1.2. GapBased Method
To enable the deployment of largescale (simulationbased) DTA models, [10] proposes a heuristic approach, which we called “gapbased” method in this study. This method updates the path flow in optimization iteration based on gap between the path generalized cost of user class and that OD shortest path generalized cost of user class :where denotes the outerloop iteration index, denotes the innerloop iteration index, and is the number of the shortest paths in iteration . This path assignment updating scheme implies a natural path flow adjustment mechanism: flows on the nonshortest paths are moved to the shortest path and the volume that moved out from a nonshortest path is proportional to the gap between the nonshortest path and shortest path over nonshortest path cost.
MSA ranking and gapbased method have two main drawbacks which are also the same for all heuristic methods that come from the flowbased traffic assignment: Using the step size guarantees the algorithm to converge but decreases the flexibility of the method to swap more users to the shortest path(s). Moreover, the final solution is obtained based on the initialization (allornothing) because the algorithm converges (by decreasing the step size) after a while even we use a not good initial assignment which is not really obvious in the large scale. Many of the methods in the literature (see, e.g., [8, 11, 12, 15, 16, 28, 30–32]) are using the average travel time of each path to calculate the indicators and to distinguish the shortest path(s) and other paths per OD. This means they do not consider users’ experienced travel time directly. Moreover, the optimizer is applied for each OD and swap users based on the different values of the paths’ average travel time. Note that users in the same departure time interval and on the same path can perceive different travel time values because of the network dynamic evolution.
These drawbacks show that we always look for the local optimum based on initialization and our step size. Here we present a probabilistic approach to help the heuristic optimization framework in order to find the global optimum for tripbased traffic assignment.
3.1.3. Probabilistic Method
The probabilistic method is a tripbased approach. This method has no step size for optimization [17]. For each user, the probability of swap to the shortest path(s) is calculated by the following formula:where is the experienced generalized cost of the user. Indeed, the pathswapping approach is applied to every user, not paths (number of users on path ). Afterwards, we use a random number uniformly distributed between 0 and 1 for taking a decision about swapping the user or not. If the random number is lower than the swap indicator, the user swaps to the shortest path; otherwise, the user keeps his/her path for the current iteration. The probabilistic process enables the algorithm to explore more the solution space without any step size. This heuristic method converges intuitively for a largescale network because the swap indicator decreases by swapping based on random numbers. Note that, at the end of each outer loop, the final assignment plan is the best assignment in terms of average gap between different innerloops solutions. It is obvious that, because of the swap indicator, the algorithm converges by random numbers to the solution with a minimum number of possible swaps, which indicates the small value of gap for each user. In other words, with every swapping decision, users try to optimize the in order to reduce the swap indicator (the possibility of swapping) for the next inner loop.
3.1.4. Hybrid Methods
Verbas et al. [33] applied an extension of the gapbased method to transit assignment. The method uses Equation (16) to determine the number of swapping users for each path and then the probabilistic method is applied to choose users from each path for swapping. This method is called “gapbased prob.” We examined this method in DTA test cases, and it fully dominates the gapbased method in different monoclass test cases [17]. Therefore, we consider this method in our numerical experiments.
The second hybrid method is similar to the probabilistic method with a different swap indicator:
This method is called “stepsize prob.” in this study. It has a good performance in largescale and congested traffic networks [20]. We developed a smart stepsize selection procedure for this method in which we specialize the step size for each OD pair and decrease the step size when the OD gap is not improved compared to the previous iteration [17]. The method is called “smart stepsize prob.” All the methods that are presented in this section will apply to multiclass multimodal DTA problem in order to evaluate their performance and also compare them with metaheuristic methods.
3.2. Metaheuristic Methods
As discussed in the previous section, in iteration , the classical (heuristic) algorithms try to improve the path flow distribution by swapping a fraction of users to the shortest path(s) from each nonshortest path. Then, one simulation is launched based on the updated path flow distribution. We have to wait until the simulation run is finished to know the new link travel times and adjust the path flow distribution to be tested in the next iteration accordingly. This is the serial process of the classical algorithms, which limits the solution space exploration and computational process. We developed two metaheuristic algorithms in [18] in order to overcome the problem of the serial approach. Two parallelized metaheuristic approaches are designed: the first derived from the simulated annealing framework and the second derived from that of the genetic algorithm.
3.2.1. Simulated Annealing Algorithm
The algorithm starts with an initial solution generated randomly; that is, users choose their path from the OD path set randomly. For solution , the total gap between the users’ travel cost and the shortest path travel cost (Equation (10)) is considered as the energy of the solution. A set of neighbor solutions are generated with respect to the current one based on the temperature of the current iteration. A neighbor solution is defined as a candidate to replace the current solution. The current phase of the iteration depends on the temperature of the process. Inspired by the physics of matter, this study distinguishes three different methods to generate a neighbor solution, gas, liquid, and solid; these methods represent the states of matter in nature. When the temperature is high ( where denotes the boiling temperature), the gas method is applied. When running the SA algorithm, by decreasing the temperature, the algorithm enters the liquid phase ( where denotes the melting temperature) and then the liquid method is applied. When the temperature is quite low , the solid method is applied. The algorithm is presented in Algorithm 2.


3.2.2. Genetic Algorithm
The genetic algorithm (GA) (first formalized by [34]) is inspired by natural selection and genetic variation. The genetic algorithm includes three main operators: selection, crossover, and mutation. The description of the operators is well presented in [35]. The GA framework is adapted to the traffic assignment problem in [18]. The algorithm is designed as a twolayer GA process to search solution space by changing the path flows in GA Inner and overcome the drawback of OD impact by considering different combinations of OD assignments in GAoperators. In other words, the heuristic algorithms plus a random method are applied in GA Inner, and the GA operators in one upper level generate different combinations of OD. GA Inner is a new layer in the optimization process compared to the classical GA. It has been added to account for correlations between OD assignments. Moreover, GA considers a set of solutions instead of a single solution at every iteration. The steps of the GA are detailed in Algorithm 3.
4. Numerical Experiments
The multiclass model is applied to a largescale multimodal network of Lyon 6th + Villeurbanne. This network that covers 25 and includes 1,883 nodes, 3,383 links, 94 origins, and 227 destinations is located in Lyon Metropolis and is described in Figure 2. All the signalized intersections in the real field have been implemented in the simulator with their actual signal timing. In this work, we use Symuvia (an open source simulator (https://github.com/Ifsttar/OpenSymuVia)) as a tripbased simulator for calculating the needed variables in the network. Symuvia is a microscopic simulator based on a Lagrangian discretization of the LWR model [38]. We get the simulation timestep to 1 second and collect the link traffic information (travel times) every 1 minute. Note that the dynamic simulator uses Equation (6) to calculate the generalized travel cost on a link at each time interval. There are buses (with/without a dedicated line), trams, and metro lines considering public transportations. In the simulation process, all these patterns, in addition to the private car and walking option, are simulated. For instance, for using a metro line, a user walks (walking time) to the station and then waits for the next train (waiting time). When the train arrives at the station, if the train has free capacity, the user can take the train (boarding time) and then take the next patterns until arriving at the destination point. All of these actions are considered for the travel time calculation. The users’ routes are determined by the DTA model and the rolling horizon technique [39], which determine the path flow distribution based on a prediction period of 30 minutes and an assignment period of 15 minutes. The multimodal timedependent shortest path algorithm is implemented based on the multiorigin Dijkstra algorithm [40] with considering hyperpaths based on [41] in order to find minimum cost path(s) in terms of travel time.
(a)
(b)
The network is loaded with 54,190 travelers of all ODs with given departure times in order to represent the 2.5 morning peak hours of the network between 6 : 30 and 9 : 00. The demand profile comes from the study of [42]. Note that, here, we fixed the departure time in order to focus on travelers’ path choices. The main reason is that this study aims to compare the performances of DTA solvers. By fixing the departure time, the input of the problem will be the same for all algorithms. Therefore, the solution space is fixed and given for the route (path and mode) choice problem, and we can compare the performances of different solution algorithms.
Figure 3(a) presents the demand profile of the numerical experiment. The figure presents the inflow rate of the trips to the traffic network. The total length of the Lyon 6th + Villeurbanne network is 496.73 kilometers. Note that each OD has a timedependent demand profile based on the real data of Lyon city [42]. The dynamic loading represents the saturation state of the network. To show a quick and synthetic overview of the network state, we plot the Macroscopic Fundamental Diagrams (MFD) of the network in Figure 3(b). This diagram represents the overall evolution of the traffic conditions in the network by plotting the total travel distance versus the total travel time [43]. Note that total travel distance is proportional to the mean network flow, while total travel time is equivalent to vehicle accumulation.
(a)
(b)
The MFD shows the network state for every 5 minutes between 7 : 00 and 9 : 00 during the full simulation. Before 7 : 00, the demand is low, and the network is almost empty. As such, we filtered out the first 30 minutes in order to present the critical part of the evolution of the network traffic conditions. For instance, point number 15 of the MFD (Figure 3(b)) shows the total travel distance and total travel time of all the travelers in the network between 8 : 10 and 8 : 15. First, the MFD curve increases from (0, 0). The traffic states remain undersaturated (point numbers 1 to 8) when demand is light, in this case from 6 : 30 to 7 : 40. Afterwards, travel production, which is equivalent to the total travel distance for a given period of time, stabilizes, while the accumulation (or total travel time) continues to increase (point numbers 9 to 13). This corresponds to the saturation level occurring from 7 : 40 to 8 : 05. The decrease in travel production and accumulation (point numbers 14 to 24) shows that the network starts to exit the saturation level at the end of the simulation period and slowly returns to the undersaturation level from 8 : 05 to 9 : 00.
The trips are divided into nine classes that are presented in Table 2. In this table, the user class’s fraction of the total demand is presented by “percentage of trips,” and, for each class, we use a specific VOT. The demand level and VOT values are extracted from the real data of the full Lyon networks [42]. The optimization goal is to find the user equilibrium based on the generalized cost functions.

The optimization process is executed for a maximum of 10 outer loops for five assignment periods. At each outerloop iteration, new shortest path(s) can be discovered; therefore, if at least a new shortest path is found per outerloop iteration, we have 11 paths per OD. In each outer loop, a maximum of 40 inner loops have been executed to find the best assignment pattern. Note that, in heuristic and SA algorithms, there is one more simulation per outer loop for keeping the best assignment pattern if the process is converged by the maximum number of the inner loops. This means that if the convergence does not happen during all inner loops, we take the best solution and simulate the network to prepare the result for the outerloop convergence test. The convergence tests of two loops are fixed for all algorithms. In the innerloop convergence test, the optimization will be stopped when we have less than 1% improvement . Improvement is related to the average gap between the current and previous inner loops. For the outer loop, we will stop the process when we do not find a new shortest path for all the ODs and also have a small value of gap or have no improvement in terms of total gap . In the next section, we will compare the performances of heuristic and metaheuristic algorithms in the multiclass multimodal DTA problem based on the numerical results.
5. Numerical Results and Discussion
The full optimization process is done for all methods. The aim of the experiments performed on the Lyon 6th + Villeurbanne network is to examine the algorithms’ convergence pattern of the solutions for a largescale network case. Thus, the UE is calculated for the considered network seven times with the heuristic and metaheuristic algorithms defined in Section 3. The indicator is used to evaluate the quality of the solution.
The results for the performance indicators of all the algorithms are presented in Table 3. GA leads to more simulations, and SA leads to only simulations compared to the 550 simulations required in heuristic algorithms. However, the computation times are significantly lower because of the parallel simulation framework. Moreover, the solutions obtained by the metaheuristic algorithms are significantly closer to the optimal UE than other methods. The metaheuristic algorithms dominate all heuristic methods in the final solution ; indeed, they are significantly better than the gapbased prob. algorithm, which is the best heuristic algorithm in this experiment. The GA and SA algorithms manage to reduce the UE of the gapbased prob. method by more than (Table 3). The violation indicator also shows that metaheuristic algorithms work much better than the heuristics methods (reduction of ).

Regarding the percentage of incomplete trips, also the metaheuristic algorithms dominate the others. The incomplete trips indicator denotes the share of trips, which is not finished by the end of the simulation in the final path flow distribution of each algorithm. A lower number of incomplete trips means a lower total travel time spent in the system over the simulation period. The SA algorithm finds the closest solution to UE (minimum ) in this study. Moreover, the final solution of SA has the best value for other quality indicators that are highlighted in bold in the SA row in Table 3.
Moreover, we measure the computation time (CT) of algorithms to evaluate the performance of them. The CTs of SA algorithm and GA are significantly better than heuristic methods. In particular, it is significantly better than the SSP method (bold values in GA row inside Table 3). Note also that the DTA process over the full simulation period requires considerable computational resources because of the network size. Therefore, the computational improvement obtained by the metaheuristic methods is high, 47 hours on average (two days) for the SA and GA, compared to the heuristic algorithms. The GA reduces the CT by one day and a half compared to SSP as the fastest heuristic method.
Regarding the computation resources, every simulation usually needs one central processing unit (core) of the processor. The heuristic algorithms, as mentioned before, are run in series, so they need only one core per iteration. The SA algorithm generates a maximum of three new path flow patterns in its liquid method per iteration. Therefore, we need a maximum of three cores for the simulation process. For the GA, according to population size, crossover rate, and the mutation rate, this value is maximum 30 new solutions from the GA in this study. Because all the experiments are conducted on a 64bit personal computer with 18 cores, the number of cores is limited to 18. If the number of new solutions in GA is bigger than 18, the algorithm executes in two successive phases: the first 18 simulations in the first phase and the remaining solutions in the second phase.
Figure 4 presents the convergence pattern for all seven algorithms in the last assignment period. The convergence pattern is drawn based on the average delay of all users. By the best heuristic methods, the final user average gap is around 30 sec, but this value for the metaheuristic method is around 10 sec. It means that the metaheuristic methods also obtain the better solutions in terms of delay. In Figure 4, the convergence pattern of all heuristic algorithms first decreases and then variates in a small range for the last iterations, where they seem to stick to the local optimum. Metaheuristic algorithms’ convergence pattern is always decreasing, which means that the algorithms are exploring and looking for the global optimum. Therefore, both metaheuristic algorithms produce not only better solutions in terms of optimality but also better (i.e., not increasing) convergence patterns than classical algorithms.
In order to consider the complexity of the problem on the performance of algorithms, a monoclass demand scenario is also designed. The volume of the demand and the traffic network are the same as those in the multiclass scenario. For monoclass users, we consider the weighted average VOT of the multiclass scenario (15.84 euros per hour), extracted from Table 2, and the same monetary cost for public transportation. The setting for the optimization framework, including the maximum inner and outerloops’ iterations, is fixed to the same values as the multiclass scenario. Table 4 presents the monoclass scenario results with the multiclass scenario. Same as the previous experiment, metaheuristic algorithms dominate heuristic algorithms. The results show that heuristic algorithms cannot guarantee the user equilibrium for the largescale network.

The results for and CT of the multiclass setting are worse than those of the monoclass setting because of the VOT of different classes that change the scale of and also require more CT to explore the solution space of the multiclass setting. However, regarding , it is even more beneficial to use metaheuristic methods over heuristic methods in a multiclass environment, as the improvements in are significantly better than those in the monoclass problem (Table 4). The same conclusion can be drawn for CT. The GA algorithm is the fastest algorithm in both scenarios. The performance comparison of solution algorithms in two demand scenarios shows that metaheuristic algorithms are more efficient when the problem’s complexity is increased.
6. Conclusion
This study develops the multiclass model by designing the hybrid generalized cost function VOT and cost function, including mode and economic classes. The model is extended for DTA, and variational inequality formulation is presented for the dynamic user equilibrium problem. This problem is also challenging from the computation point of view in the largescale application. We compared two different assignment frameworks (heuristic and metaheuristic) to execute the model on the largescale network in order to investigate which algorithm is more efficient. Seven optimization methods are discussed to calculate the user equilibrium. The model and assignment framework are applied to the Lyon 6th + Villeurbanne network, and the two optimization approaches are compared on the largescale network. The multiattribute multiclass model and optimization process have been proven to be successful when applied to our multimodal largescale network. The results show that the metaheuristic methods (SA and GA) with parallel computation work significantly better than the heuristic methods. The average user’s gap in generalized cost for metaheuristic algorithms is about 4.8 times in average better than the heuristic methods. Furthermore, the metaheuristic algorithms save a minimum of 36 hours in computation time. GA is the fastest algorithm, and SA provides the best solution for the network user equilibrium problem. Moreover, we can conclude that all heuristic methods do not guarantee the user equilibrium solution in largescale DTA problems.
As metaheuristic algorithms are new, there is room to improve them and also apply other kinds of metaheuristic algorithms. In addition to this, the authors propose the following topics for future works.
1. Consider departure time choice models. One possibility can be adding another layer of optimization for the departure time adjustment process. It can also be interesting to use the proposed parallel framework to optimize departure time and route choices simultaneously. 2. Extend the multiattribute, multiclass, tripbased model to activitybased models in order to consider travel with activity chains. 3. Compare and analyze the performance of the metaheuristic method in other largescale networks with different typologies and the other types of optimization methods. 4. Combine and investigate the metaheuristic methods with other kinds of traffic simulators and shortest path calculation methods.
Appendix
A. Literature Review on Travel Cost Definition in DTA Problems
In the literature, this classification process identifies eight categories of information about travelers:(1)Choice model: users are divided into different groups based on the behavioral model that they take into account to choose their route and mode(s). Usually, utility functions are formulated to express the choices of each class of users [44]. Random functions are used in some cases when there is less information about users’ choice model [45]. This classification considers the path (transportation mode included) choice model and users’ preferences.(2)Cost function: link cost or path cost can be different for each class of users; for example, the toll of the links can be different for each class. Also, path cost can be defined per user class [46]. In the dynamic case, the cost function is variable in each time interval per user class.(3)Economic attributes: users are categorized into different classes based on their income [47]. Considering the money budget for each class of users is another way to consider economic attributes.(4)Generalized cost function: the cost can be expressed in terms of travel time or monetary cost. The goal of the generalized cost function is to propose a function that aggregates time and monetary cost: Typically, the generalized cost function integrates travel time and travel cost . can take into account the classdependent value of time (VOT) [48]: where is VOT for class . Another approach to generalize the travel cost is to add a weighting parameter for cost summation [49]: where is a weight parameter and for class .(5)Knowledge level of the network: users can be classified by their ability to access realtime traffic information (ATIS). Each class of users attempts to optimize the travel cost based on their information from the network status [46, 50, 51].(6)Risk taking attitude: generally speaking, risk, uncertainty, and reliability are three concepts that can be defined for both sides: supplier (traffic network system) and demander (traveler) [52]. Reliability is defined as the perceived probability of favorable outcome. Risk is defined as the expected value of the unfavorable outcome. For traffic network as a supplier, uncertainty can be defined on demand. Travelers treat travel time variability as a risk in their travel choices because it introduces uncertainty for an ontime arrival at the destination [52]. Using travel time budget is an approach to classify the users in order to deal with the arrival time uncertainty. Therefore, the late arrival time will be penalized [53]. Moreover, user cost function can be developed by Value of Uncertainty (VOU) to consider the uncertainties of the travel [54].(7)Social class: classifying users into social classes does not make a difference in the equilibrium calculation. The advantage is when we have an elastic demand, and the goal is to have a better estimation of the demand in traffic assignment problems [55].(8)Vehicle class: in multimodal traffic network, the cost of travel depends not only on the traffic flow of one mode but also on the traffic loads of other modes. Vehicle classes are defined in order to calculate the travel time and cost based on mode attributes [56]. Table 5 presents a comparison between different studies in the literature and this study on multiclass traffic assignment.

Data Availability
No data were used for this study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 Research and Innovation Program (Grant agreement no. 646592 – MAGnUM project).
References
 J. G. Wardrop, “Road paper. some theoretical aspects of road traffic research,” Proceedings of the Institution of Civil Engineers, vol. 1, no. 3, pp. 325–362, 1952. View at: Publisher Site  Google Scholar
 R. Jayakrishnan, W. K. Tsai, and A. Chen, “A dynamic traffic assignment model with trafficflow relationships,” Transportation Research Part C: Emerging Technologies, vol. 3, no. 1, pp. 51–72, 1995. View at: Publisher Site  Google Scholar
 T. Schreiter, V. WageningenKessels, Y. Yuan, J. Van Lint, S. Hoogendoorn et al., “Fastlane: traffic flow modeling and multiclass dynamic traffic management,” in Proceedings of the TrailBeta Congress 2012: Mobility and LogisticsScience Meets Practice, pp. 3031, Rotterdam, The Netherlands, October 2012. View at: Google Scholar
 W. Song, K. Han, Y. Wang, T. Friesz, and E. Del Castillo, “Statistical metamodeling of dynamic network loading,” Transportation Research Procedia, vol. 23, pp. 263–282, 2017. View at: Publisher Site  Google Scholar
 X. Xu, A. Zockaie, H. S. Mahmassani et al., “Schedule consistency for daily activity chains in integrated activitybased dynamic multimodal network assignment,” Transportation Research Record: Journal of the Transportation Research Board, vol. 2664, no. 1, pp. 11–22, 2017. View at: Publisher Site  Google Scholar
 N. Alisoltani, M. Zargayouna, and L. Leclercq, “A sequential clustering method for the taxidispatching problem considering traffic dynamics,” IEEE Intelligent Transportation Systems Magazine, vol. 12, no. 4, pp. 169–181, 2020. View at: Publisher Site  Google Scholar
 R. Jayakrishnan and C. R. Rindt, “Distributed computing and simulation in a traffic research test bed,” ComputerAided Civil and Infrastructure Engineering, vol. 14, no. 6, pp. 429–443, nov 1999. View at: Publisher Site  Google Scholar
 T. L. Friesz and K. Han, “The mathematical foundations of dynamic user equilibrium,” Transportation Research Part B: Methodological, 2018. View at: Google Scholar
 S. Peeta and A. K. Ziliaskopoulos, “Foundations of dynamic traffic assignment: the past, the present and the future,” Networks and Spatial Economics, vol. 1, no. 34, pp. 233–265, 2001. View at: Publisher Site  Google Scholar
 C.C. Lu, H. S. Mahmassani, and X. Zhou, “Equivalent gap functionbased reformulation and solution algorithm for the dynamic user equilibrium problem,” Transportation Research Part B: Methodological, vol. 43, no. 3, pp. 345–364, 2009. View at: Publisher Site  Google Scholar
 G. Peque, T. Miyagi, and F. Kurauchi, “Adaptive learning algorithms for simulationbased dynamic traffic user equilibrium,” International Journal of Intelligent Transportation Systems Research, vol. 16, no. 3, pp. 215–226, 2018. View at: Publisher Site  Google Scholar
 H. Sbayti, C.C. Lu, and H. S. Mahmassani, “Efficient implementation of method of successive averages in simulationbased dynamic traffic assignment models for largescale network applications,” Transportation Research Record: Journal of the Transportation Research Board, vol. 2029, no. 1, pp. 22–30, 2007. View at: Publisher Site  Google Scholar
 I. Yang and R. Jayakrishnan, “Gradient projection method for simulationbased dynamic traffic assignment,” Transportation Research Record: Journal of the Transportation Research Board, vol. 2284, no. 1, pp. 70–80, 2012. View at: Publisher Site  Google Scholar
 M. Mahut, M. Florian, and N. Tremblay, “Comparison of assignment methods for simulationbased dynamicequilibrium traffic assignment,” in Proceeding of the Transportation Research Board 87th Annual Meeting (DVD), Washington, DC, USA, January 2008. View at: Google Scholar
 M. Florian, M. Mahut, and N. Tremblay, “Application of a simulationbased dynamic traffic assignment model,” European Journal of Operational Research, vol. 189, no. 3, pp. 1381–1392, 2008. View at: Publisher Site  Google Scholar
 H. X. Liu, X. He, and B. He, “Method of successive weighted averages (mswa) and selfregulated averaging schemes for solving stochastic user equilibrium problem,” Networks and Spatial Economics, vol. 9, no. 4, p. 485, 2009. View at: Publisher Site  Google Scholar
 M. Ameli, J. P. Lebacque, and L. Leclercq, “Cross‐comparison of convergence algorithms to solve trip‐based dynamic traffic assignment problems,” ComputerAided Civil and Infrastructure Engineering, vol. 35, no. 3, pp. 219–240, 2020. View at: Publisher Site  Google Scholar
 M. Ameli, J.P. Lebacque, and L. Leclercq, in Simulationbased Dynamic Traffic Assignment: MetaHeuristic Solution Methods with Parallel Computing, ComputerAided Civil and Infrastructure Engineering, Hoboken, NJ, USA, 2020.
 Y. Jiang, W. Y. Szeto, J. Long, and K. Han, “Multiclass dynamic traffic assignment with physical queues: intersectionmovementbased formulation and paradox,” Transportmetrica A: Transport Science, vol. 12, no. 10, pp. 878–908, 2016. View at: Publisher Site  Google Scholar
 H. Halat, A. Zockaie, H. S. Mahmassani, X. Xu, and O. Verbas, “Dynamic network equilibrium for daily activitytrip chains of heterogeneous travelers: application to largescale networks,” Transportation, vol. 43, no. 6, pp. 1041–1059, 2016. View at: Publisher Site  Google Scholar
 X. Xu, F. Fakhrmoosavi, A. Zockaie, and H. S. Mahmassani, “Estimating path travel costs for heterogeneous users on largescale networks: heuristic approach to integrated activitybased modeldynamic traffic assignment models,” Transportation Research Record: Journal of the Transportation Research Board, vol. 2667, no. 1, pp. 119–130, 2017. View at: Publisher Site  Google Scholar
 J. Wang, S. Peeta, and X. He, “Multiclass traffic assignment model for mixed traffic flow of humandriven vehicles and connected and autonomous vehicles,” Transportation Research Part B: Methodological, vol. 126, pp. 139–168, 2019. View at: Publisher Site  Google Scholar
 M. Ameli, Heuristic Methods for Calculating Dynamic Traffic Assignment, IFSTTAR and Université de Lyon, Luon, France, 2019, PhD Thesis.
 X. Zhan and S. V. Ukkusuri, “Multiclass, simultaneous route and departure time choice dynamic traffic assignment with an embedded spatial queuing model,” Transportmetrica B: Transport Dynamics, pp. 1–23, 2017. View at: Google Scholar
 M. Ameli, J. P. Lebacque, and L. Leclercq, “Daytoday multimodal dynamic traffic assignment: impacts of the learning process in case of nonunique solutions,” in Proceedings of the DTA 2018, 7th International Symposium On Dynamic Traffic Assignment, Hong Kong, June 2018. View at: Google Scholar
 M. W. Levin, M. Pool, T. Owens, N. R. Juri, and S. Travis Waller, “Improving the convergence of simulationbased dynamic traffic assignment methodologies,” Networks and Spatial Economics, vol. 15, no. 3, pp. 655–676, 2015. View at: Publisher Site  Google Scholar
 K. Nagel and G. Flötteröd, “Agentbased traffic assignment,” in The MultiAgent Transport Simulation MATSim, pp. 315–326, Ubiquity Press, London, UK, aug 2016. View at: Google Scholar
 R. Mounce and M. Carey, “On the convergence of the method of successive averages for calculating equilibrium in traffic networks,” Transportation Science, vol. 49, no. 3, pp. 535–542, aug 2015. View at: Publisher Site  Google Scholar
 H. Robbins and S. Monro, “A stochastic approximation method,” The Annals of Mathematical Statistics, vol. 22, pp. 400–407, 1951. View at: Google Scholar
 S. Boyles, J. Duthie, C. Melson, and T. Rambha, “Diverge models and dynamic traffic equilibria,” in Proceedings of the INFORMS Annual Meeting, Austin, TX, USA, November 2013. View at: Google Scholar
 C. O. Tong and S. C. Wong, “Heuristic algorithms for simulationbased dynamic traffic assignment,” Transportmetrica, vol. 6, no. 2, pp. 97–120, 2010. View at: Publisher Site  Google Scholar
 R. Mounce and M. Carey, “Route swapping in dynamic traffic networks,” Transportation Research Part B: Methodological, vol. 45, no. 1, pp. 102–111, 2011. View at: Publisher Site  Google Scholar
 Ö. Verbas, H. S. Mahmassani, and M. F. Hyland, “Gapbased transit assignment algorithm with vehicle capacity constraints: simulationbased implementation and largescale application,” Transportation Research Part B: Methodological, vol. 93, no. 116, 2016. View at: Publisher Site  Google Scholar
 J. H. Holland, Adaptation in Natural and Artificial Systems: An Introductory Analysis with Applications to Biology, Control, and Artificial Intelligence, MIT press, Cambridge, MA, USA, 1992.
 J. M. Johnson and Y. RahmatSamii, “Genetic algorithm optimization and its application to antenna design,” in Proceedings Of IEEE Antennas and Propagation Society International Symposium and URSI National Radio Science Meeting, IEEE, Seattle, WA, USA, June 1994. View at: Google Scholar
 M. Srinivas and L. M. Patnaik, “Adaptive probabilities of crossover and mutation in genetic algorithms,” IEEE Transactions on Systems, Man, and Cybernetics, vol. 24, no. 4, pp. 656–667, apr 1994. View at: Publisher Site  Google Scholar
 H. Maini, K. Mehrotra, C. Mohan, and S. Ranka, “Knowledgebased nonuniform crossover,” in Proceedings of the First IEEE Conference on Evolutionary Computation. IEEE World Congress on Computational Intelligence, Seattle, WA, USA, June 1994. View at: Google Scholar
 L. Leclercq, J. A. Laval, and E. Chevallier, “The Lagrangian coordinates and what it means for first order traffic flow models,” in Proceedings of the Transportation And Traffic Theory 2007. Papers Selected for Presentation at ISTTT17, London, UK, 2007. View at: Google Scholar
 H. S. Mahmassani, “Dynamic network traffic assignment and simulation methodology for advanced system management applications,” Networks and Spatial Economics, vol. 1, no. 3, pp. 267–292, 2001. View at: Google Scholar
 Y. Hawas, H. Mahmassani, R. Taylor, A. Ziliaskopoulos, G.L. Chang, and S. Peeta, Development of DynasmartX Software for RealTime Dynamic Traffic Assignment, University of Maryland, College Park, MD USA, 1998.
 İ. Ö. Verbas and H. S. Mahmassani, “Finding least cost hyperpaths in multimodal transit networks,” Transportation Research Record: Journal of the Transportation Research Board, vol. 2497, no. 1, pp. 95–105, 2015. View at: Publisher Site  Google Scholar
 J. Krug, A. Burianne, and L. Leclercq, “Reconstituting demand patterns of the city of Lyon by fusing multiple gis data sources,” University of Lyon, Lyon, UK, 2019, Tech. Rep. View at: Google Scholar
 N. Alisoltani, L. Leclercq, M. Zargayouna, and J. Krug, “Optimal fleet management for realtime ridesharing service considering network congestion,” in Proceedings of the 98th Annual Meeting of the Transportation Research Board, Washington, DC, USA, January 2019. View at: Google Scholar
 S.G. Li and Y.M. Su, “A multimode, multiclass dynamic network model with queues for advanced transportation information systems,” EJTIR, vol. 5, no. 1, pp. 55–80, 2005. View at: Google Scholar
 M. Patriksson, The Traffic Assignment Problem: Models and Methods, Courier Dover Publications, Mineola, NY, USA, 2015.
 H.J. Huang and Z.C. Li, “A multiclass, multicriteria logitbased traffic equilibrium assignment model under atis,” European Journal of Operational Research, vol. 176, no. 3, pp. 1464–1477, 2007. View at: Publisher Site  Google Scholar
 H. Yang and X. Zhang, “Multiclass network toll design problem with social and spatial equity constraints,” Journal of Transportation Engineering, vol. 128, no. 5, pp. 420–428, 2002. View at: Publisher Site  Google Scholar
 K. Zhang, H. S. Mahmassani, and C.C. Lu, “Dynamic pricing, heterogeneous users and perception error: probitbased bicriterion dynamic stochastic user equilibrium assignment,” Transportation Research Part C: Emerging Technologies, vol. 27, pp. 189–204, 2013. View at: Publisher Site  Google Scholar
 A. Nagurney and J. Dong, “A multiclass, multicriteria traffic network equilibrium model with elastic demand,” Transportation Research Part B: Methodological, vol. 36, no. 5, pp. 445–469, 2002. View at: Publisher Site  Google Scholar
 T.Y. Hu, C.C. Tong, T.Y. Liao, and L.W. Chen, “Dynamic route choice behaviour and simulationbased dynamic traffic assignment model for mixed traffic flows,” KSCE Journal of Civil Engineering, vol. 22, no. 2, pp. 813–822, 2018. View at: Publisher Site  Google Scholar
 B. Ran, D.H. Lee, and M. S.I. Shin, “New algorithm for a multiclass dynamic traffic assignment model,” Journal of Transportation Engineering, vol. 128, no. 4, pp. 323–335, 2002. View at: Publisher Site  Google Scholar
 J. G. Yu and R. Jayakrishnan, “Multiclass, multicriteria dynamic traffic assignment with pathdependent link cost and entropybased risk preference,” Transportation Research Record: Journal of the Transportation Research Board, vol. 2667, no. 1, pp. 108–118, 2017. View at: Publisher Site  Google Scholar
 X. Xu, A. Chen, Z. Zhou, and L. Cheng, “A multiclass meanexcess traffic equilibrium model with elastic demand,” Journal of Advanced Transportation, vol. 48, no. 3, pp. 203–222, 2014. View at: Publisher Site  Google Scholar
 M. A. AbdelAty, R. Kitamura, and P. P. Jovanis, “Using stated preference data for studying the effect of advanced traffic information on drivers' route choice,” Transportation Research Part C: Emerging Technologies, vol. 5, no. 1, pp. 39–50, 1997. View at: Publisher Site  Google Scholar
 M. Florian, J. H. Wu, and S. He, “A multiclass multimode variable demand network equilibrium model with hierarchical logit structures,” in Transportation And Network Analysis: Current Trends, pp. 119–133, Springer, Berlin, Germany, 2002. View at: Google Scholar
 W. H. K. Lam and H.J. Huang, “A combined trip distribution and assignment model for multiple user classes,” Transportation Research Part B: Methodological, vol. 26, no. 4, pp. 275–287, 1992. View at: Publisher Site  Google Scholar
 M. Meng, C. Shao, Y. D. Wong, and J. Zhang, “A multiclass, multimodal dynamic traffic assignment model with departure time,” Mathematical Problems in Engineering, vol. 2014, Article ID 812614, 9 pages, 2014. View at: Publisher Site  Google Scholar
 M. Florian and C. D. Morosan, “On uniqueness and proportionality in multiclass equilibrium assignment,” Transportation Research Part B: Methodological, vol. 70, pp. 173–185, 2014. View at: Publisher Site  Google Scholar
 S. Ryu, A. Chen, and K. Choi, “Solving the stochastic multiclass traffic assignment problem with asymmetric interactions, route overlapping, and vehicle restrictions,” Journal of Advanced Transportation, vol. 50, no. 2, pp. 255–270, 2016. View at: Publisher Site  Google Scholar
 N. Zheng and N. Geroliminis, “Modeling and optimization of multimodal urban networks with limited parking and dynamic pricing,” Transportation Research Part B: Methodological, vol. 83, pp. 36–58, 2016. View at: Publisher Site  Google Scholar
 G. R. de Andrade, Z. Chen, L. Elefteriadou, and Y. Yin, “Multiclass traffic assignment problem with flowdependent passenger car equivalent value of trucks,” Transportation Research Record: Journal of the Transportation Research Board, vol. 2667, no. 1, pp. 131–141, 2017. View at: Publisher Site  Google Scholar
 X. Pi, W. Ma, and Z. Qian, “A general formulation for multimodal dynamic traffic assignment considering multiclass vehicles, public transit and parking,” Transportation Research Part C: Emerging Technologies, vol. 104, pp. 369–389, 2019. View at: Publisher Site  Google Scholar
 R. Garcia and A. Marin, “Network equilibrium with combined modes: models and solution algorithms,” Transportation Research Part B: Methodological, vol. 39, no. 3, pp. 223–254, 2005. View at: Google Scholar
Copyright
Copyright © 2021 Mostafa Ameli 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.