Multi-Class Freeway Congestion and Emission Based on Robust Dynamic Multi-Objective Optimization

: In order to solve the problem of trafﬁc congestion and emission optimization of urban multi-class expressways, a robust dynamic nondominated sorting multi-objective genetic algorithm DFCM-RDNSGA-III based on density fuzzy c-means clustering method is proposed in this paper. Considering the three performance indicators of travel time, ramp queue and trafﬁc emissions, the ramp metering and variable speed limit control schemes of an expressway are optimized to improve the main road and ramp trafﬁc congestion, therefore achieving energy conservation and emission reduction. In the VISSIM simulation environment, a multi-on-ramp and multi-off-ramp road network is built to verify the performance of the algorithm. The results show that, compared with the existing algorithm NSGA-III, the DFCM-RDNSGA-III algorithm proposed in this paper can provide better ramp metering and variable speed limit control schemes in the process of road network peak formation and dissipation. In addition, the trafﬁc congestion of expressways can be improved and energy conservation as well as emission reduction can also be realized.


Introduction
The urban expressway system is the main framework system of the urban road system, and its operation service level directly affects the efficiency of urban road traffic. Taking Beijing, Shanghai and other big cities in China as an example, the problem of congestion on expressways is becoming increasingly serious [1]. At the same time, the on and off ramps of urban expressways are also closely related to the traffic efficiency of the main road. In the entrance and exit area of the ramps, the impact of the incoming area of the ramp on the mainline traffic can be reduced by active management and control of the incoming vehicles. Therefore, the traffic efficiency of the expressway can be improved, and the emission as well as fuel consumption can be reduced. Relevant studies show that these vehicle exhaust pollutants are one of the main factors causing urban haze [2].
Since the 1870s, abundant research results have been achieved by scholars in the field of expressway traffic control. Taking the traffic control strategy as an example, many countries already have quite scientific variable speed limit control technology. However, research in China is still lagging behind and is unable to take advantages of variable speed limit control technology. From the beginning of the 1990s, when variable speed limit control technology was put into practice, to now, despite the ramp control technology currently implemented in Shanghai, the ideal control effect still cannot be obtained. Therefore, this paper will focus on better traffic control by proposing a more practical algorithm.
Active management and control of multi-class expressway systems are studied in this paper. Through adopting effective control strategies, congestion and traffic efficiency can be improved, and vehicle exhaust emissions and fuel consumption can be reduced. In addition, considering the dynamic variation of the traffic process, it is particularly order models, the real traffic networks can be reproduced by the METANET model more accurately. For instance, the paper [7] extended a multi-class macroscopic traffic flow model, METANET, and proved that the model can lead to better fitting results because the heterogeneity of the traffic flow network structure was taken into account. Based on the above, this paper integrates the multi-class model proposed in paper [7] and a dynamic model proposed in paper [8]; therefore, a more comprehensive multi-class macroscopic traffic flow model is proposed.
Variable speed limits (VSL) and ramp metering (RM) are the most common strategies used in traffic control [9]. Based on fixed speed limits, VSL is a mainline control strategy where both frequent and incidental traffic congestion are considered. By using this strategy, dynamic speed limit suggestions can be provided to drivers in advance; therefore, the mainline flow upstream of the bottleneck can be well-managed to maximize the traffic throughput and prevent traffic congestion effectively [10]. In paper [7], a VSL-based cell transmission model was proposed to control the main roadway of a freeway, and the travel time was found to be reduced effectively. RM considers the traffic conditions and capacity on the downstream of the mainline. By using this strategy, the vehicle proportion of the ramp inflow can be determined without causing mainline traffic breakdown, therefore improving the network operation efficiency [11]. However, freeways still cannot be wellmanaged and controlled based on a single traffic control strategy due to the sophistication of expressway systems. Therefore, integrated control strategy is one of the present research directions [12]. The main idea is to regulate the expressway through coordinating more than one traffic control method. In paper [13], an integrated strategy was developed to coordinate multiple on-ramps and maximize the motorway bottleneck throughput. The limitations of the model in paper [13] include that only the throughput was taken into account, while the differences in vehicle classes were neglected. In conclusion, there are few studies on the integrated control strategy of VSL and RM on the aspect of multi-class expressways. Therefore, the research in this paper is meaningful.
In the process of optimizing practical problems, normally the original problem will be described as dynamic optimization problems (DOPs) when the environment changes dynamically [14]. In general, due to the real-time traffic flow variations on the entrance of the mainline and ramps, real-time oscillations often occur on the aspect of the traffic network states, such as density, flow, and speed, while dealing with the optimization problems of multi-class freeway congestion and emissions simultaneously. Therefore, freeway network control schemes need to be changed accordingly. The real-time variations of the demands, the network states and the control schemes were denoted in the form of parameters or performance indicators in DOPs. Moreover, the paper [15] described that historical environmental conditions and corresponding optimal solutions needed to be recorded and stored in the optimization process of DOPs. Therefore, obvious improvements in the evolution efficiency and search performance of the algorithm can be obtained through reusing these reserved solutions when detecting a similar dynamic environment. At present, robust dynamic optimization methods have been adopted by many researchers to solve DOPs, mainly including dynamic single-objective methods [16] and static multi-objective robust optimization methods [17]. However, research on robust dynamic multi-objective optimization methods is still limited. Hence, it is necessary to extend the research on robust dynamic multi-objective optimization algorithms extensively. While dealing with DOPs, the major challenge is to determine the distribution of Pareto fronts because of the irregular Pareto optimal front distribution nature. This is due to the fact that the robust dynamic solution sets, which constitute the Pareto optimal front, often change with the environment during the robust dynamic optimization process. Therefore, the research on the exploration of Pareto fronts will be reviewed below.
In paper [18], it is pointed out that clustering methods can help multi-objective optimization algorithms describe the Pareto front, detect the distribution of the solutions over the objective space, track the Pareto front continuously and guide search directions by the use of cluster centers, thereby leading to an increase in population diversity and convergence speed. The k-means method, proposed in paper [19], is one of the most widely used clustering algorithms which belongs to hard clustering. The principle of hard clustering is to divide data points into specific clusters, and the corresponding relationship between data points and clusters is unique. In paper [20], the k-means clustering algorithm was used to determine the irregular Pareto front shapes in a multi-objective optimization algorithm. By dividing the population into several cluster sets, the overall structure of the Pareto front can be discovered, which makes the population evolve in different cluster centers, thereby improving the diversity and convergence speed of the NSGA-II algorithm. However, the number of initial cluster centers in the k-means method needs to be determined by empirical methods. Moreover, random generation of initial cluster centers may lead to unstable clustering results. It is proposed in paper [21] that fuzzy clustering is one of the most widely used soft clustering algorithms which can accurately describe the mediation of the data samples. Among the fuzzy clustering methods, the fuzzy c-means clustering method, which combines traditional clustering algorithms with fuzzy theory, can provide more flexible clustering results [22]. Considering the large dimension of the objective spaces in the multi-objective optimization problems, in most cases, the solutions contained in the Pareto front are difficult to classify into a specific cluster center. Therefore, in this case, it is not reasonable to use hard clustering methods to divide the solutions into a specific cluster. Accordingly, it might be a better choice to use soft clustering methods instead to divide the samples into each cluster center according to a certain proportion. However, there are still two issues which need to be remedied. First, the fuzzy c-means clustering algorithm needs to initialize the cluster center, which will affect the clustering results. If the generated initial cluster centers are relatively poor, the algorithm may be trapped in local optima rather than obtain global optima. Second, when solving high-dimensional optimization problems, some marginal solutions may not be effectively covered by the randomly initialized cluster centers, which may reduce the diversity of solutions. In conclusion, the research on using clustering methods to determine the Pareto front is still not suitable. More in-depth research needs to be carried out due to the fact that the accuracy of the final Pareto front will be affected by the different selections of clustering algorithms and cluster centers.

