Car Following Models

Lecture Notes in Transportation Systems Engineering

Prof. Tom V. Mathew
 
PDF for Online View or PDF for Print  
(PDF for IITB students and is password protected)

Contents

1 Overview-1
2 Car following models
 2.1 Notation
 2.2 Pipe’s model
 2.3 Forbes’ model
 2.4 General Motors’ model
 2.5 Optimal velocity model
3 General motor’s car following model
 3.1 Basic Philosophy
 3.2 Follow-the-leader model
  3.2.1 Numerical Example
4 Derivation of Greenberg’s formula (under development)
5 Summary
Exercises
References
Acknowledgments
______________________________________________________________________

1 Overview-1

Longitudinal spacing of vehicles are of particular importance from the points of view of safety, capacity and level of service. The longitudinal space occupied by a vehicle depend on the physical dimensions of the vehicles as well as the gaps between vehicles. For measuring this longitudinal space, two microscopic measures are used- distance headway and distance gap. Distance headway is defined as the distance from a selected point (usually front bumper) on the lead vehicle to the corresponding point on the following vehicles. Hence, it includes the length of the lead vehicle and the gap length between the lead and the following vehicles.

2 Car following models

Car following theories describe how one vehicle follows another vehicle in an uninterrupted flow. Various models were formulated to represent how a driver reacts to the changes in the relative positions of the vehicle ahead. Models like Pipes, Forbes, General Motors and Optimal velocity model are worth discussing.

2.1 Notation

Before going in to the details, various notations used in car-following models are discussed here with the help of figure 1. The leader vehicle is denoted as n and the following vehicle as (n + 1). Two characteristics at an instant t are of importance; location and speed. Location and speed of the lead vehicle at time instant t are represented by xnt and vnt respectively. Similarly, the location and speed of the follower are denoted by xn+1t and vn+1t respectively. The following vehicle is assumed to accelerate at time t + ΔT and not at t, where ΔT is the interval of time required for a driver to react to a changing situation. The gap between the leader and the follower vehicle is therefore xnt - xn+1t.


PIC

Figure 1: Notation for car following model


2.2 Pipe’s model

The basic assumption of this model is “A good rule for following another vehicle at a safe distance is to allow yourself at least the length of a car between your vehicle and the vehicle ahead for every ten miles per hour of speed at which you are traveling” According to Pipe’s car-following model, the minimum safe distance headway increases linearly with speed. A disadvantage of this model is that at low speeds, the minimum headways proposed by the theory are considerably less than the corresponding field measurements.

2.3 Forbes’ model

In this model, the reaction time needed for the following vehicle to perceive the need to decelerate and apply the brakes is considered. That is, the time gap between the rear of the leader and the front of the follower should always be equal to or greater than the reaction time. Therefore, the minimum time headway is equal to the reaction time (minimum time gap) and the time required for the lead vehicle to traverse a distance equivalent to its length. A disadvantage of this model is that, similar to Pipe’s model, there is a wide difference in the minimum distance headway at low and high speeds.

2.4 General Motors’ model

The General Motors’ model is the most popular of the car-following theories because of the following reasons:

  1. Agreement with field data; the simulation models developed based on General motors’ car following models shows good correlation to the field data.
  2. Mathematical relation to macroscopic model; Greenberg’s logarithmic model for speed-density relationship can be derived from General motors car following model.

In car following models, the motion of individual vehicle is governed by an equation, which is analogous to the Newton’s Laws of motion. In Newtonian mechanics, acceleration can be regarded as the response of the particle to stimulus it receives in the form of force which includes both the external force as well as those arising from the interaction with all other particles in the system. This model is the widely used and will be discussed in detail later.

2.5 Optimal velocity model

The concept of this model is that each driver tries to achieve an optimal velocity based on the distance to the preceding vehicle and the speed difference between the vehicles. This was an alternative possibility explored recently in car-following models. The formulation is based on the assumption that the desired speed vndesired depends on the distance headway of the nth vehicle. i.e.vndesiredt = voptxnt) where vopt is the optimal velocity function which is a function of the instantaneous distance headway Δxnt. Therefore ant is given by

 t        opt   t    t
