enhancements_in_applied_geomechanics_mining_and_excavation_simulation_and_analysis
enhancements_in_applied_geomechanics_mining_and_excavation_simulation_and_analysis
Adam Sevi
Jose Neves
Honghua Zhao Editors
Enhancements in Applied
Geomechanics, Mining,
and Excavation Simulation
and Analysis
Proceedings of the 5th GeoChina International
Conference 2018 – Civil Infrastructures
Confronting Severe Weathers and Climate
Changes: From Failure to Sustainability, held
on July 23 to 25, 2018 in HangZhou, China
Sustainable Civil Infrastructures
Editor-in-chief
Hany Farouk Shehata, Cairo, Egypt
Advisory Board
Khalid M. ElZahaby, Giza, Egypt
Dar Hao Chen, Austin, USA
Honghua Zhao
Editors
Enhancements in Applied
Geomechanics, Mining,
and Excavation Simulation
and Analysis
Proceedings of the 5th GeoChina International
Conference 2018 – Civil Infrastructures
Confronting Severe Weathers and Climate
Changes: From Failure to Sustainability, held
on July 23 to 25, 2018 in HangZhou, China
123
Editors
Adam Sevi Honghua Zhao
Department of Civil, Environmental, Department of Engineering Mechanics
and Construction Engineering Dalian University of Technology
Norwich University Dalian, China
Northfield, VT, USA
Jose Neves
Department of Civil Engineering,
Architecture and Georesources, IST
Universidade de Lisboa
Lisbon, Portugal
This Springer imprint is published by the registered company Springer Nature Switzerland AG
The registered company address is: Gewerbestrasse 11, 6330 Cham, Switzerland
Contents
v
vi Contents
vii
Magneto-Gravity Simulation of Cone
Penetration Test in Cohesionless Soil Under
Small Gravity Fields
1 Introduction
Near-earth space exploration extends from the Moon to the Mars, after the Apollo
program. The plans of return to the Moon have been proposed in the recent decade,
following the launched lunar missions, including European SMART-1, Chinese
Chang’E, U.S. Lunar Reconnaissance Orbiter, and Google Lunar XPRIZE. Studies
about the in situ resource utilization (ISRU) have processed the lunar materials for the
construction of lunar habitat structures (Meyers and Toutanji 2007; Nakamura and
Senior 2008). The colonization of the Moon starts with the lunar outposts, whereas the
excavation and the designed underground colonies require the understanding of
geotechnical properties of the lunar regolith.
Cone penetration test (CPT) acts as an important in situ testing tool for the site
characterization, owing to its economic advantages and the continuous measurements.
However, in the field of aerospace, the gravitational environment on a planetary body
varies with its mass; e.g. 1 g indicates the Earth gravity and 1/6 g represents the lunar
© Springer International Publishing AG, part of Springer Nature 2019
A. Sevi et al. (eds.), Enhancements in Applied Geomechanics, Mining,
and Excavation Simulation and Analysis, Sustainable Civil Infrastructures,
https://doi.org/10.1007/978-3-319-95645-9_1
2 P.-Q. Mo et al.
gravity. In situ testing is important and rather essential for the further space exploration,
without taking numerous samples back to earth for investigation. Therefore the inter-
pretation of CPT data is key to obtain the local soil stratigraphy and soil properties.
DEM simulation has been conducted to investigate the effects of gravity on the cor-
relations between CPT measurements and soil properties (e.g. Jiang and Wang 2013),
while the experimental study is still rather limited.
Physical modelling of geotechnical problems under small gravity fields, especially
for gravity smaller than the Earth gravity, is of interest to space and geotechnical
engineers. The parabolic flight aircraft, microgravity rocket or the drop tower were used
to create the microgravity environments, while they were not widely applied due to
their extremely high cost and the limited testing duration (Johnson et al. 1970; Sture
et al. 1998; Colwell and Taylor 1999; Thomas et al. 2000). Although the tilting method
(Tateyama 2007; Zou et al. 2015) and the hydraulic gradient similitude method
(Zelikson 1969; Dou and Byrne 1996) were adopted along with the additional body
force or water pressure, a reliable and economical testing method for small gravity
fields is still desired for studies of geotechnical problems.
This paper aims to provide results of the preliminary tests of CPT in cohesionless
soil under small gravity fields. The magneto-gravity modelling equipment is adopted to
simulate the small gravity fields with a magnetic acceleration against the earth gravity.
The cohesionless soil is used in the magneto-gravity fields to represent the earth, lunar
and martian regolith. The effect of gravity on the penetration resistance in cohesionless
soil is investigated by varying g-level from 1/6 to 4 g. Correlations between CPT data
and soil properties are then examined to provide insights into the penetration mecha-
nisms of CPT under small gravity fields, which could contribute to the in situ resource
utilization in the near-earth space exploration.
2 Testing Methodology
2
3
(b)
Penetrometer Aluminum tube
(a) Regolith simulant
50mm
150mm
Regolith simulant
Testing samples are prepared by pouring the cohesionless soil into a container with
inner diameter of D ¼ 86 mm. Every 25 mm height of soil with vibration is validated to
ensure the uniformity of the samples, until it reaches the height of H ¼ 150 mm
(Fig. 2b). An aluminum tube is then placed on the sample surface at the location of
penetration, and additional 50 mm of soil is poured around the tube, which aims to
provide surcharge avoiding the smear of soil surface caused by the attraction between
particles in the magnetic field. The void ratio of the prepared samples is about e ¼ 1:325,
that leads to the density q ¼ 1:26 g=cm3 and the relative density Dr ¼ 22:7% .
(a) (b)
The CPT tests are conducted under different gravity fields, including 1/6 g (the
lunar gravity), 1/3, 2/3, 1, 2, 4 g. The g-level Ng , defined as:
g gM gsimu
Ng ¼ ¼ ð1Þ
g g
(a) (b)
The effect of gravity is presented in Fig. 4b, with the variation of g-level Ng ranging
from 1/6 to 4. Linear increase of the cone tip resistance is again observed for a shallow
penetration. It is obvious that qc increases with g-level under small gravity fields at a
certain depth, where the in situ stress is higher for a larger g-level. However, the
increase seems to be nonlinear, and the effect of gravity is more significant for g-level
smaller than 1 g.
In view of the observed linear relationship between the cone tip resistance and the
depth (i.e. the vertical effective stress), it is reasonable to apply the normalization of
CPT data following Bolton et al. (1993), as defined in Eqs. (2) and (3). This nor-
malization is also consistent with the definition of cone factor in sand for bearing
capacity analysis, where Nq ¼ qc =rv ¼ Q þ 1.
qc rv
Q¼ ð2Þ
rv
z
Z¼ ð3Þ
B
The results of normalised penetration resistance are provided in Fig. 5a. Note that
the initial vertical stress is a combination from both the soil gravity and the surcharge.
A comprehensive constant value of Q is obtained for the normalised penetration depth
between 6 and 14. Accordingly, Q decreases with the g-level, as also reported by
Bolton et al. (1999) with series of centrifuge tests.
1/6g
4 4 1/3g
Z=6 2/3g
6 6 1g
2g
8 8 4g
10 10
12 12
Z=14
14 14
2/3g
16 4g 1/3g 1/6g 16
2g 1g
18 18
ðqc rv Þ=r0
qc1N ¼ ð4Þ
ðrv =r0 Þ0:5
where r0 is the reference stress = 100 kPa. The corresponding results are shown in
Fig. 5b. Although the decrease of qc1N is also found, it can be concluded that Q is a
more appropriate normalization to achieve a constant value for penetration under a
small gravity field.
4.5
Normalised penetration resistance: Q
3.5
100 Q / Q1g = (Ng)-0.83 R2 = 0.995
3
Q / Q1g
2.5
2
50 Prediction of Jiang & Wang (2013)
1.5
0.5
0 0
0 1 2 3 4 5 0 1 2 3 4 5
g-level: Ng g-level: Ng
Note that this expression is obtained based on the loose sand tests for small gravity
fields (Ng \5), and comprehensive verification is suggested when it is applied for larger
magnitude of g-level. The results are also compared with the prediction of Jiang and
Wang (2013), which was deduced from two-dimensional DEM simulation of CPT
under different gravity fields. The dramatic decrease of Q=Q1g for Ng \1 is compar-
ative to the experimental data in this paper, whereas the prediction of Q was much
higher for DEM simulation of large g-level tests.
8 P.-Q. Mo et al.
(a) 0 (b) 0
2 2
Correlation from Tatsuoka et al. (1990):
Normalised penetration depth: Z=z/B
Normalised penetration depth: Z=z/B
12 12
14 14
16 16
4g 2g 1g 2/3g 1/3g 1/6g Dr=22.7%
18 18
-20 0 20 40 60 80 100 -20 0 20 40 60 80 100
Dr (%) Dr (%)
Fig. 7. Back analysis of relative density: a correlation of Tatsuoka et al. (1990); b modified
correlation
The modified results for the back analysis of Dr based on the cone tip resistance are
thus found with improved predictions, in consistence with the measured data for dif-
ferent gravity fields (see Fig. 7b).
Magneto-Gravity Simulation of Cone Penetration Test … 9
4 Conclusions
This paper presents the results of preliminary tests of cone penetration testing in
cohesionless soil under small gravity fields ranging from 1/6 to 4 g. The gravity varies
to simulate the gravitational environment on a planetary body; e.g. 1 g indicates the
Earth gravity and 1/6 g represents the Moon gravity. The effects of gravity on the
correlations between CPT data and soil properties are investigated by applying a
magnetic acceleration against the earth gravity using the magneto-gravity modelling
system. Accordingly, the cohesionless soil is made from the cement powder and Fe3O4
magnetic powder with a rational particle distribution, simulating the planetary regolith.
A miniature penetrometer with diameter of 8 mm is manufactured along with the
actuator proving a penetration speed at around 0.55 mm/s.
From the results of penetration resistance, linear increase of qc is observed for a
shallow penetration. It is obvious that the cone tip resistance increases with g-level
under small gravity fields at a certain depth, where the in situ stress is higher for a
larger g-level. Q is found to be a more appropriate normalization to achieve a constant
value for penetration under a small gravity field. Additionally, the normalised pene-
tration resistance decreases exponentially with the g-level, and a relationship between
Q=Q1g and Ng has been proposed with comparison of DEM simulation. The empirical
correlation between CPT measurements and relative density has been modified to
include the effects of gravity, and the improved predictions present relatively consistent
results with the measured data for different gravity fields. This study examined the
influence of planetary gravity to the interpretation of in situ tests, which would con-
tribute to the understanding of the stratigraphy and properties of the planetary regolith
for further in situ resource utilization and the construction of space outposts.
Acknowledgements. The author would like to acknowledge financial support from the National
Natural Science Foundation of China (no. 51323004), 111 Project (B14021) and Natural Science
Foundation of Jiangsu Province (no. BK20170279).
References
Baldi, G., Bellotti, R., Ghionna, V.N., Jamiolkowski, M., Pasqualini, E.: Interpretation of CPTs
and CPTUs; 2nd part: drained penetration of sand. In: Proceedings of the 4th International
Geotechnical Seminar, Singapore, pp. 143–156 (1986)
Bolton, M.D., Gui, M.W.: The study of relative density and boundary effects for cone penetration
tests in centrifuge. Technical Report CUED/DSOILS/TR256, University of Cambridge
(1993)
Bolton, M.D., Gui, M.W., Phillips, R.: Review of miniature soil probes for model tests. In:
Proceedings of 11th South East Asia Geotechnical Conference, pp. 85–91 (1993)
Bolton, M.D., Gui, M.W., Garnier, J., Corte, J.F., Bagge, G., Laue, J., Renzi, R.: Centrifuge cone
penetration tests in sand. Géotechnique 49(4), 543–552 (1999)
Colwell, J.E., Taylor, M.: Low-velocity microgravity impact experiments into simulated regolith.
Icarus 138(2), 241–248 (1999)
Dou, H., Byrne, P.M.: Dynamic response of single piles and soil pile interaction. Can. Geotech.
J. 33(1), 80–96 (1996)
10 P.-Q. Mo et al.
Gui, M., Bolton, M.D., Garnier, J., Corte, J.F., Bagge, G., Laue, J., Renzi, R.: Guidelines for
cone penetration tests in sand. In: International Conference on Centrifuge Modelling
(Centrifuge’98), Balkema, vol. 1, pp. 155–160 (1998)
Jamiolkowski, M., Lo Presti, D.C.F., Manassero, M.: Evaluation of relative density and shear
strength of sands from CPT and DMT. Soil Behav. Soft Ground Constr. 7(119), 201–238
(2003)
Jiang, M.J., Wang, X.X.: Numerical analysis of cone penetration tests under different gravity
fields by distinct element method. Rock Soil Mech. 34(3), 863–873 (2013)
Johnson, S.W., Pyrz, A.P., Lee, D.G.: Simulating the effect of gravitational field and atmosphere
on behavior of granular media. J. Spacecr. Rocket 7(11), 1311–1317 (1970)
Lancellotta, R.: Analisi di affidabilita in ingegneria geotechnica. Technical Report 625,
Politernico di Torino (1983)
Li, R.L., Zhou, G.Q., Chen, G.Z., et al.: Development of a new geo-mechanical magnetic sensing
soils similar material. Funct. Mater. 19(43), 2620–2623 (2012)
Meyers, C., Toutanji, H.: Analysis of lunar-habitat structure using waterless concrete and tension
glass fibers. J. Aerosp. Eng. 20(4), 220–226 (2007)
Mo, P.Q., Marshall, A.M., Yu, H.S.: Centrifuge modelling of cone penetration tests in layered
soils. Géotechnique 65(6), 468–481 (2015)
Nakamura, T., Senior, C.L.: Solar thermal power for lunar materials processing. J. Aerosp. Eng.
21(2), 91–101 (2008)
Robertson, P.K., Wride, C.E.: Evaluating cyclic liquefaction potential using the cone penetration
test. Can. Geotech. J. 35(3), 442–459 (1998)
Sture, S., Costes, N.C., Batiste, S.N., et al.: Mechanics of granular materials at low effective
stresses. J. Aerospace Eng. 11(3), 67–72 (1998)
Tatsuoka, F., Zhou, S., Sato, T., Shibuya, S.: Evaluation method of liquefaction potential and its
application. Technical report, Ministry of Education of Japan (1990)
Thomas, V.A., Prasad, N.S., Reddy, C.A.M.: Microgravity research platforms—a study. Cur.
Sci. 79(3), 336–340 (2000)
Xu, X.T.: Investigation of the end bearing performance of displacement piles in sand. Ph.D.
thesis, The University of Western Australia (2007)
Zelikson, A.: Geotechnical models using the hydraulic gradient similarity method. Géotechnique
19(4), 495–508 (1969)
Zhou, G.Q., Li, R.L., Chen, G.Z., Zhao, X.D., Wang, J.Z., Liang, H.C., Lin, C., Lai, Z.J., Chen,
X.: Geotechnical magnetic-similitude-gravity model testing method. Unpublished report
(2016)
The Effect of Several Parameters
on the Behavior of Asphalt Mixture in Libya
Abstract. In Libya, many steps have been taken to extend the life of the
pavement and improve the design of asphalt mixture to have a better connection
to Libyan climate. This study is to evaluate the effect of some parameters on the
performance of the asphalt mixture in hot arid weather as Libya. Two asphalt
concrete mixtures were designed using two different binders, Performance
Grade binder (PG70-10) and Pen Grade System (B60/70). These mixtures were
formulated in accordance with the Marshall mix design method and the
Superpave mix design method. Rutting test also tested these mixtures to eval-
uate the performance of the bituminous mixtures under certain conditional
parameters. The study clearly showed that the asphalt mixtures with binder
PG70-10 are slightly performed better than the mixtures with binder (B60/70) In
addition, the results obtained from Superpave mix design test have shown a
better performance indicator than those obtained with the Marshall Stability tests
and they were superior at all stages of tests. These results provide a guide line to
improving the production of the hot mix asphalt (HMA) and a foundation for
different paving mixes, by adopting PG binder and Superpave Gyratory com-
pactor (SGC) instead of Marshall hammer and binder B60/70.
1 Introduction
Asphalt mixture properties such as density and air voids are largely dependent on the
type of the compaction method, selection of materials and method of mix design. These
properties, in turn, have an influence over the pavement performance indicators such as
rutting, fatigue, potholes and alligator cracks (Malunga et al. 2014). Therefore, this
research aims to improve the local mix design formulation of the pavement asphalt
mixture to obtain an economical mixture that would meet the requirements for the
characteristics of the pavement in hot weather. The Marshall Mix Design Method
selected the asphalt binder based on the empirical Penetration Grading Pen Grade
System. It does not mandate for viscosity measurements and hence, it is difficult to
establish the correct HMA. Moreover, temperature can be a major contributor to several
types of distress (Salem et al. 2014). Binder and compaction method are perhaps the
most important elements for an asphalt mixture design. They influence almost all the
important properties of HMA mixtures such as stiffness, stability, durability, etc.
Therefore, the following research will adopt a new asphalt binder PG which should be
related to the climatic conditions of Libya and a new method of compaction which is
SGC. Marshall Hammer Compactor is very old method of compaction in the labora-
tory. Nevertheless, the variation in the methods of laboratory compaction is not only
the result of the methods of assessment, but is due to the technology used (Foster
1982).
2 Objectives
This research aims to show the effect of some factors on the asphalt mix design such as
the type of bitumen binder, skeleton of aggregate and the method of compaction. These
mixtures were formulated in accordance with the Marshall mix design method and the
Superpave mix design method. Rutting test also tested these mixtures to evaluate the
performance of the bituminous mixtures under certain conditional parameters. Addi-
tionally, the purpose of this paper is to show the limitation of Marshall Mix Design
Method, which is still used today and to increase the implementation of earlier
researches.
Hubbard-Field Method might be considered as the first formal design method for asphalt
mixtures. It was originally developed to design sand-asphalt mixtures and later modified
for aggregates (Roberte et al. 2002). The Marshall mix design method, initially devel-
oped by Bruce G. Marshall from the Mississippi Highway Department in 1939, and later
amended by the US Waterways Experiment Station of the US Army to include defor-
mation measurements. This is based on two criteria, a minimum stability and a minimum
range of flow values. However, it was based on a single piece of equipment to design and
control the asphalt paving mixtures which referee to as the Marshall Test. The Superpave
mix design method, for Superior Performing Asphalt Pavement System, was developed
in the 1980s as a part of the Strategic Highway Research Program (SHRP). Superpave
addresses aggregate selection and asphalt binder selection in the perspective of future
traffic and climatic requirements. It selects aggregates based on gradation distribution,
angularity, clay content, water absorption, abrasion, soundness. The asphalt binder
selection relies on the rutting analyser and the dynamic (AASHTO 1993). In Libya, mix
design of asphalt mixture is based on the Marshall method which is an empirical design
method that does not replicate the compaction, materials, and climate prevailing in
Libya. Previous studies have shown that the compaction by impact in the Marshall
method is unrealistic compared to the SGC, which simulates the field density. And the
Pen Grade system is typically conducted at a single test temperature of 25 °C, which
does not account for the entire temperature spectrum to which the asphalt-binders are
subjected to in the field (Almadwi and Assaf 2017). However, the term “Superpave”
refers to more than just the computer program. The system includes test equipment, test
methods, and criteria. Most important, it represents an improved system for specifying
component materials, asphalt mixture design and analysis, and pavement performance
prediction (Asi and Khalayleh 2011).
The Effect of Several Parameters on the Behavior of Asphalt … 13
4 Methodology
Two asphalt concrete mixtures of aggregate-asphalt binder were designed using two
different binders, PG70-10 and B60/70. These mixtures were formulated in accordance
with the Marshall mix design method and the Superpave mix design method. Rutting
test also tested these mixtures to obtain an economical mixture that would meet the
requirements for the characteristics of the pavement in hot weather. These specimens
were used for a comparison under different Libyan weather conditions. Evidently, the
test mixtures contained several asphalt bitumen contents both above and below the
optimum content of asphalt.
asphalt-binder selection criterion known as the Performance Grading system has been
developed by the Strategic Highway Research Program (SHRP) to overcome of some
the asphalt pavement distresses such as rutting, fatigue cracking, and thermal cracking
The PG concept is based on the rheological theory that an asphalt-binder’s property
should be related to the climatic conditions under which it is to be used. Engineering
properties believed to be related to the expected performance (maximum 7-day pave-
ment temperature, 1-day minimum pavement temperature, loading duration based on
truck speed, and traffic volume) are featured in the PG binder specifications to allow for
a proper selection of asphalt roads (Wahhab et al. 1997).
Table 1. Sieve analysis result and aggregate gradation for laboratory mix design
Sieve size (mm) % lower limit % upper limit % passing
28 100 100 100
20 95 100 98.7
14 67 90 85.6
10 52 75 68.1
5 34 55 44.87
2.5 24 45 27.6
1.25 16 39 17.7
0.63 9 31 12.7
0.31 6 23 9.68
0.16 4 15 7.55
0.08 3 8 6.33
The Effect of Several Parameters on the Behavior of Asphalt … 15
Table 2. Specific gravity for the compound aggregates and asphalt binders
Properties of materials in mix design
Agg. size and % of Agg. in Bulk specific Compound mix % water
asphalt binder mixture gravity specific gravity absorption
20–14 (mm) 15 2.739 2.716 0.630
14–10 (mm) 15 2.736 0.670
10–5 (mm) 26 2.733 0.700
Crushed 0–5 36 2.698 0.005
(mm)
Sand 0–5 (mm) 3 2.595 0.007
Filler (mm) 5 2.700 0.006
Bitumen (B60/70) 4.50 1.020 1.020 –
Bitumen PG 70-10 4.25 1.290 1.029 –
100
90
80
70
60
% Passing
50
40
30
20
10
0
0.01 0.1 1 10 100
Particle Diameter (mm)
Lower Limit Upper Limit Final Mix
5.25
S1 PG 70-10 S2 B 60/70
Fig. 4. Rutting values after 30,000 cycles with two different binders
7.00
RUTTING DEPTH (MM)
6.00
5.00
4.00
3.00
2.00
1.00
-
0 5000 10000 15000 20000 25000 30000
NUMBER OF CYCLES
SI PG 70-10 SII PG 70-10 SI B(60/70) SII B (60/70)
The asphalt binder that was used in the Marshall method may not be sufficient to ensure
that selected asphalt binders can satisfactorily meet the paving service temperatures
experienced in Libya. The method of selecting the binder and the bituminous content in
the mixture is still questionable. However, in this study, the conventional local spec-
ifications of selecting and blending the aggregate and binder that were used in the
Marshall mix design did not perform well and thus yielded inadequate results. On the
other hand, the Marshall mixes were designed using the SGC, PG asphalt-binders and a
proper aggregate gradation indicated that the Marshall method performed like the
Superpave method. Asphalt-binder type, asphalt-binder content, aggregate properties,
and gradation showed a significant effect on the performance of the mixes as expected.
Tables 3, 4 and 5 illustrated all the laboratory experiment results. Figures 2, 3, 4, and 5,
show the behavior of asphalt mixtures with new binder PG70-10 and a new compaction
method under the hot weather condition such as in Libya. It can be seen from Figs. 2,
3, and Tables 3, 4 that the volumetric properties values of the Superpave mixes are
properly better than that of Marshall mixtures. The flow values of Marshall mixtures
with PG70-10 are slightly less than that of Marshall mixes with B60/70. This could be
due to the differences in the binder properties and in the compaction techniques used.
The SGC rotates at a constant rate during the compaction, and this characteristic
provides around a better orientation of aggregate particles and aggregate interlock. This
process simulates closely to the field compaction. In Figs. 4, 5 and Table 5, asphalt
blends using the asphalt binder PG70-10 showed better resistance to rutting than the
mixes using the asphalt binder B60/70. This study showed that in rutting test, most of
the permanent deformation occurs in the upper four inches of HMA. The study clearly
showed that the amount of air voids in the mixture was one of the most important
properties. Although the number of air voids had a significant effect on the pavement
deformation, this parameter could not be used to predict the deformation. When air
voids were low, as a rule, there was considerable of deformation, and when air voids
were high, it was much less likely of deformation. It is therefore important that air
voids be carefully controlled during construction. In the future, the air voids will be
determined from samples compacted with the SGC. Regardless of the type of com-
paction used, it is important that the sample is compacted to about the same density as
in the field after several years of traffic.
The proper selection of the aggregates, method of compaction and the asphalt binder
can improve the pavement performance. However, the most important factors which
affects the pavement performance are the compaction technique, type of asphalt binder
and selection of aggregate. In this research, two asphalt concrete mixes were designed
using two different binders. These mixtures were formulated in accordance with the
Marshall mix design method and the Superpave mix design method and these mixtures
20 K. E. Atrash and G. J. Assaf
also tested with Rutting test. The analysis was conducted to evaluate HMA properties
such as deformation, durability, air voids, voids in mineral aggregate, voids filled with
asphalt, etc. The study in which the Marshall mixes were designed using the SGC, PG
asphalt-binders and a proper aggregate gradation (source and consensus requirements),
indicated that the Marshall method performed like the Superpave method. The type of
asphalt-binder, the asphalt-binder content, and the aggregate properties, showed a
significant effect on the performance of the mixes, as theoretically expected. A virtual
need to evaluate the mix-design methods on the hot weather locale is still mandatory. It
therefore can be concluded that the Marshall method can perform equally well if the
method of compaction and asphalt-binder selection criteria are updated in Libya.
Acknowledgements. It is with great joy and lightness of spirit that I offer my deepest regards to
my friends who have assisted and supported me in countless ways as I journeyed through the
process of undertaking, creating, and completing this paper.
References
Almadwi, F.S., Assaf, G.J.: Performance testing of paving mixes for Libya’s hot and arid
conditions, using Marshall stability and SUPERPAVE gyratory compactor methods. In:
International Congress and Exhibition “Sustainable Civil Infrastructures: Innovative Infras-
tructure Geotechnology”, pp. 313–323. Springer, Cham (2017)
American Association of State Highway & Transportation Officials: AASHTO Guide for Design
of Pavement Structures, 1993, vol. 1. AASHTO (1993)
Asi, I.M., Khalayleh, Y.: Adaptation of Superpave asphalt concrete mix design procedure to
jordan climatic and traffic conditions. Int. J. Pavement Res. Technol. 4(3), 154–161 (2011)
Button, J.W., Little, D.N., Jagadam, V., Pendleton, O.J.: Correlation of Selected Laboratory
Compaction Methods with Field Compaction. Transportation Research Record, 1454 (1994)
Cominsky, R.J., Huber, G.A., Kennedy, T.W., Anderson, M.: The Superpave Mix Design
Manual for New Construction and Overlays, No. SHRP-A-407. Strategic Highway Research
Program, Washington, DC (1994)
Foster, C.R.: Development of Marshall Procedures for Designing Asphalt Paving Mixtures,
No. IS 84. National Asphalt Pavement Association (1982)
Kandhal, P.S., Foo, K.Y., Mallick, R.B.: A Critical Review of VMA Requirements in Superpave,
No. NCAT Report No. 98-1. National Center for Asphalt Technology (1998)
Kennedy, T.W., Huber, G.A., Harrigan, E.T., Cominsky, R.J., Hughes, C.S., Von Quintus, H.,
Moulthrop, J.S.: Superior Performing Asphalt Pavements (Superpave): The Product of the
SHRP Asphalt Research Program, No. SHRP-A-410. Strategic Highway Research Program,
National Research Council, Washington, DC, USA (1994)
Lanham, M.D.: Development of Marshall Procedures for Designing Asphalt Paving Mixtures,
Information Series 84. Natl Asphalt Pavement Assoc. (1982)
Malunga, R., El Halim, A.A., Muya, M., Mwanza, A.D., Faruk, A.N., Mushota, C., Walubita, L.
F.: Theoretical review of different asphalt mix-design methods and their applicability for
developing countries like Zambia. In: Design, Analysis, and Asphalt Material Characteri-
zation for Road and Airfield Pavements, pp. 210–217 (2014). https://doi.org/10.1061/
9780784478462.026
Roberte, F.L., Mohammad, L.N., Wang, L.B.: History of hot mix asphalt mixture design in
United States. J. Mater. Civ. Eng. 14(4), 279–293 (2002)
The Effect of Several Parameters on the Behavior of Asphalt … 21
Salem, H.A., Uzelac, D., Matic, B.: Temperature zoning of Libya desert for asphalt mix design.
In: Applied Mechanics and Materials, vol. 638, pp. 1414–1426. Trans Tech Publications
(2014). https://doi.org/10.4028/www.scientific.net/AMM.638-640.1414
Wahhab, H.I.A.A., Asi, I.M., Al-Dubabe, I.A., Ali, M.F.: Development of performance-based
bitumen specifications for the Gulf countries. Constr. Build. Mater. 11(1), 15–22 (1997)
Review and Analysis on Using the Analytical
Approaches for Predicting the Pavement
Performance
1 Introduction
The root causes of the poor performing pavements or earlier deterioration of roads in
the form of potholes, undulations, fatigue cracking depends on the local factors and
may differ to the pavement to pavement. Finding an empirical solution for its main-
tenance may result again in its premature failure or the solution may be highly not
realistic for the particular condition. A number of researchers including Tapase and
Ranadive (2016, 2017) have already reported that earlier deterioration of roads and its
time-consuming rehabilitation/maintenance operation are the direct or indirect cause
behind an increase in road fatalities at around 3% compared to a preceding year in
India. The World Health Organisation (WHO), Ministry of Road Transport and
Highways (MoRTH) have already started taking actions to mitigate the fatalities and to
strengthen road safety as well as improve ease of transport across the various parts
within the country and worldwide. Yang et al. (2015) have observed number of
overlayed sections of the Arkansas Interstate roads have experienced severe premature
rutting and cracking and reported the need of pavement management tool for pavement
maintenance actions.
It is difficult to encounter the existing problem of the deteriorated road in a well-
timed and economical manner simultaneously, due to various complexities in the
maintenance procedure. The properties of existing layers are needed to be evaluated as
well as the complex interaction of different possible strategies needs to be forecasted
for planning the maintenance, rehabilitation, and reconstruction of the deteriorated
pavement. So, applying various trials at the network level with the help of analytical
tool for predicting the performance of pavement before construction of new pavement
or for maintenance or rehabilitation of deteriorated pavement will be helpful for
selecting the best possible optimum strategies for providing, evaluating and main-
taining pavements in a serviceable condition for the longer duration. In this connection,
it should be noted that versatile finite element method and the artificial neural network
holds a bright assurance. Present paper is a state of the art review focusing the
advantageous application of various methods like finite element method (FEM) and the
artificial neural network (ANN) in the civil infrastructures from failure to sustainability
and sequentially a comprehensive sensitivity analysis method is proposed as reported in
Fig. 1 to explore the influence of analytical analysis in design and planning the best
possible optimum strategies for providing, evaluating and maintaining pavements in a
serviceable condition for the longer duration.
2 Literature Review
In recent years considerable analytical research efforts have been expended in civil
engineering field due to the versatile methods like FEM, ANN and so on. The
development of the mechanistic-empirical pavement design which is applied by a
number of countries for the design of pavement is more realistic and has offered a
potential for using more improved material characterization and constitutive models of
both flexible and rigid pavements. A number of researchers have already documented
the usefulness of FEM and ANN in predicting the performance of flexible as well as the
rigid pavement. The Mechanistic Empirical Pavement Design Guide recommends
JULEA a multilayer elastic program to compute flexible pavement responses (NCHRP
2007). The recently released IRC: 37-2012 guidelines for the design of flexible
pavement recommends using the IITPAVE for layered system analysis.
Gogoi et al. (2013) studied to investigate if there is any co-relation in rutting and
fatigue distress. Tapase and Ranadive (2016), reported the usefulness of two dimen-
sional finite element analyses to study the effect of variation in thickness of different
component layers on the critical parameters. Zheng et al. (2017) have developed
genetic algorithm (GA)-back propagation (BP) neural network model based on the
experimental data for evaluating the relationship between long-term skid resistance of
epoxy asphalt mixture with the various engineering parameters for the running vehi-
cles. Xiao and Amirkhanian (2009) as stated that the ANN model had higher accuracy
then regression-based model for estimating the stiffness of rubberized asphalt concrete.
Shafabakhsh et al. (2015) have reported the effectiveness of ANN method and its
performance in evaluating the strengths of the hot mix asphalt mixtures. Ghanizadeh
and Fakhri (2017) reported a nonlinear regression equation for determining equivalent
frequencies at different depths of asphalt layer by employing a feed-forward Artificial
Neural Network (ANN). Leiva-Villacorta et al. (2017a, b) developed ANN models
capable of predicting pavement layer moduli rapidly and reliably. A significant
decrease in error by ANN predicted moduli was noticed to that of conventional iterative
approach while back-calculating the layer moduli. Rezaei-Tarahomi et al. (2017)
evaluated the sensitivity of NIKE3D-FAA software to compute critical pavement
responses respect to top-down and bottom-up cracking of concrete airport pavement
structures at different locations (interior, corner, and the edge of the slab) by varying
other NIKE3D-FAA inputs. Numerous explicit inputs to finite element analysis that
need to be considered in developing the substitute stress response prediction model are
evaluated using ANN in the study. Sensitivity evaluations were performed using a
normalized sensitivity index (NSI) as the quantitative metric which has proven the
effectiveness of finite element analysis and ANN approach. Mazari and Rodriguez
(2016) have reported the use of a hybrid gene expression programming-neural network
technique for predicting the international roughness index (IRI) which represents the
pavement roughness. The focus of the reported study is to couple genetic programming
and artificial neural network for IRI prediction on a dataset collected from the LTPP
database. The developed algorithm showed reasonable performance for prediction of
IRI using traffic parameters and structural properties of pavement. Yu et al. (2018)
study reports the application of regression model and finite element modeling for
Review and Analysis on Using the Analytical Approaches for … 25
The artificial neural network is a method which has an ability to solve the problem by
information obtained from the past experience. The configuration of ANN is done by
interconnecting the computational elements, through a learning process. There are three
basic building blocks of the artificial neural network which are as network architecture,
setting the weights and the activation function. The arrangement of neurons into layers
and the pattern of connection within and in between layer are generally called as the
architecture of the net. The neurons within a layer are found to be fully interconnected
or not interconnected. The number of layers in the net can be defined to be the number
of layers of weighted interconnected links between the particular slabs of neurons. If
two layers of interconnected weights are present, then it is found to have hidden layers.
There are various types of network architectures: feedforward, feedback, fully inter-
connected net, competitive net, etc. The method of settings the values for the weights
enables the process of learning or training. The process of modifying the weights
connections between network layers with the objective of achieving the expected
output is called training a network. The activation function is used to calculate the
output response of a neuron. The sum of the weighted input signals is applied with an
activation to obtain the response. For neurons in the same layer, same activation
functions are used. There may be linear as well as nonlinear activation functions. The
nonlinear activation functions are used in a multilayer net.
To assist decision makers in finding optimum strategies for maintaining pavement
in a serviceable condition nowadays pavement management system plays a vital role in
providing a cost-effective approach. The future performance of pavement is evaluated
in terms of the expected impact of maintenance and rehabilitation treatments that are
proposed at network level activities of PMS (AASHTO guidelines, Huang). Network
level and project level are the two main activities of PMS wherein the network level
addresses the overall budget and planning services including maintenance and reha-
bilitation and the project level is where the specific decision on maintenance strategies
and funding allocation are made. So it is learned that pavement performance if pre-
dicted at the planning stage itself or in other words at the network level then it would be
convenient to propose the cost-effective maintenance budget at the network level which
will help for assigning the most cost-effective strategies at the project level.
An analysis of the existing pavement to determine/predict the cause of deterioration
is necessary to identify and to decide the rehabilitation or maintenance alternative
available. It should be a full proof strategies which should adequately address the cause
of deterioration whether it is normal or abnormal, whether the material used is durable,
adequacy in drainage, the environmental factor, present and future traffic and the
geometrical factor.
There is a large number of maintenance, rehabilitation, and construction alternative
available for both flexible and rigid pavements, and recycling has increased that
number of options. Surface seals such as aggregate seal and slurry seal combined with
26 R. Tapase et al.
localized repairs are often used as preventive maintenance treatment for flexible
pavement, and they are also used as rehabilitation treatments on lower volume flexible
pavement when structural improvements are not required. Overlays have become more
versatile by combining them with inter-layers such as fabric, milling full or partial
widths, and recycling part of an existing layer prior to applying the overlay. Other types
of rehabilitation and reconstruction for flexible pavements include the following:
1. Cold in-place recycling followed by a new surface;
2. Hot in-place recycling with or without an overlay;
3. Reworking and stabilization of foundation materials followed a new surface;
4. Partial or full depth removal and replacement;
5. Full depth recycling; and
6. Overlay with PCC.
Another type of rehabilitation for rigid pavement may include
1. A series of maintenance treatments such as partial depth patching, full depth repair,
surface grinding and joint sealing often referred to as concrete pavement restoration
(CPR)
2. Break and seat with a concrete/asphalt overlay and
3. Bonded or unbounded PCC overlay.
Pavement design is the project level activity where detailed engineering and eco-
nomic considerations are given to alternate combinations of sub-base base, and surface
materials which will provide adequate load carrying capacity. Factors that are con-
sidered include; material, traffic, climate, maintenance drainage, and life-cycle costs.
The number of researchers have reported the usefulness of analytical tool in pavement
design and its analysis; more rigorous pavement design can be achieved by the
application of analytical tool which makes the pavement engineer check variety of
combination of material and thickness combinations for the particular location.
Pavement performance if predicted at the planning stage itself or in other words at
the network level then it would be convenient to propose the cost-effective maintenance
budget at the network level which will help for assigning the most cost-effective
strategies at the project level. Such tool will accommodate not only the thickness
design but will assist the decision makers in finding optimum strategies for providing,
evaluating and maintaining pavements in a serviceable condition for the longer
duration.
To evaluate existing structural conditions and to predict the pavement condition
after ‘n’ number of years, a methodology is proposed with the help of artificial neural
networks (ANNs), which will use the collected data from NDT tests like falling weight
deflectometer (FWD) and will interpret accordingly. The structural properties of the
existing pavements can be back-calculated using the artificial neural network. A neural
networks methodology as mentioned in Fig. 1 will be used for back-calculation model
proposed in this research. With an intention, it is proposed to develop specific models
Review and Analysis on Using the Analytical Approaches for … 27
References
Abo-Hashema, M.: Artificial neural network approach for overlay design of flexible pavements.
Int. Arab J. Inf. Technol. 6(2) (2009)
American Association of State Highway Officials AASHTO: AASHTO guide for the design of
pavement structures. AASHTO, Washington, D.C. (1993)
Diefenderfer, B.K., Al-Qadi, I.L., Diefenderfer, S.D.: Model to predict pavement temperature
profile: development and validation. J. Transp. Eng. ASCE 132(2), 162–167 (2006)
Ghanizadeh, A.R., Fakhri, M.: Quasi-static analysis of flexible pavements based on predicted
frequencies using Fast Fourier Transform and Artificial Neural Network. Int. J. Pavement Res.
Technol. (2017). https://doi.org/10.1016/j.ijprt.2017.09.002
Gogoi, R., Das, A., Chakraborty, P.: Are fatigue and rutting distress modes related? Int.
J. Pavement Res. Technol. 6(4), 269–273 (2013)
IRC: 37-2012. Guidelines for the design of flexible pavements. Indian Roads Congress, New
Delhi (2008)
Kamboozia, N., Ziari, H., Behbahani, H.: Artificial neural networks approach to predicting rut
depth of asphalt concrete by using of visco-elastic parameters. Constr. Build. Mater. 158
(2018), 873–882 (2018)
Leiva-Villacorta, F., Vargas-Nordcbeck, A., Timm, D.H.: Non-destructive evaluation of
sustainable pavement technologies using artificial neural networks. Int. J. Pavement Res.
Technol. 10(2), 139–147 (2017)
Masad, S., Little, D., Masad, E.: Analysis of flexible pavement response and performance using
isotropic and anisotropic material properties. J. Transp. Eng. ASCE 132(4), 342–349 (2006)
Mazari, M., Rodriguez, D.D.: Prediction of pavement roughness using a hybrid gene expression
programming-neural network technique. J. Traffic Transp. Eng. 3(5), 448–455 (2016)
Moroval, N., Serin, S., Terzi, S., Saltan, M.: Prediction of the pavement serviceability ratio of
rigid highway pavements by artificial neural networks. J. Adv. Technol. Sci. 2(1), 12–25
(2013)
NCHRP: Evaluation of mechanistic-empirical design procedure. National Cooperative Highway
Research Program, NCHRP Project 1-37A, National Research Council, Washington, D.C.
(2007)
28 R. Tapase et al.
Ranadive, M.S., Tapase, A.B.: Investigation of behavioral aspects of flexible pavement under
various conditions by finite element method. In: Yang, Q., Zhang, J.-M., Zheng, H., Yao, Y.
(eds.) Constitutive modeling of geomaterials. Springer, Berlin, pp. 765–770 (2013). https://
doi.org/10.1007/978-3-642-32814-5_100
Ranadive, M.S., Tapase, A.B.: Parameter sensitive analysis of flexible pavement. Int. J. Pavement
Res. Technol. (IJPRT). Elsevier, Special Issue on Sustainability on Pavement Engineering
(2016a). https://doi.org/10.1016/j.ijprt.2016.12.001
Ranadive, M.S., Tapase, A.: Pavement performance evaluation for different combinations of
temperature conditions and bituminous mixes. Innovative Infrastruct. Solutions 1(40), 1–5
(2016b). https://doi.org/10.1007/s41062-016-0040-9. (Springer)
Rezaei-Tarahomi, A., Kaya, O., Ceylan, H., Gopalakrishnan, K., Kim, S., Brill, D.R.: Sensitivity
quantification of airport concrete pavement stress responses associated with top-down and
bottom-up cracking. Int. J. Pavement Res. Technol. 10, 410–420 (2017)
Shafabakhsh, G.H., Ani, O.J., Talebsafa, M.: Artificial neural network modeling (ANN) for
predicting the performance of nano-modified hot-mix asphalt mixtures containing steel slag
aggregates. Constr. Build. Mater. 85, 136–143 (2015)
Tapase, A., Ranadive, M.: Performance evaluation of flexible pavement using finite element
method. ASCE GSP 266, Geo-China 2016: Material, Design, Construction, Maintenance, and
testing of pavement, 9–17 (2016). https://doi.org/10.1061/9780784480090.002
Thube, D.T.: Artificial neural network (ANN) based pavement deterioration models for low
volume roads in India. Int. J. Pavement Res. Technol. 5(2), 115–120 (2012). ISSN 1997-1400
Xiao, F.P., Amirkhanian, S.N.: The artificial neural network approach to estimating stiffness
behavior of rubberized asphalt concrete containing reclaimed asphalt pavement.
J. Transp. Eng. 135(8), 580–589 (2009)
Yang, H.H.: Pavement Analysis and Design, 2nd edn. Pearson Education Inc./Dorling
Kindersley Publishing Inc. (2008)
Yang, S., Braham, A., Chowdhury, N., Hossain, Z.: Linking the field and lab performance of
interstate pavements. In: Innovative Materials and Design for Sustainable Transportation
Infrastructure—Selected Papers from the International Symposium on Systematic Approaches
to Environmental Sustainability in Transportation, pp. 19–27 (2015). https://doi.org/10.1061/
9780784479278
Yu, J., Xiong, C., Zhang, X., Li, W.: The durability of innovative construction materials and
structures more accurate modulus back-calculation by reducing noise information from
in situ–measured asphalt pavement deflection basin using regression model. Constr. Build.
Mater. 158(2018), 1026–1034 (2018)
Zheng, D., Quian, Z.D., Liu, Y., Liu, C.B.: Prediction and sensitivity analysis of long-term skid
resistance of epoxy asphalt mixture based on GA-BP neural network. Constr. Build. Mater.
158(2018), 614–623 (2017)
Grouting Material Development
and Treatment of Water and Mud Inflow
Caused by TBM Tunneling in Fault Zones
Jiwen Bai1,2, Shucai Li1, Rentai Liu1(&), Xiao Feng1, Peng Jiang1,
and Bingchuan Cheng1
1
Geotechnical and Structural Engineering Research Center, Shandong
University, Ji’nan 250061, Shandong, China
[email protected], [email protected]
2
Graduate School of Engineering, Nagasaki University, Nagasaki 8528521,
Japan
Abstract. Water and mud inflow caused by TBM (Tunnel Boring Machine)
tunneling in fault zones is a serious and troublesome geological obstacle in
underground engineering. The common solving method is to disassemble the
cutterhead and tunnel shield of TBM, which can result in the unstability of the
surrounding rock in the process of disassembly. Grouting without disassembly is
a practicable method to reinforce the surrounding rock and reduce water inflow.
However, the common grouting materials are difficult to remain under the
influence of fracture flow, and the grouting pressure makes the broken sur-
rounding rock and tunnel shield much easier to consolidate together. In order to
solve the above problems, a new type of grouting material GT-1 was developed.
The GT-1 grouting material possesses the characteristics of high and early
strength. The initial setting time of grouting slurry is 92–452 s, and the retention
rate of slurry in the water flow velocity (0.4 m/s) is 95%. According to the
analysis of GT-1 grouting material and the position of water and mud inflow, the
arrangement of grouting holes, grouting parameters and safety control in the
process of grouting are designed, the grouting parameters are adjusted based on
the real-time deformation of surrounding rock. The grouting effect compre-
hensive evaluation system which includes physical exploration method, core-
drilling method and excavation method is built to examine the grouting effect.
The TBM excavated across fault zones smoothly, which possesses directing
significance for similar projects.
Keywords: Water and mud inflow TBM tunneling GT-1 grouting material
Grouting design Grouting effect comprehensive evaluation system
1 Introduction
bearing zones is normal, and sometimes leads to extreme challenges (Hassanpour et al.
2009, 2010; Zhang 1999).
In unstable fractured rock or water bearing zones, shielded hard rock TBM can be
used, which erect concrete segments to support unstable tunnel walls behind the
machine (Farrokh et al. 2012). However, the geological disasters often appear with the
investigations and evaluations in the pre-construction phase. Many scholars have
conducted a lot research on geological obstacle prevention and treatment, for example,
Dammyr et al. (2017) summarized some of the extreme challenges encountered in
Norwegian subsea road tunnels. Adverse rock mass behavior and/or sudden large water
inflow at high pressure can be challenging to handle with open-face TBM, and found
out that it can be hard to predict adverse rock mass behavior ahead of the face during
tunneling. Liu et al. (2017) described the application of a comprehensive surface
geophysical investigation of underground karst systems ahead of the tunnel face.
Surface Electrical Resistivity Tomography (S-ERT), Transient Electromagnetic
Method (TEM), Geological Drilling (Geo-D) and Three-dimensional Cross-hole
Electrical Resistivity Tomography (3D cross-hole ERT) were applied to gain a com-
prehensive interpretation. To begin with, S-ERT and TEM are adopted to detect and
delineate the underground. Khave (2014) toke Zagros Water Conveyance Tunnel in
western Iran as the research object, which crossed a vast unconfined aquifer, a series of
carefully controlled field experiments in a 100 m pilot study area along the tunnel were
studied. Suitable field models are established that may be interpreted as being asso-
ciated with either water or air-filled solution channels.
At present, most of studies are about the disaster-causing mechanism and geo-
physical investigation of underground in TBM tunneling. Whereas because of com-
plexity and agnosticism of geology and the non-back characteristics of TBM, even
shielded hard rock TBM can also be stuck by water and mud inflow in fault zones
during tunneling. Considering the above occurred disasters, the rescue of TBM caused
by water and mud inflow in Yinsong tunnel is taken as the research background. A new
type of grouting material GT-1 is developed. The grouting holes, grouting parameters
and safety control are designed. The grouting effect comprehensive evaluation system
is built to examine the grouting effect.
2 Site Specification
The case presented in this paper is the Yinsong Tunnel, the main water supply line
from Songhua River in Jilin Province in the northeastern China. The joint and fault
region crossed by the line is mainly concentrated within the area of 1 km around stake
mark K46+492–K47+492. The TBM was stuck when the TBM intersected the F23-2
fault in the position of stake mark K47+371. Figure 1 is the location and geology maps
of the studied area. The F23-2 fault is mainly composed of mud, breccia, gravel and
clay fraction, which possesses tremendous liquidity under the influence of underground
water, and the water pressure ranges from 0.5 to 0.8 MPa. The total amount of water
and mud inflow is 120 m3. However, the surrounding rock around the cutterhead and
tunnel shield of TBM is still unstable. The common solving method is to disassemble
the cutterhead and tunnel shield of TBM, which has a great risk about the unstability of
Grouting Material Development and Treatment of Water … 31
surrounding rock in the process of disassembly (Barton 2000). Grouting without dis-
assembly is a practicable method to reinforce the surrounding rock and plug water.
However, the common grouting materials are difficult to remain under the influence of
fracture flow, and the grouting pressure makes the broken surrounding rock and tunnel
shield much easier to consolidate together (Barla and Pelizza 2000; Maidl et al. 2008).
(c) (d)
W E
Tunnel Face
Fig. 1. Location and geology maps of the studied area. a Location of Jilin province in China.
b General arrangement diagram of the water supply line from Songhua River, the five-pointed
star stands for the location of the studied area—Yinsong Tunnel. c The picture of water and mud
above the shield. d The relative position of tunnel face and F23-2 fault
The most commonly used grouting material is neat cement, but it is difficult to remain
under the condition of flowing underground water (Liu et al. 2011). Based on the
requirements of stuck TBM project site and grouting, a new type of grouting material
(GT-1) was developed. The GT-1 grouting material consists of cement and admixture
(90% of sodium silicate, 5% of polycarboxylic acid water-reducer and 5% of poly-
ethylene glycol), and the main performance parameters and characteristics under
dynamic water condition were tested.
were tested according to the testing standard GBT 1346-2011. For each water and
mixture (cement and GT-1) ratio (W:C), the different volume ratio of cement (Vc) and
admixture (Va) were tested separately.
As it is shown in Table 1, the compressive strength of GT-1 grouting material
decreases with the increase of water cement rate, and increases with the increase of
volume ratio of cement and admixture. In addition, the distinction of compressive
strength between 7d and 28d is not obvious, which means that the GT-1 grouting
material possesses the characteristics of high and early strength.
450 200
400 W/C 0.8:1 175 W/C 0.8:1
100
90 Vc:Va 3:1
Slurry retention rate/%
Vc:Va 5:1
80
70
60
50
0.2 0.4 0.6 0.8 1.0 1.2
flowing velocity/m/s
Grouting area
Grouting holes
10
2m
8 5m
Deformation amount/mm
-2
2015.10.4 2015.10.8 2015.10.13 2015.10.18 2015.10.23 2015.10.28 --
Date
(a) (c)
(m)
Grouting area
Before grouting
(m)
(b) 4500
4000
After grouting
3500
Wave velocity(m/s)
Consolidated
Consolidated
area in the roof
area in the face
Cutterhead
6 Conclusions
Acknowledgements. The authors would like to acknowledge the financial support from the
National Basic Research Program of China (973 Program, Grand No. 2013CB036000), the
National Natural Science Foundation of China (Grant No. 51479106).
References
Barla, G., Pelizza, S.: TBM tunnelling in difficult ground conditions. In: ISRM International
Symposium. International Society for Rock Mechanics (2000)
Barton, N.R.: TBM Tunnelling in Jointed and Faulted Rock. CRC Press (2000)
Dammyr, Ø., Nilsen, B., Gollegger, J.: Feasibility of tunnel boring through weakness zones in
deep Norwegian subsea tunnels. Tunn. Undergr. Space Technol. 69, 133–146 (2017)
Farrokh, E., Rostami, J., Laughton, C.: Study of various models for estimation of penetration rate
of hard rock TBMs. Tunn. Undergr. Space Technol. 30, 110–123 (2012)
Frenzel, C., Käsling, H., Thuro, K.: Factors influencing disc cutter wear. Geomech Tunn. 1(1),
55–60 (2008)
Hassanpour, J., Rostami, J., Khamehchiyan, M., et al.: Developing new equations for TBM
performance prediction in carbonate-argillaceous rocks: a case history of Nowsood water
conveyance tunnel. Geomech. Geoeng. Int. J. 4(4), 287–297 (2009)
Hassanpour, J., Rostami, J., Khamehchiyan, M., et al.: TBM performance analysis in pyroclastic
rocks: a case history of Karaj water conveyance tunnel. Rock Mech. Rock Eng. 43(4), 427–
445 (2010)
Khave, G.J.: Delineating subterranean water conduits using hydraulic testing and machine
performance parameters in TBM tunnel post-grouting. Int. J. Rock Mech. Min. Sci. 70, 308–
317 (2014)
Liu, R., Li, S., Zhang, Q., et al.: Experiment and application research on a new type of dynamic
water grouting material. Chin. J. Rock Mech. Eng. 30(7), 1454–1459 (2011)
Liu, B., Liu, Z.Y., Li, S.C., et al.: Comprehensive surface geophysical investigation of karst
caves ahead of the tunnel face: a case study in the Xiaoheyan section of the Water Supply
Project from Songhua River, Jilin, China. J. Appl. Geophys. 144, 36–49 (2017)
Maidl, B., Schmid, L., Ritz, W., et al.: Hardrock Tunnel Boring Machines. Wiley (2008)
Ribacchi, R., Fazio, A.L.: Influence of rock mass parameters on the performance of a TBM in a
gneissic formation (Varzo Tunnel). Rock Mech. Rock Eng. 38(2), 105–127 (2005)
Robbins, R.J.: Hard rock tunneling machines for squeezing rock conditions. In: Golser, Hinkel,
Schubert (eds.) Three Machine Concepts. Tunnels for People. 23rd Assembly of the ITA,
Vienna, pp. 633–638
Zhang, J.J.: The application and some problems of TBM and its prospects. Chin. J. Rock Mech.
Eng. 18(3), 363–367 (1999)
Zhang, J., Fu, B.: Advances in tunnel boring machine application in China. Chin. J. Rock
Mechan. Eng. 26(2), 226–238 (2007)
Traffic Data Characterization for Road
Rehabilitation: A Case Study of the Korogwe-
Mombo Road Section in Tanzania
Abstract. Traffic loading is one of the key inputs for the structural design of
pavements. For pavement design purposes, heavy vehicles are mostly used for
the estimation of the traffic load spectra, as they cause the most structural
damage to pavements. As part of the pavement design process, heavy vehicle
volume and axle load surveys are typically carried out to assist with the accurate
estimation of the cumulative traffic loading over a pavement design period.
However, traffic volumes, axle loads and, ultimately, the cumulative traffic
loading are often not uniform due to factors, such as varying motorist population
and economic activities along the length of road. In this paper, a comparative
assessment of traffic loading estimated during the rehabilitation design in 2005
and the actual measured site-specific traffic loading in 2015, as well as the
projected future traffic loading were conducted for the Korogwe-Mombo road
section along the T2 trunk road that connects Tanzania’s business hub of Dar es
Salaam with the northern regional cities. The T2 trunk road is also used by
heavy vehicles travelling to and from the neighbouring countries of Kenya and
Uganda. This study found inter alia that the cumulative traffic loading based on
the 2015 measured site-specific traffic data is approximately 2.8 times higher
than the design traffic loading based on the 2005 traffic data, which illustrates the
importance of using the latest, most accurate and reliable traffic data during
rehabilitation design. To improve the accurate determination of traffic loading,
traffic studies should ideally be conducted over a long period (typically over one
year), but this is not practical and cost effective when traditional manual
methods are used. Hence, road agencies should consider installing portable or
permanent automated traffic and Weigh-In-Motion (WIM) monitoring systems.
1 Introduction
The primary objectives of pavement structure are to provide smooth riding quality to
vehicles and protect the weaker (lower) layers such as subgrade against the effect of
traffic loading. Although pavement deterioration may be caused by other factors such
as climatic conditions, heavy vehicle volumes and axle loads are considered to be the
primary causes of pavement deterioration. As such, heavy vehicle traffic loading has
traditionally been used as the key input for the structural design of pavements (TRH 4
1996; MOW 1999). Pavement damages caused by heavy vehicles depend not only on
the Gross Vehicle Mass (GVM), but also on the mass (weight) distribution onto the
pavement. The later in turn depends on several factors, such as the number of axles on
the vehicle, axle and wheel configuration, as well as axle load, tyre inflation pressure
and contact stress (De Beer et al. 1997; Al-Qadi and Wang 2009; Greene et al. 2009;
THM 14 2014). As part of pavement design process, heavy vehicle volume and axle
load surveys are conducted to assist with the estimation of the future traffic loading
over a pavement design period—typically 20 years for most flexible pavements (TRH
16 1991; TANROADS 2003; TRL 2004). The estimated traffic loading is used to
determine the appropriate pavement design by taking into account the available
material type and climatic conditions. However, traffic volumes, axle loads, and ulti-
mately, the cumulative traffic loading are often not uniform due to varying motorist
population and economic activities along the length of the road.
In 2005, the Tanzania National Roads Agency (TANROADS) identified a need for
the rehabilitation of the Korogwe-Mombo road section (40 km long). The Korogwe-
Mombo road section forms part of the North-East corridor (T2 trunk road), which is the
main trunk road that connects Tanzania’s East Coast (including Dar es Salaam city and
Tanga Port) with the northern regional cities of Kilimanjaro and Arusha. The T2 trunk
road is also the main route that links Dar es Salaam and Nairobi, the major trade centres
of Tanzania and Kenya, respectively. The rehabilitation design of the road section was
undertaken in 2006, based on traffic counts and axle load surveys carried out in 2005
(TANROADS 2006). The rehabilitation construction was initially planned to be
completed in 2008. However, financial constraints delayed the construction works that
only started in 2012 and were completed in 2014. Follow-up traffic counts and axle
load surveys was carried out in 2015 after the road had been opened to traffic for a
period of approximately one year. Together, the 2005 and 2015 traffic counts and axle
load surveys constitute the basic data for the assessments conducted in this paper.
The objective of the current study was to conduct a comparative assessment of the
initial traffic loading estimated during the rehabilitation design phase in 2005, the actual
traffic loading determined after construction had been completed and the road had been
opened to traffic in 2015, as well as the projected future traffic loading. The assessment
provided an indication of the adequacy of the designed pavement and its expected
future performance, and also highlighted the importance of measuring/collecting site-
specific and accurate traffic volume and loading data.
Traffic Data Characterization for Road Rehabilitation … 41
Manual traffic counts are carried out by observers at a carefully selected observation
point or counting station along a road section. The traffic count survey is usually a
classified count whereby each vehicle passing an observation point is recorded on a
prepared sheet/form according to the vehicle type, and each travel direction is recorded
separately. The manual traffic count is usually undertaken over a short period (typically
seven days), as it is not practical and cost effective to undertake manual traffic counts
over a long period (i.e. 24-hour throughout the year). To improve the reliability of the
traffic data counts, it is recommended that shorter period traffic counts be undertaken
during normal days (i.e. days on which traffic patterns are not significantly affected by
public and school holidays, or other events) (TMH 14 2014; TMH 3 2015).
Automatic traffic count an alternative to the manual traffic count can measure traffic
volumes continuously over a long period and allow for capturing affects such as the
seasonal variation of traffic volumes. The commonly used automatic traffic count
42 J. J. Komba et al.
systems can be grouped into three broad types namely: pneumatic tube, magnetic wire
loops and piezo systems (TRL 2004). It is important to recognize that each of the
available automatic traffic counting technologies has certain limitations that must be
taken into account when establishing a traffic count station. Magnetic wire loop sys-
tems are the commonly used automatic traffic count systems in Southern Africa and are
generally classified into two categories namely: intrusive loops (embedded or placed on
the road pavement typically for long period traffic counts) and non-intrusive loops
(placed on pavement surface and suitable for shorter period traffic counts) (TMH 3
2015). Figure 1 show photos of typical intrusive and non-inductive traffic count
installations.
Fig. 2. Typical portable static weigh pads and WIM bending plates
Traffic Data Characterization for Road Rehabilitation … 43
Heavy vehicle axle loads are also affected by other factors such as changes to the
legal axle load limits, the level of enforcements, as well as mechanical design and load
carrying capacity of the vehicle. Table 2 shows the legal axle load limits for Tanzania
(Government Notice No. 30 2001), alongside the legal axle load limits for South Africa
(Government Notice R225 2000) and the East African Community (EAC) Vehicle
Load Control Act (2013). It can be seen that in certain axle type/groups, the maximum
permissible axle load limits for Tanzania are slightly higher than those for South Africa
and the East African Community Vehicle Load Control Act. It is further noted that the
maximum permissible load limit for the traditional dual tyres triple axle non-steering
configuration (12 tyres) is the same as that of triple axle wide-base single (super single)
tyres [6 tyres] for Tanzania (i.e., 24 tons), which may cause accelerated pavement
damage. On the other hand, the East African Community Vehicle Load Control Act has
reduced the legal load to 22.5 tons for triple axle wide-base single tyres to compensate
for the damaging effects of the wide-base tyres, while South Africa does not encourage
the use of wide-base single tyres (i.e. use 24 tons legal load for triple axle non-steering
[12 tyres], regardless of whether they are normal or wide-base tyres).
Table 2. Legal axle load limits for Tanzania, South Africa and East African Community (EAC)
Type of axle/axle group No. of tyres Maximum permissible load on
axle/axle group (tons)
Tanzania South Africa EAC
Single steering drive operated 2 8 7.7 –
Two steering drive operated 4 14 – –
Single steering draw bar controlled 4 9 – –
Single non-steering 2 8, 10, 12 8 8
Single non-steering 4 10 9 10
Tandem non-steering 4 12 16 –
Tandem non-steering 6 15 –
Tandem non-steering 8 18 18 18
Tandem steering (dolly) 8 16 – –
Triple non-steering 10 21 – –
Triple non-steering 12 24 24 24
Triple super single tyres 6 24 – 22.5
Traditionally, dual tyres have been used to limit pavement damage by efficiently
distributing the axle loads over a larger contact area than single tyres, hence reducing
the contact stresses on the pavement. Due to economic, safety and other benefits, the
wide-base single (super single) tyres are increasingly used in the trucking industry.
However, research studies demonstrate that wide-base single tyres cause more damage
to pavements than do traditional dual tyres (Al-Qadi and Wang 2009; Greene et al.
2009; Abu Abdo 2017).
44 J. J. Komba et al.
3 Study Approach
3.1 Description of the Case Study
The road section considered in this study, an approximately 40 km-long single car-
riageway with one lane in each direction, is located between Korogwe and Mombo
towns in Tanzania’s Tanga region. The Korogwe-Mombo road section forms part of
the North-East corridor (T2 trunk road), which is the main trunk road that connects
Tanzania’s East Coast (including Dar es Salaam city and Tanga Port) with the northern
regional cities of Kilimanjaro and Arusha. The T2 trunk road is also the main road that
links Dar es Salaam and Nairobi, the major trade centres of Tanzania and Kenya
respectively. In addition, the road is used by heavy vehicles travelling to and from the
neighbouring country of Uganda.
where j is traffic growth rate (%) and n is the time in years between the determination of
traffic volume and the projection year.
Traffic Data Characterization for Road Rehabilitation … 45
For most pavement design purposes, the full axle load distribution is usually not
available. As such, the concept of “E80 per heavy vehicle” (E80/HV) or Vehicle
Equivalency Factor (VEF) is used. VEF is a factor that converts different truck loads to
an equivalent number of standard axles (i.e. 8160 kg per axle). Equation 3 was used to
process the axle load survey data of the individual heavy vehicles and determine VEF.
X
i
Vehicle Equivalency Factor ðVEF) ¼ AEF ð3Þ
2
In Eq. (3), i is the total number of axles, and AEF is the Axle Equivalency Factor
that was computed using Eq. (2). The VEF values of each heavy vehicle were sub-
sequently used to determine the average VEF for each of the heavy vehicle categories.
The average VEF was calculated separately for each direction/lane.
Figure 4 compares the actual traffic volumes obtained during the 2015 traffic counts
and the volumes projected based on the 2005 traffic counts (using Eq. 1). With the
exception of Buses, the AADT projections based on the 2005 traffic count data are
generally lower than the actual AADT determined on the basis of the actual 2015
survey. For instance, while the projected AADT for Buses is approximately 12%
higher than the actual AADT in the Korogwe-Mombo direction, the projected AADT
for MGVs, HGVs and VHGVs is lower than the actual AADT by approximately 16, 29
and 34% respectively. The difference is significantly higher for the VHGVs category.
The Origin-Destination (OD) survey performed as part of the axle loads survey indi-
cated that most of the VHGV transport cement from the Tanga cement factory to the
northern regions of Kilimanjaro and Arusha, and hence they may not have been
accounted for during the 2005 traffic counts. Additionally, the general traffic growth
trends over a 10-year period could be a contributing factor.
The comparison of the actual and projected traffic suggests that the generic traffic
growth rates obtained from the Tanzanian traffic growth baseline survey report should
be used cautiously, as they may not be realistic for some roads. It further demonstrates
the need to accurately determine site-specific traffic data for pavement structural design
Traffic Data Characterization for Road Rehabilitation … 47
purposes, as opposed to using generic traffic growth rates generally derived from the
expected growth in Gross Domestic Product (GDP) of a country. It should, however, be
mentioned that both the 2015 and the 2005 traffic counts were carried out over a short
period (seven days), hence the effects of seasonal variation of traffic may not have been
accounted for.
assumed, which is the same as the design period used during the rehabilitation design
of the road section. The commutative pavement traffic loading was computed using
traffic counts and axle load survey data for both 2005 and 2015. As indicated earlier,
the traffic growth rates used were 7.0% for Buses and 6.0% for MGVs, HGVs and
VHGVs. The following three different scenarios were considered:
• Scenario 1: Use of the traffic counts and axle load survey data for 2005 to determine
the 20-year traffic loading, with 2008 as the base year. This is similar to the
approach used during the design of the rehabilitated road (i.e., the construction was
originally planned to be completed by 2008).
• Scenario 2: Use of traffic counts and axle load survey data obtained in 2005, with
2015 as a base year. This means that the recommended traffic growth rates were
applied to the actual 2005 traffic counts to project AADT for the year 2015 and then
compute the 20-year traffic loading.
• Scenario 3: Use of the actual traffic counts and axle load survey data for the 2015
survey to determine 20-year traffic loading.
Tables 3 presents the calculated 20-year cumulative pavement traffic loading
(E80s) for Scenarios 1 to 3 above as 9.7, 14.8 and 27.2 million respectively. As
theoretically expected, due partly to traffic growth over a 10-year period, the cumu-
lative traffic loading that was calculated based on the 2015 data (Scenario 3) is sig-
nificantly higher than the estimated traffic loading based on the 2005 survey data using
Scenarios 1 and 2 (i.e. approximately 2.8 times and 1.5 times higher compared with
Scenarios 1 and 2 respectively). This was expected due to the higher AADT and VEF
values obtained during the 2015 traffic counts and axle load survey respectively. Thus,
the traffic loading that was computed based on the 2005 survey data may underestimate
the actual expected future traffic loading, and it illustrates the importance of using the
latest, most accurate and reliable traffic data during rehabilitation design.
By using the Korogwe-Mombo road section in Tanzania as a case study, this paper
presented a comparative assessment of the traffic loading estimated during the reha-
bilitation design, the currently measured site-specific traffic loading, and the projected
future traffic loading. Based on the results and discussions contained in this paper, the
following conclusions are drawn and recommendations made:
• The assessments that were conducted demonstrated the need for good quality and
reliable up-to-date traffic data to ensure an accurate determination of traffic loading
for pavement design purposes.
• The common practice to conduct traffic studies over a short period of time (usually
seven days) may cause significant errors in the prediction of design traffic loading
due to the inability to capture seasonal variation of traffic; traffic pattern changes
resulting from short- to medium-term changes in economic policies; transportation
regulations; legal axle load limits; the level of enforcements; etc.
50 J. J. Komba et al.
• The generic traffic growth rates used to determine traffic loading should be applied
cautiously. Traffic volumes may vary over time due to, for instance, fluctuation in
economic trends. Furthermore, the growth rates should not be assumed to be the
same for different heavy vehicle categories, as traffic patterns may fluctuate for a
specific heavy vehicle category.
• To improve the accurate determination of traffic loading, traffic studies should be
conducted over long periods. Road agencies should also consider investing in the
installation of permanent automated traffic and WIM-monitoring systems.
• In situations where significant time delays are expected from the rehabilitation
design up to the start of construction, it is recommended that updated traffic studies
be conducted to ascertain the traffic figures, as short- to medium-term changes in
economic activities may affect traffic patterns (such as the cement factory
influencing one direction in this paper).
Overall, this study has demonstrated the importance of periodic traffic surveys to
measure and accurately quantify the changes and growth trends in traffic patterns, both
in terms of volume counts and axle loads. In lieu of manual traffic surveys, use of more
accurate automated traffic and WIM-monitoring systems is strongly recommended to
ensure accurate traffic data characterization for optimal pavement design, rehabilitation
and planning purposes. While costly permanent traffic and WIM-monitoring systems
can be stationed on selected highway sites, portable WIM technology offers a cost
effective alternative for deployment and traffic data measurements on any desired
highway location.
Acknowledgements. The authors would like to acknowledge the Tanzania National Roads
Agency (TANROADS) and the South African Council for Scientific and Industrial Research
(CSIR) for funding this study.
References
Abu Abdo, A.M.: Effect of dual versus super single truck tyre on flexible pavement performance:
a mechanistic approach. ARPN J. Eng. Appl. Sci. 12(13), 4136–4141 (2017)
Al-Qadi, I.L., Wang, H.: Evaluation of pavement damage due to new tyre designs. Research
Report ICT-09-048, Illinois Centre for Transportation, IL USA (2009)
De Beer, M., Fisher, C., Jooste, F.J.: Determination of pneumatic tyre/pavement interface contact
stresses under moving loads and some effects on pavements with thin asphalt surfacing layers.
In: Proceedings of the 8th International Conference on Asphalt Pavements, University of
Washington, USA (1997)
Federal Highway Administration (FHWA): Traffic Monitoring Guide. Washington D.C. (2013)
Government Notice No. 30.: The Road Traffic (Maximum Weight of Vehicles). United Republic
of Tanzania (2001)
Government Notice R225: The National Road Traffic Regulation: Part IV. Republic of South
Africa, Loads on Vehicles (2000)
Greene, J., Toros, U., Kim, S., Byron, T., Choubane, B.: Impact of wide-base single tyres on
pavement damage. Research Report FL/DOT/SMO/09-528, State Materials Office, Florida,
USA (2009)
Traffic Data Characterization for Road Rehabilitation … 51
MOW: Pavement and Materials Design Manual. Ministry of Works, Dar es Salaam, Tanzania
(1999)
O’Connell, J., Mgangira, M., Komba, J., Mturi, G.: Investigating the cause of deformation along
the T2 Route between Korogwe and Mkumbara. CSIR Report Number
CSIR/BE/IE/ER/2015/0067/B. Pretoria, South Africa (2016)
TANROADS: Field Testing Manual. Ministry of Works, Dar es Salaam, Tanzania (2003)
TANROADS: Detailed design for the rehabilitation of Korogwe-Mkumbara-Same Road. Design
report: volume 1. Prepared by SMEC International, Australia in association with Engineering
Research Associated (ERA), Tanzania (2006)
TANROADS: Consultancy services for carrying out baseline traffic counts in Tanzania mainland
and establishing a comprehensive traffic census methodology for TANROADS. Prepared by
Intercontinental Consultants and Technocrats Pvt. Ltd., India in association with Data Consult
Limited, Tanzania (2009)
TMH 3: Specifications for the Provision of Traffic and Weigh-in-Motion Monitoring Service.
Committee of Transport Officials, Pretoria, South Africa (2015)
TMH 14: Traffic and Axle Load Monitoring Procedures. Committee of Transport Officials,
Pretoria, South Africa (2014)
TRH 4: Structural Design of Flexible Pavements for Interurban and Rural Roads. Department of
Transport, Pretoria, South Africa (1996)
TRH 16: Traffic Loading for Pavement and Rehabilitation Design. Department of Transport,
Pretoria, South Africa (1991)
TRL: Overseas Road Note 40: A Guide to Axle Surveys and Traffic Counts for Determining
Traffic Loading on Pavements. TRL Limited, Crowthorne, Berkshire (2004)
Effect of Width of Geosynthetic Reinforcement
within the Granular Cover on the Load
Distribution over the Tunnel Lining
1 Introduction
Buried structures have been studied for centuries as tunnel linings, pipes and culverts.
Tunnels were first constructed for underground mining in Egypt, Austria and some
other countries around 3000 B.C. Tunnels for water supply started from A.D. 50, and
for road, military and buried purposes after around A.D. 600. Civil excavations became
popular around 19th and 20th centuries (Tatiya 2005). Many tunnelling projects have
been undertaken such as waterways, subways, railways and highways. In the current
accelerated society, tunnels play a very important role in the transport purpose.
For the safety aspect, the relationship between tunnels and surrounding soil is an
important factor to expect the failure in the future. That is why many researchers have
investigated the load distribution on the buried structures, such as Valsangkar and
Britto (1978, 1979), McVay and Pappadopoulos (1986), McVay et al. (1993), Li
(2016), and Talesnick et al. (2008).
One of the reinforcements for buried structures is geosynthetic reinforcement,
which is placed above the tunnel lining or buried pipe. Geosynthetics is a generic name
representing a broad range of planer products manufactured from polymeric materials.
Several types of geosynthetics are often used in the real field applications, such as
geotextiles, geogrids, geonets, geomembranes and geocomposites, which are used in
contact with soil, rock and/or any other civil engineering-related material as an integral
part of a man-made project, structure or system (Shukla and Yin 2006). In buried
structures, the effects of geosynthetic layer have been studied by some researchers,
such as Dancygier and Yankelevsky (1996), Kawabata et al. (2003), Ahmed et al.
(2015), and Shukla and Sivakugan (2013). However, these works did not deal with the
effect of width of geosynthetic layer over the tunnel lining.
In this paper, the load reductions by using geosynthetic reinforcement within the
granular cover over the tunnel lining with surface loading are presented numerically to
supplement the very limited information that exists for tunnelling applications. In
addition, the width of geosynthetic layer is varied to investigate its effect on the load
distribution over the tunnel lining, and graphical presentations are made for their
applications as design charts.
Fig. 1. Mesh generation for the tunnel covered with unreinforced soil
Effect of Width of Geosynthetic Reinforcement within … 55
The tunnel lining was modelled using the plate element function with elastic beha-
viour as assumed. The axial stiffness EA and flexural rigidity EI based on the literature
were assumed to be 1.4 107 kN/m and 1.43 105 kN m2/m, respectively with
Poisson’s ratio, v = 0.15. The crown of the tunnel lining, having diameter, Bc = 5 m,
was placed 5 m below the ground surface. Interface elements were used to simulate the
interaction between soil and tunnel lining, and soil and geogrid.
Geogrid is one of the geosynthetic reinforcements. It was placed at a depth of 4 m
below the surface, that is, 1 m above the tunnel lining. BX1100 type geogrid having a
tensile modulus of 205 kN/m was used. Its width b was varied from Bc to 6Bc. The
centreline of the geogrid was always over the centreline of the tunnel lining, as shown
in Fig. 2.
The effect of the geogrid width for crown of tunnel lining is shown in Fig. 4. It is
observed that the pressure on the crown decreases with an increase in width of geogrid.
When geogrid width is Bc, the pressure decreases slightly from 180.8 to 179.1 kPa
which is about a 1% reduction. There is a significant decrease from 180.8 to 168.5 kPa,
about 7% reduction, when the width of geogrid increases to 2Bc. The pressure
decreases gradually with 3Bc and 4Bc, which are 7.5 and 8.5%, respectively. However,
when the width of geogrid is greater than 4Bc, there is not an obvious reduction trend
for widths of 5Bc and 6Bc, reductions in pressure are 8.4 and 8.9%, respectively.
Effect of Width of Geosynthetic Reinforcement within … 57
Fig. 4. Variation of pressure on the crown of the tunnel lining with width of geogrid
reinforcement (b/Bc)
Figure 5 shows the load distribution on the springline with different widths of
geogrid. It is noticed that the pressure decreases with an increase in width of geogrid,
but the pressure reduction was not very prominent. For example, the pressure is
54.6 kPa for unreinforced case while it is to 52.7 kPa with a geogrid width of 6Bc, thus
only a reduction of about 3.6% takes place.
Fig. 5. Variation pressure on the springline of the tunnel lining with width of geogrid
reinforcement (b/Bc)
Figure 6 presents how the pressure on the invert varies with different widths of
geogrid. There is a decrease of the pressure on the invert when the width of the geogrid
increases from Bc to 5Bc; however overall, it is relatively ineffective in reducing the
58 Y. Kou et al.
Fig. 6. Variation of pressure on the invert of the tunnel lining with width of geogrid
reinforcement (b/Bc)
pressure on the invert of the tunnel lining. Compared to the unreinforced situation, even
when the width of the geogrid is 6Bc, it only decreases when the loading pressure is
increased from 198.4 to 196.6 kPa, resulting in roughly a decrease by 0.95%.
2 Conclusions
In this paper, a numerical model has been presented for analysing the effect of width of
geosynthetic reinforcement within the granular soil cover on the pressure distribution
over the tunnel lining. Based on the results and discussion presented in the previous
section, the following general conclusions can be made:
1. The width of geogrid b equal to 4Bc is the optimum width for pressure reduction on
the crown of tunnel lining. The geogrid reinforcement of this width causes a
reduction of applied pressure over the lining by 8.5%. An increase of width from
4Bc to 6Bc does not bring significant advantage in terms of pressure reduction.
2. For a maximum reduction of pressure over the springline and invert, geogrid layer
placed horizontally above the tunnel lining causes reduction in applied pressure
over the springline and invert by 3.6 and 0.95%, respectively when the width of
geogrid is 6Bc.
3. The width of geogrid layer can change the effectiveness of the load reduction
obviously on the crown of tunnel lining with the constant buried depth, however the
buried depth of the geogrid is another factor that may affect the load reduction.
Hence, future studies may take place on investigating how the buried depth of the
geogrid layer above the tunnel lining affects the pressure reduction.
Effect of Width of Geosynthetic Reinforcement within … 59
Reference
Ahmed, M.R., Tran, V.D.H., Meguid, M.A.: On the role of geogrid reinforcement in reducing
earth pressure on the buried pipes: experiment and numerical investigations. Soil Found. 55
(3), 588–599 (2015)
Bryden, P., Naggar, H.E.I., Valsangkar, A.: Soil-structure interaction of very flexible pipes:
centrifuge and numerical investigations. Int. J. Geomech. (2015). https://doi.org/10.1061/
(ASCE)GM.1943-5622.0000442
Dancygier, A.N., Yankelevsky, D.Z.: A soft layer to control soil arching above a buried structure.
Eng. Struct. 18(5), 378–386 (1996)
Kawabata, T., Uchida, K., Hirai, T., Mohri, Y., Ling, H. I., Koyama, N.: Experiment on buried
pipe using backfill of cover with geosynthetics. In: New Pipeline Technologies, Security, and
Safety Conference; Pipeline 2003, ASCE (2003)
Li, C.: Developing an analytical method to study vertical stress due to soil arching during tunnel
construction. J. Geotech. Geol. Eng. 34, 1247–1255 (2016)
McVay, M.C., Pappadopoulos, P., Bloodmquist, D.: Long-term behaviour of buried large-span
culverts in cohesive soil. Transportation Research Record 1414, Transportation Research
Board, Washington, DC, pp. 40–46 (1993)
McVay, M.C., Pappadopoulos, P.: Long-term behaviour of buried large-span culverts.
J. Geotech. Eng. 112(4), 424–442 (1986)
Shukla, S.K., Sivakugan, N.: Load coefficient for ditch conduits covered with geosynthetic-
reinforced granular backfill. Int. J. Geomech. 13(1), 76–82 (2013)
Shukla, S.K., Yin, J.-H.: Analytical Expression for Geosynthetic Strain Due to Deflection.
Geosynth. Int. 16(5), 402–407 (2006)
Talesnick, M., Horany, H., Dancygier, A.N., Karinski, Y.S.: Measuring soil pressure on a buried
model structure for the validation of quantitative frameworks. J. Geotech. Geoenviron. Eng.
(2008). https://doi.org/10.1061/(ASCE)GM.1943-5622.0000442
Tatiya, R.: Civil Excavations and Tunnelling—A Practical Guide. Thomas Telford Ltd., London
(2005)
Valsangkar, A.J., Britto, A.M.: Centrifuge tests of flexible circular pipes subjected to surface
loading. Supplemental Rep. SR530, Transportation and Rock Research Laboratory,
Crowthorne, UK (1979)
Valsangkar, A.J., Britto, A.M.: The validity of ring compression theory in the design of flexible
buried pipes. Supplemental Rep. SR440, Transportation and Rock Research Laboratory,
Crowthorne, UK (1978)
Durability Assessment of Pavement
Foundation Materials Treated
with a Polymeric-Based Additive
1 Introduction
Australia has approximately 800,000 km of road network, which is listed as the tenth
most extensive network in the world (FactBook 2016). Of this network, 60% of these
roads are unsealed. This large network is expensive to maintain and consumes a
significant amount of high quality natural materials, which is unsustainable. Attempts
have been conducted to reduce this expense and improve sustainability by using lower
quality aggregates and improving their engineering properties through different stabi-
lization techniques.
Currently pavement designers have many types of binders to use in soil stabi-
lization. The choice to use the most suitable for the soil type is crucial to obtain
enhanced engineering properties and economic performance. This has made cement the
most common binder for soil stabilization due to its compatibility with almost all soil
types and availability (Vorobieff and Wilmot 2001). However, the development of non-
traditional stabilizers, such as polymers have gained greater attention as they demon-
strate effectiveness in the field with regards to reducing permeability, increasing
Three different soil samples from three different sites in Victoria/Australia were
selected for this study, denoted as soil I, II and III. The soil samples were collected
from the top 150 mm of an existing wearing courses of several unsealed roads that
were being stabilized with PAM. For the three soil types, plasticity indices were
determined following Australian standard AS 1289.3.3.1 (AS 2009) and their classi-
fications were determined using the unified soil classification system (USCS) (ASTM
2011). The classification and geotechnical properties of the tested soils are summarized
in Table 1.
The polymeric additive used in this study was a synthetic soluble anionic PAM,
which is produced in a granulated form. This PAM is a non-toxic water-soluble
material with a specific gravity of 0.8 and a pH value of 6.9 at 25 °C, and has a high
molecular weight typically between 12 and 15 Mg/mole.
For sample preparation, the PAM amount used was 0.002% by dry weight of the
soil, according to supplier’s recommendation. The PAM was first mixed with water in a
sealed container at a rate of 2 g per 5 l which resulted in a concentration that was twice
the recommended. The concentration of this solution was adjusted when adding water
to achieve moisture requirement of each soil sample tested. The soil-water mixture was
then mixed in a mechanical mixer for 10–15 min, and the mixture was kept in plastic
bags for 24 h in a storage room at a temperature of 20 °C ± 2 to allow for an even
moisture distribution.
The maximum dry density (MDD) and optimum moisture content (OMC) were
determined using the Australian modified proctor compaction test (AS 2003). The
samples were compacted in five layers at a compaction effort of 45 blows per layer (i.e.
4961 kN m/m3) for soil type I and 35 blows per layer (i.e. 3868 kN m/m3) for soil
types II and III. This was changed from the standard 25 blows per layer (i.e.
62 R. Georgees et al.
2703 kN m/m3) as these higher compaction efforts were found to be the optimum for
each soil type and were determined after a number of trials to ensure that laboratory
performance of the PAM-treated samples was comparable with relevant observed field
performance. Thereafter, this variation in compaction energy was applied to both
treated and untreated samples. The reader is referred to a study conducted by Georgees
et al. (2015) that explains the reason for changing the compaction effort.
The repeated load triaxial (RLT) test was conducted in accordance with the rec-
ommended Australian test procedure (Austroads 2000). The specimens were prepared
in split moulds of 100 mm in diameter and 200 mm in height and were compacted in 8
layers using mechanical compaction at a compaction effort corresponding to the soil
type. Two specimens per sample (treated and untreated) for each soil type were pre-
pared. All specimens were compacted to a target density of 100% maximum dry
density at optimum moisture content. After compaction, specimens were removed
carefully from the split moulds and left to dry back to 50% of OMC. These specimens
were then subjected to the RLT test and the values of their permanent strain and
resilient moduli were determined at different stress conditions.
Permanent deformation determination characterizes the vertical permanent strain
responses under different stress conditions to examine the stress dependent permanent
strain per single test. Specimens were loaded with three different loading stages, each
involved 10,000 cycles, at stress condition of specified deviator stresses of 250, 350
and 450 kPa and a constant 50 kPa confining stress (Austroads 2000). The loading
type used is a trapezoidal pulse with a total stress period of 3 s, with 1 s loading, and
rise and fall of 0.3 s.
Durability Assessment of Pavement Foundation … 63
The capillary rise test was conducted on soil type II only and was performed
according to Australian standard AS 1141.53 (AS 1996). Soil type II was chosen
because of its high porosity compared to the other two soil types (i.e. soils I and III).
The samples were prepared by mixing the soil samples in dry conditions with the pre-
determined optimum water content. The samples were then kept in sealed plastic bags
for 24 h in a storage room at a temperature of 20 °C ± 2 to allow even moisture
distribution. The mixture was then placed in moulds (105 mm in diameter and 116 mm
in height) and compacted in five layers using the relevant optimum compactive energy
of 3868 kN m/m3. The specimens were allowed to dry back to a constant mass, after
which they were placed in water, at room temperature, in a dish to a depth of 10 mm.
Capillary rise measurements were taken at different time intervals.
Fig. 2. RLT resilient modulus values of untreated and PAM-treated samples for soil types I, II
and III at various stress conditions
66 R. Georgees et al.
To better observe the effect of PAM additive on the resilient properties, a com-
parison was conducted between the average MR values, at a certain stress stage i.e.
deviatoric and confining stresses. Here, average MR values were determined at 100 kPa
deviatoric stress and 20 kPa confining pressure which simulates a stress level closest to
the in situ stress level for a subbase layer, as recommended by Virgil Ping et al. (2001).
The results show that the impact of PAM on the MR value is more pronounced in
treated samples of soil type III with a 55.8% increase over the untreated samples. Soil
type I also exhibits a significant increase in MR with an average increase of 35%.
However, treated samples of soil type II show only a limited increase in average MR
value of 8.8%. It is believed that soils with high porosity (i.e. high water absorption,
see Table 1) need more PAM and a longer curing time to satisfy adsorption of the PAM
onto the soil surfaces to strengthen the binding between aggregate particles.
Table 2. Difference in permanent strain between treated and untreated samples of soil types I, II
and III
Soil type Permanent strain Difference in permanent strain (%)
(%) Treated versus untreated
Untreated Treated
I 0.133 0.116 −12.8
II 0.136 0.107 −21.3
III 0.180 0.162 −10.0
Table 2 shows that soil types I and II generally exhibit less permanent strain than
soil type III. Specimens of soil type I reveal less permanent strain than those of type II.
This is consistent with the fact that soils with less clay content experience lower plastic
deformations than soils with high clay content (Puppala et al. 1999). For soil type I, the
treated samples exhibit less permanent strain than those of untreated counterparts, with
the maximum difference in permanent strain being 12.8%. It is believed that the treated
samples remain sound because the internal friction is high enough due to the bonding
action of PAM, which prevent joining of particles. However, treated samples of soil
type II exhibit far less permanent strain than their untreated counterparts. The maxi-
mum difference in permanent strain is 21.3%. Similar behaviour for the treated samples
of soil type III is also observed, with the maximum decrease in permanent strain of the
treated samples at the end of the loading cycles being 10%.
Durability Assessment of Pavement Foundation … 67
Fig. 3. Abrasion resistance test results for the treated and untreated samples of soils I, II and III
Figure 3 clearly shows that there is an increase in the mass loss for all soil types
when the materials are subjected to higher brush revolutions. For the first 250 revo-
lutions of abrasion, the mass loss is greatest for soil I and least for soil III. The mass
loss is strongly dependent on the soil type (i.e. the amount of clay fines present in the
soil); the abrasion resistance is highest for soil type III followed by soil type II and soil
type I. This is supported by the role of clay in increasing the cohesion of materials. This
trend is true for both treated and untreated soils, but less intensity was identified in the
treated samples. Figure 3 also demonstrates higher abrasion resistance for the treated
samples than the untreated samples. Abrasion resistance, at the end of 250 revolutions,
is dramatically higher in the treated samples of soils I and III with an average of 45.1
and 60.9%, respectively than their untreated counterparts. Whereas soil II only showed
a 38.1% increase in abrasion resistance when treated than untreated. On the other hand,
increasing the brush load revolution to 500 did not change the trend and the mass losses
for soils I, II and III were 48, 23.5 and 59.3%, respectively.
68 R. Georgees et al.
It is believed that the adsorption of the anionic PAM onto the clay particles occurs
by cation bridging (Theng 1982), through which the polyvalent cations bridge the
negative charged groups of clay particles and polymers together. On the other hand,
bonding between particles is increased essentially by a dual process: First, the slipping
action of PAM in the compaction stage, resulting in increased density and hence
increased contact points per unit area between the soil particles; Second, the PAM
molecules encapsulate the soil particles and upon drying an electrostatic attraction
among soil particles is developed. These dual actions resulted in stronger bonds and
less mass loss for the treated materials.
Fig. 4. Percentage water rise in treated and untreated samples of soil type II
As a matter of fact, capillary energy is a function of the liquid properties along with
soil properties. In other words, if soil properties are kept constant, increased liquid
wettability (e.g. contact angle with pore aggregate surface) and increased liquid surface
Durability Assessment of Pavement Foundation … 69
tension will increase the capillary suction in the soils. However, in the case herein,
when PAM is in contact with water, the solution’s viscosity becomes higher, resulting
in increasing the time of contact between water and soil particles (Malik et al. 1991),
thereby the absorption time for the whole sample will increase. This action can be
called “external waterproofing”, rather than internal waterproofing.
4 Conclusions
Acknowledgements. The authors would like to acknowledge the logistic support of the Earthco
projects company, particularly Mr. Mark Holding for providing the stabilizer additive for this
research.
References
AASHTO: Guide for design of pavement structures. American Association of State Highway and
Transportation Officials, AASHTO, Washington, D.C., USA (1993)
AASHTO: Standard method of test for determining the resilient modulus of soils and aggregate
materials. American Association of State Highway and Transportation Officials, AASHTO,
Washington, D.C. USA (2003)
70 R. Georgees et al.
Andrews, R., Sharp, K.: A protocol for conducting field trials for best value management of
unsealed roads. In: 24th ARRB Conference—Building on 50 Years of Road and Transport
Research (2010)
AS: Methods of sampling and testing aggregates Method 53: absorption, swell and capillary rise
of compacted materials AS1141.53. Standards Australia, Australia (1996)
AS: Methods of testing soils for engineering purposes Method 5.2.1: soil compaction and density
tests—determination of the dry density/moisture content relation of a soil using modified
compactive effort. AS 1289.5.2.1. Standards Australia, Sydney, Australia (2003)
AS: Methods of testing soils for engineering purposes Method 3.3.1: soil classification tests—
calculation of the plasticity index of a soil. AS 1289.3.3.1. Standards Australia, Sydney,
Australia (2009)
ASTM: Standard practice for classification of soils for engineering purposes (unified soil
classification system). ASTM D2487. West Conshohocken, PA, USA: American Society for
Testing Materials (2011)
Austroads: Determination of permanent deformation and resilient modulus characteristics of
unbound granular materials under drained conditions. Commentary to AG:PT/T053.
Austroads, Sydney, NSW, Australia (2000)
Austroads: Guide to Pavement Technology Part 2: Pavement Structural Design. Sydney, NSW,
Australia (2008)
Camarena, S.: Sustainable Road Maintenance and Construction Utilising New Technologies.
International Public Works Conference, Darwin, Northern Territory, Australia (2013)
Factbook: Australia. Central Intelligence Agency, In: The world factbook. Available: https://
www.cia.gov/library/publications/the-world-factbook/geos/AS.html (2016)
Georgees, R.N., et al.: Effect of the use of a polymeric stabilizing additive on unconfined
compressive strength of soils. Transp. Res. Rec. J. Transp. Res. Board, 200–208 (2015)
Jones, D.: Development of performance-based tests for nontraditional road additives.
Transp. Res. Rec. J. Transp. Res. Board, 142–153 (2007)
Malik, M., et al.: Mobility of polyacrylamide and polysaccharide polymer through soil materials.
Soil Technol. 4, 255–263 (1991)
Puppala, A., et al.: Permanent deformation characterization of subgrade soils from RLT test.
J. Mater. Civ. Eng. 11, 274–282 (1999)
Sampson, L.: Development of the mechanical wet/dry brushing test for the durability testing of
stabilized materials. Transportek, CSIR, Pretoria, SA. Confidential Report C/DPVT/3.1
(1988)
Theng, B.: Clay-polymer interactions; summary and perspectives. Clays Clay Miner. 30, 1–10
(1982)
Virgil Ping, W., et al.: Measuring resilient modulus of granular materials in flexible pavements.
Transp. Res. Rec. J. Transp. Res. Board, 81–90 (2001)
Vorobieff, G., Wilmot, T.: Australian experience on subgrade stabilisation and pavement
recycling. In: International Symposium on Subgrade Stabilisation and In Situ Pavement
Recycling Using Cement, 1st, Salamanca, Spain (2001)
Wilmot, T.D.: Applications for stabilisation in pavement construction and rehabilitation. In:
Malaysian Road Conference, 1st Kuala Lumpur, Malaysia (1994)
Analytical Method to Evaluate the Stress State
Within Vertical Backfill
Abstract. In mining industry, voids created during ore removal are generally
backfilled with mill tailings that are left over from processing of mineral ore.
Evaluating the stress state in backfill is a critical step for design of backfilled
stopes. An analytical method to evaluate the stress is proposed in this paper.
Nonuniform distribution of vertical stress is considered in this method. The
trajectory of minor principle stress in the fill from the centreline to the wall is
assumed as circular-arc. The stress state within vertical backfill can be obtained
through the equilibrium consideration of differential flat element. Several
examples carried out by analytical and numerical modeling are then presented
and compared with existing analytical solutions. The proposed analytical
method is also used to evaluate an in situ project from literature. The results by
presented analytical method are in good agreement with numerical results and
in situ data from literature. So the proposed analytical method can be used to
evaluate the stress state in backfilled stopes.
1 Introduction
Widisinghe and Sivakugan (2014, 2015) proposed the vertical stress isobars for tren-
ches containing granular soils and the vertical stress isobars for silos containing food
grains and square backfilled mine stopes containing mine fill. Other numerical mod-
eling methods such as Finite Element Methods and Discrete Element Method are also
used. Goodey et al. (2006) investigated state of stress in the stored solid and the
pressures imposed on the silo walls using a Finite Element Model. Hasan et al. (2016)
investigated the distribution of intergranular forces within uncemented mine backfills
using the Discrete Element Method.
Analytical methods serve as basis for stress calculations in designs. Analytical
methods also have the advantage of simplicity. However, analytical methods are
always suffered from some limitations due to the simplifying assumptions. The
numerical modelling results (Li et al. 2003) have shown that a uniform vertical stress
across the width is not a valid hypothesis. For instance, most current existing solutions
assume that both the horizontal and vertical stress across the width is uniform (Marston
1930; Krynine 1945; Pirapakaran and Sivakugan 2007; Ting et al. 2010).
This paper presents an analytical method to evaluate the stress state within vertical
backfilled stopes. In this method, the nonuniform distribution of vertical stress across
the width of backfill is considered. The trajectory of minor principle stress in the fill
from the centreline to the wall is assumed as circular-arc. The direction of minor
principle stress at centreline is vertical to the centreline, and the direction along the rock
mass can be determined using the geometrical relationship in the Mohr’s stress circle.
The stress state of differential flat element in backfill is analyzed through equilibrium
consideration. Then the analytical solution of stress state in backfill can be obtained.
A series of examples calculated by presented analytical method and numerical mod-
eling using ABAQUS are then presented, comparison between analytical results and
in situ data is also presented and several conclusions can be drawn.
2 Analytical Method
The stress state of soils on the edge of differential flat element can be obtained
according to the geometrical relationship in the Mohr’s circle as shown in Fig. 2. The
expressions of the shear stress and the horizontal stress:
where h ¼ 45 þ u=2; sb = shear stress on the edge of differential flat element, r1b =
major principle stress on the edge of differential flat element, r3b = minor principle
stress on the edge of differential flat element, rh = horizontal stress on differential flat
element.
where rvb = vertical stress on the edge of differential flat element. Dividing Equation
(4) by Equation (5) we can get the expression of the lateral earth pressure coefficient on
the edge of differential flat element:
rh cos2 h þ Ka sin2 h
Kb ¼ ¼ 2 ð6Þ
rvb sin h þ Ka cos2 h
The friction angle of the interface between the backfill and the rock mass is equal or
less than the backfill internal friction angle. When the wall friction angle is less than the
internal friction angle, as shown in Fig. 2b. The shear stress at the interface can be
written:
The angle between minor principle stress and vertical direction at the interface can
be obtained as:
0 qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi1
1 Ka þ ð1 Ka Þ2 4 tan2 dKa
hw ¼ arctan@ A ð8Þ
2 tan dKa
An arbitrary point in differential flat element located x apart from the wall is shown
in Fig. 1. The lateral earth pressure coefficient of this point can be obtained using the
same method:
rh cos2 w þ Ka sin2 w
Kx ¼ ¼ 2 ð9Þ
rvx sin w þ Ka cos2 w
The angle u in Eq. (10) can be substituted by d, if friction between the walls and the
soil is not fully mobilized.
Average vertical stresses acting on the differential flat element can be written as:
R 2B ZB
rvx dx rvb Kb 1 ð1 Ka Þ cos2 w
rv ¼ 0
¼ dx ð11Þ
2B B ð1 Ka Þ cos2 w þ Ka
0
Bx
cos w ¼ cos hw ð12Þ
B
rv ¼ mrvb ð13Þ
Substituting Eqs. (12), (13), and (14) into Eq. (11) respectively, the shape coeffi-
cients of circular-arc are obtained, as shown:
2 pffiffiffiffiffiffiffiffi 3
1 Ka arctan cos hwpffiffiffiffi 1Ka
¼ Kb 4 15
Ka
mcir pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ð14Þ
cos hw Ka ð1 Ka Þ
z ¼ 0; rv ¼ 0 ð16Þ
Then with Eqs. (13) and (17), the expression of vertical stress of an arbitrary point
in backfill can be derived as:
Kb rv
rvx ¼ ð18Þ
Kx m
76 Q. Chen et al.
Analytical and numerical examples are presented in this section. Numerical examples
are done by Finite Element software ABAQUS. Mohr-Coulomb plasticity model is
used to model the backfill soils and rigid body is used to model the rock mass. Young’s
modulus of the backfill soils is 300 Mpa and Poisson’s ratio is 0.2. The unit weight of
backfill is 18 kN/m3. From geotechnical consideration the angle of wall friction is
slightly less than backfill internal friction angle. Pirapakaran and Sivakugan (2007)
used the value that d ¼ 2u=3 to estimate the vertical stress in a trench. In the proposed
method, the value that d ¼ 0:9u is used.
Several analytical and numerical results are shown in Fig. 3 for H = 32 m, B = 6 m,
and u ¼ 30 : The results from previous analytical solutions in literature (Marston
1930; Krynine 1945) are also shown in Fig. 3. It can be seen from Fig. 3a–d that the
results done by proposed method fit well with numerical results. For stress along the
centerline, results by proposed method, Marston’s method and numerical method have
a good agreement for both vertical and horizontal stress. Krynine’s method failed to
estimate the stress state along the centerline. For stress state along the wall, results by
proposed method, Krynine’s method and numerical method have a good agreement for
vertical stress, and results by proposed method, Marston’s method and numerical
method have a good agreement for horizontal stress. However Marston’s method failed
to calculate the vertical stress along the wall. Figure 3e and g shows a nonuniform
distribution of vertical stress which has not been considered in previous studies. The
vertical stress near the wall is much smaller than that near the centerline. From Fig. 3f
and h we can see that the horizontal stress by numerical results across the width nearly
keeps constant, which is accord with the assumption in proposed method.
Other analytical and numerical examples for different width of backfill and different
internal friction of backfilled materials are also done. However the results are not
shown in this paper because of the limitation of paper length.
Knutsson (1981) carried out in situ stress measurements within backfill stopes at
Nasliden mine of the Boliden Company for several years. The in situ data from
Knutsson is also compared with the proposed analytical method for verification pur-
poses. Knutsson’s in situ data is reproduced after Pirapakaran (2008), which contains
vertical and horizontal stresses along the depth of the backfill. The measured stope is
6 m wide, 45 m high, and 160 m long, which can be considered as plane strain con-
dition. The unit weight of the in situ backfill is 21.6 kN/m3. The internal friction angle
is 30°. The comparisons between results by proposed analytical method and in situ data
are shown in Fig. 4. A good agreement can be obtained between measured data and
results by proposed method. At the stope bottom, calculated stresses are obviously less
than overburden pressure due to arching effect, which is verified by measured data.
Analytical Method to Evaluate the Stress … 77
(a) Vertical stress along the centerline (b) Horizontal stress along the centerline
(c) Vertical stress along the wall (d) Horizontal stress along the wall
(e) Vertical stress at different depth (f) Horizontal stress at different depth
(g) Vertical stress by numerical method (h) Horizontal stress by numerical method
4 Conclusions
An analytical method to evaluate the stress state in vertical backfill considering the
nonuniform distribution of vertical stress is proposed in this paper. The trajectory of
minor principle stress in backfill from the centreline to the wall is assumed as circular-
arc. Then the analytical solution is obtained through the equilibrium consideration of
differential flat element. Examples done by proposed analytical method and numerical
method are presented. The comparison between proposed analytical method and in situ
data from literature is also presented. Several conclusions can be drawn:
1. Results by proposed analytical method have a close match with the numerical
results by numerical method. Previous analytical method failed to match with
numerical results obtained from ABAQUS because previous approaches cannot
present the nonuniform stress state of backfill soils.
2. The stress distribution across the backfill width is nonuniform. The trajectory of the
minor principle stress from the centreline to the wall can be represented as circular-
arc. The horizontal stress component keeps same along width, while the vertical
stress component at centreline is higher than that along the wall.
3. Comparison between proposed analytical method and in situ data from literature
also shows a good agreement. So the proposed analytical method can be used to
evaluate the stress state in backfilled stopes.
Acknowledgments. The authors acknowledge the financial supports from the National Key
Basic Research Program of China (973 Program) (No. 2015CB057801), National Natural Sci-
ence Foundation of China (NSFC Grant No. 51278449, NSFC Grant No. 51238009 and NSFC
Grant No. 51338009).
Analytical Method to Evaluate the Stress … 79
References
ABAQUS User’s Manual, Version 6.1, Hibbitt, Karlsson and Sorenson Inc. (2000)
Goodey, R.J., Brown, C.J., Rotter, J.M.: Predicted patterns of filling pressures in thin-walled
square silos. Eng. Struct. 28(1), 109–119 (2006)
Hasan, A., Karrech, A., Chareyre, B.: Evaluating force distributions within virtual uncemented
mine backfill using discrete element method. Int. J. Geomech. 17(7), 06016042 (2016)
Janssen, H.A.: Versuche über getreidedruck in silozellen. Zeitschr. d. Vereines deutscher
Ingenieure 39(35), 1045–1049 (1895)
Knutsson, S.: Stresses in the hydraulic backfill from analytical calculations and in-situ
measurements. In: Conference on the Application of Rock Mechanics to Cut and Fill Mining:
01/06/1980-03/06/1980, pp. 261–268. The Institution of Mining and Metallurgy (1981)
Krynine, D.P.: Discussion on ‘Stability and stiffness of cellular cofferdams’ by Karl Terzaghi.
Trans. Am. Soc. Civ. Eng. 110(1), 1175–1178 (1945)
Li, L., Aubertin, M., Simon, R., Bussière, B., Belem, T.: Modeling arching effects in narrow
backfilled stopes with FLAC. In: Proceedings of the 3rd International Symposium on FLAC
& FLAC 3D Numerical Modelling in Geomechanics, Ontario, Canada, pp. 211–219 (2003)
Marston, A.: The theory of external loads on closed conduits in the light of the latest experiments.
In: Highway Research Board Proceedings, vol. 9. (1930)
Paik, K.H., Salgado, R.: Estimation of active earth pressure against rigid retaining walls
considering arching effects. Géotechnique 53(7), 643–653 (2003)
Pirapakaran, K., Sivakugan, N.: Arching within hydraulic fill stopes. Geotech. Geol. Eng. 25(1),
25–35 (2007)
Pirapakaran, K.: Load-deformation characteristics of minefills with particular reference to arching
and stress developments. Doctoral dissertation, James Cook University (2008)
Ting, C.H., Sivakugan, N., Read, W., Shukla, S.K.: Analytical method to determine vertical
stresses within a granular material contained in right vertical prisms. Int. J. Geomech. 12(1),
74–79 (2010)
Widisinghe, S., Sivakugan, N.: Vertical stress isobars for trenches and mine stopes containing
granular backfills. Int. J. Geomech. 14(2), 313–318 (2014)
Widisinghe, S., Sivakugan, N.: Vertical stress isobars for silos and square backfilled mine stopes.
Int. J. Geomech. 16(2), 06015003 (2015)
Numerical Simulation of Surface Subsidence
After the Collapse of a Mine
1 Introduction
Any underground work can cause surface subsidence, which could damage infras-
tructure and buildings. Longwall coal mining is not an exception. In order to mitigate
consequences of the subsidence and choose the appropriate method of mining, it is
important to predict the size and depth of the trough. Some empirical methods of
surface subsidence prediction have been developed in different countries based on the
observed local data; for example, one of them is in the Subsidence Engineers’
Handbook (NCB 1975) developed in the UK. FLAC3D is commercial software, which
has been used to predict the surface subsidence trough (Herrero et al. 2012 and Xu et al.
2013). However, due to the complex behaviour of the rock, there is minimal confidence
in predictions from numerical modelling and more research is required (Xu et al. 2013).
The purpose of this work is to increase the understanding and develop the procedure of
the numerical simulation of the surface subsidence with limited information on the
properties of the overburden.
The following procedure are followed:
– Deriving model parameters from the borehole data;
– Setting the domain of the model;
– Assigning the constitutive model;
– Assigning the parameters derived earlier;
– Setting initial and boundary conditions;
– Running the model to equilibrium;
– Assigning the special constitutive model to the excavating zone to model goaf
behavior;
– Assigning the parameters, which are inherited from the roof layer or derived from
the literature, to the goaf;
– Assigning characteristics of the volumetric behavior of the goaf;
– Running the model to equilibrium;
– Altering volumetric characteristics in the goaf to obtain the required goaf height
after the simulation (adjusting the goaf height).
The paper explains the procedure and discusses the modelling results in three parts.
The first part explains a method developed by the authors for estimating and assigning
the appropriate physical-mechanical properties to the model using a visual description
of the borehole log. This method was created after an extensive literature search and
based on the works of Balmer (1952), Deere (1968), Hoek and Brown (1980), Hansen
(1988), Palmström (1995), Hoek and Brown (1997), Palmström and Singh (2001), and
Hoek and Diederichs (2006). The second part describes the site of interest, the model
domain, the mesh density, the initial and boundary conditions. The site is at Naburn in
North Yorkshire, England. A measured subsidence profile was obtained from UKCoal.
The mine collapsed controllably after Longwall mining without stowing. The third part
discusses results, the modelled stresses in the goaf, the debris caused by a mine col-
lapse, the influence of the Geological Strength Index (GSI) and surface subsidence.
2 Physical-Mechanical Properties
1 D=2
E ¼ Ei 0:02 þ ð1Þ
1 þ eð60 þ 25DGSIÞ=11
where GSI is the Geological Strength Index introduced by Hoek and Brown (1997) and
describes the insitu state of the rock. A GSI of 100 is for a very good, undisturbed rock
mass whereas a GSI of 0 is for a very poor quality, disintegrated rock mass. D is the
disturbance factor, which is dependent upon the excavation conditions, i.e. blasting a
rock face will give the rock a disturbance factor of 1 while careful excavation will yield
a disturbance factor of 0. Ei is given by:
Ei ¼ MR rci ð2Þ
where MR is the modulus ratio classified by rock type and presented in Table 1, MR
was first proposed by Deere (1968) and later modified by Palmström and Singh (2001).
rci is the uniaxial compressive strength of the intact rock, which could be found in
Table 2.
Table 1. Guidelines for the selection of the modulus ratio. Based on Deere (1968) and
Palmström and Singh (2001)
Rock type Class Group Texture
Coarse Medium Fine Very fine
Sedimentary Clastic Conglome-rates Sandstones 200– Siltstones 350–400 Claystones
300–400 350 Greywackes 350 200–300
Breccias 230–350 Shales 150–
250a
Marls 150–
200
Non- Carbonates Crystalline Sparitic limestones Micritic limestones Dolomites
clastic limestones 400–600 600–800 800-1000 350–500
Evaporites Gypsum (350)b Anhydrite (350)b
Organic Chalk 1000+
Metamorphic Non-foliated Marble 700–1000 Hornfels 400–700 Quartzites 300–450
Metasandstone
200–300
Slightly foliated Migmatite 350–400 Amphibolites Gneiss 300–750a
400–500
Foliateda Schists 250–1100a Phyllites/mica schist Slates 400–
300–800a 600a
Igneous Plutonic Light Granitec 300–550 Dioritec 300–350
Granodioritec 400–450
Dark Gabbro 400–500 Dolerite 300–400
Norite 350–400
Hypabyssal Porphyrics (400)b Diabase 300–350 Peridotite
250–300
Volcanic Lava Rhyolite 300–500 Dacite 350–450
Andesite 300–500 Basalt 250–450
Pyroclastic Agglomerate 400– Volcanic breccia Tuff 200–400
600 (500)b
a
Highly anisotropic rocks: the value of MR will be significantly different if normal strain and/or loading occurs parallel (high
MR) or perpendicular (low MR) to a weakness plane. Uniaxial test loading direction should be equivalent to field application
b
No data available, estimated on the basis of geological logic
c
Felsic Granitoids: coarse grained or altered (high MR), fined grained (low MR)
Numerical Simulation of Surface Subsidence After the Collapse … 83
Table 2. A summary of uniaxial compressive strength based on rock type (Palmström 1995)
Rock name Uniaxial compressive Rating of the Rock name Uniaxial compressive Rating of the
strength rc, MPa factor m1)
i strength rc, MPa factor m1)
i
Low Average High Low Average High
Sedimentary rocks Metamorphic rocks
Anhydrite 120b? 13.2 Amphibolite 75 125 250 31.2
Coal 16c 21c 26c Amphibolitic 95 160 230 31?
gneiss
Claystone 2b 5b 10b 3.4 Augen gneiss 95 160 230 30?
Conglomerate 70 85 100 (20) Black shale 35 70 105
Coral chalk 3 10 18 7.2 Garnet mica 75 105 130
schist
Dolomite 60b 100b 300b 10.1 Granite gneiss 80 120 155 30?
Limestone 50a 100b 180a 8.4 Granulite 80b 150 280
Mudstone 45 95 145 Gneiss 80 130 185 29.2
Shale 36c 95c 172c Gneiss granite 65 105 140 30?
Sandstone 75 120 160 18.8 Greenschist 65 75 85
Siltstone 10b 80b 180b 9.6 Greenstone 120b 170a 280a 20?
Tuff 3b 25b 150b Greywacke 100 120 145
Igneous rocks
Andesite 75b 140b 300b 18.9 Marble 60b 130b 230b 9.3
Anorthosite 40 125 210 Mica gneiss 55 80 100 30?
Basalt 100 165 355c (17) Mica quartzite 45 85 125 25?
Diabase 227c 280c 319c 15.2 Mica schist 20 80a 170a 15?
(dolerite)
Diorite 100 140 190 27? Mylonite 65 90 120
Gabbro 190 240 285 25.8 Phyllite 21 50 80 13?
Granite 95 160 230 32.7 Quartz 70 120 175
sandstone
Granodiorite 75 105 135 20? Quartzite 75 145 245 23.7
Monzonite 85 145 230 30? Quartzitic 45 100 155
phyllite
Nepheline 125 165 200 Serpentinite 65 135 200
syenite
Norite 290c 298c 326c 21.7 Slate 120b 190b 300b 11.4
Pegmatite 39 50 62
Rhyolite 85b? (20) Talc schist 45 65 90 10?
Syenite 75 150 230 30?
Ultra basic 80b 160 360
rock
Soil materials2)
Very soft clay rc = 0.025 MPa Soft clay rc = 0.025 − 0.05 MPa Firm clay rc = 0.05 − 0.1 MPa
Stiff clay rc = 0.1 − 0.25 MPa Very stiff clay rc = 0.25 − 0.5 MPa Hard clay rc => 0.5 MPa
Silt, sand: assume rc = 0.0001 − 0.001 MPa
a
Values found by the Technical University of Norway, (NTH) Inst. for rock mechanics
b
Values given in Lama and Vutukuri, 1978
c
Values given by Bieniawski, 1984
1)
mi refers to the value of m for intact rock in the Hoek-Brown model. Values in brackets have been estimated by Hoek
et al. (1992) while those with a question mark have been assumed by Palmstrom
2)
For clays, the values of the UCS is based on ISRM (1988)
84 Y. G. Derbin et al.
where, mi refers to the value of m for intact rock in the Hoek-Brown model and is
summarized in Table 3, and
GSI 100
s ¼ exp ð5Þ
9
Then, the values for r3 must be selected. Hoek and Brown (1997) concluded that
the most consistent results are obtained when 8 equally spaced values between
0 < r3 < 0.25rci are used.
To find the tangent of the failure surface at the appropriate stress level, first the non-
linear analytical solution for Mohr’s envelope is found and then a linear regression
analysis is used to find the equation of the tangent at that point.
Balmer’s analytical solution (Balmer 1952) to Mohr’s envelope describes the
relationship between the normal and shear stresses in terms of the principal stresses as:
0 0
0 r r
0
rn ¼ r3 þ 0 1 0 3 ð6Þ
@r1 =@r3 þ 1
0 0
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
0 0
s ¼ r 1 r3 @r1 =@r3 ð7Þ
Y ¼ logðAÞ þ B X ð10Þ
The constants A and B can then be calculated using a linear regression analysis, i.e.
P P
Y =T B X=
A ¼ 10 T ð13Þ
P P P
X Y ð X Y Þ=T
B¼ P P ð14Þ
ð X 2 ð X 2 Þ=T Þ
86 Y. G. Derbin et al.
where T is the number of values in the sequence, i.e. 8, if the earlier suggestion is
followed.
Finally, the Mohr-Coulomb parameters can be deduced from the following two
equations:
0 B1 !
1 rni rtm
/ ¼ Tan AB ð15Þ
rci
0
c ¼ 0:75 s rni Tanð/Þ ð16Þ
Since A and B are known, then using rni in place of rn, (i.e. general notation of the
normal stress), in Eq. 12 will evaluate an expression for X which can be used in Eq. 10
to find a value of Y. This can be used to calculate s from Eq. 11.
It can be seen that another new parameter (rni) has been introduced and this is the
value of the normal stress at the point of interest. To determine this value we need to
turn to the work of Hoek and Brown (1980).
Hoek and Brown found the correlation (Eq. 17) between depth and vertical insitu
stress based on the collated worldwide data from researchers investigating the insitu
state of stress underground (Fig. 1).
Fig. 1. Insitu vertical stress data. after Hoek and Brown (1980) in Hudson and Harrison (2000)
The vertical stress calculated by Eq. 17 is given in MPa. This value is equivalent to
the normal stress and enables us to completely specify the mechanical properties of the
rocks underground.
Numerical Simulation of Surface Subsidence After the Collapse … 87
3 Model Development
3.1 Site of Interest
A measured subsidence profile and borehole log description were taken from both
above and under the Barnsley seam at Naburn in North Yorkshire, UK. The overburden
consists of siltstone, sandstone, mudstone, and seatearth (claystone underlying coal
seam). The bulk and shear moduli were calculated using Eqs. 1 and 2 for the elastic
stiffness and assuming a Poisson’s ratio of 0.2 for all layers. The internal friction and
cohesion were derived from Eqs. 15 and 16. The tensile strength was calculated by
Eq. 9. A density of 2700 kg/m3 was taken as an average value for these types of rock
after data collected by Shtumpf (1994) in Table 4.
E0 e
r¼ ð19Þ
ð1 e=em Þ
where E0 is the initial tangent modulus and em is the maximum strain of the goaf
material.
Since the parameters are difficult to estimate, and even sometimes impossible, the
authors of this paper go further with assumptions. Equation 18 can be rewritten con-
sidering the coefficient c = 1/em, which is used to adjust the height of the goaf after a
simulation, and E is the Young’s modulus of the roof material.
Ee
r¼ ð20Þ
ð1 ceÞ
Numerical Simulation of Surface Subsidence After the Collapse … 89
The correctness of the modelled goaf behaviour is possible to check by two facts:
the goaf height and the stresses in the goaf after the simulation. The stresses in the goaf
will be discussed in the next section.
The required goaf height at the end of the simulation can be estimated by the initial
height of the seam. The required height can be found by the multiplication between the
height of the mine and the subsidence factor: these are collected in Table 5 for different
regions of the world by Bräuner (in Bell and Donnelly 2006). The subsidence factor
depends on the region and whether or not the excavated area has been filled or packed.
The factor is used to calculate maximum possible subsidence. Table 5 shows the
subsidence factor varies from 0.33 to 0.9 for the different regions. The recommended
factor is 0.9 for mines in UK. For the case at hand with a height of the excavation of
2.8 m, the required final height of the goaf is estimated as 0.28 m. Altering the
parameter c, the simulation was repeated until the goaf height after the simulation
becomes within 5% error of the required height. For the sake of simplicity, this process
of obtaining the required goaf height by altering the volumetric characteristics of the
goaf is called ‘adjusting goaf height’.
In order to model the strain-hardening behaviour of the goaf material, the double-
yield model, which allows both shear and volumetric compression, is traditionally
implemented. In FLAC, the stress-strain curve is approximated by a table to generate a
linear piecewise curve. In the developed model, the table has 10 rows. The elastic
properties, bulk and shear moduli, and Mohr-Coulomb properties, friction and cohe-
sion, correspond to the properties of the roof material.
As it will be shown later, the double-yield model cannot simulate the goaf beha-
viour precisely enough. Instead, one of the Critical State models, namely the modified
Cam-clay model, was implemented. Derbin et al. (2016) showed that the modified
Cam-clay model predicted the goaf behaviour more accurately than the double-yield
model did. Computational application of Critical State theory includes different types
90 Y. G. Derbin et al.
of soil and soft rock (Gens, and Potts 1988). Xiao et al. (2016) successfully imple-
mented the Critical State concept to predict the behaviour of coarse granular soil
(which is a material suitable for a rockfill dam) in a true triaxial compression test. The
rockfill could best be described as a very coarse granular type of soil according to the
British Soil Classification System (BS 5930: 1981). The very coarse soils are cobbles
with sizes of 63–200 mm, boulders with sizes of 200–630 mm, and large boulders with
sizes of more the 630 mm. Singh and Singh (2011) argued that goaf consists of 22.5%
boulders and 77.5% large boulders. It should be also noted that the elastic properties
depend on strain in the modified Cam-clay model. This corresponds better to the real
behaviour of the goaf (Badr et al. 2003).
The goaf material was described by critical state parameters, i.e. lambda
(k) = 0.188, kappa (j) = 0.007 and a frictional constant (M) = 1.9, which were used
by Indraratna and Salim (2002) to model drain triaxial shearing on crushed basalt. By
changing either the specific volume at reference pressure on the normal consolidation
line or the pre-consolidation pressure, the required height of the goaf can be obtained.
Figure 3 shows curves of dependence between the obtained goaf height and specific
volume for three different pre-consolidation pressures, i.e. 1e5, 1e4, and 1e3 Pa.
4 Model Results
models agree closely with each other. This agreement means that the method of the
calculation of the Mohr-Coulomb properties from the Hoek-Brown parameters is
correct.
where cp is the plastic shear strain, %; c0 is the cohesion at the peak strength of the rock
where cp = 0, and n is the fitting parameter, which depends on rock type and its
magnitude varies from 0.29 for Sandstone to 0.34 for Mudstone (Pourhosseini and
Shabanimashcool 2014). After the evaluation of the effect of this parameter on the
subsidence profile, no effect was noticed after the goaf height was adjusted. A mean
value of 0.3 was taken in this work.
Figure 4a demonstrates that the strain-softening constitutive model predicts an
identical subsidence profile to those profiles obtained by the Mohr-Coulomb and Hoek-
Brown models. However, the earlier investigations (i.e. by Lloyd et al. 1997) showed
that the strain-softening model was capable of predicting a deeper trough than the
Mohr-Coulomb model. This can be explained by examining the zone of plasticity.
Figure 5a shows that the area of plastic deformation occurs directly above and under
the excavation. After adjusting the goaf height, the strain-softening model shows the
identical results as the Mohr-Coulomb model.
92 Y. G. Derbin et al.
Fig. 5. Plastic deformation resulting from the excavation, a GSI is constant with the depth,
b GSI gradually increased with the depth
and siltstone (which is the roof material). The layers are too thin to be clearly shown on
the diagram (Fig. 6). Figure 6a shows the constant GSI and the increasing GSI with the
depth. Figure 6b presents the diagram of the stiffness, which is calculated according to
Eq. 1 vs depth for two cases: when the GSI is constant and when the GSI changes from
25 to 85 with the depth according to Fig. 6a. In Fig. 5b, it can be seen that the model
with the increasing GSI experiences more plastic defamation above the seam than the
model with the constant GSI (Fig. 5a).
The plastic zone distribution influences the performance of the strain-softening
model. In contrast to the model with the constant GSI (Fig. 4a), Fig. 4b shows that the
strain-softening model with the increasing GSI predicts a deeper trough than the Mohr-
Coulomb model; however, this trough predicts and erroneous width if it is compared to
the field measurements. More research on computer modelling of the surface subsi-
dence is needed.
Fig. 7. Compression vertical stress within the goaf, a fictitious model (Derbin et al. 2016) and
b theoretical (after Wilson 1983)
94 Y. G. Derbin et al.
In order to investigate stress recovery in the goaf, Derbin et al. (2016) developed a
simplified fictitious subsidence model where the goaf length was equal to distance D1.
Distance D1 was calculated using the following equation suggested by Mukherjee et al.
(1994) (after Wilson 1984):
5 Conclusions
The paper describes how to simulate surface subsidence after a goaf collapse for any
seam in any part of the world using only the description of the borehole log. The
method of calculation of the Mohr-Coulomb properties out of the Hoek-Brown
parameters was introduced. The predictions of the Mohr-Coulomb and modified Hoek-
Brown constitutive models embedded in FLAC agree with each other. This proves the
correctness of the method used to calculate the properties. The strain-softening method
predicts identical deformation, magnitude, and spatial extent for a constant GSI, but the
variable GSI causes a deeper trough.
A pattern of plastic deformation around the goaf is crucial for the strain-softening
model. When plastic deformation above and under the goaf is in the same quantity, the
strain-softening effect of the model is compensated by adjusting the goaf height. All
these plastic zones act as one big goaf. The pattern of plastic deformation depends on
relationship between the GSI and depth. The GSI increasing with the depth produces a
larger failure area above the goaf, which deepens the trough of the strain-softening
model. Unfortunately, this is not sufficient to match field observations in this research,
and the trough keeps the same erroneous width. Hence, practitioners should use tra-
ditional constitutive models to predict surface subsidence with great care.
For the time being, the best solution would be utilization of the strain-softening
constitutive model keeping in mind that it can underestimate the depth and overesti-
mate the width of the subsidence trough. To improve the prediction of the sizes of the
trough, more advanced constitutive models should be implemented. The recent
research by Derbin et al. (2018) has shown that the bubble model has better results. It
predicts a deeper and narrower trough than the strain-softening model does.
The paper also discusses how to improve the modelling of the goaf behaviour.
A Critical State model, i.e. the modified Cam-clay model, was implemented into the
goaf. It was shown that if the goaf length is not sufficient to recover the primary
stresses, the modified Cam-clay model and the double-yield model predict identical
goaf behaviour. Hence, it is difficult to say how much of the model discrepancies can
be attributed to the goaf and how much can be attributed to the 700 m of overburden.
This work suggests deeper investigation of modelling the goaf behaviour by increasing
the goaf length up to distance D1 and implementing more advanced constitutive
models into the goaf.
Acknowledgements. The authors greatly acknowledge the open fund SKLGDUEK1512 of the
China University of Mining and Technology.
References
Itasca.: Manual FLAC3D Version 5.01 (2013)
Badr, S., Ozbay, U., Kieffer, S., Salamon, M.: Three-dimensional strain softening modeling of
deep longwall coal mine layouts. FLAC Numer. Model. Geomech. 233–239 (2003). https://
doi.org/10.1201/9780203883204.ch79
96 Y. G. Derbin et al.
Balmer, G.A.: General analytic solution for Mohr envelope. In: Proceedings-American Society
for Testing and Materials, pp. 1260–1271 (1952)
Bell, F.G., Donnelly, L.J.: Mining and its Impact on the Environment, p. 547. Taylor & Francis
(2006). https://doi.org/10.2113/gseegeosci.13.3.270
Bieniawski, Z.T.: Rock mechanics design in mining and tunneling. Balkema. (1984)
Bräuner, G.: Subsidence due to underground mining: theory and practices in predicting surface
deformation, US Bureau of Mines, p. 8571 (1973). https://doi.org/10.1016/0148-9062(74)
91643-x
British Standard.: Code of practice for site investigation, p. 147 (1981)
Deere, D.: Geological considerations. Rock mechanics in engineering practice, pp. 1–20 (1968)
Derbin, Y., Walker, J., Wanatowski, D., Marshall, A.: Issues related to goaf modeling. In:
Proceedings of the 19th Southeast Asian Geotechnical Conference & 2nd AGSSEA
Conference, pp. 1015–1021 (2016)
Derbin, Y., Walker, J., Wanatowski., Marshall, A.: Implementation of advanced constitutive
models for the prediction of surface subsidence after underground mineral extraction.
(Expected to be published) (2018)
Gens, A., Potts, D.: Critical state models in computational geomechanics. Eng. Comput. 5, 178–
197 (1988). https://doi.org/10.1016/0148-9062(89)90209-x
Hansen, T.H.: Rock properties. Norwegian Rock and Soil Assoc. 5, 41–44 (1988)
Herrero, C., Munoz, A., Catalina, J.C., Hadj-Hassen, F., Kuchenbecker, R., Spreckels, V., Juzwa,
J., Bennet, S., Purvis, M., Bigby, D., Moore, D.: Prediction and monitoring of subsidence
hazards above coal mines (Presidence), p. 138. Directorate-General for Research and
Innovation, Luxembourg (2012)
Hoek, E., Brown, E.T.: Underground excavations in rock, p. 527. Institution of Mining and
Metallurgy, London (1980)
Hoek, E., Brown, E.T.: Practical estimates of rock mass strength. Int. J. Rock Mech. Min. Sci. 34,
1165–1186 (1997). https://doi.org/10.1016/s0148-9062(97)00305-7
Hoek, E., Diederichs, M.S.: Empirical estimation of rock mass modulus. Int. J. Rock Mech. Min.
Sci. 43, 203–215 (2006). https://doi.org/10.1016/j.ijrmms.2005.06.005
Hoek, E., Wood, D., Shah, S.: A modifed Hoek-Brown failure criterion for jointed rock masses.
In: Rock Characterization: ISRM Symposium, Eurock’92, Chester, UK, 14–17 September
1992, pp. 209–214 (1992)
Hudson, J.A., Harrison, J.P.: Engineering rock mechanics: an introduction to the principles,
p. 444, Elsevier (2000). https://doi.org/10.1115/1.1451165
Indraratna, B., Salim, W.: Modelling of particle breakage of coarse aggregates incorporating
strength and dilatancy. Proc. Inst. Civil Eng. London 155(4), 243–252 (2002). https://doi.org/
10.1680/geng.155.4.243.38691
ISRM: Commission on standardisation laboratory and field results. Suggested methods for
determining hardness and abrasiveness of rocks. Int. J. Rock Mech. Min. Geomech. Abstr. 15,
89–97 (1988)
Lama, R., Vutukuri, V.: Handbook on mechanical properties of rocks-testing techniques and
results. 3 (1978)
Lloyd, P., Mohammad, N., Reddish, D.: Surface subsidence prediction techniques for UK
coalfields-an innovative numerical modelling approach. In: Proceeding of the 15th Mining
Congress of Turkey, Ankara, pp. 6–9 (1997)
Mukherjee, C., Sheorey, P.R., Sharma, K.G.: Numerical simulation of caved goaf behaviour in
longwall workings. Int. J. Rock Mech. Min. Sci. Geomech. 31, 35–45 (1994). https://doi.org/
10.1016/0148-9062(94)91292-0
Numerical Simulation of Surface Subsidence After the Collapse … 97
1 Introduction
investigated the bearing behavior of piled rafts in overconsolidated clay and concluded
that the interaction between piles and raft is a major influence. Sinha and Hanna [9]
have carried a 3-D numerical analysis on piled raft foundation also using finite element
method.
The numerical methods can be categorized into two main categories: continuum
and discontinuum approaches [10]. After many years’ development, finite element
method has become very mature and widely applied in geotechnical engineering field
[11, 12]. Compared with the discrete element method, its computational efficiency is a
great advantage. However, finite element method which based on continuum
mechanics could not describe the microstructure accurately of the soil and provide the
micro mechanics insights into the load transfer mechanism [13, 14]. Meanwhile, soil is
discontinuous on itself, when the finite element method is implemented, the hypothesis
brought into the simulation will affect the accuracy of the result. Cundall and Strack
first developed the discrete element method (DEM) [15, 16], which is currently an
effective simulation method to obtain inter-particle forces and trace microstructure
evolution of granular materials. Many researchers have conducted DEM simulations on
the soil behavior around pile foundations. erface undergoing loading and stress relief is
studied [17], and soil-pile interaction during pile penetration is also investigated [18].
Liu and Zhou [19] simulated the pile-soil interaction subject to uplift load. Lobo-
Guerrero and Vallejo [20] investigated the influences of pile shape and pile interaction
on the crushable behavior of granular materials around driven piles. A discrete element
modeling of a granular platform supported by piles in soft soil is performed to validate
a small scale model test [21]. Lai et al. [22] studied the “soil-arching” within geogrid-
reinforced and unreinforced pile-supported embankments with DEM. Ng and Meyers
[23] investigated the side resistance of drilled shafts in granular soils by DEM. Because
the computation cost is very high by DEM, the combined DEM and other continuous
numerical methods has been developed and applied to geomechanics problems, such as
coupled finite element and discrete element method [24, 25], coupled discrete element–
finite difference method [26] etc. However these coupled methods are more compli-
cated when performing the analysis.
At present, it is difficult to build the model at the real scale for engineering
problems with DEM method. Since the sand particle size is very small, to closely
simulate the real conditions, millions of particles are needed to build a pile-soil model
for running a DEM simulation, which is difficult to realize now based on personal
computers. In the current study, in order to obtain micro mechanics behavior and also
save computation time, a radial particle expansion method called hierarchical modeling
is proposed to study the pile-soil-cap interaction behavior for non-displacement piles.
Fig. 2. Stress-strain curves for soils of different zones under confining pressure of 100, 150, and
200 kPa
102 H. Zhao et al.
Fig. 3. Determination of the internal friction angle by the use of Mohr circles
A New Method Based on PFC3D and Hierarchical Modeling … 103
Fig. 3. (continued)
104 H. Zhao et al.
Figure 2 displayed the triaxial test results by DEM simulations for the soils of different
zones. Figure 3 plotted the Mohr coulomb strength criteria for the soils in zones 1–4
and provided the strength data. Since the middle zone is the interested area, we only
displayed the grading curve of the middle zone soils used in the triaxial test in Fig. 4
and it can be seen as gravel soil according to soil classification standard. Soils in these
four different zones have the similar friction angle. The obtained average macro
strength parameters is u = 35.5°, c = 0. This also verified that the proposed radial
particle expansion method is able to give a similar macro behavior for soils in different
zones. Table 1 summarized the parameters for the soils in different zones.
Fig. 4. Grading curve of the middle zone soil used in the triaxial test
model’s width to the pile cross section length is 22.72 greater than 20, which minimized
the boundary effect. To make the study simple, total 8 piles were built as listed in Table 3.
Model 0MM is a single pile without cap. Model LMM, MMM, SMM have different pile
cap size but with the same pile length and pile section area. Model MMM, MML, MMS
have the same pile cap size and pile section but different pile length. MBM, MSM, MMM
have different pile diameter with the same pile cap size and pile length. Pile model MMM
was used in different group studies. All the pile models are shown in Fig. 5.
Fig. 5. Model series of single pile with and without pile cap
106 H. Zhao et al.
In order to trace the change of porosity and stress of sand particles during loading, a
series of measurement spheres are installed in the model as shown in Fig. 6. Because of
the symmetry of the model, measurement spheres are installed only on the right half of
the model. More measurement spheres are installed below the pile bottom and around
the pile cap. Since the stress is the average stress over the sphere, to make the measured
stress more uniform and reliable, overlapping of measurement spheres is allowed. The
measured quantity is the average over the volume of measurement sphere. Total 132
measurement spheres are installed.
Fig. 7. The variation of cap resistance, tip resistance, shaft resistance, and total load for piles
with different pile caps
The pile cap resistance increases with the pile cap size. Pile cap resistance tends to
increase first nonlinearly to a maximum and then drops gradually to a residual state.
This can be explained that when the load applied exceeds the bearing capacity of the
soil beneath the pile cap, the soil fails finally. With the larger pile cap size, the
reduction of pile cap resistance tends to be more significant.
For pile model 0MM, the shaft friction first increases to the maximum value then
drops to a certain constant value. This indicates that the static friction between pile and
soil contributes to the shaft friction at the beginning; once the pile is moving, the static
friction turns into smaller sliding friction. Pile model SMM also has a similar variation
but with a slightly larger shaft resistance value than 0MM. Pile model MMM and
LMM have an increasing shaft resistance at the initial settlement and then drop sharply.
After the drop, the shaft resistance increases back again to a certain value. Pile model
LMM has a larger shaft resistance when the settlement is larger than 0.18 m while
MMM has a smaller value of shaft resistance.
The pile tip resistance for all the pile models fully develops at the early settlement
stage and then drops a little bit, after that gradually increases. The soil surrounding the
pile tip started to move with loading which reduces the bearing resistance of soil and
caused the first drop of pile tip resistance. Fully developed shaft friction happens earlier
than the pile tip resistance. Model MMM has the lowest pile tip resistance, while the
other three models have a similar vale of pile tip resistance.
108 H. Zhao et al.
The total load capacity for the pile models increases with the pile cap size. The
larger capped piles tend to have a larger linear portion. The slope of settlement-load
curves is more flat for the pile model with larger pile cap size. The remaining load
capacity is also increasing with the pile cap size.
(a) s=0.05m
(b) s=0.14m
(c) s=0.28m
at different settlement s = 0.05, 0.137, and 0.288 m. Pile model LMM has a very
complex force chain net. But obviously the strong force chains were developed below
the pile cap on both sides of the pile body. These force chains transfers the load mostly
downward with some portion away from the pile. There are weak chains along the
lower part of the pile above the pile bottom. Below the pile bottom, there is some force
chains also formed. These force chains are not as strong as those below the cap but
stronger than those weak force chains along the lower part of the pile. For the pile with
small pile cap or without pile cap, the force chains are weak and the major force chains
only located below the pile bottom. Arching effect was observed for LMM and MMM
because of the strong force chains formed on the top part of the model. These char-
acteristics of force chains might be able to help us to design better pile foundation
systems.
3.4 Porosity
From the DEM simulation, through the installed measurement spheres, porosity dis-
tribution can be obtained. Porosity is an important mesoscopic parameter to describe
the granular materials. Porosity distribution on the central plane of the model was
plotted in Fig. 10 at three different settlement values of 0.05, 0.137, and 0.228 m.
Three distinct porosity zones are observed from the plots. A zone with larger soil
porosity is next to the pile because of the excavation, follows a second narrow smaller
porosity zone, then a third zone with a porosity value among the above two zones. As
observed from these plots, there is a relatively large gradient change for the porosity
close to the pile. Further away from the pile, there is almost no change in porosity. With
the pile moving down, the soil is compacted below the pile cap and the porosity is
reducing, while the soil around the lower pile body is increasing. The porosity of soil
below the pile bottom is reducing due to settlement of the pile tip. This phenomenon
may be due to the dragging effects of the pile. The pile is moving down, while the
dilatancy effects made the sand particles around the pile loose. These results also
revealed that the porosity change is limited to a narrow zone close to the drilled piers.
110 H. Zhao et al.
(a) Stress σz at s= 0.05m for pile LMM, MMM, SMM, and 0MM
(b) Stress σz at s= 0.14m for pile LMM, MMM, SMM, and 0MM
(c) Stress σz at s= 0.28m for pile LMM, MMM, SMM, and 0MM
(a) s=0.05m
(b) s=0.14m
(b) s=0.28m
Fig. 10. Soil porosity distribution map with different pile cap size
112 H. Zhao et al.
4 Conclusions
In the application of discrete element method, the enormous computational cost has
always been a problem. This paper proposed a radial expansion particle generation
modeling method which can greatly reduce the computational cost and verified its
reliability. It will extend a new approach for discrete element modeling and be applied
to many cases like retaining wall and underground cavern which are in low confining
pressure conditions.
By the proposed radial expansion particle generation method, this study investi-
gated the pile cap effects. A series of DEM simulations were carried out. From the
results, the following conclusions can be obtained:
1. A radial particle expansion method is reliable to generate the model, particles in
different radial zones verified to have the same macro stress-strain curves and
strength data.
2. When the pile cap is large enough, both the pile tip resistance and pile shaft
resistance are less than the pile cap resistance, larger pile cap and larger pile
resistance. Pile cap resistance follows a variation of first increasing to a maximum
value and then reducing after the failure of soil beneath the pile cap. Shaft resistance
seems to increase with the enlarging pile cap area. All pile tip resistance increases
nonlinearly to a maximum value and then softens due to soil failure below the pile
tip. The full mobilization of pile tip resistance is later than the pile cap resistance
which is reasonable.
3. Force chains for capped pile have two strong networks on the top and at the bottom
of the capped pile foundations, while no capped pile only has one well developed
force chains network at the bottom.
4. In general, three distinct porosity zones are distributed within the soil mass. Due to
dilatancy effects, soil porosity is larger next to the pile body. During loading, soil
porosity beneath the pile cap and pile bottom is reduced to form a compaction zone.
Localized change in porosity happens in the soils close to the pile.
Overall, the interaction of pile-soil-cap is a systematic problem. Pile cap plays a
very important role in improving the loading capacity of the pile foundation. The load
transfer mechanism for the capped pile is more rational than no capped pile. This paper
does not discuss the structure failure of the pile and a rigid pile and a rigid cap
considered. Utilizing the potential resistance of pile cap may save a lot of construction
cost.
Acknowledgements. This research project was financially supported by the National Natural
Science Foundation of China (No. 11672066). The authors are greatly appreciated for this
support which made this study possible.
A New Method Based on PFC3D and Hierarchical Modeling … 113
References
1. Hongladaromp, T., Chen, N.J., Lee, S.L.: Load distribution in rectangular footings on piles.
Geotech. Eng. 4(2), 77–90 (1973)
2. Kuwabara, F.: An elastic analysis for piled raft foundations in a homogeneous soil. Soils
Found. 29(1), 82–92 (1989)
3. Clancy, P., Randolph, M.F.: An approximate analysis procedure for piled raft foundations.
Int. J. Numer. Anal. Meth. Geomech. 17, 849–869 (1993)
4. Poulos, H.G.: An approximate numerical analysis of pile–raft interaction. Int. J. Numer.
Anal. Meth. Geomech. 18, 73–92 (1994)
5. Ta, L.D., Small, J.C.: Analysis of piled raft systems in layered soils. Int. J. Num. Anal.
Methods Geomech. 20, 57–72 (1996)
6. Ta, L.D., Small, J.C.: An approximation for analysis of raft and piled raft foundations.
Comput. Geotech. 20(2), 105–123 (1997)
7. Zhang, H.H., Small, J.C.: Analysis of capped pile groups subjected to horizontal and vertical
loads. Comput. Geotech. 26, 1–21 (2000)
8. Reul, O.: Numerical study of the bearing behavior of piled rafts. Int. J. Geomech. 2(4), 59–
68 (2004). https://doi.org/10.1061/(asce)1532-3641
9. Sinha, A., Hanna, A.M.: 3D numerical mode for piled raft foundation. Int. J. Geomech. 17
(2), 04016055 (2017)
10. Duan, K., Kwok, C.Y., Tham, L.G.: Micromechanical analysis of the failure process of
brittle rock [J]. Int. J. Numer. Anal. Meth. Geomech. 39(6), 618–634 (2015)
11. Gourvenc, S.M., Powrie, W.: Three-dimensional finite element analysis of diaphragm wall
installation. Geotechnique 49(6), 801–823 (1999)
12. Yasushi, A., Osamu, K., Osamu, M., et al.: A numerical study of ground displacement and
stress during and after the installation of deep circular diaphragm walls and soil excavation.
Comput. Geotech. 35(5), 791–807 (2008)
13. Cundall, P.: A discontinuous future for numerical modelling in geomechanics? Proc. ICE
Geotech. Eng. 149(1), 41–47 (2001)
14. Potyondy, D.O., Cundall, P.A.: A bonded-particle model for rock. Int. J. Rock Mech. Min.
Sci. 41(8), 1329–1364 (2004)
15. Cundall, P.A., Strack, O.D.L.: A discrete numerical model for granular assemblies.
Geotechnique 29, 47–65 (1979)
16. Potyondy, D.O., Cundall, P.A.: A bonded-particle model for rock. Int. J. Rock Mech.
Mining Sci. 41(8), 1329–1364 (2004)
17. Peng, S.Y., Ng, C.W.W., Zheng, G.: The dilatant behaviour of sand–pile interface subjected
to loading and stress relief. Acta Geotech. 9, 425–437 (2014). https://doi.org/10.1007/
s11440-013-0216-9
18. Zhou, J., Jian, Q., Zhang, J., et al.: Coupled 3D discrete-continuum numerical modeling of
pile penetration in sand. J. Zhejiang Univ. Sci. A. 13(1), 44–55 (2012). https://doi.org/10.
1631/jzus.a1100172
19. Liu, W.-B., Zhou, J.: Numerical simulation of particle flow code for pile under uplifting
load. Chin. J. Geotech. Eng. 26(4), 516–521 (2004)
20. Lobo-Guerrero, S., Vallejo, L.E.: Influence of pile shape and pile interaction on the
crushable behavior of granular materials around driven piles: DEM analyses. Granular
Matter 9(3), 241–250 (2007)
21. Jenck, O., Dias, D., Kastner, R.: Discrete element modelling of a granular platform
supported by piles in soft soil—validation on a small scale model test and comparison to a
numerical analysis in a continuum. Comput. Geotech. 36, 917–927 (2009)
114 H. Zhao et al.
22. Lai, H.-J., Zheng, J.-J., Zhang, J., Zhang, R.-J., Cui, L.: DEM analysis of “soil”-arching
within geogrid-reinforced and unreinforced pile-supported embankments. Comput. Geotech.
61, 13–23 (2014)
23. Ng, T.-T., Meyers, R.: Side resistance of drilled shafts in granular soils investigated by
DEM. Comput. Geotech. 68, 161–168 (2015)
24. Villard, P., Chevalier, B., Le Hello, B., Combe, G.: Coupling between finite and discrete
element methods for the modelling of earth structures reinforced by geosynthetic. Comput.
Geotech. 36(5), 709–717 (2009)
25. Hamdi, P., Stead, D., Elmo, D.: Damage characterization during laboratory strength testing:
A 3D-finite-discrete element approach. Comput. Geotech. 60, 33–46 (2014)
26. Indraratna, B., Ngo, N.T., Rujikiatkamjorn, C., Sloan, S.W.: Coupled discrete element–finite
difference method for analysing the load-deformation behaviour of a single stone column in
soft soil. Comput. Geotech. 63, 267–278 (2015)
Investigation and Numerical Simulation
Analyses of the Landslides in Terrace
Formation
Abstract. In situ exploration and sampling are carried out in Sin-San Village in
Nantou County after the typhoon. Large direct shear tests in the laboratory are
performed to measure the shear strength of the material. The soil water char-
acteristic curve of the soil was derived from the pressure plate tests. The cross-
section profiles of the failed slope before and after the typhoon at this area are
obtained using digital elevation model (DEM) and 3D terrestrial scanner,
respectively. The numerical softwares, SEEP/W and SLOPE/W, are used to
analyze the slope stability under the similar rainfall by the typhoon. Different
ground water level, river level, and friction angle of the soil are considered to
compare the values of factor of safety and to reveal the failure processes of the
slope. The possible failure of the slope may be due to scouring of the toe of the
slope by the surge of rising river level and/or mud flow based on the investigated
profiles and the numerical simulations. The estimated length of scouring was
about 20 m long. The possible mechanism of failure for this length of scouring
could be due to multiple retrogressive failures after toe scouring based on the
numerical modeling.
1 Introduction
In 2009, the moderate Typhoon Morakot invaded Taiwan brought more than 2500 mm
of torrential rain in central and southern Taiwan. It was measured as 3059.5 mm in
Alishan rainfall station which has been closed to the average rainfall (3684 mm) per
year in Alishan. This astonishing rainfall has exceeded the hundred year rainstorm
frequency. This was the deadliest typhoon that ever struck Taiwan caused serious
disasters in southern and central Taiwan, leaving 461 people dead and 192 others
missing. The storm also produced huge amount of rainfall, peaking at 1623, 2361, and
2748 mm in 24, 48, and 72 h, respectively. Some mountain villages in southern Tai-
wan were covered by the soils from the landslide and debris flow without any warning.
Thus, it is very necessary to study and understand the causes of failure and failure
mechanisms under extreme weather conditions on the young deposits in these moun-
tain areas.
The infiltration of heavy rainfall made the factor of safety of slopes to decrease and
caused landslides initiated from the upper and/or downward slopes. The landslides may
then cause and trigger mudslides or debris flow. The mud flow or debris flow may
squeeze the main river channel and scour river banks and/or down slopes of the roads.
Therefore, a lot of landslides had occurred during this typhoon along the road of Tai-21
in central Taiwan. The focus of this research is to study the causes and failure
mechanism of slope failures of high terrace deposits in central Taiwan under hydro-
geology and extreme weather conditions.
Slope failures due to rainfall conditions have been studied by many researchers and
the impacts of unsaturated and saturated soil behavior on the slope stability are also
considered in many papers. The objective of this paper is to study and understand the
failure processes and mechanisms of the slope failure of high terrace deposits at Shin-
Shan Village in Nantou County during Typhoon Morakot. There are thirty-six tribu-
taries that classified as potential debris flow torrents inside this watershed. Thus, a lot of
failures and disasters had occurred at the roads and slopes within the Chen-Yu-Lan
River watershed. In situ site investigations and sampling have been performed at these
areas. Numerical modeling using softwares, SEEP/W and SLOPE/W, developed by
Geo-Studio are chosen to simulate the behavior of the slope under both extreme
weather and erosion conditions.
2 Disaster Area
The location of the study area is in Nantou County, which situates at the center of
Taiwan, and mountain areas cover most of the areas in this county, as shown in Fig. 1.
Thus, a lot of disasters, including flooding, slope failures, debris flow, and road and
bridge failures, have occurred in this county, especially inside Chen-Yu-Lan Water-
shed, during Typhoon Morakot. Investigation had been accomplished after the typhoon
to document the location and extend of various types of failures. Thus, subsequent
researches and analyses could be performed to reveal the causes of disasters. Several
major types of disasters, i.e., landslide, debris flow, slope failure due to scouring, have
been found during the in situ investigation in Nantou County. Most of the disasters
were location beside or close to the Chen-Yu-Lan River.
There are 39 slope related failures occurred in this typhoon, and 24 failures are
from upslope and the others are from downslope. This paper will focus on the slope
failures at the Shin-Shan Village which is located just beside Chen-Yu-Lan River. The
aerial photo taken before the typhoon and location of the village and watershed are
show in Fig. 1. The area within the red dashed line was eroded away during the
typhoon. There are two potential debris flow torrents, DF169 and DF170, inside this
village, and they all flow into Chen-Yu-Lan River. The extreme weather condition
induces the quick rise of river level and may cause the flow to turn and erode banks
beside the river. The downslope of this area collapsed and fallen down to the river and
caused part of the base of road and houses on this terrace to be empty, shown in Fig. 2.
Twenty three houses above the slope fell down to the slope into the river because of the
slope failure, as shown in Fig. 2. Both sides of the river banks have been eroded away
and the river channel was widen during the typhoon. Kriging interpolation method is
Investigation and Numerical Simulation Analyses … 117
Legend
Study Area
Fig. 1. Locations of study area and Chen-Yu-Lan River Watershed (Aerial Photo Taken in
2005)
used to obtain the hourly rainfall during the typhoon based on thirteen rainfall stations
inside Chen-Yu-Lan River Watershed because there is no rainfall station near the
village.
Fig. 2. Downslope of Shin-Shan Village collapsed and near fallen houses beside the river
118 S.-C. Hsu et al.
Chen-Yu-Lan River is at the central Taiwan and originates from the north peak of Yu
Mountain with an elevation of 3910 m. The river length is 42.4 km with an average
declination slope of 5% and its watershed area is about 45,000 ha. There are three
major faults, Dili, Chenyulan, and Shanshihchia Faults across this watershed which
causes the rock formation contains many fractures and discontinuities. The
Fig. 3. Geological formations of Chen-Yu-Lan River watershed (After Hsu et al. 2011)
Investigation and Numerical Simulation Analyses … 119
Chenyulanchi fault, which almost parallels Chen-Yu-Lan River, separates the Miocene
sedimentary rocks of the Western Foothills from the east side of the Oligocene
metamorphic rocks of the Shuehshan Range (Fig. 3). Differential uplifting along this
fault has generated great topographic relief and abundant fractures that resulted in
frequent landslides and debris flows even prior to the Chi-Chi earthquake (Lin et al.
2002). Slates and meta-sandstones are the dominant lithologies in the metamorphic
terrane. Based on the relative amount of slate and meta-sandstone, the metamorphic
strata are divided into four formations: in ascending order, the Shihpachuangchi, the
Tachien Meta-Sandstone, the Paileng Meta-Sandstone, and the Shuichangliu. In the
sedimentary terrane, sandstone and shale predominate and include the Nankang,
Nanchuang, and Kueichulin Formations (Lin et al. 2002). The slope failure at Sin-San
Village is located in the recent deposit and can be classified as terrace deposit. This
terrace deposit could be formed by previous flood, debris flow, and/or sliding materials.
Fig. 4. Volumetric water content versus matric suction and fitted VG model
Before
After
Elevation, m
Distance, m
The maximum grain size of this terrace deposit is 30 cm based on the in situ sampling.
Thus, the equal-weighted method is used to reduce the maximum grain size of the
tested samples to 10 cm and adjust the grain-size distribution accordingly in order to
perform large direct shear test. Modified and standard compaction tests are made to
obtain the compaction energy similar to the in situ soil condition. Then, a large direct
shear tests (50 50 cm) are used to acquire the shear strength of the terrace deposit
under saturated and unsaturated conditions. The shear strength of the sample with 5%
water content is larger than the saturated sample. According to the in situ large direct
shear test, 1.5 1.5 m (Chu et al. 1989) performed in central Taichung, the peak and
residual friction angles were 54.3° and 44.9°, and apparent cohesions were 14.7 and
0 kN/m2, respectively. The tested results of samples with 5% water content are used to
compare and combine with the in situ large direct shear test and plotted in Fig. 6. The
peak and residual friction angles are about 54° and 46°, respectively. The peak apparent
cohesion is about 25 kPa. Considering the shear strength may decrease as water content
increases during rainfall. Therefore, the friction angle and apparent cohesion of the
terrace deposit used for later numerical modeling are chosen as 46° and 15°,
respectively.
Fig. 6. Results of direct shear tests from this study and Chiu et al. (1989)
122 S.-C. Hsu et al.
6 Numerical Modeling
Rainfall Intensity
Fig. 8. The schematic profiles and boundary conditions before and after rainfall infiltration
The slope is modeled using SEEP/W under the rainfall intensity in Fig. 7, and then uses
SLOPE/W to calculate the factor of safety at each time. The factor of safety for the
slope before rainfall is 2.61. The lowest factor of safety for the slope during rainfall
without considering change of river level is 1.23 after 60 h of raining, i.e. the slope is
still stable. According to the observation from the residents in that village, the river
level goes up very quick and started scouring the toe of the slope during the typhoon.
Thus, the river level is also assumed to go up as the groundwater level rose due to
infiltration. The tail groundwater level rises from 330 to 343 m after 60 h of rainfall,
and the head groundwater level rises from 335 to 350 m. The lowest factor of safety for
the slope considering rising of river water level is down to 1.187, as shown in Fig. 9,
after 60 h of rain. The lowest factor of safety for the slope will be down to 0.915 if the
analysis considering rise of river water level and scouring effect at the toe pf the slope.
Thus, retrogressive failure will be initiated since the factor of safety is less than 1. At
least three subsequent sliding failures could occur after the first slide. Figure 10 shows
slope profile and the factor of safety after the toe erosion and the first slide. The factor
of safety is 0.937. The slope is still unstable and will slide again. The obtained slope
profile after four sliding failure is similar to the actual failed profile after the typhoon.
Therefore, the major cause of the slope failure could be due to scouring and erosion of
the river on the toe of the slope.
124 S.-C. Hsu et al.
Rise of River
Water Level
Fig. 10. Profile of lowest factor of safety after toe erosion and first slide
8 Conclusions
Failure of high terrace deposit during Typhoon Morakot beside Chen-Yu-Lan River,
Nantou County, Taiwan is chosen for this study. The following conclusions can be
drawn based on the study and numerical modeling:
1. Part of the village is eroded away based on the aerial photos taken before and after
the typhoon.
2. Soil water characteristic curve is obtained from the pressure plate tests and fitted
with VG model.
3. The slope could remain stable if the river water level had not gone up and the toe of
the slope was not eroded away during Typhoon Morakot.
4. The factor of safety will decrease if the river water level goes up.
Investigation and Numerical Simulation Analyses … 125
5. The slope will become unstable if the toe of the slope is eroded away due to the
river flow.
6. Retrogressive failure will occur if the river water level goes up and the toe of the
slope is eroded away.
Acknowledgements. The authors would like to thank National Science Foundation in Taiwan
for financially supporting this research. Also thanks to Nantou Branch of Directorate General of
Highways of Taiwan for providing the valuable information and photos.
References
Chu, B.-L., Hwang, C.L., Cheng, S.Y., Pan, G.M.: Study of in situ direct shear tests in gravel
formation in Taiwan. In: The Third Geotechnical Conference, Taiwan, pp. 695–706 (1989)
Geo-Slope: User’s Manual for SEEP/W and SLOPE/W: Version 4. Geo-Slope International,
Canada (2013)
Hsu, S.-C., Lai, B.J., Hsu, W., Lai, J.: Modeling of gravel properties and anchors using discrete
element method. Adv. Mater. Res. 189–193, 1726–1731 (2011)
Lin, P.-S., Lin, J.-Y., Hung, J.-C., Yang, M.-D.: Assessing debris-flow hazard in a watershed in
Taiwan. Eng. Geol. 66, 295–313 (2002)
Wu, C.-C., Wang, C.-H., Chan, Y.-T.: Infiltration behavior of gravel-rich soil. J. Chin. Soil Water
Conserv. 43(1), 65–74 (2012)
Performance of Three Atmospheric Density
Models on Precise Orbit Determination
for Haiyang-2A Satellite Using DORIS Data
Abstract. DORIS has become a matured space geodetic technique after more
than ten years of development. This geodetic technique is mainly applied to
determine the orbit for the low Earth orbit (LEO) satellite. There are a number of
non-conservative forces acting on Haiyang-2A (HY-2A) satellite with altitude of
about 970 km, in which the atmospheric drag is the most dominant and
uncertainty in the precise orbit determination (POD) with the dynamic method.
In order to achieve POD for HY-2A, MSIS-86, Jacchia 1971 and DTM87
models were evaluated in this study. The precise orbits of HY-2A from DORIS
data were compared with the precise orbit ephemeris computed by the Centre
National d’Etudes Spatiales (CNES). Tests demonstrated that the relative opti-
mal atmospheric density model was the empirical MSIS-86 model for HY-2A
satellite with the corresponding drag coefficient of 2.0. The RMS of orbit dif-
ference between the derived orbits using MSIS-86 and the CNES orbits was
0.842 cm in the radial direction, and 3.899 cm in three dimensions (3D). This
study will provide valuable references for the LEO satellites with similar altitude
and surface, especially for the other HY-2 series satellite of China.
1 Introduction
Haiyang-2A Satellite (HY-2A) was launched on August 16, 2011. It is the first satellite
designed to monitor the global marine dynamic environment with microwave sensors
and detect ocean surface wind field, height and temperature in China. It is also the first
mission of China’s Earth observation satellite HY-2 series (including HY-2A, HY-2B,
HY-2C, and HY-2D) and China’s first radar altimetry satellite. Precision orbit deter-
mination (POD) is an essential objective for HY-2A to provide fundamental references
for the analysis and application of the on-board equipment. In order to achieve the
essential goal, HY-2A satellite carries three kinds of instrument for POD activities, that
is, a Doppler Orbit Determination and Radio positioning Integrated by Satellite
(DORIS) receiver, a dual frequency Global Positioning System (GPS) receiver and a
Satellite Laser Ranging (SLR) retro-reflector array (Guo et al. 2013a; Kong et al. 2014).
DORIS is an excellent geodetic technique supporting precise orbit determination
(POD) of the low Earth orbit satellite (LEO) and geodesy (Zelensky et al. 2010). The
current DGXX generation DORIS receiver onboard HY-2A and Jason-2 can simulta-
neously track 7 DORIS beacons with 2 7 channels (Mercier et al. 2010), and range-
rate and phase measurements on both channels with RINEX format can be provided at
the same epoch. The evenly distributed beacons and the capability of tracking 7 DORIS
beacons can improve the geometry strength and increase the observation quantity
significantly from the former generation DORIS receiver.
There are a number of non-conservative forces acting on a LEO satellite, in which
the atmospheric drag is the most dominant and also contains the most uncertainty. The
atmospheric drag is directly proportional to atmospheric density. There are many
factors affecting the variability of atmospheric density, and some of these factors have
not been well modeled, such as atmospheric heating and the solar and geomagnetic
activity levels. The Earth’s upper atmospheric density changes with the solar activity,
time, season, latitude and geomagnetic index, and the mechanisms are too complex to
be accurately modeled. Atmospheric drag is one dominant constrain on the accurate of
orbit determination. It is very important to make a full research on the atmospheric
density model and atmospheric drag coefficient to achieve POD for the LEO satellite
using DORIS data. At present, significant researches have taken place to determine the
proper modeling for atmospheric density. Fattig (2011) demonstrated that the RMS
differs very little when the Jacchia family of atmospheric density models, namely CIRA
1972, Jacchia (1971), and Jacchia-Robert (Robert 1971) were used to determine the
precise orbits of GRACE-A and GRACE-B, and the test showed that the Jacchia 1971
was the optimal model in the Jacchia family. Apart from the well-known and frequently
applied Jacchia model, a variety of other density models of the upper atmosphere exits,
such as MSIS (Mass Spectrometer and Incoherent Scatter) and DTM (French Density
Temperature Model) families. Li et al. (2010) found that the impact of MSIS-86 (Hedin
1987) atmosphere density model on the accuracy of the orbit varied largely with
different altitude, and the orbital accuracy in the radial direction is about 1.5 cm, and
the one in along track direction is about 3.0 cm using DORIS technique. Mance (2010)
precisely determined orbits of GEOSAT Follow-On (GFO), Starlette, Stella, and Geo-
Forschungs-Zentrum-1 (GFZ-1) satellites using MSIS-86 atmospheric density model.
Lechtenberg (2010) had tested the performance of different atmospheric density models
(such as Jacchia, MSIS and DTM families) using GPS technique for the POD of
CHAMP and using SLR technique for the POD of GRACE, and found that the Jacchia
based models outperformed the other two density model families. Lechtenberg (2015)
estimated atmospheric densities using SLR data and made a full research on the density
correction factors and standard deviations for the ANDE Castor satellite.
In order to realize precise orbit determination for the HY-2A satellite, several
researches have presented their strategies using different geodetic techniques and dif-
ferent atmospheric density models. Zhao et al. (2013) and Wang et al. (2014) have used
the DTM94 (Berger et al. 1998) atmospheric model and have achieved the radial
accuracy of better than 3 and 4 cm using SLR geodetic technique, respectively. Guo
et al. (2013b) and Lin et al. (2014) have used the same model and have got 1–2 cm and
128 Q. Kong et al.
3 cm radial accuracy using GPS technique, respectively. Kong et al. (2014) have
applied DTM87 model (Barlie et al. 1978) and about 1.5 cm accuracy has been gained
in radial direction using DORIS technique. Gao et al. (2015) have used MSIS-86
atmospheric density model and accuracy about 1.1 cm has been obtained in radial
direction with DORIS data. Both DTM and MSIS atmospheric density family models
were applied based on SLR, GPS and DORIS data for HY-2A, respectively.
Many studies and even some international standards have promoted one model over
another, while it is necessary to select the best model for a particular mission. The goal
of this research is to invest the performance of MSIS-86, Jacchia 1971 and DTM87
density models, and thus to find a relatively accurate coefficient and density model to
enhance atmospheric drag calculations for the POD of HY-2A satellite.
The observed atmospheric drag forces act in the opposite direction of the satellite
velocity. Accurately modeling atmospheric forces is difficult because of the difficulty of
availability of the physical properties, complexity of the interaction of neutral gas with
different spacecraft surfaces, and variation of attitude of spacecraft. Usually satellite
surface is divided into a number of limited basis units to calculate the atmospheric drag
on per unit area. The widely accepted equation for atmospheric drag is
1 A
D€rd ¼ CD qv2r ^vr ð1Þ
2 ms
where the negative sign indicates that the acceleration is anti-parallel to the unit relative
velocity vector (^vr ), which provides the direction of the acceleration. CD is introduced
to model the actual momentum transfer dependent upon the interaction between the
satellite and the atmosphere. Typically, the nominal value of CD ranges from 2.0 to 2.5
(Montenbruck and Gill 2000). A is the cross sectional area of the orbiting body and is
normal to the velocity vector. ms is the mass of the orbiting body, and vr is the velocity
of the satellite.
Considering the amount of tracked DORIS data and strength of geometric structure of
beacons, we should apply the dynamic method to determine the orbit for HY-2A
satellite. Six prior independent parameters at the initial epoch are needed for the orbit
integration for the dynamic method. These six parameters can be the three-dimensional
position and velocity vectors or six orbital elements. HY-2A is disturbed by the con-
servative forces and the non-conservative forces as this mission travels around the
Earth. The conservative forces include the Earth’s gravitation, n-body perturbation,
Performance of Three Atmospheric Density Models on Precise … 129
oceanic and solid earth tides, et al., and the non-conservative forces include the solar
and earth radiations, atmospheric drag and relativistic effect, and other small perturbing
factors. The multi-step COWELL II numerical integration was used during the process
of orbit determination for HY-2A (Balmino et al. 1990), and the step size for orbit
integration is 10 s and the output interval of orbit is set to 60 s. During the orbit
determination, the empirical force model was applied to compensate the other orbit
dynamic perturbations which cannot be modeled and estimated (Kang et al. 2006; Liu
2013). In the process of parameter estimation using the DORIS Doppler tracking data,
all parameters are solved together, such as the coordinate, velocity, solar radiation
coefficients and empirical acceleration coefficients. The initial state vectors and
empirical acceleration coefficients in the cross track are estimated, and the coefficients
of solar radiation are estimated every 6 h. The least squares estimation method is
applied to solve the unknown parameters. The cutoff elevation angle is set to 10°. Due
to TAI time system used in DORIS system, the constant offset (leap second) of 35 s is
introduced to be in consistent with UTC time. Table 1 summarizes the applied dynamic
models.
In this paper, we applied DORIS range rate data to compute the precise orbit of
HY-2A. During the 3-day arc from 8 Sep., 2012 to 10 Sep., 2012, HY-2A was tracked
by 48 DORIS stations. Figure 1 displays the geographical location of DORIS beacon
stations and the 3-day trajectory of HY-2A.
130 Q. Kong et al.
Fig. 1. Distribution of DORIS tracking stations on the global and the trajectory of HY-2A from
September 8 to 10, 2012
The coefficient of drag Cd is ultimately difficult to define for the complex shape of HY-2A
satellite. The determination of the Cd is related to several parameters such as the
accommodation coefficient, momentum accommodation coefficient, the temperature of
the atmosphere, the temperature of the satellites surface, altitude, the velocity of the
satellite, solar minimum and solar maximum, as well as other parameters (Mehta et al.
2014; Pilinski et al. 2011). Drag coefficients for satellites in the upper atmosphere are
typically approximated from 2.0 to 2.5. The empirical drag coefficients are 2.3 and 2.5,
while the coefficient also has the close relation to the altitude of the satellite. The popular
density models, such as Jacchia 1971, MSIS-86 and DTM87 atmospheric density model,
were applied to compute the orbit for HY-2A. In order to search for the optimal coeffi-
cients for this mission, we have estimated the orbits using MSIS-86, Jacchia 1971 and
DTM87 model with the coefficients of 2.0, 2.1, 2.2, 2.3, 2.4 and 2.5 for HY-2A satellite
with 970 km altitude, respectively.
The statistics of RMS differences between the CNES orbit and the orbit derived
from MSIS-86 with different atmospheric drag coefficients are showed in Table 2.
Table 2 indicates that the smallest RMS value is 0.842 cm in radial direction
computed using coefficients of 2.0, 2.2, 2.3, respectively, and the biggest RMS value is
computed using MSIS-86 with coefficients 2.4 and 2.5. The smallest and the biggest
RMS values in along track direction are 2.682 and 2.893 cm achieved from MSIS-86
Performance of Three Atmospheric Density Models on Precise … 131
Table 2. Statistics of orbit difference between the CNES orbits and the derived orbits using
MSIS-86 model with different coefficients (cm)
Direction Radial Along track Cross track 3D
Mean RMS Mean RMS STD RMS RMS
2.0 −0.118 0.842 0.690 2.682 2.701 2.702 3.899
2.1 −0.118 0.843 0.692 2.686 2.702 2.703 3.903
2.2 −0.118 0.842 0.695 2.688 2.704 2.705 3.905
2.3 −0.120 0.842 1.050 2.884 2.694 2.696 4.037
2.4 −0.120 0.843 1.049 2.887 2.695 2.697 4.040
2.5 −0.119 0.843 1.053 2.893 2.697 2.699 4.045
with coefficient 2.0 and 2.5 respectively. The smallest RMS in 3D is 3.899 cm
achieved using this density model with coefficient 2.0. From Table 2, we can notice
that the relative optimal accuracy of the orbit is derived using the empirical atmospheric
density model MIS-86 with drag coefficient 2.0, which has the closest agreement with
the precise CNES orbit, and the differences between both orbits are listed in Fig. 2.
Therefore, 2.0 is the relative optimal coefficient for MSIS-86 model for HY-2A
satellite. At the same time, we can know that there is little change with the RMS in
radial direction and cross track direction from the coefficient 2.0 to 2.5, and there are
large changes in along track direction. The main reasons are that the atmospheric drag
acts on the satellite conversely along the track direction, and it is in proportion to the
atmospheric density, and the change of drag coefficient will bring change of atmo-
spheric drag.
1.82 Radial
1.32
0.82
0.32
-0.18
-0.68
-1.18
-1.68
-2.18
Orbit differences (cm)
Epoch
Fig. 2. The orbit difference between the derived orbit using MSIS-86 model with coefficient 2.0
and the CNES orbit
132 Q. Kong et al.
The statistics of differences between the CNES orbits and the orbits derived from
Jacchia 1971 with different atmospheric drag coefficients are listed in Table 3.
Table 3. Statistics of orbit difference between the CNES orbits and the derived ones using
Jacchia 1971 model with different coefficients (cm)
Direction Radial Along track Cross track 3D
Mean RMS Mean RMS Mean RMS RMS
2.0 −0.120 0.842 1.146 2.933 0.102 2.693 4.070
2.1 −0.121 0.842 1.148 2.936 0.102 2.694 4.072
2.2 −0.121 0.842 1.149 2.939 0.102 2.695 4.075
2.3 −0.121 0.842 1.150 2.942 0.100 2.696 4.079
2.4 −0.120 0.843 1.049 2.887 0.102 2.697 4.040
2.5 −0.120 0.842 1.048 2.880 0.102 2.700 4.037
Table 3 indicates that the smallest RMS value is 0.842 cm in radial direction
computed using coefficients of 2.0, 2.1, 2.2, 2.3 and 2.5, respectively, and the biggest
RMS is 0.843 cm derived from the Jacchia 1971 with coefficient 2.4. The smallest
RMS value in along track direction is 2.880 cm derived from this model with coeffi-
cient 2.5, and the biggest one is 2.942 cm achieved with coefficient 2.3. The smallest
RMS in 3D is 4.037 cm achieved using this density model with coefficient 2.5. From
Table 3, we can see that the relative optimal accuracy of the orbit is derived using the
empirical atmospheric density model Jacchia 1971 with drag coefficient 2.5. Therefore,
2.5 is the relative optimal coefficient for Jacchia 1971 model for HY-2A satellite. At the
same time, we can know that there are the smallest changes with the RMS in radial
direction and cross track direction, and there is the largest change in along track
direction.
The statistics of differences between orbits derived using DTM87 with different
atmospheric drag coefficients with respect to the CNES orbits are listed in Table 4.
Table 4. Statistics of orbit differences between the CNES orbits and the orbits derived using
DTM87 model with different coefficients (cm)
Direction Radial Along track Cross track 3D
Mean RMS Mean RMS Mean RMS RMS
2.0 −0.120 0.845 1.157 2.951 0.109 2.692 4.083
2.1 −0.121 0.845 1.157 2.951 0.109 2.692 4.083
2.2 −0.121 0.846 1.159 2.956 0.108 2.694 4.088
2.3 −0.121 0.845 1.159 2.959 0.107 2.695 4.090
2.4 −0.120 0.846 1.162 2.964 0.107 2.695 4.094
2.5 −0.120 0.846 1.162 2.967 0.107 2.697 4.098
Performance of Three Atmospheric Density Models on Precise … 133
Table 4 shows that the smallest RMS value is 0.845 cm in radial direction com-
puted using coefficients 2.0, 2.1 and 2.3, respectively, and the biggest one is 0.846 cm
derived from DTM87 with coefficients 2.2, 2.4 and 2.5, respectively. The smallest
RMS value in along track direction is 2.951 cm derived with Jacchia 1971 with
coefficient 2.0 and 2.2, respectively, and the biggest one is 2.967 cm achieved using
this model with coefficient 2.5. The smallest RMS in 3D is 4.083 cm computed using
this density model with coefficients 2.0 and 2.1, respectively. From Table 4, we can see
that the relative optimal accuracy of the orbit is derived using the empirical atmospheric
density model DTM87 with drag coefficient 2.0. At the same time, we can know that
the RMS in radial direction and cross track direction vary little between the results
derived using different coefficients, and vary largely in along track direction.
5 Conclusions
The work was performed to determine the relative optimal drag coefficients and drag
density model for the HY-2A satellite by using DORIS data. In order to obtain the
optimal solution, we compared the performance of different coefficients of atmospheric
drag using different atmospheric density models. We computed orbits using MSIS-86,
Jacchia 1971 and DTM87 atmospheric models applying DORIS technique, respec-
tively, and analyzed the orbit difference RMS with respect to CNES orbits. The results
indicate that the empirical atmospheric density model MSIS-86 performed significantly
better than the other two density models. Among the coefficients of 2.0, 2.1, 2.2, 2.3,
2.4 and 2.5, we found that the relative optimal atmospheric coefficient is 2.0 for MSIS-
86 density model, and 2.5 for Jacchia 1971 density model. The smallest RMS differ-
ence with respect to the CNES orbit is 0.842 cm in radial direction, and the biggest
RMS is 2.702 cm in cross track direction computed using MSIS-86 with coefficient
2.0.
There are mainly three possible reasons for the RMS of orbit difference better than
1 cm derived from DORIS data with respect to the CNES orbit. Firstly, the DORIS
beacon stations are distributed evenly; secondly, there are large amount of measure-
ments of DORIS; finally there are more tracking stations. From the comparison
between different atmospheric models, we can know from Fig. 2 that at the beginning
and end of the arc, there are diverse phenomena, and this is mainly because the
uncertainties are significantly higher at the beginning of the arc than the remnant part,
and the integrator is not stable, and the uncertainties accumulate with the time goes by.
It should be recognized that DORIS, as another space geodetic technique besides
GPS, SLR and VLBI, can be used to achieve high precise orbit. In our future work, we
will put more research on the improvement of atmospheric density model and inves-
tigation of the relationship between atmospheric drag and the latitude, longitude,
altitude, the plate face and so on, and realize the POD for the other satellites of HY-2
series.
Acknowledgements. This work was supported by the National Natural Science Foundation of
China (Grant No. 41704015, 41774001), A Project of Shandong Province Higher Educational
Science and Technology Program (J17KA077), the Public Benefit Scientific Research Project of
134 Q. Kong et al.
China (Grant No. 201412001), International Science and Technology Cooperation Program of
China (Grant No. 2009DFB00130), the Basic Science and Technology Research Project of China
(Grant No. 2015FY310200), and Open Fund of Engineering Laboratory of Spatial Information
Technology of Highway Geological Disaster Early Warning in Hunan Province (Changsha
University of Science & Technology) (Grant No. kfj150605).
References
Balmino, G., et al.: Numerical integration techniques revisited. J. Geod. (1990)
Berger, C., et al.: Improvement of the empirical thermospheric model DTM: DTM94—a
comparative review of various temporal variations and prospects in space geodesy
applications. J. Geod. 72, 161–178 (1998)
CIRA.: COSPAR international reference atmospheric 1972. In: COSPAR Committee for CIRA.
Akademie Verlag, Berlin (1972)
Fattig, E.: Comparison of precision orbit derived density estimates for CHAMP and GRACE
satellites. Master’s thesis, University of Kansas, Kansas (2011)
Gao, F., et al.: Analysis of HY2A precise orbit determination using DORIS. Adv. Space Res.
(2015). https://doi.org/10.1016/j.asr.2014.11.032 (Springer)
Goad, C.C., Goodman, L.: A modified Hopfield tropospheric refraction correction model. In:
American Geophysical Union Annual Fall Meeting, San Francisco, USA (1974)
Guo, J., et al.: On precise orbit determination of HY-2 with space geodetic techniques. Acta
Geophys. (2013a). https://doi.org/10.2478/s11600-012-0095-8 (Springer)
Guo, J., et al.: Centimeter level orbit determination for HY2A using GPS data (in Chinese).
Geomatics Inf. Sci. Wuhan Univ. (2013b)
Hedin, A.E.: A revised thermospheric model based on mass spectrometer and incoherent scatter
data MSIS-83. J. Geophys. Res. (1983). https://doi.org/10.1029/ja088ia12p10170 (Springer)
Hedin, A.E.: MSIS-86 thermospheric model. J. Geophys. Res. (1987). https://doi.org/10.1029/
JA092iA05p04649
Jacchia, L.G.: Revised static models of the thermosphere and exosphere with empirical
temperature profiles. Smithsonian Astrophysical Observatory, Special Report (1971)
Kang, Z., et al.: Precise orbit determination for the GRACE mission using only GPS data.
J. Geod. 80, 322–331 (2006). https://doi.org/10.1007/s00190-006-0073-5
Knocke, P.C., et al.: Earth radiation pressure effects on satellites. In: Proceedings of the
AIAA/AAS Astrodynamics Conference, Minneapolis, Minnesota (1988)
Kong, Q., et al.: Precise orbit determination and accuracy analysis of HY-2A satellite using
DORIS Doppler data. Acta Geod. Geophys. (2014). https://doi.org/10.1007/s40328-014-
0066-4 (Springer)
Kong, Q.: Centimeter level precise orbit determination of HY-2a satellite using DORIS and SLR
techniques. Doctoral thesis, Shandong University of Science and Technology, Qingdao
(2015)
Lechtenberg, T.F.: Derivation and observability of upper atmospheric density variations utilizing
precision orbit ephemerides. Master’s thesis, University of Kansas, Kansas (2010)
Lechtenberg, T.F.: Density model corrections derived from orbit data to characterize upper
atmospheric density variations. Doctoral thesis, University of Kansas, Kansas (2015)
Li, P.: High accuracy orbit determination with DORIS system. J. Spacecr. TT&C Technol.
(2010). 1674-5620(2010)03-0058-07 (Springer)
Lin, M., et al.: Precise orbit determination technology based on dual-frequency GPS solution for
HY-2 satellite. Eng. Sci. (2014)
Performance of Three Atmospheric Density Models on Precise … 135
Liu, J.: Satellite orbit determination method with DORIS system. Master’s thesis. National
University of Defense Technology, Changsha, China (2013)
Locke, T.C.: Further developments in orbit ephemeris derived neutral density. Master’s thesis,
University of Kansas, Kansas (2012)
Lyard, F., et al.: Modeling the global ocean tides: modern insights from FES2004. Ocean Dynam.
(2006). https://doi.org/10.1007/s10236-006-0086-x (Springer)
Mance, S.: An investigation into satellite drag modeling performance. Master’s thesis, University
of Kansas, Kansas (2010)
McCarthy, D.D., Petit, G.: Equations of motion for an artificial earth satellite. In: IERS
Conventions 2003. IERS Technical Note 32, Paris, France (2003)
Mehta, P.M., et al.: Comparing physical drag coefficients computed using different gas-surface
interaction models. J. Spacecr. Rockets. (2014). https://doi.org/10.2514/1.a32566 (Springer)
Mercier, F., et al.: Jason-2 DORIS phase measurement processing. Adv. Space Res. (2010).
https://doi.org/10.1016/j.asr.2009.12.002 (Springer)
Montenbruck, O., Gill, E.: Satellite Tracking and Observation Models, Satellite Orbits. Springer,
Berlin, Heidelberg, German (2000)
Pavlis, N.K., et al.: The development and evaluation of the earth gravitational model 2008
(EGM2008). J. Geophys. Res. (2012). https://doi.org/10.1029/2011jb008916 (Springer)
Petit, G., Luzum, B.: IERS conventions (2010). Bureau International des Poids et Mesures,
Sevres, France (2010)
Pilinski, M.D. et al.: Drag coefficients of satellites with concave geometries: comparing models
and observations. J. Spacecr. Rockets (2011). https://doi.org/10.2514/1.50915 (Springer)
Rim, H.J.: TOPEX orbit determination using GPS tracking system. Doctoral thesis, University of
Texas at Austin, Austin, USA (1992)
Robert, C.E.: Densities based upon Jacchia’s 1971 models. Celest. Mech. 4, 368–377 (1971)
Schlegel, K., et al.: Thermospheric density structures over the polar regions observed with
CHAMP. Ann. Geophys. (2005)
Standish, E.M.: JPL planetary and lunar ephemerides, DE405/LE405. Jet Propulsion Laboratory,
InterOffice Memorandum, IOM (1998)
Wang, X., et al.: Precise orbit determination technology based on SLR solution for HY-2 satellite
(in Chinese). Eng. Sci. (2014)
Willis, P., et al.: Recent improvements in DORIS data processing at IGN in view of ITRF2008,
the ignwd08 solution. In: Geodesy for Planet Earth International Association of Geodesy
Symposia (2012). https://doi.org/10.1007/978-3-642-20338-1_6 (Springer)
Willis, P. et al.: DPOD2008, A DORIS-oriented terrestrial reference frame for precise orbit
determination. In: Rizos, C., Willis, P. (eds.) International Association of Geodesy Symposia,
p. 143. Springer, Switzerland (2013)
Zelensky, N.P., et al.: DORIS/SLR POD modeling improvements for Jason-1 and Jason-2. Adv.
Space Res. (2010). https://doi.org/10.1016/j.asr.2010.05.008 (Springer)
In Situ Test of Traffic-Load-Induced
Settlement of Alluvial Silt Subsoil Treated
by Unslaked Lime
Abstract. The soft wet alluvial silt is widely distributed material in the world.
In order to improve the bearing capacity and decrease the traffic-load-induced
settlement of silt subsoil, the shallow subsoil always treats with the unslaked
lime. However, the mitigating effect of this ground treatment method on traffic-
load-induced settlement of alluvial silt subsoil is inconclusive. Therefore, with
the developed falling weight simulation equipment of traffic load, in situ tests
are carried out on the natural and unslaked lime treated alluvial silt ground in the
Yellow River delta of China to study traffic-load-induced settlement. Further-
more, Chai-Miur cumulative deformation model of soil is employed to
numerically simulate the long-term traffic-load-induced cumulative settlement.
In situ test results indicate that because unslaked lime treatment enhances wave
impedance of the reinforced soil layer, the wheel-load-induced dynamic stress
and excess pore water pressure in the substratum decreases. The decrease of
excess pore water pressure reduces the cumulative settlement of unslaked lime
treated subsoil. For short-term cumulative settlement, there are differences
between tested and calculated results, but not much. The calculated results imply
that after opening to traffic for 10 years, compared with the natural ground, the
cumulative settlement of the unslaked lime treated subsoil reduces by about
21%, and the change of transverse slope of pavement induced by cumulative
settlement decreases by 1/3. In situ test and numerical calculation results
demonstrate that shallow layer treatment with unslaked lime can effectively
mitigate the cumulative settlement of alluvial silt subsoil.
1 Introduction
Alluvium is loose soil, which has been eroded and reshaped by water in some form,
and typically made up of a variety of materials, including fine particles of silt and clay
and larger particles of sand and gravel. Alluvium is widely distributed in the world. In
China, the Yellow River delta is the youngest large river delta due to frequent river
diversion in history, and in the Yellow River delta, the alluvium is mainly low liquid
limit silt. The alluvial silt is newly formed under consolidated deposit and has poor
engineering performances with low liquid limit and plasticity index, small cohesion,
low strength, intensive capillarity, poor graduation and water stability and high
liquefaction potential. Furthermore, because the alluvial silt is liquefiable and the
groundwater level is very high in the region, the traffic-load-induced cumulative set-
tlement of the road with low embankment is significant after opening to traffic, and this
can cause serious pavement diseases. In order to improve the bearing capacity of the
soft wet subsoil in the region and reduce its cumulative settlement, shallow layer
treatment with unslaked lime is usually employed. The improving effect of unslaked
lime treatment on the bearing capacity of soft wet subsoil is undeniable (Kamon 1992;
Narasimha Rao and Rajasekaran 1996). However, at present the mitigating effect of
hard crust formed by unslaked lime treated subsoil on the traffic-load-induced cumu-
lative settlement is inconclusive. Note that in this study, the pavement means the hard
layered structure that forms a road carriageway, the embankment means the compacted
soil layer below the pavement of a road, and the subsoil means the natural soil below a
road embankment; the term subgrade includes all layers above the natural ground
surface.
Previous studies on the cumulative settlement of subsoil focus on dynamic defor-
mation test of soil in the laboratory and numerical simulation. In the early 1950s, Seed
and his co-workers (Seed et al. 1955; Seed and McNeill 1956) studied the settlement
behaviour of road under repeated loads by normal compression tests. Monismith et al.
(1975) analyzed the characteristics of permanent deformation of subsoil due to repeated
loadings based on repeated-load triaxial compression tests. Yildirim and Ersan (2007)
studied consolidation settlements of soft clay by undrained cyclic simple shear tests in
the laboratory. Shahu and Yudhbir (2008) studied cumulative plastic strain of a quasi-
saturated compacted silty clay under cyclic load by monotonic and cyclic undrained
triaxial tests. Cyclic triaxial tests were conducted by Liu and Xiao (2010) to study the
behaviour of silt subsoil under various physical states and stress conditions. Chai and
Miura (2002) modified the Li-Selig model (Li and Selig 1996) and calculated the
permanent settlement of road with a low embankment on soft subsoil. Akira et al.
(2003) analyzed the traffic-load-induced settlement of low embankment road on silty-
clay. Abdelkrim et al. (2003) adopted a general structure analysis approach to predict
the traffic-load-induced residual settlement. Fujiwara and coworkers (Fujiwara et al.
1985; Fujiwara and Ue 1990) researched the effect of preloading on the cumulative
settlement of clay under cyclic loadings, and they concluded that the amount of set-
tlement that occurs after construction depends strongly on the soil over consolidation
ratio, degree of consolidation at the time of unloading, static loading magnitude, and
repeated loading magnitude. Chai and Miura (2002) calculated the permanent settle-
ment of the improved subsoil with soil-cement columns. Indraratna et al. (2010) carried
out a large-scale triaxial testing on clay-drain complex subjected to cyclic loading
representing a typical track environment to predict the behaviour of the soft estuarine
subsoil with short vertical drains. However, currently there has been no studies on the
cumulative settlement of unslaked lime treated subsoil.
Though the laboratory tests can evaluate the traffic-load-induced cumulative
deformation of subsoil, they are not accurate because the stress and boundary condi-
tions are more complex in real practices. The numerical method, although it has the
advantage of analyzing the subsoil settlement under complicated stress conditions, has
its shortcomings, such as long calculating time and uncontrollable cumulative calcu-
lation errors, which limit its practical use in engineering. In addition, the road
138 Q. Jin et al.
settlement observed in the field after opening to traffic is the sum of the settlement
induced by the weight of subgrade and traffic-load-induced settlement. However, there
is no effective method to separate them and analyze individually traffic-load-induced
cumulative settlement. Therefore, it is significant to conduct in situ tests for directly
simulating the traffic-load-induced settlement of subsoil. In this paper, with the
developed Falling-Weight Traffic Load Simulation Equipment (FWTLE), in suit tests
are performed to study the short-term traffic-load-induced cumulative settlement of
unslaked lime treated alluvial silt subsoil in the Yellow River delta. And by the
numerical calculations, the long-term cumulative settlement is studied. Comparing the
cumulative settlement of the natural ground with that of unslaked lime treated subsoil
demonstrates the validity of unslaked lime treatment method.
2 In Suit Tests
In order to improve the bearing capacity of the subsoil and decrease the traffic-load-
induced cumulative settlement, the design requires that the 60 cm thick unsaturated
topsoil is treated with 6% unslaked lime powder. The physical mechanical properties of
unslaked lime treated soil are shown in Table 2.
time of wheel load (i.e. vehicle speed). Therefore, the size of square loading plate used
in situ test is determined as 1.2 m 1.2 m, and its thickness is 3 mm. And the radius
of its equivalent area circle is 0.68 m, which can meet the test precision requirement in
the case of the zero-fill embankment. The falling weight in the test is the hammer used
in the standard penetration test (SPT) for the engineering geological investigation and it
is 62.5 kg in weight. Before in situ test, by adjusting the height of falling weight and
the stiffness of air spring, different vertical stress under the loading plate can be
obtained to meet the requirements of different engineering cases. The amplitude and
duration of vertical dynamic stress response obtained in tests should be approximately
consistent with the computed results in the cases of different heights of the embank-
ment and different vehicle speeds.
Because the cumulative settlement of the zero-fill road embankment is significantly
larger, in this study, the falling distance and the stiffness of air spring were adjusted to
simulate the case of the zero-fill embankment.
cyclic loading, the parameters such as stress, displacement and pore pressure were
tested. Because the vehicle-load-induced excess pore pressure need a long time to
dissipate, the settlement and pore pressure after cyclic loading were also tested. The
total cumulative settlement was the sum of the settlement during and after loading.
The subsoil settlement was tested with the laser displacement sensors. Four
dynamic soil pressure sensors were installed at different characteristic positions
between loading plate and subsoil surface. Two dynamic pore pressure sensors were
installed at 0.7 and 1.2 m underground, respectively, below the center of loading plate.
And in the same way, two dynamic soil pressure sensors were installed at 0.6 and
1.0 m underground, respectively (Fig. 3).
60
60
70
40
50
Unit: cm
20
30
Vertical stress /kPa
20
10
Fig. 4. Vertical stress response on the subsoil surface under the center of loading plate
Compared with the natural silt subsoil, it can be seen from Fig. 4b that the stress
amplitude on the unslaked lime treated subsoil surface increases. This is because after
mixing the soft wet soil with the reasonable amount of unslaked lime, a series of
reactions take place: water absorption, exothermic action and expansive action; ion
exchange action; carbonation (chemical cementation reaction); pozzolanic action
(chemical gelation reaction) and crystallization action. These reactions can make the
moisture in soil reduce, slit particles coagulate to form larger aggregates. With the
gradual hardening of subsoil, soil particles bonded together and the physical and
mechanical properties of the subsoil are improved. With the subsoil treated by unslaked
lime, the wave impedance of shallow subsoil increase, and the hard curst effect causes
the increase of the dynamic stress amplitude on the subsoil surface.
Figure 5 shows that the variation curves of vertical dynamic stress amplitude of the
subsoil under the center of the loading plate with depth. It is can be seen that compared
with the natural ground, the vertical stress amplitude in the ground treated by unslaked
lime more quickly decreases, and this makes the dynamic stress in the substratum of
hard crust significantly reduce. The decrease of the stress level in the substratum can
induce the reduction of the cumulative deformation.
In Situ Test of Traffic-Load-Induced Settlement … 143
32
16
0
0.0 0.2 0.4 0.6 0.8 1.0
Depth /m
Figure 7 shows the variations of the excess pore water pressures in the natural and
unslaked lime treated ground with the load numbers N during the loading process. It
can be seen that, compared with the natural ground, the excess pore water pressure in
the soft substratum of the unslaked lime treated ground is obviously lower. This is
because with the ground treated by unslaked lime, the dynamic stress level in the
substratum decreases (as shown in Fig. 5).
At the initial stage of loading, the pore water pressure in the soft substratum of the
unslaked lime treated ground develops faster than that in the natural ground. Never-
theless, after certain load numbers, the pore water pressure in unslaked lime treated
ground does not increase but gradually declines to a stable state at last. By contrast, the
pore water pressure in natural ground continuously increases with the load numbers.
144 Q. Jin et al.
16
Depth: 0.7m
12
p /kPa
8
4 Natural ground
Unslaked lime treated subsoil
0
0 1500 3000 4500
N
(a) 0.7 m depth under ground
16
Depth: 1.2m
12
p /kPa
4
Natural ground
Unslaked lime treated subsoil
0
0 1500 3000 4500
N
(b) 1.2 m depth under ground
Fig. 7. Variation curves of excess pore water pressure with load numbers
The difference above is mainly because of the different drainage conditions of the
unslaked lime soil treated and natural grounds during the loading process. Because the
permeability coefficient of unslaked lime treated soil is smaller than the natural soil, the
soil-lime hard crust can block the dissipation of pore water pressure in the substratum.
Therefore, at the beginning of cyclic loading, the pore water pressure in unslaked lime
treated ground increases faster than the one in the natural ground. But as the pore water
pressure increases, the strength of the substratum decreases, and this induces the
cracking of the soil-lime hard crust under the traffic load (observed in tests). The cracks
provide good drainage channels. As mud comes out through the channels, the pore
water pressure in the substratum of hard crust does not increase, but decreases.
However, for the natural ground, the flexibility of top soil is larger, so its main failure
mode is not local cracking but the plastic deformation. Therefore, the seepage paths
change little in the process of cyclic loading, and this makes the pore water pressure
continuously increase.
Figure 8 shows the variations of excess pore water pressure in natural and unslaked
lime treated grounds after terminating loading. The excess pore water pressure sharply
declines firstly, and then stabilizes for some time, finally declines to hydrostatic
pressure. This phased variation implies that after terminating loading, the seepage path
of water and dissipation process of pore water pressure are complex.
After terminating loading, compared with the natural ground, the dissipation of
excess pore water pressure in unslaked lime treated ground needs more time. The
dissipation of excess pore water pressure in natural ground needs seven hours, but the
one in unslaked lime treated ground takes more than twenty hours. This is because after
In Situ Test of Traffic-Load-Induced Settlement … 145
16
12
p /kPa
8
Depth: 0.7m
4 Natural ground
Unslaked lime treated subsoil
0
0 5 10 15 20 25
t /h
(a) 0.7 m depth under ground
15
12
p /kPa
6 Depth: 1.2m
3 Natural ground
Unslaked lime treated subsoil
0
0 5 10 15 20 25
t /h
(b) 1.2 m depth under ground
terminating cyclic loading, the mud hardens in the cracks of soil-lime hard crust
induced by cyclic loading, and this makes the hard crust becomes impervious. How-
ever, the practical traffic load is cyclic and continual, the cracks in hard crust are not
blocked. So the delay effect of dissipating the pore water pressure does not occur.
3.0
Natural ground
Settlement /cm
2.4 Unslaked lime treated subsoil
1.8
1.2
0.6
0.0
0 1500 3000 4500
N
Figure 10 shows the variations of cumulative settlement with time after terminating
loading. With the dissipation of excess pore pressure, the cumulative settlement
gradually increases, and tends to be stable finally. Furthermore, after terminating
loading, the cumulative settlement of unslaked lime treated ground has a little differ-
ence with that of natural ground. But the cumulative settlement of unslaked lime treated
subsoil takes more time to reach steady state. This is because, after terminating loading,
the hardening of mud in the soil-lime hard crust blocks the dissipation channels of the
excess pore water pressure (as seen in Fig. 8).
0.8
Settlement /cm
0.6
0.4
Natural ground
Unslaked lime treated subsoil
0.2
0.0
0 4 8 12 16 20 24
t /h
Fig. 10. Variation curves of cumulative settlement with time after terminating loading
1
qf ¼ ccu cos /cu =ð1 sin /cu Þ þ ð1 þ K0 Þrcz sin /cu =ð1 sin /cu Þ ð2Þ
2
where K0 is the static soil pressure coefficient at rest; rcz is overlying soil pressure.
a, b, m and n in Eq. (1) reflect the combined effect of stress state, physical state and
types of soil. Series of triaxial tests, the constants in Eq. (1) were obtained: a = 0.64,
b = 0.10, m = 1.70 and n = 1.00.
3.1 Comparisons and Analyses of the In Situ Test and Numerical Results
Before calculating long-term cumulative settlement of subsoil, in order to prove the
validity and feasibility of calculation method, the cumulative settlement of natural
ground was numerically calculated by simulating in situ test.
Firstly, the static module of the finite difference program Flac3D was used to
calculate initial static deviatoric stress qs in the subsoil. Secondly, the dynamic module
of Flac3D was employed to simulate the dynamic response of subsoil. In the dynamic
calculation, the vertical dynamic stresses obtained by soil pressure sensors (Fig. 3a)
under the loading plate in the tests were loaded on the subsoil surface. The dynamic
deviatoric stress qd can be obtained from the peaks of the dynamic stresses in three
orthogonal directions. Finally, qs and qd were taken into Eq. (1), and cumulative
deformation of subsoil can be calculated. The static strength of soil qf is calculated by
Eq. (2). The cumulative settlement on the subsoil surface can be obtained by inte-
grating the vertical cumulative deformation of soil along depth.
Figure 11 shows the comparisons of cumulative settlements from in situ tests and
numerical simulations. It can be seen that the variation trend of testing cumulative
settlement with the load numbers is basically consistent with calculated results.
However, the calculated settlement develops faster than the tested results at the early
stage of cyclic loading. There are many reasons for the differences between the tested
and calculated settlements. The main reason is that the subsoil is partially drained at the
site, i.e. the development and dissipation of excess pore water pressure in the subsoil
148 Q. Jin et al.
2.5
Settlement (cm)
2.0
1.5
1.0 Tested results
Calculated results
0.5
0.0
0 1000 2000 3000 4000 5000
Loading numbers
35 kN 65 kN
D A
1.8 m
C B
35 kN 65 kN
4.0 m
simplified as four point loads on the pavement shown in Fig. 12. The loading mode of
wheel load proposed by Huang (1993) is adopted:
(
p
F ¼ Fmax sin2 t 0 t Ts
Ts ð3Þ
F¼0 t [ Ts
where t is time; Fmax is the wheel load peak, 65 kN for the front wheel and 35 kN for
the rear wheel; Ts is the duration time of single vehicle load, and has an inverse
relationship with the vehicle speed. Herein Ts is taken as 0.031 s, representing the
equivalent vehicle speed of 120 km/h (Huang 1993).
The settlement calculation considers the variation of annual traffic volume. The
traffic volumes at eight flyovers along the Xinhe-Xinzhuangzi highway were investi-
gated. The predicted cumulative traffic volume of trucks is:
73N1
N¼ ð1 þ cÞt 1 ð3Þ
c
where N1 = 13463 is the annual average daily traffic at the early stage after opening to
traffic; r = 5.775% is the average annual growth rate; t is time in year.
The physical and mechanical parameters of pavement, embankment and subsoil are
shown in Tables 1, 2 and 3. In this study, for the unslaked lime treated ground, the
cumulative settlement is only from the deformation of the natural soil layer, and
the deformation of soil-lime hard crust was ignored in the calculation. Figure 13 shows
16
12 Time
εp/(10-3)
1 year
8 2 years
10 years
4
0
0 4 8 12 16
Depth/m
the variations of cumulative deformations with the depth (one year, two years and ten
years after opening to traffic). It can be seen that the cumulative deformation sharply
reduces within 5 m underground, and the decreasing rate becomes slow beyond 5 m
underground. This illustrates that the traffic-load-induced settlement is mainly from the
cumulative deformations of silts within 5 m underground.
Figure 14 shows the variations of the cumulative settlement of subsoil with time
after opening to traffic. It can be seen that the cumulative settlement rapidly develops in
the initial stage and then gradually gets slow. Compared with the natural ground, the
cumulative settlement of unslaked lime treated ground evidently reduces. For example,
the cumulative settlement in ten years reduces by 21.4%. This demonstrates that the
unslaked lime treatment obviously mitigates the traffic-load-induced settlement, and
this is consistent with the results of in situ tests.
6
Settlement /cm
4
Natural ground
2 Unslaked lime treated subsoil
0
0 2 4 6 8 10
t /year
8
Carriageway
Hard shoulder
Overtaking lane
Settlement /cm
6 Left marginal Earth shoulder
strip
Median
strip
4
Natural ground
Unslaked lime treated subsoil
2
0 2 4 6 8 10 12 14 16
Distance from road midline /m
Fig. 15. Transverse distribution curves of cumulative settlement after opening to traffic for
10 years
4 Conclusions
In order to prove the mitigating effect of unslaked lime treatment on the traffic-load-
induced settlement of soft wet alluvial silt subsoil, in situ tests were conducted to study
the short-term settlements of natural and treated grounds in the Yellow River delta with
the developed FWTLE. In addition, based on the Chai-Miura cumulative deformation
model, the long-term cumulative settlements of the subsoils were analyzed numerically.
The following main conclusions are drawn out:
(1) Compared with the natural ground, the wave impendence of unslaked lime treated
ground increases. This makes traffic-load-induced dynamic stress in substratum
significantly reduce. Moreover, the excess pore water pressure in the substratum of
unslaked lime treated subsoil is less than that in the natural ground.
(2) In comparison with the natural ground, the cumulative settlement of the unslaked
lime treated ground obviously decreases. This proves the mitigating effect of
unslaked lime treatment on the cumulative settlement of alluvial silt subsoil.
(3) The cumulative settlement of subsoil aggravates the development of pavement
diseases. Unslaked lime treatment of subsoil can effectively decrease the harms of
traffic loads on the road with low embankment.
Acknowledgements. This work was supported by the Chinese Natural Science Foundations
(Nos. 51279094, 51078222 and 50708056), the Natural Science Foundations of Shandong
Province, China (No. ZR2011EEM012) and the Independent Innovation Foundation of Shan-
dong University (IIFSDU) (No. 2012HW003).
References
Abdelkrim, M., Bonnet, G., Buhan, P.D.: A computational procedure for predicting the long term
residual settlement of a platform induced by repeated traffic loading. Comput. Geotech. 30(6),
463–476 (2003)
Akira, S., Lawalenna, S., Norihiko, M.: Partially-drained cyclic behavior and its application to
the settlement of a low embankment road on silty-clay. Jpn. Geotech. Soc. 43(1), 33–46
(2003)
152 Q. Jin et al.
Chai, J.C., Miura, N.: Traffic-load-induced permanent deformation of road on soft subsoil.
J. Geotech. Geoenviron. Eng. 128(11), 907–916 (2002)
Cui, X.Z.: Traffic-induced settlement of subgrade of low liquid limit silt in Yellow River delta.
China Civil Eng. J. 45(1), 154–162 (2012). (in Chinese)
Fujiwara, H., Ue, S., Yasuhara, K.: Consolidation of alluvial clay under repeated loading. Soils
Found. 25(3), 19–30 (1985)
Fujiwara, H., Ue, S.: Effect of preloading on post-construction consolidation settlement of soft
clay subjected to repeated loading. Soils Found. 30(1), 76–86 (1990)
Huang, Y.H.: Pavement Analysis and Design. Pearson Education, Delhi (1993)
Indraratna, B., Rujikiatkamjorn, C., Ewers, B.: Class A prediction of the behavior of soft
estuarine soil foundation stabilized by short vertical drains beneath a rail track. J. Geotech.
Geoenviron. Eng. 136(5), 686–696 (2010)
Kamon, M.: Recent developments of soil improvement. In: Proceedings of International
Symposium on Soil Improvement and Pile Foundation, Nanjing, China, vol. I, pp. 1–16
(1992)
Li, D., Selig, E.T.: Cumulative plastic deformation for fine-grained subgrade soils. J. Geotech.
Eng. 122(12), 1006–1013 (1996)
Liu, J.K., Xiao, J.H.: Experimental study on the stability of railroad silt subgrade with increasing
train speed. J. Geotech. Geoenviron. Eng. 136(6), 833–841 (2010)
Monismith, C.L., Ogawa, N., Freeme, C.R.: Cumulative deformation characteristics of subsoil
due to repeated loading. Transp. Res. Rec. 537, 1–17 (1975)
Narasimha, R.S., Rajasekaran, G.: Reaction products formed in lime-stabilized marine clays.
J. Geotech. Eng. 122, 329–336 (1996)
Seed, H.B., Chan, C.K., Monismith, C.L.: Effects of repeated loading on the strength and
deformation of compacted clay. Highw. Res. Board Proc. 34, 541–558 (1955)
Seed, H.B., McNeill, R.L.: Soil deformation in normal compression and repeated loading test.
Highw. Res. Board Bull. 141, 44–53 (1956)
Shahu, J.T., Yudhbir, Hayashi S.: Cumulative plastic strain and threshold stress of a Quasi-
saturated compacted silty clay. Lowland Technol. Int. 10(2), 10–20 (2008)
Shen, Z.J.: Earth pressure of clay based on effective consolidation stress theory. Chi-
nese J. Geotech. Eng. 22(3), 353–356 (2000). (in Chinese)
Wang, X.: Test on dynamic stress of roadbed and pavement under heavy loads. J. Vib. Shock 26
(2), 169–173 (2007). (in Chinese)
Yildirim, H., Erşan, H.: Settlements under consecutive series of cyclic loading. Soil Dyn. Earthq.
Eng. 27, 577–585 (2007)
Experimental Study on Shear Strength
Behavior of Glass Fiber-Reinforced Sand
1 Introduction
stress for fiber-reinforced soil (Falorca and Pinto 2011). The effect of initial soil density
is not significant on the shear strength of fiber-reinforced soil reinforced with com-
paratively high fiber content (Li and Zornberg 2013). Addition of polypropylene fibers
to cemented soil increases peak strength, residual shear strength and energy absorption
capacity but it reduces the initial stiffness and brittleness index (Hamidi and Hoores-
fand 2013). Tensile strength of fiber-reinforced soil increases with an increase in dry
density and decreases with an increase in water content (Li et al. 2014). Fiber-
reinforcement reduces the amount of particle breakage of sand, more significantly for
well graded sand (Pino and Baudet 2015). Dry fiber reinforced sand provides same
shear strength that of heavily compacted unreinforced moist sand (Eldesouky et al.
2016). Claria and Vettorelo (2016) found that the addition of fibers increases the shear
strength and ductility but decreases the initial stiffness of sand. The shear strength was
increased indefinitely with fiber length.
The objective of the present study is to investigate the influence of glass fibers on
the shear strength characteristics of a locally available sand at varying relative density
for use as a construction material.
The experimental study was conducted with sand collected from the deposits of the
nearby Brahmaputra River. The sand can be designated as SP as per the Unified Soil
Classification Systems (USCS). The specific gravity of sand was 2.69 with a uniformity
coefficient of 1.59 and coefficient of curvature of 0.97. The maximum and minimum
void ratios of the sand were 0.96 and 0.63, respectively. Glass fiber of 0.15 mm
diameter, 10 mm length, 1.5 GN/m2 tensile strength, and 2.57 specific gravity was
used as reinforcement. Varying fiber contents (0.5, 1, 2, 3, and 4% by dry weight of
soil) were mixed with the sand along with 2% water content addition for producing a
homogeneous mix.
Prior to triaxial tests, compacted sand-fiber specimens of 38 mm diameter and
76 mm height were prepared on the base plate of the triaxial cell by using a split mould
at three different relative densities (Dr = 35, 65 and 85%) of the sand. Unreinforced
sand specimens were also prepared in a similar manner. Consolidated drained triaxial
tests were conducted on all specimens as per ASTM D4767-11 (2011) under different
confinement pressures (r3) ranging from 100 to 400 kPa which reflect the ranges
typical of relevant geoengineering applications. Shearing was done with a strain rate of
0.12 mm/min.
1200
fc = 4% Dr = 35%
fc = 3% L = 10 mm
1000 σ3 = 100 kPa
fc = 2%
400
200
0
0 5 10 15 20
Axial strain (%)
1200
fc = 4% Dr = 65%
fc = 3% L = 10 mm
1000 σ3 = 100 kPa
fc = 2%
Deviator stress (kPa)
fc = 1%
800
fc = 0.5%
No fiber
600
400
200
0
0 5 10 15 20
Axial strain (%)
relative densities. There is also an increase in failure strain indicating that the ductility
of the reinforced sand has improved. The reinforcing action of the fibers is primarily
due to local stitching of the potential slip planes through interlocking resistance.
The stress-strain plots of specimen with 35% relative density (Fig. 1) have not
shown any clear peak at any fiber content. However for reinforced specimens com-
pacted at 85% relative density (Fig. 3), the stress-strain plots show clear peak even at
low fiber content (0.5%). The reinforced specimen with 65% relative density has shown
an intermediate behavior (Fig. 2). In addition, the plots of the reinforced sand in the
post-yield stage have exhibited strain-softening behavior and this behavior is more
prominent with increasing specimen density. However, the residual strength is still
greater than that of the sand alone at any density level. This indicates that the glass
156 S. K. Patel and B. Singh
1200
Dr = 85% fc = 4%
L = 10 mm fc = 3%
1000 σ3 = 100 kPa fc = 2%
fc = 1%
Deviator stress (kPa) 800 fc = 0.5%
No fiber
600
400
200
0
0 5 10 15 20
Axial strain (%)
fiber-reinforced sand can bear greater load even at larger settlement compared to parent
soil.
The initial stiffness of specimens with 65 and 85% relative density is decreased
with fiber addition while it increases for 35% relative density. With addition of fibers,
peak appears at strain higher than 5% and it increases with fiber content. The peak
becomes more expressive at higher fiber as the density increases (Fig. 3). Post-peak
strength reduction occurs for all fiber contents and it reaches its residual strength
around 15% axial strain for specimen with 85% density, and at 20% axial strain for 35
and 65% relative density. At this stage, the resistance action of the fibers is lost partially
through pullout. The improvement in strength reaches optimum at 3% fiber content and
4% fiber-reinforced sand shows only marginal strength gain. Similar trends were found
at all confinement pressures, wherein the peaks in the stress-strain responses are more
prominent when the fiber content is higher.
Figure 4 shows a comparison of the stress-strain behavior of reinforced sand with
3% fiber at different relative densities under low (r3 = 100 kPa) and high confinement
pressure (r3 = 400 kPa). The initial stiffness of the reinforced sand increases with
either the increase of soil relative density or confinement pressure keeping the other
constant indicating that the specimen density has major role in stiffness improvement.
However, the peak strength is found to be occurring either at a lower axial strain with
increase in relative density at the same confinement pressure, or at a higher axial strain
with increase in confinement pressure at the same relative density. The strength at
higher strain level comes closer to each other irrespective of the initial specimen
density.
Experimental Study on Shear Strength Behavior of Glass … 157
1800
Dr = 85%, σ3 = 400 kPa
1500 Dr = 65%, σ3 = 400 kPa
Dr = 35%, σ3 = 400 kPa
Deviator stress (kPa)
1200
Dr = 85%, σ3 = 100 kPa
900 Dr = 65%, σ3 = 100 kPa
Dr = 35%, σ3 = 100 kPa
600
fc = 3%
300
0
0 5 10 15 20
Axial strain (%)
Fig. 4. Effect of relative density and confinement pressure on stress-strain behavior of sand
specimens of 65 and 85% relative density, friction angle marginally decreases as fiber
content increases. Improvement of cohesion is due to the fact that at the time of
loading, presence of fibers reduces the movement of soil particles by sharing the stress
coming on them and also provides some extra internal confinement to the specimen.
The sharing of load by fibers provides increased tensile strength to the sand matrix and
the extra confinement effect of fiber increases the specimen load carrying capacity in
terms of increased shear strength.
Experimental Study on Shear Strength Behavior of Glass … 159
1000
fc = 4% Dr = 85%
fc = 3% L = 10 mm
800 fc = 2%
fc = 1%
600 fc = 0.5%
τ (kPa) No fiber
400
200
0
0 200 400 600 800 1000
σn (kPa)
Fig. 6. Typical failure modes of unreinforced and reinforced sand specimens (fc = 3%):
a unreinforced with 35% density, b reinforced with 35% density, c unreinforced with 85%
density and d reinforced with 85% density
remains like a column and seems to penetrate in lower half section under loading. This
indicates that the fiber reinforcement contains the specimen failure in a relatively
smaller region.
4 Conclusions
With glass fiber reinforcement, the shear strength or the load carrying capacity of the
sand is found to have improved significantly. The post-failure residual strength of the
fiber-reinforced sand is also higher than that of the unreinforced sand. As the failure
strain increases, it can be concluded that through introduction of glass fibers, any
chance of sudden failure of the reinforced sand is significantly reduced. The sand with
3% optimum fiber content can be used for pavement subgrade and foundation beds. For
improving the specimen stiffness, it is better to increase specimen density. Inclusion of
fiber contains the specimen failure in relatively small zone by restricting the bulging of
sand specimen with 35% relative density and shear failure of dense sand specimen with
85% relative density.
References
ASTM D4767-11: Standard Test Method for Consolidated Undrained Triaxial Compression Test
for Cohesive Soils, ASTM International, West Conshohocken, PA (2011)
Benson, C.H., Khire, M.V.: Reinforcing sand with strips of reclaimed high-density polyethylene.
J. Geotech. Eng. ASCE 120(5), 838–855 (1994)
Claria, J.J., Vettorelo, P.V.: Mechanical behavior of loose sand reinforced with synthetic fibers.
Soil Mech. Found. Eng. 53(1), 12–18 (2016). https://doi.org/10.1007/s11204-016-9357-9
Experimental Study on Shear Strength Behavior of Glass … 161
Consoli, N.C., Festugato, L., Heineck, K.S.: Strain-hardening behavior of fiber-reinforced sand in
view of filament geometry. Geosynth. Int. 16(2), 109–115 (2009). https://doi.org/10.1680/
gein.2009.16.2.109
Eldesouky, H., Morsy, M.M., Mansour, M.F.: Fiber-reinforced sand strength and dilation
characteristics. Ain Shams Eng. J. 7(2), 517–526 (2016). https://doi.org/10.1016/j.asej.2015.
06.003
Falorca, I.M.C.F.G., Pinto, M.I.M.: Effect of short randomly distributed polypropylene
microfibers on shear strength behavior of soils. Geosynth. Int. 18(1), 2–11 (2011). https://
doi.org/10.1680/gein.2011.18.1.2
Hamidi, A., Hooresfand, M.: Effect of fiber reinforcement on triaxial shear behavior of cement
treated sand. Geotext. Geomembr. 36, 1–9 (2013). https://doi.org/10.1016/j.geotexmem.2012.
10.005
Li, C., Zornberg, J.G.: Mobilization of reinforcement force in fiber-reinforced soil. J. Geotech.
Geoenviron. Eng. ASCE 139(1), 107–115 (2013). https://doi.org/10.1061/(ASCE)GT.1943-
5606.0000745
Li, J., Tang, C., Wang, D., Pei, X., Shi, B.: Effect of discrete fiber reinforcement on soil tensile
strength. J. Rock. Mech. Geotech. Eng. 6(2), 133–137 (2014). https://doi.org/10.1016/j.jrmge.
2014.01.003
Lirer, S., Flora, A., Consoli, N.C.: On the strength of fiber-reinforced soils. Soils Found. 51(4),
601–609 (2011). https://doi.org/10.3208/sandf.51.601
Michalowski, R.L., Cermak, J.: Triaxial compression of sand reinforced with fibers. J. Geotech.
Geoenviron. Eng. ASCE 129(2), 125–136 (2003). https://doi.org/10.1061/(ASCE)1090-0241
(2003)129:2(125)
Pino, L.F.M., Baudet, B.A.: The effect of particle size distribution on the mechanics of fiber-
reinforced sand under one-dimensional compression. Geotext. Geomembr. 43(3), 250–258
(2015). https://doi.org/10.1016/j.geotexmem.2015.02.004
Tang, C.S., Shi, B., Zhao, L.Z.: Interfacial shear strength of fiber-reinforced soil. Geotext.
Geomembr. 28(1), 54–62 (2010). https://doi.org/10.1016/j.geotexmem.2009.10.001
Tingle, S., Webster, S.L. Santoni, R.L. Discrete fiber reinforcement of sands for expedient road
construction. Technical Report GL-99-3, U.S. Army Corps of Engineers Waterways
Experiment Station (1999)
Yetimoglu, T., Salbas, O.: A study on shear strength of sands reinforced with randomly
distributed discrete fibers. Geotext. Geomembr. 21(2), 103–110 (2002). https://doi.org/10.
1016/S0266-1144(03)00003-7
Comparative Analysis on the Effect of Asphalt
Film Aging Test and Actual Production,
Transportation and Paving on Asphalt Aging
Degree
1 Introduction
Asphalt aging can reduce adhesion, water stability performance, and thermal crack
resistance of asphalt mixture [1, 2]. During mixing, transporting, and placing, hot mix
asphalt mixtures need to be covered and kept at high temperature. Maintaining at high
temperature ensures the good workability of hot mix asphalt mixture. Meanwhile, it
also causes the complex physical and chemical changes of asphalt mixture that leads to
aging [3].
Usually, thin film oven test (TFOT) and rotation thin film oven test (RTFOT) were
employed to simulate asphalt aging of construction. However, there are not many
studies devoted to understanding asphalt aging in asphalt mixtures in China. Based on
ASTM D2127 [4] and D5404 [5], the paper investigated asphalt aging under different
construction condition by asphalt separation centrifuge extractor and rotary evaporator.
In addition, this paper discussed the relationship between lab aging and actual onsite
aging. Methodologies to reduce asphalt aging and improve asphalt pavement durability
are proposed in the paper.
2 Test Method
To extract aged asphalt from mixture, asphalt separation centrifuge extractor and rotary
evaporator were utilized in this study. To investigate the effects of test instruments on
aged asphalt test results, basic engineering indexes of origin asphalt, extracted asphalt
from asphalt trichloroethylene solution, and extracted asphalt from asphalt mortar
which was mixture of origin asphalt and mineral filler have been measured and
comparisons were made.
3 Test Instrument
Digital asphalt penetration tester (Type WSY-026), automatic asphalt soft point tester
(Type SYD-2806-E), low temperature ductility tester (Type LYY-9B), RTFOT tester
(Type SYD-0610), asphalt separation centrifuge extractor (Type 20-1100), and rotary
evaporator (Type BUCHI) were employed to measure engineering indexes.
The origin asphalt used in the study is 70# (Penetration 60 * 80) asphalt. The asphalt
solution was prepared by trichloroethylene to dissolve the origin asphalt, and the
asphalt was extracted by the rotary evaporator. The mineral filler for origin asphalt
weight ratio to asphaltic mortar was 1.6. The dissolved mineral filler in the solution was
extracted by asphalt separation centrifuge extractor and rotary evaporator. Test results
were shown in Table 1.
The proportion of extracted filler was up to 99.7% from Table 1, which indicated
that there was a little mineral filler missed during separation centrifuge and rotary
evaporation process.
The basic engineering indexes [6] of asphalt were tested after three treatments, and
the results were shown in Table 2.
The results in Table 2 showed that there were some differences in engineering
indexes of origin asphalt and the asphalt extracted from trichloroethylene asphalt
solution: the penetration of extracted asphalt decreased 0.05 mm, soft point decreased
0.6 °C, and no change in ductility. The asphalt extracted from asphalt mortar was harder
than origin asphalt: its penetration increased 0.09 mm, soft point increased 0.3 °C,
164 Q. Zhang and Z. Pu
and ductility kept stable. Minor variability in the numerical value met the standard
permissible error requirement of repetition, reproducibility, which meant asphalt sepa-
ration centrifuge extractor and rotary evaporator had no effect on the three basic engi-
neering indexes. The studies of Fujie Zhou [7] and Hongsheng Li [8] also confirmed
asphalt separation centrifuge extractor and rotary evaporator effectiveness for study on
extracted asphalt.
5 Test Results
Ten different 70# origin asphalt were included in the study, and their engineering
indexes were shown in Table 3. RTFOT test was used to simulate short time aging
[3, 6], and results were shown in Table 4.
Meanwhile, the mixing plant selected these origin asphalts to produce asphalt
mixtures with the same aggregate gradation. After placing, aged asphalts have been
extracted from mixtures by asphalt separation centrifuge extractor and rotary evapo-
rator. The corresponding engineering indexes are shown in Table 5.
It view of Tables 3 and 5, it is found that extracted asphalt from mixture were quite
different in engineering indexes compared with origin asphalt as penetration decreased
significantly, soft point increased drastically, and ductility decreased significantly.
These showed that extracted asphalt from mixture became stiffer and more brittle that
lead to lower thermal-crack resistance. It indicates that extracted asphalts show signs of
significant aging.
By comparing Tables 4 and 5, it was found that 10 asphalt samples extracted from
mixture had smaller penetration, higher soft point temperature, and lower ductility as
compared with sample after RTFOT. It is concluded that mixture samples after RTFOT
show signs of significant aging.
After RTFOT, penetration and soft point of asphalt samples extracted from mix-
tures were compared with those of origin asphalt samples, as shown in Table 6.
In terms of retained penetration ratio in Table 6, the minimum value among
samples after RTFOT was 62%, but only two asphalt samples were extracted from
mixtures that had higher retained penetration value of more than 62%. The minimum
value among samples extracted from mixture was 23%. Most of samples became more
consist and stiffer. For softening point, it was found that the maximum and minimum
increased values among samples after RTFOT were 17, and 8% respectively. Note that
166 Q. Zhang and Z. Pu
Table 6. Comparisons of asphalt aging between lab simulation and field samples
Serial Retained penetration ratio (%) Soft point increment (%)
number Sample after Sample extracted Sample after Sample extracted
RTFOT/origin from RTFOT/origin from
sample mixture/origin sample mixture/origin
sample sample
1# 68 49 11 41
2# 62 23 11 61
3# 62 57 8 17
4# 65 49 17 19
5# 63 50 15 24
6# 62 53 10 14
7# 62 62 15 17
8# 62 26 16 51
9# 64 54 13 20
10# 68 66 11 15
Average 63.8 48.9 12.7 27.9
maximum and minimum increases among samples extracted from mixture were 61, and
14%, respectively. Softening point of samples extracted from mixtures had obvious
increased, which mean reducing anti-aging ability, and more negative impact on
engineering projects.
In this study, the AC-20 mixture is made of skeleton suspended hot mix asphalt and
limestone was used as gravel. Before mixing, 70# asphalt was usually heated to
155 * 165 °C, and the temperature of mineral aggregate was higher than asphalt
temperature 10 * 30 °C for drying and heating. During mixing, asphalt contacted
with aggregate particles in thin film state and aged most rapidly. The mixing time was
limited in 3 min, but it played important role on aging of asphalt [10]. The average
temperature of asphalt mixture is about 160 °C, and the temperature coefficient of
variation is below 1.6%, so we think the discharge temperature of the mixture is
uniform and stable. The time hot mix asphalt mixture produced by mixing factory, then
shipped to the construction site is shown in Table 7. The spread speed of asphalt
mixture is controlled at 2 m/min to ensure uniform spread. The average temperature of
paved asphalt mixture is 145 °C. Compaction is divided into three stages which are
initial compaction, repeat compaction and final compaction. As can be seen from
Table 7, the transportation time of the Nos. 2, 8 and 9 samples exceeds 2 h, and the
deterioration shown in Table 5 is more serious than that of other samples.
In the process of storage, transporting and placing of hot mix asphalt, the asphalt
mix is kept in the high temperature and thus asphalt aging was continuous. Although
the rotary evaporator method would cause a little degree of aging during the process of
asphalt extraction, the influence is relatively limited as comparing to the origin asphalt
[11–13]. So most aging of asphalt extracted from field samples occurred in the stages of
mixing, transporting, and placing.
6 Conclusions
It was found that there was low influence on basic engineering index properties of
extracted asphalt from mixture by using asphalt separation centrifuge extractor and
rotary evaporator.
Through the analyses of the test results, it is concluded that the degrees of aging are
much higher for asphalt samples from the field than those samples from lab aging
simulation. Because of imitated field samples, it could not represent all field conditions.
Thus, it is suggested that additional studies for asphaltic aging and their associated
components should be investigated thoroughly.
It was recommended to strictly control the heating time of asphalt and to minimize
any delays on mixing, transporting and placing. It is because the higher the heating
temperature of asphalt and mineral aggregate, the greater impacts on the aging of
asphalt. Thus, it is advisable to strictly control the mixing temperature to a certain
temperature interval. In addition, heating and drying temperature of aggregate should
be reduced as much as possible and the mixing time should be limited. It further
suggests of the benefits of utilizing warm mix techniques.
References
1. Ma, Q., Rong, G., Juan, Z., et al.: Analysis on the influence of asphalt aging on the
properties of asphalt and asphalt mixture [J]. Highway traffic technology 3, 76–78 (2012)
2. Lu J.: Study on aging behavior and recycling technique of asphalt pavement [D]. Chang An
University (2008)
3. Shi, Z., Yin, X., Zhang, W., et al.: Effect of asphalt aging methods on asphalt performance
[J]. Build. Mat. World 33(5), 45–47 (2014)
4. ASTM: Standard Practice for Recovery of Asphalt from Solution using the Rotary
Evaporator D5404. West Conshohocken, PA (2011)
5. ASTM: Standard Test Methods for Quantitative Extraction of Bitumen from Bituminous
Paving Mixtures D2127. West Conshohocken, PA (2011)
6. JTG E20-2011.: The test criterion for highway project pitch and pitch compounds [S]
7. Zhou, F.: RAS binder characterization and blending with virgin binders. Annual Meeting of
Transportation Research Board, Washington, DC
8. Li, H.S., Wu, Y.B., Guo, Y.F.: Validation of reclaimed shingles asphalt binder extraction
and recovery methods. In: Geo-Hubei 2014 International Conference on Sustainable
Infrastructure: Advanced Characterization of Asphalt and Concrete Materials, pp. 17–23.
Yichang, Hubei, China, American Society of Civil Engineers (ASCE) (2014)
9. He, H., Chen, H.: Evaluation of two experimental methods for the aging performance of
asphalt [J]. West. Explor. Eng. 1, 249–250 (2006)
168 Q. Zhang and Z. Pu
10. He, Z., Zumin, O., Nan, L., et al.: Effect of heating temperature and heating time on anti-
aging performance of asphalt mixture [J]. J. Chongqing Jiaotong University 31(2), 239–242
(2014)
11. Tian, X., Mo, Y., Zheng, J.: Effect of extracting reclaiming process on aging degree of
asphalt [J]. J. Transp. Eng. 5(2), 23–26 (2005)
12. Zheng, C., Xin, M., Dan, L., et al.: Contrast test of asphalt reclaimed method for aging
degree [J]. Highw. Eng. 32(6), 165–168 (2007)
13. Wang, H., Wang, F.: Discussion on the effects of asphalt aging and pavement performance
[J]. China Water Transp. 08(01), 96–97 (2008)
Response of Swelling Clays to Superstructure
Vertical Loads
1 Introduction
The estimated annual cost related to the problems associated with the expansive
behavior of the soil is billions of dollars globally (Jone and Jones 1987). Furthermore,
designing sustainable structure able to resist this behavior is very costly. Therefore, a
realistic estimation of the swelling pressure attracts the attention of researchers and
practitioners alike in order to obtain realistic figures. In order to assess the uplift
movement of the swelling soil, the swelling pressure shall be addressed. Swelling
pressure can be defined as the pressure required keeping a soil element at a constant
volume (Caterina Di Maio 2001). A proper evaluation of swelling pressure is con-
sidered as a vital parameter in designing sustainable structures. According to
Prabhakara Rao (2013), the swelling of the soils depends on many factors such as type
of clay, clay minerals, moisture content and surcharge weight. At a laboratory scale, the
determination of swelling pressure can be carried out through different methodologies
which are summarized in Table 1.
In Table 1 the constant volume swell pressure test is conducted with water added
frequently to a swelling soil with increments of applied pressure to maintain constant
volume. In the different pressure test samples are subjected to a series of loading and
unloading until reducing the sample to its original volume. The double oedometer
method includes two samples subjected to oedometer tests. One tested for swell under
wetting and the second sample tested at the natural water content. In the swell con-
solidation method both swelling and consolidation tests are carried out together. From
logical point of view, if an element is subjected to lateral confinement then its chance to
expand vertically is much higher than a freely swelling matrix in 3D. Hence, the major
drawback of the conventional oedometer method is the constraint condition of samples
with the ability of upward movement only. In order to simulate the actual field con-
dition more closely where a soil element can swell freely in all directions, the 3D-Swell
apparatus have been introduced by numerous researchers (Aravind Pedarla et al. 2016;
Response of Swelling Clays to Superstructure Vertical Loads 171
Mutaz et al. 2011a, b and Al-Shamrani et al. 2010) to evaluate the vertical swelling as
well as the volumetric strains. The 3D-Swell setup is comprised of 70 mm compaction
mold with 140 mm height. A vertical dial gauge was placed at the top and used to
determine vertical elongation of the specimen at select time periods and a pie tape was
used to measure lateral expansion of the specimen at same time periods (Al-Shamrani
et al. 2010). Once a sample is prepared under certain water content and density con-
ditions, then it will be covered by a membrane and placed inside the bottom of the mold
while porous stones were placed at the top and bottom of the tested sample. Before
running the test, the vertical dial gauge reading as well as the initial diameter of the
sample were recorded and noted as reference readings. Then, water is introduced to the
sample from top and bottom with vertical and volumetric expansion monitored.
This study is aimed at providing a guide for the practice and design engineers to
help them adjust the figures quoted by laboratories or literature (Dafalla and
Al-Shamrani 2012) and amends them as appropriate to reflect realistic field conditions.
A research study was conducted on samples from Al-Ghatt province in Saudi Arabia.
A series of 3D-Swell tests were carried out at untreated and treated samples with
different percentages of lime and cement at different conditions of moisture content for
the sake of measuring the vertical and volumetric strains and the ratio between these
strains. In addition, oedometer tests were conducted in order to measure the vertical
swelling as well as swelling pressure at lateral restraint condition. A comparison was
done between the 3D-Swell and oedometer test results and a realistic design criterion
was presented.
2 Experimental Program
Materials for this study were collected from Al-Ghatt province located 270 km
Northwest of Riyadh in Saudi Arabia. The severity nature of Al-Ghatt soil has been
addressed by several researchers (Dafalla and Al-Shamrani 2008; Mutaz et al. 2011a, b).
Referring to Al-Shamrani et al. (2010), Al-Ghatt soil has been characterized as high
expansive clay with presence of expansive minerals pertaining to smectite group. Basic
soil tests were carried out for Al-Ghatt clay such as; gradation, hydrometer, Atterberg
limits, specific gravity and compaction test. According to the gradation test, the per-
centage of fine material is 86%. From Atterberg limits, the liquid limit, plastic limit and
plasticity index are 60, 30 and 30 respectively. The specific gravity value is 2.8. Al-
Ghatt clayey soil is classified as CH in accordance with the unified soil classification
system (USCS). Compaction tests were conducted in order to determine the optimum
moisture content and its relevant maximum dry density. In order to study the rela-
tionship between vertical and volumetric strain for virgin and treated clays, Al-Ghat
clay was stabilized with different percentage of lime (4 and 8%) and cement (3 and 6%)
at different conditions of moisture; namely OMC (optimum moisture content), DOMC
(dry of optimum moisture content) and WOMC (wet of optimum moisture content).
Compaction curves of untreated and treated clays are depicted in Fig. 1.
172 M. A. Dafalla et al.
15.0
Fig. 2. Influence of lime addition on the vertical swell in an oedometer test for Al-Ghatt clay
untreated clay is reduced by adding lime. Figure 2 shows a graph of typical reduction
in vertical swell of Al-Ghatt clay by 4% lime addition. The swelling pressure required
to bring the vertical swell to the original height is assumed at three stress levels. Stress
level 1 is for untreated Al-Ghatt clay and stress level 2 is for Al-Ghatt clay treated with
4% lime and stress level 3 is for Al-Ghatt clay treated with 8% lime. It can be
demonstrated that stress level 1 is greater than stress level 2 and stress level 2 is greater
than stress level 3. No quantitative figures are given in this stage of the research study.
Figure 3 presents the axial/volumetric strain ratio as obtained in a 3D free swell test for
the three stress levels defined for Al-Ghatt clay.
0.6
Axial volumetric strain
0.5
0.4 optimum moisture
state
ratio
0.3
0.2 dry of optimum
0.1
0 wet of optimum
0 1 2 3 4
Stress level
Fig. 3. Presents the axial/volumetric strain ratio as obtained in a 3D free swell test for the three
stress levels defined for Al-Ghatt clay
174 M. A. Dafalla et al.
The percentages of vertical swell obtained from the triaxial swell tests are found to be
considerably lower than the corresponding values measured in the oedometer tests.
Al-Shamrani and Al-Mhaidib (2000). When comparing the swelling pressure obtained
in an oedometer using one dimensional ASTM method of expanding and loading, none
realistic swell pressure is obtained due to the lateral confinement. Practicing engineers
often use a factor of one third to estimate actual value. However, this study is aimed at
refining this factor by suggesting 3D swell simple test. The authors see that the ratio of
vertical strain to volumetric strain can be utilized to give better estimation.
The results of vertical and volumetric strains for all samples as well as the ratio
between both strains are presented under Table 2. It is shown from Table 2 that, the
effectiveness of treatment is simultaneously increased with the increment of stabilizer
dose. Moreover, significant reduction in both vertical and volumetric strains was
observed with the treatment for DOMC and WOMC conditions. Therefore, stabilizer
dose and compaction condition play a vital role in the effectiveness of the treatment.
This is similarly valid for cement treated clays but only lime was used to demonstrate
the suggested approach for this procedure. The data for clay stabilized with cement
indicates a similar trend of that shown for lime stabilized clay. Al-Mhaidib (1998)
suggested that triaxial tests provides lower values of vertical swell than swell measured
from one dimensional oedometer tests. His work indicated that for all swelling pressure
values, the ratio of triaxial swell to the oedometer swell is constant and estimated at
0.333. This is not likely to be the case when the mineralogy of the clay or clay content
is variable. The current study suggests that the axial to volumetric strain ratio decreases
with decrease of the swell pressure. Figure 2 does not suggest quantitative values but
the decrease with decrease is clearly shown. The constant volume method and the one
dimensional oedometer method were compared for identical sample and results
obtained showed that the oedometer method gives higher swelling pressure values than
the constant volume method, El Fatih and Muawia (1984).
The test results suggest that the ratio of axial to volumetric strain can be used in
estimation of swelling pressure provided that a correction factor is applied. The authors
of this work suggest considering a correction factor to predict the swelling pressure.
Assuming a linear elastic behaviour, swell deformation measured can be used to correct
swelling pressure values using a correction factor coupled with the ratio of axial to
volumetric strain as obtained from a simple 3D swell test. This correction factor can be
of local validity and can be established for specific region or clay type. This shall be
taken as a rough guide in preliminary design works and must be confirmed by detailed
testing.
Response of Swelling Clays to Superstructure Vertical Loads 175
Fig. 4. Relationship between swell pressure and free swell. After Kaybali and Demir (2011)
4 Conclusion
The ratio of axial to volumetric strain obtained in a 3D swell test can be useful in
predicting actual swelling pressure when only one dimensional oedometer testing is
performed. Expanded and loaded test method tends to overestimate the swelling
pressure due to the confinement of sample during the test. It is suggested to use the 3D
swell test as an additional tool to better predict the swelling pressure obtained from one
dimensional oedometer test. This study suggests an equation for predicting the actual
swelling pressure. The actual swelling pressure can be obtained by multiplying the
oedometer swelling pressure by the 3D swell test ratio of axial to volumetric strain and
a regional correction factor. Regional and area correction factor can be obtained from
testing a series of samples.
References
Al-Mhaidib, A.: Prediction of swelling potential of an expansive shale. In: Proceedings of the
Second International Conference on Unsaturated Soils, 27–30 Aug 1998, Beijing, China
(1998)
Al-Shamrani, M., Al-Mhaidib, A.: Swelling behavior under oedometric and triaxial loading
conditions. Geo-Denver 2000, 5–8 Aug 2000, Denver, Colorado, United States (2000)
Caterina Di Maio: Swelling pressure of clayey soils: the influence of stress state and pore liquid
composition. Rivista Italiana di Geotecnica 3/2001 (2001)
Dafalla, M.A., Al-Shamrani, M.A.: Expansive soil properties in a semiarid region. Res.
J. Environ. Earth Sci. 4, 930–938 (2012)
Dafalla, M.A., Al-Shamrani, M.A.: Performance-based solutions for foundations on expansive
soils-Al-Ghatt region, Saudi Arabia. In: Geo-Chiangmai 2008, Chiangmai, Thailand (2008)
El Fatih, M.A., Muawia, A.D.E.: Comparison of two methods for the measurement of swelling
pressure. In: Proceedings of the Fifth International Conference on Expansive Soils, Adelaide,
South Australia, 21–23 May 1984
Response of Swelling Clays to Superstructure Vertical Loads 177
Jone, D.E., Jones, K.A.: Treating expansive soils. Civil Eng. ASCE 43, 11 (1987)
Kayabali, K., Demir, S.: Measurement of swelling pressure: direct method versus indirect
methods. Can. Geotech. J. 48(3), 354–364 (2011)
Mutaz, E., Dafalla, M.A., Shamrani, M.A.: Stabilization of Al-Ghatt clay shale by using a
mixture of lime and cement. In: Geo-Hunan Conference, China, American Society of Civil
Engineers—Geotechnical Special Publication GSP 223 (2011a)
Mutaz, E., Shamrani, M.A, Puppala, A.J., Dafalla, M.: Evaluation of chemical stabilization of a
highly expansive clayey soil. Transp. Res. Board J. (Washington DC, USA, Low-Volume
Roads 2011) 2, 148–157 (2011b)
Pedarla, A., Puppala, A.J., Laureano, R.H.: Evaluation of swell behavior of expansive clays from
internal specific surface and pore size distribution. J. Geotech. Geoenviron. Eng. 142(2)
(2016)
Prabhakara Rao, V.V.N.: Swelling pressure of soil using a predictive tool. Int. J. Comput. Appl.
(0975–8887) 63(19) (2013)
Shamrani, M.A, Mutaz, E., Puppala, A.J., Dafalla, M.A.: Characterization of problematic
expansive soils from mineralogical and swell characterization studies. American Society of
Civil Engineers—Geotechnical Special Publication GSP 199, Geoflorida, Florida, No. 199,
pp. 793–802 (2010)
Author Index
A J
Aldar, Dilip, 22 Jiang, Peng, 29
Al-Shamrani, M.A., 169 Ji, Shunying, 98
Assaf, Gabriel J., 11 Jin, Qing, 136
Atrash, Khlifa El, 11
K
B Komba, Julius J., 39
Bai, Jiwen, 29 Kong, Qiaoli, 126
Kou, Yan, 52
C
Cheng, Bingchuan, 29 L
Chen, Qizhi, 71 Liang, Luju, 71
Cui, Xinzhuang, 136 Li, Shucai, 29
Liu, Ming-Hung, 115
D Liu, Rentai, 29
Dafalla, M.A., 169 Lu, Tu, 136
Derbin, Y.G., 80
M
E Maina, James W., 39
Evans, Robert, 60 Malisa, John T., 39
Marshall, A. M., 80
F Mataka, Mussa, 39
Fan, Xiaozhen, 71 Mohyeddin, Alireza, 52
Feng, Xiao, 29 Mo, Pin-Qiang, 1
Mutaz, E., 169
G
Gao, Feng, 1 P
Georgees, Romel, 60 Patel, Suchit Kumar, 153
Guo, Jinyun, 126 Pu, Zhichao, 162
H Q
Han, Litao, 126 Qiu, Peng, 98
Hassan, Rayya, 60 Quah, Tai-Seong, 115
Hsu, Sung-Chi, 115
Huang, Yishuo, 115
S Wanatowski, D., 80
Shen, Yi, 126 Wang, Zhongxiao, 136
Shukla, Sanjay Kumar, 52
Singh, Baleshwar, 153 X
Su, Junwei, 136 Xu, Changjie, 71
T Z
Tapase, Anand, 22 Zhang, Jin, 98
Tapase, Rajashree, 22 Zhang, Lei, 136
Zhang, Qingqing, 162
W Zhao, Honghua, 98
Walker, J., 80 Zhou, Guoqing, 1
Walubita, Lubinda F., 39