Robust Dynamic Multi-Objective Optimization Problem in Freeway Traffic Congestion and Emission
Normally the optimization problem in the traffic field has the characteristics of high dimension, dynamic variation, and certain stability of the implemented optimization scheme. In this paper, the traffic congestion and emission optimization problems of multiclass expressways are described as robust dynamic multi-objective optimization problems. The traffic efficiency and environmental effects of expressways are measured based on the three indexes of travel times, traffic emissions and ramp queues. All the abbreviations can be found in Table A1 and relevant notations are concluded in Table A2.

Robust Dynamic Performance Indicators
When solving robust dynamic multi-objective optimization problems, the objective function can be designed to ensure the robustness of the solutions in a continuously changing environment. In this paper, the performance indicators are designed as follows in the robust dynamic multi-objective optimization problems of multi-class freeway congestion and emission.

Robust Dynamic Travel Times
The length of temporal windows is defined as T, and the number of the temporal windows is denoted as A. The robustness of travel times in continuous variable temporal windows with each length T and total number A can be demonstrated by this performance indicator. It can be denoted as the average travel time spent in the predesigned temporal windows, where two parts are included: the travel time spent in the mainstream and the waiting time spent at the ramps. Therefore, the objective function can be defined as: where c denotes the vehicle classes, i denotes the indexes for freeway sections and the freeway on-ramps are indicated by o. The time horizon is divided into K time steps, and k = 0, 1, 2, . . . , K indicates the temporal step counter, which relates to the time kT with sample time interval T (h). In addition, l o,i,c (k) is the queue length of class c in section i of on-ramp o in time interval (kT, (k + 1)T)). ρ i,c (k) is the traffic density of class c in section i in time interval (kT, (k + 1)T)). Additionally, L i and λ i represent the length of section i and the number of lanes in section i, respectively.

Robust Dynamic Traffic Emissions
To calculate robust dynamic traffic emissions, the length of temporal windows is denoted as T, and the number of them considered by decision maker is described as A analogously. This indicator is denoted as the average emissions generated by the total vehicles in A continuous variable temporal windows with the length T. This objective function can be calculated as: where A denotes the total number of temporal windows and k is the time step counter. The freeway stretch is divided into i sections, and the vehicle classes are indexed by c. In addition, the index y ∈ {CO, HC, NO x , FC} denotes the set of all the emission categories. J t y,i,c (k) indicates the emissions generated by class c travelling in section i in time interval (kT, (k + 1)T), in which the superscript t represents the normal travelling state. J s y,i,c (k) denotes the emissions generated by class c waiting in section i in time interval (kT, (k + 1)T), in which the superscript s denotes the stopping state. J y,on,o,c (k) represents the emissions generated by class c at on-ramp o in time interval (kT, (k + 1)T), in which the on-ramp is denoted as on and indexed by the notation o.

Robust Dynamic Ramp Queues
Similarly, the temporal window length is denoted as T, and the number of these temporal windows is described as A. Robust dynamic ramp queues can be denoted as the average number of vehicles waiting at the on-ramps in the continuous variable temporal windows. This performance indicator can be expressed as: where k denotes the time step counter, i denotes the on-ramp i of the freeway and the total number of temporal windows is described as A. The length of the temporal window or the sample time interval can be denoted as T (h). Additionally, l i,c (k) denotes the queue length of vehicle class c at on-ramp i in time interval (kT, (k + 1)T). d i,c (k) indicates the allowed traffic volume of class c entering mainstream at on-ramp i in time interval (kT, (k + 1)T), whereas r i,c (k) represents the actual traffic volume of class c entering the mainstream at on-ramp i in time interval (kT, (k + 1)T)).

Constraint
The detection of environmental variations, traffic safety requirements, and the range of variable values to meet actual demands are the three main constraints of this robust dynamic optimization problem of multi-class, multi-on-ramp, and multi-off-ramp freeway congestion and emission.

Environmental Variation Detection
It is judged whether the entering traffic volume at the on-ramps exceeds the variable range compared with the former time instant. This can also be expressed as the variations of r i,c (k), the allowed on-ramp traffic volume of vehicle class c entering the mainstream, which can be computed as: where r i,c (k) denotes the allowed on-ramp traffic volume of class c entering section i at time instant kT, which can be obtained as the ramp control flow of the last time in the optimization algorithm. r i,c (k) indicates the demanded on-ramp traffic volume of class c entering section i at time instant kT, and µ i,c (k) is the metering rate, which can also be described as the portion of the flow allowed to enter the mainstream under ramp metering.
In addition, three on-ramps are set in this paper; therefore i = 1, 2, 3, and σ denotes the threshold of the detection, which is set equal to 0.1 in this paper. It is also judged whether the entering traffic volume in the mainstream exceeds the variable range compared with the last time, which can be expressed as: where q i,c (k) denotes the traffic volume of class c entering section i in time interval (kT, (k + 1)T). Analogously, σ denotes the threshold of the detection. Variations are regarded to occur in the environment if the difference between two adjacent times exceeds the setting threshold. Additionally, the threshold σ is suitable to be set equal to 0.1 by multiple trials.

