Abstract

The vehicle nonstrict priority give-way behavior (VNPGWB) is a common part of traffic interaction between motorized and nonmotorized vehicles in many countries. This study proposes a mixed-flow cellular automaton model to simulate the passing of vehicles in front of bicycles at crosswalks. The mixed-flow model combines a vehicle model with a bicycle model, using nonstrict priority give-way and strict give-way two driving behaviors defined as relating to the decision point rule and the launching rule, respectively. Simulation results showed that as the vehicle and bicycle inflow rates increased, a critical inflow rate divided vehicle and bicycle traffic flow into free flow and saturated flow conditions. The values of vehicle saturation flow decreased from 0.34 to 0.05, and the values of bicycle saturation flow decreased from 0.54 to 0.44, indicating that the mixed traffic flow has a negative effect on vehicle and bicycle saturated flow. Results also showed that VNPGWB effectively improves vehicle saturation flow over that of the strict give way. The advantage of VNPGWB is more significant when vehicles and bicycles are in saturation traffic flow.

1. Introduction

In many countries in Europe and Asia, the bicycle is still an important mode of transportation. For example, in China, more than 38% of commuters choose bicycles as their main travel mode [1, 2]. It is noteworthy that electric bicycles have grown rapidly in recent years, with the number of e-bikes exceeding 250 million in China [3, 4]. In addition, shared bicycles also have gained tremendous popularity as a result of their convenience [5, 6]. However, due to their vulnerability in a collision, cyclists are subject to higher safety risks than drivers of vehicles [7, 8]. Collisions between bicycles (including e-bikes) and vehicles tend to cause very severe injuries and fatalities. Some studies show that 90% of all cyclist fatalities are caused by collision with vehicles [6, 911]. In fact, vehicles and bicycles competing for priority are a common reason for conflicts and crashes at crosswalks [12]. Most traffic managers believe that drivers should strictly give way to bicycles passing through a crosswalk [3, 13, 14]. However, actual investigations have found that drivers in many countries and regions, such as Norway, Finland, Germany, and China, do not strictly adhere to this rule [8, 1519]. That is, during the real traffic conditions, vehicle drivers may not always assume they must comply with the rule of giving way, and instead, they compete for priority [17, 18]. For example, the decision a driver makes to pass through the crosswalk or not depends on the bicycle’s position and speed. Since it is based on the driver’s subject assessment, the behavior is defined as vehicle nonstrict priority give-way behavior (VNPGWB).

In the process of vehicles passing in front of bicycles at crosswalks, drivers have only two choices: passing without stopping or stopping to give way [17]. Drivers should generally adapt their speeds to avoid endangering cyclists at crosswalks, and if necessary, drivers should stop to give way to cyclists [12]. This rule ensures the orderly passing of vehicles and bicycles at the crosswalk. However, when the vehicle and bicycle are in a condition of saturated traffic flow, the vehicle flow is usually blocked by the bicycle flow, which has a significant impact on the traffic capacity of the vehicle flow, especially near unsignalized intersections [15, 17]. Specifically, there is a great impact on traffic flow when vehicles are waiting for bicycles at crosswalks. In this case, nonstrict priority allows vehicles to leave the conflict area earlier, mitigating their effects on subsequent vehicles and thereby improving the capacity of the road and intersection. Therefore, despite the negative impact on traffic safety, many countries acquiesce to nonstrict priority in the practical management of traffic, administering no punishment for violators as long as no crash occurs [17, 19].

This study uses a cellular automaton (CA) model to simulate a vehicle passing in front of bicycles at crosswalks engaged in VNPGWB. In order to understand VNPGWB, the key focus is to combine a mixed-flow model (i.e., integrating a vehicle model and a bicycle model) with driving behaviors. Driving behaviors in this study are primarily divided into the decision point rule and the launching rule. By adjusting the proportion of driving behaviors (strict give way and nonstrict priority give way) in these two rules, the real traffic conditions for VNPGWB can be simulated.

