Experimental Study of Behaviour of Circu
Experimental Study of Behaviour of Circu
Experimental Study of Behaviour of Circu
An International Journal
To cite this article: Danny Useche-Infante, Gonzalo Aiassa Martinez, Pedro Arrúa & Marcelo
Eberhardt (2019): Experimental study of behaviour of circular footing on geogrid-reinforced sand,
Geomechanics and Geoengineering, DOI: 10.1080/17486025.2019.1683621
Article views: 25
The aim of present experimental investigation is to method was used to determine the particle size distribu-
evaluate the effect of different parameters of the sand- tion of the sand and the results are shown in Figure 1.
geogrid configuration, on the behaviour of a circular Standard triaxial compression tests on this sand at relative
footing supported on reinforced sand. The parameters densities of 30%, 60% and 90% revealed internal friction
considered in the model tests are: top layer spacing (u), angles of 34.3º, 40.5º and 45.4º, respectively.
diameter of reinforcement (D), vertical spacing between
reinforcement (h), number of layers (N), embedment
depth of footing (Df), type of the reinforcement and the Geosynthetic specimens
effects of wraparound ends of the geogrid (see Figure 3). Two types of geosynthetics were used in the experimen-
Finally, with the results obtained from test, two regres- tal programme: uniaxial and biaxial geogrid. Figure 2
sion models are proposed for prediction of the bearing shows the geogrid used as reinforcement in this study.
capacity of a circular foundation supported on sand The uniaxial geogrid is made of polyvinyl alcohol
reinforced with geogrid. The expressions obtained (PVA). The biaxial geogrid is fabricated of polypropy-
were compared with different analytical and multiple lene (PP) yarns. In Table 1, the mechanical strength
regression models reported in the technical literature. parameters of the two types of geosynthetics used in
the tests are shown.
Figure 2. Geosynthetics used in the tests (Dimensions in mm), (a) uniaxial geogrid, (b) biaxial geogrid.
Table 1. Mechanical properties of geogrid (provided from the suppliers' technical brochures).
Property Reinforcement
Type of geosynthetic material Uniaxial geogrid Biaxial geogrid
Type of polymer Polyvinyl (PVA) Polypropylene (PP)
Nominal mass per unit area (g/m2) 240 200
Modulus to def. 5% LDa 630 ≥ 360
(kN/m) CDb - ≥ 360
Modulus to def. 2% (kN/m) LDa 700 ≥ 400
CDb - ≥ 400
Mesh opening (mm) LDa 20 40
CDb 30 40
Thickness (mm) LDa 1.0 1.5
CDb 1.0 1.0
a
Longitudinal direction.
b
Cross machine direction.
rough circular-steel footing with diameter, B, of 100 mm to Chummar (1972), the soil failure surface extends
and thickness of 15 mm, is placed over a geogrid- around of 2B on each side from the edge of the footing
reinforced sand. A hydraulic jack with 50-ton capacity and its depth is about 1.1B from the base of the founda-
welded against a reaction frame was used to apply the tion. Therefore, in the tests carried out, it is guaranteed
load to the footing. The diameter of the foundation and that the fault surface develops freely. Besides, the B/D50
size of the tank were chosen in such a way that the edge ratio was 64.1, accordingly, the size of the soil particles
effect was minimum. The diameter of the tank was has no significant influence on the results (see Kusakabe
6B and the height of the soil stratum 4.5B. According (1995)).
4 D. USECHE-INFANTE ET AL.
Figure 3. Schematic diagram of the test set-up, (a) Series I to VI, (b) Series VII.
The tank was fixed and supported by two steel beams, density of 30%, the sand was placed by dry pluviation
and connected to four steel columns which were placed in method, at a constant height of fall to achieve uniformity
the floor of the lab firmly by anchor bolts. The rate of the of sand index properties. The density was obtained by
vertical movement of the foundation was kept to calibrating the height of fall and the size of the hole of the
a constant value of 1 mm/min for the entire test. The spreader. The calculated fall height was 0.8 m, which
amount of load applied was measured by using a load remained constant during the deposition. At the required
cell. The displacements were measured by a digital dial. depth according to the geometric configuration, the sand
Data were recorded by means of a data logger that was raining was stopped to place the geosynthetic sheet. The
connected to the computer for processing. A schematic density of the sample was controlled with small alumi-
diagram of the test set-up is shown in Figure 3 and the test nium cans of known volumes placed in different parts of
apparatus is shown in Figure 4. the tank. Differences in densities below 2% were mea-
sured at different locations in the test tank. For relative
densities of 60% and 90%, the soil layers were prepared
Sample preparation with dynamic compaction using a 150 mm in diameter
vibratory poker to achieve the desired relative density.
The test tank was filled with sand in five layers, four layers Proper care was given to achieve the same dry unit weight
of 100 mm thick and a top layer of 50 mm thick; main- in all the five lifts, in order to have a homogenous sand
taining a total thickness of the sand bed of 450 mm. In bed. For the preparation of reinforced sand sample, geo-
order to study the influence of relative density of founda- grid is placed at desired depth from bottom of footing
tion soil on the performance of geogrid-reinforced foun- after levelling the surface to make it horizontal. Model
dation beds (series I), the model tests were performed at footing was then placed so that the centre of the plate
relative densities of 30%, 60%, and 90%. In the other coincides with the centre of the loading jack to avoid
series of tests, the sand beds were maintained at fixed eccentric loading. (In this paragraph the following is
relative density of 60%. For samples with a relative added: At the required depth according to the geometric
GEOMECHANICS AND GEOENGINEERING 5
configuration, the sand raining was stopped to place the to find the sand-geogrid configuration that provides the
geosynthetic sheet. The density of the sample was con- greatest increase in the bearing capacity. In each series
trolled with small aluminium cans of known volumes a parameter was studied while the other variables were
placed in different parts of the tank. Differences in den- kept constant. Table 2 summarises the test programme
sities below 2% were measured at different locations in and test variables.
the test tank.) Under series I, samples were prepared with different
relative densities. In this series of tests, the effect of sand
density on the improvement obtained in reinforced soils
Experimental testing programme was determined.
Seven series of tests were performed to analyse the effect Test series II were conducted by varying the depth of
produced in the footing when geogrid layers are the top geogrid layer (u). The objective of these tests is
included in the sand. Tests were carried out in order to determine the depth of the geogrid layer with which
6 D. USECHE-INFANTE ET AL.
Figure 5. Variation of bearing pressure with footing settlement for different relative densities test, series I (N = 1, u/B = 0.25, D/B = 3, Df/B = 0),
(a) Uniaxial geogrid and (b) Biaxial geogrid.
8 D. USECHE-INFANTE ET AL.
Figure 6. Variation of bearing pressure with settlement for different depths of the top geogrid layer, series II (Dr = 60%, N = 1, D/B = 3,
Df/B = 0), (a) Uniaxial geogrid and (b) Biaxial geogrid.
Figure 7. Variation of bearing pressure with settlement for different diameters of geogrid layer, series III (Dr= 60%, N= 1, u/B = 0.25, Df/B = 0),
(a) Uniaxial geogrid and (b) Biaxial geogrid.
Figure 8. Bearing pressure vs. Footing settlement for different spacings of geogrid layers, series IV (Dr = 60%, N = 2, u/B = 0.25, D/B =
4, Df/B = 0), (a) Uniaxial geogrid and (b) Biaxial geogrid.
dense sands and the second type of behaviour occurs in Effect of relative density variation
medium-dense sands (Vesic 1973). When three or more
Bearing pressure versus settlement responses for different
layers of geosynthetic are used, it can be seen that the
values of relative density are shown in Figure 5. The pres-
typical behaviour for dense sands is altered by the develop-
sure – settlement responses, presented a visible improve-
ment of the tensile force in the geogrid layers.
ment in the bearing capacity of the soil as the relative
GEOMECHANICS AND GEOENGINEERING 9
Figure 9. Bearing pressure vs. Footing settlement for different numbers of reinforcement layers, series V (Dr= 60%, u/B = 0.25, D/B = 4,
h/B = 0.25, Df/B = 0), (a) Uniaxial and (b) Biaxial geogrid.
Figure 10. Bearing pressure vs. Footing settlement for different embedments of a footing, series VI (Dr = 60%, N = 1, u/B = 0.25, D/B = 4), (a)
Uniaxial geogrid and (b) Biaxial geogrid.
Figure 11. Bearing pressure vs. Footing settlement for geogrid reinforcement with wraparound ends, series VII (Dr = 60%, N = 1, u/B = 0.25,
D/B = 4), (a) Uniaxial and (b) Biaxial geogrid.
density increases. Variation of the bearing capacity ratio It could be observed that as the relative density increases
(BCR) with relative density of soil Dr, is shown in Figure 12. the bearing capacity ratio also increases, which is in good
10 D. USECHE-INFANTE ET AL.
Figure 12. Variations of BCR with different relative densities, (a) for ultimate bearing capacity and (b) for settlement ratio, s/B, of 0.03
and 0.10.
agreement with the finding by Durga-Prasad et al. (2016). that develops in it has less effect on the stress influence
The BCR presented values between 1.13 and 2.13 for the zone. In the results of the series II, a greater increase in
series I, as can be seen in Table 2, this shows a significant the bearing capacity was obtained when u = 0.25B,
enhancement in the bearing capacity of the soil for the however, to find the optimum value of the depth of
different densities considered. The inclusion of the geosyn- the first layer, it is necessary to carry out tests with u <
thetic generates a stiffening effect in the soil, which pro- 0.25B.
duces a greater improvement in the reinforced samples. (In
this paragraph the citation of the document by Fragaszy
and Lawton (1985) was changed to that of Durga-Prasad Effect of the diameter of reinforcement
et al. (2016). There was an error in the reference. The work
In order to determine the size of the reinforcement that
reported by Durga-Prasad et al. (2016) presents the
produces the optimal anchorage between the geogrid and
improvement in load-carrying capacity of footing on rein-
the soil, in the series III was studied the effect of the
forced sand was higher for dense sand beds and the experi-
diameter of the geogrid on the footing behaviour. The
mental conditions are similar to those presented in this
stress–settlement curves for different values of geogrid
paper.)
diameter ratio (D/B) are shown in Figure 7. This figure
shows that the bearing capacity of the soil begins to grow
for values of D/B ≥ 3, with the bearing capacity values
Effect of reinforcement’s top spacing
being higher for the samples with D/B = 4. For D/B = 5
Figure 6 shows that the bearing capacity of the soil the increase in bearing capacity does not change signifi-
decreases as the depth of the geogrid layer increases, to cantly in relation to that achieved with D/B = 3 and 4.
a value from which there is no increase in the soil This behaviour illustrates that a sufficient diameter must
carrying capacity (for u/B > 0.75). This can occur be provided to the geogrid to fully develop the pull-out
because when the depth of the geogrid layer is increased capacity in the reinforcement. Figure 14 shows the varia-
too much, the failure occurs between the base of the tions of BCR with the geogrid diameter ratio. The BCR
foundation and the reinforcement layer (Chen 2007); increases with increasing geogrid diameter, the highest
therefore, the tensile force in the geosynthetic is not values of increase in the bearing capacity are reached for
developed. The variation of BCR with u/B is shown in the samples with D/B ≥ 3. For values of D/B > 5, the size
Figure 13. The results indicate that the BCR values at the of the geogrid layer is not cost-effective. Therefore, the
ultimate loads and for s/B = 0.03 and 0.10 generally best performance of reinforced soil occurs when D/B is in
decreased as top layer spacing increased. Abu-Farsakh the range between 3 and 5. This is consistent with the
et al. (2013) reported the same observation. This beha- recommended size of geogrid reported by El-Sawwaf
viour occurs because there is a direct relationship y Nazir (2010), Chakraborty and Kumar (2014) and
between the location of the geogrid compared to the Tavangar and Shooshpasha (2016), who found optimal
size of the stress influence zone. The average vertical values between 3 and 5B. Similar to results obtained in the
stress increase below the centre of a circularly loaded series II, the relation between reinforcement diameter
area has greater influence at depths close to the base of and size of the stress influence zone becomes visible.
the foundation (see Saikia (2012)). Therefore, when the The length of this zone of influence for depths close to
depth of geosynthetic location is increased, the tension the base of the foundation is around 4B (see Chummar
GEOMECHANICS AND GEOENGINEERING 11
Figure 13. Variations of BCR with u/B ratio, (a) for ultimate bearing capacity and (b) for settlement ratio, s/B, of 0.03 and 0.10.
Figure 14. Variations of BCR with D/B ratio, (a) for ultimate bearing capacity and (b) for settlement ratio, s/B, of 0.03 and 0.10.
(1972)). For this reason, the samples in which the dia- of the geogrids in relation to the size of the stress
meter of the geogrid crossed the stress influence zone had influence zone.
a higher increase in capacity.
Figure 15. Variations of BCR with h/B ratio, (a) for ultimate bearing capacity and (b) for settlement ratio, s/B, of 0.03 and 0.10.
Figure 16. Variations of BCR with N, (a) for ultimate bearing capacity and (b) for settlement ratio, s/B, of 0.03 and 0.10.
number of geogrid layers, for N > 3 the growth rate in Effect of the wraparound of the geogrid
the BCR was much lower. Therefore, when three layers reinforcement ends
of geogrid were used, the optimum increase in bearing
The behaviour of circular foundations resting on sand
capacity was obtained.
beds reinforced with a single layer of the geogrid with
wraparound ends was studied. This technique can be
used when there is not enough space to place the geo-
Effect of footing depth
synthetic layer with the optimum size. The variable
The effect of footing depth ratio (Df/B) of single layer of taken into account in this series is the folded length, l/
geogrid reinforcement was studied by applying an over- B. Samples were prepared with l/B = 0, 0.3, 0.6 and 1.0
load of Df/B = 0.25, 0.5, 0.75 and 1.0. Figure 10 indicates (see Figure 3 and Table 2). In Figure 11, it is observed
that the bearing capacity of the unreinforced and rein- that the wraparound ends cause an increase in the
forced soil is increased when the depth of foundation is bearing capacity of the reinforced soil. Figure 18 shows
increased. Similar results were obtained by Sitharam and the variation of the BCR against l/B. The results of the
Sireesh (2004), Gu (2011) and Kazi et al. (2016). On the tests indicated that the geogrid with wraparound ends
other hand, Figure 17 shows the variations of BCR with the (with l/B = 0.3), produces an increase in the BCR that
footing depth ratio. The figure shows that as the footing depends on the type of reinforcement, in this case the
depth increases, the BCR decreases gradually. The decrease improvement was 10% and 5% for the samples with
in the BCR is due to the fact that a greater increase in the uniaxial and biaxial geogrid, respectively. For l/B > 0.6,
bearing capacity of the soil was obtained in the unrein- the enhancement is insignificant. The improvement
forced samples compared to the reinforced samples when obtained is attributed to the confinement effect pro-
the footing depth was increased. The trend presented by duced when the edges of the geogrid are folded
BCR vs. Df/B is similar to that found by Gu (2011). (Shukla 2012, Kazi et al. 2015).
GEOMECHANICS AND GEOENGINEERING 13
Figure 17. Variations of BCR with Df/B, (a) for ultimate bearing capacity and (b) for settlement ratio, s/B, of 0.03 and 0.10.
Figure 18. Variations of BCR with l/B, (a) for ultimate bearing capacity and (b) for settlement ratio, s/B, of 0.03 and 0.10.
Effect of the type of geosynthetic material relation to the main footing can have an effect on the
increase of the bearing capacity. Investigations carried
Two different types of geogrid with different tensile mod-
out with strips and rectangular footings show important
uli were used in the model footing tests, uniaxial and
increases in the bearing capacity contributed by uniaxial
biaxial geogrid (see Table 2). According to the results
geogrid (Patra et al. 2005, El-Sawwaf and Nazir 2010).
obtained in the tests, the improvement in the bearing
capacity of the soil depended on the type of geogrid
used. As seen in Figure 5 to 18; the sand reinforced by
Regression models
the biaxial geogrid, which has a higher modulus and
resists greater tensile stresses in the transverse ribs than Basic model
the uniaxial geogrid, reached higher values of bearing
The following basic regression model was selected to
capacity than sand reinforced by the uniaxial geogrid.
perform the data analysis:
This difference can be attributed to the resistance pro-
vided by the transverse fibres and the interaction between y ¼ 1 x1 þ 2 x2 þ ::: þ p xp (2)
the soil particles and the geosynthetic fibres, which is
presented differently in the two types of geogrid. The where y represents the estimated value for y’, x1, x2, . . ., xp
interaction between soil and the two types of geogrid are independent variables and ξ1, ξ2, . . ., ξp are unknown
used in this research has been studied by the authors in parameters and p is the number of independent variables.
previous investigations (Useche-Infante et al. 2016). The ξi coefficients are selected in such a way that the sum
Besides, a study conducted by Liu et al. (2009) showed of squared residuals Σ(y – y’)2 is again minimised.
the contribution of geogrid transverse bars to the shear To evaluate the degree of adjustment of the regres-
strength at soil–geogrid interface. The direction of geo- sion line to the data, the coefficient of determination R2
grid placement (especially of the uniaxial geogrid) in was used, which is defined as:
14 D. USECHE-INFANTE ET AL.
The value of R2 will never decrease by adding some where wi0 is the positive square root of C-1iiEs2, where
explanatory variables; an alternative to evaluate the C-1ii is the ith diagonal element of the matrix (x´x)−1.
model in these cases is to make an adjustment to the The decision rule for the ‘t’ test is:
coefficient R2. The adjusted coefficient of determination
(R2adj) is calculated as follows: ● Reject hypothesis if tcal > t(1 – α/2, n – p) o, tcal < – t
(1 – α/2, n – p).
2
2
n 1 ● Accept hypothesis if – t(1 – α/2, n – p) ≤ tcal ≤ t(1 –
R adj ¼ 1 1 R (4)
n p α/2, n – p).
where n and p are the number of independent variables where t(1 – α/2, n – p) is a value chosen from the
and the number of observations, respectively. t table for the appropriate level of significance, α.
The efficiency of the model was measured with the
standard error. This value is an unbiased estimate of the
variance. The prediction obtained with the model will Linear and non-linear model
be better when the value of Es is smaller. Es2 is given by:
The results of the load tests showed that the soil type,
ðy y0 Þ2
P
2 the geosynthetic type and the geometric parameters of
Es ¼ (5)
ðn pÞ soil-geogrid system influence the stress-settlement
behaviour of the foundation. The following parameters
were included in the regression models: relative density
Significance of regression coefficients (Dr), settlement of the footing ratio (s/B), depth of the
first layer ratio (u/B), vertical separation between layers
In order to determine that if the assumed equation is
ratio (h/B), diameter of the reinforcement ratio (D/B),
statistically significant or not, F and t statistics were
number of layers (N) and tensile modulus of reinforce-
calculated. The significance of the regression coefficient
ment (J). The linear regression model used to investigate
as a whole can be assessed by the F test, while the
the effects of all variables is expressed as follows:
significance of the partial regression coefficient can be
determined by t test.
Dr s u D
With the F test, the variance contributed by the BCR ¼ 1 þ 2 þ 3 þ 4 þ 5
100 B B B
regression and the error variance are compared. This h J
test allows evaluating the suitability of the model. The þ 6 þ 7 ðN Þ þ 8
B 100kN=m
null hypothesis is that the partial regression coefficients, (8)
ξ1, ξ2, . . ., ξp are equal to zero. From the definition of R2;
F can be expressed as: After performing the procedure described in the previous
sections and obtaining the values of the coefficients, the
½R2 =ðp 1Þ following linear regression equation was obtained:
Fcal ¼ (6)
½ð1 R2 Þ=ðn pÞ
Dr s u
BCR ¼ 1:76 þ 0:862 0:454 0:833
The decision rule for the ‘F’ test is: 100 B B
D h
þ :::::: þ 0:191 0:539 þ
● Reject hypothesis if Fcal > F(1 – α, p – 1, n – p). B B
● Accept hypothesis if Fcal ≤ F(1 – α, p – 1, n – p).
J
0:409ðN Þ þ 0:609
100kN=m
where F(1 – α, p – 1, n – p) is a value chosen from the
(9)
F table for the appropriate level of significance, α. If
the null hypothesis in the case of the F test is rejected, Likewise, the analysis of variance and the F statistic are
the contribution of individual variables to explain the presented in Table 3, while the t statistic is presented in
variation of dependent variable is carried out with the Table 4. The values of R2 and R2adj for the linear model
t test. If the results of t statistics show that any regression were 0.848 and 0.839, respectively.
coefficient is not statistically significant, a new equation Based on Figure 12–17, the following non-linear
must be proposed without the insignificant coefficients. regression model was chosen as the one that best repre-
The statistic t is calculated as follows: sents the behaviour of the data:
GEOMECHANICS AND GEOENGINEERING 15
Table 3. Analysis of variance and F statistic (linear model). Table 5. Analysis of variance and F statistic (non-linear model).
Sum of squares Mean square F Pr>F Sum of squares Mean square F Pr>F
Regression 50.3524 7.1932 101.80 <0.0001 Regression 14.9501 2.1357 111.30 <0.0001
Error 9.0448 0.0707 Error 2.4562 0.0192
Total 59.3972 Total 17.4063
Figure 19. BCR versus N for linear and non-linear models and different analytical methodologies, (a) Uniaxial geogrid, (b) Biaxial
geogrid.
and biaxial geogrid are presented. The values measured The additional bearing capacity provided by the lat-
in the laboratory and the values estimated with the eral forces produced by the equivalent depth can be
linear and non-linear models presented a more adjusted evaluated by the effect known as ‘wideslab’, according
behaviour to the models of Kumar and Saran (2003), to which, the foundation with equivalent depth has
Sharma et al. (2009) and Latha et al. (2013). Some a resistant base equal to the original foundation. The
differences between the results measured in the labora- ‘deep footing’ effect is explicit or implicit in the methods
tory and the linear and non-linear regression models are of Huang and Menq (1997) and Wayne et al. (1998).
evident with the methodologies of Huang and Menq According to the behaviour shown by the BCR with the
(1997) and Wayne et al. (1998). Latha et al. (2013) different parameters investigated, the laboratory results
obtained similar results and attributed the difference and the regression equations are in the range of values
presented with the models to the effect called ‘deep estimated by the different analytical and empirical
footing’, according to which, the behaviour of the rein- methodologies. Therefore, the two regression models
forced soil is very similar to the behaviour of the unrein- developed (Equations (9) and (11)) present an accepta-
forced soil with a depth of foundation equal to the depth ble approximation and allow estimating the bearing
of the reinforced zone. capacity of the reinforced soil.
GEOMECHANICS AND GEOENGINEERING 17
Shahin, H.M., et al., 2018. Effective use of geosynthetics to Tavangar, Y. and Shooshpasha, I., 2016. Experimental and
increase bearing capacity of shallow foundations. Canadian numerical study of bearing capacity and effect of specimen
Geotechnical Journal, 54 (12). doi:10.1139/cgj-2016-0505. size on uniform sand with medium density, reinforced with
Sharma, R., et al., 2009. Analytical modeling of geogrid rein- nonwoven geotextile. Arabian Journal for Science and
forced soil foundation. Geotextiles and Geomembranes, 27, Engineering, 41 (10), 4127–4137. doi:10.1007/s13369-016-
63–72. 2101-y.
Shukla, S.K., 2012. Handbook of geosynthetic engineering. Useche-Infante, D.J., et al., 2016. Shear strength behavior of
London: ICE. different geosynthetic reinforced soil structure from direct
Sitharam, T.G. and Sireesh, S., 2004. Model studies of shear test. International Journal of Geosynthetics and
embedded circular footing on geogrid-reinforced sand Ground Engineering, 2 (17), 1–16. doi:10.1007/s40891-
beds. Ground Improvement, 8 (2), 69–75. doi:10.1680/ 016-0058-2.
grim.2004.8.2.69. Vesic, A.S., 1973. Analysis of ultimate loads of shallow
Sitharam, T.G. and Sireesh, S., 2006. Effects of base geogrid on foundations. Journal of the Soil Mechanics and
geocell-reinforced foundation beds. Geomechanics and Foundations Division, ASCE, 99 (SM1), 45–73.
Geoengineering: an International Journal, 1 (3), 207–216. Wayne, M.H., Han, J., and Akins, K., 1998. The design of
doi:10.1080/17486020600900596. geosynthetic reinforced foundations. Proceedings of ASCEs
Skempton, A.W., 1951. The bearing capacity of clay. 1998 annual convention & exposition ASCE GSP, 76, 1–18.
Proceedings, Building Research Congress, 1, 180–189. Yadu, L. and Tripathi, R.K., 2013. Effect of the length of
Srinivasan, V. and Ghosh, P., 2013. Experimental investiga- geogrid layers in the bearing capacity ratio of geogrid
tion on interaction problem of two nearby circular footings reinforced granular fill-soft subgrade soil system. 2nd
on layered cohesionless soil. Geomechanics and Conference of Transportation Research Group of India
Geoengineering: An International Journal, 8 (2), 97–106. (2nd CTRG), 104, 225–234.