Traffic Safety Requirements
It is required that the variable range of the two variable speed limits should be constrained within 20 km/h, so that traffic incidents can be avoided though preventing the two continuous variations from being too large. It can be expressed as: where VSL can be described as the value of the variable speed limit and VSL c,i (k) denotes the variable speed limit of vehicle class c in section i at time kT.
It is also required that the difference between two adjacent sections should be constrained within 20 km/h, so that traffic breakdowns, such as rear-end collisions, can be avoided through preventing the deviation from being overlarge. It can be defined analogously as:

Range of Variables
The range of variables is required to meet the actual demands; therefore in this paper, the VSL of cars is set between 60 and 90 km/h and the VSL of trucks is set between 50 and 80 km/h, which can be expressed respectively as: VSL 1 ∈ (60, 90) (km/h) and VSL 2 ∈ (50, 80) (km/h). Additionally, both the ramp metering rates are set in the range (0, 1).

The Robust Dynamic Nondominated Sorting Multi-Objective Genetic Algorithm Based on Density Fuzzy C-Means Method DFCM-RDNSGA-III
The nondominated sorting multi-objective genetic algorithm NSGA-III proposed in paper [23] is mainly used to solve traditional static multi-objective optimization problems. This paper needs to solve the robust dynamic multi-objective optimization control problems of the freeway system. Therefore, the original NSGA-III algorithm needs to be improved to meet the requirements of providing robust solutions in a dynamic environment. At the same time, due to the difficulty in determining the Pareto fronts in the optimization process, a clustering algorithm is used in this paper to obtain the specific distribution of Pareto fronts. The clustering centers are used to describe the Pareto fronts and detect the distribution of the solutions in the target space. Through continuously tracking the Pareto fronts, the search is guided forward to the optimal Pareto fronts, and the population diversity and convergence speed can also be increased. Therefore, this paper improves the nondominated sorting multi-objective genetic algorithm NSGA-III in paper [23], and a robust dynamic nondominated sorting multi-objective genetic algorithm, RDNSGA-III, is proposed based on the concept of robust optimization and dynamic optimization.
Based on the original NSGA-III algorithm, the environment detection operator is introduced to store the environment information and the set of optimal solutions corresponding to the current environment. Additionally, the robust dynamic objective function is set to meet the needs of robust dynamic optimization. Secondly, in order to solve the problem of the difficulty of determining the Pareto fronts in the optimization process, a density fuzzy c-means clustering algorithm DFCM is proposed to determine the distribution of Pareto fronts. This algorithm can detect the distribution of solutions in the target space, track the Pareto front continuously, and guide the search direction of the algorithm. Finally, combining the robust dynamic optimization algorithm given in this section with density fuzzy c-means clustering algorithm DFCM, a robust dynamic nondominated sorting multi-objective genetic algorithm based on the density fuzzy c-means clustering algorithm DFCM-RDNSGA-III is proposed. The details of the proposed algorithm are shown in the following.

Robust Dynamic Nondominated Sorting Multi-Objective Genetic Algorithm RDNSGA-III
In the robust dynamic nondominated sorting multi-objective genetic algorithm RDNSGA-III proposed in this section, there are two main improvements: 1. In terms of dynamic characteristics, an environment detection operator is introduced based on the NSGA-III algorithm to test whether the external environment has changed. If environmental variations are detected, optimization will be restarted to obtain the optimal solutions which satisfy the current environment. Moreover, in the process of the algorithm operating, the environment information and the optimal solution sets corresponding to the current environment will be stored. It is ensured that repeated optimization under a similar environment can be avoided during the dynamic optimizing process; 2. In terms of robustness, the robust dynamic objective function is introduced to measure the robustness of the solutions in the nondominated sorting, i.e., the robust dynamic objective function is used to sort the solutions. By minimizing the average value of the objective function in multiple continuous temporal windows, the solution with stronger robustness is not only applicable to the current dynamic environment, but also applicable to multiple continuous dynamic environments. The robust dynamic objective function setting has been described in detail in Formula (1) to Formula (3) above.
The framework of the RDNSGA-III algorithm is described as follows: Step 1: Initialize the environment detection parameter, set the environmental detection counter t = 1, the maximum environmental detection times t max , and bestpop f inal = ∅ to store the current environment information and corresponding optimal solutions.
Step 2: If the environment detection counter t = 1, or if the environment detection counter t > 1 and the environment has changed, go to step 3. Otherwise, copy the current population P 1 , and go to step 12.
Step 3: Initialize algorithm parameters, including the maximum iterations Gen max , the number of the population pop, the current generation gen = 0, initialize the population P 0 = {x(1), . . . , x(pop)}, nondominated solution set S t = ∅, archive set D t = ∅, and reference point set Z s = ∅. Step 4: Recombine, crossover and mutate the population P t to generate offspring population Q t , Q t = Recombination(P t ) + Mutation(P t ).
Step 5: Combine archive set D t and offspring population Q t to generate the combined population R t , R t = D t ∪ Q t .
Step 8: Generate the next-generation population P t+1 . If the number of solutions in S t is exactly equal to N, i.e., |S t | = N, the next generation of parent population P t+1 is generated directly, P t+1 = S t , and t = t + 1, return to step 2; otherwise, structure P t+1 with F 1 , F 2 , . . . , F l−1 , i.e., P t+1 = ∪ l−1 j=1 F j . The remaining K solutions need to be selected from the layer F l according to the niche count, i.e., the number of solutions that the population P t+1 needs to select from F l is K = N − |P t+1 |.
Step 9: Generate reference points; the reference point set Z r = Normalize( f n , S t , Z s , Z r ).
Step 10: Associate the solutions in the solution set S t with the reference points.
Step 11: Calculate the niche count of each reference point j in the reference point set Z r , select elements according to the niche count to construct the population P t+1 , and select K solutions from F l to join the population P t+1 .
Step 12: Reserve the existing population, randomly take out half of the solutions and store them in the archive set D t+1 . Let gen = gen + 1; if gen < Gen max , go to step 4. if gen = Gen max , output the final solution set P f inal , and go to step 13.
Step 13: Save the current environment information and the final optimal solution P f inal to the optimal solution set bestpop, and add it to bestpop f inal (t). Let t = t + 1; if t < t max , go to step 2. If t = t max , output the final optimal solution set.