The rest of this paper is organized as follows: Section 2 reviews the literature; Section 3 describes the proposed mixed method in detail; Section 4 measures numerical simulation; and Section 5 concludes the paper.

2. Literature Review

Few previous studies explored VNPGWB, Räsänen [16] was the first to examine the behavior of road users changed based on their knowledge of priority regulations, such as whether a turning vehicle needs to give way to a cyclist from the same or opposite direction. Results showed that the effect of priority regulations on road user behavior depended on the characteristics of the bicycle crossings. Most subsequent studies found that drivers choose a give-way behavior depending on their understanding of priority order. For example, Lin et al. [17] developed a microdriving force model, which included safety driving force and efficiency driving force, for right-turning drivers, who constitute the dominant group confronted with the nonstrict priority passing situation. Silvano et al. [7] presented a modeling framework to describe driver-cyclist interactions when they were approaching a conflict zone. In the framework, driver yielding, or give-way, behavior is modeled as a function of several explanatory variables. Ma et al. [8] established a three-layered mathematical model, including a decision layer, operation layer, and constraint layer, to simulate the variation in trajectories of right-turn vehicles. Bai et al. [19] estimated the capacity of left-turn vehicles under nonstrict priority. Results showed that the model was valid for estimating the capacity of an exclusive left lane with a permitted phase under nonstrict priority. In summary, most studies have assumed that all drivers follow VNPGWB; however, drivers who choose to strictly give way coexist in mixed traffic flow with those who do not. Further, a given driver may make different decisions under free flow and saturated flow conditions.

When passing bicycles at crosswalks, most drivers constantly adjust their vehicles’ speed. In order to simulate this behavior, this study will establish a microscopic simulation model. The most popular microscopic simulation methods used in previous studies are car-following models and cellular automaton (CA) models [2022]. Because the car-following model primarily simulates the interaction between a front and rear vehicle, the model is only applicable to one-dimensional movement and cannot represent the full mixed traffic flow. The CA traffic flow model, in contrast, can make full use of computer operations to flexibly change its rules according to various traffic conditions. In fact, the movements of vehicles and bicycles are discrete, and the CA model can use discrete space-time and state variables to regulate evolution rules and thus to describe nonlinear behavior. Additionally, the CA model can simulate the gradual change in vehicle and bicycle behavior as conditions change from free flow to saturated flow [2325]. It can be found, through long-term simulation, whether there is a phase transformation from free flow to saturation flow. In recent years, a large number of studies have used CA models to simulate mixed traffic flow. For example, Meng et al. [26] proposed a single-lane CA model to simulate mixed traffic with motorcycles and investigated the relationship between motorcycle lane changing behavior and density of traffic flow. Zhao [24] established a mixed bicycle traffic model, comprised of two bicycle types of bicycles, that shows bicycle traffic characteristics on eight physically separated bicycle paths in China. Ren et al. [21] improved a cellular automaton model by incorporating social forces that can describe interactions between pedestrians, making it useful for modeling the bidirectional pedestrian flow at crosswalks. Lu et al. [25] proposed a simulation model to represent vehicles yielding to pedestrians at crosswalks and demonstrated the relationship between saturated flow and yield behavior. The CA model can demonstrably be applied to simulate the behavior of vehicle drivers passing bicycles at crosswalks.

3. Model

The crosswalk is one of the most serious conflict areas between vehicles and bicycles. Obviously, cyclists are vulnerable road users in this conflict area. In order to avoid a crash, a vehicle driver’s decision process begins from upstream of the conflict area, where the driver confirms whether there is a potential conflict or not. The decision process of giving way behavior can be divided into two situations, as shown in Figure 1. Situation 1 (Figures 1(a)–1(c)): the white car (red box) has observed the bicycle (red box) at the crosswalk and chooses to pass without giving way at all; Situation 2 (Figures 1(d)–1(f)): the blue car (red box) has observed the bicycle (red box) at the crosswalk and stops to give way. Based on field observation, a vehicle’s passing through bicycles at the crosswalk can be divided into the following steps: (1) the vehicle reaches the decision point; (2) the vehicle stops outside the conflict area; (3) the vehicle passes through the conflict area. It is worth noting that VNPGWB only occurs in Steps 1 and 2. The schematic diagram of these steps is depicted in (Figures 2(a)–2(c)). When a vehicle needs to pass through a crosswalk, the driver’s decision process begins at a certain distance upstream of the conflict area, a distance defined as the decision point (Figure 2(a)). When a vehicle stops to give way to a bicycle, the driver usually adjusts speed and stops at the boundary of the conflict area. Once an acceptable gap exists, the driver will finish the traverse. The stop position in this study is defined as the launching point (Figure 2(b)). However, if the driver chooses to pass without giving way, Step 2 will not occur.

3.1. Definition and Delimitation of Cellular Space

A vehicle lane and a crosswalk were chosen as the research area for this study. First, the cellular space is defined as a two-dimensional matrix, and the nth cellular unit in the cellular space position is . Then, because of the different sizes of vehicles and bicycles, a finer cellular unit size is defined in order to comply with the actual speed and minimum vehicle space for each. One bicycle occupies two cellular units and one vehicle occupies 6 × 4 cellular units. The vehicle lane widths are  m and the bicycle lane widths are  m; thus, the overlap area between the vehicle lane and the bicycle lane is  m2. Finally, the conflict area in the cellular space position is shown in Figure 3.

3.2. Vehicle Model

According to the vehicle’s position in its lane, the vehicle model can be divided into three rules: the driving rule, the decision point rule, and the launching rule. When following the driving rule, the driver selects an acceptable distance in understanding the interaction rules with consideration of safety and other factors. Then, the decision point rule is used: when drivers observe the bicycle lanes, different drivers (strict give way and nonstrict priority give way) make different decisions about whether to give way. If the vehicle slows down and stops at the boundary of the conflict area, it will be launched through the area when there is an acceptable gap. The launching rule is used to determine the acceptable gap for the different types of drivers making different decisions.

3.2.1. Vehicle Rule 1: Driving Rule

This study assumes that the vehicle is not affected by lane changes and reverse lanes. All vehicles move at an expected maximum velocity, and they adjust their velocity in order to avoid collisions with vehicles in front of them. For a process of , the driving rule evolution is as follows:Step 1: acceleration; Step 2: slowing down; Step 3: randomization with probability ; if , then Step 4: movement; ,where and represent the position and velocity, respectively, of vehicle n; represents vehicle acceleration; is the number of cellular spaces between vehicle n and vehicle n+ 1; is the length of the vehicle; and is a random number between 0 and 1, where is the randomization with probability.

This model adopts an open boundary: when vehicles’ road position is updated, the positions of the head vehicle and the tail vehicle are monitored as and , at the moment of . If , vehicles at the velocity of will enter the cellular with an inflow rate of . At the point of exiting the road, if , vehicles will leave the road.

3.2.2. Vehicle Rule 2: Decision Point Rule

The decision point is defined as the place where the vehicle driver needs to decide whether or not to give way when encountering a bicycle. In order to mirror the process of a vehicle driver’s decision, a deceleration restriction is introduced in the decision point rule and the decision time is a time step. is the acceleration and D is the deceleration. It is important to note that the position of the decision point is not a fixed point, but varies with different drivers and is impacted by actual conditions. At the instant of t, the decision point is given bywhere and represent the position and velocity of vehicle n, respectively; is the conflict area boundary; the summation represents the distance from deceleration to stopping; and means the time from deceleration to stop.

When a vehicle is at the decision point, the driver adjusts to a suitable velocity to make sure the cyclist safely traverses the conflict area. Suitable velocity is defined as a velocity that ensures the vehicle not only can pass through the conflict area safely but also can decelerate to a timely stop. In order to capture the stochasticity of different drivers at the decision point, a binary variable () is presented aswhere represent the driver choosing to give way to the cyclist; represents the driver choosing not to give way, that is, selecting nonstrict priority give way; and represents the proportion of nonstrict priority give-way drivers. Thus, the suitable velocity is presented at the instant of t:where “Case I” is the vehicle’s maximum velocity that will avoid a crash with the bicycle; “Case II” represents the choices of different drivers to give way or not; and D is the deceleration of vehicles.