an = [1∕τ][V  (Δx n)- vn]
(1)

where 1
τ is called as sensitivity coefficient. In short, the driving strategy of nth vehicle is that, it tries to maintain a safe speed which in turn depends on the relative position, rather than relative speed.

3 General motor’s car following model

3.1 Basic Philosophy

The basic philosophy of car following model is from Newtonian mechanics, where the acceleration may be regarded as the response of a matter to the stimulus it receives in the form of the force it receives from the interaction with other particles in the system. Hence, the basic philosophy of car-following theories can be summarized by the following equation

[Response]n α [Stimulus]n
(2)

for the nth vehicle (n=1, 2, ...). Each driver can respond to the surrounding traffic conditions only by accelerating or decelerating the vehicle. As mentioned earlier, different theories on car-following have arisen because of the difference in views regarding the nature of the stimulus. The stimulus may be composed of the speed of the vehicle, relative speeds, distance headway etc, and hence, it is not a single variable, but a function and can be represented as,

 t
an = fsti(vn,Δxn, Δvn)
(3)

where fsti is the stimulus function that depends on the speed of the current vehicle, relative position and speed with the front vehicle.

3.2 Follow-the-leader model

The car following model proposed by General motors is based on follow-the leader concept. This is based on two assumptions; (a) higher the speed of the vehicle, higher will be the spacing between the vehicles and (b) to avoid collision, driver must maintain a safe distance with the vehicle ahead.

Let Δxn+1t is the gap available for (n + 1)th vehicle, and let Δxsafe is the safe distance, vn+1t and vnt are the velocities, the gap required is given by,

Δxtn+1 = Δxsafe + τvtn+1
(4)

where τ is a sensitivity coefficient. The above equation can be written as

xn - xtn+1 = Δxsafe + τvtn+1
(5)

Differentiating the above equation with respect to time, we get

 vtn - vtn+1 = τ.atn+1
      1
atn+1 =--[vtn - vtn+1]
      τ
General Motors has proposed various forms of sensitivity coefficient term resulting in five generations of models. The most general model has the form,
       [           ]
at   =  αl,m-(vtn+1)m- [vt- vt  ]
 n+1    (xtn - xtn+1)l  n   n+1
(6)

where l is a distance headway exponent and can take values from +4 to -1, m is a speed exponent and can take values from -2 to +2, and α is a sensitivity coefficient. These parameters are to be calibrated using field data. This equation is the core of traffic simulation models.

In computer, implementation of the simulation models, three things need to be remembered:

  1. A driver will react to the change in speed of the front vehicle after a time gap called the reaction time during which the follower perceives the change in speed and react to it.
  2. The vehicle position, speed and acceleration will be updated at certain time intervals depending on the accuracy required. Lower the time interval, higher the accuracy.
  3. Vehicle position and speed is governed by Newton’s laws of motion, and the acceleration is governed by the car following model.

Therefore, the governing equations of a traffic flow can be developed as below. Let ΔT is the reaction time, and Δt is the updation time, the governing equations can be written as,

  vtn  =   vt-nδt+ atn-δt× δt                              (7)
  xt  =   xt-δt+ vt-δt× δt+ 1at-δtδt2                    (8)
   n      [n     n        ]2 n
 t           αl,m (vtn-+δ1t)m      t-ΔT    t- ΔT
an+1  =    --t--ΔT---t-ΔT-l- (vn    - vn+1 )             (9)
           (xn   - xn+1  )
The equation 7 is a simulation version of the Newton’s simple law of motion v = u + at and equation 8 is the simulation version of the Newton’s another equation s = ut + 1
2at2. The acceleration of the follower vehicle depends upon the relative velocity of the leader and the follower vehicle, sensitivity coefficient and the gap between the vehicles.

3.2.1 Numerical Example