The Robust Dynamic Nondominated Sorting Multi-Objective Genetic Algorithm Based on Density Fuzzy C-Means Method DFCM-RDNSGA-III
In the paper [24], it is pointed out that the main challenge of solving the Pareto front distribution problem is the unknown distribution of Pareto optimal fronts in multi-objective optimization algorithms. The paper [25] pointed out that clustering algorithms can help multi-objective optimization algorithms detect the distribution of solutions in the objective space and guide the search to the Pareto fronts. However, the randomness of initial cluster generation and the difficulty of determining the number of initial cluster centers occur when using fuzzy c-means clustering algorithm. Therefore, a density-based clustering method (DBSCAN) [26] is adopted in this paper to conduct density clustering of the initial cluster centers when initializing the cluster centers.
The reasons to adopt this method can be explained from three aspects. Firstly, the spatial characteristics of data points can be reflected, and the rationality of cluster center distribution can be ensured when using density clustering, therefore improving the quality of the initial cluster centers of the fuzzy c-means clustering algorithm. Secondly, the number of cluster centers does not need to be set in advance because cluster centers can be generated automatically according to the specific density distribution of the population, which can solve the problem of setting the number of cluster centers by trial and error in fuzzy c-means clustering algorithm. In addition, through the fuzzy c-means clustering algorithm, the given initial cluster center set is continuously optimized, and the position of the cluster center is adjusted; therefore the shape information of more real Pareto fronts can be obtained.
Based on the RDNSGA-III algorithm in Section 4.1 and the DFCM algorithm proposed in paper [26], a robust dynamic nondominated sorting multi-objective genetic algorithm based on density fuzzy c-means method DFCM-RDNSGA-III is proposed in this section. The major improvement of this algorithm is that, in the following description of the DFCM-RDNSGA-III algorithm, specifically in Steps 9 and 10, the difficulty of determining Pareto fronts during the algorithm optimization process is solved based on density fuzzy c-means method DFCM. Specific distribution of the Pareto fronts can be acquired through adopting this clustering method. By describing the Pareto fronts through cluster centers, the distribution of the solutions over the objective space is detected, and the Pareto fronts can be tracked continuously [27]. Therefore, the search direction can be guided towards the Pareto fronts, leading to an increase in population diversity and convergence speed. Additionally, the details of the improvement are presented in sub-algorithm 1 and subalgorithm 2 in this section.
The flowchart of the proposed algorithm DFCM-NSGA-III is given in Figure 1. The framework of the proposed algorithm DFCM-NSGA-III will be described first, then subalgorithm 1 and sub-algorithm 2 used in the proposed algorithm will be discussed later.
tering algorithm, the given initial cluster center set is continuously optimized, and the position of the cluster center is adjusted; therefore the shape information of more real Pareto fronts can be obtained.
Based on the RDNSGA-III algorithm in Section 4.1 and the DFCM algorithm proposed in paper [26], a robust dynamic nondominated sorting multi-objective genetic algorithm based on density fuzzy c-means method DFCM-RDNSGA-III is proposed in this section. The major improvement of this algorithm is that, in the following description of the DFCM-RDNSGA-III algorithm, specifically in Steps 9 and 10, the difficulty of determining Pareto fronts during the algorithm optimization process is solved based on density fuzzy c-means method DFCM. Specific distribution of the Pareto fronts can be acquired through adopting this clustering method. By describing the Pareto fronts through cluster centers, the distribution of the solutions over the objective space is detected, and the Pareto fronts can be tracked continuously [27]. Therefore, the search direction can be guided towards the Pareto fronts, leading to an increase in population diversity and convergence speed. Additionally, the details of the improvement are presented in sub-algorithm 1 and sub-algorithm 2 in this section.
The flowchart of the proposed algorithm DFCM-NSGA-III is given in Figure 1. The framework of the proposed algorithm DFCM-NSGA-III will be described first, then subalgorithm 1 and sub-algorithm 2 used in the proposed algorithm will be discussed later. Specifically, the execution steps of the DFCM-RDNSGA-III algorithm are as follows: Step 1: Initialize the environment detection parameter, set the environmental detection counter t = 1, the maximum environmental detection times t max , and set bestpop f inal = ∅ to store the current environment information and corresponding optimal solutions.
Step 2: Initialize the optimization process of the algorithm, or if t > 1 and the environment has changed, go to Step 3. Otherwise, copy the current population P 1 , and go to Step 11.
Step 3: Initialize algorithm parameters, including the maximum iterations Gen max , the number of population pop, the current generation gen = 0, the population the nondominated solution set S t = ∅, the archive set D t = ∅, and the reference point set Z s = ∅.
Step 4: Recombine, crossover and mutate the population P t to generate the offspring population Q t , Q t = Recombination(P t ) + Mutation(P t ).
Step 5: Randomly select half of the solutions as the set D past from the current archive set {D t }. Combine the parent population P t and offspring population Q t to generate the combined population R t , R t = P t ∪ Q t ∪ D past .
Step 7: Generate nondominated solution set S t .
Step 8: Generate the next generation population P t+1 . If the number of solutions in S t is exactly equal to N, i.e., |S t | = N, the next generation of parent population P t+1 is generated directly, P t+1 = S t , and t = t + 1, go to Step 2; otherwise, construct P t+1 with F 1 , F 2 , . . . , F l−1 , i.e., P t+1 = ∪ l−1 j=1 F j . The remaining K solutions need to be selected from the layer F l according to the niche count, i.e., the number of solutions that the population P t+1 needs to be selected from F l is K = N − |P t+1 |.
Step 9: Generate reference points according to the density fuzzy c-means clustering algorithm DFCM and construct the reference point set Z r . The specific process is described in sub-algorithm 1.
Step 10: Associate the solutions in the solution set S t with the reference points. Calculate the niche count of each reference point j in the reference point set Z r , select elements according to the niche count to construct the population P t+1 , and select K solutions from F l as set {rest} to join the population P t+1 . The specific process is shown in sub-algorithm 2.
Step 11: Reserve the existing population, randomly take out half of the solutions and store them in the archive set D t+1 . Let gen = gen + 1. If gen < Gen max , return to step 4. If gen = Gen max , output the final solution set P f inal , and go to Step 12.
Step 12: Save the current environment information and the final optimal solution P f inal to the optimal solution set bestpop, and add it to bestpop f inal (t). Let t = t + 1; if t < t max , go to Step 2. If t = t max , output the final optimal solution set.
(1) Sub-algorithm 1: Calculate the cluster centers of F l and generate reference point set Z r according to the algorithm DFCM.
Step 1: Initialize the radius of neighborhood ε, and the minimum number of samples MinPts.
Step 2: Generate the initial cluster center set {Center 1 }. Calculate the initial membership degree µ ij of all solutions to the initial cluster center set {Center 1 } in the current situation, H denotes the number of cluster centers; m(m > 1) denotes the fuzzy allocation matrix index used to control the degree of fuzzy overlap. x i refers to the ith solution and c j denotes the jth cluster center.
Step 3: Calculate the initial adaptiveness function J of fuzzy c-means clustering algorithm, where G denotes the number of all solutions in the F l and H denotes the number of cluster centers.
Step 4: If the adaptiveness function J does not reach the acceptable extent or if the maximum number of iterations is satisfied, go to Step 5; otherwise, end the procedure and output the clustering result.
Step 5: Update the cluster centers to form a new cluster center set c j = Step 6: Update the membership degree µ ij and adaptiveness function J according to the new cluster centers and repeat Steps 3-5 until the output conditions are satisfied. Finally, output the final clustering results {Center} = {c 1 , c 2 , . . . , c H }, i.e., the reference point set Z r .
Step 7: Generate the solution set {xc i } of the corresponding solutions included in the ith cluster center c i according to the cluster center set {Center}. Define the set of all the cluster centers and the solutions included as {xc}.
(2) Sub-algorithm 2: Solution selection method based on the density fuzzy c-means clustering algorithm DFCM.
Sub-algorithm 2 provides the solution selection method based on the density fuzzy c-means clustering algorithm DFCM. The purpose is to select K solutions from F l , K = N − |P t+1 |, and join them to the set {rest} to construct the next generation population. The specific process is described as follows.
Step 1: Calculate the distance of all the solutions {x 1 , x 2 , . . . , x G } in F l and each reference point in reference point set Z r = {c 1 , c 2 , . . . , c H }. Distance = x i − c j 2 , x i denotes the ith solution in F l , while c j denotes the jth reference point in the reference point set Z r .
Step 2: Initialize the number of all the solutions of F l distributed by c j , nc j = 0. x i being dominated by c j can be denoted as x j i = argmin Distance x i , c j . If x i is dominated by c j , set nc j = nc j + 1.
Step 3: For j = 1 : H reference points, if nc j > 0, select the solution x j i min among all the solutions that has the closest distance to the reference point, and join it to set For DFCM-RDNSGA-III, calculating the cluster centers in step 9 requires O N 2 (subalgorithm 1). The density determination of the reference points requires O(N). In step 10 of DFCM-RDNSGA-III, the solution selection based on the density fuzzy c-means clustering method requires O MN 2 (sub-algorithm 2). Therefore, the computational complexity of DFCM-RDNSGA-III in the worst case is O(N 2 log M−2 N + MN 2 ), which is a little larger comparable to NSGA-III and RDNSGA-III.