3.2.3. Vehicle Rule 3: Launching Rule

The vehicle stops at the launching point and waits for an acceptable gap. Once an acceptable gap exists, the driver will take the opportunity to finish the traverse. The acceptable gap is defined as the time interval in which a vehicle successfully traverses the conflict area. In order to capture the stochasticity of different drivers at the launching point, another binary variable () is presented aswhere represents the driver who needs to give way and waits for bicycles to traverse the conflict area; represents the driver who chooses nonstrict priority to give way, for whom the acceptable gap is equal to the launching time of the vehicles; and represents the proportion of VNPGWB drivers. Thus, the acceptable gap is given by

In order to prevent vehicles from remaining static for a long period of time, a threshold vehicle waiting time T is defined. If the waiting time of a vehicle is longer than T, the vehicle will launch into the conflict area, and bicycles will be forced to stop.

3.3. Bicycle Model

There is great flexibility in the movement of bicycles, including both lateral and vertical movement. This study chose a new cellular type of bicycle model to simulate one-way bicycle travel in order to explore the interference between bicycles and vehicles in conflict areas. This model consists of two steps, lateral movement and vertical movement, as shown in Figure 3. Both the two steps adopt parallel rules. A bicycle can move to the left, forward, and right: , , and describe the number of empty cellular units on the left front, front, and right front, respectively; and , describe the number of empty cell on the perpendicular left and right, respectively. The bicycle cellular coordinate is represented by , where means that there is no bicycle occupying this cell; and means that there is a bicycle occupying this cell; represents the lateral velocity of the bicycle; represents the vertical velocity of the bicycle, and means the bicycle moves to the left; means the bicycle moves to the right; is the lateral acceleration of the bicycle; and is the vertical acceleration of the bicycle.

3.3.1. Bicycle Rule 1: Lateral Movement

For the process of , the lateral movement rules are as follows:Step 1: acceleration; Step 2: slow down; , Step 3: randomization with probability ; if rand()>, then Step 4: movement;

This model adopts an open boundary: when the vehicles’ road position is updated, the positions of the head vehicle and the tail vehicle are monitored as and , at the moment of . If , bicycles at the speed of will enter the cellular with an inflow rate of . At the exit point of the road, if , the bicycles will exit the road.

3.3.2. Bicycle Rule 2: Vertical Movement

In this study, there are two main situations in which bicycles can move perpendicularly. First, if there is no space in the front, a cyclist can choose vertical movement. Second, if a bicycle moves laterally, vertical movement can be chosen if the number of empty cells on either side is more than or equal to the number in front. The two situations are illustrated in detail as follows:

Situation 1. If  =  =  = 0, the cyclist chooses vertical movement in order to pass the crosswalk as soon as possible, and . If , then , and .

Situation 2. If , , and are not all zero, then the cyclist chooses the largest space in front as far as possible, and . If there are two or more choices with the same maximum number of cells, then , and  = .

3.4. Algorithm for the Mixed-Flow CA Model

In order to simulate the mixed-flow CA model, the abovementioned methods are integrated into one executable algorithm, which is processed in a parallel computing setting using the Python programming environment. The algorithm (Algorithm 1) calculates the process of . In each step, the following three rules are conducted in order from first to last, and all are applied to every vehicle. When the algorithm ends, the updates are applied to all vehicles and bicycles in parallel. The complete framework of this algorithm’s four steps is shown in Figure 4.

Step 1: driving rule.
 Step 1.1: input the current vehicle position ; current vehicle velocity ; a randomization with probability .
 Step 1.2: if the vehicle position is at decision point (equation (1)), then go to Step 2.
 Step 1.3: if the vehicle position +> boundary, then vehicles will leave the road.
 Step 1.4: the vehicle adjusts velocity according to the distance to the front vehicle; go to Step 4.
Step 2: decision point rule.
 Step 2.1: input the current bicycle positions and current bicycle velocities and .
 Step 2.2: assume that there exists a time , . If and , then the vehicle will crash with bicycles. This vehicle needs to stop at the conflict area boundary and go to Step 3.
 Step 2.3: if the driver chooses to give way, then go to Step 2.2; if the driver chooses nonstrict priority give way, then go to Step 4. The vehicle velocity is given by equation (4).
Step 3: launching rule.
 Step 3.1: input the current bicycle positions, current bicycle velocities and , and waiting time .
 Step 3.2: assume that there exists another time and the acceptable gap is given by equation (6). If the driver chooses to give way and , then the vehicle continues to wait; for waiting time, add 1; if the driver chooses nonstrict priority give way, then go to Step 4. The vehicle velocity .
 Step 3.3: if the waiting time is longer than T, the vehicle will launch into the conflict area, and bicycles will be forced to stop.
Step 4: movement.
 Update the position of vehicles and bicycles based on their velocities and ,. For the vehicle, , . For the bicycle, ,.

4. Numerical Simulation

4.1. Simulation Parameter Setting

The data inputs in the following simulations are set as follows. The cellular size is 1 m 1 m; the vehicle lane length () is 100 cell units, which therefore corresponds to 100 m; vehicle lane width () is 4 cells, which corresponds to 4 m; the crosswalk length () is 100 cell units, and it corresponds to 100 m; crosswalk width () is 6 cell units, and it corresponds to 6 m and the size of the conflict area is 6  4 cell, and it corresponds to 6  4 m2. The simulation parameters are shown in Table 1. In order to obtain vehicle and bicycle traffic flow data, the virtual detectors are set at the boundaries of the vehicle and crosswalk. When a vehicle or bicycle leaves the lane, the detector’s counter adds one. Moreover, and qb are the vehicle and bicycle flow, respectively, and the unit is veh/(time  lane). The simulation steps are 100,000 s, and the preceding 20,000 s are abandoned.

4.2. Model Validation

To verify the reliability of the mixed-flow CA model, the traffic flux-density relationship is the fundamental diagram that describes driver-cyclist interactions. With settings of  = 0.1,  = 0.1, and T = 30 s, Figures 5(a) and 6(a) show the microscopic fundamental relationships between flow and inflow rate for the vehicle and bicycle. When the bicycle’s inflow rate  = 0, there is no bicycle passing the crosswalk, and the vehicle flow shows an upward trend and remains stable at 0.34 (Figure 5(b)). When bicycle’s inflow rate  = 1, the maximal point of vehicle flow is 0.05 (i.e., saturation flow), which is obviously lower than  = 0. With more and more bicycles entering the mixed-flow cellular space, drivers must slow down and stop outside the conflict area to avoid crashes. Thus, only a few vehicles cross the conflict area. However, this situation does not occur for bicycle flow , because bicycles have priority in passing over the crosswalk. The saturation flow rate of bicycles drops only slightly when bicycles are waiting for vehicles to cross the conflict area. Thus, in the proposed mixed-flow CA model, the transition phase from the free flow to saturation is observed (Figure 7(a)). As shown in Figure 7(a), conflicts between vehicles and bicycles in the traffic system result in a drop in saturated flows, demonstrating that the model can reveal the interactions between vehicles and bicycles in the mixed traffic system.

4.3. Vehicle and Bicycle Flow Transition Phase

As can be seen in Figures 5(a) and 6(a), each curve has an obvious turning point (i.e., critical inflow rate) (), which divides the flow into two regions: the free flow and saturated flow. This section mainly analyzes the transition phase between the vehicle flow and bicycle flow .