Let a leader vehicle is moving with zero acceleration for two seconds from time zero. Then he accelerates by 1 m∕s2 for 2 seconds, then decelerates by 1m∕s2for 2 seconds. The initial speed is 16 m/s and initial location is 28 m from datum. A vehicle is following this vehicle with initial speed 16 m/s, and position zero. Simulate the behavior of the following vehicle using General Motors’ Car following model (acceleration, speed and position) for 7.5 seconds. Assume the parameters l=1, m=0 , sensitivity coefficient (αl,m) = 13, reaction time as 1 second and scan interval as 0.5 seconds.

Solution The first column shows the time in seconds. Column 2, 3, and 4 shows the acceleration, velocity and distance of the leader vehicle. Column 5,6, and 7 shows the acceleration, velocity and distance of the follower vehicle. Column 8 gives the difference in velocities between the leader and follower vehicle denoted as dv. Column 9 gives the difference in displacement between the leader and follower vehicle denoted as dx. Note that the values are assumed to be the state at the beginning of that time interval. At time t=0, leader vehicle has a velocity of 16 m/s and located at a distance of 28 m from a datum. The follower vehicle is also having the same velocity of 16 m/s and located at the datum. Since the velocity is same for both, dv = 0. At time t = 0, the leader vehicle is having acceleration zero, and hence has the same speed. The location of the leader vehicle can be found out from equation  as, x = 28+16×0.5 = 36 m. Similarly, the follower vehicle is not accelerating and is maintaining the same speed. The location of the follower vehicle is, x = 0+16×0.5 = 8 m. Therefore, dx = 36-8 =28m. These steps are repeated till t = 1.5 seconds. At time t = 2 seconds, leader vehicle accelerates at the rate of 1 m∕s2 and continues to accelerate for 2 seconds. After that it decelerates for a period of two seconds. At t= 2.5 seconds, velocity of leader vehicle changes to 16.5 m/s. Thus dv becomes 0.5 m/s at 2.5 seconds. dx also changes since the position of leader changes. Since the reaction time is 1 second, the follower will react to the leader’s change in acceleration at 2.0 seconds only after 3 seconds. Therefore, at t=3.5 seconds, the follower responds to the leaders change in acceleration given by equation i.e., a = 13×0.5
28.23 = 0.23 m∕s2. That is the current acceleration of the follower vehicle depends on dv and reaction time Δ of 1 second. The follower will change the speed at the next time interval. i.e., at time t = 4 seconds. The speed of the follower vehicle at t = 4 seconds is given by equation as v= 16+0.231×0.5 = 16.12 The location of the follower vehicle at t = 4 seconds is given by equation as x = 56+16×0.5+12×0.231×0.52 = 64.03 These steps are followed for all the cells of the table.



Table 1: Car-following example









t a(t) v(t) x(t) a(t) v(t) x(t) dv dx
(1) (2) (3) (4) (5) (6) (7) (8) (9)









t a(t) v(t) x(t) a(t) v(t) x(t) dv dx









0.00 0.00 16.00 28.00 0.00 16.00 0.00 0.00 28.00









0.50 0.00 16.00 36.00 0.00 16.00 8.00 0.00 28.00









1.00 0.00 16.00 44.00 0.00 16.00 16.00 0.00 28.00









1.50 0.00 16.00 52.00 0.00 16.00 24.00 0.00 28.00









2.00 1.00 16.00 60.00 0.00 16.00 32.00 0.00 28.00









2.50 1.00 16.50 68.13 0.00 16.00 40.00 0.50 28.13









3.00 1.00 17.00 76.50 0.00 16.00 48.00 1.00 28.50









3.50 1.00 17.50 85.13 0.23 16.00 56.00 1.50 29.13









4.00 -1.0018.00 94.00 0.46 16.12 64.03 1.88 29.97









4.50 -1.0017.50102.88 0.67 16.34 72.14 1.16 30.73









5.00 -1.0017.00111.50 0.82 16.68 80.40 0.32 31.10









5.50 -1.0016.50119.88 0.49 17.09 88.84 -0.5931.03









6.00 0.00 16.00128.00 0.13 17.33 97.45 -1.3330.55









6.50 0.00 16.00136.00-0.2517.40106.13-1.4029.87









7.00 0.00 16.00144.00-0.5717.28114.80-1.2829.20









7.50 0.00 16.00152.00-0.6116.99123.36-0.9928.64









8.00 0.00 16.00160.00-0.5716.69131.78-0.6928.22









8.50 0.00 16.00168.00-0.4516.40140.06-0.4027.94









9.00 0.00 16.00176.00-0.3216.18148.20-0.1827.80









9.50 0.00 16.00184.00-0.1916.02156.25-0.0227.75









10.00 0.00 16.00192.00-0.0815.93164.24 0.07 27.76









10.50 0.00 16.00200.00-0.0115.88172.19 0.12 27.81









11.00 0.00 16.00208.00 0.03 15.88180.13 0.12 27.87









11.50 0.00 16.00216.00 0.05 15.90188.08 0.10 27.92









12.00 0.00 16.00224.00 0.06 15.92196.03 0.08 27.97









12.50 0.00 16.00232.00 0.05 15.95204.00 0.05 28.00









13.00 0.00 16.00240.00 0.04 15.98211.98 0.02 28.02









13.50 0.00 16.00248.00 0.02 15.99219.98 0.01 28.02









14.00 0.00 16.00256.00 0.01 16.00227.98 0.00 28.02









14.50 0.00 16.00264.00 0.00 16.01235.98-0.0128.02









15.00 0.00 16.00272.00 0.00 16.01243.98-0.0128.02









15.50 0.00 16.00280.00 0.00 16.01251.99-0.0128.01









16.00 0.00 16.00288.00-0.0116.01260.00-0.0128.00









16.50 0.00 16.00296.00 0.00 16.01268.00-0.0128.00









17.00 0.00 16.00304.00 0.00 16.00276.00 0.00 28.00









17.50 0.00 16.00312.00 0.00 16.00284.00 0.00 28.00









18.00 0.00 16.00320.00 0.00 16.00292.00 0.00 28.00









18.50 0.00 16.00328.00 0.00 16.00300.00 0.00 28.00









19.00 0.00 16.00336.00 0.00 16.00308.00 0.00 28.00









19.50 0.00 16.00344.00 0.00 16.00316.00 0.00 28.00









20.00 0.00 16.00352.00 0.00 16.00324.00 0.00 28.00









20.50 0.00 16.00360.00 0.00 16.00332.00 0.00 28.00











PIC

Figure 2: Velocity vz Time



PIC

Figure 3: Acceleration vz Time


The earliest car-following models considered the difference in speeds between the leader and the follower as the stimulus. It was assumed that every driver tends to move with the same speed as that of the corresponding leading vehicle so that

 t   1  t+1   t
an = τ(vn  - vn+1)
(10)

where τ is a parameter that sets the time scale of the model and 1
τ can be considered as a measure of the sensitivity of the driver. According to such models, the driving strategy is to follow the leader and, therefore, such car-following models are collectively referred to as the follow the leader model. Efforts to develop this stimulus function led to five generations of car-following models, and the most general model is expressed mathematically as follows.

 t       αl,m [vtn-+δ1t]m    t-ΔT    t-ΔT
an+1 = [xt-ΔT---xt-ΔT]l(vn    - vn+1 )
         n      n+1
(11)

where l is a distance headway exponent and can take values from +4 to -1, m is a speed exponent and can take values from -2 to +2, and α is a sensitivity coefficient. These parameters are to be calibrated using field data.

4 Derivation of Greenberg’s formula (under development)

The Greenberg’s stream model can be derived from General Motors car follwing model. The GM model is given as:

 t   --αl,m-[vnt-+δ1t]m---  t-ΔT    t-ΔT
an = [xt-nΔT - xtn-+Δ1T]l(vn    - vn+1 )
(12)

Putting l = 1 and m = 0 the model becomes:

 t  αl,m(vt-n-ΔT--vtn-+Δ1T-)
an =  (xtn-ΔT - xt-ΔT )
                n+1
(13)

At steady state condtion, the equation becomes:

a =  α (vn-1 --vn)
 n   (xn-1 - xn)
(14)

Now integrating both sides with respect to tm we get:

∫       ∫
  a dt =   α-(vn-1---vn)dt
   n        (xn- 1 - xn)
(15)

Now, put u = xn-1 - xn, then

du= v    - v , or du = (v - v )dt.
dt   n- 1   n          n-1   n
(16)

Therefore,

 ∫         ∫  du
∴   andt = α  --.                          (17)
              u
∴ vn ==   ααllooggu(x+n-c1 - xn) + c                  ((1189))

     =   αlogsn +c.                           (20)
where, sn is the spacing. Now, under steady state conditions, there is no identity for vehicles. Hence, the above expression will be as follows:
v  =  α logs+ c                          (21)
           1
   =  α log--+ c.                         (22)
           k
To find the constant c, apply the boundary condition, namely, at jam density, the speed is zero. That is:
  0  =  α log 1kj +1 c
∴ c  =  - αlogk-
               j
∴ v  =  α log 1-- α log 1-
             k        kj
     =  α log 1-+ α log k
             k         j
  v  =  α log kj.                             (23)
             k
Now, to find α we do the following:
               k
q = kv = α k log-j                          (24)
               k
Differentiating with respect to density and equating to zero:
     ddqk- =  α k [[ddkd-log kkj]+dα (log]kkj) (   k )
         =  α k  ---logkj - --logk  +α  log-j
                [dk ]    ( dk  )           k
         =  α k  --1 + α  log kj
                  k          k
                 kj
         =  α log k - α = 0.
     kj
∴ log k   =  1,
     kj               kj
  or k   =  e, or, k0 = e .                             (25)
Noted that we denote k corresponding to the maximum flow as k0. Now, substituting the above result, we get the value of alpha as
v0 =   αloge = α.
                                         (26)
Now the Greenberg’s stream model is:
v = v0log kjk
(27)

5 Summary

Microscopic traffic flow modeling focuses on the minute aspects of traffic stream like vehicle to vehicle interaction and individual vehicle behavior. They help to analyze very small changes in the traffic stream over time and space. Car following model is one such model where in the stimulus-response concept is employed. Optimal models and simulation models were briefly discussed.

Exercises

  1. Discuss the concepts and model formulations of Generalised GM model, Gipps’ model, and Wiedemann 74 car-following models.
  2. Discuss the evolution of the five generations of General Motors car following model highlighting the drawback and advantages of each generation of the model.
  3. A line of vehicles are in car following mode and all vehicles are travelling at 18 m/s with distance headway of 20 m. After 1.2 seconds, the lead vehicle suddenly decelerates at a rate of 1.2 m∕s2 until it stops completely. simulate the behaviour of first following vehicle using the GM fifth car following model for the first 2.5 seconds. Tabulate the results. Assume headway exponent 1.2, speed exponent 1.6, sensitivity coefficient 0.8, reaction time 0.6 seconds, and scan interval 0.3 seconds.
  4. Simulate the following vehicle behaviour for the following data using Widemann 74 model. (a) For the case of stand still distance 3.5m, additive part of safety distance 1.5, and multiplicative part of safety distance 0.8. (b) For the case of stand still distance 3.5m, additive part of safety distance 1.5, and multiplicative part of safety distance 0.8. Comment on the following vehicle behaviour for the above two cases.
  5. A car is travelling with a speed of 16 m/sec at time t=0. Another car follows the first at a distance of 28 m with same velocity. If the first car accelerated by 1 m/sec2 from t=1 to 2 and decelerate by 1 m/sec2 from t=2 to 3, find the speed, acceleration and spacing of the follower at time t=3.0 sec. Assume the reaction time is 1 sec, vehicle dynamics are updated every 0.5 seconds, and the car following model is given by Eq. 28. (Use of a tabular form is encouraged).
                 [                    ]
