C PDF
C PDF
C PDF
Appendix C
INTEGRATION Description
FORTRAN Code
ANOVA SAS File and Output
Simulation Output
112
Appendix C
Appendix C
This appendix presents a brief description of the INTEGRATION model that was used in
this thesis. The most important simulation features will be presented as well as the input
and output files of the model. Also included is the FORTRAN code used for the safety
model and the calculation of the accident risk.
C.1 Background
The INTEGRATION traffic simulation model was developed in the mid-1980s by Dr.
Michel Van Aerde. It is an integrated simulation and traffic assignment model (Van
Aerde, 1985). INTEGRATION is a dynamic and microscopic traffic simulation model.
The traffic flow in the model is represented as a series of individual vehicles that each
follow pre-specified macroscopic traffic flow relationships.
The input data are required to run the model and these are divided into fundamental and
advanced data. Table C-1 below provides a listing of the fundamental input data files and
their content in brief.
113
Appendix C
On-Screen Graphics/text
Simulation Run time Errors (runerr.out file)
Labeled Output Statistics (file10.dat file)
Summary Statistics (summary.out file)
Of greater interest in this thesis was the summary statistics file that provided a summary
of important aggregate network level measures of performance. The summary.out file is
automatically produced whenever a simulation run is completed. Data are categorized
according to the vehicle type (5 vehicle types) with a cumulative total across all vehicle
types. A sample summary output is attached for explanation purposes in the following
section.
The following is a summary output file after the sample network explained in
chapter 4 was simulated using INTEGRATION. Accident risk results are highlighted in
order to stress the output statistics produced by the model in terms of its safety model.
1 1 1 1 1 1 5 - vehicles
2 30 36 45 60 90 261 - veh-secs
3 0 0 0 0 0 4 - veh-km
4 0 0 0 0 0 0 - veh-stops
5 0 0 0 0 0 0 - wrng turn
6 0 0 0 0 0 0 - fuel (l)
7 0 0 0 0 0 0 - HC (g)
8 5 2 1 1 1 12 - CO (g)
9 0 0 0 0 0 1 - NO (g)
10 0 0 0 1 3 6 - acc* 10e6
11 0 0 0 0 0 0 - tot delay
12 0 0 0 0 0 0 - stp delay
13 0 0 0 0 0 0 - ac/dc del
14 0 0 0 0 0 0 - dollars
15 1 1 1 1 1 5 - persons
16 30 36 45 60 90 261 - per-hrs
17 0 0 0 0 0 4 - per-km
18 15 7 2 0 0 26 - acc-noise
19 0 0 0 0 0 0
20 0 0 0 0 0 0
114
Appendix C
C.5 Discussion
For the development of this thesis and for study purposes, INTEGRATION 2.20e
for Windows was used. The model included the safety model input, which was
furthermore analyzed for the purposes of this thesis. Input data (charts) as well as the
115
Appendix C
The FORTAN code below consists of a main program (main) and three subroutines
(acc_rat, acc_dam and acc_inj) as they were discussed in chapter 5.
program main
c
c - Macroscopic Safety Model 1.3 - Copyright M.Van Aerde(1999)
c - for MMDI Evaluation Study purposes only
c
c 1. rate1(i) - crashes/ million veh-sec type 1 to 14 plus total (15)
c 2. rate2(i) - crashes/ million veh-km type 1-14 plus total
c 3. dam(i) - property damage levels: 1 to 4
c a. No damage
c b. Minor damage
c c. Moderate damage
c d. Major damage
c 4. inj(i) - injury levels: 1 to 5
c a. No injury
c b. Minor injury
c c. Non-incapacitating injury
c d. Incapacitating injury
c e. Fatal injury
c
real rate1(15),rate2(15),dam(4),inj(5)
open(unit=9,file=acc_mod.out,status=unknown)
write(9,*) accident rates
c - speeds from 20 to 120 km/h
do 1100 i = 20,120,20
xspd = i
call acc_rat(xspd,rate1,rate2)
irate = 1
write(9,(f6.1,i5,15f9.5)) xspd,irate,(rate1(j),j=1,15)
irate = 2
write(9,(f6.1,i5,15f9.5)) xspd,irate,(rate2(j),j=1,15)
c - damage probabilities per crash
call acc_dam(xspd,dam)
itype = 1
write(9,(f6.1,i5,15f9.5)) xspd,itype,(dam(j),j=1,4)
c - injury probabilities per crash
call acc_inj(xspd,inj)
itype = 2
write(9,(f6.1,i5,15f9.5)) xspd,itype,(inj(j),j=1,5)
1100 continue
end
116
Appendix C
c ----------------------------------------------------------------
subroutine acc_rat(xspd,rate1,rate2)
real rate1(15),rate2(15),spd
real c1(15),c2(15)
data c1/-0.021501436,-0.023559071,-0.01709544,-0.019589495,
* -0.041110596,-0.023089944,-0.026041357,-0.032079684,
* -0.020100423,-0.020957402,-0.019926275,-0.019976246,
* -0.014536839,-0.017781349,-0.020170247/
data c2/-4.903643188,-5.013290246,-4.483434319,-4.283242043,
* -6.600147706,-5.340775704,-6.694593649,-5.641286798,
* -5.154935427,-4.715896567,-4.589321771,-5.073630589,
* -7.133251805,-4.848826862,-2.419018353/
c - spd is in mph
spd = xspd/1.6
c - model only valid from 20 to 120 km/h
if (spd.lt. 20) spd = 20
if (spd.gt. 70) spd = 70
do 1100 i =1,15
p1 = c1(i)*spd + c2(i)
c - crash rate per million veh-sec
rate1(i) = exp(p1)
c - accident rate per million veh-km
rate2(i) = rate1(i) / (spd*1.6) * 3600
c write(9,(a,i5,6f10.5)) p1 c1 c2 spd ,
c * i,p1,c1(i),c2(i),spd,rate1(i),rate2(i)
1100 continue
end
c ----------------------------------------------------------------
subroutine acc_dam(xspd,dam)
real dam(4)
real c(16)
spd = xspd/1.6
if (spd.lt. 20) spd = 20
if (spd.gt. 70) spd = 70
s1= spd
s2= spd**2.0
s3= spd**3.0
117
Appendix C
end
c ============================================================
subroutine acc_inj(xspd,inj)
real inj(5)
real c(20)
spd = xspd/1.6
if (spd.lt. 25) spd = 25
if (spd.gt. 70) spd = 70
s1= spd
s2= spd**2.0
s3= spd**3.0
end
c ============================================================
118
Appendix C
'$7$PPGL
,QILOH
H?KHVKDP?PPGL?SKRHQL[?JSV?ZLQB?DQRYD?UXUBVXPRXW
,QSXWVHTIODJFDUGD\SHUGLULVHT6SG6WRSV'HOD\)XHO+&&212[
&UDVK,QMXU\)DWDO0LQRU0RG0DMRU
SURFSULQW
YDUVHTIODJFDUGD\SHUGLULVHT6SG6WRSV'HOD\)XHO+&&212[
&UDVK,QMXU\)DWDO0LQRU0RG0DMRU
UXQ
SURFXQLYDULDWHQRUPDO
YDU6SG6WRSV'HOD\)XHO+&&212[&UDVK,QMXU\)DWDO0LQRU0RG
0DMRU
UXQ
SURFJOP
FODVVIODJGLUSHU
PRGHO6SG6WRSV'HOD\)XHO+&&212[&UDVK,QMXU\)DWDO0LQRU0RG
0DMRU IODJ_GLU_SHUVV
UXQ
SURFVRUW
E\GLUSHUIODJ
SURFPHDQV
E\GLUSHUIODJ
119
Appendix C
Extremes
Lowest Obs Highest Obs
23.99( 106) 55.28( 80)
24.72( 97) 55.77( 216)
25.15( 284) 55.85( 160)
25.18( 132) 56.23( 126)
25.22( 114) 56.47( 227)
Extremes
Lowest Obs Highest Obs
3.3( 126) 10.9( 89)
3.78( 282) 10.99( 132)
3.87( 216) 11.21( 79)
4.16( 84) 11.56( 96)
4.46( 160) 11.56( 97)
Extremes
Lowest Obs Highest Obs
18.08( 59) 43.41( 284)
19.74( 126) 43.83( 132)
19.87( 227) 44.41( 114)
20.13( 216) 44.48( 97)
20.19( 160) 46.56( 106)
120
Appendix C
Extremes
Lowest Obs Highest Obs
8.73( 59) 20.6( 284)
9.51( 126) 20.78( 132)
9.55( 227) 21.07( 97)
9.69( 80) 21.08( 114)
9.74( 216) 22.09( 106)
Extremes
Lowest Obs Highest Obs
0.08( 59) 0.14( 114)
0.09( 280) 0.14( 132)
0.09( 274) 0.14( 269)
0.09( 256) 0.14( 284)
0.09( 227) 0.15( 106)
Extremes
Lowest Obs Highest Obs
0.99( 80) 1.49( 96)
0.99( 59) 1.49( 122)
1.03( 244) 1.49( 132)
1.05( 286) 1.55( 114)
1.05( 256) 1.59( 106)
121
Appendix C
Extremes
Lowest Obs Highest Obs
7.58( 59) 21.68( 284)
8.25( 216) 21.89( 132)
8.33( 126) 21.98( 114)
8.44( 160) 22.5( 97)
8.44( 125) 23.23( 106)
Extremes
Lowest Obs Highest Obs
5.59( 59) 11.83( 284)
5.93( 80) 11.94( 132)
5.97( 227) 12( 97)
6.09( 126) 12.2( 114)
6.1( 244) 12.7( 106)
122
Appendix C
123
Appendix C
Dependent Variable: HC
Source DF Sum of Squares Mean Square F Value Pr > F
Model 11 1.05661239 0.09605567 8.28 0.0001
Error 289 3.35357300 0.01160406
Corrected Total 300 4.41018538
R-Square C.V. Root MSE HC Mean
0.239585 8.870992 0.10772214 1.21431894
Dependent Variable: CO
Source DF Sum of Squares Mean Square F Value Pr > F
Model 11 13.73945748 1.24904159 5.83 0.0001
Error 289 61.89824352 0.21418077
Corrected Total 300 75.63770100
R-Square C.V. Root MSE CO Mean
0.181648 2.981031 0.46279669 15.52471761
124
Appendix C
125
Appendix C
126
Appendix C
127
Appendix C
128
Appendix C
129
Appendix C
130
Appendix C
131
Appendix C
132
Appendix C
The INTEGRATION summary files are presented below. Two summary files refer to
before signal coordination and after signal coordination conditions for the road network
of Phoenix, Arizona (described in chapter 6).
C.8.1 Summary1.out
Before Scenario
1 123961 0 0 0 0 123961 - vehicles
2 67204328 0 0 0 0 67204328 - veh-secs
3 493410 0 0 0 0 493410 - veh-km
4 528153 0 0 0 0 528153 - veh-stops
5 0 0 0 0 0 0 - wrng turn
6 82441 0 0 0 0 82441 - fuel (l)
7 1250482 0 0 0 0 1250482 - HC (g)
8 8520539 0 0 0 0 8520539 - CO (g)
9 205691 0 0 0 0 205691 - NO (g)
10 1281017 0 0 0 0 1281017 - acc* 10e6
11 38096212 0 0 0 0 38096212 - tot delay
12 22248634 0 0 0 0 22248634 - stp delay
13 15847670 0 0 0 0 15847670 - ac/dc del
14 0 0 0 0 0 0 - dollars
15 123961 0 0 0 0 123961 - persons
16 67204328 0 0 0 0 67204328 - per-hrs
17 493410 0 0 0 0 493410 - per-km
18 -2147483648 0 0 0 0 -2147483648 - acc-noise
19 0 0 0 0 0 0
20 0 0 0 0 0 0
----------------------------------------------------------------------------------------------
133
Appendix C
C.8.2 Summary2.out
After Scenario
1 123961 0 0 0 0 123961 -vehicles
2 62416856 0 0 0 0 62416856veh-secs
3 493079 0 0 0 0 493079 -veh-km
4 527437 0 0 0 0 527437 -veh-stops
5 0 0 0 0 0 0 -wrng turn
6 80540 0 0 0 0 80540 -fuel (l)
7 1218181 0 0 0 0 1218181-HC (g)
8 8325895 0 0 0 0 8325895-CO (g)
9 204130 0 0 0 0 204130 -NO (g)
10 1219053 0 0 0 0 1219053-acc* 10e6
11 35034584 0 0 0 0 35034584-tot delay
12 19800048 0 0 0 0 19800048-stp delay
13 15234769 0 0 0 0 15234769-ac/dc del
14 0 0 0 0 0 0 -dollars
15 123961 0 0 0 0 123961 -persons
16 62416856 0 0 0 0 62416856-per-hrs
17 493079 0 0 0 0 493079 -per-km
18 -2147483648 0 0 0 0 -2147483648-acc-noise
19 0 0 0 0 0 0
20 0 0 0 0 0 0
-------------------------------------------------------------------------------------------------------------------------
1 1 0 0 0 0 1 -vehicles
2 503.52 0 0 0 0 503.52 -veh-secs
3 3.978 0 0 0 0 3.978 -veh-km
4 4.255 0 0 0 0 4.255 -veh-stops
5 0 0 0 0 0 0 -wrng turn
6 0.65 0 0 0 0 0.65 -fuel (l)
7 9.827 0 0 0 0 9.827 -HC (g)
8 67.165 0 0 0 0 67.165 -CO (g)
9 1.647 0 0 0 0 1.647 -NO (g)
10 9.834 0 0 0 0 9.834 -acc* 10e6
11 282.626 0 0 0 0 282.626 -tot delay
12 159.728 0 0 0 0 159.728 -stp delay
13 122.9 0 0 0 0 122.9 -ac/dc del
14 0 0 0 0 0 0 -dollars
15 1 0 0 0 0 1 -persons
16 503.52 0 0 0 0 503.52 -per-hrs
17 3.978 0 0 0 0 3.978 -per-km
18 115317.063 0 0 0 0 115317.063-acc-noise
19 0 0 0 0 0 0
20 0 0 0 0 0 0
134