Simulation
The simulation network and parameters are introduced in this section, then the traffic control strategy is described, and finally the experimental results are compared and analyzed.

Simulation Network
The multi-on-ramp and multi-off-ramp expressway network in paper [6] is adopted in this experiment. This freeway network has an 18 km main road, and there are three on-ramps (O 1 , O 2 , O 3 ) and three off-ramps (O 4 , O 5 , O 6 ). Two lanes are subsumed in the main road, which are equally divided into 18 sections. Each section is 1 km. Additionally, the on-ramp and off-ramp are all single-lane ramps, and the length of both on-ramps and off-ramps is 1 km. The specific network diagram is as shown in Figure 2. control strategy is described, and finally the experimental results are compared and analyzed.

Simulation Network
The multi-on-ramp and multi-off-ramp expressway network in paper [6] is adopted in this experiment. This freeway network has an 18 km main road, and there are three onramps ( 1 , 2 , 3 ) and three off-ramps ( 4 , 5 , 6 ). Two lanes are subsumed in the main road, which are equally divided into 18 sections. Each section is 1 km. Additionally, the on-ramp and off-ramp are all single-lane ramps, and the length of both on-ramps and offramps is 1 km. The specific network diagram is as shown in Figure 2:

Network Overview and Model Parameters
In order to simulate the peak hour curve of traffic inflow, the time span of the simulation time needs to be enough long to cover the whole process. Therefore, 2 h is the appropriate value. As for the sampling period, and control period, the experiment needs to collect enough data for calculation. On the one hand, it needs to quickly respond to changes in the current traffic flow. On the other hand, it cannot issue control commands very frequently, as this could result in traffic accidents. Based on the above considerations, the classical value is adopted, which can be seen in the paper [8]. To sum up, the simulation time is 2 h, the sampling period is 10 s and the control period is 5 min. The implemented control variables are adjusted according to the detection of environmental change and the robustness of the solution in the proposed DFCM-RDNSGA-III algorithm. The proportion of vehicles in the simulation road network is 80% cars and 20% trucks, assuming that the driver compliance rate for the control variables is 100%. The metering rate of the off-ramp is 10% of its adjacent mainline flow. The control parameter settings are shown in Table 1.

Network Overview and Model Parameters
In order to simulate the peak hour curve of traffic inflow, the time span of the simulation time needs to be enough long to cover the whole process. Therefore, 2 h is the appropriate value. As for the sampling period, and control period, the experiment needs to collect enough data for calculation. On the one hand, it needs to quickly respond to changes in the current traffic flow. On the other hand, it cannot issue control commands very frequently, as this could result in traffic accidents. Based on the above considerations, the classical value is adopted, which can be seen in the paper [8]. To sum up, the simulation time is 2 h, the sampling period T is 10 s and the control period is 5 min. The implemented control variables are adjusted according to the detection of environmental change and the robustness of the solution in the proposed DFCM-RDNSGA-III algorithm. The proportion of vehicles in the simulation road network is 80% cars and 20% trucks, assuming that the driver compliance rate for the control variables is 100%. The metering rate of the off-ramp is 10% of its adjacent mainline flow. The control parameter settings are shown in Table 1. Table 1. Multi-ramp expressway network simulation parameters.

Parameter Name Parameter Value
Simulation time 2 h Sampling period (T) 10 s Control period (KT) 5 min Table 2 is the actual value of the model simulation parameters. The multi-on-ramp and multi-off-ramp expressway network in paper [6] is adopted in this experiment, and all the actual values in Table 2 are set according to this expressway network. In Table 2, τ 1 denotes the car lag time parameters, η 1 denotes the car correction coefficient, v 1 denotes the car expectation constant, k 1 denotes the car ramp influence coefficient, τ 2 denotes the truck lag time parameters, η 2 denotes the truck correction coefficient, v 2 denotes the truck expectation constant, and k 2 denotes the truck ramp influence coefficient. In the freeway network, the real-time variations of the mainline car inflow and the truck inflow are shown in Figure 3. The real-time variations of the on-ramps (O 1 , O 2 , O 3 ) car and truck inflow are shown in Figure 4. inflow are shown in Figure 3. The real-time variations of the on-ramps ( 1 , 2 , 3 ) car and truck inflow are shown in Figure 4.