For <, the vehicle flow is free flow, and the vehicle flow only depends on its own inflow rate . Conversely, for >, the vehicle flow is saturation flow, the vehicle flow is independent of its own rate, and it reaches its saturation flow , as shown in Figure 5(a). However, with an increase of bicycle inflow rate , both the critical value of the inflow rate and the saturation flow value diminish until a minimum value is reached. To illustrate the relationship between the inflow rate and the vehicle flow , Figure 5(b) shows four curves of different bicycle inflow rates . From this diagram, it is evident that the values of vehicle saturation flow decrease from 0.34 to 0.05, and the critical values of inflow rate decrease from 0.38 to 0.08. The drop ratio of is about 85%. Similarly, for <, the bicycle flow is free flow, and for >, the bicycle flow is saturation flow. The bicycle flow is independent of its own rate and reaches its saturation flow value , as shown in Figure 6(b). The values of bicycle saturation flow decrease from 0.54 to 0.44, and the critical values of inflow rate decrease from 0.64 to 0.48. The drop ratio of is about 19%.

It is noteworthy that the vehicle (bicycle) critical inflow rate () gradually decreases with the bicycle (vehicle) inflow rate (). Moreover, the collective effect of vehicles and bicycles only appears when and surpass their critical value. The critical inflow rate () is calculated and presented in the phase diagram in Figure 7. This transition phase can be classified into four zones; for example, line1 and line3 are the boundaries of Zone I (III) and Zone II (IV), which correspond to the critical value of the vehicle inflow rate . In Zone I and Zone III, the vehicle is in free flow, whereas in Zone II and Zone IV, it is in saturation flow. As the bicycle inflow rate () increases, the vehicle critical inflow rate () first decreases and then remains stable. Similarly, the bicycle critical inflow rate shows the same trend (line2 and line4). With the gradual interaction between the two traffic flows, the vehicle’s critical inflow rate and bicycle’s critical inflow rate reach equilibrium at cross point o( ).

In a word, the mixed-flow CA model effectively illustrates the transition phase from free flow to saturation for both vehicles and bicycles. It is interesting that the collective effect of the vehicle flow and bicycle saturation flow only appears when and surpass their critical inflow rate.

4.4. Comparison between Give-Way and Nonstrict Priority Give-Way Behavior

Demonstrating the effects of VNPGWB in mixed traffic flow, the proportion of VNPGWB drivers and has increased from 0.1 to 0.9, respectively. Figure 8(a) shows the relationships between flow and inflow rate for the vehicle and bicycle in the case of  = 0.9,  = 0.9, and T = 30 s. The vehicle saturation flow increases from 0.34 to 0.36 when  = 0, while the vehicle saturation flow increases from 0.08 to 0.12 when  = 1. When the bicycle inflow rate is low, the VNPGWB driver will quickly pass through the crosswalk, whereas the strict give-way driver will slow down to pass through the crosswalk due to caution. Thus, the vehicle saturation flow resulting from nonstrict priority behavior is slightly higher than from strict give-way behavior. With the further increase of the bicycle traffic flow rate, the VNPGWB driver can pass through the crosswalk more easily than can the strict give-way driver. In this way, VNPGWB can indeed improve the vehicle saturation flow, as the vehicle saturation flow of nonstrict priority is 1.5 times larger than that of strict give way. Figure 8(b) shows the effect of VNPGWB on vehicle saturation flow in the decision point rule and launching rule. The vehicle saturation flow decreases nonlinearly with an increasing proportion of VNPGWB drivers, with the advantage of VNPGWB being more significant in the decision point rule than in launching rule. As shown in Figure 7(b), VNPGWB also affects the vehicle and bicycle critical inflow rates. For  = 0.9 and  = 0.9, the boundary of line1 and line3 moves up and boundary of line2 and line4 moves down. In other words, the proportion of free flow vehicle traffic increases with the increase of VNPGWB.

