2.1. Roundabout Flow Control
Consider a typical roundabout shown in
Figure 2, which has four approaches indexed by
. The incoming flows on the four approaches are denoted by
, assumed constant during each time period
, and known at the beginning. In real world,
can be predicted using real-time loop detector data, while the length of
can be adjusted according to the changing rate of
. The queue length on approach
,
, is a state variable of
. In each time period,
starts from 0 and ends at
(the time period length). Merge-in flows,
, are control inputs, which represent the flows merging into the circular flow within the roundabout from approach
. The queue length dynamics can be described as
Apart from the non-negative constraint, it is also reasonable to have an upper bound
for queue length
In applications,
can be decided by the critical queue length that may cause spillback on approach
. Therefore, the state constraints are
To model the constraints of control inputs
, the roundabout segments are introduced and described. As seen in
Figure 2, the circular roundabout is separated by four approaches into four arc-shaped segments,
, and
. Since these segments are likely to have similar road design as they belong to the same roundabout, it is assumed that all segments have the same service capacity
(i.e., the upper bound of segment flow
). Note that
may consist of vehicles from every approach. For example, vehicles entering from approach 1 heading to approach 3 will drive through
and
, and thereby contribute to both
and
. Similarly, vehicles entering from approach 2 heading to approach 1 will also drive through
, thereby contributing to
. Denote flows from approach
heading to approach
as
, which can be represented as
, where
is the proportion of the flow heading to approach
in the merge-in flow from approach
. To describe the relationships between
and
, define the service tables of four segments as follows (note that this study assumes that vehicles will drive through each segment no more than once when passing through the roundabout). The element in row
and column
of the service table for
,
, is an indicator variable of whether
will pass through
. For instance,
, indicating the U-turn flow from approach 2 will pass through
in the roundabout; while
means vehicles coming from approach 3 heading to approach 2 will not drive through
. Therefore, the flow on
is
Then, the control input constraints can be represented as
The control objective is
where
is the estimated waiting time for vehicles at the end of the queue on approach
, and
is the number of incoming vehicles or vehicles joining the queue at the end on approach
. Therefore, the integration can be interpreted as the total estimated waiting time of vehicles arriving from
to
. Compared to the vehicle-level delay representation, Equation (6) provides an informative and computationally tractable estimation for the total system delay. The complexity of the problem does not increase with the number of vehicles, enabling its use in real-time applications. The optimal control problem defined by (1)~(6) falls in the domain of nonlinear optimal control, which generally does not have analytical solutions. Therefore, a direct collection method [
16] is applied to solve it, which converts the optimal control problem into a nonlinear optimization problem by discretizing state variables
and control inputs
into vectors
.
2.2. Vehicle Merge-in Decision Control
Given the optimal merge-in flows from the last subsection, the intersection controller needs to coordinate vehicles coming from four approaches on when they should merge to achieve the desired merge-in flows. Therefore, a probabilistic merge-in decision control model is developed to determine when the vehicles at the head of queues can start to merge in. In this way, the flow-level control inputs are converted into an optimal vehicle passing sequence, which can be connected to vehicle-level operations.
Let the spare capacity
of each segment
denote the difference between the capacity
and the flow on
, that is
. Define
and
as the upstream segment of approaches
and 1, respectively. Apart from the spare capacity from approach
’s upstream segment, there is additional
capacity available for vehicles on approach
, which can be described as follows:
Note that
can also be interpreted as the flow that diverges out at approach
(diverge-out flow). Therefore, the total available capacity for the merge-in flow of approach
(suppose its upstream segment is
is
. Let
denote the probability of the leading vehicle on approach
choosing to merge into the circular flow within the roundabout whenever there is available capacity. Then, the probability
should satisfy the following condition to achieve the desired merge-in flows:
where
is the upstream segment of approach
. With this condition met, the actual merge-in flow from approach
would be
. Following (3) and (7), it is proved that
, which guarantees that our decision control is applicable at all times.
Since the merge-in decision control is probabilistic, the actual merge-in flows given by this model are different from the desired merge-in flows. Therefore, the optimal system performance given by the flow control model is an upper bound to the actual system performance, which is illustrated in our numerical studies as well.
2.3. Vehicle Platoon Control
Once vehicles choose to merge into the roundabout based on the probabilistic decision control model, they will start the merge-in procedure and eventually cruise with the circular platoon until they diverge to the approaches they are heading to. Note that this subsection focuses on longitudinal control therefore does not consider non-holonomic vehicle dynamics or low-level vehicle control such as throttle/steering control. The merge-in procedure consists of a self-adjusting stage and a virtual platoon control stage.
Figure 3 illustrates an example of the two stages using the merge-in procedure of vehicle
. First,
waits at position
(marked by a red line) when it is at the head of the queue on approach 2. When the decision control model allows it to merge into a spare slot (labeled the matching slot of
) in the platoon,
enters the self-adjusting stage. The objective of this stage is to have vehicle
drive through a critical position
(marked by a blue line on approach 2) to enter the virtual platoon control stage at the desired speed
after some desired time
. The desired speed is also the desired cruising speed of the platoon. The desired time
is estimated as the time it takes for the matching slot to travel to a virtual critical position
(marked by the yellow line on the right part of the roundabout in
Figure 3) from the beginning of the self-adjusting stage. The virtual critical position
is a projection of
on the roundabout circle, that is
, where
is a merging point of the roundabout and the on-ramp from approach 2. When vehicle
travels from
to
, a virtual vehicle
, which is the projection of
, also travels from
to
. At the same time,
communicates with vehicles
and
, making
behave like a vehicle communicating with its surrounding vehicles in the circular platoon. This stage is labeled as the virtual platoon control stage. Note that when determining the length of
,
should be ensured, where
is the diverging point of the off-ramp to approach 2. It is assumed that once a diverging vehicle passes the diverging point and drives to the off-ramp, it is no longer a part of the circular platoon (i.e., its movement will not affect surrounding vehicles in the platoon anymore). This guarantees that during the virtual platoon control stage, there will be no real diverging vehicle at the slot where the virtual vehicle projected from a real merge-in vehicle is located. Otherwise, they may have a conflicting influence on the surrounding vehicles since they are at the same slot in the platoon.
The control of a vehicle from approach
in the self-adjusting stage can be formulated as
where
are the maximum speed, maximum deceleration, and maximum acceleration, respectively. As shown in
Figure 4, in a speed-time graph, the formulation seeks for trajectories starting from
ending at
, whose time integration equals to
. The formulation can have multiple feasible solutions if
, where
are the blue (right)-shaded and green (left)-shaded areas, respectively, in
Figure 4. The blue-shaded area stands for the maximum distance a vehicle can travel within
(the vehicle accelerates at
, all the way to
, then maintains the speed until an instant after which it can keep decelerating at
to arrive at
with desired final speed
). Similarly, the green-shaded area is the lower bound of the distance with the given constraints (9b-d). One feasible strategy is to follow the speed trajectory indicated by the yellow (dash-dot)/red (solid) lines when
is greater than or equal to/less than
. Specifically, the yellow dash-dot speed profile indicates that the vehicle accelerates at
, all the way to a specific speed
, then maintains the speed until an instant after which it can keep decelerating at
to arrive at
with the desired final speed
. The red solid line shows that the vehicle stops at
until moment
, then accelerates at
until the speed reaches
, and maintains this speed up to the end
. Both
and
can be derived from the constraint (9a).
The merge-in procedure will be seamless if the transition of two stages is as precise as described by (9a–d). However, it is possible that when the matching slot arrives at the virtual critical position, the vehicle is not exactly at the critical position at the desired speed because a critical parameter
is estimated and possibly inaccurate when solving (9a–d). For example, in
Figure 3, the matching slot of vehicle
arrives at the diverging point with the diverging vehicle
earlier than
. Additionally, it is possible that, at
,
still needs a little more time to get to the off-ramp. Then,
will have passed the critical point when
leaves the platoon. Note that the start time of the virtual platoon control stage is when the available slot arrives at the diverging point, not when the vehicle arrives at the critical point or the blue lines in
Figure 3. Therefore, the virtual platoon control should ensure stability to mitigate the initial speed and position errors introduced by virtual vehicles when they enter the virtual platoon control stage.
To ensure stability, the circular platoon in the roundabout is investigated, which consists of real vehicles, virtual vehicles projected from the vehicles on on-ramps at the second merge-in stage, and some spare vehicle slots. Note that here we assume a circular roundabout for simplicity; however, the model is not restricted only to ideal circular roundabouts. The platoon control model can be extended to non-circular roundabouts by mapping the roadways to circular roundabouts with the same road length. Suppose that when the roundabout serves at its maximum capacity, there are
homogeneous vehicles cruising in the roundabout with a desired angular speed
(where
is the radius of the roundabout). It implies that the roundabout has
vehicle slots labeled from 1 to
. When there are spare vehicle slots, the circular platoon can be converted into one (if there is only one spare vehicle slot) or more linear vehicle platoons using the spare slots to separate them. Then, existing linear CACC platoon control strategies [
17] can be applied.
For a circular platoon with no spare vehicle slots, a symmetric bidirectional control architecture is developed (where predecessor and follower position errors influence the vehicle controller symmetrically). As shown in
Figure 3, an angular coordinate
is used to describe the position and angular speed of vehicle
(i.e., vehicle matched with vehicle slot
). Each vehicle is modeled as a double integrator. The control inputs for vehicle
in the platoon are assumed to depend only on its speed error
and the relative headway errors between itself and its immediate neighbors (i.e., its predecessor and follower). Denote the positions of the preceding and following vehicles of vehicle
as
and
, respectively. The desired angular headway is
. The vehicle dynamics can be described as
where
are positive constants. To facilitate analysis, a coordinate change is considered using the initial position of slot 1,
, as a reference point:
To describe the dynamics of the circular platoon, define
. Substituting (11a–b) into (10a–b), the following equation is obtained:
Equation (12) indicates that the circular platoon to be controlled is a linear system. The following lemma [
18] illustrates the non-negativity of
’s eigenvalues.
Lemma 1 [
18]
. Given an circulant matrix with the formthe normalized eigenvectors are the Fourier modes, namely,where and
is the imaginary unit. The corresponding eigenvalues are then given by Note that
is a
circulant matrix with
. According to Lemma 1, the
eigenvalue of
is
where
. Note that
. Therefore,
Theorem 1 [
19]
. The equilibrium states for a general linear time-invariant system are stable if and only if each eigenvalue of satisfies and if is defective, where is the real part of .
From Theorem 1 (Theorem 5.4 of [
19]), the eigenvalues of the system can be described by (12) as follows:
With
’s eigenvalues
, we have
Recall that (18) has shown the non-negativity of . Therefore, for , and when , . According to Theorem 1, the equilibrium states of (12) are stable.