an+1(t) = 15×  un(t-- 1)---un+1(t--1)
               xn(t - 1) - xn+1(t- 1)
    (28)

  6. In a simulation experiment on a single lane road, one vehicle is travelling at 18 m∕s. After 1.5 seconds, the vehicle suddenly accelerates at a rate of 1.5 m∕s2 for the next 1.8 seconds. Simulate the behaviour of subsequent vehicle with an initial speed of 16 m/s using GM fifth car following model for the first 3 seconds if the initial distance headway is 20 m. Tabulate the results. Assume headway exponent 1.2, speed exponent 1.5, sensitivity coefficient 0.8, reaction time 0.6 seconds, and update interval of 0.3 seconds.
  7. A line of vehicles are in car following mode and all vehicles are travelling at 15 m/s with distance headway of 20 m. After 1.2 seconds, the lead vehicle suddenly decelerates at a rate of 1.2 m∕s2 until it stops completely. simulate the behaviour of first following vehicle using the GM fifth car following model for the first 2.5 seconds. Tabulate the results. Assume headway exponent 1.2, speed exponent 1.6, sensitivity coefficient 0.6, reaction time 0.6 seconds, and scan interval 0.3 seconds.
  8. In a simulation experiment on a single lane road, one vehicle is travelling at 16 m∕s. After 0.6 seconds, the vehicle suddenly accelerates at a rate of 1.2 m∕s2 for the next 0.9 seconds. Simulate the behaviour of subsequent vehicle with an initial speed of 16 m/s using GM fifth car following model for the first 2.1 seconds if the initial distance headway is 25 m. Tabulate the results. Assume headway exponent 1.2, speed exponent 1.4, sensitivity coefficient 0.6, reaction time 0.6 seconds, and update interval of 0.3 seconds.
  9. A line of vehicles are in car following mode and all vehicles are travelling at 15 m/s with distance headway of 25 m. After 1 second, the lead vehicle suddenly decelerates at a rate of 1.2m∕s2 until it stops completely. Simulate the behaviour of first following vehicle using the GM fifth car following model for the first 3 seconds. Tabulate the results. Assume headway exponent 1.0, speed exponent 1.5, sensitivity coefficient 0.5, reaction time 0.5 seconds, and scan interval 0.25 seconds.

Program Codes

GMV CF Model

Instructions

Modify the input file and use. To compile in linux use gcc prog.c -lm -o prog.exe and to run .\prog.exe.

Input File
GMV_Sensitivity_Coefficient_a 26.0  
GMV_Distance_Exponent_______l  2.0  
GMV_Speed_Exponent__________m  1.0  
Reaction_Time______________DT  1.0  
Update_Interval____________dt  0.5  
Initial_Speed_Leader_______vl 16.0  
Initial_Position_Leader____xl 28.0  
Initial_Speed_Follower_____vf 16.0  
Initial_Position_Follower__xf  0.0  
No_of_Update_Intervals______n   50  
Leader_Accelration_Values_for_each_n_update_intervals_below  
 0.0    0.0    0.0    0.0    1.0    1.0    1.0    1.0   -1.0  -1.0  
-1.0   -1.0    0.0    0.0    0.0    0.0    0.0    0.0    0.0   0.0  
 0.0    0.0    0.0    0.0    0.0    0.0    0.0    0.0    0.0   0.0  
 0.0    0.0    0.0    0.0    0.0    0.0    0.0    0.0    0.0   0.0  
 0.0    0.0    0.0    0.0    0.0    0.0    0.0    0.0    0.0   0.0

