Variable Cell Transmission Model for Mixed Traffic Flow with Connected Automated Vehicles and Human-Driven Vehicles

,


Introduction
In recent years, the rapid development of automatic driving and communication technology enabled the large-scale application of CAVs. e intelligent transportation system with CAVs has attracted the increasing attention of scholars [1]. Nevertheless, the popularity of CAVs takes a specific time; there will be a mixed traffic flow composed of CAVs and HDVs for a long time in the future. erefore, it is of great significance to study the traffic characteristics and propagation mechanism of mixed traffic flow. In real-time management and control of traffic flow, many studies have been applied modeling and simulation of mixed traffic flow with CAVs. ese studies can be divided into two categories: one focuses on theoretical analysis [2][3][4] and the other on simulation analysis [5][6][7]. is work mainly focuses on the second part, i.e., simulation analysis.
Against such a background, to gain the basic characteristics of general traffic flow, the homogeneous HDVs traffic flow environment is studied at the macro level at first. Daganzo [8,9] discretized LWR (Lighthill-Whitham Richards) model to obtain the Cell Transmission Model (CTM). Compared with higher order traffic flow models [10], CTM is a first-order model, which is widely applied in traffic simulation [11], traffic distribution [12], and traffic signal control [13] because of the low requirement of calculation. en, many scholars have focused on making corresponding improvements on the CTM for different study scenarios, such as city roads, expressways, intersections, and air transportation. To predict the traffic flow accurately, Tiaprasert et al. [14] and Szeto et al. [15] proposed a multiclass CTM and a real-time traffic control method, respectively. Carlos and Ferrara [16] studied the optimal speed control design of loop traffic based on the variable-length CTM (VLM) to optimize the travel time. Zhang et al. [17] studied the impact of traffic demand of the target road network on the overall progress of congestion evacuation based on the CTM and the macroscopic fundamental diagram (MFD). Results showed that the model could improve the computational efficiency of the dynamic traffic assignment model. Shirke et al. [18] proposed an arterial CTM (ACTM) to capture the congestion and overflow of turning lanes in the arterial traffic network. Compared with the classical CTM, ACTM can improve the accuracy of queue estimation by 20%∼80%. Dong et al. [19] developed an improved CTM to predict the duration of traffic congestion and estimate the spatial diffusion law of congestion. Moreover, an asymmetric CTM based on the upper limit of ramp flow rate and queue length was proposed to achieve the optimal control of on-ramp flow by Gomes and Horowitz [20]. For unsignalized intersections, Flötteröd and Rohde [21] studied an incremental node model (INM) to solve the problem of priority allocation based on the incremental transfer principle of the CTM. For signalized intersections, Srivastava [22] and Gao et al. [23] estimated the capacity of intersections based on the CTM with a new demand function and the Enhanced-CTM, respectively. Furthermore, according to the CTM and the fundamental diagram model, Lo [24] proposed a dynamic traffic control equation to predict traffic conditions and demand patterns. In traffic signal optimization, Lin and Wang [25] developed an improved mixed-integer linear programming model to obtain the optimal signal sequence based on CTM, which is used to capture the signal priority of vehicles at intersections.
However, as previously mentioned, the studies are based on HDVs to expand the CTM, without considering the impact of CAVs. At present, wireless communication and automatic driving technology promote the development of CAVs. e application of CAVs will significantly improve traffic efficiency and safety [26][27][28]. e large-scale application needs the support of related infrastructure such as intelligent roads. However, the upgrading of smart road infrastructure is a long process. Consequently, the transition period from homogeneous HDVs to homogeneous CAVs is long [29]. To analyze the characteristics of the mixed traffic flow composed of CAVs and HDVs has attracted a great deal of attention. According to the traffic flow theory, Levin and Boyles [30] established a multilane CTM to study the optimization algorithm for lane direction decision-making of CAVs. e results showed that the proposed heuristic algorithm could reduce the travel time of the traffic system, but they only focused on the traffic flow characteristics in the homogeneous CAVs environment. Based on the fundamental diagram of mixed traffic flow, Qin and Wang [31] proposed a CTM for mixed traffic flow with different MPR of CAVs. Results showed that the degree of queue length and dissipation time decreased with the MPR of CAVs. Whereas they did not consider the limitation of the CTM (e.g., the equal).
To sum up, the current researches inspire us to focus on the following three gaps: (1) most of the existing research uses microscopic models to capture the driving behavior of vehicles and describe the propagation mechanism of traffic flow. However, there are many parameters of the microscopic model that must be calibrated by collecting many data. In addition, compared with the meso and macro models, the micro model has low efficiency with much calculation. erefore, using the macro CTM will reduce the complexity of the simulation model and improve the simulation efficiency. (2) Most CTM models only study the characteristics of homogeneous traffic flow (e.g., HDVs or CAVs). With the development of CAVs, there will be a mixed traffic flow composed of CAVs and HDVs in the future. (3) In some studies of mixed traffic flow, they did not break through the limitation of equal cell length of the CTM. What is more, its basic assumption is that the traffic flow obeys the triangle fundamental diagram, and the hysteresis phenomenon is not considered. It is impossible to accurately describe the relationship between traffic volume and density under congestion, and the results of the study are far from the actual situation. Moreover, based on the three-phase traffic flow theory, the volume-density-speed of traffic flow is not an ideal triangular relationship.
ere is a substable region, which presents a converse 'λ' shape. When traffic density is greater than the critical density, there will be 'ghost jam'. erefore, it is more realistic to propose a VCTM to describe the characteristics of mixed traffic flow considering the hysteresis phenomenon.
To fill in these gaps, this study proposes a VCTM for mixed traffic flow with CAVs and HDVs. e main contributions of this work can be summarized as follows: (1) e VCTM is proposed based on "λ" shape fundamental diagram of mixed traffic flow to consider the cell density and cell length simultaneously (2) A microscopic simulation is designed to verify the effectiveness of the VCTM in SUMO based on IDM car-following model (3) Some congestion indexes, such as travel speed, delay, and congestion scale, are selected to discuss the characteristics of mixed traffic flow based on the VCTM e remainder of this paper is organized as follows. Section 2 proposes the VCTM based on the classical CTM and the CTM for mixed traffic flow. e validity of the VCTM is verified by SUMO simulation, and the applicability of CTM and VCTM is presented in Section 3. Section 4 analyzes the relationship between the time-space change of traffic flow and the MPR of CAVs. Finally, the conclusions and future works are summarized in Section 5.
where q is the traffic flow, q max is the maximum traffic flow, w is the reverse wave speed, v f is the free-flow speed, k is the traffic density, k j is the congestion density, and k 1 and k 2 are the minimum and maximum traffic density, respectively. Figure 2, the number of vehicles in cell i at time t can be expressed as equation (2), and f i (t) is the number of vehicles entering cell i at time t, which can be calculated by

Basic Link. From
where n i (t) is the number of vehicles in cell i at time t, Q i (t) is the maximum flow of cell i at time t, andN i (t) is the maximum number of vehicles that cell i can hold at time t. ese two capability parameters can describe the number of vehicles that can be sent and received by cell i at time t, which is defined in equations (4) and (5). en, (2) can be simplified to equation (6).
erefore, the traffic flow characteristics of the basic link of CTM can be described by equations (2)-(6).

Merging Link.
After satisfying the transport constraints between cells in Figure 3, considering the receiving vehicles of cell C limit the transmission vehicles of cells A and B, the diverging model can be divided into two cases.
(2) When S A (t) + S B (t) > R C (t), considering the priority of merging, the merging coefficients α A (t) and β A (t) need to be introduced here.

Diverging Link.
e receiving vehicles of cells F and G limit the transmission vehicles to two cells when cell E diverges, as shown in Figure 4. e diverging coefficient β F (t) and β G (t) is introduced to consider a diverging priority.
erefore, equations (2), (6), and (9) can be used to describe the traffic flow characteristics of the CTM diverging link.

Car-Following Modes of Mixed Traffic Flow.
e mixed traffic flow is composed of HDVs and CAVs. ere are three car-following modes based on vehicle-to-vehicle (V2V) communication conditions [32]. As shown in Figure 5(a), if the current and the following vehicle are CAVs, it can be described as a cooperative adaptive cruise control (CACC) mode. As shown in Figure 5  upstream Journal of Advanced Transportation HDV, the following car is a CAV. Because the HDV lacks network communication equipment, the CAV cannot drive cooperatively with it, resulting in CACC vehicles degradation in the mixed traffic flow [28]. Moreover, it can be described as adaptive cruise control (ACC) mode. In addition, if the current vehicle is an HDV or CAV, and the following vehicle is HDV, as shown in Figure 5(c).
In Figure 5(c), it can be described as a human-driven mode (HDV). erefore, if the platoon intensity is not considered [33], that is, the random distribution of CAVs and HDVs. e MPR of CAVs in mixed traffic flow is p, the proportion of different car-following modes in mixed traffic flow can be obtained by equation (15).where p CACC , p ACC , and p HDV are the proportion of CACC, ACC, and HDV carfollowing mode, respectively.

