Next Article in Journal
Optimizing Multi-Layer Perovskite Solar Cell Dynamic Models with Hysteresis Consideration Using Artificial Rabbits Optimization
Next Article in Special Issue
Optimizing Propellant Distribution for Interorbital Transfers
Previous Article in Journal
Emotional Health and Climate-Change-Related Stressor Extraction from Social Media: A Case Study Using Hurricane Harvey
Previous Article in Special Issue
Drone-Based Decentralized Truck Platooning with UWB Sensing and Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Microscopic Traffic Model Incorporating Vehicle Vibrations Due to Pavement Condition

by
Faryal Ali
1,*,
Zawar Hussain Khan
1,
Khurram Shehzad Khattak
2,
Thomas Aaron Gulliver
1 and
Ahmed B. Altamimi
3
1
Department of Electrical and Computer Engineering, University of Victoria, Victoria, BC V8W 2Y2, Canada
2
Department of Computer Systems Engineering, University of Engineering and Technology, Peshawar 25000, Pakistan
3
College of Computer Science and Engineering, University of Ha’il, Ha’il 55476, Saudi Arabia
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(24), 4911; https://doi.org/10.3390/math11244911
Submission received: 7 September 2023 / Revised: 4 December 2023 / Accepted: 6 December 2023 / Published: 9 December 2023
(This article belongs to the Special Issue Modeling and Simulation in Engineering, 3rd Edition)

Abstract