To further verify that the proposed model can simulate VNPGWB in the real world, the space-time trajectory diagrams for the vehicle lane are shown in Figures 9 and 10. Obviously, the microscopic model proposed in this study is primarily used to describe different microlevel driving behaviors. When the vehicle and bicycle are in free flow (Figure 9), most vehicles pass through the crosswalk normally, while some drivers choose to slow down and temporarily stop to give way to bicycles at cellular spaces between 40 and 60. As mentioned in Section 3.2.2., the decision point is not a fixed point, which is also depicted in the space-time trajectory diagrams. When the vehicle and bicycle are in free flow, only  = 0.1 produces congestion at cellular space 60, and it dissipates very quickly (Figures 9(a) and 9(b)). On the other hand, when vehicles and bicycles are in saturation flow, most vehicles need to slow down and stop outside the conflict area, as shown in Figure 10. In this case, the congestion cannot dissipate immediately. The in the launching rule has a more significant effect on vehicle saturation flow than does in the decision point rule. Due to bicycle saturation flow, the vehicles form a long queue and waiting for an acceptable gap. For  = 0.1, the stopped vehicle can only be launched when its waiting time exceeds the waiting threshold (Figures 10(a) and 10(c)). For  = 0.9, some of the space trajectories continue to pass into the conflict area, primarily due to nonstrict priority drivers sometimes choosing to follow the front vehicle through the crosswalk when the bicycle flow is high.

Overall, the space-time trajectory diagrams show that the proposed model in this study effectively simulates the interference between vehicles and bicycles. It is further demonstrated that VNPGWB can improve vehicle flow over strict give way. Especially, when the vehicle and bicycle are in saturation flow, this advantage of VNPGWB is significant. These conclusions can provide support for acquiescence to nonstrict priority behavior for practical management in those countries (e.g., China).

5. Conclusions

The crosswalk is one of the most serious conflict areas between vehicles and bicycles. With the aim of improving traffic flow without compromising safety, this study proposed a new mixed-flow cellular automaton model (CA) to simulate the vehicle nonstrict priority give-way behavior (VNPGWB) in the crosswalk conflict area. To consider the driving behaviors in the appropriate proportions, the proposed model simulates VNPGWB in real traffic conditions. This model was divided into three rules: the driving rule, the decision point rule, and the launching rule; VNPGWB occurs in the decision point rule and the launching point rule. The main results of the simulation model are as follows: (1)The mixed-flow CA model effectively illustrated the transition phase from free flow to saturation flow for both vehicles and bicycles. As the vehicle and bicycle inflow rate increased, there is a critical inflow rate in each curve which divided the traffic flow between free flow and saturated flow. Moreover, the collective effect of vehicles and bicycles appeared only when the inflow rate surpassed its critical value.(2)Phase diagrams showed that the mixed traffic flow had a negative effect on the saturated flow of both vehicles and bicycles. When the bicycle inflow rate increased from 0 to 1, vehicle saturation flow values decreased from 0.34 to 0.05, and the critical values of inflow rate decreased from 0.38 to 0.08. When vehicle inflow rate increased from 0 to 1, bicycle saturation flow values decreased from 0.54 to 0.44, and the critical values of inflow rate decreased from 0.64 to 0.48.(3)The simulation results showed that the vehicle saturation flow of nonstrict priority behavior is 1.5 times larger than that of the strict give way, indicating that VNPGWB can improve the vehicle saturation flow. The space-time trajectory diagrams confirmed these results and thereby demonstrate the reliability of the simulation model. The advantage of VNPGWB over strict give way is most significant when the vehicle and bicycle are in saturation flow.

The simulation results, however, still require calibration of observed data to confirm the model’s accuracy. Using field data to verify the model’s association of parameters is recommended in the future study. In addition, the lack of variation in the size and speed of the simulated vehicles and bicycles in this study was a limitation that will be improved upon in future models. Moreover, the proposed model can also consider machine learning methods such as logistic regression and random forest to predict driving behavior. Finally, similar studies can be conducted for the purpose of implementing traffic control at crosswalks with devices such as signal lights or stop signs. Nonetheless, this model provides insights that can be beneficial to traffic engineers and managers.

Data Availability

The simulation result data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by Jiangsu Provincial Key Research and Development Program (Social Development) Project (BE2019713).