e CTM for Mixed Traffic Flow.
In order to facilitate the derivation of the fundamental diagram, the lane-changing behavior of vehicles is not considered. e relationship between traffic volume, density, and speed in mixed traffic flow can be deduced based on related research [15,17,34].
where M is the type of car-following modes; in this study, the M is 3, which means CACC, ACC, and HDV. p m is the proportion of the mth car-following mode of mixed traffic flow, v is the equilibrium speed, T m is the headway of vehicles in mth car-following mode, and d m is the minimum headway of vehicles in mth car-following mode. Moreover, when the velocity is free-flow speed, the q max corresponding critical density k c can be deduced in equations (14) and (15) based on equation (11).
e maximum traffic volume is where the speed is 0, the congestion density is erefore, the fundamental diagram of mixed traffic flow can be obtained by equations (14)- (16). From Figure 6, the reverse wave speed and speed of the mixed traffic flow can be calculated by en, equation (19) can be obtained by taking equations (15)-(17) into equation (18).

Analysis of Fundamental Diagram of Mixed Traffic
Flow. e volume, density, and speed of mixed traffic flow are related to the MPR of CAVs based on equation (19). In particular, refer to [3,35], the parameters setting of different car-following modes are shown in Table 1. In addition, the MPR are set as 0, 0.2, 0.4, 0.6, 0.8, and 1, respectively. en, the fundamental diagram of mixed traffic flow can be obtained, as shown in Figure 7. Figure 7 represents that traffic capacity increases with the MPR of CAVs, and the growth is more visible when the MPR is within 0.6 to 1. is indicates that the traffic capacity of mixed traffic flow is positively correlated with the MPR of CAVs. Further analysis shows that when the MPR increases, the number of CACC and ACC increases gradually, and many vehicles can realize V2V real-time communication. In this scenario, when the leader vehicles suddenly change their driving behavior, the following vehicle can take reasonable countermeasures quickly. erefore, the traffic capacity of mixed traffic flow can be improved by maintaining a small headway.
Moreover, when the MPR is 1, the mixed traffic flow can be regarded as homogeneous CAVs traffic flow. In this case, CAVs realize cooperative driving through intelligent sensor devices. Currently, the cell transmission capacity and road traffic capacity reach the maximum.
is means that the large-scale use of CAVs is conducive to improving traffic efficiency, thus reducing delay and congestion, and improving traffic flow performance.