:
A microscopic traffic flow model is developed that incorporates vehicle vibrations due to pavement condition. The Intelligent Driver (ID) model employs a fixed exponent so traffic behavior is the same regardless of the road condition. Thus, it ignores the underlying physics. To address this limitation, the proposed model employs the Pavement Condition Index (PCI) in describing traffic behavior. The performance of both models is evaluated on a 3000 m circular road using the Euler numerical discretization technique. The results show that the performance of the proposed model varies with the pavement condition (PCI), as expected. Furthermore, the traffic flow increases with vehicle speed. The oscillations in speed and density with the proposed model decrease as the PCI increases, and are larger when the speed is higher. Consequently, the results with the proposed model align more closely with reality as they are based on the PCI, and so are a more accurate representation of traffic behavior.

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 55 % when the road condition is poor compared to when it is excellent, and average emissions increase by 2.49 %. Road safety is a primary concern worldwide as road accidents cause 1.3 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 0 and 100 [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 v d , distance to align with leading conditions s , and the difference in speed v with the leading vehicle [20]. Driver response is a function of the ratio of average speed v to desired speed v d , and is expressed as [20]
d v d t = a 1 v v d δ H s 2
where a is the maximum acceleration and δ is a fixed acceleration exponent. H is the desired distance headway during traffic alignment to leading conditions and is given by [20]
H = J + T v + v v 2 a d
where d is the deceleration or minimum acceleration, J is the jam spacing as illustrated in Figure 2, and T is the time required by a vehicle to adjust its speed to the speed of the leading vehicle [5].   H 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 12 times with speeds of 35 km/h ( 9.72 m/s), 45 km/h ( 12.50 m/s), and 55 km/h ( 15.27 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
δ = 0.0169 P C I + 4.068
δ = 0.0265 P C I + 5.037
δ = 0.0251 P C I + 5.209
for speeds of approximately 9.72 m/s, 12.50 m/s, and 15.27 m/s, respectively. The PCI ranges from 0 to 100 where 0 corresponds to a poor road condition and 100 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 9.72 m/s, 12.50 m/s, and 15.27 m/s, respectively
d v d t = a 1 v v d ( 0.0169 P C I + 4.068 ) H s 2   ( v d = 9.72   m / s )
d v d t = a 1 v v d ( 0.0265 P C I + 5.037 ) H s 2   ( v d = 12.50   m / s )
d v d t = a 1 v v d ( 0.0251 P C I + 5.209 ) H s 2   ( v d = 15.27   m / s )
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 P C I = 100 . A poor road condition can result in congestion due to the reduction in vehicle speed. In this case, P C I = 0 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 D = 1 / s e [29] where s e is the distance headway at equilibrium. In this case, v = 0 so substituting (2) in (1) gives for the ID model
a 1 v v d δ J + T v s e 2   = 0
and rearranging we obtain
s e = J + T v 1 v v d δ 1 2
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
s e = J + T v 1 v v d ( 0.0169 P C I + 4.068 ) 1 2   ( v d = 9.72   m / s )
s e = J + T v 1 v v d ( 0.0265 P C I + 5.037 ) 1 2   ( v d = 12.50   m / s )
s e = J + T v 1 v v d ( 0.0251 P C I + 5.209 ) 1 2   ( v d = 15.27   m / s )
The product of density and speed is traffic flow [8,30] so that
F = v s e
and substituting (10) in (14) gives the flow for the ID model as
F = v J + T v 1 v v d δ 1 2
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
F = v J + T v 1 v v d ( 0.0169 P C I + 4.068 ) 1 2   ( v d = 9.72   m / s )
F = v J + T v 1 v v d ( 0.0265 P C I + 5.037 ) 1 2   ( v d = 12.50   m / s )
F = v J + T v 1 v v d ( 0.0251 P C I + 5.209 ) 1 2   ( v d = 15.27   m / s )
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 v e associated with s e . The corresponding change in distance headway, denoted by a , is also small as is the change in velocity denoted by b . The distance headway can then be expressed as
s = s e + a ,
and
v = v e s e + b .
The temporal change in velocity during traffic alignment over the distance headway is [32]
a t = d a d t = b l b F ,
where the subscripts F and l denote the following and leading vehicles, respectively. Given the minor variations in v e s e , the adjustments in headway are negligible. Consequently, b t during alignment can be expressed as [31]
b t = d b d t = f s a F + ( f v + f v ) b F f v b l ,
where f v , f v   , and f s denote the partial derivatives w.r.t. velocity, change in velocity, and distance headway, respectively, which are
f v = f v ,   f v = f v   and   f s = f s .
Using Fourier–Ansatz to express (21) and (22) gives
a t = a ^ e γ t + i k ,
b t = b ^ e γ t + i k ,
so (23) and (24) can be written as
a t b t = a ^ b ^ e γ t + i k ,
where γ = α + i ω corresponds to the traffic oscillations during alignment and i = 1 . The real part α corresponds to the amplitude change and ω = 2 π T is the oscillation frequency with oscillation period T . The parameter k denotes driver delay [31], while a ^ and b ^ are the changes in velocity and distance headway, respectively.
Substituting (25) in (21) and (22), gives
a t = b b e i k ,
b t = f s a e i k + f v + f v b e i k f v b .
Model stability requires that the real components of the eigenvalues are negative. The eigenvalues are the solution of
j λ 0 0 λ = 0 .
The Jacobian matrix is
j = j 11 j 12 j 21 j 22 ,
where j 11 and j 21 are the gradients of (26) and (27) w.r.t. a and j 12 and j 22 are the gradients of (26) and (27) w.r.t. b . We have
j = e i k 0 e i k 1 f s f v + f v f v e i k .
and substituting this in (28) gives
λ 1 e i k f s λ f v f v + f v e i k = 0 ,
so that
λ 2 + f v f v + f v e i k λ + f s 1 e i k = 0 .
Setting M k = f v f v + f v e i k and N k = f s 1 e i k , (24) becomes
λ 2 + M k λ + N k = 0 .
Thus, the eigenvalues from (32) are
λ 1,2 = M k 2   1 ± 1 4 N k M 2 k ,
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, k 0 , leading to minimal delay between flow changes (traffic waves) [31].
Approximating M ( k ) and N ( k ) using Taylor series for a small delay, i.e., k 0 , gives
M k = f v i f v k ,
N k = i f s k + f s 2 k 2 .
From [31], at equilibrium
f s = v e s e f v .
where v e s e is the equilibrium speed gradient relative to the distance headway. Then, (35) becomes
N k = i v e s e f v v e s e 2 f v .
Let
M k = x 1 + x 2 k , N k = y 1 k + y 2 k 2 ,
where
x 1 = f v , x 2 = i f v , y 1 = i v e s e f v = i v e s e x 1 , y 2 = v e s e 2 f v = v e s e 2 x 1 .
Considering a Taylor series expansion, the square root in (33) can be approximated as
1 4 N k M 2 k = 1 2 N k M 2 k 2 N 2 k M 4 k ,
which gives
λ 2 = N k M 2 k N 2 k M 3 k .
Using (38)
λ 2 = y 1 x 1 k + y 1 x 2 x 1 2 y 2 x 1 y 1 2 x 1 3 k 2 ,
and then from (39), we obtain
λ 2 = i v e s e k + v e s e f v 2 f v f v 2 v e s e k 2 .
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
v e s e   0   a n d   f v < 0 .
Then, 2 f v f v 2 v e s e is the string stability criterion [27] which can be expressed as
v e s e f v 2 f v .
From (44) and (45), the product of 2 f v f v 2 v e s e and v e s e f v indicates that λ 2 has a negative real part. Further, at equilibrium
f v = a δ v e ( s e ) δ 1 v d δ 2 T J + v e ( s e ) T s e 2 ,
f v = v e ( s e ) s e a d J + v e ( s e ) T s e .
Using (46) and (47), the criterion for string stability from (44) is
v e s e a δ s e 2 v e s e δ 1 + 2 T J v d δ + 2 v e s e T 2 v d δ 2 s e 2 v d δ + v e s e a d s e + T v e s e s e 2 d
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 9.72 m/s, 12.50 m/s, and 15.27 m/s are then
v e s e a 0.0169 P C I + 4.068 s e 2 v e s e 0.0169 P C I + 4.068 1 + 2 T J v d 0.0169 P C I + 4.068 + 2 v e s e T 2 v d ( 0.0169 P C I + 4.068 ) 2 s e 2 v d ( 0.0169 P C I + 4.0679 ) + v e s e a d s e + T v e s e s e 2 d
v e s e a 0.0265 P C I + 5.037 s e 2 v e s e 0.0265 P C I + 5.037 1 + 2 T J v d 0.0265 P C I + 5.037 + 2 v e s e T 2 v d ( 0.0265 P C I + 5.037 ) 2 s e 2 v d ( 0.0265 P C I + 5.037 ) + v e s e a d s e + T v e s e s e 2 d
v e s e a 0.0251 P C I + 5.209 s e 2 v e s e 0.0251 P C I + 5.209 1 + 2 T J v d 0.0251 P C I + 5.209 + 2 v e s e T 2 v d 0.0251 P C I + 5.209 2 s e 2 v d δ + v e s e a d s e + T v e s e s e 2 d
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.

4. The Euler Technique

The Euler technique is used to evaluate the proposed and ID models. It is a simple but effective method to solve systems of differential equations and is widely used in traffic simulators such as SUMO [34] and AIMSUN [35]. This technique divides time into discrete steps and the vehicle position, speed, and acceleration are approximated using the model at each time step. The change in distance w.r.t. time results in a change in speed given by
d s d t = v ,
and temporal changes in speed lead to changes in acceleration. Denote the right-hand side of (1), (6), (7), and (8) by Υ , then
d v d t = Υ .
For the Euler technique, the position and speed for the ID and proposed models is
s f x + 1 = s f x + Δ t × v f x
v f x + 1 = v f x + Δ t × Υ f x
where x is the current time step and x + 1 is the next time step. s f x , v f x , and Υ f x are the position, speed, and acceleration, respectively, of the following vehicle in the xth time interval where
t = x Δ t
and Δ t is the duration of a time step.

5. Performance Evaluation

In this section, the performance of the proposed model and ID models is evaluated on a circular road of length 3000 m. The Euler scheme is employed with time step Δ t = 0.50 s. The proposed model is simulated for 400 s and the ID model for 150 s. Based on (3), (4), and (5) the desired speed v d for the proposed model is set to 9.72 m/s, 12.50 m/s, and 15.27 m/s. The desired speed for the ID model is 20 m/s [22]. The jam spacing is set to 2.0 m [31], the maximum acceleration is 0.73 m/s², and the minimum acceleration is 1.67 m/s² [20]. The acceleration exponent δ is typically 1 or greater and is often set to 4 [20]. Thus, here δ = 1 ,   4 and 20 . The PCI values considered are P C I = 0 ,   50 and 100 . The maximum normalized density is set to 1 / J = 0.50 and the critical density is 0.25 [36]. The maximum flow is obtained at the critical density with speed v d . Thus, the speed is normalized by v d and the flow is normalized by 0.25 × v d . The simulation parameters are summarized in Table 1.
Figure 3 gives the normalized flow for the proposed model with v d = 9.72 and P C I = 0 ,   50 , and 100 . When P C I = 0 , the flow at 19.0 s is 0.0010 . It is 0.0021 at 106.5 s, increasing to 0.0032 at 293.0 s and 0.0035 at 400 s. When P C I = 50 , the flow at 19.5 s is 0.0010 . It is 0.0020 at 162.5 s, increasing to 0.0031 at 288.0 s and 0.0038 at 400 s. When P C I = 100 , the flow at 21.5 s is 0.0010 . It is 0.0025 at 172.0 s, increasing to 0.0033 at 285.5 s and 0.0044 at 400 s.
Figure 4 gives the normalized flow for the proposed model with v d = 12.50 m/s and P C I = 0 ,   50 , and 100 . When P C I = 0 , the flow at 20.0 s is 0.0012 . It is 0.0029 at 156.0 s, increasing to 0.0042 at 268.5 s and 0.0054 at 400 s. When P C I = 50 , the flow at 22.5 s is 0.0013 . It is 0.0033 at 167.0 s, increasing to 0.0045 at 276.0 s and 0.0063 at 400 s. When P C I = 100 , the flow at 23.5 s is 0.0010 . It is 0.0037 at 179.0 s, increasing to 0.0062 at 313.5 s and 0.0110 at 400 s.
Figure 5 gives the normalized flow for the proposed model with v d = 15.27 m/s and P C I = 0 ,   50 , and 100 . When P C I = 0 , the flow at 22.0 s is 0.0013 . It is 0.0050 at 217.5 s, increasing to 0.0079 at 307.5 s and 0.0140 at 400 s. When P C I = 50 , the flow at 25.0 s is 0.0013 . It is 0.0052 at 220.5 s, increasing to 0.0110 at 324.0 s and 0.0270 at 400 s. When P C I = 100 , the flow at 27.5 s is 0.0010 , increasing to 0.0070 at 234.5 s and 0.0770 at 400 s.
Figure 6 gives the normalized flow for the ID model with δ = 1 ,   4 , and 20 and v d = 20 m/s. When δ = 1 , the flow at 31.5 s is 0.0013 , increasing to 0.0024 at 94.0 s and 0.0063 at 150 s. When δ = 4 , the flow at 33.0 s is 0.0017 . It is 0.0020 at 80.0 s, increasing to 0.0039 at 116.0 s and 0.0088 at 150 s. When δ = 20 , at 28.0 s the flow is 0.0012 . It is 0.0032 at 100.5 s, increasing to 0.0060 at 130.0 s and 0.0095 at 150 s.
Figure 7 gives the normalized speed with v d = 9.72 m/s and P C I = 0 ,   50 and 100 for the proposed model. When P C I = 0 , the speed is 0.39 from 0.5 s to 15.0 s, decreasing to 0.23 at 15.5 s, and then increasing to 0.46 at 20.0 s. The speed oscillates between 0.16 and 0.64 from 236.5 s to 399.0 s as indicated in Figure 7a. The speed when P C I = 50 is similar to that when P C I = 0 . It is 0.39 from 0.5 s to 15.0 s, decreasing to 0.23 at 15.5 s, and then increasing to 0.46 at 21.0 s. The speed oscillates between 0.21 and 0.58 from 263.0 s to 399.5 s as indicated in Figure 7b. When P C I = 100 , the speed is 0.40 from 0.5 s to 15.0 s, decreasing to 0.23 at 15.5 s, and then increasing to 0.46 at 19.5 s. The speed oscillates between 0.32 and 0.48 from 315.0 s to 399.0 s as indicated in Figure 7c. For all PCI values, there are road segments where the speed is constant such as between 1947.8 m and 426.7 m at 397.5 s when P C I = 0 , between 1932.3 m and 228.4 m at 392.0 s when P C I = 50 , and between 1962.4 m and 428.1 m at 393.5 s when P C I = 100 .
Figure 8 gives the normalized speed with v d = 12.50 m/s and P C I = 0 ,   50 , and 100 for the proposed model. When P C I = 0 , the speed from 0.5 s to 15.0 s is 0.15 , decreasing to 0.09 at 15.5 s, and then increasing to 0.18 at 19.5 s. The speed oscillates between 0.06 and 0.27 from 255.0 s to 399.0 s as indicated in Figure 8a. Similarly, when P C I = 50 the speed is 0.15 from 0.5 s to 15.0 s, decreasing to 0.09 at 15.5 s, and then increasing to 0.18 at 20.5 s. The speed oscillates between 0.06 and 0.26 from 258.0 s to 399.0 s as indicated in Figure 8b. The speed is also similar when P C I = 100 . It is 0.15 from 0.5 s to 15.0 s, decreasing to 0.09 at 15.5 s and then increasing to 0.18 at 19.5 s. The speed oscillates between 0.09 and 0.22 from 294.5 s to 399.5 s as indicated in Figure 8c. For all PCI values, there are road segments where the speed is constant such as between 2022.5 m and 237.5 m at 397.0 s when P C I = 0 , between 1922.7 m and 245.5 m at 395.0 s when P C I = 50 , and between 1953.75 m and 230.0 m at 390.0 s when P C I = 100 .
Figure 9 gives the normalized speed with v d = 15.27 m/s and P C I = 0 ,   50 , and 100 for the proposed model. When P C I = 0 , the speed from 0.5 s to 15.0 s is 0.11 , decreasing to 0.06 at 15.5 s, and then increasing to 0.13 at 19.5 s. The speed oscillates between 0.04 and 0.20 from 241.0 s to 399.0 s as indicated in Figure 9a. The speed behavior is similar when P C I = 50 . It is 0.11 from 0.5 s to 15.0 s, decreasing to 0.06 at 15.5 s, and then increasing to 0.14 at 20.5 s. The speed oscillates between 0.04 and 0.20 from 267.5 s to 399.0 s as indicated in Figure 9b. Similar speed behavior also occurs when P C I = 100 . It is 0.11 from 0.5 s to 15.0 s, decreasing to 0.06 at 15.5 s, and then increasing to 0.13 at 20.5 s. The speed oscillates between 0.05 and 0.18 from 278.5 s to 399.5 as indicated in Figure 9c. For all PCI values, there are road segments where the speed is constant such as between 1945.4 m and 244.3 m at 391.0 s when P C I = 0 , between 1911.8 m and 239.7 m at 388.5 s when P C I = 50 , and between 1817.1 m and 371.4 m at 384.0 s when P C I = 100 .
Figure 10 gives the normalized speed for the ID model with v d = 15.27 m/s and δ = 1 , 4 , and 20 . When δ = 1 , the speed is 0.10 until 15.0 s. It is 0.05 at 15.5 s and then increases to 0.11 at 21.0 s. The speed oscillates between 0.09 and 0.10 from 118.0 s to 149.5 s as indicated in Figure 10a. When δ = 4 , the speed is 0.09 until 15.0 s, decreasing to 0.05 at 15.5 s, and then increasing to 0.11 at 19.5 s. The speed oscillates between 0.08 and 0.11 from 122.5 s to 149.5 s as indicated in Figure 10b. When δ = 20 , the speed is 0.1 until 15.0 s, decreasing to 0.05 at 15.5 s, and then increasing to 0.11 at 20.5 s. The speed oscillates between 0.08 and 0.11 from 118.5 s to 150.0 s as indicated in Figure 10c. For all the values of δ , there are road segments where the speed is constant such as between 2590.0 m and 322.0 m at 148.0 s when δ = 1 , between 2498.0 m and 292.0 m at 147.0 s when δ = 4 , and between 2594.0 m and 374.0 m at 149.0 s when δ = 20 .
Figure 11 gives the normalized density for the proposed model with v d = 9.72 m/s and P C I = 0 ,   50 and 100 . When P C I = 0 , the density is 0.17 until 31.0 s. It is 0.20 at 32.0 s, decreasing to 0.13 at 32.5 s. The density oscillates between 0.13 and 0.21 from 279.5 s to 398.5 s as indicated in Figure 11a. When P C I = 50 , the density is 0.16 until 31.0 s. It is 0.20 at 32.0 s, decreasing to 0.13 at 33.0 s. The density oscillates between 0.14 and 0.19 from 277.5 s to 399.5 s, and then it varies between 0.16 and 0.23 as indicated in Figure 11b. When P C I = 100 , the density is 0.16 until 31.0 s. It is 0.20 at 32.0 s, decreasing to 0.13 at 32.5 s. It oscillates between 0.14 and 0.18 from 313.5 s to 398.5 s, and then varies between 0.16 and 0.18 as indicated in Figure 11c.
Figure 12 gives the normalized density for the proposed model with v d = 12.50 m/s and P C I = 0 ,   50 and 100 . When P C I = 0 , the density is 0.17 until 31.0 s, increasing to 0.20 at 32.0 s and then decreasing to 0.13 at 32.5 s. The density oscillates between 0.13 and 0.21 from 271.5 s to 398.0 s, and then it varies between 0.17 and 0.28 as indicated in Figure 12a. When P C I = 50 , the density is 0.17 until 31.0 s, increasing to 0.20 at 32.0 s, and then decreasing to 0.13 at 32.5 s. The density oscillates between 0.13 and 0.23 from 269.0 s to 399.0 s, and then it varies between 0.17 and 0.27 as indicated in Figure 12b. When P C I = 100 , the density is 0.16 until 31.0 s, increasing to 0.20 at 31.5 s, and then decreasing to 0.13 at 32.5 s. From 275.5 s to 398.0 s, the density oscillates between 0.14 and 0.19 , and then it varies between 0.16 and 0.21 as indicated in Figure 12c.
Figure 13 gives the normalized density for the proposed model with v d = 15.27 m/s and P C I = 0 ,   50 and 100 . When P C I = 0 , the density is 0.17 until 31.0 s, increasing to 0.20 at 32.0 s, and then decreasing to 0.13 at 33.0 s. The density oscillates between 0.13 and 0.21 from 266.5 s to 398.0 s, and then it varies between 0.17 and 0.28 as indicated in Figure 13a. When P C I = 50 , the density is at 0.17 until 31.0 s, increasing to 0.20 at 32.0 s, and then decreasing to 0.13 at 33.0 s. The density oscillates between 0.13 and 0.22 from 269.5 s to 398.0 s as indicated in Figure 13b. When P C I = 100 , the density is 0.16 until 31.0 s, increasing to 0.20 at 32.0 s, and then decreasing to 0.13 at 33.0 s. The density oscillates between 0.13 and 0.20 from 288.5 s to 398.5 s as indicated in Figure 13c.
Figure 14 gives the normalized density for the ID model with v d = 15.27 m/s and δ = 1 , 4 , and 20 . When δ = 1 , the density is 0.15 until 31.0 s, increasing to 0.19 at 32.0 s, and then decreasing to 0.12 at 32.5 s. The density oscillates between 0.14 and 0.16 from 123.0 s to 149.0 s, and then it varies between 0.15 and 0.16 as indicated in Figure 14a. When δ = 4 , the density is 0.16 until 31.0 s, increasing to 0.20 at 32.0 s, and then decreasing to 0.13 at 32.5 s. It oscillates between 0.15 and 0.17 from 102.5 s to 149.0 s as indicated in Figure 14b. When δ = 20 , the density is 0.16 until 31.0 s, increasing to 0.20 at 32.0 s, and then decreasing to 0.13 at 33.5 s. The density oscillates between 0.15 and 0.17 from 96.5 s to 149.0 s and then it varies between 0.16 and 0.18 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.

6. Conclusions

A microscopic traffic flow model was developed based on pavement condition. The Pavement Condition Index (PCI) was used to characterize traffic behavior. The performance of the proposed model was evaluated and compared with that of the Intelligent Driver (ID) model. The results obtained demonstrate that the proposed model provides realistic traffic flow dynamics. In particular, the traffic flow under excellent pavement conditions ( P C I = 100 ) is high while the flow under poor pavement conditions ( P C I = 0 ) is low, as expected. Conversely, the ID model has a fixed acceleration exponent which does not reflect the relationship between flow and road condition. Furthermore, the oscillations in speed and density with the proposed model vary according to the pavement condition. They are negligible when the PCI is high, which is expected traffic behavior. In contrast, the ID model produces unrealistic speed and density oscillations based on δ. The results given indicate that the proposed model can be used in traffic simulators for realistic and effective traffic prediction.
The proposed model is a deterministic rather than a probabilistic system. Future research can integrate random variables to provide a probabilistic framework. This will allow the model to deal with the uncertainties and variability in complex traffic environments. Furthermore, it can be implemented for road networks to examine challenging traffic situations and propose solutions. Future research can also consider additional parameters. While PCI is a key factor in traffic flow dynamics, it is important to include other factors such as road emergencies to increase the applicability and improve the accuracy and effectiveness in real-world scenarios. This will contribute to the development of more comprehensive and robust models for traffic flow analysis and management.

Author Contributions

Conceptualization, F.A. and Z.H.K.; methodology, F.A.; software, F.A.; validation, F.A., Z.H.K., K.S.K., A.B.A. and T.A.G.; formal analysis, F.A., Z.H.K., K.S.K. and T.A.G.; investigation, F.A., Z.H.K., K.S.K. and T.A.G.; writing—original draft, F.A.; writing—review and editing, F.A., Z.H.K., K.S.K., A.B.A. and T.A.G.; visualization, F.A., Z.H.K., K.S.K., A.B.A. and T.A.G.; funding acquisition, A.B.A., Z.H.K., K.S.K. and T.A.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Scientific Research Deanship at the University of Ha’il, Ha’il, Saudi Arabia through project RG23-186.

Data Availability Statement

No data was generated in this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Setyawan, A.; Kusdiantoro, I.; Syafi’i. The effect of pavement condition on vehicle speeds and motor vehicles emissions. Procedia Eng. 2015, 125, 424–430. [Google Scholar] [CrossRef]
  2. Kumar, R.; Suman, S.K.; Prakash, G. Evaluation of pavement condition index using artificial neural network approach. Transp. Dev. Econ. 2021, 7, 20. [Google Scholar] [CrossRef]
  3. Road Traffic Injuries, World Health Organization. Available online: https://www.who.int/news-room/fact-sheets/detail/road-traffic-injuries (accessed on 23 September 2022).
  4. Accidents Caused by Poor Road Quality and Conditions in DC, Gilman & Bedigian. Available online: https://www.gilmanbedigian.com/accidents-caused-by-poor-road-quality-and-conditions-in-dc/ (accessed on 23 September 2022).
  5. Ali, F.; Khan, Z.H.; Khan, F.A.; Khattak, K.S.; Gulliver, T.A. A new driver model based on driver response. Appl. Sci. 2022, 12, 5390. [Google Scholar] [CrossRef]
  6. Khan, Z.H. Traffic Modelling for Intelligent Transportation Systems. Ph.D. Dissertation, University of Victoria, Victoria, BC, Canada, 2016. [Google Scholar]
  7. Henein, C.M.; White, T. Microscopic information processing and communication in crowd dynamics. Phys. A Stat. Mech. Appl. 2010, 389, 4636–4653. [Google Scholar] [CrossRef]
  8. Ali, F.; Khan, Z.H.; Khattak, K.S.; Gulliver, T.A.; Khan, A.N. A microscopic heterogeneous traffic flow model considering distance headway. Mathematics 2023, 11, 184. [Google Scholar] [CrossRef]
  9. Khan, Z.H.; Gulliver, T.A. A macroscopic traffic model based on anticipation. Arab. J. Sci. Eng. 2019, 44, 5151–5163. [Google Scholar] [CrossRef]
  10. Gazis, D.C.; Herman, R.; Rothery, R.W. Nonlinear follow-the-leader models of traffic flow. Oper. Res. 1961, 9, 545–567. [Google Scholar] [CrossRef]
  11. Newell, G.F. Nonlinear effects in the dynamics of car following. Oper. Res. 1961, 9, 209–229. [Google Scholar] [CrossRef]
  12. Wiedemann, R. Simulation des Straßenverkehrsflusses. Ph.D. Dissertation, Instituts für Verkehrswesen, Universität Karlsruhe, Karlsruhe, Germany, 1974. [Google Scholar]
  13. Fritzsche, H.T.; Ag, D. A model for traffic simulation. Traffic Eng. Control 1994, 35, 317–321. [Google Scholar]
  14. Askari, A.; Farias, D.A.; Kurzhanskiy, A.A.; Varaiya, P. Measuring impact of adaptive and cooperative adaptive cruise control on throughput of signalized intersections. arXiv 2017, arXiv:1611.08973. [Google Scholar]
  15. Wiedemann, R. Modelling of RTI-Elements on multi-lane roads. In Proceedings of the DRIVE Conference, Brussels, Belgium, 4–6 February 1991. [Google Scholar]
  16. Kovács, T.; Bolla, K.; Gil, R.A.; Fábián, C.; Kovács, L. Parameters of the intelligent driver model in signalized intersections. Teh. Vjesn. 2016, 23, 1469–1474. [Google Scholar]
  17. Bando, M.; Hasebe, K.; Nakayama, A.; Shibata, A.; Sugiyama, Y. Dynamical model of traffic congestion and numerical simulation. Phys. Rev. E 1995, 51, 1035–1042. [Google Scholar] [CrossRef] [PubMed]
  18. Helbing, D.; Tilch, B. Generalized force model of traffic dynamics. Phys. Rev. E 1998, 58, 133–138. [Google Scholar] [CrossRef]
  19. Gipps, P.G. A behavioural car-following model for computer simulation. Transp. Res. Part B 1981, 15, 105–111. [Google Scholar] [CrossRef]
  20. Treiber, M.; Hennecke, A.; Helbing, D. Congested traffic states in empirical observations and microscopic simulations. Phys. Rev. E 2000, 62, 1805–1824. [Google Scholar] [CrossRef]
  21. Cao, Z.; Lu, L.; Chen, C.; Chen, X.U. Modeling and simulating urban traffic flow mixed with regular and connected vehicles. IEEE Access 2021, 9, 10392–10399. [Google Scholar] [CrossRef]
  22. Dahui, W.; Ziqiang, W.; Ying, F. Hysteresis phenomena of the intelligent driver model for traffic flow. Phys. Rev. E 2007, 76, 2–8. [Google Scholar] [CrossRef] [PubMed]
  23. Kesting, A.; Treiber, M.; Helbing, D. Enhanced intelligent driver model to access the impact of driving strategies on traffic capacity. Philos. Trans. R. Soc. Math. Phys. Eng. Sci. 2010, 368, 4585–4605. [Google Scholar] [CrossRef]
  24. Schakel, W.J.; Van Arem, B.; Netten, B.D. Effects of cooperative adaptive cruise control on traffic flow stability. In Proceedings of the International IEEE Conference on Intelligent Transportation Systems, Funchal, Portugal, 19–22 September 2010; pp. 759–764. [Google Scholar]
  25. Milanés, V.; Shladover, S.E.; Spring, J.; Nowakowski, C.; Kawazoe, H.; Nakamura, M. Cooperative adaptive cruise control in real traffic situations. IEEE Trans. Intell. Transp. Syst. 2013, 15, 296–305. [Google Scholar] [CrossRef]
  26. Wang, M.; Treiber, M.; Daamen, W.; Hoogendoorn, S.P.; van Arem, B. Modelling supported driving as an optimal control cycle: Framework and model characteristics. Procedia-Soc. Behav. Sci. 2013, 80, 491–511. [Google Scholar] [CrossRef]
  27. Derbel, O.; Peter, T.; Zebiri, H.; Mourllion, B.; Basset, M. Modified intelligent driver model for driver safety and traffic stability improvement. IFAC Proc. 2013, 46, 744–749. [Google Scholar] [CrossRef]
  28. Sohail, A.M.; Khattak, K.S.; Iqbal, A.; Khan, Z.H.; Ahmad, A. Cloud-based detection of road bottlenecks using OBD-II telematics. In Proceedings of the International Multitopic Conference, Islamabad, Pakistan, 29–30 November 2019. [Google Scholar]
  29. Kessels, F. Traffic Flow Modelling: Introduction to Traffic Flow Theory through a Genealogy of Models; Springer: Cham, Switzerland, 2019. [Google Scholar]
  30. Ali, F.; Khan, Z.H.; Khattak, K.S.; Gulliver, T.A. A microscopic traffic flow model characterization for weather conditions. Appl. Sci. 2022, 12, 12981. [Google Scholar] [CrossRef]
  31. Treiber, M.; Kesting, A. Traffic Flow Dynamics: Data, Models and Simulation; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  32. Hallerbach, S.; Xia, Y. Simulation-based identification of critical scenarios for cooperative and automated vehicles. SAE Int. J. Connect. Autom. Veh. 2018, 1, 93–106. [Google Scholar] [CrossRef]
  33. Feng, S.; Zhang, Y.; Li, S.E.; Cao, Z.; Liu, H.X.; Li, L. String stability for vehicular platoon control: Definitions and analysis methods. Ann. Rev. Control 2019, 47, 81–97. [Google Scholar] [CrossRef]
  34. Krajzewicz, D.; Erdmann, J.; Behrisch, M.; Bieker, L. Recent development and applications of SUMO-simulation of urban mobility. Int. J. Adv. Syst. Meas. 2012, 5, 128–243. [Google Scholar]
  35. Casas, J.; Ferrer, J.L.; Garcia, D.; Perarnau, J.; Torday, A. Traffic simulation with Aimsun. In Fundamentals of Traffic Simulation; International Series in Operations Research & Management Science, 145; Springer: New York, NY, USA, 2010; pp. 173–232. [Google Scholar]
  36. Ahmed, A.; Ngoduy, D.; Adnan, M.; Baig, M.A.U. On the fundamental diagram and driving behavior modeling of heterogeneous traffic flow using UAV-based data. Transp. Res. Part A Policy Pract. 2021, 148, 100–115. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the methodology.
Figure 1. Flowchart of the methodology.
Mathematics 11 04911 g001
Figure 2. Intelligent Driver (ID) model parameters.
Figure 2. Intelligent Driver (ID) model parameters.
Mathematics 11 04911 g002
Figure 3. Normalized flow for the proposed model with v d = 9.72 m/s over a 3000 m circular road.
Figure 3. Normalized flow for the proposed model with v d = 9.72 m/s over a 3000 m circular road.
Mathematics 11 04911 g003
Figure 4. Normalized flow for the proposed model with v d = 12.50 m/s over a 3000 m circular road.
Figure 4. Normalized flow for the proposed model with v d = 12.50 m/s over a 3000 m circular road.
Mathematics 11 04911 g004
Figure 5. Normalized flow for the proposed model with v d = 15.27 m/s over a 3000 m circular road.
Figure 5. Normalized flow for the proposed model with v d = 15.27 m/s over a 3000 m circular road.
Mathematics 11 04911 g005
Figure 6. Normalized flow for the ID model with δ = 1 ,   4 , and 20 over a 3000 m circular road.
Figure 6. Normalized flow for the ID model with δ = 1 ,   4 , and 20 over a 3000 m circular road.
Mathematics 11 04911 g006
Figure 7. Normalized speed for the proposed model with v d = 9.72 m/s over a 3000 m circular road: (a) P C I = 0 ; (b) P C I = 50 ; (c) P C I = 100 .
Figure 7. Normalized speed for the proposed model with v d = 9.72 m/s over a 3000 m circular road: (a) P C I = 0 ; (b) P C I = 50 ; (c) P C I = 100 .
Mathematics 11 04911 g007
Figure 8. Normalized speed for the proposed model with v d = 12.50 m/s over a 3000 m circular road: (a) P C I = 0 ; (b) P C I = 50 ; (c) P C I = 100 .
Figure 8. Normalized speed for the proposed model with v d = 12.50 m/s over a 3000 m circular road: (a) P C I = 0 ; (b) P C I = 50 ; (c) P C I = 100 .
Mathematics 11 04911 g008aMathematics 11 04911 g008b
Figure 9. Normalized speed for the proposed model with v d = 15.27 m/s over a 3000 m circular road: (a) P C I = 0 ; (b) P C I = 50 ; (c) P C I = 100 .
Figure 9. Normalized speed for the proposed model with v d = 15.27 m/s over a 3000 m circular road: (a) P C I = 0 ; (b) P C I = 50 ; (c) P C I = 100 .
Mathematics 11 04911 g009aMathematics 11 04911 g009b
Figure 10. Normalized speed for the ID model over a 3000 m circular road: (a) δ = 1 ; (b) δ = 4 ; (c) δ = 20 .
Figure 10. Normalized speed for the ID model over a 3000 m circular road: (a) δ = 1 ; (b) δ = 4 ; (c) δ = 20 .
Mathematics 11 04911 g010
Figure 11. Normalized density for the proposed model with v d = 9.72 m/s over a 3000 m circular road: (a) P C I = 0 ; (b) P C I = 50 ; (c) P C I = 100 .
Figure 11. Normalized density for the proposed model with v d = 9.72 m/s over a 3000 m circular road: (a) P C I = 0 ; (b) P C I = 50 ; (c) P C I = 100 .
Mathematics 11 04911 g011aMathematics 11 04911 g011b
Figure 12. Normalized density for the proposed model with v d = 12.50 m/s over a 3000 m circular road: (a) P C I = 0 ; (b) P C I = 50 ; (c) P C I = 100 .
Figure 12. Normalized density for the proposed model with v d = 12.50 m/s over a 3000 m circular road: (a) P C I = 0 ; (b) P C I = 50 ; (c) P C I = 100 .
Mathematics 11 04911 g012aMathematics 11 04911 g012b
Figure 13. Normalized density for the proposed model with v d = 15.27 m/s over a 3000 m circular road: (a) P C I = 0 ; (b) P C I = 50 ; (c) P C I = 100 .
Figure 13. Normalized density for the proposed model with v d = 15.27 m/s over a 3000 m circular road: (a) P C I = 0 ; (b) P C I = 50 ; (c) P C I = 100 .
Mathematics 11 04911 g013
Figure 14. Normalized density behavior for the ID model over a 3000 m circular road: (a) δ = 1 ; (b) δ = 4 ; (c) δ = 20 .
Figure 14. Normalized density behavior for the ID model over a 3000 m circular road: (a) δ = 1 ; (b) δ = 4 ; (c) δ = 20 .
Mathematics 11 04911 g014aMathematics 11 04911 g014b
Table 1. Simulation Parameters.
Table 1. Simulation Parameters.
ParameterValues
Desired speed for the proposed model, v d 9.72 m/s, 12.50 m/s and 15.27 m/s
Desired speed for the ID model, v d 20 m/s
Time headway for ID and proposed models, T 2.0 s
Critical density 0.25
Jam spacing, J 2.0 m
Maximum density 1 / J = 0.50
Maximum acceleration, a 0.73 m/s2
Vehicle length,   L 5.0 m
Acceleration exponent for the ID model, δ 1 ,   4 and 20
Pavement Condition Index, P C I 0 ,   50 and 100
Minimum acceleration, d 1.67 m/s2
Time step, t 0.50 s
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ali, F.; Khan, Z.H.; Khattak, K.S.; Gulliver, T.A.; Altamimi, A.B. A Microscopic Traffic Model Incorporating Vehicle Vibrations Due to Pavement Condition. Mathematics 2023, 11, 4911. https://doi.org/10.3390/math11244911

AMA Style

Ali F, Khan ZH, Khattak KS, Gulliver TA, Altamimi AB. A Microscopic Traffic Model Incorporating Vehicle Vibrations Due to Pavement Condition. Mathematics. 2023; 11(24):4911. https://doi.org/10.3390/math11244911

Chicago/Turabian Style

Ali, Faryal, Zawar Hussain Khan, Khurram Shehzad Khattak, Thomas Aaron Gulliver, and Ahmed B. Altamimi. 2023. "A Microscopic Traffic Model Incorporating Vehicle Vibrations Due to Pavement Condition" Mathematics 11, no. 24: 4911. https://doi.org/10.3390/math11244911

APA Style

Ali, F., Khan, Z. H., Khattak, K. S., Gulliver, T. A., & Altamimi, A. B. (2023). A Microscopic Traffic Model Incorporating Vehicle Vibrations Due to Pavement Condition. Mathematics, 11(24), 4911. https://doi.org/10.3390/math11244911

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop