
The remainder of the paper is organized as follows:
Section II formulates the problem that is solved in this paper.
In Section III the method for solving the stated problem is
presented, followed by Section IV where simulation results
illustrate the coordination algorithm. Section V concludes the
work and provides some possible extensions.
II. PROBLEM FORMULATION
We consider Nafully automated vehicles on a road
network with cross-intersection, path merges and path splits.
The road network is assumed to be fully in a confined area,
such that non-controlled traffic participants (e.g. manually
operated vehicles, pedestrians, bicyclists etc.) are absent. We
further assume that the paths of all vehicles, i.e., their routes
through the road network are known, that no vehicle reverses,
and that overtakes are prohibited.
A. Vehicle modelling
The motion of the vehicles along their path is described
by
˙pi(t) = vi(t)(1)
˙xi(t) = fi(pi(t), xi(t), ui(t)) (2)
0≤hi(pi(t), xi(t), ui(t)).(3)
where pi(t)∈Ris the position, xi(t)∈Rnthe vehicle
state, ui(t)∈Rmthe control input, with i∈ {1, . . . , Na}.
The state is subdivided as xi(t)=(vi(t), zi(t)), with the
speed along the path vi(t)∈Rand zi(t)∈Rn−1collecting
possible other states. The functions fiand hi, both assumed
smooth, describes the dynamics and constraints that capture,
e.g., actuator and speed limits respectively.
B. Vehicle modelling in the spatial domain
For confined site optimization, it is beneficial to optimize
the trajectories of the vehicles over their full paths. However,
the time it takes a vehicle to traverse a path is dependent
on the solution, and not known a-priori. Consequently, it
is inappropriate to plan the vehicle’s motion with time
as the independent variable. Due to this, the problem is
reformulated in the spatial domain, using that dpi
dt=vi(t)
and dt= dpi/vi(t). The formulation of the vehicle dynamics
(1) in the spatial domain is
dti
dpi
=1
vi(pi)(4)
dxi
dpi
=1
vi(pi)fi(pi, xi(pi), ui(pi)) (5)
0≤h(pi, xi, ui).(6)
where the position piis the independent variable.
C. Conflict zone modelling
Aconflict zone (CZ) is described by the entry and exit
position [pin
i, pout
i]on the path of each vehicle. From the
known positions, the time of entry and exit of vehicle iis
tin
i=ti(pin
i),tout
i=ti(pout
i), respectively. In this paper,
two types of conflict zones are considered, as depicted in
Figure 1: the “intersection-like” and the “merge-split”. We let
I={I1, I2, ..., Ir0}denote the set of all intersections in the
confined site, with r0being the total number of intersection
CZs, and let Qr={qr,1, qr,2, ..., qr,l}denote the set of
vehicles that cross an intersection Ir. In the intersection-
like CZ, it is desired to only have one vehicle inside the
CZ, i.e., not allowing the vehicle jto enter the CZ before
vehicle i6=jexits the CZ, or vice-versa. The order in
which the vehicles cross the intersection Iris denoted OI
r=
sr,1, sr,2, ..., sr,|Qr|, where sr,1, sr,2, ... are vehicle indices
and we let OI=OI
1,...,OI
r. A sufficient condition
for collision avoidance for the r-th intersection CZ can be
formulated as
tsr,i (pout
sr,i )≤tsr,i+1 (pin
sr,i+1 ), i ∈I[1,|Qr|−1],(7)
where tis determined from (4). In the merge-split CZ case,
let M={M1, M2, ..., Mw0}denote a set of all merge-
split zones, with w0being the total number of merge-split
CZs in the site and Zw={zw,1, zw,2, ..., zw,h}denote
the set of vehicles that cross the merge-split CZ Mw.
For efficiency, it is desirable to have several vehicles in
the zone at the same time, instead of blocking the whole
zone. This requires having rear-end collision constraints
once the vehicles have entered the CZ safely. In this case,
the order in which the vehicles enter the zone is denoted
as OM
w=sw,1, sw,2, ..., sw,|Zw|, and we let OM=
OM
1,...,OM
w. The collision avoidance requirement for
this w-th CZ is described with the following constraints:
tsw,i (pin
sw,i )+∆t≤tsw,i+1 (pin
sw,i+1 +c)(8a)
tsw,i ,ki+ ∆t≤tsw,i+1 (psw,i,ki−pin
sw,i +pin
sw,i+1 +c),
kin
sw,i ≤ki≤kout
sw,i (8b)
tsw,i (pout
sw,i )+∆t≤tsw,i+1 (pout
sw,i+1 +c),(8c)
i∈I[1,|Zw|−1].
That is, while in the CZ, the vehicles must be separated
by at least a time-period ∆tiand a distance cior by ∆tj
and cj, depending on if vehicle jis in front of vehicle i
or vice versa. This is equivalent to the standard offset and
time-headway formulation often used in automotive adaptive
cruise controllers.
D. Discretization
The independent variable is discretized as pi=
(pi,1, . . . , pi,Ni), where pi,Niindicates the end position for
vehicle i, and the input is approximated using zero order
hold such that u(p) = ui,k, p ∈[pi,k, pi,k+1[. The equations
(4), (5) are (numerically) integrated on this grid, giving the
“discretized” state transition relation
ti,k+1
xi,k+1=F(xi,k , ui,k, pi,k , pi,k+1)(9)
where Fdenotes the integration of (4), (5) from pi,k to
pi,k+1.