Algorithm Parameter Setting
Two existing multi-objective optimization algorithms, NSGA-III and RDNSGA-III, are used in this paper for comparison of results with the proposed DFCM-NSGA-III algorithm. The parameter setting of the three algorithms are shown in Table 3, where the vacancy value indicates that the type of parameter is not involved in the corresponding algorithm. Among them, the NSGA-III algorithm is a static multi-objective optimization algorithm that optimizes the variable speed limit and the ramp control strategy without involving the environmental detection parameters. Through several tests, it is shown that a satisfactory solution with better robustness can be obtained if the number of continuously changed temporal windows is set to be 3. To ensure safety in actual traffic driving, the control period implemented by each variable speed limit is 5 min. inflow are shown in Figure 3. The real-time variations of the on-ramps ( 1 , 2 , 3 ) car and truck inflow are shown in Figure 4.

Algorithm Parameter Setting
Two existing multi-objective optimization algorithms, NSGA-III and RDNSGA-III, are used in this paper for comparison of results with the proposed DFCM-NSGA-III algorithm. The parameter setting of the three algorithms are shown in Table 3, where the vacancy value indicates that the type of parameter is not involved in the corresponding algorithm. Among them, the NSGA-III algorithm is a static multi-objective optimization algorithm that optimizes the variable speed limit and the ramp control strategy without involving the environmental detection parameters. Through several tests, it is shown that a satisfactory solution with better robustness can be obtained if the number of continuously changed temporal windows is set to be 3. To ensure safety in actual traffic driving, the control period implemented by each variable speed limit is 5 min.

Algorithm Parameter Setting
Two existing multi-objective optimization algorithms, NSGA-III and RDNSGA-III, are used in this paper for comparison of results with the proposed DFCM-NSGA-III algorithm. The parameter setting of the three algorithms are shown in Table 3, where the vacancy value indicates that the type of parameter is not involved in the corresponding algorithm. Among them, the NSGA-III algorithm is a static multi-objective optimization algorithm that optimizes the variable speed limit and the ramp control strategy without involving the environmental detection parameters. Through several tests, it is shown that a satisfactory solution with better robustness can be obtained if the number of continuously changed temporal windows is set to be 3. To ensure safety in actual traffic driving, the control period implemented by each variable speed limit is 5 min.

Control Strategy
The variable speed limit of cars is set as (60, 90), the variable speed limit of trucks is set as (50, 80) and the ramp metering rate is set in the range (0, 1). A specific schematic diagram of control strategy implementation is shown in Figure 5, where VSL represents the section of variable speed limit control and ramp control indicates the on-ramp implementing the ramp control strategy.

Control Strategy
The variable speed limit of cars is set as (60, 90), the variable speed limit of trucks is set as (50, 80) and the ramp metering rate is set in the range (0, 1). A specific schematic diagram of control strategy implementation is shown in Figure 5, where VSL represents the section of variable speed limit control and ramp control indicates the on-ramp implementing the ramp control strategy. By adopting the robust dynamic multi-objective optimization algorithm, a set of optimal solutions will finally be obtained, therefore leading to the problem of choosing the optimal solution. Three performance indicators are considered in this paper: robust dynamic travel time, robust dynamic ramp queue and robust dynamic traffic emissions, which correspond to the patency degree of the main road, the patency degree of the ramps and the environmental influence. The optimal solution is selected by the weight method.

DFCM Clustering Analysis
Fuzzy c-means clustering is an unsupervised clustering method, so the following two indicators can be considered in judging its clustering effect: 1. The silhouette coefficient (SC) [28] is a way to evaluate the clustering effect, which can evaluate clustering results for different algorithms based on the same original data. The SC of clustering results is in the range (−1, 1). The greater the value, the closer the same samples are and the further the different samples are, therefore leading to a better clustering effect.
2. The Davies-Bouldin Index (DBI) [29], also known as the classification and suitability indicator, was proposed by David Davies and Donald Bouldin to evaluate the effect of the clustering algorithms. This index measures the mean value of maximum similarity of each cluster. A smaller DBI value represents that the clustering results are close within the cluster interior, and the different clusters are farther separated; that is, the smaller the distance is within the class, and the greater the distance is between the classes.
The two indicators are calculated as follows: where ( ) denotes the average distance from sample to other samples in the same cluster. The smaller it is, the more sample should be clustered into this cluster. It can be By adopting the robust dynamic multi-objective optimization algorithm, a set of optimal solutions will finally be obtained, therefore leading to the problem of choosing the optimal solution. Three performance indicators are considered in this paper: robust dynamic travel time, robust dynamic ramp queue and robust dynamic traffic emissions, which correspond to the patency degree of the main road, the patency degree of the ramps and the environmental influence. The optimal solution is selected by the weight method.