C-code
#include<math.h>  
#include<stdio.h>  
#include<stdlib.h>  
#define MAX_DUR 1000  
int main()  
{  
   char text[10][100];  
   int curTime=0;         /* current time in seconds */  
   int maxTime=0;         /* maximum simulation time in seconds */  
   int DelT=1;            /* reaction time in no of update intervals = DT/dt */  
   float delX,delV;       /* relative speed and relative distance */  
   float DT=0.0;          /* reaction time in seconds */  
   float dt=0.0;          /* updated interval in seconds */  
   float gm5_a=0.0;  
   float gm5_l=0.0;  
   float gm5_m=0.0;  
   float vFol[MAX_DUR];  
   float xFol[MAX_DUR];  
   float aFol[MAX_DUR];  
   float vLdr[MAX_DUR];  
   float xLdr[MAX_DUR];  
   float aLdr[MAX_DUR];  
   float clockTime[MAX_DUR]={0.0}; /* simulaion clock time in seconds */  
   FILE* finp=fopen("input.dat","r");  
 
   /* Read all input data from the file */  
   fscanf(finp,"%s %f",text[0],&gm5_a);  
   fscanf(finp,"%s %f",text[1],&gm5_l);  
   fscanf(finp,"%s %f",text[2],&gm5_m);  
   fscanf(finp,"%s %f",text[3],&DT);  
   fscanf(finp,"%s %f",text[4],&dt);  
   fscanf(finp,"%s %f",text[5],&vLdr[0]);  
   fscanf(finp,"%s %f",text[6],&xLdr[0]);  
   fscanf(finp,"%s %f",text[7],&vFol[0]);  
   fscanf(finp,"%s %f",text[8],&xFol[0]);  
   fscanf(finp,"%s %d",text[9],&maxTime);  
   fscanf(finp,"%s",text[11]);  
   for(curTime=0;curTime<maxTime;curTime++)  
      fscanf(finp,"%f",&aLdr[curTime]);  
 
   /* Compuation of the states of leader and follower */  
   for(curTime=1;curTime<maxTime;curTime++)  
   {  
      clockTime[curTime]=clockTime[curTime-1]+dt;  
      vLdr[curTime]=vLdr[curTime-1]+aLdr[curTime-1]*dt;  
      xLdr[curTime]=xLdr[curTime-1]+vLdr[curTime-1]*dt+0.5*aLdr[curTime-1]*dt*dt;  
      if(clockTime[curTime]>DT) DelT=(int)(DT/dt);  
      delX=xLdr[curTime-DelT]-xFol[curTime-DelT];  
      delV=vLdr[curTime-DelT]-vFol[curTime-DelT];  
      aFol[curTime]=gm5_a*pow(vFol[curTime-1],gm5_m)/pow(delX,gm5_l)*delV;  
      vFol[curTime]=vFol[curTime-1]+aFol[curTime-1]*dt;  
      xFol[curTime]=xFol[curTime-1]+vFol[curTime]*dt+0.5*aFol[curTime-1]*dt*dt;  
   }  
 
   /* Printing all results */  
   printf("%s\t%.1f",text[0],gm5_a);  
   printf("\n%s\t%.1f",text[1],gm5_l);  
   printf("\n%s\t%.1f",text[2],gm5_m);  
   printf("\n%s\t%.1f",text[3],DT);  
   printf("\n%s\t%.2f",text[4],dt);  
   printf("\n%s\t%.2f",text[5],vFol[0]);  
   printf("\n%s\t%.2f",text[6],xFol[0]);  
   printf("\n%s\t%.2f",text[7],vLdr[0]);  
   printf("\n%s\t%.2f",text[8],xLdr[0]);  
   printf("\n%s\t%d",text[9],maxTime);  
   printf("\nNo\tTime\taLdr\tvLdr\txLdr\taFol\tvFol\txFol");  
   for(curTime=0;curTime<maxTime;curTime++)  
   {  
      printf("\n%d\t%.2f\t%f\t%f\t%f\t%f\t%f\t%f",\  
         curTime+1,clockTime[curTime],aLdr[curTime],vLdr[curTime],\  
         xLdr[curTime],aFol[curTime],vFol[curTime],xFol[curTime]);  
   }  
   return 0;  
}  
/* End of File */

References

  1. L R Kadiyali. Traffic Engineering and Transportation Planning. Khanna Publishers, New Delhi, 1987.
  2. L J Pignataro. Traffic Engineering: Theory and practice. Prentice-Hall, Englewoods Cliffs,N.J., 1973.

Acknowledgments

I wish to thank several of my students and staff of NPTEL for their contribution in this lecture. I also appreciate your constructive feedback which may be sent to tvm@civil.iitb.ac.in

Prof. Tom V. Mathew
Department of Civil Engineering
Indian Institute of Technology Bombay, India

_________________________________________________________________________

Thursday 31 August 2023 12:13:06 AM IST