1. Introduction
Pavement condition significantly impacts traffic behavior. Pavement deterioration causes traffic accidents, congestion, pollution, and time delays [
1]. Moreover, poor roads impact the smooth flow of traffic, resulting in rider discomfort and increased vehicle operating costs [
2]. Congestion lowers vehicle speeds so emissions are increased [
1]. Further, vehicle speed is reduced by an average of
% when the road condition is poor compared to when it is excellent, and average emissions increase by
%. Road safety is a primary concern worldwide as road accidents cause
million fatalities each year [
3]. It is also dangerous as uneven pavement and potholes, damaged concrete, cracks, and exposed rebar can cause drivers to lose control, resulting in severe accidents [
4]. Efficient traffic forecasting and control are essential to alleviate traffic problems such as congestion and improve road infrastructure [
5]. This requires a practical model for traffic prediction.
Traffic models are typically microscopic or macroscopic, and mesoscopic. Macroscopic models focus on speed and density to describe traffic flow [
6], while microscopic models focus on individual vehicles and drivers [
7]. They incorporate speed, position, and distance and time headway [
8,
9].
Gazis, Herman, and Rothery proposed a microscopic model commonly known as the GHR model [
10]. This model characterizes driver response considering the speed and distance of leading vehicles. However, driver behavior in changing conditions is ignored as speed adjustments are based on a constant and not traffic physics. Newell [
11] characterized vehicle behavior in dense traffic and showed that velocity (speed) is impacted by the distance headway. An increase in this headway results in higher speeds and lower density. However, high speeds can produce large acceleration, which is neither safe nor realistic.
Wiedemann [
12] and Fritzsche [
13] developed similar models based on driver behavior under varying conditions. Their results are employed in the PTV VISSIM and PARAMICS simulators, respectively [
14]. However, the traffic states have different equations so their models are complex. Wiedemann [
15] created an improved model using simulation results for traffic on motorways. However, this model is not stable for a large number of vehicles [
16].
An improvement to the Newell model was given in [
17], but it neglects speed differences, resulting in acceleration which is very high. Moreover, driver behavior is based on a constant and so traffic physics is ignored. It was shown in [
18] that speed differences can be used to accurately characterize speed and time headway in dense traffic. However, average and slow driver behavior are not considered so the results only pertain to aggressive drivers. The model in [
19] is widely used because it produces realistic traffic behavior. As a consequence, it is employed in the AIMSUN simulator [
14]. However, this model cannot differentiate between acceleration and deceleration and is limited to a small range of parameters [
20].
The Intelligent Driver (ID) model was developed in [
20] based on driver reaction. This model considers desired velocity (speed) and distance headway to characterize driver behavior [
21,
22,
23]. Unlike the Gipps model, the ID model provides realistic acceleration and deceleration [
14]. As a consequence, it is widely utilized in Adaptive Cruise Control (ACC) and cooperative ACC [
24,
25,
26]. The ID model is also employed in Simulation of Urban MObility (SUMO) and PTV VISSIM [
27]. However, it uses a fixed exponent to characterize traffic. This means that driver behavior is not based on traffic conditions. This is unrealistic as real-world traffic dynamics are influenced by various factors including pavement condition, and this affects driver behavior.
This study introduces a microscopic traffic model that incorporates the pavement condition to accurately represent traffic behavior. The pavement condition is evaluated using the Pavement Condition Index (PCI), which is an indicator of pavement condition and quality, and thus affects driver behavior and traffic flow. It ranges between
and
[
2]. Incorporating the PCI results in a model that provides a more comprehensive and accurate representation of traffic behavior. A flowchart of the methodology employed in this research is given in
Figure 1. First, field experiments to determine the impact of vehicle vibrations on the PCI were conducted on the Grand Trunk highway in Peshawar, located in the Khyber Pakhtunkhwa province of Pakistan. This road section spans 7 km and extends from the Chamkani Bus Rapid Transit (BRT) station to Pabbi. Then, the proposed and ID models are implemented using the Euler technique in MATLAB. The results obtained indicate that the proposed model is more suitable for evaluating traffic behavior.
The rest of this paper is organized as follows. In
Section 2, traffic flow models are introduced and their stability is analyzed in
Section 3.
Section 4 outlines the Euler technique and the performance is evaluated in
Section 5. The results of this paper are summarized in
Section 6.
2. Traffic Models
The ID model is used for microscopic traffic characterization and incorporates factors such as the desired speed
, distance to align with leading conditions
, and the difference in speed
with the leading vehicle [
20]. Driver response is a function of the ratio of average speed
to desired speed
, and is expressed as [
20]
where
is the maximum acceleration and
is a fixed acceleration exponent.
is the desired distance headway during traffic alignment to leading conditions and is given by [
20]
where
is the deceleration or minimum acceleration,
is the jam spacing as illustrated in
Figure 2, and
is the time required by a vehicle to adjust its speed to the speed of the leading vehicle [
5].
indicates driver desire to maintain a safe distance from the leading vehicle. This is crucial for ensuring safety on the road and preventing collisions. The ID model employs (1) and (2) for traffic by incorporating driver response and distance headway for the alignment of traffic [
5].
The ID model characterizes driver response to traffic conditions based on a fixed value . Thus, driver behavior does not vary based on these conditions, so it is unrelated to traffic physics and results in inadequate and unrealistic traffic characterization.
An acceleration exponent based on the PCI is proposed for the realistic characterization of traffic. Then, is a function of vehicle vibrations which are mechanical oscillations. These vibrations are largely generated by the interaction between the road surface and tires, and thus are a major contributor to passenger fatigue and discomfort.
Field experiments were conducted by driving a test vehicle over the road segment in Peshawar, Pakistan, between 12 AM and 2 AM. One lane in each direction was traversed
times with speeds of
km/h (
m/s), 45 km/h (
m/s), and 55 km/h (
m/s). Thus, for a given speed, a lane was traversed four times. These speeds were selected to represent typical traffic observed on the road segment. Data were collected using an On-Board Diagnostic-II scanner connected to a smartphone with the BotlnckDectr [
28] mobile app. This allowed for the recording of various parameters including GPS location, in-vehicle noise, vibration, and time [
28]. During the experiments, the smartphone was positioned on the vehicle dashboard. The data were transmitted to the Amazon Web Services (AWS) cloud. It was then analyzed to obtain the
PCI of the road segment. The relationships between PCI and vehicle vibrations obtained are
for speeds of approximately
m/s,
m/s, and
m/s, respectively. The PCI ranges from
to
where
corresponds to a poor road condition and
to an excellent road condition. Thus,
and PCI are linearly related. As the pavement condition degrades, the oscillations increase, which reduces passenger comfort, i.e., a higher PCI corresponds to lower vibrations. Substituting (3), (4) and (5) in (1) gives the proposed model for speeds of
m/s,
m/s, and
m/s, respectively
An excellent road condition is required to avoid traffic congestion and accidents and efficiently align to forward vehicles. In this case, there is free flow traffic which corresponds to . A poor road condition can result in congestion due to the reduction in vehicle speed. In this case, and vehicle acceleration and deceleration are large so the emissions are high. With the proposed model, alignment is according to the PCI and is more realistic compared with fixed .
The traffic density can be expressed as
[
29] where
is the distance headway at equilibrium. In this case,
so substituting (2) in (1) gives for the ID model
and rearranging we obtain
Thus, the fixed
in the ID model results in a constant distance headway between vehicles at equilibrium regardless of the traffic conditions. In contrast, in the proposed model the distance headway is based on the PCI. The distance headway at equilibrium is obtained by substituting (3), (4), and (5) in (10) which gives
The product of density and speed is traffic flow [
8,
30] so that
and substituting (10) in (14) gives the flow for the ID model as
This is unrealistic as it relies on a fixed exponent. The proposed model considers the PCI to determine traffic flow and so is more realistic. The traffic flow can be expressed as
The proposed model indicates that when the road condition is poor, the vehicle vibrations are large and the flow is small, whereas when the road condition is excellent, the vehicle vibrations are small and the flow is large. Further, the proposed model can predict traffic behavior in real-time to help ACC systems better anticipate and adapt to changes in traffic conditions. An ACC system guided by the proposed model can adjust the vehicle speed and following distance in response to the observed traffic density. When the density is high, the ACC can reduce the speed and maintain a safe distance to ensure safety and a smooth traffic flow. Conversely, when the density is low, the ACC can increase the speed while maintaining a safe distance to improve efficiency.
3. Stability Analysis
This section presents an analysis of the stability of traffic models considering an infinitely long road. Identical vehicles are assumed with a constant equilibrium distance headway [
31]. Therefore, drivers adjust to forward conditions with minimal acceleration, so there are only small changes in the equilibrium velocity
associated with
. The corresponding change in distance headway, denoted by
, is also small as is the change in velocity denoted by
. The distance headway can then be expressed as
and
The temporal change in velocity during traffic alignment over the distance headway is [
32]
where the subscripts
and
denote the following and leading vehicles, respectively. Given the minor variations in
, the adjustments in headway are negligible. Consequently,
during alignment can be expressed as [
31]
where
,
, and
denote the partial derivatives w.r.t. velocity, change in velocity, and distance headway, respectively, which are
Using Fourier–Ansatz to express (21) and (22) gives
so (23) and (24) can be written as
where
corresponds to the traffic oscillations during alignment and
. The real part
corresponds to the amplitude change and
is the oscillation frequency with oscillation period
. The parameter
denotes driver delay [
31], while
and
are the changes in velocity and distance headway, respectively.
Substituting (25) in (21) and (22), gives
Model stability requires that the real components of the eigenvalues are negative. The eigenvalues are the solution of
The Jacobian matrix is
where
and
are the gradients of (26) and (27) w.r.t.
and
and
are the gradients of (26) and (27) w.r.t.
. We have
and substituting this in (28) gives
so that
Setting
and
, (24) becomes
Thus, the eigenvalues from (32) are
A model is string stable [
31] if the real components of the eigenvalues are negative. Under this condition, traffic oscillations diminish over time and the flow becomes stable and smooth [
33]. Conversely, a model is considered unstable if traffic oscillations increase and are large as in congestion. In this case, acceleration is high unlike when there is string stability [
29]. As the model becomes unstable,
, leading to minimal delay between flow changes (traffic waves) [
31].
Approximating
and
using Taylor series for a small delay, i.e.,
, gives
From [
31], at equilibrium
where
is the equilibrium speed gradient relative to the distance headway. Then, (35) becomes
Let
where
Considering a Taylor series expansion, the square root in (33) can be approximated as
which gives
Using (38)
and then from (39), we obtain
The real part of (43) represents the rate at which the traffic oscillation amplitude changes, signifying growth or decay. When this real part is negative, the traffic flow is string-stable, since
Then,
is the string stability criterion [
27] which can be expressed as
From (44) and (45), the product of
and
indicates that
has a negative real part. Further, at equilibrium
Using (46) and (47), the criterion for string stability from (44) is
Thus, the velocity with the ID model is determined by
. A higher value improves stability but may lead to optimistic performance in congestion. Consequently, increasing
for stability reasons ignores traffic physics and can produce unrealistic results [
5]. Changes in velocity during traffic alignment are influenced by driver response and thus pavement condition. Hence, more realistic behavior is obtained using (3), (4), and (5) for
according to the speed. The stability criteria for the proposed model with speeds
m/s,
m/s, and
m/s are then
respectively. When the pavement is in good condition, vehicles can more easily adjust to changes in traffic ensuring string stability. Conversely, pavement in poor condition results in greater adjustments to changes in traffic which may not result in a smooth flow.
5. Performance Evaluation
In this section, the performance of the proposed model and ID models is evaluated on a circular road of length
m. The Euler scheme is employed with time step
s. The proposed model is simulated for
s and the ID model for
s. Based on (3), (4), and (5) the desired speed
for the proposed model is set to
m/s,
m/s, and
m/s. The desired speed for the ID model is
m/s [
22]. The jam spacing is set to
m [
31], the maximum acceleration is
m/s², and the minimum acceleration is
m/s² [
20]. The acceleration exponent
is typically
or greater and is often set to 4 [
20]. Thus, here
and
. The PCI values considered are
and
. The maximum normalized density is set to
and the critical density is
[
36]. The maximum flow is obtained at the critical density with speed
. Thus, the speed is normalized by
and the flow is normalized by
. The simulation parameters are summarized in
Table 1.
Figure 3 gives the normalized flow for the proposed model with
and
and
. When
the flow at
s is
. It is
at
s, increasing to
at
s and
at
s. When
, the flow at
s is
. It is
at
s, increasing to
at
s and
at
s. When
, the flow at
s is
. It is
at
s, increasing to
at
s and
at
s.
Figure 4 gives the normalized flow for the proposed model with
m/s and
and
. When
, the flow at
s is
. It is
at
s, increasing to
at
s and
at
s. When
the flow at
s is
. It is
at
s, increasing to
at
s and
at
s. When
, the flow at
s is
. It is
at
s, increasing to
at
s and
at
s.
Figure 5 gives the normalized flow for the proposed model with
m/s and
and
. When
, the flow at
s is
. It is
at
s, increasing to
at
s and
at
s. When
, the flow at
s is
. It is
at
s, increasing to
at
s and
at
s. When
, the flow at
s is
, increasing to
at
s and
at
s.
Figure 6 gives the normalized flow for the ID model with
and
and
m/s. When
, the flow at
s is
, increasing to
at
s and
at
s. When
, the flow at
s is
. It is
at
s, increasing to
at
s and
at
s. When
, at
s the flow is
. It is
at
s, increasing to
at
s and
at
s.
Figure 7 gives the normalized speed with
m/s and
and
for the proposed model. When
, the speed is
from
s to
s, decreasing to
at
s, and then increasing to
at
s. The speed oscillates between
and
from
s to
s as indicated in
Figure 7a. The speed when
is similar to that when
. It is
from
s to
s, decreasing to
at
s, and then increasing to
at
s. The speed oscillates between
and
from
s to
s as indicated in
Figure 7b. When
, the speed is
from
s to
s, decreasing to
at
s, and then increasing to
at
s. The speed oscillates between
and
from
s to
s as indicated in
Figure 7c. For all PCI values, there are road segments where the speed is constant such as between
m and
m at
s when
, between
m and
m at
s when
, and between
m and
m at
s when
.
Figure 8 gives the normalized speed with
m/s and
and
for the proposed model. When
, the speed from
s to
s is
, decreasing to
at
s, and then increasing to
at
s. The speed oscillates between
and
from
s to
s as indicated in
Figure 8a. Similarly, when
the speed is
from
s to
s, decreasing to
at
s, and then increasing to
at
s. The speed oscillates between
and
from
s to
s as indicated in
Figure 8b. The speed is also similar when
. It is
from
s to
s, decreasing to
at
s and then increasing to 0.18 at
s. The speed oscillates between
and
from
s to
s as indicated in
Figure 8c. For all PCI values, there are road segments where the speed is constant such as between
m and
m at
s when
, between
m and
m at
s when
, and between
m and
m at
s when
.
Figure 9 gives the normalized speed with
m/s and
and
for the proposed model. When
, the speed from
s to
s is
, decreasing to
at
s, and then increasing to
at
s. The speed oscillates between
and
from
s to
s as indicated in
Figure 9a. The speed behavior is similar when
. It is
from
s to
s, decreasing to
at
s, and then increasing to
at
s. The speed oscillates between
and
from
s to
s as indicated in
Figure 9b. Similar speed behavior also occurs when
. It is
from
s to
s, decreasing to
at
s, and then increasing to
at
s. The speed oscillates between
and
from
s to
as indicated in
Figure 9c. For all PCI values, there are road segments where the speed is constant such as between
m and
m at
s when
, between
m and
m at
s when
, and between
m and
m at
s when
.
Figure 10 gives the normalized speed for the ID model with
m/s and
and
. When
, the speed is
until
s. It is
at
s and then increases to
at
s. The speed oscillates between
and
from
s to
s as indicated in
Figure 10a. When
, the speed is
until
s, decreasing to
at
s, and then increasing to
at
s. The speed oscillates between
and
from
s to
s as indicated in
Figure 10b. When
, the speed is
until
s, decreasing to
at
s, and then increasing to
at
s. The speed oscillates between
and
from
s to
s as indicated in
Figure 10c. For all the values of
, there are road segments where the speed is constant such as between
m and
m at
s when
, between
m and
m at
s when
, and between
m and
m at
s when
.
Figure 11 gives the normalized density for the proposed model with
m/s and
and
. When
, the density is
until
s. It is
at
s, decreasing to
at
s. The density oscillates between
and
from
s to
s as indicated in
Figure 11a. When
, the density is
until
s. It is
at
s, decreasing to
at
s. The density oscillates between
and
from
s to
s, and then it varies between
and
as indicated in
Figure 11b. When
, the density is
until
s. It is
at
s, decreasing to
at
s. It oscillates between
and
from
s to
s, and then varies between
and
as indicated in
Figure 11c.
Figure 12 gives the normalized density for the proposed model with
m/s and
and
. When
, the density is
until
s, increasing to
at
s and then decreasing to
at
s. The density oscillates between
and
from
s to
s, and then it varies between
and
as indicated in
Figure 12a. When
, the density is
until
s, increasing to
at
s, and then decreasing to
at
s. The density oscillates between
and
from
s to
s, and then it varies between
and
as indicated in
Figure 12b. When
, the density is
until
s, increasing to
at
s, and then decreasing to
at
s. From
s to
s, the density oscillates between
and
, and then it varies between
and
as indicated in
Figure 12c.
Figure 13 gives the normalized density for the proposed model with
m/s and
and
. When
, the density is
until
s, increasing to
at
s, and then decreasing to
at
s. The density oscillates between
and
from
s to
s, and then it varies between
and
as indicated in
Figure 13a. When
the density is at
until
s, increasing to
at
s, and then decreasing to
at
s. The density oscillates between
and
from
s to
s as indicated in
Figure 13b. When
, the density is
until
s, increasing to
at
s, and then decreasing to
at
s. The density oscillates between
and
from
s to
s as indicated in
Figure 13c.
Figure 14 gives the normalized density for the ID model with
m/s and
and
When
, the density is
until
s, increasing to
at
s, and then decreasing to
at
s. The density oscillates between
and
from
s to
s, and then it varies between
and
as indicated in
Figure 14a. When
, the density is
until
s, increasing to
at
s, and then decreasing to
at
s. It oscillates between
and
from
s to
s as indicated in
Figure 14b. When
, the density is
until
s, increasing to
at
s, and then decreasing to
at
s. The density oscillates between
and
from
s to
s and then it varies between
and
as indicated in
Figure 14c.
The results for the proposed model indicate that pavement condition influences traffic flow as expected. In particular, the flow increases with speed as shown in
Figure 3,
Figure 4 and
Figure 5. The flow with the ID model increases with
, which is not based on traffic physics. Furthermore, the oscillations in speed and density with the proposed model vary with the PCI and decrease over time as the PCI increases. These results are more realistic as they are based on real parameters such as the PCI. Conversely, the oscillations in speed and density with the ID model are the result of an arbitrary fixed parameter, and they increase over time as
increases with no justification. This is an inadequate and unrealistic traffic characterization.