DFCM Clustering Analysis
Fuzzy c-means clustering is an unsupervised clustering method, so the following two indicators can be considered in judging its clustering effect: 1. The silhouette coefficient (SC) [28] is a way to evaluate the clustering effect, which can evaluate clustering results for different algorithms based on the same original data. The SC of clustering results is in the range (−1, 1). The greater the value, the closer the same samples are and the further the different samples are, therefore leading to a better clustering effect.
2. The Davies-Bouldin Index (DBI) [29], also known as the classification and suitability indicator, was proposed by David Davies and Donald Bouldin to evaluate the effect of the clustering algorithms. This index measures the mean value of maximum similarity of each cluster. A smaller DBI value represents that the clustering results are close within the cluster interior, and the different clusters are farther separated; that is, the smaller the distance is within the class, and the greater the distance is between the classes.
The two indicators are calculated as follows: where a(i) denotes the average distance from sample i to other samples in the same cluster. The smaller it is, the more sample i should be clustered into this cluster. It can be regarded as the dissimilarity of sample i within the cluster. Similarly, b(i) denotes the average distance from sample i to all samples in other cluster centers c j , which can be described as the dissimilarity between sample i and other cluster centers c j .
where N denotes the number of clusters. Variables i and j denote the cluster class i and cluster class j. S i and S j indicate the average distance from the data of cluster class i and cluster class j to the cluster center, respectively, which represent the dispersion degree of samples in cluster class i and cluster class j. w i − w j 2 denotes the distance between cluster class i and cluster class j.
In this paper, DFCM-RDNSGA-III, a robust dynamic nondominated sorting genetic algorithm based on density fuzzy c-means clustering is proposed. The whole generated population in the iterative process of the algorithm is used as the test set. The clustering effect of the DFCM algorithm is compared with that of the basic fuzzy c-means clustering algorithm. The number of samples in the data set is 100, the number of iterations is 100, and the data dimension is 3. The initial number of cluster centers is set according to the number of cluster centers generated by density clustering. The two algorithms are clustered 10 times, where EX is the mean of 10 times, DX is the variance of 10 times, DFCM denotes the density fuzzy c-means clustering algorithm, and FCM denotes the basic fuzzy c-means clustering algorithm. The comparison of indicators is shown in Table 4. It can be seen from Table 4 that the density fuzzy c-means clustering algorithm DFCM has a better clustering effect than the basic fuzzy c-means clustering algorithm FCM in terms of SC and DBI among the 10 clustering results. From the 10 running results, the average score of the DFCM algorithm in SC increased by 14%, which shows that the distance between different samples in DFCM is further and the clustering effect is better. Compared with FCM, DFCM improves 8% on DBI, which shows that DFCM has a closer distance within clusters and further distance between clusters; that is, DFCM generates a better clustering effect. Secondly, DFCM is relatively stable. Comparing the variance of 10 results, it can be found that the variance of DFCM in the SC and DBI is increased by 73% and 84%, respectively, which indicates that DFCM can be more stable in the final clustering results, and it is difficult to have a large difference in clustering results. Therefore, the density fuzzy c-means clustering algorithm DFCM can deal with the overall structure distribution of Pareto fronts more stably and accurately.

Results
Each experiment was run for 5 times and the results of the three indicators are the average results. The results of three objective functions, including travel time, ramp queue, and traffic emission of the three algorithms, NSGA-III, RDNSGA-III and DFCM-RDNSGA-III, are compared in Table 5. Among these three algorithms, a static optimization framework is adopted in algorithm NSGA-III to statically optimize the variable speed limit and metering rate in the whole simulation, while a dynamic optimization framework is adopted to dynamically optimize these variables according to the environmental changes in the RDNSGA-III algorithm and the DFCM-RDNSGA-III algorithm.
It can be concluded from Table 5 that: (1) The travel time and traffic emission indicators performed the worst on the algorithm NSGA-III with a static optimization framework. Compared to NSGA-III, the two indicators decreased by 14.64% and 16.72%, respectively, in the algorithm RDNSGA-III, as well as 19.04% and 21.64%, respectively, in the algorithm DFCM-RDNSGA-III. (2) Although NSGA-III performs poorly on the travel time and traffic emission indicators, it performs the best on the aspect of the ramp queue indicator. Compared to NSGA-III, the ramp queue indicator increased by 20.51% and 30.36%, respectively, in the algorithm RDNSGA-III and DFCM-RDNSGA-III. The above comparison results show that, under the control scheme of NSGA-III, which has a static optimization framework, the ramp-entering vehicles are not limited when the main road capacity is relatively tight. Therefore, the pressure on the main road will be aggravated due to the excessive ramp-entering, which increases the total travel time and traffic emissions but decreases the ramp queue. For the algorithms under robust dynamic optimization frameworks, including RDNSGA-III and DFCM-RDNSGA-III, although the ramp queue of the DFCM-RDNSGA-III algorithm is improved by 9% compared with the former, the travel time and traffic emission are both improved by more than 5%. The reason is that the inflow of part of the ramp vehicles to the main road is sacrificed in the DFCM-RDNSGA-III algorithm, contributing to the improved traffic efficiency of the main road and the reduced environmental pollution caused by gas emissions.
The VSL and RM control schemes during the peak forming period and peak subsiding period of the RDNSGA-III and DFCM-RDNSGA-III algorithms are shown in Tables 6 and 7, while the niche detection threshold is set as ε = 0.1. The environment detection is carried out every 5 min. When it is detected that the environment has not changed, i.e., the threshold before and after the change is set within the predefined range, the optimization does not need to be performed again and the current control scheme continues to be used. This also means that the situation where the control scheme has not changed represents that the optimization need not to be restarted. Therefore, the robustness of the control scheme of the RDNSGA-III and DFCM-RDNSGA-III algorithms can be compared under the same environment detection threshold through the change of control schemes. A static framework is adopted in the basic algorithm NSGA-III to optimize the VSL and RM, so the change of the control scheme is not involved in it.
The variation in traffic control schemes of the RDNSGA-III algorithm and the DFCM-RDNSGA-III algorithm according to the change of the environment in the processes of peak hour and dissipation hour are shown in Tables 6 and 7. It can be seen that the average ramp metering rate of the RDNSGA-III algorithm at the three on-ramps is 68% within 0 min to 5 min of the peak formation process, and the ramp metering rate of the DFCM-RDNSGA-III algorithm is 48%. In the early stages of peak formation, the DFCM-RDNSGA-III algorithm begins to limit the ramp entry, so as to ensure the patency of the main road in the process of peak formation. After that, the ramp metering rate of the RDNSGA-III algorithm is still maintained at more than 60% within 5-15 min. At 15-20 min, it starts to limit the ramp entry, and the average ramp metering rate is set to 56%, while the ramp metering rate of the DFCM-RDNSGA-III algorithm is always set at about 50% during this period. Therefore, in the early peak formation process, the inflow of the DFCM-RDNSGA-III algorithm is maintained within the scope of environmental detection by limiting the ramp entry. While the ramp inflow in the RDNSGA-III algorithm is not restricted well in the early stage, the flow changes at the last and next time exceed the detection threshold. Therefore, within 20 min of peak formation, the control strategy is continuously adjusted before it starts to maintain stability. Similarly, RDNSGA-III has a certain conservative setting of VSL within 60 min to 75 min of peak dissipation. At the beginning of the peak dissipation process, the VSL is still controlled below 70 km/h. However, the DFCM-RDNSGA-III algorithm starts to maintain the VSL above 70 km/h after 60 min to 65 min of early peak dissipation process, therefore ensuring the rapid outflow of the main road in the process of peak dissipation. In general, the DFCM-RDNSGA-III algorithm has better control optimization stability than the RDNSGA-III algorithm in the process of peak formation. The control strategy can be better dealt with in the process of peak formation, and solutions with better robustness can be generated under the same environmental sensitivity.  Figures 6 and 7 show the heat map of car speed and car flow in each road section of the three algorithms, in which the ordinate denotes the main road section number link 1 to link 18, and the abscissa denotes the time interval. In Figure 6, the speed heat map, it can be noticed through comparison that in the process of peak formation, the speed of each road section in DFCM-RDNSGA-III can be maintained at a relatively high level compared with the other two algorithms, which can better help control the ramp metering rate. Therefore, the traffic efficiency and vehicle speed of the main road can be improved. As can be seen from the flow heat map in Figure 7, the overall process of traffic flow is smoother from link1 to link18 in the DFCM-RDNSGA-IIII algorithm. The flow variation in each road section is within 500 veh, and the flow difference between each road section is small. However, the traffic flow of the NSGA-III algorithm and the RDNSGA-III algorithm have quite obvious and abrupt variations in the ramp merging section, especially at the ramp merging part of link 10. The traffic flow of NSGA-III algorithm suddenly changes from 2000 veh to 3000 veh, and the traffic flow of RDNSGA-III algorithm also increases by more than 500 veh. In the ramp merging section, the traffic flow surge indicates that in the on-ramp area, there is too much traffic flow at the on-ramp entrance, which leads to traffic congestion due to a significant variation in the traffic flow of the front and rear sections.  Figures 8 and 9 show the heat map of truck speed and truck flow in each road section of the three algorithms. As can be seen in Figure 8, in the second half of the road, that is, after link 10, the overall speed of the trucks in the DFCM-RDNSGA-III algorithm is slightly better than the other two algorithms. The main reason is that the DFCM-RDNSGA-III algorithm performs relatively better in the control of the last two ramps, which ensures the speed of the trucks on the main road. As can be seen in Figure 9, the performance of the heat map of the truck flow on the aspect of the three algorithms is basically the same as that of the car flow. The flow of the NSGA-III algorithm changes greatly before and after link 10. Compared with the RDNSGA-III algorithm, the DFCM-RDNSGA-III algorithm is smoother in the flow transmission and has a gradual increase process. In general, the control effect of the DFCM-RDNSGA-III algorithm is better than the other two algorithms.
rithm have quite obvious and abrupt variations in the ramp merging section, especially at the ramp merging part of link 10. The traffic flow of NSGA-III algorithm suddenly changes from 2000 veh to 3000 veh, and the traffic flow of RDNSGA-III algorithm also increases by more than 500 veh. In the ramp merging section, the traffic flow surge indicates that in the on-ramp area, there is too much traffic flow at the on-ramp entrance, which leads to traffic congestion due to a significant variation in the traffic flow of the front and rear sections.   Figures 8 and 9 show the heat map of truck speed and truck flow in each road section of the three algorithms. As can be seen in Figure 8, in the second half of the road, that is, after link 10, the overall speed of the trucks in the DFCM-RDNSGA-III algorithm is slightly better than the other two algorithms. The main reason is that the DFCM-RDNSGA-III algorithm performs relatively better in the control of the last two ramps, which ensures the speed of the trucks on the main road. As can be seen in Figure 9, the rithm have quite obvious and abrupt variations in the ramp merging section, especially at the ramp merging part of link 10. The traffic flow of NSGA-III algorithm suddenly changes from 2000 veh to 3000 veh, and the traffic flow of RDNSGA-III algorithm also increases by more than 500 veh. In the ramp merging section, the traffic flow surge indicates that in the on-ramp area, there is too much traffic flow at the on-ramp entrance, which leads to traffic congestion due to a significant variation in the traffic flow of the front and rear sections.   Figures 8 and 9 show the heat map of truck speed and truck flow in each road section of the three algorithms. As can be seen in Figure 8, in the second half of the road, that is, after link 10, the overall speed of the trucks in the DFCM-RDNSGA-III algorithm is slightly better than the other two algorithms. The main reason is that the DFCM-RDNSGA-III algorithm performs relatively better in the control of the last two ramps, which ensures the speed of the trucks on the main road. As can be seen in Figure 9, the performance of the heat map of the truck flow on the aspect of the three algorithms is basically the same as that of the car flow. The flow of the NSGA-III algorithm changes greatly before and after link 10. Compared with the RDNSGA-III algorithm, the DFCM-RDNSGA-III algorithm is smoother in the flow transmission and has a gradual increase process. In general, the control effect of the DFCM-RDNSGA-III algorithm is better than the other two algorithms.