VCTM Based on Mixed Traffic Flow.
e actual traffic link length is not equal, and the classical CTM is used to divide the road link into isometric cells, which may not accurately describe the characteristics of mixed traffic flow. Furthermore, the division standard is the product of the simulation step length and the free-flow velocity. When the length of the road section cannot be divided equally by the cells, as shown in Figure 8, the method of reducing the simulation step length is usually adopted. However, this method will reduce the simulation efficiency and is not conducive to CTM's simulation of large-scale transportation networks. To solve this problem, the VCTM is proposed. e VCTM should meet the requirements, x i ≥ v f Δt, that is, the length of the cell is not less than the distance the vehicle passes through in time at free-flow speed. Based on the assumption of uniform distribution of vehicles in the same cell, that is, the traffic density in the cell is evenly distributed [8,9,36], the sending capacity S i (t) and receiving capacity R i (t) of cell i are defined by Equations (22)-(23).
Both sides of Equations (22) and (23) are divided by Δt: where q max i+1 (t) is the maximum flow of cell i + 1 at time t, and k j i (t) is the congestion density of cell i at time t, j is the jam state.   In CTM, the cell length is fixed and equal, and the number of cell vehicles is used to represent the cell state. However, the cell length changes in VCTM, and the traffic density is used to describe the cell state.
where k i (t + 1) is the traffic density of cell i at time t.
Based on the three-phase traffic flow theory [37], there is a substable region in the actual traffic flow. When the density is greater than the critical density, 'ghost jam' will occur. erefore, the traditional triangle fundamental diagram is improved to consider the influence of the hysteresis phenomenon, as shown in Figure 9. It can capture the actual traffic flow characteristics more accurately. Figure 9 shows that when the free-flow direction changes to the congested flow, the traffic flow density is k 1 , and when the congestion density changes to the free flow, the traffic flow density is k 2 , and k 1 < k 2 .
For better understanding, we introduce a parameter φ i (t), which is defined as the traffic state of cell i at time t. When φ i (t) � 0, it means that cell i is free-flow at time t, and when φ i (t) � 1, it means that it is congested flow at time t.
It is easy to know that the sending capacity of a cell is determined by the traffic density of the upstream cell and the outflow capacity of the downstream cell. e receiving capacity of the cell is mainly determined by the traffic state of the cell. erefore, the expressions of cell transmitting and receiving ability can be obtained by Equations (29)- (30).
en, equations (24), (28)- (30) can be used to describe the traffic flow characteristics of the VCTM in the basic link, and the merging and diverging link models remain unchanged.

Model Validation
e existing vehicle trajectory data are HDVs, such as NGSIM [38] and pNEUMA [39] databases. e trajectory data of CAVs are lacking in these databases. e mixed traffic flow studied in this work is composed of HDVs and CAVs, and these databases cannot be used to calibrate the proposed model. Instead, the driving behavior of CAVs is easy to implement in micro-simulation software [40]; thus, this section verifies the validity of the proposed model based on micro-simulation software.

VCTM Simulation.
Due to the frequent and serious consequences that arise from expressway traffic crashes, we focus on capturing the characteristics of expressway traffic flow [41]. e basic link of the expressway with a total length of 4500 m is selected as the simulation segment. e segment is divided into 30 cells with different lengths based on the method proposed by Hu et al. [42], as shown in Figure 10. e simulation time step is 3 s, the simulation time is 1200 s, the free-flow velocity is 33.3 m/s, and the input flow is 1200 veh/h. ere is a traffic accident in the 18th cell at 300 s, which leads the output capacity to be reduced to 0 veh/h. erefore, this accident caused a large area of traffic congestion. en, the accident treatment is completed at 600 s, and the congested vehicles begin to evacuate. In addition, the MPR of CAVs are set as 0, 0.2, 0.4, 0.6, 0.8, and 1, respectively. e macro numerical simulation experiment is carried out by MATLAB, and the change of traffic density of each cell with time and space is obtained, as shown in Figure 11. Figure 11 shows that after the traffic accident at the 18th cell at 300 s, the traffic density of the cell increases rapidly; the former cell starts to queue and waits until the traffic accident is cleared. From 600 s, the congestion begins to dissipate, the vehicle returns to a normal state, and the cell density decreases gradually. From Figure 11, the queue propagates to the 8th cell and the queue dissipation time needs 912 s when p is 0. When p is 1, the queue propagates only to the 10th cell and the queue dissipation time only takes 741 s. is indicates that the larger the MPR of CAVs is, Journal of Advanced Transportation the fewer congestion cells and the shorter the queue dissipation time are. Moreover, the 23.08% reduction of queuing dissipation time of homogeneous CAVs traffic flow than that of homogeneous HDVs. is is because when the MPR is low, there are many HDVs in mixed traffic flow, which cannot realize V2V communication. e headway in the VCTM is adopted with a more significant value to prevent collision between two adjacent vehicles. e more significant headway leads to the delay in transmitting more cells, so the queue dissipation time is longer. When the MPR is high, the number of vehicles accommodated by the cell increases with minor headway. en, the number of delay cells decreases, and the queuing dissipation time is shorter. is indicates that the capacity of road congestion dissipation increases with the MPR of CAVs.

Microsimulation.
To verify the effectiveness of the macro VCTM, the micro-simulation of mixed traffic flow in Section 3.1 is carried out. en, the results are compared with the numerical simulation results of the VCTM.
More recently, to accurately describe the car-following behavior of mixed traffic flow, many scholars have used different car-following models to study the characteristics of traffic flow from multiple perspectives at the microlevel, such as intelligent driver model (IDM) [43], adaptive cruise control (ACC) model [44], and cooperative adaptive cruise control (CACC) model [45]. On the basis, the stability and safety of mixed traffic flow are analyzed. In terms of stability, Wang et al. [46] use the CACC model and ACC model to describe the car-following behavior of CAVs with connected communication and connected degradation, respectively, while for HDV, the IDM is used to reflect the car-following behavior. en, it is found that the degradation of the CAVs leads to the decreasing of flow stability remarkably. In terms of safety, Yao et al. [28] also used the CACC model, ACC model, and IDM to describe the car-following behavior in mixed traffic flow. Furthermore, further analyzed the safety of mixed traffic flow by using indicators such as standard deviation of vehicle speed (SD), time exposed rear-end crash risk (TER), exposed time-to-collision (TET), and integrated time-to-collision (TIT).
To sum up, no matter to analysis the stability or safety of mixed traffic flow, most of the existing studies use different car-following models to describe the car-following behavior between CAVs and HDVs. Whereas, there are significant differences in the car-following mechanism and expression between different car-following models, in practical application. In addition, the differences will make human drivers feel strange and inadaptable. What is more, according to the research [47,48], the unified car-following model with deep human imitation is more conducive to simulating the real traffic flow. In other words, using the same car-following model to describe the car-following behavior of CAVs and HDVs will be more matched with the driving habits of human drivers. erefore, in this study, IDM with different MPR is used to describe the car-following behavior of CAVs and HDVs. Its equation will be presented as where a is the maximum acceleration, b is the comfortable deceleration, v is the speed of the rear vehicle, v f is the freeflow speed, Δv is the speed difference between the front and the rear vehicle, δ is the acceleration index, and its value is 4, S 0 is the minimum distance of the vehicle, and its value is 2 m, T m is the safe headway of vehicles in m-th car-following mode (the value as shown in Table 1), and h m is the headway of vehicles in m-th car-following mode, L is the length of the vehicle, and its value is 5 m. e basic link of the expressway with a total length of 4500 m is selected as the simulation segment, which is divided into 30 cells with different lengths, as shown in Figure 10. e simulation parameter setting is consistent with Section 3.1.
e microscopic simulation experiment is carried out by SUMO, and the change of traffic density of each cell with time and space is obtained, as shown in Figure 12. Figures 11 and 12 show that the results of macro VCTM simulation are consistent with those of micro-simulation. With the increase of MPR of CAVs, the queue dissipation time decreases gradually. e ability to alleviate congestion is also improved.

CTM Simulation.
In order to compare the applicability of VCTM and CTM models, the equal length simulation section is selected, which is divided into 45 links with equal lengths, as shown in Figure 13. e simulation parameter setting is consistent with Section 3.1. e numerical simulation experiment is carried out by MATLAB, and the change of traffic density of each cell with time and space is obtained, as shown in Figure 14.
Comparing Figure 14 with Figures 11 and 12, the overall trend is roughly the same; that is, with the increase of MPR of CAVs, the queue dissipation time gradually decreases.

Comparison of Results.
To quantitatively describe the difference between CTM, VCTM, and microscopic simulation, the queue dissipation time is selected as the analysis index to calculate mean absolute error (MAE) and mean absolute percentage error (MAPE). e results are shown in Table 2. Table 2 shows that with the increase of the MPR of CAVs, the congestion dissipation time decreases gradually. According to the error results in Table 2, compared with the microscopic simulation results, the errors of the two models are minimal. e maximum value of MAE of the CTM was 34 s, and the maximum value of MAPE was 4.18%. Moreover, the maximum value of MAE of VCTM was no more than 22 s, and the maximum value of MAPE of VCTM was 2.82%. erefore, compared with CTM, the error of VCTM is smaller, the average error is reduced by 0.85%, and the relative error is reduced about two times. is means that the VCTM can accurately describe the characteristics of mixed traffic flow in the actual situation.   In general, the simulation results of VCTM are closer to reality in the same simulation environment. erefore, the VCTM can reflect the physical phenomena of queuing propagation and dissipation well. is is because the CTM needs to divide the research links into equal length cells for simulation experiment and does not consider the hysteresis phenomenon. On the contrary, the VCTM is proposed based on the improved three-phase traffic flow theory, which considers the substable region of the actual traffic flow in congestion.
erefore, the VCTM can analyze the whole evolution process of normal, congestion, and dissipation of unequal length segments. To sum up, the VCTM not only increases the flexibility of road traffic simulation but also make the simulation results more accurate and fit with the actual situation.

VCTM Simulation for Mixed Traffic Flow.
e expressway with a basic link, merging link, and a diverging link is selected as the study road segment. e total length is 1300 m, and it is divided into nine cells with unequal lengths, as shown in Figure 15. Among them, the simulation step is 3 s, the free-flow speed is 33.3 m/s, the input flow is 1200 veh/ h, the main road merging coefficient is 0.6, the ramp entrance is 0.4, the main road diverging coefficient is 0.7, and the ramp exit is 0.3. In addition, the traffic accident is set at 700 m section, and the duration is 300 to 600 s. To study the influence of the MPR of CAVs on physical phenomena, such as mixed traffic flow congestion and queuing, the MPR is set as 0, 0.2, 0.4, 0.6, 0.8, and 1.0, respectively. e numerical simulation experiment is carried out by MATLAB, and the change of cell flow with time on the road segments is obtained, as shown in Figure 16.
As shown in Figure 16, the traffic accidents occurred at 300 s, the traffic flow of each cell decreased rapidly, and the traffic volume reached the highest when the congestion began to dissipate. When p � 1, the maximum volume is 2975 veh/h, and the road capacity recovery time is 438 s. When p � 0, the maximum volume is only 2105 veh/h, and the road capacity recovery time is 501 s. e homogeneous CAVs traffic flow capacity can reach 1.41 times of the homogeneous HDVs traffic flow, and the congestion dissipation time can be reduced by 25%. e results show that with the increase of the MPR of CAVs, the road capacity increases gradually; that is, CAVs are conducive to improving the efficiency of the traffic system.
In addition, compared with the basic link, the traffic volume of the merging and diverging link has obvious changes. is indicates that the location of the ramp entrance and exit has a specific impact on the capacity of the expressway. e changing length of cells can make the road traffic layout more reasonable. erefore, the VCTM can more accurately describe the mixed traffic flow characteristics of the expressway.
To understand the congestion dissipation mechanism of mixed traffic flow well, the VCTM simulation results are further analyzed.

Congestion Analysis.
To further study the influence of the MPR of CAVs on the congestion dissipation capacity of the expressway, the VCTM is used to simulate the mixed traffic flow of Section 4.1. en, the average travel speed, congestion delay, and congestion scale are combined to analyze and evaluate the road congestion dissipation capacity [49][50][51][52].

Average Travel Speed.
e average travel speed refers to the average travel speed of all vehicles in the traffic network. It is an important index to evaluate road traffic smoothness. e calculation formula is as follows based on the cells.
where v min is the minimum speed corresponding to congestion density (generally 13 km/h). Figure 17 shows that the evolution law of traffic congestion and queue dissipation after the traffic accident of 300 s to 600 s at 700 m of the simulation section. At 300 s, vehicles in the fifth cellular began to queue and wait. e congestion spread forward the sequential cell. e average travel speed of the vehicles in the whole simulation road suddenly dropped below 5 m/s. When the traffic accident was cleared up at 600 s, the congestion began to dissipate, and the vehicles in each cell gradually recovered to normal driving speed. In addition, when the road is the homogeneous HDVs traffic flow, the vehicle speed recovery time needs 804 s after congestion dissipation. While the speed recovery time of the homogeneous CAVs traffic flow only needs 738 s. erefore, the recovery time of speed is inversely related to the MPR of CAVs; that is, the greater MPR of CAVs, the shorter the recovery time of vehicle speed. is indicates that the CAVs traffic flow has a strong recovery ability after congestion evacuation, which is conducive to alleviating the traffic congestion and improving the efficiency of the traffic system.

Congestion Delay.
Individual congestion delay refers to the difference between the number of vehicles in cell i and the number of vehicles leaving at time T. e overall congestion delay reflects the operation efficiency of the whole traffic system and the state of the road link at the beginning and end congestion. We assume the cell number is I, and the overall congestion delay is presented as Figure 18 shows that the road congestion and delay in the whole simulation period, where the cell number is 9. At 300 s, the congestion begins to spread from the fifth cell, and the delay of the whole road increases rapidly. e congestion begins to dissipate at 600 s, and the delay gradually decreases. In addition, the delay of traffic congestion gradually decreases with the increase of the MPR of CAVs. erefore, the congestion delay is inversely proportional to the MPR of CAVs; that is, the higher the MPR is, the lower the congestion delay is, and the shorter the dissipation time is. is suggests that the large-scale application of CAVs is helpful to reduce traffic congestion and delay.

Congestion Scale.
Congestion scale is the number of congested road segments to the total number of road segments in the traffic network. In the VCTM, it is the ratio of the number of congested cells to the total number of road cells. erefore, the congestion scale is an important indicator reflecting the congestion situation and traffic capacity of the whole traffic network. Its calculation formula is as follows: where M C is the number of congestion cells, and M T is the total number of road cells. Figure 19 shows the degree of traffic congestion in the simulation period. At 300 s, the number of congestion cells began to increase. e congestion began to dissipate at 600 s, and the number of congestion cells gradually decreased. In addition, when the traffic system is homogeneous HDVs traffic flow, the maximum congestion degree reaches 70%, while that of homogeneous CAVs traffic flow is only 50%. erefore, the traffic congestion degree is inversely related to the MPR of CAVs. e greater the MPR, the lower the congestion degree. is indicates that the CAVs traffic flow has a strong ability to improve the traffic congestion after the traffic accident.

Results.
e results of congestion analysis show that the congestion delay of mixed traffic flow decreases with the increase of the MPR of CAVs. e reason may be that the vehicle can make safe and efficient decisions on accelerating and decelerating in the CAVs environments [53]. erefore, CAVs can reduce the congestion degree of the whole traffic system. In addition, the traffic capacity will not increase infinitely with the increase of the MPR of CAVs, which is related to the free-flow speed, the minimum safe distance, and the expected headway [29].

Conclusions and Future Work
is paper proposes a VCTM to analyze the characteristics of mixed traffic flow with CAVs and HDVs. Based on the simulation, the following conclusions can be drawn.
(1) e mixed traffic flow is proportional to the MPR of CAVs, that is, with the increase of the MPR, the capacity of the mixed traffic flow increases gradually. e closer the MPR is to 1, the more significant the improvement effect on the traffic capacity is.
(2) Compared with the classical CTM, the VCTM can more accurately reflect the changes of traffic volume and density with time in different types of road sections (merging, diverging, and basic road links). e VCTM can clearly describe the physical phenomena of congestion and dissipation in each road segment.
(3) e results of congestion analysis show that the congestion delay of mixed traffic flow decreases with the increase of the MPR of CAVs.
is paper studies the characteristics of mixed traffic flow based on VCTM and draws some exciting conclusions. However, there are two regrets in this paper. First of all, we did not consider the influence of different headways, traffic demand, and ramp entrance and exit location on the characteristics of mixed traffic flow. In addition, we have not been extended the VCTM to the research of multiple lanes. erefore, more deep research can be conducted from the following aspects in the future: (1) different factors will be introduced for sensitivity analysis, and the influence of different trajectory planning strategies of CAV on the characteristics of mixed traffic flow will also be studied. (2) e model is further extended to study multiple lanes of the expressway.

Data Availability
No data were used to support this study.

Conflicts of Interest
e authors declare that they have no conflicts of interest.