Conclusions
In this paper, the robust dynamic traffic multi-objective optimization problem of a multi-class expressway is studied. By considering three traffic indicators-travel time, ramp queue and traffic emissions-at the same time, the implementation of variable speed limits and ramp metering strategies on an urban expressway is explored in this paper. Based on the multi-class macroscopic traffic flow model multi-class METANET and the multi-class emission and fuel consumption model multi-class VT-macro, a multi-on-ramp and multi-off-ramp expressway model is proposed. The multi-on-ramp and multi-offramp expressway system optimization problem is described as a robust dynamic multiobjective optimization problem. The DFCM-RDNSGA-III algorithm is proposed to solve

Conclusions
In this paper, the robust dynamic traffic multi-objective optimization problem of a multi-class expressway is studied. By considering three traffic indicators-travel time, ramp queue and traffic emissions-at the same time, the implementation of variable speed limits and ramp metering strategies on an urban expressway is explored in this paper. Based on the multi-class macroscopic traffic flow model multi-class METANET and the multi-class emission and fuel consumption model multi-class VT-macro, a multi-on-ramp and multi-off-ramp expressway model is proposed. The multi-on-ramp and multi-offramp expressway system optimization problem is described as a robust dynamic multiobjective optimization problem. The DFCM-RDNSGA-III algorithm is proposed to solve the variable speed limit and ramp metering scheme under this problem, and the algorithm effect is verified based on VISSIM simulation. In general, compared with NSGA-III and RDNSGA-III, DFCM-RDNSGA-III, a robust dynamic genetic algorithm based on density fuzzy c-means clustering method, can better deal with the robust dynamic optimization problem which considers environmental factors and traffic congestion.
In future research, the robust dynamic multi-objective optimization algorithm needs to be improved in the following aspects: the high evolution pressure when the dimension is high, the insufficient factors considered in the evolution process to adapt to the complexity of high-dimensional problems, the irregular Pareto front exploration, and the balance between dynamic capabilities and robustness. Moreover, the existing traffic flow model does not match the actual traffic operation completely. Factors affecting the traffic condition must be studied from more aspects, and a traffic flow model better matching the actual road network must be designed.

Conflicts of Interest:
We declare there are no conflict of interest regarding the publication of this paper. We have no financial and personal relationships with other people or organizations that could inappropriately influence our work.