(Limin Zhang, Yu Wang, Gang Wang, Li Dianqing) Geo
(Limin Zhang, Yu Wang, Gang Wang, Li Dianqing) Geo
(Limin Zhang, Yu Wang, Gang Wang, Li Dianqing) Geo
Zhang
Wang
Wang
Li
Geotechnical Safety and Risk IV contains the contributions presented at the
4th International Symposium on Geotechnical Safety and Risk (4th ISGSR,
Hong Kong, 4-6 December 2013), which was organised under the auspices of
Geotechnical Safety and Risk IV
the Geotechnical Safety Network (GEOSNet), TC304 on Engineering Practice
of Risk Assessment and Management and TC205 on Safety and Serviceability
an informa business
GEOTECHNICAL SAFETY AND RISK IV
Editors
L.M. Zhang
The Hong Kong University of Science and Technology, HKSAR
Y. Wang
City University of Hong Kong, HKSAR
G. Wang
The Hong Kong University of Science and Technology, HKSAR
D.Q. Li
Wuhan University, China
All rights reserved. No part of this publication or the information contained herein may be reproduced,
stored in a retrieval system, or transmitted in any form or by any means, electronic, mechanical, by
photocopying, recording or otherwise, without written prior permission from the publisher.
Although all care is taken to ensure integrity and the quality of this publication and the information
herein, no responsibility is assumed by the publishers nor the author for any damage to the property or
persons as a result of operation or use of this publication and/or the information contained herein.
Table of contents
Preface xi
Organisation xiii
Acknowledgements xv
2 Keynote lectures
Advances in geotechnical risk and reliability for offshore applications 29
R.B. Gilbert, S. Lacasse & F. Nadim
Homogenization of geomaterials using the random finite element method 43
D.V. Griffiths, J. Paiboon, J. Huang & G.A. Fenton
Selecting optimal probability models for geotechnical reliability analysis 53
H.W. Huang, W.W. Su & J. Zhang
Robust design of geotechnical systems—a new design perspective 69
C.H. Juang, L. Wang & S. Atamturktur
Integrated geo risk management: Crossing boundaries 79
M.Th. van Staveren
Is landslide risk quantifiable and manageable? 101
H.N. Wong
5 Geohazards
Integrating seismic hazard analyses with geotechnical site characterization for liquefaction
potential assessment in Kaohsiung area 201
D. Huang & J.P. Wang
Landslide considerations for low cost remedial works for the Karakorum Highway, Pakistan 207
A.D. Mackay
Landslide mitigation considerations for low cost highway remedial works, East Timor 215
A.D. Mackay
The landslide of Kirf: A chain of governance failures 221
S. Van Baars, M. Sosson, S. Jung & R. Becker
A risk-based active fault classification 225
J.P. Wang
Introducing non-stationary earthquake process concept: Including an analytical model and
a case study in Central Taiwan 229
Y. Xu & J.P. Wang
Modelling techniques of submarine landslide in centrifuge 235
C. Zhao, W.Z. Zhang, J.H. Zhang, K.Z. Wang & Z.S. Xiong
vi
vii
viii
ix
Preface
The 4th International Symposium on Geotechnical Safety and Risk (4th ISGSR) was organised by the
Hong Kong University of Science and Technology under the auspices of the Geotechnical Safety Net-
work (GEOSNet; Chair, Daniel Straub; Co-chair Limin Zhang), Technical Committee TC304 on Engi-
neering Practice of Risk Assessment and Management (Chair, K.K. Phoon) and Technical Committee
TC205 on Safety and Serviceability in Geotechnical Design (Chair, Brian Simpson) of the International
Society of Soil Mechanics and Geotechnical Engineering (ISSMGE). The Symposium was also sup-
ported by Hong Kong Geotechnical Society, the Geotechnical Division of the Hong Kong Institution of
Engineers, Chinese Institution of Soil Mechanics and Geotechnical Engineering, the Engineering Risk
and Insurance Branch of China Civil Engineering Society, and American Society of Civil Engineers—
Hong Kong Section.
The 4th ISGSR was a continuation of a series of symposiums and workshops on geotechnical risk and
reliability starting with LSD2000 in Melbourne, Australia, IWS2002 in Tokyo and Kamakura, Japan,
LSD2003 in Cambridge, USA, Georisk2004 in Bangalore, India, Taipei2006 in Taipei, 1st ISGSR in
Shanghai, China in 2007, 2nd ISGSR in Gifu, Japan in 2009 and 3rd ISGSR in Munich, Germany in
2011.
Safety, reliability, and risk assessment and management have attracted growing interests of the geotech-
nical community in recent years due to the frequent occurrences of natural and man-made disasters and
the needs for safe and cost-effective design, construction and operations of infrastructures. At the same
time there is an increasing expectation of the general public that requires the engineering community to
provide quantitative information concerning risks posed by geotechnical hazards. The 4th ISGSR pro-
vided an excellent opportunity to better understand the geotechnical safety and risk management issues
in engineering practices and research. The proceedings contain seven invited keynotes and 69 accepted
papers from 28 countries and regions. Each accepted paper in the conference proceedings was subject to
review by two peers. These papers cover six themes: (1) geotechnical uncertainty and variability, (2) geo-
hazards such as landslides, earthquakes and climate changes, (3) reliability and risk analysis, (4) reliability-
based design and limit-state design in geotechnical engineering, (5) risk assessment and management in
geotechnical engineering and infrastructural projects, and (6) practical applications.
One of the highlights of this symposium was the 3rd Wilson Tang Lecture. The lecture was inaugurated
during the 2nd ISGSR in Gifu to recognize and honor the significant contributions of the late Professor
Wilson Tang, who was one of the founding researchers in geotechnical reliability and risk. The first lec-
ture was given by Prof. T. H. Wu of the Ohio State University and the second lecture by Prof. Y. Honjo of
Gifu University. The 3rd lecture was given by Prof. Suzanne Lacasse of Norwegian Geotechnical Institute
during the 4th ISGSR.
The credit for the proceedings goes to the authors and reviewers. The publication of the proceedings
was financially supported by the National Basic Research Program of China (Grant No. 2011CB013500)
and the National Natural Science Foundation of China’s Oversea Collaborative Research Program (Grant
No. 51129902).
Limin Zhang
Chairman of the Regional Organising Committee
August 2013, The Hong Kong University of Science and Technology, HKSAR
xi
Organisation
Members
Ivan S.K. Au Jinhui Li
K.T. Chau Jack Pappin
Tony C.T. Cheung L.G. Tham
John Endicott Jun Yang
Albert Ho Kelvin K.V. Yuen
ISSMGE-TC304 Members
K.K. Phoon (Chair) B.K. Low
Gordon Fenton (Vice Chair) Dagang Lu
Jianye Ching (Secretary) Farrokh Nadim
Y. Ashkey Laszlo Nagy
Gregory Baecher Shin-ichi Nishimura
Hansgeorg Balthaus Ali Noorzad
Robert Berkelaar Lars Olsson
Laura Caldeira Lance Roberts
Paul Cools Adrian Rodriguez-Marek
Federica Cotecchia Nick Sartain
Guillermo Espin Bernd Schuppener
P. Ganne Timo Schweckendiek
Anthony TC Goh V.I. Sheinin
Vaughan Griffiths Kenichi Soga
Ken Ho Abdul-Hamid Soubra
Hongwei Huang M. Suzuki
Mark Jaksa Pavlos Tyrologou
Rafael Jiménez Martin van Staveren
Leena Korkiala-Tanttu Yu Wang
Kishor Kumar Jianfeng Xue
Giuseppe Lanzo Limin Zhang
Dianqing Li
xiii
xiv
Acknowledgements
MANUSCRIPT REVIEWERS
The editors are grateful to the following people who helped to review the manuscripts and hence assisted
in improving the overall technical standard and presentation of the papers in these proceedings:
xv
ABSTRACT: The paper is in homage to Professor Wilson Tang for his inspiration to fellow engineers
in the area of geotechnical engineering. The role of statistics, probability and reliability in geotechni-
cal engineering is first outlined. Examples of solutions based on Wilson Tang’s pioneering work are
presented: uncertainties in soil parameters; Bayesian updating applications; reliability of tailings dam;
model uncertainty and calibration of safety factor. Two aspects of special interest to Wilson Tang are
also briefly discussed: improving the cost-effectiveness of site investigations and the reliability of offshore
structures.
Table 2. Probability distributions for different soil system. Both statistical variability (i.e. second-
properties (adapted from Lacasse & Nadim 1996). moment) and spatial variability (i.e. spatial cor-
relation) of soil properties affect the reliability of
Soil property Soil type PDF geotechnical systems.
Cone resistance Sand, clay N/LN The number of studies making use of random
Undrained shear Clay (triaxial tests) LN field simulation, finite elements and Monte Carlo
strength Clayey silt N simulation is still limited. The importance of the
Normalized undrained Clay N/LN results so far, however, should be a stimulus for the
shear strength transposition of results to practice.
Plastic limit Clay N Popescu et al. (2005) investigated the differen-
Submerged unit weight Clay, silt, sand N tial settlements and bearing capacity of a rigid
Friction angle Sand N strip foundation on an overconsolidated clay layer.
Void ratio, porosity Clay, silt, sand N The undrained strength of the clay was modelled
Overconsolidation ratio Clay N/LN as a non-normal random field. The deformation
modulus was assumed to be perfectly correlated
to undrained shear strength. The settlements (uni-
– There generally exists a critical correlation
form and differential settlements) were computed
distance which corresponds to a minimum
with non-linear finite elements in a Monte Carlo
reliability.
simulation framework. Anisotropy in spatial cor-
– Phenomena governed by highly non-linear
relation was addressed, with the horizontal scale of
behavior laws are affected the most by spatial
fluctuation exceeding the vertical scale of fluctua-
variations.
tion by a factor of 10.
Variance reduction alone cannot convey a Figure 3a shows the contours of maximum
comprehensive picture of the implications of spa- shear strain for a uniform soil deposit with und-
tial variability on the behavior of a geotechnical rained strength of 100 kPa and for a normalized
Figure 3. Results of investigation on homogeneous and spatially random foundation soil (Popescu et al. 2005).
Figure 4. Left: Single footing and two footings founded on a spatially heterogeneous soil; Right: 3D finite element
mesh of spatially heterogeneous soil volume supporting two footings (Fenton & Griffiths 2005).
location and shape of the failure surface is strongly (DSS) tests were used (Fig. 6). No outliers were
related to the presence of weaker soil zones (shown evident from visual inspection. The data from
in lighter colors) and is, in both cases, markedly depths 0 to 5 m were excluded for this analysis. The
different from the shapes assumed in earth pres- following model for the total Coefficient of Varia-
sure theory. tion (COVtot) was used:
Griffiths et al. (2013) brings new developments
on soil variability and random finite element Vto2t
COV COV 2ω + COV
Vm2 COV
Vse2 (1)
analysis.
where COVω is the coefficient of variation of
3.4 Application to Troll clay inherent variability, representative of aleatory
uncertainty; COVm is the coefficient of varia-
Uzielli et al. (2006a) did an uncertainty-based
tion of measurement error; and COVSE is the
geotechnical characterization of the Troll clay,
coefficient of variation of statistical estimation
a site offshore Norway for the world’s largest
uncertainty.
gravity structure. Second-moment statistics were
Figures 7 (CAUC data) and 8 (DSS data) and
obtained from laboratory and in situ tests. Bayesian
Table 3 present the second-moment estimates of
updating combined the values of undrained shear
strength resulting from triaxial compression tests
and piezocone tests. Some of the results are pre-
sented herein.
The authors believe that the approach followed
would have been close to what Wilson Tang would
have done himself, if he had been asked to interpret
and calculate the uncertainty in the Troll data.
The characterization consisted of four steps:
(a) visual inspection of data by soil unit and pre-
liminary second-moment data analysis; (b) iden-
tification of a deterministic trend function and
decomposition; (c) identification of a suitable
uncertainty model; and (d) quantification of the
uncertainty (mean, variance, standard deviation or
coefficient of variation). The Kendall’s tau statistic
test (Uzielli et al. 2006b) was run to check whether
or not the data were statistically independent.
A trend function was obtained by regression analy-
sis (Ang & Tang 2007). An uncertainty model was
used to merge the different uncertainty compo-
nents to estimate the total uncertainty.
μt (kPa)* 28.4 80.5 25.8 69.2 Figure 10. Statistics of cone resistance from 5 CPTU
COVω 0.09 0.06 0.21 0.11 tests: (a) measured data. (b) coefficients of variation,
COVm 0.20 0.20 0.20 0.20 (c) mean and standard deviation (Uzielli et al. 2006a).
COVSE 0.02 0.01 0.04 0.02
COVtot 0.22 0.21 0.29 0.23
10
11
4.1.1 Bayesian network for earthquake-triggered Table 5. Nodes and their possible states in the Bayesian
landslide risk assessment network in Figure 2 (after Nadim & Liu 2013a).
Nadim & Liu (2013a) provided a brief review of # of
Bayesian networks. Figure 13 presents graphically Node states Possible discrete states
a simple Bayesian network with five nodes and five
arcs. The nodes are: Magnitude (M), Distance (D), Magnitude (Mw) 6 4.0–4.5–5.0–5.5–
Seismic severity (S), Landslide severity (L), and 6.0–6.5–7.0
Building damage (B). These nodes are connected Distance (km) 6 22–25–28–31–34–37–40
PGA (g) 6 0–0.08–0.16–0.24–
0.32–0.40–0.48
Landslide 2 Happens; does not happen
Building 3 No damage; some damage;
damage collapse
Alarm 2 On; off
Measure 2 Yes; no
Decision 4 Passive; active; no action;
warning on
Cost measure, ---
Figure 13. Simple Bayesian network (Nadim & Liu cost, utilities
2013a).
12
Figure 15. Discrete probabilities of distance to the where Dn is the Newmark displacement (cm).
seismic source (Nadim & Liu 2013a). The calculated probabilities of slope failure for
different ranges of PGA are listed in Table 7. As
mentioned above, countermeasures made to land-
slide can reduce risk. The probability of slope fail-
ure when active actions were used are also listed
Table 7.
For a building subjected to a multi-hazard situa-
tion involving additive load effects (e.g. earthquake
and landslide), the damage was increased. For the
PGA
(10−2 g) 0–8 8–16 16–24 24–32 32–40 40–48
13
Building No damage 0.4 0.1 0.52 0.1 Figure 19. Sensitivity analysis of the risk as a function
damage Some damage 0.3 0.1 0.43 0.1 of the probability of slope failure for different mitigation
actions—horizontal arrows indicate range where type of
Collapse 0.3 0.8 0.05 0.8
mitigation measure is the optimum (after Nadim & Liu
2013a).
14
15
Table 10. Results of Bayesian updating of risk associated with talus landslide along
the Shuifu-Maliuwan highway.
16
Figure 25. Cross-section of tailings dam in Romania 5.2 Event tree analysis
(Corser, P. 2009. Personal comm. MWH Americas Inc.
Bucharest, Romania). To establish whether the dam provides acceptable
safety against “uncontrolled” release of tailings and
To establish whether the dam provides acceptable water during its life, an event tree analysis was done.
safety against “uncontrolled” release of tailings A workshop was organized to develop the event trees
and water during its life, an event tree approach was and reach a consensus when quantifying the hazards.
used to do the hazard analyses. This technique iden- The analysis involved breaking down the complex
tified potential failure mechanisms and followed system into its fundamental components, and deter-
how a series of events leading to non-performance mining the potential “failure” mechanisms leading
of a dam might unfold. The probability of each sce- to non-performance of the dam and the physical
nario, given a triggering event, was quantified. processes that could cause such mechanisms.
The event tree hazard analyses considered the The key factors considered in the analyses
dam at different stages of its life and estimated included: dam configuration (Starter Dam, dam
the probability of non-performance. A non- during construction and Completed Dam), and
satisfactory performance of the dam was defined triggers, including earthquake shaking, extreme
as an uncontrolled release of tailings and water rainfall or snowmelt, natural terrain landslide in
from the dam over a period of time. The release the valley or failure of the waste stockpile into the
could be due to a breach of the dam or overtop- tailings reservoir.
ping without breach of the dam. The analyses Acts of war or sabotage, impact by meteorites
looked at critical scenarios, including all potential or other extreme events of this type were not con-
modes of non-performance under extreme triggers sidered, as they would result in so low probabilities
such as a rare, unusually strong earthquake and of non-performance that they are not realistic to
extreme rainfall in a 24-hour period. consider.
The non-performance modes considered
included:
5.1 Design considerations
1. Foundation failure, due to, e.g. excess pore
The most significant requirements that influenced pressures or weak layer in foundation leading to
the probabilities in the hazard analyses include: cracking, instability and breach of the dam.
– Operational freeboard at all times of one meter 2. Dam slope instability downstream or upstream,
above storage level for maximum reclaim pond due to e.g., construction pore pressure in core
and 2 PMP (probable maximum precipitation); of Starter Dam, excessive pore pressures caused
the requirement leads to a storage volume capac- by static or earthquake loads or instability due
ity of two 1/10,000-yr rainfall within the same to inertia forces.
24 hours. 3. Unravelling of downstream toe and slope, due
– Gentle slopes for the Starter Dam (≈2H:1V to e.g. overtopping or excessive leakage through
upstream and ≈2H:1V downstream). or under the dam. This can be caused by a slide
– Gentle downstream slopes for the Completed into the reservoir, dam crest settlement due to
Dam (3H:1V). deformations of the Starter Dam, piping, inter-
– Good quality rockfill for the Starter Dam con- nal erosion and sinkhole formation, or exces-
struction and the Completed Dam. sive deformations (slumping) of the top vertical
– “Well drained” tailings beach at the upstream part of the Completed Dam during earthquake
face of the dam, where equipment can move shaking.
in for repairs, in case of movement or partial 4. Dam abutment failure followed by breach, due to
breach. e.g. slide close to and/or under part of the dam.
– Secondary Containment Dam (SCD) with about 5. Liquefaction of the tailings.
50,000 m3 containment capacity after 16 years. Figure 26 presents some of the configurations
– Diversion channels along the sides of the valley to and examples of the non-performance modes
divert excess rainfall runoff away from the TMF analyzed. Overtopping without breach of the
pond to minimize the risk of overtopping. dam, including under-capacity or damage of the
17
Secondary Containment Dam was also considered, Table 11. Total probabilities of non-performance.
not as a separate non-performance, but as one of
the events in the sequence of events in the trees. P [non-
Configuration performance]
Different conditions can affect the probabil-
ity of a hazard occurring or severity of a conse- Starter dam (t = 1.5 yr, internal erosion) 1.3 × 10−6 /yr
quence, for example construction deficiencies or Completed dam (t = 16 yrs) 1.3 × 10−6/yr
inadequate response of the field control team at Intermediate stage (t = 4 yrs) 6.5 × 10−7/yr
the site when warning signals may appear. The
Intermediate stage (t = 9–12 yrs) 1.3 × 10−6/yr
analyses also looked into construction deficien-
cies, e.g. inadequate filters leading to uncontrolled
internal erosion, inadequate drainage, very weak
construction layers or zones in the embankment, considered successively. The total probability of
inadequate types of material(s) in the embankment non-performance is the sum of all contributing
fill, or insufficient quality control and unforeseen probabilities to the non-performance for each
construction schedule changes. These conditions of the dam configurations. Table 11 presents the
were also integrated in the event trees as separate total probabilities for each configuration of the
events during the course of the construction of the dam (all triggers included). The probabilities were
Starter Dam and Completed Dam. presented as a function of the release of tailings
and water associated with the non-performance of
the dam. The highest annual probability of non-
5.3 Probability of non-performance
performance was 10−6.
At the event tree workshop, the critical times in the The highest probabilities of non-performance
life of the TMF were defined: during construction were associated with earthquake shaking of the
of the Starter Dam, during the downstream con- completed dam and the static liquefaction of the
struction stages, during the centerline construc- tailings at time 9 to 12 years after the start of
tion of the dam, and/or in the early years after the construction. The non-performance scenarios
Completed Dam is built. A matrix of dam config- would result in some material damage and some
uration versus time was prepared. The modes seen contamination, but only in the vicinity downstream
as most critical and susceptible to lead to the high- of the dam. For the Starter Dam, no reasonable
est probabilities of non-performance were listed. expected scenario lead to a significant release of
As part of the mode screening, the following con- tailings and water because of the limited quantity
siderations were subjected to a consensus decision: of water available and the reserve capacity provided
extreme and critical precipitation (rainfall, flood (2 PMP’s). Internal erosion may cause, with an
and snowmelt), likelihood of failure of the waste annual probability of 10−6, a small escape of tailings
stockpile, critical situations after construction of and water. The escape would cause only modest con-
the dams, and geo-environmental considerations. tamination of the immediate vicinity downstream.
Event trees were developed for each dam con- Essentially all material released could be contained
figuration and trigger, with each non-perform- by the Secondary Containment Dam.
ance mechanism looked at separately. In some The analyses showed (1) no plausible events
cases, two non-performance mechanisms were result in an annual probability of non-performance
18
Figure 27. Annual probability of dam failure by inter- Figure 28. Physical impact of TMF dam breach at
nal erosion for different dams in the USA (Von Thun, Roşia Montanǎ (Corser, P. 2009. Personal comm. MWH
1985; Vick 2002). Americas, Inc. Bucharest, Romania).
19
20
21
Figure 29. Cost reduction with increased number of borings (Lacasse & Nadim 1998 based on Tang 1987).
22
23
24
25
26
R.B. Gilbert
The University of Texas at Austin, Austin, Texas, USA
ABSTRACT: This paper describes recent advances in geotechnical reliability and risk for offshore
applications. The topics addressed include spatial variability, model uncertainty, hazard characterization,
reliability-based design, system reliability and risk management. Conclusions from the evolution of reli-
ability and risk approaches include that practical implementation is key, that assessment is best consid-
ered in the context of decision making, and that collaboration of multiple disciplines and stakeholders is
important to managing risk effectively.
29
Consequently, it is not feasible to gather 100 for this model and its calibration are provided in
percent knowledge of the geotechnical proper- Cheon (2011) and Cheon & Gilbert (2013).
ties at the location of or along every foundation The model shows that the influence of the spa-
element. tial variability relative to the mean decreases with
Recent advances have been made in developing depth (Fig. 1), possibly reflecting the increasing
realistic models of spatial variability to account overburden stress damping variations in mineral-
for it in designing foundations and optimizing site ogy or depositional history. The effect of spatial
investigation programs (e.g., Keaveny et al., 1990, averaging in reducing variability for the depth-
Gambino & Gilbert 1999 and Valdez-Llamas et al., averaged strength increases with averaging length
2003). (Fig. 1).
An example of a random field model for the The horizontal correlation distance2 obtained
design capacity of deep foundations is shown in was between 2 and 6 km, and is therefore hundreds
Figures 1 to 4. The geologic setting is normally to of times greater than the vertical correlation dis-
slightly overconsolidated marine clays in 1,500 to tance (Fig. 2). Note that the correlation distance
3,000 m deep water in the Gulf of Mexico. The is much greater for the design undrained shear
random field model represents spatial variations in strength compared to that for individual measure-
the design shear strength. The design strength is the ments of undrained shear strength since the design
strength selected by a designer for the purposes of profile implicitly averages out small-scale varia-
foundation design based on all available laboratory tions (either real or due to measurement methods)
and field test data and geologic information at a and reflects larger-scale variations. Both the hori-
given location. The available data for this geologic zontal and vertical correlation distances are greater
setting included over 100 design profiles of und- for the depth-averaged versus the point strength
rained shear strength from site investigations with (Fig. 2). The horizontal correlation structure is
soil borings, jumbo piston cores, field vane tests best modelled as anisotropic, with a longer hori-
and Cone Penetration Tests (CPT). These design zontal correlation distance moving away from the
profiles are located as close as 100’s of meters to as continental shelf (in the direction of depositional
far as 1,000’s of kilometers from one another. flow) compared to moving along the continental
The three-dimensional random field model con- shelf (Fig. 3).
sists of two cross-correlated models for the design This model of spatial variability can be used to
undrained shear strength: one for the design strength support design decisions. An example application
at a particular depth below the sea floor (to calculate
end bearing) and one for the depth-averaged design
2. Correlation distance was defined here as the separation
strength from the sea floor to that depth (to calcu- distance at which the correlation coefficient is equal to
lated side shear). The model incorporates means 0.37 for an exponentially decreasing correlation coef-
and standard deviations that increase with depth, an ficient with separation distance. This correlation dis-
anisotropic spatial correlation structure, and hori- tance is one-half the scale of fluctuation defined by
zontal correlations that increase with depth. Details Vanmarcke (1983).
30
3 MODEL UNCERTAINTY
31
32
33
34
5 RELIABILITY-BASED DESIGN
35
To evaluate the required resistance factor, standard deviation, coefficient of variation and
Lacasse et al. (2013a, 2013b and 2013c) calculated Probability Density Function (PDF). The model
the annual probability of failure for piles on off- uncertainty was obtained by comparing the pre-
shore jackets designed with the API method and dicted to the measured axial pile capacity from
with the newer CPT-based methods. The goal relevant and reliable pile model tests. The NGI
was to make a recommendation on the appropri- database of “super pile” load tests NGI (2000;
ate resistance factor and minimum pile penetra- 2001) was used.
tion depth to use for the design of the piles on an The calibration used (1) the results of the deter-
offshore jacket. Table 1 lists the axial pile capacity ministic analyses giving the ultimate axial pile
methods considered. capacity with the characteristic strength parameters
The reliability analyses of the axial pile capacity (Qult char); (2) the probabilistic analyses giving the
methods included a statistical analysis of the soil PDF of the ultimate axial pile capacity (Qult mean);
parameters; statistical analysis of the model uncer- and (3) the results of the probabilistic analyses giv-
tainty for the different pile capacity calculation ing the annual probability of failure, Pf.
methods used; statistical analysis of the static Figure 15 is a simplification in two dimensions
(permanent) and environmental loads on the top of the overlap of the probabilistic ultimate pile
of the piles; deterministic analysis of the ultimate capacity (Qult) and probabilistic environmental
axial pile capacity, Qult; probabilistic analyses of load (Penv). The probability density function for the
axial pile capacity to obtain the PDF of the ulti- Penv was taken as the same for Pf1 and Pf2 in the cal-
mate capacity, Qult; calculation of the annual prob- culations. The calibration of the resistance factor
ability of failure by combining the statistics of the was coordinated with the definition of characteris-
loads and the probabilistic description of Qult; and tic design load and the characteristic soil strength
calibration of the safety factors (load and resist- profile used for the calculation of axial pile capac-
ance factor) for each pile capacity design method, ity. The calibration details are described in Lacasse
for a target annual probability of failure of 10−4. et al., 2013a.
Three sites, where jackets are currently under Table 2 presents the results of the calibration
design, were analyzed. For Jacket A, the soil con- of the resistance factor for the case study jackets
ditions are characterized by mainly clay layers with to achieve an annual probability of failure of 10–4.
intermittent thin sand and silt layers. For Jacket The resistance factor was obtained based on the
B, the soil consists of mainly dense to very dense axial pile capacity calculated with the characteris-
sand layers, with rather thin clay layers in between. tic undrained shear strength (Qult char). The load fac-
For Jacket C, the soil profile consists of alternat- tors were maintained at the recommended values
ing very dense sand and very stiff clay units. The in the design guidance, although the load factor at
parameters were estimated with statistical analyses the design point was smaller.
of the soil data, combined with well-documented For a given pile length, the calibrated resistance
correlations and experience (bias factors). factor varied with the pile design method. The fac-
An extended study of the model uncertainty tors reflect the varying influence of the uncertainty
was carried out for the different axial pile capac- in the soil parameters and of the model uncer-
ity calculation methods (Lacasse et al., 2013c). The tainties for the different methods. The results are
model uncertainty was expressed as a bias (mean), generally consistent, where the axial pile capacity
36
Site C
Site A (clay) Site B (sand) (clay & sand)
Method 90-m pile 26-m pile 40-m pile
37
38
39
40
41
42
D.V. Griffiths
Colorado School of Mines, Golden, CO, USA
University of Newcastle, NSW, Australia
Jumpol Paiboon
Colorado School of Mines, Golden, CO, USA
Jinsong Huang
University of Newcastle, NSW, Australia
Gordon A. Fenton
Dalhousie University, Halifax, NS, Canada
ABSTRACT: The homogenized stiffness of geomaterials that are highly variable at the micro-scale has
long been of interest to geotechnical engineers. The purpose of this study is to investigate the influence
of porosity and void size on the homogenized or effective properties of geomaterials. A Random Finite
Element Method (RFEM) has been developed enabling the generation of spatially random voids of given
porosity and size within a block of geomaterial. Following Monte-Carlo simulations, the mean and stand-
ard deviation of the effective property can be estimated leading to a probabilistic interpretation involving
deformations. The probabilistic approach represents a rational methodology for guiding engineers in the
risk management process. The influence of block size and the Representative Volume Elements (RVE) are
discussed, in addition to the influence of anisotropy on the effective Young’s modulus.
43
44
3 CONTROLLING POROSITY
Thereafter, any element assigned a random field
The random field generator in the RFEM model
value in the range |Z| > zn/2 is treated as intact material
known as the Local Average Subdivision method
with a Young’s modulus and Poisson’s ratio given by
(LAS) (Fenton & Vanmarcke 1990) is used in this
E0 = 1 and υ0 = 0.3, while any element where |Z| ≤ zn/2
paper to model spatially varying voids properties.
is treated as a void element with Young’s modulus
The target mean porosity n is obtained by using the
and Poisson’s ratio given by E0 = 0.01 and υ0 = 0.3
standard normal distribution shown in Figure 3.
(100 times smaller than the surrounding intact
A single value of the random variable Z is initially
material). As can be seen in Figure 4, for the case
assigned to each element of the finite element
when n = 0.2, the results show a small influence of
mesh. Once the standard normal random field
the arbitrarily selected Young’s modulus of the void
values have been assigned, cumulative distribution
elements. In the current work, a void stiffness one
tables Φ (suitably digitized in the software) are then
hundred times less than the surrounding intact mate-
used to estimate the value of the standard normal
rial gave reasonable (and stable) results. The nature
variable zn/2 for which
of random fields is that the mean porosity is under
the user’s control, but the porosity of each individual
Φ ( )
n2 Φ ( 0) n2 (1) simulation processed by the Monte-Carlo method
will vary from one simulation to the next.
where Φ is the cumulative normal distribution
function, and n is the target porosity as shown in
Figure 3. 4 CONTROLLING OF VOID SIZE
ρ ( 2 τ θ) (2)
45
δx δy δ
εx = , εy = , εz = z (5)
L L L
Q
E= (6)
Lδ z
Figure 5. Typical simulations showing generation of δx
voids at (a) low and (b) high spatial correlation lengths υx = (7)
θ (n = 0.2 in both cases). δz
δy
of ρ varies from 0 to 1. Points close together are υy = (8)
δz
strongly correlated and therefore likely to belong
to the same void. In the limiting case of θ → 0, the
random field value changes rapidly from point to where E = the effective elastic Young’s modulus,
point delivering numerous small voids. At the other Q = stress loading at the top side, υx and υx = the
extreme as θ → ∞, the random on each simulation effective Poisson’s ratios based on the displacement
becomes increasingly uniform with some simula- in the x- and y-directions respectively.
tions representing entirely intact material and In each simulation, the effective Young’s modu-
other consisting entirely of voids. For example as lus is normalized as E/E0 by dividing by the intact
shown in Figure 5, the models show typical simula- Young’s modulus E0. In the current study, fol-
tions of different void clustering for two materials lowing some numerical experiments as shown in
with the same mean porosity. Figure 6, it was decided that 1000 simulations for
each parametric combination would deliver rea-
sonably repeatable results. In this study, we have
5 MONTE-CARLO SIMULATIONS expressed the spatial correlation length in dimen-
sionless form
A “Monte-Carlo” process is combined with the
RFEM and repeated until stable output statistics θ
Θ= (9)
are achieved. The primary outputs from each elas- L
tic analysis are the vertical and horizontal defor-
mations of the block δz, δx and δy. Although all where L is the width of the loaded element
simulation use the same θ and n, the spatial loca- (L = 50).
tion of the voids will different each time. In some
cases, the voids may be located just below the top
of the block leading to a relatively high δz. While
in others, the voids may be buried in the middle of
the block leading to a relatively low δz. Following
each simulation, the computed displacements δz, δx
and δy are converted into the “effective” values of
Young’s modulus and Poisson’s ratio as follows
Based on Hooke’s law,
1
εx (σ x − υ (σ y + σ z ))
E
1
εy (σ y − υ (σ z + σ x )) (3) Figure 6. Sensitivity of the mean effective Young’s
E modulus as a function of the number of simulations for
1 n = 0.2 and Θ = 0.4. It was decided that 1000 simulations
εz (σ z − υ (σ x + σ y ))
E would deliver reasonably repeatability.
46
1 0.2 0.2
bigger, but at a slower rate for higher values of Θ.
2 0.2 0.7
In both Figures 8, it is noted that the influence of
3 0.7 0.2
Θ on block statistics is greater than that of n. The
4 0.7 0.7
RVE depends more on spatial correlation length
than porosity.
7 RESULTS OF RFEM
47
μE
The result obtained from Equations 7 and 8 for μK = (10)
the effective Poisson’s ratio were in good agreement 3(1 − 2 μυ )
as expected for the range of n and Θ considered.
μE
In the isotropic material model, the two Poisson’s μS = (11)
ratios are essentially identical after Monte-Carlo 2(1 + μυ )
simulation; however the results are based on an
average to account for any small differences. The where μK = the mean effective bulk modulus,
plots shown in Figures 11 and 12 give the mean μS = the mean effective shear modulus.
48
Figure 15. Analysis of tied freedom in four “cubic element test” models with voids: (a) isotropic model, (b) aniso-
tropic model along x-axis, (c) anisotropic model along y-axis and (d) anisotropic model along z-axis.
49
ACKNOWLEDGEMENT
REFERENCES
50
51
ABSTRACT: The cohesion and friction angle of soil are two common random variables encountered
in geotechnical reliability analysis. Its joint distribution can be fully specified with marginal distributions
and a copula function. In this study, a probabilistic framework is suggested to construct, calibrate and
rank probability models for failure probability estimation. It is shown that calibrating the marginal dis-
tributions and copula function may miss the optimal model parameters. The optimal probability model
is problem specific, and the commonly adopted models based on the Gaussian copula function may not
be optimal. When copula function is used to specify the joint distribution of cohesion and friction angle,
many existing algorithms such as the first order reliability method cannot be directly used for failure prob-
ability calculation. Monte Carlo simulation can be applied to general performance functions and hence is
superior to the existing one-dimensional integration method. Both the marginal distribution and copula
function could significantly affect the failure probability calculation. If the probability model is arbitrarily
chosen, the failure probability may be severely overestimated or underestimated. In general, the effect of
probability model on failure probability estimation is more obvious when the failure probability is small.
The effect of probability model on failure probability calculation also depends on the deterministic model
and the available data. Increasing the amount of the measured data is likely to provide more constraint
to the marginal distributions and hence reduce its effect on failure probability estimation. The effect of
copula function on failure probability calculation decreases as there is less correlation between the cohe-
sion and the friction angle.
53
2 PROBABILITY MODELS FOR SHEAR where μc and μϕ are the mean values of c and ϕ,
STRENGTH PARAMETERS respectively, and σc and σϕ are the standard devia-
tions of c and ϕ, respectively.
Let c and ϕ denote the cohesion and friction angle, Comparing Eq. (3) with (4), we can see that the
respectively. Let F1(c) and F2(ϕ) denote the Cumu- correlation coefficient is a function of the copula
lative Distribution Function (CDF) of c and ϕ, model. Previously, the Gaussian copula is often
54
{ }
−θ θ −1−
× ⎡⎣ F1 (c ) ⎦⎤ + ⎣⎡ F2 (ϕ ) ⎤⎦ −1 θ
f1 ( ) f2 ( ) (5) i i
where d1 and d are the first and second elements
2
of di, respectively.
Based on the principle of maximum likelihood,
where f1(c) and F1(c) in this case are respectively the
the optimal values of Θ can be obtained by maxi-
PDF and CDF of a normal variable with a mean
mizing the likelihood function, or equivalently, by
of μc and a standard deviation of σc, and f2(ϕ)
maximizing the logarithm of the likelihood func-
and F2(ϕ) are respectively the PDF and CDF of a
tion as follows
normal variable with a mean of μϕ and a standard
deviation of σϕ. As an example, f1( c) and F1(c) can n
be written as follows L ∑
i =1
f ( i
|Θ ) (10)
1 ⎡ ( x − )2 ⎤
f1 (c ) = p ⎢−
exp c
⎥ (6) If Eq. (10) is used to find the optimal values of
2π σ c ⎢⎣ 2σ c2 ⎥⎦ Θ, the marginal distributions and the copula func-
tion are calibrated simultaneously. Thus, if there is
c any interaction between the optimal marginal dis-
F1 (c ) = ∫ f ( z ) dz (7) tributions and copula function, such interaction
−∞ 1
−
can be readily considered. Nevertheless, estimating
Sklar’s theory implies that the structure of a the values of Θ using Eq. (10) requires solving the
probability model for cohesion and friction angle optimization problem in a five-dimensional space,
can be constructed in two steps: (1) determining which could be non-trivial. To reduce the computa-
the marginal distributions of c and ϕ, and tional challenge involved in the optimization work,
(2) determining the copula function between c and one can assume there is no interaction between
ϕ. After the structure of the probability model is the optimal marginal distributions and the copula
determined, the next step is how to calibrate the function. In such a case, the marginal distributions
model parameters of the probability model, as and copula function can be calibrated separately
described in the following section. (Joe & Xu 1996). For the shear strength probability
model studied here, the optimal values of μc, σc, μϕ,
σϕ, and θ can be obtained by respectively maximiz-
3 MAXIMUM LIKELIHOOD ing the following three likelihood functions
CALIBRATION OF PROBABILITY
MODELS
( )
n
μc c ∏ f1 d1i | μc ,σ c (11)
In the probability model for shear strength param- i=1
eters, the parameters to be calibrated include μc, μϕ,
( )
n
σc, σϕ and θ. Let d = {d1, d2} denote a measurement of l μϕ ϕ ∏ i
f2 d | μϕ , σ ϕ (12)
c and ϕ. Let d1, d2, d3, …, dn denote n measurements i=1
55
56
Table 3. Calibrated parameters of 16 probability models based on the given data (model parameters calibrated
simultaneously).
Model Marginal
no. Copula distributions* AIC θ μc σc μϕ σϕ
*n,n: both c and ϕ follow the normal distribution; l,n: c and ϕ follow the lognormal and normal distributions, respectively;
n,l: c and ϕ follow the normal and lognormal distributions, respectively; l,l: both c and ϕ follow the lognormal distributions.
**For the FGM copula, the value of θ is bounded between −1 and 1 (Trivedi et al. 2005). This constraint is considered
in the maximum likelihood method.
57
Model Marginal
no. Copula distribution AIC θ μc σc μϕ σϕ
and 124.1 kPa, respectively. However, when the 4. For a given copula structure, the model assu-ming
marginal distribution of c is changed to be log- c and ϕ are lognormally distributed is most sup-
normal, the mean and standard deviation of c is ported by the data. Comparing models with dif-
changed to be 196.7 kPa and 351.9 kPa, respec- ferent copula functions, Model 16 with the FGM
tively. The optimal model parameters are sensi- copula function has the smallest AIC, and hence
tive to the assumed marginal distribution. is most supported by the data. Models based on
2. When the marginal distributions of cohesion the commonly used Gaussian copula function
and friction angle are both normal distributions are not the optimal probabilistic models.
but when the copula function is assumed to be
Gaussian and Clayton, respectively, the mean
values of c and ϕ are both changed, indicating 7 IMPACT OF PROBABILITY
that the optimal marginal distribution is also MODELS ON FAILURE
affected by the copula function. Hence, it may PROBABILITY CALCULATION
not be appropriate to calibrate the marginal
distribution and copula function separately. To 7.1 Example 1
verify this conclusion, Table 4 shows the cali- To investigate the effect of probability model on
bration results of the 16 models when marginal failure probability calculation, consider an infinite
distributions and the copula function are cali- slope as shown in Figure 1. Its factor of safety can
brated separately based on Eqs. (11)–(13). For be calculated as follows:
the same model the AIC value in Table 4 is
generally larger than that in Table 3, indicating
c γ H s 2 α tan φ
that the model calibrated ignoring the interac- Fs = (20)
tion between marginal distributions and copula γ H i α cos α
function is less supported by the data.
3. When the copula function is Gaussian, the four where H = thickness of soil layer; α = slope angle;
probability models with different marginal dis- γ = unit weight of the soil. In this example, H = 5 m
tributions have different AIC values, indicat- and γ = 17 kN/m3.
ing that the validity of each probability model Monte Carlo simulation is first used to calcu-
is affected by the marginal distributions. Also, late the failure probability of the slope when dif-
when the marginal distributions are the same, ferent probability models are adopted and as the
the AIC values of the model vary with the slope angle varies. In the Monte Carlo simulation,
copula function, indicating the validity of each 1,000,000 samples are generated. The obtained
model is also affected by the copula function. failure probabilities are summarized in Table 5.
58
Table 5. Failure probability calculated based on different probability models (example 1).
pf
Model Marginal 1-D Monte Carlo 1-D Monte Carlo 1-D Monte Carlo
no. Copula distribution integration simulation integration simulation integration simulation
59
Figure 3. Effect of copula function on failure probability calculation (example 1) (a) models 1, 5, 9, 13;
(b) models 2, 6, 10, 14; (c) models 3, 7, 11, 15; (d) models 4, 8, 12, 16.
60
Figure 4. Effect of marginal distributions on the distribution of samples (example 1): (a) model 1; (b) model 2;
(c) model 3; (d) model 4.
61
Qu 0 B ⋅ Nγ c ⋅ Nc 2 D f ⋅ Nq (21)
Figure 6. Strip footing analyzed in example 2.
2
⎛ ⎛ π φ ⎞⎞
Nq eπ φ
⋅ ⎜ tan ⎜ + ⎟ ⎟ (22)
⎝ ⎝ 4 2⎠⎠
bottom. In this example, Df = 0.5 m, and B = 2 m,
Nγ 1.8 ( Nq − 1) ⋅ tan φ (23) γ1 = γ2 = 19.5 kN/m3.
As the performance function here is complex,
the one-dimensional integration method can
Nq − 1 hardly be applied. However, the Monte Carlo sim-
Nc = (24)
tan φ ulation algorithm is still applicable. Table 6 sum-
marizes the failure probability calculated using
where Df = embedded depth of foundation 16 different probability models (N = 1,000,000).
(m); B = foundation width (m); γ1 = unit weight Similar to Figures 2, 3 and Figures 7, 8 show the
of the soil below the foundation bottom and effect of marginal distribution and copula function
γ2 = unit weight of the soil above the foundation on failure probability calculation in this example,
62
pf
Model Marginal
no. Copula distribution q = 50 kPa q = 80 kPa q = 110 kPa
*0 denotes that no failure sample is found during the conditional sampling with 1,000,000 samples.
Figure 7. Effect of marginal distributions on failure probability calculation (example 2): (a) models 1, 2, 3, 4;
(b) models 5, 6, 7, 8; (c) models 9, 10, 11, 12; (d) models 13, 14, 15, 16. (failure probability of 0 in Table 6 is plotted
as 1.0 × 10−6).
respectively. The phenomena observed in Figures 7 Figures 9(a) and (b) show the calculated fail-
and 8 are generally similar to those observed in ure probabilities based on 16 models for Exam-
Figures 2 and 3, i.e. the effect of marginal distribu- ple 1 and Example 2, respectively. The difference
tions and copula functions on failure probability between calculated failure probabilities based on
calculation is more obvious when the failure prob- the 16 probability models in Example 2 is even
ability is small. larger than that observed in Example 1, indicating
63
7.3 Example 3
In the above two examples, the probability mod-
els have important effect on the failure probability.
The third example shows a case where such an
effect is less obvious. In this example, the data are
adopted from Xu & Yang (1998) about the shear
strength parameters of clay in Shanghai, as shown
in Table 7. In this example, there are 81 measure-
ments, compared to only 25 measurements in
Example 1. The same 16 models are considered in
this example, and the calibration results are shown
in Table 8.
Among the 16 models considered, Model 8 has
the smallest AIC value and hence is the optimal
one, in which both c and ϕ are assumed to be log-
normal and the copula function is Clayton. The
optimal model in this example is different from
that found in Example 1, indicating the optimal
probability model is problem specific.
Using the calibrated samples, the failure prob-
ability of the infinite slope in Example 1 is calcu-
lated again. Figure 10 shows the effect of marginal
Figure 9. Failure probabilities calculated based on distribution on failure probability calculation. The
16 models: (a) example 1; (b) example 2; (c) example 3. difference between failure probabilities calculated
(Failure probability of 0 in Table 6 is plotted as 1.0 × 10−6). based on different probability models is generally
64
Model Marginal
no. Copula distribution AIC θ μc σc μϕ σϕ
no more than two orders of magnitude, which are mum difference between failure probabilities calcu-
less than those observed in Example 1, indicating lated based on different copula functions occurs in
the effect of marginal distribution is less obvious in Figure 11(d), which is less than one order of mag-
this example. Figure 11 shows the effect of copula nitude and is also significantly than that observed
function on failure probability calculation. This is in Example 1. Note the copula function measures
possibly because there are more measurements in the dependence relationship between the cohesion
this example, which provide more constraint on and the friction angle. To find out the reason why
the shape of the marginal distributions. The maxi- the effect of copula is less obvious in example, the
65
Figure 11. Effect of copula function on failure probability calculation (example 3) (a) models 1, 5, 9, 13; (b) models
2, 6, 10, 14; (c) models 3, 7, 11, 15; (d) models 4, 8, 12, 16.
66
67
68
ABSTRACT: This paper presents a new geotechnical design concept, called Robust Geotechnical Design
(RGD). The new design methodology seeks to achieve design robustness and economics, in addition to
meeting safety requirements. Here, a design is considered robust if the variation in the system response
(i.e., failure probability) is insensitive to the statistical characterization of noise factors, such as uncertain
geotechnical parameters. It can be shown that when safety requirements are met, cost and robustness are
conflicting objectives with a trade-off, meaning that a single best design may be unattainable. In such a
case, a multi-objective optimization considering cost and robustness can be performed to obtain a suite of
best designs within the solution space in which safety is guaranteed. To this end, the suite of best designs
forms the Pareto Front, which is a useful tool to aid in the design decision making within the RGD meth-
odology. The new design methodology is illustrated with examples of spread foundation designs and rock
slope designs, through which the significance of the RGD methodology is demonstrated.
69
70
71
72
73
74
75
5 CONCLUDING REMARKS
76
77
78
ABSTRACT: Any engineering and construction project involves dealing with risk, which has often a
geotechnical origin. Therefore, the objective of applying Geotechnical Risk Management (GeoRM) is
contributing to the engineering, construction, and maintenance of successful projects, despite the inher-
ently uncertain character of its ground conditions. However, for realizing these projects in an effective
and cost-efficient way, GeoRM has to be well-integrated within Project Risk Management (ProjectRM).
As ISSMGE TC304 Task Force 3 revealed, a complete GeoRM application that is integrated with Pro-
jectRM, seems not yet common practice in a lot of project phases, in a lot of projects, in a lot of coun-
tries. This requires crossing often hidden boundaries between risk concepts, disciplines, project phases,
scopes, project types and even boundaries between industries and countries. Finally and perhaps first
of all, we have to cross our own mental boundaries. Within this paper these boundaries are highlighted
and crossed, for catching benefits of learning, effectiveness, and cost-efficiency. The promise of well-
integrated GeoRM in ProjectRM, by crossing the mentioned boundaries, is substantially contributing to
successful construction projects from a geotechnical engineering and construction perspective, with fewer
cost overruns, less delays, less non-conformances of safety and quality, and therefore happier project
participants and stakeholders.
79
80
81
82
83
84
85
86
87
88
Done
1 Generate and protect value Make all geotechnical risks in each project phase
explicit, with risk effects remediation measures
2 Participate in decision making Make a geotechnical risk file from the start of the
in all project phases project, use it for decisionmaking
3 Make geotechnical uncertainty Include geotechnical sensitivity analyses with
explicit margins in project reports
4 Work systematically, structured, Include GeoRM explicitly in the project planning
and in time and reserve resources for it
5 Use all available information Work from a general level to a detailed level, from
using geological maps to geotechnical monitoring
6 Work transparently together with Indicate and communicate any dependencies of
all stakeholders geotechnics with other disciplines in the project
7 Include the role of the human Make differences in organizational culture of all
factor involved project parties visible and feasible
8 Use experiences and lessons for Use all available and relevant project evaluations, risk
continuous improvement checklists, and experiences from professionals
89
90
91
92
93
Mechanical industry How to increase productivity, despite the unique character of construction
projects and the related risks.
IT-industry How to use simulation tools for optimizing construction processes under risky
conditions.
Space industry How to integrate risk management in systems engineering in construction projects.
How to develop risk management as part of continuous organizational
improvement programs.
Insurance industry How construction insurance claims can be avoided or reduced, by using qualified
and competent people.
Chemical and nuclear industry How to develop clear risk benchmarks
How to combine risk and safety management
Financial services industry How to consider risk in a holistic and integral way by applying Enterprise Risk
Management (ERM).
Consulting industry How to apply proven project management techniques within ProjectRM and
GeoRM processes.
94
95
96
97
98
99
H.N. Wong
Civil Engineering and Development Department, Geotechnical Engineering Office, Government of the
Hong Kong Special Administrative Region, Hong Kong, China
ABSTRACT: Quantitative Risk Assessment (QRA) enables decision-makers to exercise informed and
rational judgement in evaluating and managing risk. It is commonly acknowledged that Hong Kong has
showcased not only the state-of-the-art in, but also the practicability of, quantification and management
of landslide risk. The author would caution that this is a simplistic view. With reference to Hong Kong’s
experience, it is highlighted that there are different categories of landslide problems, each with its own
factors that combine to give rise to risk. As to whether the risk is quantifiable and manageable, different
categories have their particular circumstances and the answers are not the same. It is credible to manage
landslide risk without quantification, although risk quantification can undeniably aid risk management.
Landslide risk quantification may not be practicable in all cases. Even if landslide risk is quantifiable, it is
important not to lose sight of the degree of uncertainty involved.
101
Table 1. The slope safety system in Hong Kong (based on Malone 1998).
Policing
– Checking new slope works
– Slope maintenance audits
– Recommending safety clearance
of vulnerable squatters and unauthorised
structures threatened by hillsides
– Exercising geotechnical control through
input in land use planning
– Safety screening studies and
recommending statutory repair orders
for private slopes
Works projects
– Retrofitting substandard government
man-made slopes
– Natural terrain landslide mitigation
and boulder stabilisation works
Research and setting standards
Education and information
– Slope maintenance campaigns
– Risk awareness programmes
and personal precaution campaigns
– Information services
– Landslip warning and emergency services
102
103
104
105
106
Table 2. Breakdown of landslides affecting different facilities in Hong Kong since 1994.
*A ‘block’ is a multi-storey building, which may comprise up to several dozen of flats/units. A ‘house’ is typically within
3 storeys, which comprises several flats/units.
107
108
109
110
111
W. Du & G. Wang
Department of Civil and Environmental Engineering, The Hong Kong University of Science and Technology,
Clear Water Bay, Hong Kong SAR, China
115
2
⎛ a ⎞ ⎛ a ⎞
ln D = −1.56 − 4.58 ⎜ c ⎟ − 20.84 ⎜ c ⎟
⎝ PGA
G ⎠ ⎝ PGA
G ⎠
3 4
⎛ a ⎞ ⎛ a ⎞
+ 44.75 ⎜ c ⎟ − 30.5 ⎜ c ⎟
⎝ PGA
G ⎠ ⎝ PGA
G ⎠
− 0.64 ln( PG ) + 1.55 ln(l ( )
σ lln D = 0.411 + 0.52 ( ac /P
PGA
P ) (4)
2
⎛ a ⎞ ⎛ a ⎞
ln( ) = −0.74 − 4.93 ⎜ c ⎟ − 19.91⎜ c ⎟
Figure 1. Illustration of Newmark displacement ⎝ PGA
G ⎠ ⎝ PGA
G ⎠
3 4
with critical acceleration ac = 0.2 g. (A) Earthquake ⎛ a ⎞ ⎛ a ⎞
acceleration-time history. (B) Velocity of sliding block ver- + 43.75 ⎜ c ⎟ − 30.12 ⎜ c ⎟
sus time. (C) Displacement of sliding block versus time.
⎝ PGA
G ⎠ ⎝ PGA
G ⎠
− 1.3 ln( PGA
G ) + 1.04
04 l ( ) + 0.67 ln( Ia )
σ ln D = 0.2 + 0.79 ( ac /PGA )
Ground Motion Prediction Equations (GMPEs) (5)
developed in the past four decades for IMs such
as PGA, PGV, Ia, etc. However, the epistemic 5. [PGA, M] RS09 model (Rathje & Saygili,
and aleatory uncertainties of the Newmark dis- 2009):
placement models have not yet been thoroughly
2
investigated. ⎛ a ⎞ ⎛ a ⎞
This study aims at quantifying the epistemic ln( ) = 4.89 − 4.85 ⎜ c ⎟ − 19.64 ⎜ c ⎟
⎝ PGA
G ⎠ ⎝ PGA
G ⎠
uncertainty and aleatory variability of Newmark 3 4
displacement models. Ten recently developed ⎛ a ⎞ ⎛ a ⎞
+ 42.49 ⎜ c ⎟ − 29.06 ⎜ c ⎟
Newmark displacement models are used in this ⎝ PGA
G ⎠ ⎝ PGA
G ⎠
study and listed as follows: + 0.72 ln( G ) + 0.89( − 7)
w
1. [PGA] RS08 model (Saygili & Rathje, 2008):
7 ( ac /P
PGA ) 0.54 ( ac /P ) 2
σ lln D = 0.73 + 0.79 PGA PGA
P
2 (6)
⎛ a ⎞ ⎛ a ⎞
ln D = 5.52 − 4.43 ⎜ c ⎟ − 20.39 ⎜ c ⎟
⎝ PGA
G ⎠ ⎝ PGA
G ⎠
3 4
6. [PGA] J07 model (Jibson, 2007):
⎛ a ⎞ ⎛ a ⎞
+ 42.61⎜ c ⎟ − 28.74 ⎜ c ⎟
⎝ PGA
G ⎠ ⎝ PGA
G ⎠ ⎛ a ⎞
+ 0.72 ln( PGA
G ) log10 ( D ) .215 + 2.
.215 log10 ⎜1 − c ⎟
⎝ PGAG ⎠
σ lln D = 1.13 (2) ⎛ ac ⎞
− 1. log10 ⎜ ⎟
⎝ PGA
G ⎠
2. [PGA, Ia] RS08 model (Saygili & Rathje, σ log10 D = 0.51 (7)
2008):
116
117
118
119
120
ABSTRACT: A reliability analysis of 640 m long retaining wall for an embedded highway has been
carried out. The soil profile at the site mainly consists of alternation of three cohesive and cohesionless
soil layers. The soil investigations have been conducted by the borings (SPT N-values) of about 100 m
interval. It was found that the most dominant source of uncertainty to control the performance of the
retaining wall is the thickness of the soil layers which are different from a section to another. The layers
thicknesses are modeled as correlated random processes whose thicknesses are known at the location of
the borings. The thickness of a layer at a certain location is estimated by Co-Kriging method and the
conditional simulation based on it with quantified uncertainty. The design calculation model error is
also evaluated by comparing the calculated results with laboratory model experiment results and actual
site measurements. The response surfaces of the performance of the retaining wall as a function of layer
thicknesses and soil property variations are obtained by extensive design calculations of the retaining
wall. Finally, the reliability evaluation of the retaining wall along the excavation has been carried out by
Monte Carlo Simulation (MCS). The final result takes into account the spatial variability of soil proper-
ties both soil mechanical property variation and the layer thicknesses, the statistical estimation error due
to sparse soil investigations, and design calculation model errors. The quantitative contributions of each
component are also presented.
121
122
4 UNCERTAINTY ANALYSIS
123
124
Table 5. Model error in the elasto-plastic method. Statistical estimation NB 15.0 6.43 0.43
Error Γ2 NDg 19.0 1.83 0.10
Item Mean COV n (Normal) NDc 12.0 1.13 0.09
Spatial variability NB 15.0 7.69 0.51
Case histories Max disp. 0.70 0.29 58 Λ2 NDg 19.0 2.20 0.12
Max BM 0.69 0.43 42 (Normal) NDc 12.0 1.35 0.11
Model Max disp. 1.08 0.67 38 Model error Disp. 1.08 0.67 0.62
experiments Max BM 0.73 0.44 24
δME BM 0.73 0.44 0.60
(log normal)
(note 1 error = (observed)/(calculated); (note 2 BM:
bending moment, disp.: displacement).
125
figures, contribution of each uncertain source are properties are quite limited on the reliability of
illustrated for case θh = 100 (m). the retaining wall.
The following observations can be made for
these cases.
REFERENCES
1. Pf for both the maximum horizontal displace-
ment and the bending moment takes the maxi- Honjo, Y.,Y. Otake & Kato H. 2012. A simpli-fied scheme
mum value between distance 300 to 450 (m). In to evaluate spatial variability and statistical estima-
this section, the thickness of Dg layer is thicker, tion error of local average of geotechnical parameters
which makes the earth pressure lager. in reliability analysis. J. Geotechnical Engineering
2. The accuracy of the layer thickness prediction (JSCE) 68(1):41–55. (in Japanese).
Honjo, Y. & Otake Y. 2012. Verification of statistical esti-
is less for θh = 50 (m) than for 100 (m). It is mation error evaluation theory of local averages of
also less in the area that are distanced from the geotechnical parameters. J. Geotechnical Engineering
observations are made. However, the difference (JSCE) 68(3):475–490. (in Japanese).
is rather small for θh = 50 (m) and for 100 (m). Honjo, Y. & Otake Y. General and local estimation of
3. The uncertainty sources that are critical are the local average of geotechnical parameters in reliabil-
model error and the statistical estimation error of ity analysis. Proc. APSSRA, Phoon, K.K., Beer, M.,
the layer thicknesses. The spatial variability and Quek, S.T., Pang, S.D. (editors), 236–241.
the statistical estimation error of SPT N-value Honjo, Y. & Otake Y. 2013. A Simple Method to Assess
of each layer have only negligible impacts. the Effects of Soil Spatial Variability on the Perform-
ance of a Shallow Foundation. ASCE Geotechnical
Special Publication No. 229, Foundation Engineering in
the Face of Uncertainty (honoring Fred H. Kulhawy),
6 CONCLUSION J.L. Whithiam, K.K. Phoon and M.H. Hussein (eds.),
385–404.
The sources of uncertainty concerning retaining Jounel, A.G. & Huijbregts C.J. 1978. Mining Geostatis-
wall design have been identified and quantified tak- tics, Academic Press.
ing 640 m long embedded highway retaining wall as Otake, Y. & Honjo Y. 2012. Reliability based design on long
an example. The reliability analysis scheme proposed irrigation channel considering the soil investigation
by Honjo, Otake and Kato (2012) has been applied locations, Proceedings of Geo Congres, 2836–2845.
Otake, Y. 2012. Development of geotechnicl reliability
to the reliability evaluation of the retaining wall. based design and chllenges in geotechnical design,
The result of the evaluation has shown that the PhD dissertation submitted to Gifu university, 215 pp.
two major sources of uncertainty are the statisti- (in Japanese).
cal estimation error of the layer thicknesses and Vanmarcke, E.H. 1977. Probabilistic modeling of soil
the model error of the elasto-plastic method for profiles. J. of Geotechnical Engineering, ASCE 103
retaining wall design. The impacts of the spatial (GT11): 1227–1246.
variability of soil properties within each layer,
and the statistical estimation error of the soil
126
K. Huang
Department of Civil and Architectural Engineering, Shenzhen Research Institute, City University of Hong Kong,
Hong Kong, China
Z. Cao
State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University,
Wuhan, P.R. China
Y. Wang
Department of Civil and Architectural Engineering, Shenzhen Research Institute, City University of Hong Kong,
Hong Kong, China
ABSTRACT: This paper aims to identify underground soil stratum and classify the soil types using
Bayesian approaches and Cone Penetration Tests (CPT). The uncertainty in the CPT-based soil classifica-
tion using the Robertson chart is modeled explicitly in the Bayesian approaches. The proposed Bayesian
approaches contain two major components: a Bayesian model class selection approach to identify the
most probable number of underground soil layers and a Bayesian system identification approach to esti-
mate the most probable layer thicknesses and determine the soil types simultaneously. Equations are
derived for the Bayesian approaches, and the proposed approaches are validated using a set of simulated
CPT data. It has been shown that the proposed approaches correctly identify the underground soil strati-
fication and classify the soil type of each layer.
1 INTRODUCTION
127
128
129
130
7 SIMULATED EXAMPLE
131
REFERENCES
132
J.H. Li
Centre for Offshore Foundation Systems, The University of Western Australia, Perth, Australia
Harbin Institute of Technology Shenzhen Graduate School, Shenzhen, China
L.M. Zhang
The Hong Kong University of Science and Technology, Hong Kong, China
Y. Zhao
China Merchants Shekou Industrial Zone Co., Ltd., Shenzhen, China
ABSTRACT: Cracks in soils are three-dimensional (3D) and provide important preferential path-
ways for rainfall infiltration. The geometrical properties of 3D cracks are crucial parameters for analysis
of seepage in discontinuous soils or rocks. However 3D crack planes in soils are still largely unknown
because the cracks in soils are prone to disturbance and sensitive to moisture content. This paper presents
a method to characterize 3D crack planes in soils by employing a nondestructive Computer Tomography
(CT) test. The traces of cracks are first obtained from sections along three directions in the CT test. Then
three criteria are proposed to determine a crack plane based on these crack traces. Finally the random
crack network in the soil is established by assembling the planes. The probability distribution and statisti-
cal parameters can be obtained based on the obtained crack planes.
133
3 TRACES IN THREE-DIMENSIONAL
SPACE
134
135
136
137
138
ABSTRACT: Linear regression is often used to obtain the sample mean value that can be regarded as
the population mean. A common misuse of linear regression is on the standard error of the mean, which
decreases with sample number and has nothing to do with the real variance of regression parameters.
Derivations are reported to calculate the “real” variance and numerical tests are used to verify these
derivations. Further, linear regression is used to analyze the tri-axial test data. The results demonstrate
that an equation transformed from traditional mean stress relation can lead to mean and variance values
with higher accuracy in linear regression, especially under the condition of small sample size (smaller than
1000).
139
140
141
Var( ) = μX 2σ b 2 + σ a2 + σ ε2
Var( ) b2 Var( (9) Var( ) = μX 2σ b 2 + σ a2 + 2 ρaabσ aσ b μX
Var( ) b2 Var(
(11)
where b is the mean of b calculated by Eq. 3, which
is considered to be of high accuracy; σε is normally Eq. 11 can be solved by quadratic regression,
regarded as a negligible value. So there are two taking μX as independent variable and [var(Y) −
unknowns σa and σb. b2var(X)] as dependent variable.
If the sample is divided into n groups, and To solve for var(a) and var(b), the following
Eq. 9 is applied to each group, n equations can be equation can be also used in LR instead of Eq. 9
obtained and used to solve the two unknowns σa
and σb. Var( ) b2 Var( ) 1
Obviously, the error of var(X) will affect the = σ b2 + σ a2 (12)
μX 2 μX 2
accuracy of the solved σa and σb. To reduce
the value of var(X), n groups of data should be
The choice of Eq. 9 or 12 depends on the accu-
selected according to the X value, but not selected
racy indicated by the standard errors while apply-
randomly. For example, σa and σb can be solved by
ing these two equations.
using the following procedure:
1. dividing the sample into m groups accord- 4.2 Numerical test of the estimation of variance
ing to the X value with an interval of
[max(X) – min(X)]/m. Taking a, b as normal distributions and X as a uni-
form distribution with input settings (Tests 4 & 5
–2 each group data, a pair of [μX , var(Y) −
2
2. For
b var(X)] can be obtained and a total m pairs of in Table 1): a = 10, σa = 2, b = 2, σb = 0.4, xmin = 0,
data can be obtained. xmax = 100, the results of numerical test are listed
3. Taking μX2 as independent variable and [var(Y) − in Table 5. Referring to Table 5, Eq. 9 can lead to
– accurate var(b), but unreliable var(a), which has
b 2var(X)] as dependent variable, σa2 and σb2 can
be solved using linear regression (Eqs. 3 and 4). great difference with the input value and is of high
standard error. Whereas, Eq. 12 can lead to values
The above procedure is valid for general cases, with higher accuracy for var(a) and similar accu-
such as for X following a uniform distribution. If X racy for var(b) comparing to Eq. 9.
is m groups of repeated data and m groups of data If X is taking as m groups of repeated data,
are selected according to X value, var(X) is 0 and the results of numerical tests are listed in Table 6.
Eq. 9 can be simplified as (if σε is negligible), According to Table 6, Eq. 12 can lead to values
with higher accuracy comparing to Eq. 10. Using
Var(Y ) = μX 2 b
2
+ σ a2 (10) X as m groups of repeated data instead of X with
n* m♥ μ (Eq. 4) Se2 (Eq. 6) μ (Eq. 3) Se2 (Eq. 5) μ (Eq. 3) Se2 (Eq. 5) μ (Eq. 4) Se2 (Eq. 6)
142
σa2
σb
2
σa2 σb2
n* m♥ μ (Eq. 4) Se2 (Eq. 6) μ (Eq. 3) Se2 (Eq. 5) μ (Eq. 3) Se2 (Eq. 5) μ (Eq. 4) Se2 (Eq. 6)
143
8 40 N (50, 100) N (0.36, 0.002) 50, 100, 150, 200 ∼N (142, 842) ∼N (2.03, 0.029)
9 400 N (50, 100) N (0.36, 0.002) 50, 100, 150, 200 ∼N (142, 842) ∼N (2.03, 0.029)
10 4000 N (50, 100) N (0.36, 0.002) 50, 100, 150, 200 ∼N (142, 842) ∼N (2.03, 0.029)
11 800 N (10, 4) N (0.73, 0.006) 100, 200, 300, 400 ∼N (39.4, 68.2) ∼N (3.89, 0.236)
12 8000 N (10, 4) N (0.73, 0.006) 100, 200, 300, 400 ∼N (39.4, 68.2) ∼N (3.89, 0.236)
13 800 N (100, 400) N (0.73, 0.006) 100, ..., 2000 ∼N (394, 6820) ∼N (3.89, 0.236)
14 8000 N (100, 400) N (0.73, 0.006) 100, ..., 2000 ∼N (394, 6820) ∼N (3.89, 0.236)
8 N (86.2, 575) N (0.35, 0.003) N (151, 21.8) N (1.94, 0.105) N (144, 523) N (2.02, 0.066)
9 N (81.7, 606) N (0.36, 0.0021) N (143, 818) N (2.01, 0.044) N (143, 838) N (2.01, 0.030)
10 N (82.5, 617) N (0.36, 0.0020) N (144, 883) N (2.02, 0.041) N (144, 864) N (2.02, 0.043)
11 N (−14.1, −39.6) N (0.615, 0.0040) N (41.5, −834) N (3.88, 0.265) N (41.3, −30.4) N (3.88, 0.250)
12 N (−12.8, 216) N (0.614, 0.0039) N (42.7, 356) N (3.88, 0.242) N (42.4, 270) N (3.88, 0.242)
13 N (82.2, 7490) N (0.606, 0.0018) N (427, −12840) N (3.85, 0.268) N (399, 9104) N (3.89, 0.256)
14 N(69.4, 1458) N(0.609, 0.0021) N (389, 25607) N (3.90, 0.236) N (390, 8864) N (3.90, 0.259)
6 CONCLUSIONS ACKNOWLEDGEMENTS
In this paper, theoretical and numerical efforts This study is supported by the State Key
are conducted to study the capacity of the lin- Development Program of Basic Research of China
ear regression method. The study demonstrates (Project No. 2012CB026104), a Chinese National
that: Key Technology R&D Program of the Ministry of
Science and Technology Project (2012BAH10B01),
1. A common misuse of linear regression is to
and the National Natural Science Foundation of
take the standard errors of the mean (SEM,
China (Grant No. 51109003).
Eqs. 5 and 6) as the “real” standard deviations
of regression parameters a, b. In fact, SEMs
are inversely proportional to the sample size REFERENCES
number, have nothing to do with the real var(a)
and var(b), but indicate the accuracy of the Chen, L.H. & Chen, Z.Y. 2007. Effect of nonlinear
mean value a, b (Eqs. 3 and 4). strength of rockfill on slope stability of high earth-rock
2. If a is relatively small comparing to the mean dam. Rock and Soil Mechanics 9(28): 1807–1810.
of independent variable, μX, Eq. 8 can improve Mann, C.J. 1987. Misuses of linear regression in earth
the accuracy of regression analysis comparing sciences. In W.B. Size (ed.), Use and Abuse of Statisti-
to Eq. 1, especially under the condition of small cal Methods in the Earth Sciences: 74–106. New York:
Oxford University Press.
sample size (smaller than 1000). Weisberg S. 2005. Applied linear regression. Hoboken:
3. In the statistical analysis of shear strength based John Wiley & Sons, Inc.
on triaxial test data, Eqs. 8 and 15 are recom- Yan, X. & Su, X.G. 2009. Linear regression analysis-
mended because they lead to mean and variance theory and computing. Singapore: World Scientific
values with higher accuracy. Publishing Co. Pte. Ltd.
144
ABSTRACT: Uncertainty quantification is the process of determining the effect of input uncertainties
on response metrics of interest. Many approximate methods have been developed so far for the purpose,
and among these methods the polynomial chaos expansion is considered a technique with strong math-
ematical basis and ability to produce functional representation of stochastic variability. The approach has
proven to be an efficient methodology to study several stochastic problems, considered the original form
where optimal convergence for orthogonal Hermite polynomials is only achieved for gaussian stochastic
processes. Regarding the developments in the stochastic response surface methodology and some heuris-
tic and optimisation concepts, results for a design example with correlated nonnormal random variables
are presented in the light of the quality of the approximate metamodels. Considered truncated full and
sparse polynomial chaos expansions, the efficiency and accuracy provided by different schemes of experi-
mental design are analysed and the convergence process is lastly discussed.
145
146
147
ρϕwϕw ρϕwγw ρϕwϕf ρϕwq ρϕwcf ρϕwγf Table 4. Summary description of results for the bearing
ργwϕw ργwγw ργwϕf ργwq ργwcf ργwγf capacity model.
ρϕfϕw ρϕfγw ρϕfϕf ρϕfq ρϕfcf ρϕfγf
ρqϕw ρqγw ρqϕf ρqq ρqcf ρqγf Reliability Reliability
ρcfϕw ρcfγw ρcfϕf ρcfq ρcfcf ρcfγf index index
βFORM βMCS
ργfϕw ργfγw ργfϕf ργfq ργfcf ργfγf
1.0 0.5 0.9 0.0 0.0 0.5 2.484 2.529
0.5 1.0 0.5 0.0 0.0 0.9 Output Output standard
0.9 0.5 1.0 0.0 0.0 0.5 mean deviation Output Ouput
0.0 0.0 0.0 1.0 0.0 0.0 (kN/m) (kN/m) skewness kurtosis
0.0 0.0 0.0 0.0 1.0 0.0
0.5 0.9 0.5 0.0 0.0 1.0 1267 899.1 1.802 6.391
ρ-coefficient of correlation. MCS and statistics results from 106 simulation steps.
148
149
REFERENCES
relative errors considered quite acceptable, but
the sparse polynomial chaos expansions may Armani, U., Khatir, Z., Khan, A., Toropov, V.V.,
present improved solutions. So, the fourth order Polynkin, A., Thompson, H., Kapur, N. & Noakes C.J.
presents higher accuracy for the first order reli- 2011. Control of physical consistency in metamodel
ability method, and the fifth order may improve building by genetic programming. In: Proc. of the
in general the results from Monte Carlo simula- Second International Conference on Soft Computing
tion. The statistical description of the output for Technology in Civil, Structural and Environmental
the bearing capacity metamodels is quite accurate Engineering, Chania, Crete, Greece: Civil-Comp Pro-
considered the reference model, with exception to ceedings 97, Paper 43, 1–18. Civil-Comp Press.
the case studies for a number of 150 design sam- Blatman, G. & Sudret, B. 2010a. An adaptive algorithm
to build up sparse polynomial chaos expansions for
pling points, where some differences in kurtosis stochastic finite element analysis. Probabilistic Engi-
are more important. neering Mechanics 25(2): 183–197.
Blatman, G. & Sudret, B. 2010b. Efficient computation
5 CONCLUSION of global sensitivity indices using sparse polynomial
chaos expansions. Reliability Engineering and System
Safety 95(11): 1216–1229.
In the field of abstraction, metamodelling is the Can, B. & Heavey, C. 2011. Comparison of experimental
construction of a cluster of concepts within a cer- designs for simulation-based symbolic regression
tain domain, highlighting properties of the model of manufacturing systems. Computers & Industrial
itself. Due to the computational complexity of Engineering 61(3): 447–462.
the analysis of current systems, metamodelling Haldar, A. & Mahadevan, S. 2000. Reliability assessment
has become an indispensable tool for exploring using stochastic finite element analysis. John Wiley &
the design space and performing sensitivity analy- Sons.
sis, particularly useful for the implementation of Hofwing, M., Strömberg N. & Tapankov M. 2011.
Optimal polynomial regression models by using a
numerical-based reliability techniques, to be used genetic algorithm. In: Proc. of the Second International
for design purposes and risk analysis. Therefore, Conference on Soft Computing Technology in Civil,
there is great interest in the development of tech- Structural and Environmental Engineering, Chania,
niques that facilitate the construction of approxi- Crete, Greece: Civil-Comp Proceedings 97, Paper 39,
mation models by using criteria of efficiency and 1–20. Civil-Comp Press.
accuracy. Huang, S.P., Liang B. & Phoon, K.K. 2009. Geotechnical
Regarding these issues, results for a design probabilistic analysis by collocation-based stochastic
example with correlated nonnormal random varia- response surface method: an Excel add-in implemen-
bles are presented. For the purpose it is considered tation. Georisk: Assessment and Management of Risk
for Engineered Systems and Geohazards 3(2): 75–86.
an inclined eccentric loading problem referred to Kiureghian, A.D. & Liu, P.L. 1986. Structural reliability
a bearing capacity nonlinear calculation model for under incomplete probability information. Journal of
a concrete gravity retaining structure, moreover a Engineering Mechanics 112(1): 85–104.
geotechnical reliability-based analysis problem. Li, D., Chen, Y., Lu, W. & Zhou, C. 2011. Stochastic
Considered the important nonlinear behaviour response surface method for reliability analysis of
of the geotechnical bearing capacity model, results rock slopes involving correlated non-normal variables.
show that the bearing capacity metamodels built Computers and Geotechnics 38(1): 58–68.
by multi-dimensional Hermite polynomial chaos
expansion may provide a general approximation
150
Y. Hu & J. Ching
Department of Civil Engineering, National Taiwan University, Taipei, Taiwan
ABSTRACT: A previous study by Ching et al. (2013) showed that the mobilized shear strength of
a spatially variable soil is governed by two factors: (a) line averaging effect along potential slip curves;
and (b) the emergent feature of the critical slip curve, i.e., the location of the critical slip curve depends
on the random field realization. In the study done by Ching et al. (2013), the latter factor plays a major
role as the slip curves are allowed to propagate freely without constraints. For scenarios where slip
curves are constrained, e.g., the slip curve in a vertical retaining wall typically passes through the toe,
it is not clear whether this latter factor will be fully suppressed. It is therefore the goal of this study to
understand the behavior of the mobilized shear strength for spatially variable soils when slip curves are
constrained.
153
⎛1 x2 H2 ⎞
Figure 2. Process of simulating a random sample of Pa Pa ax F ( x ) ma γ H 2 − τ LLA
f ⋅ , 0⎟
using the LEM.
x ⎝2 x ⎠
(2)
Finite Element Analysis (RFEA) results, it will The x position at which F(x) is maximized is
be clear that the actual effect of this assumption denoted by x*, which is the daylight position of the
is minimal. It is also assumed that the potential critical slip plane.
slip planes pass through the toe. This is barely an Active lateral forces also can be simulated by
assumption because the actual slip curves obtained RFEA. The first step of the RFEA is a geostatic
in RFEA always pass through the toe as well. For step that builds up the in situ stress field, and in the
clarity, only three of these planes are shown. second step a horizontal displacement is applied to
The sectioned τf profiles along the three potential the wall (away from the soil mass) to reach an active
slip planes are plotted in Figure 2b. The average of failure state. Then, the lateral force is calculated as
each profile (the dashed lines in Fig. 2b) is equal to the integration of the soil-wall contact forces in the
the line average for each section. Let us denote the FEA. Figure 3 shows a comparison between the
line average of τf(x,z) along a potential slip plane Pa samples simulated by RFEA and the LEM for
as τfLA. This line average can be directly simulated 100 cases. These 100 cases use wall heights in the
by the Fourier series method developed in Jha & range of [5 m, 20 m]. The random field realizations
Ching (2013) (Eq. 19 in Jha & Ching (2013)). for these 100 cases have fairly wide ranges of μ, σ,
Let us denote the line averages over the three δx, and δz (δx and δz are the horizontal and vertical
potential slip planes shown in Figure 2a by τfLA(xA), Scale Of Fluctuation (SOF), respectively). Each
τfLA(xB), and τfLA(xC). There are an infinite number case is analyzed by both RFEA and the LEM to
of potential slip planes. Each potential slip plane has determine the Pa values under the same random
a line average τfLA, and the resulting continuous line field realizations. It is clear that the Pa values simu-
average process τfLA(x) forms a One-Dimensional lated by both methods generally follow the 1:1 line.
(1D) random field. This 1D random field is illus- As mentioned earlier, although the assumption in
trated in Figure 2c. For each potential slip plane with
daylight position x, force equilibrium can be used to
derive the lateral force F required to support the fail-
ing soil mass. The force equilibrium of the wedge
formed by the potential slip plane can be written as
W N ⋅ cos( β ) + τ LLA
f ⋅ L p ⋅ sin( β )
F τ LA
f Lp ⋅ cos( β ) = N ⋅ sin( β )
1 x
H x ⋅γ = N ⋅ + τ LLA
f ⋅H
2 Lp
H
F τ LA
f x=N⋅
Lp
2
1 Lp 1 x2 + H 2
F γ H 2 − τ LLA
f ⋅ = γ H 2 − τ LLA
f ⋅ (1) Figure 3. Comparison between the Pa values simulated
2 x 2 x by RFEA and the LEM.
154
155
156
157
⎛ 1 ⎞
Pa H 2 − 2τ mf × H , 0 ⎟ (6)
⎝ 2 ⎠
158
159
T.M.H. Le
Department of Civil and Transport Engineering, Norwegian University of Science and Technology—NTNU,
Norway
ABSTRACT: Spatial variabilities of the soil shear strength parameters (friction angle and cohesion
coefficient) have long been recognised as an important factor influencing the reliability of geo-structures
including slopes. However, these two parameters are frequently considered separately in research studies
even though, in natural soils, both parameters are likely to vary simultaneously with existence of cross-
correlation between them. This study stochastically investigates the reliability of a slope constructed in
soil having spatially varying both friction angle and cohesion coefficient, and compares that with the
scenarios where each soil parameter varies individually. The finite element method is merged with the ran-
dom field theory to probabilistically evaluate the factor of safety and probability of failure of the slope via
Monte-Carlo simulations. A simple procedure to create perfect cross-correlation is discussed. The results
show that the variabilities of both friction angle and cohesion coefficient even without cross-correlation
can elevate the probability of failure relative to the cases where each of them varies individually. If a per-
fectly positive cross-correlation exists, the probability of failure increases significantly due to increasing
occurrences of local failures.
161
162
163
164
3.2 Factor of safety and probability of failure Figure 7. Variations of the μ(FoS) and COV(FoS)
with COV(Shear strength) for four scenarios
The variations of μ(FoS), COV(FoS) and pf (μ(tanφ′) = 0.375, μ(c′) = 15 kPa, COV(Shear
with the degree of variability of the strength strength) = COV(tanφ′) = COV(c′), θv = θ, θh = 2θ).
165
166
167
T. Namikawa
Shibaura Institute of Technology, Tokyo, Japan
ABSTRACT: The strength of cement-treated soil varies greatly owing to variability of the in-situ soil,
variability of mixing effectiveness, and other factors. This paper presents the numerical experiments that
investigate the effect of the spatial correlation on the behavior of cement-treated soil columns with dif-
ferent sizes. Finite element analyses were performed to simulate the unconfined compression behavior of
the columns with three different heights. The spatial distribution of the material properties was provided
by means of Monte Carlo simulation. The simulations were carried out with the samples in the presence
of different degrees of spatial autocorrelation. The numerical results provide good understanding of the
effects of the column height on the compression behavior of the cement-treated soil columns with spatial
variability.
169
170
171
172
173
174
M. Wasiul Bari
Department of Civil Engineering, Curtin University, Perth, WA, Australia
M.A. Shams
Housing and Building Research Center, Cairo, Egypt
M.A. Shahin
Department of Civil Engineering, Curtin University, Perth, WA, Australia
ABSTRACT: A stochastic approach that investigates the effects of soil spatial variability on stabili-
zation of soft clay via Prefabricated Vertical Drains (PVDs) is presented and discussed. The approach
integrates the local average subdivision of random field theory with the Monte Carlo finite element tech-
nique. A special feature of the current study is the investigation of impact of spatial variability of soil
permeability and volume compressibility in the smear zone as compared to that of the undisturbed zone,
in conjunction with uncoupled 3D finite element analysis. The results of this study indicate that the spatial
variability of soil properties has a significant impact on soil consolidation by PVDs; however, it is also
found that the spatial variability of soil properties in the smear zone has a dominating impact on soil
consolidation by PVDs over that of the undisturbed zone.
175
176
177
178
179
( )
ard deviation of the lognormally distributed t that σ ln t (U ) = lln ⎜1 + 2 s ⎟ (12)
achieves Us are calculated from the suite of 2000
s
⎜⎝ μt (U ) ⎟⎠
s
realizations using the following transformation
functions: Following the procedure set out above, prob-
nsim abilities of required time t to achieve Us that is less
1 i
μt (U ) =
s
nsim
∑ ti (U s ) (8) than or equal to ts can be estimated for any com-
bination of υ and θ, and the stochastic behavior
i i =1
of soil consolidation by PVDs can be investigated.
180
181
Figure 7. Effect of υu and υs on P[t(U90) ≤ ts(U90)] for θ ks θ mvs = 0.25 , indicate that μU is more or less
θu = θs = 0.5. independent of θu. Similar to θu, the influence of
θs (θ ku and θ mv are fixed at 0.25) on μU is also mar-
u
ginal as can be seen in Figure 4b.
182
ts(U90)] remains marginal. The effect of υs at fixed turbed and smear) that are equal to the υ and θ
values of υ ku = 50% and υ mvu = 10% is shown in of the smear zone does not significantly affect the
Figure 7b, which shows that varying the values of final results. Overall, the results obtained from
υs has a considerable impact on the estimated val- this research highlight valuable insights into the
ues of P[t(U90) ≤ ts(U90)]. At any certain consolida- impact of soil spatial variability on soil improve-
tion time, P[t(U90) ≤ ts(U90)] decreases significantly ment by PVDs and clearly demonstrate the
with the increase of υs. The overall observation potential of stochastic analyses in routine design
that can be derived from comparing the results in practice.
Figure 7 is that the effect of υs on P[t(U90) ≤ ts(U90)]
is dominant.
Figure 8 investigates the effects of θ on REFERENCES
P[t(U90) ≤ ts(U90)] at fixed values of υ ku υ ks = 50%
and υ mv υ mv = 10%. In Figure 8a, the influence Abuel-Naga, H.M., Pender, M.J. & Bergado, D.T. 2012.
of increasing θu on P[t(U90) ≤ ts(U90)] is shown
u s
Design curves of prefabricated vertical drains includ-
at θ ks θ mv = 0.25, and the results yield almost ing smear and transition zones effects. Geotextiles and
s
identical curves indicating that varying the values Geomembranes 32: 1–9.
Badaoui, M., Nour, A., Slimani, A. & Berrah, M.K.
of θu has little or no impact on the probabilistic 2007. Consolidation statistics investigation via thin
behavior of degree of consolidation. On the other layer method analysis. Transport in Porous Media
hand, the estimated P[t(U90) ≤ ts(U90)] for different 67(1): 69–91.
values of θs is plotted in Figure 8b at a fixed value Bari, M.W., Shahin, M.A. & Nikraz, H.R. 2012a. Effects
of θ ku θ mv = 0.25. It can be seen that unlike θu, of soil spatial variability on axisymmetric versus plane
θs has a considerable impact on the estimated val-
u
strain analyses of ground improvement by prefabri-
ues of P[t(U90) ≤ ts(U90)]. The comparison between cated vertical drains. International Journal of Geotech-
Figures 8a and b reveals that, the effect of θs on nical Engineering 6(2): 139–147.
P[t(U90) ≤ ts(U90)] is more significant than θu. It is Bari, M.W., Shahin, M.A. & Nikraz, H.R. 2012b.
Probabilistic analysis of soil consolidation via prefab-
interesting to know that the deterministic solution ricated vertical drains. International Journal of Geome-
yields P[t(U90) ≤ ts(U90)] < 50% for all combinations chanics, ASCE: (in press).
of values of υu, υs, θu, and θs, as can be seen in Beacher, G.B. & Christian, J.T. 2003. Reliability and
Figures 7 and 8. Statistics in Geotechnical Engineering, John Wiley &
Sons: Chichester, England.
Bo, M.W., Chu, J., Low, B.K. & Choa, V. 2003. Soil
5 CONCLUSIONS Improvement: Prefabricated Vertical Drain Techniques,
Thomson Learning: Singapore.
This paper has used the random field theory and Chu, J. 2004. Practical considerations for using vertical
drains in soil improvement projects. Geotextiles and
finite element modeling to investigate the rela- Geomembranes 22(1): 101–117.
tive significance of soil spatially variability in the Dagan, G. 1989. Flow and transport in porous media,
smear and undisturbed zones in soil improve- Springer: New York.
ment by prefabricated vertical drains. Both the Elkateb, T., Chalaturnyk, R. & Robertson, P.K. 2003.
coefficient of permeability, k, and coefficient of An overview of soil heterogeneity: Quantification and
volume compressibility, mv, were treated as inde- implications on geotechnical field problems. Canadian
pendent random variables and uncoupled 3D Geotechnical Journal 40(1): 1–15.
finite element analysis was applied. The effect of Fenton, G.A. & Griffiths, D.V. 2008. Risk assessment in
coefficient of variation, υ, and spatial correlation geotechnical engineering, Wiley: New York.
Fenton, G.A. & Vanmarcke, E.H. 1990. Simulation of
or scale of fluctuation, θ, of the undisturbed zone random fields via local average subdivision. Journal of
on the estimated mean and standard deviation Engineering Mechanics 116(8): 1733–1749.
of the degree of consolidation was found to be Freeze, R.A. 1977. Probabilistic one-dimensional consol-
marginal. On the other hand, the estimated statis- idation. Journal of Geotechnical Engineering Division
tics and probability associated with the degree of 103(GT7): 725–742.
183
184
Yong-Hua Cao
CCCC-Tianjin Port Engineering Institute Ltd., Tianjin, China
ABSTRACT: In the process of reliability analysis in geotechnical engineering, the variance reduction
function is the key parameter to transit the “point property” into “spatial average property”. In the paper,
improvement on the method to calculate the no correlation distance put forward by Yan Shu-wang and
Zhu Hong-xia was made, which can greatly minimize the computation efforts. The procedure of calcula-
tion of reduction function was analyzed, and general regularity of determining reduction function through
calculation of no correlation distance was obtained, which was for the first time put forward by us.
185
186
2 τ
Γ2 ( ) = ( bh
1) (7) ρ (τ ) = e −τ δu
cos ,
b2 h2 δu
Replace b in Eq. (7) with δu, and it can be δ u ⎛ δ u ⎞ 2 −h /δu h
Γ2 ( ) = − e sin (10)
obtained that: h ⎝ h ⎠ δu
2 2
Γ ( )= 2 2 ( bh
1) (9)
−e bh ⎡2ω bsin(
⎣ b siin((ω hh)) (ω
2
c ω h )⎤⎦
b2 ) cos( }
b h (11)
Figure 2. Reduction process of variance for correlation Figure 3. Reduction process of variance for correlation
function in exponential type. function in exponential type with cosine.
187
⎛δ u⎞ ⎡−0.495 − e − h / 0.505δ ⎛ h ⎞⎤
2
δu h
⎜⎝0.1sin 5.05δ − 0.99 cos 5.05δ ⎟⎠⎥
2
0.1 1/0.505b Γ ( )= + 0.095
⎝ h ⎠ ⎢⎣
u
5δu
h u u ⎦
⎛δ u⎞ ⎡−0.455 − e − h / 0.545δ ⎛ h ⎞⎤
2
δu h
⎜⎝0.3 sin 1.82δ − 0.91cos 1.82δ ⎟⎠⎥
2
0.3 1/0.545b Γ ( )= + 4.5δu 0.106
⎝ h ⎠ ⎢⎣
u
h u u ⎦
⎛δ u⎞ ⎡−0.375 − e − h / 0.625δ ⎛ h ⎞⎤
2
δu h
⎜⎝0.5 sin 1.25δ − 0.75 cos 1.25δ ⎟⎠⎥
2
0.5 1/0.625b Γ ( )= + 0.119
⎝ h ⎠ ⎢⎣
u
4δu
h u u ⎦
δu ⎛δ u⎞
2
⎡ − h / 0.745δ ⎛ h h ⎞⎤
⎜⎝0.7 sin 1.064δ − 0.51cos 1.0643δ ⎟⎠⎥
2
0.7 1/0.745b Γ ( )= + ⎢0.255 − e
u
3δu 0.189
h ⎝ h⎠ ⎣ u u ⎦
2
2 δu ⎛ δ u ⎞ e − h/δ sin h
1 1/b Γ ( )= − u
2δu 0.251
h ⎝ h⎠ δu
δu ⎛δ u⎞
2
⎡ − h / 1.345δ ⎛ h h ⎞⎤
⎜⎝1.3 sin 1.035δ + 0.69 cos 1.035δ ⎟⎠⎥
2
1.3 1/1.345b Γ ( )= + ⎢0.345 − e
u
3δu 0.215
h ⎝ h⎠ ⎣ u u ⎦
δu ⎛δ u⎞
2
⎡ − h//1.625δ ⎛ h h ⎞⎤
⎜⎝1.5 sin 1.083δ + 1.25 cos 1.083δ ⎟⎠⎥
2
1.5 1/1.625b Γ ( )= + ⎢0.625 − e
u
5δu 0.096
h ⎝ h⎠ ⎣ u u ⎦
⎛ δ u ⎞ ⎡0.945 − e − h / 1.945δ ⎛ h ⎞⎤
2
δu h
⎜⎝ 1.7 sin 1.144δ + 1.89 cos 1.144δ ⎟⎠ ⎥
2
1.7 1/1.945b Γ ( )= + 6.5δu 0.076
⎝ h ⎠ ⎢⎣
u
h u u ⎦
δu ⎛ δu ⎞
2
⎡ − h / 2.5δ ⎛ h h ⎞⎤
⎜ 2.sin 1.25δ + 3 cos 1.25δ ⎠⎟ ⎥
2
2 1/2.5b Γ ( )= + ⎢1.5 − e
u
9.5δu 0.057
h ⎝ h⎠ ⎣ ⎝ u u ⎦
188
δu ⎛ δu ⎞ 2 ⎡
Γ2 ( ) = + α 1 − e − h / α 2δ u
h ⎝ h⎠ ⎣
⎛ h h ⎞⎤
× ⎜ α 3 sin + α 5 cos
α 4δ u ⎟⎠ ⎥⎦
(12)
⎝ α 4δ u
4 CASE STUDY
2 bh
Exponential type ρ(τ) = e−b|τ| 2 2
( bh e 1) (where, h = L ) 0.095 1
b h
1
1. Type of exponential with ⎡
bh − e
( bh )
si bh ⎤⎦ (where, h = L ) 0.251 1
b h ⎣
2 2
Cosine ρ(τ) = e−b|τ| ⋅ cos(bτ)
δu/h* 1
δu⎛δ ⎞
2
⎡α − e − h / α δ ⎛ h h ⎞⎤
2. Type of exponential with + u ⎢1 ⎜⎝α 3 sin α δ + α 5 cos α δ ⎟⎠⎥ (where, ω/b =α,
2 u
189
Table 4. Vertical correlation distance of Tianjin Port. 4.2 Variance reduction function
Correlation distance (m) Figure 7 describes laws for the reduction function
Cone penetration with ω/b. Results of this area can be employed to
test number Silt Silt clay Silty clay prove whether the law is ubiquitous. Reduction
functions of variance corresponding to the ratio
1 0.139 0.533 0.400 ω/b are plotted in Figure 7.
2 0.330 1.730 0.358 It can be known from Figure 7 that the law men-
3 0.209 1.220 0.431 tioned above is ubiquitous. That is with increas-
4 0.182 1.067 0.400 ing of ω/b, the no correlation distance decreases
5 0.190 0.273 0.483 when the ratio is smaller than 1, and the reduc-
6 0.412 0.800 0.431 tion coefficient increases; conversely, with increas-
7 0.067 1.077 0.297 ing of ω/b, the no correlation distance increases
8 0.270 1.239 0.325 when the ratio is larger than 1, and the reduction
9 0.984 1.239 0.308 coefficient decreases. The no correlation distance
10 0.175 0.400 0.492
is the largest and the reduction coefficient is the
11 0.543 0.554 0.335
smallest when the ratio ω/b is 1. To conservatively
12 0.176 1.103 0.224
consider in practical project, the smallest ratio ω/b
13 0.649 0.690 0.163
will be taken to determine the reduction function
14 0.130 0.825 0.400
of the variance. In addition, the value of variance
15 0.341 0.226 0.400
16 0.113 0.588 0.431
17 0.701 0.865 0.400
18 0.181 0.550 0.400
19 0.452 1.000 0.457
20 0.511 1.000 0.400
21 0.923 0.650 0.398
22 0.195 0.550 0.500
23 0.320 0.667 0.276
24 0.088 0.738 0.328
25 0.202 1.077 0.400
26 0.400 0.500 0.363
27 0.667 0.320 0.363
28 0.690 0.550 0.328
29 0.546 0.054 0.384
30 0.267 0.780 0.328
31 0.495 0.832 0.303
32 0.100 0.780 0.328 Figure 7. The change law of variance reduction
33 0.100 1.429 0.260 function.
34 0.100 0.667 0.355
35 0.100 0.920 0.328
36 0.100 1.107 0.328 Table 5. Variance reduction function of typical layers
in Tianjin Port.
37 0.098 1.176 0.303
38 0.183 1.440 0.270 No Reduction
39 0.175 0.976 0.303 Name of ω/b correlation function of
40 0.180 0.148 0.270 the soil (min) distance the variance
Minimum value 0.040 0.054 0.163
Maximum value 0.984 1.730 0.500 Silt 0.17 4.75 δu 0.374
Average 0.317 0.808 0.356 Silt clay 0.21 4.75 δu 0.312
Standard deviation 0.237 0.368 0.072 Silty clay 0.38 4.50 δu 0.378
190
191
ABSTRACT: Natural slopes are highly heterogeneous. When inverse segregation takes place in a land-
slide soil deposit, the coefficient of permeability of the soil statistically decreases from the surface to the
bottom of the deposit. Such material spatial variability can be characterized using non-stationary random
fields. The objectives of this paper are to develop algorithms to generate non-stationary random fields for
a slope with the mean coefficient of saturated permeability decreasing with depth and to carry out analy-
sis of infiltration in an anisotropic non-stationary slope. Basic anisotropic stationary random fields with
orthogonal major and minor scales of fluctuation are considered first. Non-stationary random fields are
then formulated by adding a trend to the stationary random fields. The Cholesky decomposition method
is applied to produce correlated random numbers following an exponential correlation structure. After
generating the anisotropically heterogeneous non-stationary random fields, the steady-state pore-water
pressures in each random field are obtained by solving the Richards equation. The effects of coefficient
of variation of the saturated permeability and the ratio between the major and minor scales of fluctuation
on the pressure bound are investigated.
193
194
Y( x, y ) x, y ) × f ( y )
X( x, (8)
195
∂q ∂q
+ =0 (12)
∂x ∂y
196
Case Coefficient
no. of variation Scale of fluctuation
1 1.0 θh = 10 m, θv = 0.1 m
2 0.5 θh = 10 m, θv = 0.1 m
3 0.1 θh = 10 m, θv = 0.1 m
4 1.0 θh = 10 m, θv = 1.0 m
5 0.5 θh = 10 m, θv = 1.0 m
6 1.0 θh = 10 m, θv = 0.5 m
197
198
ABSTRACT: Probabilistic Seismic Hazard Analysis (PSHA) has been widely used to estimate the
ground motion intensity that should be accounted for in structure design. Here, an in-depth PSHA study
for Kaohsiung, the metropolitan area in South Taiwan, is carried out. On the other hand, soil liquefaction
and subsequent structural damage are usually reported in a catastrophic earthquake event, such as numer-
ous significant incidents of liquefaction in inland alluvial areas during the 1999 Chi-Chi earthquake.
A quantitative framework to evaluate earthquake-induced soil liquefaction potential is proposed herein,
including earthquake hazard assessment and geotechnical liquefaction evaluation. This framework is fur-
ther demonstrated with a case study for a thoroughly-investigated site in Kaohsiung. Results show that
the liquefaction potential evaluated with existing methods is to some extent dependent on the earthquake
hazard in terms of the occurrence probability of large-size earthquakes at the study site. Ground improve-
ment is suggested considering severe consequences of ground failure.
201
202
Pr( y* λ ,,t)
t) 1 e− λ t (4)
NS NM N D
λ( y*)
*) ∑ i ∑ ∑ Pr ⎡⎣Y y* | m j , d k ⎤⎦
i =1 j =1 k =1
× Pr ⎡⎣M m j ⎤⎦ × Pr [ D dk ] (3)
203
Maximum
Source zone a-value b-value magnitude
Figure 6. PGA exceedance probabilities within four
A 2.670 0.828 6.5 time spans at the study site, given the annual rate of seis-
B 3.790 1.074 6.5 mic hazard calculated from PSHA (Fig. 5).
C 3.579 1.014 7.1
D 4.802 1.334 7.3
E 4.420 1.087 7.3 Given PGA > 0.31 g with an annual rate of 0.0022
F 4.720 1.117 7.3 through PSHA, the occurrence probability for
G 3.580 0.876 6.5 such an event is about 20% as far as a time span
H 3.071 0.756 7.3 of 100 years is concerned. On the other hand, the
I 4.450 1.137 6.5 occurrence probability is reduced to 2% when the
J 3.210 0.979 6.5 time of interest is decreased to 10 years.
K 4.490 1.101 6.5
L 3.510 0.918 7.3
M 3.307 0.803 6.5 4 LIQUEFACTION POTENTIAL
N 5.345 1.134 8.0 ASSESSED OVER TIME
O 5.440 1.107 8.3
P 4.920 1.058 7.8 Accordingly, at this study site we found that L = 15
Q 4.170 0.898 7.8 was on the basis of PGA > 0.31 g, and we have
R 3.640 0.913 7.8 also characterized a 20% occurrence probability
S 5.104 1.036 8.0 for PGA > 0.31 g in 100 years, we can therefore
T 4.110 0.900 7.8 combine the two pieces of underlying information
to evaluate time-variant soil liquefaction with a
new index referred to as “soil liquefaction index
per
p~
unit time,” denoted as LT t* . In this instance,
LT 100 yrs is equal to 3.0 per 100 years given a
unique scenario PGA > 0.31 g. Because there are
infinite scenarios other than PGA > 0.31 g that can
cause different levels of soil liquefaction p
potential,
~
Figure 7 shows the relationships between LT t* and
PGA given four periods of time, basically a combi-
nation of Figures 3 and 6. As a result, the integra-
tion of the curve is the resulting soil liquefaction
index per unit time covering any possible scenario
with
~
a small PGA to a large PGA. For example,
LT 100 yrs is about 1.18 at the study site.
The relationship between LT t* and t* at the
site with repetitive calculations from t* = 1 year
to t* = 1,000 years is nearly a perfect polynomial
function obtained with regression analysis on the
Figure 5. PSHA hazard curves at the study site with the 1,000 points:
use of four PGA attenuation relationships to account for ~
such an epistemic uncertainty. T t 0.000003 t*2 + 0.0112 t* + 0.0887 (5)
204
205
206
A.D. Mackay
Nishimatsu Construction Company Ltd., Hong Kong, China
Formerly SMEC Asia Ltd., Hong Kong, China
ABSTRACT: The Pakistan section of the Karokorum Highway runs from Attobad in the south, crossing
Pakistan’s north-western frontier mountainous region to the Khunjerab Mountain Pass, in the north.
This section of the highway is over 800 long and rises to a maximum elevation of about 4,700 m above
sea level, making it the world’s highest paved international highway. Along the majority its length the
alignment follows the route of the Indus River, situated between the Himalaya and Karakorum mountain
ranges. During 8 October 2005 a severe earthquake struck the region triggering landslides impacting the
highway and destabilizing the natural terrain and cut slopes. This paper provides an overview of the insta-
bility of the region with particular attention specific inspections along the highway, carried out prior to
2010 impacted by the 2005 earthquake event. Details of the landslide assessment approach, particularly
detailed site reconnaissance survey, are summarised for specific locations along the highway running
through the Jhelan Valley. The findings of the landslide characterization and stabilization measures based
on the landslide assessment are presented. The proposed stabilization measures included river erosion
protection, debris removal, slope cut back and slope support.
1 INTRODUCTION
207
208
3 GROUND CONDITIONS
209
210
5 MITIGATION MEASURES
Table 2. Proposed rock slope cut design, exceeding 7–10 m vertical height.
Medium jointed rock slightly 1:3 Locally—30 mm; Locally—50 mm dia. 12/2.6–3
weathered as required 6–9 m long slotted
PVC
Intensively jointed rock, 1:2 50 mm; N/A; 25 mm 4 m c/c, 50 mm dia., 12/3
unfavourable dia., 4–6 m length 6–9 m long slotted
discontinuities PVC
Soft rock, weak, highly 1:1 70 mm; N/A; 25 mm 4 m c/c, 50 mm dia., 7/3
fractured and weathered, dia., 6 m length 9–12 m long slotted
shale, partly disintegrated PVC
N/A—Not applicable.
211
N/A—Not applicable.
6 CONCLUSIONS
Figure 10. Slope stabilization measures with slope cut
back. The detailed inspections assessing the ground
conditions along specific landslide locations along
the alignment, considering influences affecting the
rock dowel and/or soil nail and retaining walls stability, such as the tectonic effects triggering
with sub-surface drainage measures, cut back in landslides, rainfall, river erosion and the exisitngf
combination with slope support measures. River g steep natural gradients, were considered applica-
erosion protection measures are presented in ble in assessing preliminary landslide mitigation
Plate 5. proposals.
212
213
A.D. Mackay
Nishimatsu Construction Company Ltd., Hong Kong, China
Formerly SMEC Asia Ltd., Hong Kong, China
ABSTRACT: East Timor has an extensive road network often prone to land-sliding; typically this
causes disruption to its transportation system which is the island’s life line. The instability impacting
the roads is mainly a combination of steep terrain; high annual precipitation, typical of a sub-tropical
climates, and complex ground conditions affected by the tectonic uplift and deep weathering zones com-
mon throughout the island. Other influences affecting the susceptibility include erosion sensitivity and
roads being placed on poorly prepared and un-compacted ground. Remedial measures against landslides
have considered optimization of road alignment corridors; the avoidance of areas of concern, were pos-
sible, and a reduction in the reliance of traditional slope stabilization measures. Traditional stabiliza-
tion measures have included stone pitching and gabion boxes, which have proven insufficient to support
modern transportation frequency and load impacts. Up to date landslide stability evaluation has allowed
assessment of landslide susceptibility for different ground conditions. Recommendations for suitable,
regular maintenance with particular attention to effective road drainage networks have been based on
this. Notwithstanding effective resource, with suitable training and allocation, is required to carry out
this maintenance, and allocation using suitable construction materials, sourced from suitable locations as
needed. This paper provides an overview of the road network throughout East Timor and factors affect-
ing the landslide susceptibility, such as precipitation, terrain and the geological conditions. Details of
some of the landslides, representing different types of susceptibility, impacting the roads are presented,
with suggested slope stability mitigation measures for each representative landslide.
1 INTRODUCTION 2 SETTING
East Timor’s road network has an extensive and 2.1 Geographical features impacting instability
large proportion of paved, surfaced roads com-
Timor Island is located south of the Banda Arc
pared with many other countries with similar capi-
Archipelago, which extends eastwards from Bali
tal expenditure. Notwithstanding there has been
Island. East Timor occupies the north eastern
a recent increase in the transportation frequency
portion of Timor Island. Refer to Figure 1 for the
and loading, which has impacted and/or been
location.
effected by the instability occurring along the road
network. The instability is typically a result of a
combination of poor road base preparation, com-
plex geology, including tectonic uplifting and deep
weathering zones, and high precipitation. To miti-
gate this instability initiatives have been taken to
carry out remedial measures and implement more
efficient road maintenance. These measures have
included more robust cut and fill support along
the road alignments, with particular attention to
effective surface and sub-surface drainage systems,
optimizing road corridors and avoiding zones sus-
ceptible to landslides were possible. This paper
provides an overview of the road network, geog-
raphy, climate, terrain and geology impacting the
landslide susceptibility. Details of proposed stabi-
lization initiatives for different and representative
landslides encountered are presented. Figure 1. East Timor location (Google Earth, 2011).
215
216
217
Geological Facies/lithology
Formation age (remarks)
218
Zone Description
Zone Description
Plate 6. Mudflow, away from road, Cablaci Limestone. The zoning generally corresponded with the
source geological conditions, the weathering inten-
sity, natural slope gradient and precipitation. An
4 INTERPRETATION approximate zoning corresponding with the source
geological formation encountered is summarized
The stability affecting the roads could generally in Table 4, as Zone 1 comprises flat recently depos-
be zoned based on the general scale of ited material, such as alluvium, this has not been
instability associated with ground conditions, included.
were information is available. A generic zoning Generally instability increases with younger,
presenting basic landslide susceptibility, is uplifted geological formations; hence the inclu-
presented in Table 3. sion of the Viqueque, and Barbanoro Formations
219
5 MITIGATION MEASURES
ACKNOWLEDGEMENTS
REFERENCES
220
S. Jung
Dr. Jung and Lang Engineering, Germany
R. Becker
Holmalux, Luxembourg
ABSTRACT: In January 2011 a landslide occurred in the German village of Kirf, which is between
the town Saarburg and the border of Luxembourg. This slope failure caused a destruction of a federal
main road over 50 meters. Stability calculations were made by a geotechnical consulting company,
commissioned from the Federal Road Administration, in order to design the repair works and also to find
out what really had happened. These calculations showed remarkably high safety factors. Nevertheless
the consulting company blamed the families of the four houses down the slope for the failure, because
they had made gabion walls in the toe of the slope. A research team of the University of Luxembourg
discovered that this hilly area originally was called “On the slide” and has never been stabile according to
their own calculations. So, the real question is why so many apparently uncorrelated mistakes were made
in a single case.
221
222
223
224
J.P. Wang
Department of Civil and Environmental Engineering, Hong Kong University of Science and Technology,
Kowloon, Hong Kong, China
ABSTRACT: Engineering risk is referred to as occurrence probability multiplied with the consequence
when the event occurs. Few engineering decisions are now risk-based, using an alternative with the lowest
risk as the solution to a problem. Under such a risk framework, this study develops a classification system
for active faults. In addition to population at risk, this classification also considers the earthquake mag-
nitudes and earthquake energy release. Applying such a classification to six active faults in Taiwan, we
suggested that the Sanchiao Fault and Changhua Fault should be Category III active faults, the highest
level in this classification system.
Engineering risk is referred to occurrence The risk score of an active fault in this study is
probability multiplied with the consequence when referred to the relative level of risk among the
the event occurs (Wang et al., 2012a). As for risk- faults analyzed. Note that the original score is
based decision making, the decision with the lowest normalized, so that the resulting risk score is
risk will be adopted. Few engineering decisions are dimensionless. In other words, after normalizing
now risk-based, including site selections for criti- the units of the underlying parameters should not
cal structures (Wang et al., 2012b), and earthquake affect the classification.
early warning (Wang et al., 2012a). With the concept of risk, the risk score (RS) of
Understandably, probability is used to quan- an active fault in this study can be first expressed
tify the level of uncertainty or randomness. For as follows:
example, the probability in a random process that
“Point One” occurs while tossing a fair dice is 1/6, RS P (M > m*) CS (1)
which reflects the randomness of this stochastic
process. Since the earthquake is hardly predictable where Pr( *) is the probability of magnitude
given our limited, imperfect understanding (Geller exceedance for M > m* when an earthquake occurs
et al., 1997), probabilistic analyses are gradu- at an active fault; CS denotes the consequence
ally applied to evaluating earthquake potentials when such an event occurs. This expression is the
in terms of probability. For example, a couple underlying governing equation of this study.
of probabilistic studies (Wang et al., 2013a) were
performed and suggested that indeed the earth-
2.1 Earthquake energy release
quake’s temporal occurrence follows the statistical
Poisson model that has been commonly used in Earthquake energy release exponentially increases
earthquake analyses (Wang et al., 2012c, 2013b). with magnitude. For example, an M 7.0 earthquake
Another instance is the analysis evaluating large- can release energy 32 times as large as an M 6.0
earthquake probabilities from studying the pool event (Keller, 1996). Figure 1 shows the energy
of annual maximum earthquakes in Taiwan since release ratio on the basis of M 6.0 earthquakes. It
1900 (Wang et al., 2011). is around a 1,000-time difference in energy release
This paper presents a risk-based classicification between M 6.0 and M 8.0 earthquakes.
system for active faults. In addition to Population
At Risk (PAR), the best-estimate earthquake
2.2 Earthquake energy release and consequence
magnitude and earthquake energy release are
accounted for in the classicification developed. Relatively speaking, the relationship between earth-
The method was then applied to six active faults in quake energy release and the consequence caused by
Taiwan, which is considered one of the earthquake earthquake is more linear than that between earth-
“hot zones” in the world. quake magnitude and the consequence. Therefore,
225
226
227
228
ABSTRACT: From observation and experience, earthquake recurrence induced by the same fault should
be influenced by “memory”, or such a stochastic process should be non-stationary rather than being
stationary. For instance, the same fault triggering the recent Tohoku earthquake in Japan is less likely to
trigger another one in the coming few years, but the probability will be increasing with time. However, the
commonly used Poisson process considers the earthquake as a stationary stochastic process, resulting in the
identical probability regardless whether the next earthquake occurs in 2020∼2030 or 2220∼2230. As a result,
based on the mechanism of faulting, this paper introduces a new analytical model to properly take the earth-
quake memory effect into account. A case study in Central Taiwan is given, which is considered to provide a
more realistic result compared to that from the stationary stochastic process suggested by Poisson model.
229
3 EARTHQUAKE PROBABILITY
INDUCED BY THE MEISHAN FAULT
WITH THE POISSON PROCESS
=
( × t ) e − λ × Δt
n
(1)
stresses in horizontal and vertical directions, respectively.
n!
stress accumulation (Δτ ) and stress accumulation
where n = the number of the event; t1 = the starting (τ t) at time t on the fault plane can be formulated
date; Δt = time interval; and λ = the mean annual as follows:
rate, respectively.
Given the occurrence probability of at least one T t
exceedance in a specific time window can be esti- τt ∫ T = 0 Δτ (3)
mated as follows:
Applying the boundary condition derived
Pr( ; t1 ≤ T t1 + Δt; λ ) beforehand, the relationship between the shear
stress and the strength at return period R can be
= 1 − e − λ × Δt (2) expressed as follows:
4.1 Mechanism of faulting and faulting 4.2 Model verification on a deterministic basis
probability
On a deterministic basis, the earthquake probabil-
Figure 2 shows a systematic diagram of the under- ity is zero before the return period is due. Since the
lying mechanism of faulting. The model is dictated stress accumulation in every 162 years in this case is
by the law that when the shear stress exceeds the equal to the strength S, the unit stress accumulation
strength of the rock along the fault, the faulting is per year is S /162. Therefore, the stress accumula-
initiated and leads to an earthquake (Reid, 1911). tion from 1906 to 2020 is equal to ( ) /162,
Therefore, the first boundary condition in this and ( ) /162 till 2040. As a result, the earth-
model is that the accumulated shear stress equals the quake probability is zero in 2020∼2040 because
strength of the rock when the return period is due. at any moment during this time period the stress
Given the ground subject to a stress accumula- is less than the strength. Applying the calcula-
tion with time, the relationship between transient tion to evaluating the earthquake probability in
230
E ⎡∫ Δτ ⎤
T R
E [τ R ] = S (6)
⎣ TT= 0 ⎦
⎡t T ⎤
E ⎡ ∫ Δτ ⎤
T t R
E [τ t ] = E ⎢ ∫ Δτ ⎥
⎣T o ⎦ ⎣R T o ⎦
t
= ×S (7)
R
VAR ⎡ ∫ Δτ ⎤⎥ = VAR[τ R ] σ
T R *
⎣⎢ T = 0 ⎦
⎡t T R ⎤
⇒ VAR ⎡ ∫ Δτ ⎤ = VAR ⎢ ∫
T t
Δτ ⎥
⎣⎢ T = 0 ⎦⎥ ⎣ R T =0 ⎦
2
⎛t⎞
= ×σ* (8) Figure 3. Earthquake probability with respect to differ-
⎝ R⎠ ent coefficient of variation.
231
232
233
C. Zhao, W.Z. Zhang, J.H. Zhang, K.Z. Wang & Z.S. Xiong
State Key Laboratory of Hydroscience and Engineering, Tsinghua University, Beijing, China
ABSTRACT: Gas hydrate is ice-like compound consisting of mostly methane and water. It is a
promising energy source but very unstable and easy to dissociate when it is disturbed. If gas hydrate in
the sediments under the sea dissociates, it may trigger a submarine slide and render high risk to human
beings and facilities in the sea. In this paper, modeling techniques were developed to simulate changes
of strength and pore pressure in sediments caused by the dissociation of gas hydrates in a geotechnical
centrifuge in order to assess the risk potential associated with submarine slides. The soft toothpaste,
Kaolin paste and Kaolin paste underlying sand were used to simulate sediments containing gas hydrate.
Centrifuge modeling tests were carried out on submarine landslides using the materials. The effectiveness
of the methods is discussed.
235
236
237
238
G.X. Chen
School of Civil Engineering, Tongji University, Shanghai, China
School of Civil Engineering and Architecture, Henan University of Technology, Zhengzhou, China
ABSTRACT: Subway construction schedule management is one of the major tasks for the government
departments, owner and construction units. As a new schedule management method, critical chain method
is effective at shortening the project duration, optimizing the resources scheduling and economizing the
project investment. Based on the analysis of basic theory and assumption of the critical chain, and com-
bined with the characters of subway station construction, steps applying critical chain to subway project
were illustrated in this paper, and methods to calculate the buffer size and activity duration were also
introduced. Finally, the critical chain method is systematically used in one subway station of Zhengzhou
Metro 1 route and the performance of critical chain method is satisfactory.
241
242
243
A1000 A Building 23 22
demolition
A1010 B Pipeline 30 30
moving
Figure 1. Activity-on-node network of this project. A1020 C Maintenance 175 175
construction
A1030 D Dewatering 10 10
construction
A1040 E Foundation pit 25 25
construction
A1050 F Main 140 140
construction
A1060 G Gateway 70 70
construction
A1070 H Affiliated 95 95
construction
Figure 2. Initial scheduling.
4.2 Activity duration and buffer size calculation relocation construction can be carried out together
with building demolition. Then, adjust their logical
According to the critical chain by Goldratt, the relations as SS; adjust logical relations of pipeline
critical chain activity duration is half the amount relocation and maintenance construction as FS;
of planning duration, and the buffer size is half the adjust logical relations of dewatering construction
amount of safety time related. and foundation pit construction as FS-6; adjust
logical relations of main construction and founda-
1 tion pit construction as SS+24, and other logical
d′ d (1)
2 relations remain invariant. The logical relationship
of the activities and the critical chain activity dura-
n
⎛1 ⎞ tion after adjustment are shown in Table 3.
Buffer size = ∑ Δti (2)
⎝
i =1 2
⎠
4.4 Determining the critical chain
where: d ′ is critical chain activity duration, Δti and the buffer
stands for safety time, n is the number of the activi-
Take the necessary time of the activities as the
ties. Then the critical chain activity duration and
critical chain duration, then get PB and FB of the
buffer size are shown in Table 2.
project using the cut & past method and insert it
into the corresponding chain, which makes rea-
4.3 Adjusting the initial scheduling
sonable use of the safety time and ensures that the
After analyzing the actual situation of the project project can be implemented in accordance with
and optimizing the collective resources and the logi- the schedule, and the project duration is shortened
cal relationship of the activities, we find that pipeline as well. According to the calculation result, the
244
ACKNOWLEDGMENT
245
ABSTRACT: Rainfall-induced slope failures are one of the most frequent hazards in hilly terrains.
This paper presents a method that adopts a distributed cell model for predicting rainfall-induced slope
failures in a large area and the runout distances of detached materials. The model consists of five compo-
nents; namely, a digital terrain model, a spatial rainfall distribution model, a rainfall infiltration analysis
model, a slope reliability analysis model and a runout distance prediction model. A cluster of unstable
cells located together are viewed as an entity and the size effect of the combined detached materials is
considered. The materials move down the slopes along the steepest paths on the digital elevation map. The
method predicts the locations of rainfall-induced slope failures and their depositional zones reasonably
well. The material movement traces assessed by this model also offer information for the assessment of
landslide and rockfall risks. The presented method shows promise for use as a module in a real-time warn-
ing system for rainfall-induced slope failures.
2 STUDY AREA
247
248
249
REFERENCES
250
ABSTRACT: The failure of a dike system, as a characteristic of long linear systems, is ruled by its
weakest link. Meanwhile the most serious hazard in a dike section is its breaching, which can develop
due to different failure mechanisms. The objective of this paper is to identify the most common fail-
ure mechanisms and evaluate them from the point of view of risk assessment. As illustrated via four
case studies, the most frequent problem in the design and risk assessment is the inappropriate address-
ing of potential threats. The results of risk analysis are not reliable if relevant failure mechanisms are
neglected, the structural behaviour of the dike-section is not understood or inadequate maintenance is
provided.
251
Investigation period
(1564–2010) (1134–2006)
Not known 58.1% Erosion of inner slope prot. & crest 67%
Overtopping 32.4% Ice drift 11%
Human activity 2.3% Erosion or instability of 6%
Slope failure 2.0% outer slope prot.
Subsoil failure 1.8% Sliding in. slope 5%
Other known 1.2% External (human 4%
Structural failure 1.3% and animal)
Wave erosion 0.6% Sliding out. slope 3%
Liq. of shore line 2%
Piping 1%
Micro instability 0.5%
Horizontal shear 0.5%
Total 2858 cases 1735 cases
structures). Based on their results, the most initiations are hardly observable without a moni-
frequently observed cause of dike breaching was toring system. More likely they are only noticed
the external erosion due to overtopping. Erosion in the developed phase. By then the opportunity
of the outer slope, distresses caused by human of implementing rescue measures is limited. In
or animal activities and sliding of the inner slope Figure 2 the most common failure modes related
during flood event are the most serious threats in to external and internal erosion are summarized.
addition to overtopping. Moreover potential factors that may accelerate
erosion are highlighted.
Additionally non-erosion related failure modes
2.2 Most important failure modes
like slope failures must be evaluated both under
In the historical overview it was revealed that the flood conditions (the down-stream side) and under
most probable failure scenario is the overtopping rapid draw-down conditions (the upstream side).
of dike section. However the external erosion Local slope failures caused by seepage in the inner-
might not be governed by the elevation difference slope may be referred to as micro-instabilities in
between the flood water level and the dike crest, but the literature.
by waves. When the erosion starts in the upstream Other failures might be governed by settlement,
side (e.g. the Upper Jones Tract levee breach), the horizontal shearing (sliding), earthquakes, drift-
failure of the riverside slope will become the most ing ice or collision with floating objects. Non-
likely cause of breaching. erosion related failure modes are summarised in
Although the most probable failure mechanism Figure 3.
is the external erosion of the section it must be
stated that, from risk management point of view,
2.3 Failure modes vs. risk
their assessment may not be the most challenging
one. Nowadays flood water levels can be estimated The first issue in the risk management is whether a
with acceptable accuracy; therefore the overtop- threat is identified or not. One of the most common
ping of the section is quite predictable. Moreover, reasons for dike breaching is that a failure mode is
the observation of the external erosion type of fail- considered as “unlikely”, hence not included in the
ures is quite straightforward. Thus the warning of design or risk analysis. Once the failure unexpect-
the population at risk and their evacuation from edly develops, adequate countermeasures may not
the flood zone seem to be more manageable than be well prepared.
in the case of a less ductile failure. Other possible failure scenarios occur when the
The failure modes of dike sections originated potential failure modes are identified but the dike
from internal erosion are more deceitful since their construction does not meet the design criteria.
252
253
254
255
256
ABSTRACT: Pressure transducers are increasingly being used within soil mass or at soil-structure inter-
face to monitor stresses at the point of installation. Calibration of these transducers provides unique rela-
tionship between applied pressure and measured stress. Routinely, calibration is performed by applying
fluid pressure on diaphragm of transducer, however interpretation of data using fluid calibration results
lead to measurement errors when transducer is used in soil for stress measurement. This paper presents
comparison of calibration results of pressure transducer using fluid and soil calibration techniques. The
bias factor in the earth pressure measurement using fluid calibration results are evaluated, which can be
used with fluid calibration results to precisely arrive at the values of earth pressures on the retaining wall,
in the absence of in-soil calibration test results.
257
3 EVALUATION OF CALIBRATION
FACTORS
258
each loading were obtained once the fluid pressure Figure 3. Particle size distribution of sand.
had stabilized, typically in 30 seconds.
The relationship between applied pressure and
acquired data (measured in terms of micro strain) Table 2. Physical and mechanical properties of sand.
and measurement non-linearity for both the trans-
Property Value
ducers was obtained and presented in Table 1.
Gs 2.65
3.2 In-soil calibration Cu 1.42
Cc 0.93
Effect of sand layer thickness on the performance γdmin (kN/m3) 14.58 (ASTM D4254-00)
of Type 2 transducer is observed. Same modi- γdmax(kN/m3) 17.10 (Pluviator)
fied triaxial setup with additional plastic tube of Cohesion, kPa (c′)* 0
2.5 mm thickness and 100 mm Internal Diameter Friction angle (φ′)* 39°
(ID) placed tightly on the brass pedestal was
used with the transducer fixed flush on it. In the *Obtained from direct shear test at 68% relative density.
present study silica sand (Indian Standard sand
of Grade III) was used. The typical particle size
distribution curve of the sand used in the study (Type 1 and Type 2) were compared, discussed in
is shown in Figure 3. Some of the physical and the following section.
mechanical properties of the sand used in the
present study are presented in Table 2.
Considering diameter of EPC (DEPC) as refer- 4 RESULTS AND DISCUSSION
ence, sand layer thicknesses of 10 mm (0.25DEPC),
20 mm (0.5DEPC), 40 mm (DEPC), 60 mm (1.5DEPC) Calibration performances of Type 1 and Type
and 100 mm (2.5DEPC) were placed above the EPC 2 transducers in the range of 0–180 kPa for first
to obtain relation between applied pressure and loading-unloading cycle were obtained. Typical
measured strain. Greased polyethylene sheets of fluid calibration results are presented in Figure 4.
60 μm thickness were pasted to inner surface of The fluid calibration results of first loading-un-
plastic tube in order to reduce friction between loading cycle indicated maximum non-linearity of
sand particles and plastic tube surface. A rigid 0.6% of FS (full scale) and 2.64% of FS for Type 1
wooden block followed by thick rubber pad was and Type 2 transducers, respectively.
placed on sand layer for proper transfer and uni- The calibration factors obtained from fluid cali-
form distribution of pressure between load cell bration of these two pressure cells are presented in
of triaxial apparatus and transducer, as shown in Table 3.
Figure 1. In case of in-soil calibration, the relation
Air pluviation was adopted for preparation of between applied pressure and measured strain
sand layer of 16 kN/m3 unit weight. Incremental is non-uniform for sand layer thickness of 0.25
displacements were applied manually to the triaxial DEPC and 0.5 DEPC. As the pedestal body was more
base, thereby monitoring applied stress on wooden rigid than diaphragm of EPC, arching of sand
block through load cell and induced strains in the might have caused non-uniform transfer of pres-
transducer. Strains were measured upto 50 kPa sure, in line with observations by Ingram (1968).
pressure. Results obtained from both fluid calibra- With increase in the sand layer thickness effect of
tion and in-soil calibration of both pressure sensors arching would have been reduced, thus allowing
259
260
261
ABSTRACT: This paper presents an example of S-PSA, Seismic Probabilistic Safety Assessment,
applied to an existing estuary dam. Severe damage probabilities of respective parts constituting the dam,
which are estimated by reliability analysis with using geotechnical analysis (FEM) and statistical analysis
of uncertainties, are described. Finally, the seismic retrofit planning based on risk is discussed.
Seismic design with considering occurrence prob- The target facility is an existing estuary dam to
ability of larger earthquakes than the design one supply city, industrial and agricultural water due
has been an important issue in Japan, because sev- to prevent running up of seawater to the river. The
eral large earthquakes caused tremendous damages facility length is 192.3 m and it consists of six gates
to human and social economy, such as Kobe quake (25 m in width, 7.3 m in height) and seven hanging
in 1995 and Tohoku quake in 2011. up structures of the gates. Figure 1 shows the exter-
Therefore, the authors have conducted a study nal appearance of the dam. The target structure in
on application of S-PSA, Seismic Probabilistic this study is the hanging up structures of the gates.
Safety Assessment (IAEA 1993, 2009), to com- This structure is built on the liquefiable ground
mon infrastructures except nuclear plants. In this during earthquake with 35 steel pipe piles (600 mm
paper, an example of S-PSA applied to an exist- in diameter, 9 mm in thickness, 25 m in length).
ing estuary dam is presented. Dynamic effective
stress FEM analysis is adopted in order to evaluate
the seismic response of the dam built on liquefi- 3 PROCEDURE FOR S-PSA
able ground during earthquake by using response
surface method (Honjo 2011) in this S-PSA. In Figure 2 shows the procedure for the S-PSA in this
the reliability analysis, following uncertainties are study. At first, the uncertainty of seismic intensity
considered. is set as a probabilistic density function from the
Response uncertainties: Earthquake intensity published probabilistic seismic hazard curve at the
based on probabilistic hazard curve at the site, site, which corresponds to residual service term of
reproducibility of the FEM analysis adopted in the target facility. Then, Response Surfaces (RS)
this study with respect to experimental responses
of the structure built on liquefied ground, and
reproducibility of the assumed response surface
with respect to the seismic responses estimated by
FEM analysis.
Resistance uncertainties: Strength estimation of
the structural members, bearing capacity of pile
foundation, spatial variation and transformation
error of geotechnical parameters.
Details of the S-PSA with using FEM analysis
as well as severe damage probabilities of the respec-
tive parts constituting the dam are introduced and
the seismic retrofit planning based on risk is dis-
cussed in this paper. Figure 1. External appearance of the target dam.
263
4 UNCERTAINTY OF SEISMIC
INTENSITY Figure 4. Uncertainty of seismic intensity at the site.
264
265
Verification section
Figure 8. RS of maximum shear force of gate pier.
River axis Dam axis
266
267
268
269
ABSTRACT: Analytical probabilistic analysis and Monte Carlo simulation based on elasto-plastic
Finite Element Method (FEM) on dry soil mix columns are presented. It is shown that analytical method
is over conservative because it ignores the supports from adjacent columns. Probabilistic FEM analysis
can provide more accurate predictions, and thus lead to more economic designs. Probabilistic FEM analy-
ses show that the effects of adjacent columns can be destructive when applied load is close to the strength.
The reliability of the system of columns is analyzed by setting residual strength to zero. Results show that
close spacing has more safety margin than loose spacing.
271
272
273
W
= 1.5 μCu (4) Table 2. Probability of column collapse.
A
Column Probabilistic
If the Tresca failure criterion is used, the yield strength (kPa) FEM 1D estimate
stress is
100 1/2000 0.35
W 150 0.001 0.35
σ 0 75 μCu (5)
200 0.0015 0.35
2A
250 0.0065 0.35
The factor of safety is 1.333. The probability of 300 0.0165 0.35
column collapse is 500 0.0665 0.35
⎛
P 2Cu
⎝
W⎞
A⎠
P Cu( 0.75 μCu )
Table 3. Probability of column collapse.
= Φ⎜
(
⎛ ln 0.75 μC − μl
u
) Cu
⎞
⎟ Column Spacing FEM
⎜⎝ σ llnCu ⎟⎠ strength (kPa) (m) analysis 1D estimate
274
ACKNOWLEDGEMENT
REFERENCES
275
ABSTRACT: On the basis of the interval analysis, the non-probabilistic reliability analysis method
makes the reliability of slope stability and risk analysis more complete. Through the comparison
between the safety factor, probability and reliability indexes and the non-probabilistic reliability
index of dangerous sliding surface, the non-probabilistic reliability index’s advantage is proved
and a practical evaluation method is supplied for the reliability and risk assessment of geotechnical
engineering.
277
But σ z = σ R 2 + σ S 2 (5) l
in G (δ1,δ 2 , ..., δ n )
δ i ∈Δ l
In the process of calculating the reliability Zu x G (δ1,δ 2 , ..., δ n )
index, we compare formula (4) with formula (5), δ i ∈Δ l
Z ( f Li + σ tan
)Li = cL t ϕ L i τ Li
Figure 2. Normal distribution scheme estimates (small Figure 3. The ABAQUS model of the section gravel
sample). pile roadbed.
278
n
∑ τ f Li
i =1
K= n
= 181.02/89.32 = 2.03
∑ τ Li
i =1
n
∑ μz
i =1
β= = 22.93/6.66 = 3.44
n
Figure 4. Shear stress. ∑σ z 2
i =1
Division of grid, make each quadrilateral ele- For the first embankment layer τ f and τ along
ment approximation for 2 m length, pile and the sliding surface of the integral unit length, the
surrounding transverse length of 0.5 m units, the safety coefficient and probabilistic reliability index
soil elements between the pile lateral side length is can be calculated:
about 2.7 m. The maximum stress and minimum
stress of each unit are chosen from Figures 5 and 6. n
After calculation, the reliable indicators and the ∑ τ f Li
i =1
coefficients of safety of every unit in most likely K= n
= 752.16/666.84 = 1.13
damage surface of slope are listed in Table 3.
The Table 3 shows that the safety factor and prob- ∑ τ Li
i =1
ability reliability of embankment soil are large and
the shear of muddy silty clay layer unit is large, and
n
the safety factor through calculation of every unit
is close to 1 and the probabilistic reliability index is ∑ μz
i =1
close to zero, which tends to be damaged, and the β= = 85.56/42.88 = 2
n
probability reliable index of pile is not very large.
∑ σ z2
i =1
3.2.2 The safety coefficient and reliable index
of overall reliability in every soil layer
For embankment layers τ f and τ along the sliding 3.2.3 The safety coefficient and reliability index
surface of the integral unit length, the coefficient of embankment slope’s overall reliability
of safety and probabilistic reliability index can be We can get the coefficient of safety and probabil-
calculated: ity reliability index through integrating along the
279
μc μϕ μz σz Probability Safety
Soil layer Unit (KPa) (degree) (KN/m) (KN/m) index β coefficient K
n
∑ μz
i =1
β= = 177.26/44.91 = 3.95,
n
∑σ z 2
i =1
280
sliding surface of the integral unit length, the non- ∑ (Zu − Zl )dLLi
i =1
probabilistic reliability index can be calculated:
281
REFERENCES
Figure 7. The relationship curve of k, β, η.
Junqi Zhang. Non-probabilistic reliability analysis of
slope stability. Hebei University of Technology, 2011.
we can not get the result, because it is influenced Minghua Zhao, Chong Jiang & Wengui Cao. 2008.
by parameter’s distribution. If the coefficient Non-probabilistic reliability analysis based on the the
of variation is bigger, the probability reliability interval theory retaining wall stability.The Journal of
index will decrease. So we can say that the non- Geotechnical Engineering 30(4): 467–472.
probability reliability theory is good for the slope Ruijiao Xin & Qiang Li. 2011. Structure interval fuzzy
engineering. stochastic finite element reliability analysis. Ha’erbin
University of Technology, 2011.
Shengfei Yu, Zhengzhou Chen & Mingrui Zhang. Slope
4 CONCLUSION stability analysis based on the range of uncertainty
analysis method. Journal of Engineering Geology,
1004-9665/2012/20(2)-0228.
Through the analysis of engineering examples, the Wengui Cao & Yongjie Zhang. 2007. Non-probabilistic
relationship between the three parameters is shown fuzzy reliability analysis of slope stability based on
in Figure 7. interval method. The Journal of Civil Engineering,
Figure 7 shows that the non-probabilistic reli- 2007.11.
ability indexes are similar to the safety coefficient Yong Liu & Huansheng Mu. 2011. Non-probabilistic
values. There are no clear provisions on the non- reliability analysis of highway roadbed slope stability
probabilistic reliability index, so the selecting on in soft ground conditions. Road Engineering 1002-
the value target of probabilistic reliability index 4786 17-0096-02.
can refer to the target value of safety coefficient. In Yingying Zhang. 2011. Non-probabilistic reliability
analysis based on particle swarm optimization and
addition, the quantity of soil parameters is small, its application in Slope Engineering. Zhongnan
so parameters distribution can only be assumed University, 2011.
to be normal distribution, and the coefficient of Zhengzhong Wang, Shengyuan Mou & Jun Liu. 2009.
variation is not accurate, when calculating the Calculated safety factor of slope aspect centrifugation.
probability of slope stability reliability. Combined Rock and Soil Mechanics, 2009.9.
282
ABSTRACT: Metamodelling is actually a very interesting research area for the replacement of
simulation models when the trade-off between efficiency and accuracy is appropriate. Several approaches
have been proposed, and localised metamodelling techniques applied in reliability-based analysis may
involve the search for the most probable failure point. Response surface methodology has been suggested
for the approximation of functions, namely implicit limit state functions. The response surface is typically
constructed from a polynomial function fitted at a number of points, but no clear guidance has been
given for the effect of varying their location. Results for a geotechnical design example with correlated
nonnormal random variables are presented in the light of the quality of the design point approximation.
In particular, interpolation methodologies in conjunction with appropriate experimental design are
suggested for extracting nonlinear response surfaces whose roughness varies substantially over the input
domain, so that sequential algorithms may adapt to nonlinear features.
283
284
be applied beyond structural engineering; where the sum of B1 and B2 is the foundation width B;
further check the convergence process for H1 is the wall height; H2 is the foundation height; γc
the coordinates of the design point; is the unit concrete weight; ϕw is the friction angle of
[Step 9] If required, repeat steps 7 and 8 until the soil on the active and passive sides of the wall;
convergence to the selected tolerance γw is the unit soil weight on the active and passive
criterion, then express the final response sides of the wall; cf is the cohesion of the founda-
surface, and estimate the reliability index tion soil; ϕf is the friction angle of the foundation
β and the corresponding coordinates of soil; γf is the unit weight of the foundation soil; and
the design point; Monte Carlo simulation q is the variable load at ground surface. Other con-
based on counting is now efficient since sidered parameters, namely for the earth pressure
the evaluation of the response surface coefficients, are the soil-wall interface friction angle
function requires a reduced computa- on the active side of the wall δ w ϕ w; the soil-
tional effort and input random variables wall interface friction angle on the passive side of
are correlated in the standard normal the wall δwp = 0; and the soil-foundation interface
space; it is further noted that transforma- friction angle δf = ϕf; considered either a 50% reduc-
tions to the space of uncorrelated ran- tion for passive earth pressures.
dom variables may be considered. Table 1 summarises the description of basic
Thus, the efficiency and accuracy provided by input variables, with different distribution types.
different schemes of experimental design are ana-
lysed in the next sections with a design example for Table 1. Summary description of basic input variables.
a concrete gravity retaining structure, considered
correlated nonnormal random variables. Basic input
variables Distribution mv cv Statistics
3 DESIGN EXAMPLE B1 (m) Deterministic – – 1.75a
1.85b
The design example is referred to the concrete grav- 1.95c
ity retaining structure on a relatively homogeneous B2 (m) Deterministic – – 1.00
c − ϕ soil shown in Figure 1, wherein groundwater H1 (m) Deterministic – – 7.00
H2 (m) Deterministic – – 1.00
γc (kN/m3) Deterministic – – 24.00
ϕw (º) Lognormal 33.00 0.10 μ1 − 3.4915
σ1 − 0.0998
γw (kN/m3) Normal 18.80 0.05 –
cf (kN/m2) Lognormal 14.00 0.40 μ1 − 2.5648
σ1 − 0.3853
0.05a μ1 − 3.4645
σ1 − 0.0500
ϕf (º) Lognormal 32.00 0.10b μ1 − 3.4608
σ1 − 0.0998
0.15c μ1 − 3.4546
σ1 − 0.1492
γf (kN/m3) Normal 17.80 0.05 –
q (kN/m2) Weibull – – sh−3.25
sc−10.00
285
Correlation matrix for the bearing capacity model Coefficient Reliability Reliability
of variation index index
ρϕ ϕ ρϕwγw ρϕwϕf ρϕwq ρϕwcf ρϕwγf Model ϕf βFORM βMCS
w w
ργwϕw ργwγw ργwϕf ργwq ργwcf ργwγf
ρϕfϕw ρϕfγw ρϕfϕf ρϕfq ρϕfcf ρϕfγf Bearing 0.05 2.4842 2.5288
capacity 0.10 2.2248 2.2787
ρqϕw ρqγw ρqϕf ρqq ρqcf ρqγf
ρcfϕw ρcfγw ρcfϕf ρcf q ρcfcf ρcfγf 0.15 2.0410 2.0952
ργfϕw ργfγw ργf ϕf ργf q ργfcf ργfγf MCS results from 106 simulation steps.
1.0 0.5 0.9 0.0 0.0 0.5
0.5 1.0 0.5 0.0 0.0 0.9
0.9 0.5 1.0 0.0 0.0 0.5 variation for the friction angle of the foundation
0.0 0.0 0.0 1.0 0.0 0.0 soil. It is noted that the differences between the reli-
0.0 0.0 0.0 0.0 1.0 0.0 ability index obtained from both methodologies are
0.5 0.9 0.5 0.0 0.0 1.0 somewhat notorious due to the important nonlinear
behaviour presented by the bearing capacity model.
ρ—coefficient of correlation. The bearing capacity metamodels are built
according to different schemes of experimental
design of axial type, for a selected distance from
Correlation is the property in that, on average,
the center point of 0.01, 1, 2, 3 and 4. Grouped
two variables are linearly associated with one
these models, a comparative study about the con-
another, and the value of one provides information
vergence process is then summarised in Tables 4
on the probable value of the other. The strength of
and 5, for two convergence criteria and three dif-
the association is measured by a correlation coef-
ferent reference models corresponding to coeffi-
ficient in the range [−1,1], for the case that both
cient of variation ϕf 0.05, 0.10 and 0.15, considered
variables vary inversely or together exactly, respec-
the first order reliability method.
tively. The considered coefficients of correlation
According to Haldar & Mahadevan (2000), the
between the random variables are either presented
number of design sampling points should be kept
in Table 2, noted that the correlation matrices are
to a minimum to increase efficiency but must be
idealised with the main purpose of creating high
at least equal to the number of coefficients needed
correlation between some of the random variables.
to define the performance function, therefore effi-
It is further noted that the correlation between the
ciency is related to the number of design sampling
cohesion and the friction angle of the foundation
points and iterations, in other words the number
soil depends particularly on the material, and the
of design of experiments. Haldar & Mahadevan
question whether the two parameters are corre-
(2000) refer also that the efficient location of
lated is still not clearly decided in the literature, as
design sampling points around the center point
stated by Fenton & Griffiths (2003). According to
is essential for the accurate construction of the
Marques (2012), references for negative correlation
performance function, and increased accuracy is
are common but regarding a sensitivity analysis for
described in relation with an increased number of
the design example, the influence of negative cor-
design sampling points, with the discussion of dif-
relation on reliability is favourable and quite con-
ferent layouts. These issues are focused hereafter.
siderable when compared to the effects of positive
Thus, two convergence criteria are considered in
correlation. Thus, the hypothesis of uncorrelated
order to compare the accuracy of the models, given
random variables is probably cautious and further-
by the relative errors obtained by the first order
more, the characteristic values of soil properties are
reliability method, and also the efficiency from the
often calculated separately, neglecting the effects of
required number of iterations. It is noted yet that the
correlation, so both procedures are concordant.
percentage errors are currently referenced in the lit-
erature to Monte Carlo simulation based on count-
4 RESULTS AND DISCUSSION ing, although the reference is in fact the first order
reliability method, used in the iterative algorithm.
Considered the implementation of the foremost From the analysis of results, it is concluded that
detailed methodology, Table 3 presents the sum- the efficiency and accuracy of each model depend not
mary description of results for the bearing capacity only on the number of design sampling points and
model, included the reliability index obtained by the layout, but also on the distance from the center point
first order reliability method and the Monte Carlo that is considered to set the axial points, even whether
simulation, considered the different coefficients of the number of design sampling points and layout are
286
similar. When this distance is reduced in order to be reference models. The methodology is based on the
close to zero, the search for the design point in non- local approximation of the coordinates of the most
linear analysis is improved. It is also concluded that probable failure point, the design point, unlike the
as this distance becomes larger, the different models stochastic response surface methodology which
may present inaccurate results and important conver- may provide a general approximation statistically
gence problems, in particular for the tighter tolerance representative of the performance function.
criterion and the higher coefficient of variation for
the friction angle of the foundation soil.
The response surface methodology itera- 5 CONCLUSION
tive algorithm generates nonstochastic local
metamodels, in that they present statistical param- The design of civil engineering structures relies
eters more or less nonrepresentative of their heavily on models, abstract and simplified
287
metamodels built by the response surface meth- where the sum of B1 and B2 is the foundation
odology iterative algorithm show that improved width B; H1 is the wall height; H2 is the founda-
techniques for the search of the design point are tion height; γc is the unit concrete weight; ϕw is the
required, in particular if a tighter tolerance crite- friction angle of the soil on the active and passive
rion is considered. It is noted that the efficiency sides of the wall, expressed in degrees; γw is the unit
and accuracy of the iterative algorithm are very soil weight on the active and passive sides of the
sensitive to the distance from the center point that wall; cf is the cohesion of the foundation soil; ϕf is
is considered to set the axial points, so when this the friction angle of the foundation soil, expressed
distance is reduced in order to be close to zero, the in degrees; γf is the unit weight of the foundation
search for the design point in nonlinear analysis is soil; q is the variable load at ground surface; Nc,
improved. Nq, and Nγ are the bearing resistance factors; sc, sq,
Going back in history, several cases of instabil- and sγ are the foundation shape factors; ic, iq, and
ity have been reported in the literature regarding iγ are the load inclination factors; Be is the effec-
the use of reduced distances in schemes of experi- tive foundation width; Evb is the resultant vertical
mental design, but results for the selected design action for bearing capacity; and Ehb is the result-
example show a possibility for progression when ant horizontal action for bearing capacity.
considered the detailed response surface method-
ology iterative algorithm. It is noted lastly that the
axial type schemes are not always the most appro- REFERENCES
priate option.
Allaix, D.L., Carbone, V.I. & Mancini G. 2011. Advanced
response surface method for structural reliability. In:
APPENDIX 1 EQUATION (3) Proc. of the Thirteenth International Conference on
Civil, Structural and Environmental Engineering Com-
puting, Chania, Crete, Greece: Civil-Comp Proceedings
M (
f B B H H c w w cf f f q ) (3) 96, Paper 88, 1–14. Civil-Comp Press.
288
289
H. Nagao
Naigai Engineering Corp., Kyoto, Japan
K. Fujisawa
Graduate School of Agriculture, Kyoto University, Kyoto, Japan
ABSTRACT: The risk of earth fills during heavy rains is evaluated in this research. The rainfall intensity
and the variability of the numerous soil parameters of the earth-fills related to soil erosion are dealt with as
probabilistic parameters. When the peak overflow head on a spillway bed becomes greater than the design
overflow head, overflow occurs. The peak overflow head, is determined within 72 hours of the rainfall,
and the various rainfall patterns are tested by the Monte Carlo method. Furthermore, the statistical values
of the parameters for the soil erosion are assumed, and the variability of the discharge hydrographs are
derived from the values. The discharge hydrographs are applied for the flood simulation as input waves. The
submerged area and the head of the flood discharge in the downstream area can be predicted. Considering
disaster loss, the risk to the downstream area of the earth-fill dams can be evaluated.
291
Figure 1. Cumulative probability distribution of rain- Figure 2. Cumulative distribution of maximum annual
fall intensity after 6 hours. continuous precipitation.
292
Qin f prA / ..6 (9) Once the event of the overflows happens, the flood
simulation must be done as a next step to estimate
the damage in the downstream area. The discharge
where Qin = inflow to the reservoir (m3/s), fp is
hydrographs are determined to be applied for
the peak runoff coefficient, r is the quasi-rainfall
the flood simulation as input waves. An example
intensity (mm/h) generated in chapter 2 and A is
of section of the embankments analyzed in this
the area of the basin (km2). Uniform random num-
study is shown in Figure 4. When the overflow
bers are used for fp in the 0.7 to 0.8 range (JSIDRE
occurs, the embankment is supposed to be failed
2002). The discharge equation for a rectangular
as in Figure 5, which shows a vertical section of an
weir as used in this study is;
embankment. This section of the flood way is to be
assumed to have a trapezoid shape. The overflow
Qout CBS h2 3
CB (10) hydrograph from broken section is determined by
the runoff discharge Qdis. is formulated as follows;
where Qout = discharge (m3/s), C = discharge
coefficient, Bs = width of spillway, and h = static
Qdiis = ggA3 ( ∂A / ∂h ) (14 )
or piezometric head on a weir referred to the weir
crest. The storage of water in the water reservoir U Qdiis / A (15)
Vr is estimated as follows;
dH / dt = − E /(1 n p ) (16 )
Vr Aw h (11)
E = α (τ − τ c ) (17 )
where Aw is area of water reservoir (km2) and h is τ ρC f U 2
( )
overflow head (m). The decreasing rate of the stor-
age V with the runoff is: Cf 2
/l {R exp((κ Ar − 1)}/
)}// s] ( )
Vr dt = Qin Qout
dV (12)
where A is the cross-sectional area of broken
The overflow head h is determined from section (m2), U is the flow velocity (m/s), H is the
Equation (12), and the maximum h within the overflow head on the bottom of the failed section
72 hours, is defined as the peak overflow head on (m), E is the erosion rate of embankment material
the spillway. When hp becomes greater than the (m/s), np is the porosity of the embankment, α is
design overflow head hd, the overflow occurs. Then, the erosion rate coefficient (m/s/Pa), τ is the shear
the probability of overflow is defined by Equation stress (Pa), τc is the critical shear stress (Pa), ρ is
(13) as the times of hp < hd in the iterations of the the water density (kg/m3), κ is Karman constant,
Monte Carlo simulation (Rubinstein 1981). R is the hydraulic radius (m), Ar is the constants
(=8.5) and κs is roughness height (m).To consider
uncertain material values, the variabilities of the
Pf ⎡ hd < hp ⎤ (13)
⎣ ⎦ several soil parameters of the earth-fills are dealt
293
(21)
294
( )
dU i 1 1 1
dt
=−
Ωi
∑ E ⋅ nij ΔΓ ij + Si (25) gh* =
2
ghhL + ghR +
4
( L − R )⋅n (30)
j =1
⎧(E ) ⋅ n ( ≥ )
⎪ 5 RELIABILITY-BASED DESIGN
⎪ SR (E ) n − SL (E ) n
⎪
E n = ⎨ + SRSL ⎡⎣( U ) ( U )⎤⎦ The probability of submerge in a cell i per year is
⎪ SR − SL
(SL < < ) defined as:
⎪
⎪( ) ⋅ n ( ≤ )
NE
⎩ Pi ∑ Pf ,ii jjp ⋅ Psi jjp (31)
(26) jp =1
where UL and UR = reconstructions of U on the left in which the subscript jp corresponds to the event,
and right sides, respectively, and SL and SR = wave which is the combination of the failure dams, NE
speed estimates. is number of events, Pf is probability of overflow,
and Ps is the probability of submerge under the
⎧ ⎛ q ⋅ n − ghL ,⎞ condition that the overflow happens.
⎪min ⎜ *L ⎟ if both sides are wet The expected total cost within the lifetime
⎪⎪ ⎝u ghh* ⎠ period = t is given by the following equations.
SL = ⎨
⎪ q L ⋅ n − ghL if the right
i side is dry M
⎪ CT C + ∑ C f ,i E [ ni ] (32)
⎪⎩ q R n − 2 ghR if the left f side is dry
i
(27)
E [ni ]
⎧ ⎛ q ⋅ n − ghR ,⎞ ⎧ tI
⎪max ⎜ *R
⎪⎪ ⎝u ghh*
⎟ if both sides are wet
⎠
⎪∑ ⎡P (1 − PO ,i )
= ⎨ k =1 ⎣ O ,i
k −1
+ (tl { } (
k ) PI i ⎤⎦ )
SR = ⎨ ⎪
⎪ q L ⋅ n − ghL if the right
i side is dry ⎩tl ⋅ PI ,i (Improved)
⎪
⎪⎩ q R n − 2 ghR if the left f side is dry (33)
(28) in which M is number of cells, CT is the expected
total cost n is the frequency of overflows within a
lifetime span of t (years), PO and PI are the prob-
abilities of overflow a year corresponding to the
original and the improved states of the embank-
ment, respectively, C0 is the cost of the improve-
ment, and Cf is the damage loss due to flooding. PO
and PI are derived from Equation (31).
In this research, the improvement of the spill-
way is considered, and the improvement brings
about a drastic increase in the discharge ability of
the spillway. In Equation (32), improvement cost
C0 is zero for the original state of the embank-
ment, and it is assumed that when embankments
are failed due to overflows, that they will be
restored to the same level as the improved state of
Figure 7. Boundary condition and input cells. the embankment.
295
Pf (A ∩ B ) ( ) (
0 Pf A ∩ B ≠ 0, Pf A ∩ B = 0 )
Pf (C ∩ D) 0 Pf (C D) = 0, Pf (C ∩ D) = 0 (34)
296
297
ABSTRACT: The authors proposed a reliability analysis method for long continuous structures
considering investigation sites to counter liquefaction failure. This method is characterized by continu-
ously calculating the failure probability caused by liquefaction at an arbitrary point and quantifying the
uncertainty concerning the investigation site. The purpose of this study is to verify its effectiveness based
on application to an actual river dike that suffered liquefaction failure due to the 2011 Earthquake off the
Pacific coast of Tohoku.
299
Figure 2. Temporal history of strong movement. Figure 3. Collapse of river dike caused by liquefaction.
parameters for the dike body and alluvium. The 2 COLLAPSE RISK ASSESSMENT USING
river body soil is mainly sandy silt (Bc), but sandy PL-VALUE
earth filling has possibly accumulated under the Bc
in the former riverbed area. The liquefaction risk at an arbitrary point is
The alluvium layers under the dike body are evaluated using the liquefaction index i.e. PL-value,
roughly classified into two types, silt or sandy silt as proposed by SFHB2002 (JRA 2002). In fact,
(Ac) (which is a non-liquefiable layer) and medium we calculated the correlation between the PL-value
or fine sand (As) (which is liquefiable), and the total and collapse risk as follows:
thickness of the alluvium layer is between 10 and
1. Calculation of PL-value at all investigated sites.
25 m. Furthermore, diluvium, which is regarded as
2. Development of empirical fragility curve in the
engineering bedrock, has accumulated under these
form of a normal distribution function using a
alluvial strata.
two-value regression analysis based on informa-
The soil investigations mainly consist of a
tion about whether the site collapsed or not.
Standard Penetrating Test (SPT) and determination
of the particle size. There were 83 SPTs in total,
but only ten particle size determinations. 2.1 Calculation of PL-value
2.1.1 Calculation procedures
1.3 Profiles in strong movement and collapse PL-value at the investigation sites was basically
status of target river dike based on SFHB2002, but there were improvements
to the original method in order to take the long
Figure 2 shows the observed time history of accel- duration and local site of shaking amplification
eration at the site located on bedrock near the into consideration. The calculation procedures are
target river dike. The maximum acceleration of described based on notations in the improvement
this earthquake wave was moderate by Japanese points as follows.
standards, but it was characterized by a very long Liquefaction evaluation along the depth was
duration. calculated as follows.
The total duration was around 200 seconds,
and during this period, the main shaking occurred FL R /L
L (1)
between 90 and 145 seconds (i.e. for 55 s).
In addition, four points along the target river R cW RL (2)
dike interval collapsed after the shaking (Fig. 1).
A summary of the collapsed status is presented in where FL is the resistance ratio of liquefaction, R
Figure 3. All of the collapsed area was along the is the dynamic shear resistance ratio, L is the shear
former riverbed or natural levee, and was consid- stress load during strong movement, cw is the adjusted
ered to have been caused by liquefaction of earth coefficient of strong movement duration, and RL is
filling Bs, which resulted in leveling of ground on the dynamic triaxial resistance. RL is defined in Equa-
the flood plain. tion (3), and is correlated with SPT N-value.
300
Figure 4. Results of ground response analysis in shaking and PL-value at all investigational sites.
301
2
20
⎛ 10 − 0.5x ⎞
Var [ PL ] = ∑ ⎜⎝ L i ⎟⎠ Var[R L ] (9)
i=1 i
N
L ∏ [ FFRR(l PLi )]α i [ − FR(lln PLi )]1−α i (10)
i =1
∂ ln L ∂ ln L
= =0 (12)
∂c ∂ξ
302
303
Figure 9. Estimated PL-value, failure probability and contribution of uncertainties (Zone 2, STA 35-45).
but it was considered that additional investigations Otake, Y. & Honjo, Y. 2012. Reliability based design
should be conducted. on long irrigation channel considering the soil
investigation locations. Proceedings of Geo Congress
2012, 2836–2845.
Otake, Y., Honjo, Y. & Koike, K. 2012. Evaluation of
4 CONCLUSION soil liquefaction risk on linearly structure considering
the soil investigation locations, Japanese Geotechnical
Conducting countermeasures for a long continuous Journal 7(1): 283–293. (in Japanese).
structure such as this river dike is time-consuming Schnabel, P.B., Lysmer, J. & Seed, H.B. 1972.
and costly. The information provided here is very SHAKE-A computer program for earthquake response
useful for planning such a construction and deter- analysis of horizontally layered sites, Report EERC
mining its priorities. 72-12, Berkeley California, Earthquake Engineering
It should be pointed out that the owner could Research Center, University of California.
obtain more information on soil properties by Yoshida, N., Ohya, Y., Sawada, S. & Nakamura, S. 2009.
Simplified procedure for evaluating liquefaction
conducting additional soil investigations. The potential under ocean trench type long period earth-
results from this method provide very valuable quake, Journal of JAEE 9(3): 28–47. (in Japanese).
information in this regard.
REFERENCES
304
ABSTRACT: This study investigated factors that control the reliability of Augered Cast-In-Place (ACIP)
piles in predominately cohesionless soils under axial compression at the Serviceability Limit State (SLS).
A simple probabilistic hyperbolic model was used to account the uncertainty in the load-displacement
relationship using correlated bivariate curve-fitting parameters. Contrary to previous studies, the curve-
fitting parameters were found to be dependent on pile slenderness ratio (D/B) and the effect of D/B and
other pertinent variables (e.g., uncertainty in capacity, displacement) on SLS reliability was investigated
using a First-Order Reliability Method (FORM). The D/B ratio had a considerable effect on founda-
tion reliability, illustrating the importance of the dependence between the load-displacement behavior
(i.e. curve-fitting parameters) and pile geometry and stiffness. In general, the uncertainty in the capacity
model had a larger effect on reliability than that of the allowable displacement; the reliability index was
found to approach an upper bound limit regardless of the level of uncertainty in allowable displacement
and the pile capacity model.
305
Q y
= (1)
QSTC k1 k2 y
3 DATABASE
Navg QSTC
Variable D (m) B (mm) D/B (bl/0.3 m) (kN)
Figure 1. The dependence between slenderness ratio,
Minimum 7.5 300 20.0 4 367 D/B, and model parameters, (a) k1 and (b) k2 and the
Maximum 29.0 800 68.5 54 5300 corresponding Kendall tau correlation coefficients and
p-values.
306
Figure 2. Correlation between model parameters (a) k1 Figure 3. Empirical, lognormal, and normal marginal
and k2 and (b) k1,t and k1,t and the corresponding Kendall cumulative distributions for the transformed hyperbolic
tau correlation coefficients and p-values. model parameters: (a) k1,t, and (b) k2,t.
307
6 TRANSLATIONAL MODEL
FOR BIVARIATE PROBABILITY
DISTRIBUTIONS
and where λ1,t, ζ1,t and λ2,t, ζ2,t are the approximate 7 RBD FOR THE SERVICEABILITY
lognormal mean and standard deviation of k1,t and LIMIT STATE USING A FIRST-ORDER
k2,t, respectively, and ρ is the standard product- RELIABILITY METHOD
normal correlation coefficient for two normally
distributed variables. The second moment statistics The SLS is reached when foundation displacement,
in Equation 4 were calculated as: y, equals or exceeds allowable settlement, ya. This
study followed the approach outlined in Phoon &
Kulhawy (2008), where the SLS can be evaluated
ζ i ,t ( σ i2,t / ki2,t ) (5a) using a performance function, P:
The correlated, lognormal marginal distribu- Failure is defined as P ≤ 0, and the probability
tions of k1,t and k2,t were thus simulated using: of exceeding the SLS, pf, is:
k2,t e
(ζ 2,t X 2 + λ2,t ) (6b)
By combining Equations 1, 7, and 8, and defin-
ing a deterministic mean global factor of safety,
FS, the probability of failure is:
In order to adequately reproduce the uncer-
tainty in the observed load-displacement curves, ⎛ ya 1 Q′ ⎞
k1,t and k2,t must be back-transformed into k1 and p f = Pr ⎜ < ⎟ (9)
k2. This study calculated k1 and k2 using determin- ⎝ 1
k k y
2 a FS
F Q′ p ⎠
istic values of D/B because the uncertainty asso-
ciated with pile geometry could not be evaluated where Q′ and Q′p correspond to the applied load
from the database. and predicted pile capacity, respectively. In order
308
β = − Φ −1 ( ) (10)
309
310
Pham Quang Tu
Delft University of Technology, The Netherlands
Hanoi Water Resources University, Vietnam
ABSTRACT: The province of Thai Binh, an area with a large population and high socio-economic
status, is located at the downstream part of the Red River Delta. The area is protected by dike systems
whose safety has been upgraded significantly. While overtopping was an important failure mode in the
past, after raising crest levels, now piping and uplift phenomena are considered the dominant threats, hav-
ing led to several dike breaches in the Red river dike system. In order to address this issue and the large
uncertainties in commonly used uplift models, a series of field tests had been designed and carried out
on the Thaibinh formation (a regional top soil layer) for determining the critical uplift pressure. A special
structure is installed to increase water head artificially, by excavation through the top layers, and ground
water flow is concentrated by plastic sheet piles. Phenomena are monitored and measured in an excava-
tion, for thinning the thickness of impervious layers, during the test. The results cover critical head, uplift
displacement, time series and visual observations. In this paper, the authors will perform a multivariate
analysis for the field test results in order to develop statistical model for uplift mechanism. Furthermore, a
proposed model of uplift has also been generated including a calibration coefficient which represents for
uncertainties of both inherent and epistemic sources. The proposed model is also compared to test results
and statistical model to give a better understanding of this phenomenon. The findings of the research
might contribute to improving dike assessments in the Red River Delta.
311
312
313
ρsat − ρw
ico = (4)
ρw
314
⎛ ⎞
ρ d + θ ⎜ C d + d K o ⋅ tanφ ⎟
2
hcpro (7)
⎝ 3 ⎠
In this case, groundwater level is equal to the Figure 10. Critical gradient of uplift in different models.
level of hinterland and parameters of foreland’s
blanket layers are the same as that of hinterland; 3.2 Discussion
hcpro c is critical uplift head of the proposed model; From the observation during the tests, it can be
θ is the calibration coefficient. seen that uplift commences first at the weakest spot
Critical gradient of uplift can be calculated from in the impervious layers, then develops widely on
Equation 7, then we have: area if water pressure till increases continuously.
The weak spots mentioned here are cracks or het-
⎛ ⎞
= ρ ′ + θ ⎜ C + d K o ⋅ tanφ ⎟
hc pro 2 erogeneous structures in soils and the drought
icpro = (8) cracks are also supposed being a main disturbed
d ⎝ 3 ⎠
factor of the top layers in the depth from 0.8 to
Calibration coefficient can also be figured out by: 1.2 m, see (Bui 2009).
After crack, heave is also visualized in all test
3( cpro − ρ ′ ) locations and the moment of heave start is defined
θ= (9) as end of crack developing process. Time depend-
3C + d ⋅ K otanφ ent in test series may also be influenced to its
results. Unfortunately, all of these issues are not
By fitting icpro to uplift critical gradient from included in this study.
test’s data series, we could elaborated the values The statistical model fitted quite well with data
of θ. series of test, see Figure 8. However, it has not
shown the physical relationship between parameters
3 RESULTS AND DISCUSSION and its physical behaviours. It is clear from statisti-
cal analysis that, correlated coefficient of this model
3.1 Results has a negative trend (ρ = −5.3%) comparing to posi-
tive trend of two other models (ρ = 8.9–16.9%). It
Both uplift and heave are observed during the seem unreasonably that different heads increase
increasing of water level indicating for flood water lead to decline of critical gradient, following the
level. Crack was initiated first in weak spot of the statistical model. So that, in practice it might be
blanket, then appears in many spots with sand or incorrect if we apply the statistical model without
silt flowing out in a certain area. understanding about its physical behaviours.
There are total of fifty eight tests performs on It also can be seen from Figure 10 that the pro-
Thaibinh formation for uplift and heave phenom- posed model has the smallest variation comparing
ena. The plotting of the final results of critical gra- to others models. Calibration coefficient has been
dient of uplift could be found in Figure 10. assumed to be normal distributed and parameters
The critical uplift gradients are found to be nor- are mean and standard deviation by μθ = 0.5519
mal distributed of mean and standard deviation and σθ = 0.746 respectively. In practice, it is recom-
equivalent to 0.836 and 0.105 respectively. From mended that θ should be chosen by 0.5 for engi-
the theoretical Equation 4, the so-called theoreti- neering purposes.
cal critical gradient has a value of 0.758. Whereas, In term of comparison between statistical to pro-
in the proposed model, critical gradient could be posed model, the probability if actual critical gradi-
formulated in range of normal distribution with ent is larger than that of theoretical model could be
parameters N(0.8301;0.0497). It is figured out that calculated, for instance P (i icpro ) P (i .758)
statistical model fits quite well to normal distribu- The final results show that exceeded probability is
tion with parameters of N(0.8323;0.066). P (icpro ico ) = 92. % for the proposed model and
315
REFERENCES
316
ABSTRACT: Rainfall is one of the most common factors triggering landslides, since infiltration of
water into the soil has a significant impact on pore water pressure buildup that affects slope stability. In
this study, the influence of the wetting front development on the reliability of an infinite slope is analyzed.
The failure condition of the slope is expressed in terms of the factor of safety. Rainfall infiltration is
simulated by a time-dependent model, based on the Green and Ampt assumptions. The vertical variability
of the saturated hydraulic conductivity and soil strength parameters are modeled as random fields. The
reliability of the slope is evaluated by Monte Carlo simulation. A numerical example demonstrates the
influence of the vertical variability on the analysis results.
317
318
l
2.4 Pore water pressure distribution
in which Kl is the saturated hydraulic conductivity The pore water pressure varies with time and space
of each wetted computational layer. within one unit slice. In this study, emphasis is on
We assume that the rainfall supply is larger than the pore water pressure within the wetted zone and
the infiltration capacity. Therefore the cumulative below the water table.
infiltration will be governed by the infiltration Within the wetted zone, the spatial variability of
capacity and can be obtained by the saturated hydraulic conductivity may introduce
F z ⋅ Δθ (6) large hydraulic gradients that will result in vari-
able pore water pressure. Employing the plug flow
in which Δθ is the change in moisture content, assumption given by the Green and Ampt infiltra-
given as follows: tion model, we assume that the effective flow rate
of the wetted zone f computed by Equations (4)
Δθ η − θ0 (7) and (9) equals the flow rate at each wetted layer
(Liu et. al. 2008), i.e.
where η is the porosity of the soil and θ0 is the
initial moisture content (see Fig. 2). Noting that f1 f2 = … = fl f (12)
fin ddF /ddt , we can obtain the cumulative infiltra-
tion time as follows: where fl is the flow rate at wetted layer l, computed
by
t z Δθ
∫0 dt = ∫0 fin
dz (8)
Δhl
fl K l ⋅ il Kl ⋅ (13)
Δz
The second phase that follows the rainfall event
is governed by plug flow. It is assumed that the wherein il is the hydraulic gradient of layer l; Δhl
suction at the bottom of the wetted zone equals the change of hydraulic head within the lth layer.
the one at the top [see Fig. 3(b)]. Therefore, the Combining Equations (12) and (13), we obtain:
hydraulic gradient equals unity and the plug flow f
rate fp equals the harmonic mean of the saturated Δhl Δz (14)
hydraulic conductivities of the layers correspond- kl
ing to the wetted zone:
At the bottom of lth wetted layer, the hydraulic
fp Ke head hl is obtained by summing the incremental
(9) heads Δhk of each computational layer k ≤ l
Therefore, the time needed for the wetting front
to reach the bottom of each layer is obtained as hl h − ∑ k l Δhk (15)
follows:
in which h0 is the initial hydraulic head at the top of
⎧ z Δθ
the wetted zone, i.e. h0 = 0 during Phase I [Fig. 3(a)]
⎪ ∫0 f dz , t T
and h0 = −z0 − S during Phase II [Fig. 3(b)], where
⎪ in
t=⎨ (10) z0 is the initial elevation head. Since the hydrau-
⎪T +
z Δθ
⎪⎩ ∫z L f p dz , t T lic head hl consists of the pressure head ψl and the
elevation zl, the pressure head at the bottom of zl th
layer is evaluated as:
wherein T is the duration of the rainfall event; L
is the length of wetting front exactly after rainfall ψ l = hl zl (16)
319
320
321
REFERENCES
Figure 8. Influence of the coefficient of variation of
the saturated hydraulic conductivity on the factor of Chen L. & Young M.H. 2006. Green-Ampt infiltration
safety for r = 0.5 m. model for sloping surface. Water Resources Research
42(7): W07420.
Corominas, J., Moya, J. & Hürlimann, M. 2002. Landslide
rainfall triggers in the spanish eastern pyrenees. Proc.
4th EGS Plinius Conference. Mallorca, Spain.
Der Kiureghian, A. & Ke, J.-B. 1988. The stochastic finite
element method in structural reliability. Probabilistic
Engineering Mechanics 3(2): 83–91.
Freeze, A.R. & Cherry, J. A. 1979. Groundwater. Upper
Saddle River, NJ: Prentice Hall.
Gelhar L.W. 1986. Stochastic subsurface hydrology from
theory to application, Water Resources Research 22(9):
135–145.
Green, W. H. & Ampt, G. A. 1911. Studies on soil physics,
Figure 9. Influence of the coefficient of variation of part I—the flow air and water through soils. Journal
the saturated hydraulic conductivity on the probability of Agricultural Science 4(1): 1–24.
of failure for r = 0.5 m. Griffiths, D.V., Huang, J. & Fenton, G.A. 2011.
Probabilistic infinite slope analysis. Computers and
Geotechnics 38(4): 577–584.
which will cause the wetting front to move slower Liu J. T., Zhang J. B. & Feng J. 2008. Green–Ampt Model
towards the water table and therefore failure will for Layered Soils with Nonuniform Initial Water Con-
be governed by shallow slip surfaces. Moreover, tent Under Unsteady Infiltration. Soil Science Society
a large coefficient of variation leads to steeper of America Journal 72(4): 1041–1047.
changes between small and large values of the Muntohar, A.S. & Liao H. 2010. Rainfall infiltration:
hydraulic conductivity within the wetted zone; infinite slope model for landslides triggering by rain-
storm. Natural Hazards 54(3): 967–984.
hence the pore water pressure at layers above the
Phoon, K.K. & Kulhawy, F.H. 1999. Characterization
wetted zone increases, which leads to a larger prob- of geotechnical variability. Canadian Geotechnical
ability of shallow slope failure and smaller corre- Journal 36(4): 612–624.
sponding mean factors of safety. Polemio, M. 1997. Rainfall and Senerchia landslides,
southern Italy. Proc. 2nd COBRAE, Rio de Janeiro.
Rawls, W.J., Brakensiek, D.L. & Saxton, K.E. 1982. Esti-
4 CONCLUSION mation of soil water properties. Transactions of the
ASAE 25(5): 1316–1320 & 1328.
Santoso, A.M., Phoon, K., & Quek, S. 2011. Effect of
In this paper, we developed a simplified time-
1D Infiltration Assumption on Stability of Spatially
dependent model to study the reliability of slopes Variable Slope. Geo-Risk 2011: Risk Assessment and
subjected short-term and intense rainfall events. Management. 704–711.
The model incorporates infinite slope stability Wu, T.H. & Abdel-Latif, M.A. 2000. Prediction and
analysis, one-dimensional infiltration analysis mapping of landslide hazard. Canadian Geotechnical
based on the Green and Ampt assumptions as Journal 37(4): 579–90.
well as a random field modeling of the effective Zêzere, J.L. Trigo, R.M. & Fragoso, M., Oliveira, S.C. &
friction angle and the saturated hydraulic con- Garcia, R.A.C. 2008. Rainfall-triggered landslides in
ductivity of the soil. An example demonstrated the Lisbon region over 2006 relationships with the
north Atlantic Oscillation. Natural Hazards and Earth
the influence of the correlation length and coef-
System Science 8(3): 483–499.
ficient of variation of the hydraulic conductivity
on the development of the factor of safety and
probability of failure with time. It is shown that
a small correlation length and high coefficient
322
ABSTRACT: Internal erosion and transport of fine particles are among the possible consequences of
the infiltration of rainwater into a slope. The transport of fine particles may lead to local variation of
soil porosity and hydraulic properties, which affects the stability of a soil slope under rainfall infiltration.
In this study, the governing equation for transport of fine particles based on the conservation of mass
of moving particles is coupled with the governing equation of transient seepage analysis for unsaturated
soils. An internal erosion law is used to describe the relationship between the rate of eroded fine particles
and the hydraulic gradient. The effect of internal erosion on change of porosity and saturated perme-
ability is also considered. The influence of transport of fine particles on slope stability under infiltra-
tion is investigated by a finite-element analysis. The numerical results show that the rainfall infiltration
can induce reduction of the density of fine particles in the shallow depth of soil slope. This can further
increase the permeability and water seepage in the slope and the slope stability is reduced consequently.
A parametric study is also conducted to investigate the effect of saturated coefficient of permeability on
water seepage and slope stability.
323
324
n3 (1 − n0 )2
ks = × ks 0 (12)
(1 − n ) 2
n03
325
326
327
REFERENCES
328
A.H. Albatal
Arab Academy for Science, Technology and Maritime Transport, Cairo, Egypt
H.H. Mohammad
Zagazig University, Zagazig, Egypt
ABSTRACT: The appropriate information or data is the keystone of any successful design. Site
investigation is one of the first steps of any construction projects. The purpose of a soil subsurface inves-
tigation is to provide data concerning the engineering properties of the soil for the proper design and
safe construction of a project. The site investigation phase of any geotechnical design plays a vital role to
provide the geotechnical engineer by the most appropriate data to ensure that the design data represent
the investigated soil. Inadequate soil investigation may contribute to either a significantly over designed
foundations, that is not cost-effective, or an under designed foundations, which may lead to potential
failures. Insufficient geotechnical investigation is one of first sources of projects’ delays, disputes, claims,
and projects’ cost overruns. This paper aims to focus on the impact of varying the scope of the site investi-
gation process, on the financial risk of construction projects. Another goal is to compare the cost of extra
site investigation with the repairing or reconstruction cost result from improper site investigation. The
results of limited site investigation scope are clearly shown in a recycling factory at the 6th of October city
governorate, Egypt, which suffer from cost overrun and long delays. The apparent savings due to making
inadequate site investigation leads to cost overrun by about 64.2% of the project cost.
331
332
333
334
335
336
Z. Cao
State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan,
P.R. China
Y. Wang
Department of Civil and Architectural Engineering, Shenzhen Research Institute, City University of Hong Kong,
Hong Kong, China
ABSTRACT: This paper presents a Markov Chain Monte Carlo Simulation (MCMCS)-based approach
for probabilistic characterization of the undrained Young’s modulus, Eu, of clay, which utilizes both
the prior knowledge and project-specific SPT data to generate a large number of equivalent samples
of Eu for its probabilistic characterization. The proposed approach combines the prior knowledge and
project-specific SPT data systematically under a Bayesian framework and allows general choices of
realistic prior knowledge (e.g. an arbitrary histogram type of prior distribution). Equations are derived
for the proposed approach, and a sensitivity study is performed to explore the effects of prior knowledge
on probabilistic characterization of soil properties. It is shown that the proposed equivalent sample
approach integrates the information provided by different types of prior knowledge with project-specific
information in a rational manner and improves significantly probabilistic characterization of soil
properties by incorporating consistent prior knowledge.
337
338
339
340
341
342
343
A.D. Mackay
Nishimatsu Construction Company Limited
Formerly SMEC Asia Ltd.
N.R. Wightman
Sinclair Knight Merz (HK) Limited
Formerly SMEC Asia Ltd.
ABSTRACT: Due to the current resource demand from the People’s Republic of China and the
availability of these resources in the Gobi desert, Mongolia, there is a need to construct a railway to
provide an efficient haulage system from the Gobi desert mines to the People’s Republic of China. The
proposed railway currently runs for about 225 kilometres within Mongolia from Ukhuu Khudag to the
People’s Republic of China border crossing facilities at Gashun Sukhait. As the Oyo Tolgoi open cast
mine will also come into operation in the near future provision for haulage of these mined resources needs
to be made; the alignment passes in close proximity to this location accordingly. As part of the feasibil-
ity design a preliminary site investigation was carried out, including a limited desk study, a walkover
survey, a topographic survey and a basic ground investigation, including boreholes and trial pit explo-
ration stations. The ground investigation stations were located according to geotechnical construction
risks assessed from the desk study, topographic survey and walkover survey information; at this pre-
liminary stage these included stream and river crossings, areas with considerable volumes of cut and fill,
material suitability for ballast and embankment construction and areas with near surface groundwater.
The ground conditions along the alignment included superficial deposits; alluvium, cemented alluvium,
talluvium and aeolian (wind-blown) deposits; and solid geology; sedimentary (coal measures, shale and
sandstone) and igneous (tuff, basalt, andesite and granite) rock. This paper outlines the methods adopted
for carrying out the site investigation, which included considerations for the past environmental and geo-
logical history when assessing the physical characteristics relevant to construction.
345
346
3 GROUND CONDITIONS
347
3.4 Findings
The major geological units, including areas of
exposed solid ground beneath the superficial
deposits, generally corresponded with the terrain
evaluation as follows:
Aeolian deposits—generally covering the entire
area with localized increases in thickness in zones
of “small hillocks” and in areas of flatter land
Talluvium—typically standing at shallow gradi-
ents beneath the rock outcrops
Alluvium and lake deposits—extensive areas of
flat ground typically in the vicinity of major rivers
Sedimentary rock (sandstone, shale and coal
measures)—typically undulating terrain
Igneous rock (dyke intrusions; basalt and
andestite)—associated with abrupt changes in
topography and surface exposure
Igneous rock (batholitic intrusions; granite)—
“tor” landscape and surface exposure
The general findings of the SI, are presented in
Table 2 and Figure 7.
4 INTERPRETATION
4.1 Weathering
The relative effects of the extent of weathering in
Figure 6. Terrain evaluation (a) and Plates 1–3 (b–d). different environments are presented in Figure 8
(Fookes 1997).
As shown weathering in desert environments,
The BHs were advanced using undisturbed such as the Gobi desert, is minimal. Relatively com-
sampling and Standard Penetration Testing tech- petent ground with an increased bearing capacity
niques and bulk samples and in-situ testing within can therefore be expected near surface. As presented
the pits. Subsequent laboratory testing included in Plate 5 the most recent deposits, referenced Ap-La
particle size distribution, atterberg limit and Q111-N (Table 2), located further to the north-west,
aggregate testing. The GI specification followed showed the weaker superficial deposits to have a
international British and Eurocode standards (BS minimal thickness and be present over more compe-
1999) as required. tent rock, as shown by the difficulty in excavation.
348
Recent (66 Ma to now) Alluvium/lake deposits (gravel, clay and sand (key Ap-La Q111-N); talluvium;
alluvium/lake deposits (gravel, clay and sand (key L-aN2) and alluvium/lake
deposits (cemented granular sand and clay). (key P3).
Cretaceous (66–145 Ma) Lake-alluvial deposits (key LaK2) and coal measures (“aleurolite”—siltstone
and coal measures)
Jurassic (145–200 Ma) Intrusive rocks (basalt, andesite and granite)
Carboniferous (300–350 Ma) Shale with basalt and andesite intrusions
Devonian (350–417 Ma) Sandstone with tuff and basalt and andesite intrusions
349
5 CONCLUSIONS
ACKNOWLEDGEMENTS
REFERENCES
Baynes, F.J., Fookes, P.G. & Kennedy, J.F. 2005. The total
engineering geology approach applied to railways in
the Pilbara, Western Australia. Bulletin of Engineering
Geology and the Environment 64: 67–94.
British Standards (BS), 1999. Code of practice for site
investigations. London: British Standards Institu-
tion (BSI).
Fookes, P.G. 1997. The 1st glossop lecture, geology
for engineers: the geological model, prediction and
performance. Quaterly Journal of Engineering Geology
30(4): 293–424.
Fookes, P.G., Hearn, G.J. & Tollin, B. 2005. Aligning oil
Figure 9. Geotechnical properties. and gas pipelines through the Makharov Mountains,
Sakhalin Islands, Russia. In M. Sweeney (ed.), Int.
Conf. on Terrain and Geo-hazards Challenges Facing
On-shore oil and gas pipelines, London, 2–4 June
uncemented was typically present at depth beneath 2004.
the duri-crust. Refer to Plate 7 on Figure 10. Mackay, A.D. 2010. The role of the civil engineer in the
The thickness of the duricrust and its vicinity to changing demand for space in HK. In HKIE Civil
ephemeral stream course was an important consid- Division Conference 2010: infrastructure solutions for
eration for estimating the bearing capacity. tomorrow, Hong Kong, 12–14 April 2010.
350
ABSTRACT: Monitoring is an essential element of modern tunneling construction. The most common
monitoring method is measuring displacements, for example convergence of the tunnel opening or surface
settlements. Measurement outcomes can be used to update the knowledge on material properties of the soil
or other parameters that enter numerical models of the structural behavior of the tunnel. In probability
theory, this process can be formalized in the concept of Bayesian updating. In this paper, we apply the
Bayesian concept to update the numerical model of a tunnel in soft soil conditional on settlement meas-
urements. The tunnel is constructed by means of the conventional tunneling method and modeled with 2D
finite elements applying the stress reduction method. We assume that settlement measurements are taken
at full excavation and utilize the measurements to update the material properties of the soil as well as the
the relaxation factor of the stress reduction method. Updating is performed by means of BUS, a recently
proposed method for Bayesian updating of mechanical models with structural reliability methods.
351
352
L f( )
f Z) = (4)
∫R n L f ( )dx
353
354
355
Eoed, Eud and the relaxation factor β. In this case, that models the 3D arching effect of the stress
which implies higher information content of the distribution.
measurement, the posterior mean of β is further
decreased and the weak zone around the tunnel
with low values of Eoed and Eud is increased. REFERENCES
Figure 6 demonstrates the influence of the prior
knowledge on the relaxation factor β. Therein, the Au, S.K. & Beck, J.L. 2001. Estimation of small failure
prior and posterior PDFs of β are plotted for two probabilities in high dimensions by subset simulation.
different assumed prior coefficients of variation Probabilist. Eng. Mech. 16(4): 262–277.
(10% and 20%) and for σ∈ = 1 mm. It is observed Der Kiureghian, A. 2005. First- and second-order reli-
ability methods. In E. Nikolaidis, D.M. Ghiocel,
that as the prior knowledge on beta decreases, i.e. and S. Singhal (Eds.), Engineering Design Reliability
as its prior coefficient of variation increases, the Handbook. Boca Raton, FL: CRC Press.
influence of the measurement becomes higher. Der Kiureghian, A. & Liu, P.L. 1986. Structural reliabil-
Comparing the posterior PDFs for the two cases, ity under incomplete probability information. J. Eng.
we see that the same measurement information Mech.-ASCE 112(1): 85–104.
leads to much lower values of β when a larger Gelman, A. (2004). Bayesian data analysis. Boca Raton,
prior coefficient of variation is assumed. This FL: Chapman & Hall/CRC.
result further highlights the influence of the 2D Ghanem, R. & Spanos, P. 1991. Stochastic Finite
modeling of the arching effect on the surface set- Elements—A Spectral Approach. Berlin: Springer.
Gilks, W.R., S. Richarson & D.J. Spiegelhalter (Eds.)
tlements. Moreover, it shows how the confidence 1998. Markov chain Monte Carlo in practice. Boca
on the prior assumption can influence the updat- Raton, FL: Chapman & Hall/CRC.
ing results that may provide a basis for further risk Möler, S. 2006. Tunnel induced settlements and structural
and reliability assessment. forces in linings. Ph.D. thesis, Institut für Geotechnik,
Universität Stuttgart.
Panet, M. & Guenot, A. 1982. Analysis of convergence
5 CONCLUSION behind the face of a tunnel. In Proc. International
Conference on Tunnelling ’82, London. IMM.
In this paper, we performed Bayesian updating Papaioannou, I., Betz W., Zwirglmaier K. & Straub D. 2013.
MCMC algorithms for subset simulation. Manuscript,
of the parameters of a 2D numerical model of a Engineering Risk Analysis Group, TU München.
tunnel in soft soil, conditional on settlement meas- Papaioannou, I. & Straub D. 2012. Reliability updating
urements. We applied BUS, a recently proposed in geotechnical engineering including spatial variabil-
method for Bayesian updating with structural reli- ity of soil. Comput. Geotech. 42: 44–51.
ability methods, combined with subset simulation, Ranjan, R., Betz W., Papaioannou I. & Straub D. 2013.
an adaptive Monte Carlo method that is able to A two-step approach for reliability assessment of a
handle efficiently problems with a large number tunnel in soft soil. In Proc. 3rd International Confer-
of random variables. The results demonstrate the ence on Computational Methods in Tunnelling and Sub-
influence of the accuracy of the measurement surface Engineering EUROTUN: 2013, Bochum.
SOFiSTiK AG. 2012. SOFiSTiK analysis programs ver-
device as well as the prior knowledge of the uncer- sion 2012. Oberschleißheim: SOFiSTiK AG.
tain parameters on their posterior distributions. It Straub, D. 2011. Reliability updating with equality infor-
was shown that the highest impact of the measure- mation. Probabilist. Eng. Mech. 26(2): 254–258.
ment fell on the stiffness moduli and the relaxa- Straub, D. & Papaioannou I. (under review 2013). Bayesian
tion parameter of the stress reduction method updating of mechanical models. J. Eng. Mech.-ASCE.
356
Satoko Ryuo
Railway Technical Research Institute, Tokyo, Japan
Yoshiharu Asaka
Shimizu Corporation, Tokyo, Japan
Kenichi Soga
University of Cambridge, Cambridge, UK
ABSTRACT: In typical conventional foundation design, the inherent variability of soil properties,
model uncertainty and construction variability are not modeled explicitly. A main drawback of this is that
the effect of each variability on the probability of an unfavorable event cannot be evaluated quantitatively.
In this paper, a method to evaluate the uncertainty-reduction effect on the performance of a vertically-
loaded pile foundation by monitoring the pile performance (such as pile load testing or placing sensors in
piles) is proposed. The effectiveness of the proposed method is examined based on the investigation of a
120-pile foundation placed on three different ground profiles. The computed results show the capability
of evaluating the uncertainty-reduction effect on the performance of a pile foundation by monitoring.
357
Weight density 0–0.15 Harr (1984), Lacasse and Nadim (1997), Phoon
and Kulhawy (1999a)
Undrained shear strength 0.1–0.70 Harr (1984), Cherubini (1997), Lacasse and Nadim
(UU, UC) (1997), Phoon and Kulhawy (1999a, 1999b),
Duncan (2000)
Internal friction angle 0.02–0.15 Harr (1984), Phoon and Kulhawy (1999a, 1999b)
SPT N-value 0.15–0.5 Harr (1984), Barker et al. (1991), Phoon and
Kulhawy (1999a)
Shear modulus or 0.15–0.65 Phoon and Kulhawy (1999a)
Young’s modulus
DMT Horizontal stress index 0.2–0.6 Phoon and Kulhawy (1999a)
Compression index 0.1–0.37 Harr (1984), Duncan (2000)
358
of the ultimate base resistance is also 0.42. These 1. First, using the borehole data, the thickness of soil
values of uncertainty approximately coincide with layers and the soil parameter values are estimated
the variation coefficients of resistance calculated by the Kriging method at all pile locations.
from field tests in previous studies (ex. Barker et al. 2. At borehole locations, the ultimate shaft resist-
(1991), Okahara et al. (1993), Kay (1993), Archi- ance and ultimate base resistance of a “virtual”
tectural Institute of Japan (AIC) (2001) ). Hence pile are calculated. Also, based on the thickness
it is concluded that the developed procedure for of soil layers and the soil parameter values eval-
calculating uncertainty is reasonable. uated from Step (1) above, the ultimate shaft
For a pile that is directly monitored by instru- resistance and ultimate base resistance of each
mentation or if a pile load test was carried out, the pile in the actual foundation are calculated by
uncertainty of the “monitored” pile is governed by using each borehole data.
the accuracy of measurement devices, the parame- 3. A pile or piles that will be monitored are ran-
ters used in conversion of measurement and physi- domly selected as a given monitoring scenario.
cal quantity such as Young’s modulus of concrete The mean value of the ultimate shaft resist-
or steel. In this paper, for monitored piles, the vari- ance and base resistance obtained from Step
ation coefficients of ultimate shaft resistance and (2) has a variability characterized by a normal
ultimate base resistance are assumed as 0.05 (see distribution and the 95% confidence bounds
Table 3). These values are smaller than the values is set as 0.05 for monitored piles and 0.42
assigned for non-monitored piles, which are given for “virtual” piles. Then, 1000 data sets of
in Table 2. shaft resistance and base resistance are ran-
domly produced for both monitored piles and
“virtual” piles.
2.2 Uncerainty of a whole pile group foundation
4. From each of 1000 data sets, the ultimate shaft
The uncertainty of a whole pile group foundation resistance and base resistance for each pile in the
that consists of normal and monitored piles can be actual foundation are estimated by the Kriging
evaluated by the following steps. method.
359
The estimated value is called “average of varia- The results for the cases with one monitored pile are
tion co-efficient” in this paper. As a separate study, shown in Figures 3 to 5 for the three different soil
other evaluation criteria for selecting the moni- profiles. In these figures, the effectiveness of moni-
toring point(s) to reduce the uncertainty of piled toring for the whole piled foundation is defined
foundation are considered and the results will be by the difference between the average of variation
reported elsewhere in the future. coefficient for whole pile foundation calculated
in one monitored pile situation and that without
monitoring. The differential value is computed at
each monitored pile location and then a contour
3 CASE STUDY
plot showing the spatial distribution of the differ-
ential value is generated. The location that gives the
3.1 Condition for calculation
maximum value is considered to be the ideal place
A piled foundation which has 120 single piles to monitor. Also the maximum value itself gives
placed at 7-meter intervals is considered as a hypo- the value of monitoring; that is, the larger the value
thetical case study. The plan area is composed of is, the uncertainty reduces more by monitoring.
7-span-49 meter × 14-span-98 meter as shown in Figure 3 shows the result for Case 1. In the case
Figure 1. Each pile is cast-in-place pile with diam- of horizontally-layered conditions, the pile located
eter of 0.8 m and its length coincides with the thick- away from the borehole, especially the one near the
ness of clay as shown in Tables 4 to 6. The ultimate center of the piled foundation, will be ideal in terms
resistance is given by the shaft resistance and the base of reducing the uncertainty of the foundation.
resistance considering the soil conditions. Three soil The result for Case 2 is shown in Figure 4. It indi-
profiles are considered; (i) Case 1: two horizontally cates that a pile located away from the boreholes
layered conditions, (ii) Case 2: two layers with an and on the left hand side of the ground that has a
inclined interface, in which the slope angle changes moderate interface slope angle would be effective
at the centerline of the foundation, and (iii) Case 3: one for monitoring to reduce the uncertainty.
360
cu of clay (kPa) 50 50 55 55 50 50 45 45
qb at toe (kPa) 12000 11700 11700 12300 12000 12300 12000 12000
Thickness of clay = pile length (m) 20 19.5 19.5 20.5 20 20.5 20 20
Diameter of model pile, D (m) 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8
Adhesion factor of shaft resistance, α 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6
Ultimate shaft resistance, USR1 (kN) 1508.0 1470.3 1617.3 1700.2 1508.0 1545.7 1357.2 1357.2
Ultimate base resistance, UBR (kN) 6031.9 5881.1 5881.1 6182.6 6031.9 6182.6 6031.9 6031.9
USR1 + UBR(kN) 7539.8 7351.3 7498.3 7882.9 7539.8 7728.3 7389.0 7389.0
cu of clay (kPa) 50 50 55 55 50 50 45 45
qb at toe (kPa) 10200 9900 10800 12300 12000 13800 15000 15000
Thickness of clay = pile length (m) 17 16.5 18 20.5 20 23 25 25
Diameter of model pile, D (m) 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8
Adhesion factor of shaft resistance, α 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6
Ultimate shaft resistance, USR1 (kN) 1281.8 1244.1 1492.9 1700.2 1508 1734.2 1696.5 1696.5
Ultimate base resistance, UBR (kN) 5127.1 4976.3 5428.7 6182.6 6031.9 6936.6 7539.8 7539.8
USR1 + UBR(kN) 6408.8 6220.3 6921.6 7882.9 7539.8 8670.8 9236.3 9236.3
cu of clay (kPa) 50 50 55 55 50 50 45 45
qb at toe (kPa) 12000 11700 11700 12300 12000 17700 18000 18300
Thickness of clay = Pile length (m) 20 19.5 19.5 20.5 20 29.5 30 30.5
Diameter of model pile, D (m) 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8
Adhesion factor of shaft resistance, α 0.6 0.6 0.6 0.6 0.6 0.6 0.6 0.6
Ultimate shaft resistance 1, USR1 (kN) 1508 1470.3 1617.3 1700.2 1508 2224.2 2035.8 2069.7
Ultimate base resistance, UBR (kN) 6031.9 5881.1 5881.1 6182.6 6031.9 8897 9047.8 9198.6
USR1 + UBR (kN) 7539.8 7351.3 7498.3 7882.9 7539.8 11121.2 11083.5 11268.3
In Case 3, the result is given in Figure 5. Results new information about ground profile is found
show that monitoring in the space enclosed by to be similar to the existing ground profile, this
BH6, BH7 and BH8 where the clayey soil layer is reduces the uncertainty in the ground profile. This
relatively thicker will not be effective in reducing would increase the precision of fitting function
the uncertainty. The most effective location for to the ground profile, which in turn reduces the
monitoring is at the left hand side of the ground uncertainty of the whole foundation.
in which the soil profile is uniform. Another effec-
tive monitoring area is the space enclosed by BH4,
3.3 Computed results for multi-monitored pile
BH5 and BH6 where the thickness of the clayey
soil layer is variable. Based on the contour plots, the location(s) of
The efficient locations obtained from the calcu- the monitoring are determined and the effect of
lations are consistent with typical rules of thumb. monitoring number in reducing the uncertainty
According to Case 2 and 3, the effective area for was examined. The number of monitoring was
monitoring is located at the ground that has a increased up to ten and a relationship between
moderate interface slope angle and that has uni- monitoring number and the effect is shown in
form soil profile in large part. This is because when Figure 6. The average of variation coefficient for
361
4 CONCLUSIONS
REFERENCES
362
363
T. Schweckendiek
Deltares, Unit Geo-Engineering, Delft, The Netherlands
Department of Hydraulic Engineering, Delft University of Technology, Delft, The Netherlands
ABSTRACT: Backward erosion piping is an important internal erosion failure mechanism in deltaic
areas. Piping is very sensitive to ground conditions, especially to the permeability of aquifers under levees.
The uncertainties in these ground properties can be reduced by monitoring the hydraulic head in the piping-
sensitive aquifer in order to update the probability distributions of the basic random variables influencing
the geo-hydrological conditions such as the hydraulic resistance of the aquifer. This paper demonstrates how
Bayesian posterior analysis can be used to incorporate head monitoring information in order to update the
reliability with respect to uplift, heave and piping for river levees. The paper starts by recapping background
on uplift, heave and piping reliability, before we discuss how the equality type of information involved can
be used in posterior analysis by following the approach proposed by Straub (2011). The workings of the
theory are illustrated by an example considering uplift (of a blanket layer) only using a simplified ground-
water model for the sake of illustration. Sensitivity analysis on the example shows that the effect of reliability
updating can be considerable for large prior uncertainties and for observations which are rather unlikely
given the prior distributions of the random variables (i.e. far from the expected value of the measurement).
365
366
inequality domain:
where X is the vector of random variables and fX(x)
{ }
its (prior) Probability Density Function (PDF).
Fragility curves are a frequently used alternative εe ⎡ fe (φm − φ ( ))⎤⎦ < (10)
⎣ m
way of expressing the uncertainty in the aggre-
gated resistance R with respect to a dominant load where u is a standard Normal distributed random
variable. In other words, R contains all variables variable, Φ −1 is the inverse cumulative standard
in X except the dominant load variable, which for Normal distribution function, fem is the PDF of
river levees is the water level h), providing the con- the measurement error and c is a positive constant
ditional probability of failure with respect to h: to be chosen such that the term in square brackets
is always between zero and one. For the deriva-
P ( Fi | h ) ∫g ( f ( r))dr (6) tion and further details of this method we refer to
i h) < 0 R
Straub (2011).
As explained in the previous section, failure This transformation allows us to compute the
can only occur if uplift precedes heave and heave updated probability of failure using standard reli-
precedes piping, because piping can only develop ability analysis methods by evaluating the follow-
if the groundwater can flow upwards through the ing expression:
ruptured blanket (i.e. uplift) sufficiently to move
sand particles vertically (i.e. heave). In reliability P (F ∩ ε e ) ∫F ∩ε fX ϕ u )du
d dx
terms this relationship can be described by means P (F | ε ) = = e
(11)
P (ε e )
of a parallel system, for which the failure set is ∫ε e
fX ϕ u )du dx
given by:
where ϕ(u) is the standard normal PDF. Notice
F {gu < } ∩ {gh < } ∩ {g p < } (7) that this way the posterior probabilities of
failure for the separate failure mechanisms
Note that the prior probability of failure and the uplift, heave and piping can be computed as well
fragility curves as defined above can be determined as the posterior (parallel system) probability
per failure mode individually (i.e. Fu, Fh and Fp) as of failure for the whole failure mode (i.e.
well as for the “system” failure event F. P(Fu∩Fh ∩Fp|ε)).
3 RELIABILITY UPDATING
4 EXAMPLE: UPLIFT
This section provides the theory and definitions
4.1 Limit state
used in this study for reliability updating with infor-
mation from pore pressure monitoring. We apply In the following we provide a numerical elabora-
Bayes’ rule (Bayes, 1763) to obtain the posterior tion for uplift failure as defined in section 2.2. The
(updated) probability of failure, here expressed in limit state function for this example is obtained by
367
⎛ γ ssat ⎞
hc,u mu d /λ hp (13)
⎝ γw ⎠
368
5 CONCLUSION
369
370
ABSTRACT: The capability of predicting Blasting Induced Damage (BID) in the surrounding rock mass
is important to excavation of high rock slopes. A tensile-compress damage model linked to the computer
code LS-DYNA through the user subroutine interface was used by introducing the consideration of com-
press damage and modifying the method of definition of the elastic constants based on the existed blasting
damage models. To verify the rationality of the proposed model, based on the case study of the excavation
of high rock slope at the Xiluodu Hydropower Station in Sichuan province of China, spatial distribution
of blasting induced damage zone at the 640 m berm were determined with cross-hole sonic test. Numerical
simulations of BID at the 640 m berm were implemented with tensile-compress model and other four
existing blasting damage models. The damage zone obtained by the modified tensile-compressive damage
model agreed with observations better than other four existing blasting damage models.
373
374
E E ( − D) (12)
3 SONIC WAVE TEST AND BLASTING
G G ( − D) (13) INDUCED DAMAGE DISTRIBUTION
The above method of determining the elastic To fully understand the damage distribution in
constants could not meet the mathematic relation the remaining rock mass near the berm, the cross-
between the elastic constants. Normally, the elastic hole sonic wave test was employed to determine
constants have the following relations: the extent of damage zone surrounding the berm
at the 640 m elevation of the Xiluodu high rock
E 2(1 μ )G (14) slope. The arrangement of the testing holes was
E 3(1 2 μ )K (15) shown in Figure 1. Three groups measuring results
were obtained in this test. One group of them was
The method of determining the elastic constants in the left and other two groups were in the right.
was modified in this study, the degraded bulk mod- The change of acoustic velocities and the damage
ulus and Poisson’s ratio were definited as the same depths are listed in Table 1.
as the origin method. According to the results of sonic wave test, the
acoustic velocities in the remaining rock mass near
K K ( − D) (16) the berm range from 1500 m/s to 5800 m/s. The
acoustic velocities decreases form the outer flank
⎛ 16 ⎞ to the inside of the berm. The minimum acoustic
μ μ exp − βCd (17)
⎝ 9 ⎠ velocity was found at the outer flank of the berm,
which means the rock mass in this part is the most
Other elastic constants were determined by the serious damaged.
mathematic relation: The outline of EDZ for rock mass near the
berm was showed in Figure 2. In the vertical
3(1 − 2 ) direction the damage depth ranges from 1.0 m to
G K (18)
2(1 + μ ) 4.0 m and decreases nonlinearity from the outer
flank to inside of the berm. In the horizon direction
9KG
E= (19)
3K G
Based on the above analysis, a tensile-compress Figure 1. Schematic illustration of sonic wave test holes
damage model was established. LS-DYNA is a at the EL.640 m berm.
375
the damage radius is between 1.2 m and 3.0 m, and tensile-compress damage model for the rock mass
decreased nonlinearity from the top to the bottom of basalt were listed in Table 2. The material con-
of slope. The maximum vertical damage depth and stants such as elastic modulus, Poisson’s ratio et al
the maximum horizontal damage extent are found were determined according to the geological data
to be at the outer flank of the berm. and field experiments. To simulate borehole blast-
It can be seen that damage degree and extent of ing, a radial pressure was applied to the elements on
rock mass near the berm is more serious and much the borehole wall. The amplitude and the duration
larger than other parts of slope. There are three of the pressure wave can be determined according
induced factors for this: firstly, the rock mass near to the method proposed by Lu et al. (2011), the
the outer flank of the berm was affected by buffer density of explosive is 1000 kg/m3 and the detona-
blasting of former bench most seriously; secondly, tion velocity is 3200 m/s. To prevent the effects of
there are two free surfaces at the top of slope, reflected waves at the edge of the region, the trans-
which the enlarge degree of stress redistribution mitting boundaries were adopted in numerical sim-
and induce a more larger damage extent; thirdly, in ulation. Damping in the model is done through the
the current bench blasting excavation, the combina- Rayleigh classical approach by making the damp-
tion of effect the shock wave and detonation gas is ing matrix equal a linear combination of the mass
another main damage induced factor. To guarantee and stiffness matrix.
the safety of the slope, the reinforcement measure
should be carried out on the top of slope in time. 4.2 Numerical results
Figure 3 shows the damage distribution character-
istics of the former bench and presplitting blasting
4 NUMERICAL SIMULATION OF
of the current bench.
BLASTING INDUCED DAMAGE ZONE
It can be seen that the damage depth caused by
the former excavation bench is about 0.6 m. The
4.1 Numerical model
damage degree decreases from the outer flank to
In the present study, 3D numerical simulation of the inner side of the berm. After the presplitting
blasting induced damage was implemented with the blasting of underlayer excavation, a columnar
computer code LS-DYNA. The height of blasting damage zone was formed, the degree of dam-
excavation layer is 15 m depth. Diameter of the blast age decreases with the distance to the blast hole.
holes is 110 mm. The excavation region can be divided The rock mass close to the blast hole which was
into many same sections. One section contains five crushed by the shock wave was critical destroyed
presplit holes, two buffering holes and one produc- for the damage scalar equals 1.0, then the shock
tion hole. Considering the limit of mesh number and wave changed into stress wave and the damage
computation time, one section was selected to estab- degree decreases to about between 0.2 and 0.8. The
lish the numerical model. The number of nodes is maximum damage extent is located in the top face
566016 and the element number is 541400. of the blast hole.
Assuming to be isotropic and homogenous To verify the accuracy of the tensile-compress
material, the damage constants used in the blasting damage model, other four exiting blasting
376
377
378
ABSTRACT: Reliability analysis of slope stability has received considerable attention. This paper aims
to propose a Non-Intrusive Stochastic Finite Element Method (NISFEM) for slope reliability analy-
sis with an implicit performance function. This procedure is formulated on the basis of the Stochastic
Response Surface Method (SRSM) and the deterministic finite element method. The probabilistic analy-
sis of slope stability is decoupled with the deterministic finite element analysis. A global sensitivity analy-
sis based on Sobol’s indices is carried out to identify the key factors influencing the factor of safety of a
slope. A practical rock slope example is presented to demonstrate the capability of the proposed method.
The results indicate that the proposed method can effectively evaluate the slope reliability with an implicit
performance function with a sufficient accuracy. The internal friction angles of the fault f42-9 and the class
III2 rock mass have a significant effect on the stability of the Jinping left abutment slope.
379
1. Identify the random variables and estimate their + ⋅ ⋅ ⋅ + ∑ ∑ ∑ ⋅ ⋅ ⋅ ∑ ai1 i2 , …,, in Γn ξi , ξi2 , …, ξin
i1 = 1 i2 = 1 i3 = 1 in = 1
means, Coefficients Of Variation (COV), distri-
bution types and correlations between them for (1)
the considered slope reliability problem.
2. Construct the FEA model of slope stability with where FS(ξ) is the vector of the factor of safety;
the mean values of the input random variables a = (a0, ai1, ai1i2, ⋅⋅⋅) are the unknown coefficients
using the software SIGMA/W and SLOPE/W. to be estimated, in which n is the number of ran-
dom variables representing the uncertainties of
the model inputs, and i1 + i2 + ⋅⋅⋅ + in ≤ n; ξi = (ξi1,
ξi2, ⋅⋅⋅, ξin) is the vector of independent stand-
ard normal variables resulting from the trans-
formation of the input uncertain parameters,
Xi; Γn(ξi1, ξi2, ⋅⋅⋅, ξin) is the multivariate Hermite
polynomials of degree n. (Li et al. 2011, Mollon
et al. 2011).
4. Generate and select Nc sets of collocation
points in the independent standard normal U
space, ξ = (ξ1, ξ2, …, ξNc) based on the linearly
independent principle. The linearly independent
Probabilistic Collocation Method (PCM) can
yield accurate results with much less computa-
tional cost. The number of runs of slope stability
model equals to the number of unknown coef-
ficients in the PCE of order p, Nc, Nc = (n + p)!/
(n! × p!). The reader is referred to Mao et al.
(2012) for details.
5. Map the selected collocation points ξ in the
independent standard normal U space into the
collocation points X in the physical space. If
the input random variables are correlated non-
normal variables, the Nataf transformation
(Nataf 1962) can be used to transform ξ into X.
Then, the Nc sets of collocation points X can be
obtained and taken as input parameters for the
deterministic FEA model of slope stability.
6. Replace the corresponding mean values of
input random variables in the source FEM-FS
*.xml file generated in step (2) with the Nc sets
of the collocation points X. Then Nc different
new FEM-FS *.xml input files can be gener-
Figure 1. Flowchart of the NISFEM for slope reliabil- ated. Thus, no programming effort is required
ity analysis. to modify the existing finite element codes in
380
381
Table 2. Statistics of random variables in the Jinping as well as the topography of the computational
left abutment slope stability model. domain, the Finite Element Method (FEM) (Zou
et al. 1995, Kim & Lee 1997) is used to calculate the
Random factor of safety of the left abutment slope because
Materials variables Mean COV
this approach can account for the influence of the
Lamprophyre dike X c1 (kPa) 20 0.25 deformation of rock masses on the slope stability.
φ1 (°) 16.7 0.15 Theoretically, a Three-Dimensional (3D) model
Fault f42-9 c2 (kPa) 20 0.3 for the slope stability problem considering spatial
φ2 (°) 16.7 0.2 variability of rock masses should be employed.
Class III2 rock mass c3 (kPa) 900 0.15 For simplicity, a finite element model based on
φ3 (°) 45.57 0.08 2D plane strain analysis is built in the SIGMA/W
Class IV1 rock mass c4 (kPa) 600 0.18 with the means of random variables, as shown in
φ4 (°) 34.99 0.10 Figure 3. The analyzed domain is 420 m in height
Class IV2 rock mass c5 (kPa) 400 0.2 from the elevation 1650 m to the elevation 2070
φ5 (°) 30.96 0.12 m and 384.7 m in length of the direction perpen-
dicular to the river from the centerline of rive to
the within slope. The gravity loads induced by the
self-weight of rock masses are applied. For illustra-
friction angles of five types of materials, respec- tive purposes, a conventional elastic and perfectly
tively. The means of the shear strength parameters plastic model based on the Mohr-Coulomb failure
are determined based on field test, laboratory test criterion is adopted to represent the stress-strain
supplemented with engineering judgment, such as behavior of rock masses and structural planes. The
direct shear test and triaxial test (Wu et al. 2008). finite element model consists of three-node trian-
The Coefficients of Variation (COV) of the shear gular and four-node quadrilateral hybrid elements.
strength parameters are adopted from the litera- There are totally 14363 zones and 14113 grid points
ture (Phoon & Kulhawy 1999, Baecher & Christian in the mesh. With regard to the boundary condi-
2003, Tang et al. 2012). Additionally, all random tions, the horizontal movement on the vertical left
variables are assumed to be characterized statisti- boundary of the mesh is restrained, whereas the
cally by a lognormal distribution. The slope sta- base of the mesh is not allowed to move in both the
bility is evaluated under two working conditions, horizontal and the vertical directions. The initial
namely natural condition and rainfall condition. stress field is obtained using the SIGMA/W. Then,
For the natural condition, the groundwater table is it is imported into the SLOPE/W for slope stability
below the slip surface of the slope. If the ground- analysis. Based on the geological investigation, it
water table is above the slip surface of the slope, can be expected that a potential deep sliding may
the resulting load should be taken into account, occur along the lamprophyre dike X and fault f42-9,
which corresponds to the rainfall condition. The and shear off superficial rock masses in accord-
pore-water pressure ratio Ru = 0.1, a ratio of the ance to the path highlighted with a red dashed line
pore-water pressure to the weight of overburden in Figure 3. For the left abutment slope with the
rock masses is used to account for the effect of prescribed slip surface, the factors of slope safety
rainfall on the stability of the left abutment slope. obtained from the FEM are 1.083 and 0.981 for
In order to model the structures and the gen- the natural and rainfall conditions, respectively,
eralization of geomechanical model of the slope which are consistent with 1.110 and 1.014 using
382
383
384
ABSTRACT: A novel bond contact model for rocks was proposed, based on a series of laboratory tests
on bonded granules idealized by two glued aluminum rods. The model was then implemented into a two-
dimensional distinct element code to numerically carry out the uniaxial and biaxial compression tests on
Lac du Bonnet granite. The results were compared to the experimental data and the DEM simulations
using an existing Bonded-Particle Model (BPM). The results demonstrate that the proposed model can
reproduce the strength envelope of Lac du Bonnet granite adequately. The slope of the strength envelope
obtained from our simulation is higher than that obtained with the BPM due to higher peak shear and
rolling strength in the proposed model, and that the ratio of shear and rolling failure to tension failure is
also higher.
385
⎧ Rc
Figure 1. Schematic illustration of the mechanical ⎪μb R R
responses of the bond contact model: (a) normal contact ⎪ c t
⎪ ⎡ ⎛ R Rt ⎞ s ⎤
f
model; (b) tangential contact model; (c) rolling contact Rs = ⎨ ⋅ (Fnb Rt ) ⋅ ⎢1 + gs ⋅ ⎜ ln c ⎥ , Fn ≤ Rc ()
model (Jiang et al. 2006, 2012a). ⎪ ⎢⎣ ⎝ Fnb Rt ⎟⎠ ⎦⎥
⎪F , Fn > Rc (2)
⎪ n b
to the residual strength when the rolling bond ⎩
strength Rr is reached.
⎧1 Rc
⎪ 6 β br R R ⋅ ( Fnb Rt )
2.2 Experimental setup ⎪ c t
⎪ ⎡ ⎛ R Rt ⎞ r ⎤
f
In order to validate the contact law of bonded gran- Rr = ⎨ ⋅ ⎢1 + gr ⋅ ⎜ ln c ⎥, ≤ (3)
⎝ Fnb Rt ⎟⎠ ⎥⎦
n c
ules as well as verify the mechanical responses of ⎪ ⎢⎣
the model, Delenne et al. (2004) first presented an ⎪1
experimental investigation on mechanical behavior ⎪ Fn β br , n > c ( 4)
⎩6
of cemented granules by performing simple tests
on a pair of aluminum rods glued together with
where Rt is tensile strength; μb is friction coeffi-
epoxy resin. Taking the normal force and the bond
cient of cement; βb is rolling resistance coefficient
material into account, Jiang et al. (2012a, b) con-
of cement; r = 2r1r2/(r1 + r2) is the common radius
ducted a series of simple and complex loading tests
of two contact particles with radii of r1 and r2;
on a pair of aluminum rods glued with a rock-like
material, cement, in addition to epoxy resin.
The analogous cemented sample adopted in
the experiment is shown in Figure 2. The alumi-
num rods glued together with cement are 12 mm
in diameter and 50 mm in length. The bond of this
sample is 0.6 mm thick and 3 mm wide. Five dif-
ferent loading paths (tension, compression, shear
under different normal forces, rolling under dif-
ferent normal forces and shear-rolling tests under
different normal forces) were applied to charac- Figure 2. Analogous cemented sample used in experi-
terize the mechanical behavior of the analogous ments (Jiang et al. 2012a, b).
386
⎧ Rc
⎪ μb R R ⋅ ( Fnb Rt )
⎪ c t
⎪⎪ ⎡ ⎛ R Rt ⎞ s ⎤
f
Rs = ⎨ ⋅ ⎢1 + gs ⋅ ⎜ ln c ⎥ + μ p Fnp Fn Rccbr (7 )
⎝ Fnb Rt ⎟⎠ ⎥⎦
np
⎪ ⎢⎣
⎪
⎪ μ p Fnp
n b Fnnb , Fn Rccbr (8)
⎪⎩
⎧1 Rc
⎪ 6 β b r R R ⋅ ( Fnb Rt )
⎪ c t
Figure 3. Strength envelopes for cemented samples with ⎪ ⎡1 + gr ⎤
bond thickness of 0.6 mm in terms of: (a) shear strength ⎪ ⎢ f ⎥ 1
Rr = ⎨ ⋅ ⎢ ⎛ Rc Rt ⎞ r ⎥ + β p Fnp , 0 < n ≤ cbr (9 )
vs. normal force; (b) rolling strength vs. normal force. ⎪ ⎢⋅ ⎜⎝ ln F ⎟ ⎥ 6
⎪ ⎣ nb Rt ⎠ ⎦
⎪1 1
⎪ 6 β p np + 6 ⋅ β b rFnnb
and fs, gs, fr, gr are fitting parameters that control , Fn > Rccbr (10 )
the envelope shapes, whose values are 0.986, 2.15, ⎩
0.761 and 3.055 respectively.
387
388
389
390
ABSTRACT: By employ physical experiment and numerical simulation, this study focusses on cracking
risk, the cracking stability and the reinforcement effectiveness of the high and steep left hand bank of the
Jinping dam. The study results relating to the major slope at Jinping showed that: (1) The bending and
toppling cracks in the unreinforced zone of the slope as well as strong relaxation and tension cracks were
extended simultaneously. The rock cracking and slope failure controlling factors could be directly seen in
the structurally weak rock mass. (2) The controlling failure mode for the Jinping steep and high left bank
slope consists of the integral catastrophic instability mode, whereby the slope energy is totally dissipated
and the slope destabilizes at a limit state. The slope cracking factor of safety and the integral stability fac-
tor of safety satisfy the stability requirements. (3) The reinforcement has been effective in keeping cracks
closed in the rock, maintaining the slope stability, and increasing the slope’s inherent safety factor.
391
2.1 Experimental method and model 2.2 Numerical method and model
The physical experiments were performed at FEM (Finite Element Method) was used to simu-
the State Key Laboratory of Hydroscience and late the stability of the high and steep left hand
Hydraulic Engineering at Tsinghua University. It bank of the Jinping dam in this paper. Generally,
consists of a rotational structure system, a defor- strength reduction FEM is often applied to analyze
mation monitoring, cracking and data collection the slope stability numerically recent years. But the
system and an analysis system (Zhou 2008). The strength reduction FEM is not suitable to com-
slope model is built on a rotating support system, pare with the physical experiments because it is not
and can be rotated to any angle of θ around the axis easy to operate experiments with reduced mate-
OY by hydraulic lift (as shown in Fig. 2). Rotation rial properties or parameters. In order to compare
of the test platform at a velocity of 0.5°/min. When the numerical modeling method with the physical
the angle reaches 5°, the excavation of the abut- experiments, the same technique is used in numeri-
ment groove in the model and recording of the cal model, just like what had been done in physical
slope deformation of every measuring point is then modeling. That is, rotating the numerical model to
conducted. The testing platform is continually some angle and calculating the stress condition of
rotated at a velocity of 0.5°/min, until the cracks the whole model until slip surfaces occur.
are initiated on the slope, and then coalescence is Also the following equation is utilized to cal-
propagated until ultimate failure occurs. culate the slope stability factor. This traditional
The simulation scale and faults are shown in method can be used to compare with the method
Figure 1. The model range extends 180 m upstream, developed in this paper.
270 m downstream, and is 600 m wide, taking in
left bank. The simulated elevations and heights are Fanti sliding
K= (1)
1600~2250 m and 660 m, respectively, and sufficient Fsliding
392
(ci Fi + fi Ni ) + Pi
K′ =
Pi ( i i ) Wi sin α i
[ci Fi + fi Wi cos α i + Pi − i ( i −1 − α i ))] + Pi
=
Pi − ( i − i − ) + Wi sin α i
(2) Figure 4. Numerical model of Jinping left bank slope.
393
394
395
396
ABSTRACT: Reference to the principle of employing centrifugal force field to simulate gravity field, the
geomechanics magnetic model test is a new test method that simulates gravity field with electromagnetic
force field to study geomechanical engineering problems. Using electromagnetic method generates a mag-
netic field where the magnetic flux density gradient keeps a fixed nonzero value, where the mixture of fer-
romagnetic material and geological mechanics prototype material can get a uniform electromagnetic force.
The FEM software for 3D electromagnetic field numerical simulation is used to study the influence of the
electromagnetic parameters and magnetic factors to the magnetic field intensity, which provides reason-
able design parameters for the production of the experimental machine. The results of simulation indicate
that the geomechanics magnetic model method can effectively compensate for the landslide weight loss
caused by the reduced-scale and nicely reflect the characteristics of the deformation and failure.
397
398
∂B
∂α
(
= n × 0.044 T ⋅ m −1 ) (3)
3 IMPLEMENTATION OF UNIFORM
GRADIENT MAGNETIC FIELD
FOR GEOMECHANICS MAGNETIC
MODEL TEST
399
magnetic flux density changes greatly at both ends, Table 1. Material parameters.
i.e., the gradient is bigger corresponding to the
Figure 4. In the height of 1.4 m–1.8 m, namely the Name Value Unit
distance from the test area bottom is 0.4 m–0.8 m,
the change of magnetic flux density gradient is Young’s modulus 10e6 Pa
gentle, between 3.1 T/m–3.7 T/m, where can be Poisson’s ratio 0.3 1
taken for the best area for test. Density 1800 kg/m3
Cohesion 42e3 Pa
Angle of internal friction 17 Degree
4 NUMERICAL SIMULATION OF
GEOMECHANICS MAGNETIC MODEL
TEST OF A TYPICAL LANDSLIDE
400
401
402
W. Wang
Key Laboratory of Ministry of Education for Geomechanics and Embankment Engineering, Hohai University,
Jiangsu, Nanjing, P.R. China
Geotechnical Research Institute, Hohai University, Jiangsu, Nanjing, P.R. China
J.F. Shao
Key Laboratory of Ministry of Education for Geomechanics and Embankment Engineering, Hohai University,
Jiangsu, Nanjing, P.R. China
Geotechnical Research Institute, Hohai University, Jiangsu, Nanjing, P.R. China
Laboratory of Mechanics of Lille, UMR 8107 CNRS, Cité Scientifique, France
ABSTRACT: A discrete thermodynamic approach is presented for the modeling of coupled anisotropic
viscoplastic and damage behaviour in cohesive-frictional geomaterials. The idea is to extend the aniso-
tropic coupled elastoplastic damage model to the anisotropic one using a discrete approach. The visco-
plastic strain is induced by frictional sliding along weakness planes randomly distributed in the elastic
matrix. The evolution of induced damage is controlled by the evolution of weakness planes in connection
with the propagation of microcracks. Perzyna’s viscoplastic theory is applied to formulate the macro-
scopic viscoplastic strain. The modified Coulomb-type plastic viscoplastic potential function and damage
evolution criterion proposed by Mazars are given for each family of weakness planes. The effective elastic
modulus of damaged material is determined by the damage variable associated with each family of weak
sliding planes. The proposed model is applied to a typical geomaterial (e.g. sandstone) for the modeling
of time-dependent behaviour. The comparison between numerical prediction and experimental data is
shown. The validity of the proposed model is verified by the good fitting of strain-time curves obtained
by the multi-step triaxial creep tests for sandstones.
403
404
microcracks in some preferred orientations. For which is seen as the cumulated equivalent plastic
the description of this anisotropic damage, it is distortion in weakness sliding planes. Based on
proposed to replace the isotropic damage variable experimental data obtained on typical geomateri-
d by a damage distribution function ω ( ). In addi- als, the following particular form is proposed:
tion, it is shown that the fourth order tensor T( )
( )
is directly related to the degradation of the shear − b1γ p r
modulus (Walpole 1981). Thus, the term d K can αpr α fp − α fp α op e (13)
be generalized in the following integral form:
Physically, the hardening function defines the
1 current mobilized frictional coefficient along weak-
dK ζ ∫ ω ( n )T( n )dS
4π S ness sliding planes. α op and α fp are respectively
the initial and asymptotic values of the hardening
∫S 2 T( n )dS / 4π 2 K / 5 (8)
function. The parameter b1 controls its kinetics.
In order to define a non-associated plastic flow
It follows, from the case of isotropic damage rule, the following expression is used as plastic
distribution, i.e. ω ( ) d , that ζ = 5 /2 . The effec- potential g :
pr
tive elastic stiffness tensor C can be rewritten into
( )=σ
the form:
gp r r
t
r
n dγ
d, p ,r r
t + ( − ζ d )η p,rσ nr = 0 (14)
1
C(ω ( )) Cs μ sκ ∫ 2 ω ( )T( )dS (9) In which η is a coefficient related to plastic
pr
4π S
volumetric strain for undamaged material. In order
to describe the transition from plastic compress-
with κ 5κ /2 , and in discrete form: ibility to dilatancy observed in geomaterials; η p r is
defined as a function of plastic hardening variable
κ 15 i i γ p r and expressed as follows:
C(ω ( )) Cs μs ∑ω T
15 i =1
= 3k s J
⎛ κ 15 ⎞
μ s I − ∑ ω i Ti ⎟ : K (10)
ηp r (
η fp − η fp ηop e ) − b2γ p r
(15)
⎝ 15 i =1 ⎠
The current value of η is usually called as
pr
405
The damage evolution is physically related to The same function is proposed for the viscoplas-
propagation of microcracks and defeats in various tic potential:
orientations. Generally, the propagation is not uni-
form in space orientation. However, an isotropic
damage is adopted in this simplified version, and
gvp r ( r
t
r
n dγ
d, p,r
)=σ r
t + ( − ζ d )η p,rσ nr = 0 (23)
406
3 NUMERICAL SIMULATION
407
408
ABSTRACT: A completely new constitutive model is proposed for describing the viscoplastic behavior
of cohesive geological material—granite, which essentially possesses the morphology of granular texture.
Based on micromechanics consideration, the local inelastic deformation is attributed to the sliding along
specified direction on certain crystallographic (weakness) planes within each granular. The corresponding
macro deformation is obtained with the classical homogenization approach. The special features
of geomaterial, e.g. pressure sensitivity and plastic dilatancy, are taken into account by introducing
Mohr-Coulomb type yield criterion and non-associated plastic potential. Comparisons between model’s
predictions and test data have been carried out to check the validity of proposed model.
409
410
qα =
1
b
( − (−bλ )) α
(6) 3 NUMERICAL SIMULATION
( ( ))
3.1 Representation of the model
Rα = Q ∑ hαβ bλ β (7)
The constitutive equations can be combined
where hαβ is the interaction matrix allowing the intro- with any commercial finite element software, e.g.
duction of cross influence of β th slip system on the ABAQUS, without any difficulty. From micro-
hardening of α th slip system. If there is no cross hard- scopic point of view, each integration point is
ening or latent hardening, hαβ is an identity matrix. regarded as a polycrystal, which associates with a
large number of randomly distributed grains. The
stereographic projection of those 40 grains is as
2.4 Determination of slip rate λ α Figure 3 shows.
For time-dependent (viscoplastic) model, the slip
rate is explicitly expressed as the function of σ and 3.2 Parameters calibration
Rα as follows
The predictive capacity of the modified KBW
polycrystalline model for geomaterials is now
⎧ ⎛ 1 ⎞
⎪1 ⎜⎛ f α ⎞p ⎟ α
+ − ⎟ if f > 0
(
λ α σ ,Rα ) ⎪ ⎜⎜ α
= ⎨ η ⎝ τ c + Rα
⎜
1⎟
⎠
1
⎟⎠
⎪ ⎝
⎪0 if f α ≤ 0
⎩
(8)
411
E η τc Q
(MPa) v uf vf (s) p (MPa) hαβ b (MPa)
ACKNOWLEDGEMENT
REFERENCES
412
413
ABSTRACT: Deep mixing with lime-cement columns is a ground improvement method used to improve
the strength and deformation properties of soft cohesive soils. Due to the complex manufacturing proc-
ess, the variability in the strength and deformation properties is normally high. A rational approach to
include variability in the design process is by introducing Reliability-Based Design (RBD). This paper
presents a reliability-based design approach for Serviceability Limit State (SLS) design of soil improved
by lime-cement columns using the First-Order Reliability Method (FORM). The paper further presents
the impact of uncertainties, reliability indices and area replacement ratios on the relationship between the
characteristic value and the design value with respect to the column modulus of elasticity.
417
418
(6)
4 RESULTS
where μ E is the mean value of Ecol , α E is the
evaluated sensitivity factor, σ E is the standard
col
col col
Reliability index (β) 1.0 and 1.5 Adopted value and Eurocode 0, Honjo et al. (2010)
Maximum allowed settlement (δmax) 15 cm Adopted value
Characteristic value of column Varying Parameter is altered with different values of COV
VE ,TOT
modulus of elasticity (Ek) considering acceptable value of β
Characteristic value of clay 500 and 1500 kPa Adopted values. Normally consolidated clay
oedometer modulus ( soil )
Unit weight embankment (γemb) 20 kN/m3 TK Geo (2011). In the present study γemb is considered
and treated deterministically
Height embankment (hemb) 4m Adopted value
COVE,TOT 5–60% (altered Adopted range
in 5% intervals)
COVM,soil 30% Adopted value
Column length (L) 10 m Adopted value
Lime-cement column area 0.1, 0.2 and 0.4 Adopted values
replacement ratio (a)
419
420
421
422
J.Y. Ching
Department of Civil Engineering, National Taiwan University, Taipei, Taiwan
K.K. Phoon
Department of Civil and Environmental Engineering, National University of Singapore, Singapore
ABSTRACT: Reliability-Based Design (RBD) is known to provide a rational basis for incorporating
uncertainties in the design environment explicitly into geotechnical design. In principle, RBD should
be able to translate site investigation efforts into design savings. However, the constant partial factor
approach prohibits such translation. This study shows that the combination of two recent developments
made by the authors allows such translation. The first development is the multivariate probability distri-
bution of soil parameters. Such a multivariate distribution provides a solid basis that can translate site
investigation efforts into smaller parameter variability, e.g. smaller Coefficient Of Variation (COV). The
second development is called the Quantile Value Method (QVM). The QVM adopts constant quantile val-
ues as the design values, so the resulting partial factors scale with the parameter COV—when parameter
COV is small, partial factors are closer to one. The combination of the above two novel developments is
a cost-effective framework that is able to translate site investigation efforts into smaller parameter COVs
and in turn translate into partial factors that are closer to one.
423
424
Statistics
Characteristic
Variable Description Distribution Mean COV value
425
Table 3. Calibrated γsu for the constant PF method and calibrated η for constant QVM.
Calibrated γsu
426
427
T1 T2 T3 T4
P0 115.8
(0.35) (3.33 m)
P1 111.3 107.3 111.3 107.3
(0.32) (3.21 m) (0.27) (3.04 m) (0.32) (3.21 m) (0.27) (3.03 m)
P2 111.4 105.7 102.4 100.4
(0.31) (3.19 m) (0.25) (2.94 m) (0.25) (3.02 m) (0.21) (2.87 m)
P3 111.3 105.5 99.9 99.0
(0.31) (3.16 m) (0.24) (2.91 m) (0.24) (2.98 m) (0.20) (2.81 m)
P4 112.0 105.3 99.5 98.1
(0.31) (3.17 m) (0.24) (2.88 m) (0.23) (2.93 m) (0.19) (2.84 m)
Table 6. Final design dimensions (deviation in B to the rigorous RBD)/βA for the
constant PF method.
T1 T2 T3 T4
P0 3.61 m
(0.28 m)/3.72
P1 3.68 m 3.75 m 3.68 m 3.75 m
(0.47 m)/4.14 (0.71 m)/4.75 (0.47 m)/4.10 (0.72 m)/4.75
P2 3.68 m 3.78 m 3.84 m 3.88 m
(0.49 m)/4.16 (0.84 m)/4.75 (0.82 m)/≥4.75* (1.01 m)/≥4.75*
P3 3.68 m 3.78 m 3.88 m 3.90 m
(0.52 m)/4.16 (0.87 m)/≥4.75* (0.90 m)/≥4.75* (1.09 m)/≥4.75*
P4 3.67 m 3.78 m 3.89 m 3.92 m
(0.49 m)/4.16 (0.90 m)/≥4.75* (0.96 m)/≥4.75* (1.08 m)/≥4.75*
Table 7. Final design dimensions (deviation in B to the rigorous RBD)/βA for the
QVM.
T1 T2 T3 T4
P0 3.21 m
(−0.12 m)/2.96
P1 3.14 m 3.02 m 3.14 m 3.01 m
(−0.07 m)/3.05 (−0.02 m)/3.16 (−0.07 m)/3.04 (−0.02 m)/3.16
P2 3.10 m 2.94 m 3.02 m 2.89 m
(−0.08 m)/3.04 (−0.00 m)/3.20 (0.00 m)/3.21 (0.02 m)/3.26
P3 3.10 m 2.92 m 2.99 m 2.84 m
(−0.06 m)/3.08 (0.01 m)/3.22 (0.01 m)/3.21 (0.03 m)/3.28
P4 3.08 m 2.91 m 2.97 m 2.84 m
(−0.09 m)/3.01 (0.03 m)/3.27 (0.03 m)/3.29 (0.00 m)/3.21
3.3 Constant PF
design savings rationally using the rigorous RBD.
However, practicing engineers may not be com- Table 6 shows the required width B based on the
fortable with the rigorous RBD. The design table calibrated partial factor γsu = 0.331 given in Table 3.
in Table 3 developed by the two simplified RBD To find the required B for each information sce-
methods is more likely to be used in practice. nario, one only needs to solve Eq. (4) for B with
428
429
M. Dithinde
Department of Civil Engineering, University of Botswana, Botswana
J.V. Retief
Department of Civil Engineering, University of Stellenbosch, South Africa
431
432
433
434
the angle of friction obtained from the SPT meas- Nonetheless, the relatively lower β values depicted
urement is not directly used. Instead the key pile by NC are mainly due to the data set from D-NC
design parameters in the form of bearing capac- which has the lowest β. It appears that D-NC sig-
ity factor (Nq), earth pressure coefficient (ks) and nificantly diverges from the rest of the dataset and
pile-soil interface friction (δ) are obtained from the therefore warrants further investigations in future
derived angle of friction on the basis of empirical studies.
correlation and thus introducing some additional A comparison of all Bored (B) versus all Driven
uncertainties. (D) piles irrespective of soil type indicates very lit-
With piles in none-cohesive soils, there is a tle difference in terms of β values. This suggests
significant difference in β values between bored that pile installation method has little influence on
and driven piles. In this regard β values for bored β values. Therefore from reliability perspective, pile
piles are higher than that for driven piles. Again design should be classified on basis of soil only.
this attributed to the higher uncertainty depicted This further implies that resistance partial fac-
by driven piles relative to bored piles presented in tors and model factors should be differentiation
Table 3. Dithinde & Retief (2013) further explains on basis of soil properties. Further differentiation
that this implies that the densification of the soil into pile installation method is just a refinement
surrounding the pile emanating from the pile driv- and elaboration.
ing process is not well captured in the selection of Generally the β values for all pile classes are
the soil design parameters and hence current prac- above the target β of 3.0 for the reference reliability
tice is conservative in selecting design parameters class R2 for which reliability procedures are speci-
for driven piles. fied in SANS 10160. When taking into account the
Comparing β values for all piles in Cohesive redundancy due to group and system effects, the
materials (C) versus all piles in Non-Cohesive (NC), β values will increase appreciably. Therefore for
shows distinctively different level of reliability sug- pile groups (which is the common practice), SANS
gesting the two broad classes based on soil type 10160-5 yields β values that are significantly higher
should be treated separately in calibration studies. than the target β of 3.0 for class R2. Conversely,
435
436
J.P. Li
Department of Geotechnical Engineering, Tongji University, Shanghai, China
J. Zhang
Department of Geotechnical Engineering, Tongji University, Shanghai, China
ABSTRACT: As part of the worldwide efforts for implementing reliability-based design in geotechnical
engineering, Shanghai also calibrated the resistance factors for design of shallow foundations using
reliability theory. This paper introduces how the resistance factors are determined in Shanghai. A typical
shallow foundation is assumed to be constructed on 142 sites in Shanghai. The mean values of cohesion
and friction angle are determined based on the site-specific data. The Coefficients Of Variation (COV)
of the cohesion and friction angle are determined based on regional experience considering variance
reduction due to spatial correlation of soil properties. The first order reliability method is then used to
determine the partial factors for each site. It is found that if the notional Factor Of Safety (FOS) is 2.5,
the reliability indexes of the 142 foundations are in the range of 2.3–4.2 with a mean of 3.35. If the target
reliability index is 3.35, the partial factor for cohesion for the 142 sites have a mean of 2.76 and a standard
deviation of 0.32, and the partial factor for friction angle has a mean of 1.17 and a standard deviation
of 0.07. It is then recommended that the partial factors for cohesion and friction angle be 2.7 and 1.2,
respectively. The adopted partial factors produce design bearing capacity predictions consistent with local
experience.
437
Substituting Eq. (10) into Eq. (9) yields Let c* denote the value of c at the design point.
The partial factor for c can be determined using
Ru (cc φ ) − SL SD = 0 (11) the following relationship
438
Substituting Eq. (18) into Eqs. (16) and (17), the 5 CALIBRATION RESULTS
mean values of live and dead loads can be calcu-
lated as follows 5.1 Current level of reliability and target
reliability index
Rnu Figure 1 plots the calculated reliability indexes and
SnD = (19)
( + )K mean values of cohesion for the 142 sites. The reli-
ability indexes are in the range of 2.3–4.2 with a
ρRnu mean of 3.35. For comparison, Phoon et al. (2003)
SnL = (20)
( + )K reported that the reliability index of existing drilled
shafts under compression to support transmission
line structures in North America is in the range
Given the values of cn, ϕn, ρ and K, Eqs. (19)
of 2.4–3.6. Figure 2 plots the calculated reliability
and (20) can be used to determine the values of SnD
indexes and mean values of friction angle for the
and SnL based on which the reliability index cor-
142 sites. Comparing Figures 1 and 2, there are less
responding to a specific K can be found out using
scatter in the relationship between reliability index
Eq. (14) as the limit state function. If one knows
and the friction angle, probably because the reli-
the target reliability index, one can adjust the val-
ability indexes of the foundation are more affected
ues of SnD and SnL in Eq. (14) subjected to the con-
by the friction angle.
straints set by Eqs. (19) and (20) until the target
The target reliability index is often determined
reliability index is satisfied, which can then be used
with reference to the reliability index of exist-
to determine the partial factors corresponding to
ing buildings (e.g., Allen 1975; Phoon 2003) or
the target reliability index.
To calibrate the resistance factors for design,
consider a strip foundation with B = 3 m, D = 1 m.
In Shanghai, the typical design FOS is K = 2.5. As
the soil data may vary with site, the foundation is
assumed to rest on 142 different sites in Shanghai
to obtain a set of robust partial factors that can
be generally used in Shanghai. At a site, while the
mean values of soil properties are relatively easier
to determine, the Coefficients Of Variation (COV)
are much hard to estimate based on limited site-
specific data. Thus, the mean values of the cohe-
sion and friction angle take the site-specific values,
and the COV values of the cohesion and friction
angle take the regional values. In Shanghai, the Figure 1. Relationship between cohesion and reliability
COV of cohesion of the clay at shallow depths index for the 142 sites.
439
440
ϕ (°) MB MD Mc
0 0 1 3.14
2 0.03 1.12 3.32
4 0.06 1.25 3.51
6 0.1 1.39 3.71
8 0.14 1.55 3.93
10 0.18 1.73 4.17
12 0.23 1.94 4.42
14 0.29 2.17 4.69
Figure 5. Relationship between γc and μφ for the 142 sites. 16 0.36 2.43 5
18 0.43 2.72 5.31
20 0.51 3.06 5.66
22 0.61 3.44 6.04
24 0.8 3.87 6.45
26 1.1 4.37 6.9
28 1.4 4.93 7.4
30 1.9 5.59 7.95
32 2.6 6.35 8.55
34 3.4 7.21 9.22
36 4.2 8.25 9.97
38 5 9.44 10.8
40 5.8 10.84 11.73
β γc γϕ
441
442
S.G. Paikowsky
University of Massachusetts Lowell (UML), Lowell, Massachusetts, USA
GeoDynamica, Inc., Newton, Massachusetts, USA
C.M. Marchionda
Nobis Engineering, Inc., Lowell, Massachusetts, USA
Previously UML, MA, USA
S. Amatya
Golder Associates Pty Ltd., Richmond, Victoria, Australia
Previously UML, MA, USA
M.C. Canniff
GeoDynamica, Inc., Newton, Massachusetts, USA
A.S. Budge
Minnesota State University Mankato, Mankato, Minnesota, USA
ABSTRACT: Driven piles are the most common foundation solution used in bridge construction
(Paikowsky et al., 2004). Their safe use requires to reliable verification of their capacity and integrity.
Dynamic analyses of driven piles are methods attempting to obtain the static capacity of a pile, utilizing its
behavior during driving. Dynamic equations (aka pile driving formulas) are the earliest and simplest forms
of dynamic analyses. The development and the examination of such equation tailored for given demands are
presented. After establishing the bridge pile construction practices of MnDOT, a database of driven pile case
histories relevant to these practices was built. The databases were utilized to investigate previous MnDOT
(and other) dynamic formulas and use object oriented programming for linear regression to develop a new
formula that was then calibrated for LRFD methodology and evaluated for its performance. A standalone
control database and a database of dynamic measurements enabled independent evaluation of the for-
mula leading to the final stage in which the equation was adjusted for field application and past experience.
Though developed for MnDOT, its use is universal for driven piles within the provided limitations.
443
444
2.3 MnDOT/LT-Pipe-Piles
Figures 3 and 4 present graphically some of the is presented along with the number of piles it is
important features of MnDOT/LT 2008 Pipe Piles based upon. For example, 12 and 65 case histories
database along with a comparison of data reflect- of the database are related to 12.00 and 12.75 inch
ing the MnDOT foundation practices as presented diameter piles, respectively. The mean failure load
earlier in section 1. The information in Figure 3 of these cases was 388 +/− 173 kips and 372 +/−
provides mean failure load (+/− 1 standard devia- 195 kips (1 SD) for the 12 and 12.75 inch diam-
tion) for each pile type/size category (by pipe pile eter piles, respectively. The mean factored load of
diameter) comprising the database, along with the 12 and 12.75 inch piles by the MnDOT is 155
the number of cases related to that information. and 120 kips, based on 1055 and 116 piles, 12.0
In addition, the mean LRFD factored (design) and 12.75 inch diameter pipe piles, respectively.
load for the MnDOT for the applicable pile cases These data suggests that the mean safety margin
of the MnDOT 12 and 12.75 inch diameter piles
is 2.503 +/− 1.116 and 3.100 +/− 1.625 in compari-
son with the database information (not including
the load factor) or the covering of approximately
1.3 S.D range (lower value of resistance is 1.29 and
1.35 standard deviations from the mean, hence
approximately taken as 1.3) translates to a target
reliability of β = 1.3 and a probability of failure
pf = 9.8%. The information in Figure 4 presents the
distribution of the case histories in the database
based on the pile sizes in comparison with the dis-
tribution of use of the same pile by MnDOT. To
be relevant, the frequency of use of the pipe piles
by the MnDOT, presented in Figure 4, reflect-
ing the use of the particular pile type out of the
pipe piles only and not out of all driven piles. For
example, 77 pile cases of 12 and 12.75 inch diam-
eter piles are available in the database (46.1% of
all cases), while the MnDOT uses these diameter
piles in 50.0% of the projects where pipe piles are
being used. The major difficulty of the database
Figure 3. Range of pile capacity based on static load as presented in Figure 4 is evidently related to the
test (mean +/− 1 S.D.) and MnDOT mean factored design 16 inch diameter closed ended pipe piles. This type
loads sorted by pipe pile type and pipe pile diameter. of piles are used by the MnDOT in 22.9% of all
445
446
Table 2. Statistical parameters and LRFD calibration for the developed equation (2) for H piles and pipe piles, EOD
condition only.
Resistance factor φ
Mean bias β = 2.33, pf = 1%, φ/λ
# of measured/ Stand. Coeff. Best fit line Coeff. of redundant Efficiency
Pile cases calculated dev. of var. equation determination factor
type (n) (mλ) (σλ) (COVλ) (least square) (r2) FOSM MC3 Recom. (%)
H 125 1.0163 0.3599 0.3542 Ru = 0.880 * Rs 0.896 0.495 0.542 0.55 54.1
Pipe 99 1.1089 0.5955 0.5370 Ru = 0.805 * Rs 0.849 0.364 0.385 0.35 31.6
447
Table 3. Dynamic equation predictions for H piles and pipe piles EOD condition only.
Notes: 1Searched coefficient for the equation Ru = Coeff. Eh log(10 N ), 2MnDOT energy range contains
hammers with rated energies between 42.4 and 75.4 k-ft.
analysis under the various sub-categorizations to good performance of the proposed format and
the EOD cases only. The obtained results summa- obtained coefficients.
rized in Table 3 suggest the following: c. When restricting the EOD data to diesel ham-
mers only and a blow count of equal or greater
a. For both pile types under all EOD data selec- to 4 BPI (with or without the outliers) the rec-
tion criteria (with or without the outliers), the ommended coefficients are 32.1 to 33.5 for the
recommended coefficient varied between 34.2 H piles and 30.0 to 30.5 for the pipe piles. Both
to 36.7 reaffirming the coefficient of 35 recom- subsets contain significant number of cases
mended for the general equation as appeared in (38 H piles and 38 pipe piles when eliminating
equation 2. outliers).
b. All regressions resulted with a coefficient of d. When further restricting the conditions
determination greater than 0.85 suggesting described in (c) above by looking at the energy
448
factor (%)
Notes: Calculated capacity using MnDOT new dynamic equation Ru = Coeff. Eh log(10 N ), Rs is the static capacity of the pile examined by Davisson’s failure criterion,
Efficiency
MnDOT practice, the subsets decrease to 13/12
H pile cases and 16/14 pipe pile cases, with and
54.1
50.6
55.8
62.2
48.2
57.0
φ/λ
without outliers, respectively. These are mar-
ginal size sets that result with coefficients vary-
Recom
β = 2.33, pf = 1%, redundant
ing between 31.2 to 34.4 for H piles and 33.1 to
0.55
0.60
0.60
0.65
0.55
0.60
33.3 for pipe piles.
e. Close examination of the most restrictive sub-
Resistance factor φ
sets described in (d) above (i.e. 13 H piles and
0.542
0.632
0.628
0.674
0.553
0.620
MC3
16 pipe piles before removing the outliers) show
that in both subsets a relatively (to the subset
size) large group of cases are of different piles
FOSM
of the same size tested at the same site (e.g. 7
0.495
0.578
0.566
0.598
0.508
0.558
out of the 16 pipe piles are 14” diameter piles
from Deer Island project in Massachusetts and
6 of the H piles are 12 × 53 from site no. 37
Coeff. of
determ.
in Canada). As such, the data are too biased
0.896
0.896
0.907
0.935
0.870
0.924
as not only the set is marginal in size, but
(r2)
about 50% of the cases are related to the same
project. The statistics and coefficient obtained Statistical parameters and resistance factors of the new detailed MnDOT dynamic equation for H piles.
Ru = 0.888 Rs
Ru = 0.754 Rs
Ru = 0.812 Rs
Ru = 0.873 Rs
Ru = 0.759 Rs
Ru = 0.889 Rs
from that subset should, therefore, cautiously
(least sq.)
be applied.
line eq.
Best fit
f. As a lower coefficient means a more conserva-
tive capacity evaluation, the above discussion
and the observations presented in (c) should
serve as the guideline for the new MnDOT
(COVλ)
of var.
0.3542
0.3542
0.3176
0.2762
0.3968
0.3135
dynamic equation that suits better to MnDOT
Coeff.
0.4531
0.3297
Stand.
hammers) is therefore:
Ru 30 Eh * log(10 * N ) (3)
Mean bias
measured/
calculated
1.0163
1.1856
1.0760
1.0458
1.1419
1.0518
(mλ)
EOD, diesel,
EOD, diesel,
dynamic equation
EOD only
EOD only
Condition
30
30
125
39
38
13
12
(n)
2
3
449
Notes: Calculated capacity using MnDOT new dynamic equation Ru = Coeff. Eh log(10 N ), Rs is the static capacity of the pile examined by Davisson’s failure criterion,
factor (%)
Efficiency
evaluation compared to that of equation (2) that
was developed for that situation specifically. The
31.6
34.8
42.0
49.6
45.3
81.3
φ/λ
greater mean bias obtained when using equation
(3) allows, therefore, to select a consistent resist-
ance factor of φ = 0.60 to be used for all the cases
Recom
β = 2.33, pf = 1%, redundant
0.35
0.45
0.40
0.45
0.85
0.90
when applying equation (3). This conclusion was
further examined and reaffirmed against an inde-
Resistance factor φ
0.905
1.012
MC3
0.766
0.825
are presented in Table 5. This is done to examine
the applicability of using the detailed equation (3)
under all driving conditions.
Coeff. of
determ.
0.974
0.988
Ru = 0.878 Rs
Ru = 0.895 Rs
Ru = 0.805 Rs
0.5955
0.5370
0.4284
0.3472
0.1818
0.1151
0.2051
0.1273
Stand.
1.1284
1.1065
(mλ)
4.1 Overview
EOD, diesel, BC ≥ 4 BPI
EOD, diesel,
30
30
(n)
99
41
38
16
14
2
3
450
451
5 EQUATION ADJUSTMENT
AND FINAL FORMULATION
5.1 Overview
The equations presented above were found to
best predict pile capacity when examined with
separate data not related to its development.
Several additional developments were required in
order to finalize the equation format to be used
by the MnDOT. These developments include the
Figure 10. CAPWAP vs. new MnDOT formula following:
(coefficient 35) all CIP piles.
Statistics
Figure No. of Coeff. of
Relations no Category cases n Mean S.D. COV Best fit line determ. (r2)
Notes: Ru is the calculated capacity using each of the dynamic formulae. Rc is the Capacity determined by
CAPWAP.
452
453
454
ABSTRACT: By designing mainly axially loaded offshore foundation piles in North Sea conditions,
two commonly used design methods, the β-method and the ICP-method, were used. Thereby the required
minimum embedded pile length may vary a lot. The actual safety which is achieved with a certain pile
length and the corresponding prescribed partial safety factors is not known. For general assumed pile and
soil conditions it could be determined that a reliability index between β = 3.5 and β = 4.1 is established.
Further it could be shown that the prescribed partial safety factors should be decreased for the β-method,
where an increase for the ICP-method should be considered. By performing reliability based design also a
more accurate determination of the required pile length for a specified safety is possible.
One of the main reasons for the increase of The tension bearing capacity of mainly axially
renewable energy in Germany over the last years loaded piles consists basically of the mobilized
is the continued extension of wind energy. For friction between the pile outer shaft area and
the next years additionally several offshore wind the surrounding soil. Additionally if an open
farms are planned to be built in the North Sea. ended pile is used two different condition states
Since many projects at the German North Sea sites “plugged” or “unplugged” have to be considered.
are going to be realized outside of the 12 sea mile In the unplugged case also the friction resistance
border, water depths greater than 30 m are often between the inner pile shaft area and the inner soil
faced. For such water depths jacket and tripod can be taken into account. By assuming a plugged
supporting structures with mainly axially loaded condition the effective weight of the inner soil plug
foundation piles are mostly used. can be added to the bearing capacity.
Thereby the pile resistance can be calculated by Generally it can be observed that in almost all
applying the well known β-method (also known as design cases the tension limit state is the control-
API-method) recommended in the guideline of the ling one with regard to the required pile length.
American Petroleum Institute (API 2007). In the Therefore only tensile capacity is considered
course of publishing errata and supplements for here. The pile resistance for the unplugged or
the existing guideline, new simplified CPT-based plugged condition can be computed by applying
design methods were introduced following the Equation 1:
recommendation of the corresponding research
groups, like the ICP-method introduced by Jardine Rt Ao ∫ ft ( z ) z Ai ∫ ft ( z ) dz G′
G p] (1)
et al. (2005).
For many design cases the ICP-method leads to
considerably differing resistances compared with where Ao = outer pile shaft area; ft(z) = skin fric-
the API-method. This affects also the deterministic tion for tension loading; Ai = inner pile shaft area;
design, where this deviation leads to a strong vari- and G′p = effective weight of the inner soil plug.
ation of the required embedded pile length. For Following the API-method the skin friction of
practical applications it is of interest which method offshore foundation piles in sand for tension con-
is more reliable for determining the pile length and ditions can be computed by Equation 2. Thereby
which method should be used, respectively? In the friction is basically determined by multiplying
addition it is also unknown which safety level— the effective vertical stress with a β-value. Further
measured by the reliability index β—is finally a limitation of the skin friction is included. Both,
achieved by applying one of these methods with the shaft friction factor as also the limiting fric-
the corresponding partial safety factors in a deter- tion value, only depend on the relative density of
ministic design. the soil and may vary with depth. The additional
455
ft ( z ) = 2 3 β ′v ≤ 2 3 ft ,max (2)
456
Vk Lγ R Rk (5)
457
ICP-method, all pile tests 1.01 0.28 0.28 Unit weight γ′ in kN/m3 10 1 N
API-method, all pile tests 1.15 0.58 0.50 Inherent wqc in MN/m2 0 6 N
ICP-method, only related tests 1.16 0.19 0.16 variability qc
API-method, only related tests 1.26 0.14 0.11 Measurement eqc 1 0.15 N
error qc
Transformation eϕ in ° 0 2.8 N
error ϕ
Jardine et al. (2005) obtained a model error for Model mAPI 1.26 0.14 N
both methods where the following pile tests were error API
considered: Closed-ended and open-ended piles, Model mICP 1.16 0.19 N
tension and compression test, steel and concrete error ICP
piles and loose to very dense soil state. Pile diameter D in m 2 – D
Altogether 81 pile tests were evaluated. The sta- Pile wall t in m 0.05 – D
tistical values are shown in Table 2. Therein the thickness
model error is defined as the ratio of the measured Pile length L in m 20–60 – D
to the calculated resistance. Axial loading V in MN 5.4 1.89 G
Achmus & Müller (2010) compared measured *N = Normal distributed; D = deterministic value; and
with calculated results for the API and the ICP- G = Gumbel distributed.
method only with pile tests which are closely
related to the boundary conditions in the North
Sea (see section 3). By this it can be seen how relia- The vertical autocorrelation length for the cone
ble these methods are to a special boundary sector. tip resistance is indicated to be between 0.1 m and
In detail the following pile tests were considered: 2.2 m. Hence, the value of θ = 0.6 m was chosen for
Open-ended piles, tension test, steel piles, dense to all simulations.
very dense soil state and slenderness ratios between The buoyant unit weight γ ′, the measurement
10 and 40. error for the cone tip resistance eqc and the transfor-
By taking into account only pile tests within mation error for the internal friction angle eϕ were
a range of 95% confidence interval, only 6 pile modeled as uncorrelated and normally distributed
test for the ICP-method and only 4 pile tests for random fields with typical values for the mean and
the API-method could be found. The mean and standard deviation. The stochastic formulation for
standard deviations for these related tests are sum- the load was chosen according the recommenda-
marized in Table 2 and should be seen as rough tion by Holicky et al. (2007) for a 50-year extreme
approximations of the real values. environmental event.
As it can be seen the overall model error for All applied variables with their corresponding
the ICP-method is smaller in mean and standard stochastic moments are summarized in Table 3.
deviation compared to the overall model error for
the API-method. By comparing the model error
6 OBSERVED RESULTS
only for the related tests, it can be said that in aver-
age both methods underpredict the resistance. The
6.1 Impact of uncertainties
standard deviation seems also to be closer related
to each other as in the overall case. For the computed data sets first the COVs of the
resistance distributions are evaluated. Addition-
ally the corresponding parts which arise from
5.3 Stochastic variables and simulation
the model error and soil uncertainties were esti-
For each performed simulation with a certain pile mated, as depicted in Figure 4 for very dense soil
length 6 million realizations within a plain Monte- condition.
Carlo Simulation (MCS) were computed. For both methods the total variation of the
The inherent variability of the cone tip resist- resistance is much lower than of the load. For the
ance wqc was assumed to be constant with depth. API-method almost 11% and for the ICP-method
By applying a standard deviation of 6 MN/m2 about 17% were obtained. Further the model error
the COV is varying between 0.75–0.20 for dense is almost exclusively responsible for the variation of
and 0.40–0.10 for very dense soil conditions. the resistance, especially for the API-method. Also
According to Phoon & Kulhawy (1999) the COVs it can be noticed that the deviation arising from
were thereby in a typical range of 0.81 and 0.10. soil uncertainties is decreasing with embedded pile
458
6.2 Safety
The safety of a system can be described in terms of
a reliability index β corresponding to Equation 7.
In the following, β is also termed safety.
Figure 5. Evaluated safety with increasing pile length.
In the Eurocode 0 three reliability classes with
the corresponding safeties were proposed. Further
it is remarked that the partial safety factors are
construed for the second reliability class and there- Table 4. Required deterministic pile lengths with
fore a safety of β = 3.8 should be established in a corresponding reliability index.
design system. However, offshore pile foundations
may also be associated to the first reliability class Dr = 0.75 Dr = 0.93
which demands a safety of β = 3.3. Design method API ICP API ICP
β = Θ −1
−
( − ) (7) Deterministic length 51.85 45.91 45.19 30.24
Reliability index 4.1 3.5 4.0 3.5
459
Dr = 0.75 Dr = 0.93
Design
method API ICP API ICP
γL γR
β = 3.3 1.59 (1.52)* 1.86 1.55 2.60 (1.87)*
β = 3.8 1.96 (1.82)* 2.29 1.84 3.12 (2.31)*
460
461
A.A. Henriques
University of Porto—LABEST, Porto, Portugal
ABSTRACT: This paper aims to show reliability concepts for axial pile foundations’ design and safety
aspects, revealing that reliability analysis do not imply a significant increased effort when compared to
the traditional design methodologies. In fact, calculations using reliability concepts provide very useful
information, such as the uncertainty that mostly influence the probability of failure. Therefore, the most
used reliability methods, FORM and MCS, are here presented and used to evaluate the safety of the bear-
ing capacity of an axial bored pile considering the geotechnical uncertainties and different load values.
For the application example presented, the physical uncertainties of actions, the inherent soil variability
and model error (pile’s resistance) were determined by experimental in situ tests (SPT) or collected from
bibliography.
463
Zero I II III IV
Geotechnical parameters
Calculation method (deterministic)
Design parameter (statistical basis)
Variability of parameters:
– mean and SD
– PDF
Costs/consequences
Type of analysis: Global SF Partial SF e.g. FORM e.g. MCS Risk
parameters—mean and Standard Deviation in calculations (Phoon & Honjo 2005, Phoon
(SD) or coefficient of variation (COV = SD/ 2008, Zhang et al. 2011).
mean); A full account of the reliability methods devel-
– RA level III: full probabilistic analyses, it takes opment and evolution can be found in Manohar &
into account all the probabilistic characteristics Gupta (2005) and Baecher & Christian (2003).
of the RV—mean, SD and Probability Distribu-
tion (PDF) and when the problem is complex
simulations methods are used; 2.3 RA-based safety evaluation
– RA level IV: risk analyses, where the conse- In order to carry the reliability-based safety evalu-
quences of failure are also taken into account and ation, the following procedure is used for both
the risk is used as a measure of the reliability. FORM- and MCS-based RA:
– definition of the significant failure modes and
formulation of their functions → g(Xi):
2.2 RA methods
One of the first attempts to solve this kind of M R−E g(X i ) (1)
problems involving probabilities and statistics,
with non-normal variables and complex integral
where M is the safety margin, R denotes the
calculations, was the First Order Second Moment
resistance, E denotes the action, g is the perform-
(FOSM) proposed by Cornell (Cornell 1969) that
ance function, and Xi are the random variables;
first introduced the concept of reliability index.
– identification of the random and deterministic
FOSM ignores the shape of the PDF, using only
variables → Xi;
the mean and SD of the RV, and the calculation
– description and characterisation of the RV,
model is linearized using Taylor’s expansion.
namely the statistical parameters—mean, SD,
Next, other authors (e.g. Hasofer & Lind 1974)
COV and PDF type—as well as identification of
continued to develop the method, and FORM
the dependencies among them (if existent);
was proposed, solving the invariance problem
– selection of the target reliability index → βT and
of FOSM and being nowadays one of the basic
correspondent probability of failure → pf.
tools for RA level II. These are called first and
second order methods and are approximation The target depends on many factors, such as
techniques. the type of structure, the social tolerance for non-
For more accurate calculations (RA level III) compliance, among others. Its selection can be
the simulation methods are used, these methods based on previous similar construction projects
can be applied to RV with non-normal distribu- that met predefined requirements or on recom-
tion and complex performance functions (e.g. mendations in design codes—βT for ultimate limit
requiring nonlinear functions or finite element state design should be 3.3 to 4.3 (CEN 2002a)
methods). RA level III use all the statistical infor- or 1.3 to 4.3 (ISO 2394, 1998). These values cor-
mation about the RV and the most known and respond to pf of approximately 10−1 and 10−5, as
method of reference is the MCS. To MCS it is shown in Figure 1.
possible to add reduction of variance techniques, For safety evaluation of the bearing capacity,
such as Stratified sampling, Latin hypercube or the problem main characteristics are fixed, such
Importance sampling, with the goal of reducing as dimensions and materials, and probabilities of
calculation time, by selecting only points near fail- failure (reliabilities) are calculated for different
ure, or selecting specific groups to be represented load values.
464
465
466
3.5 Results
Recall that the safety evaluation approach assumes
a fixed length and diameter of the pile while it cal-
culates the probabilities for different load values.
Safety evaluation approach was carried out for the
following load values E = Gk + Qk = [400; 600; 700;
750; 800; 850; 900; 1000; 1200; 1400] kN. FORM
Model
– tip 1.12 0.706 Lognormal Figure 5. Reliability-based safety evaluation of vertical
– side 1.07 0.492 Lognormal bearing capacity of application example.
Soil
– NSPT,tip 10.26.1.91 z 4.6* Normal
– NSPT,side 4.6** Normal and MCS results are presented in Figure 5. MCS
Actions†
using n1 = 200,000 and n2 = 1,000,000 are presented
– permanent 1.0 0.10 Normal
to see that the number of simulations stabilizes
– variable 0.6 0.21 Gumbel
from n > 200,000. Figure 5 also depicts a light line
*Will be reduced taking into account the influence zone on marking the hypothetical/considered design load
the pile tip (3 × diameter) as averaging over the thickness. (total of 800 kN) and the load test result (1350 kN).
**Will be reduced taking into account the length of the Moreover, Figure 6 presents the obtained values
pile as averaging over the thickness. †Considered Gk = Qk. for α factors (FORM).
467
REFERENCES
Figure 6. FORM sensitivity factors for the application
example.
Baecher, G.B., & Christian, J.T. 2003. Reliability and
statistics in geotechnical engineering. Chichester, UK:
John Wiley & Sons.
CEN. 2002a. European Committee for Standardization.
It is possible to conclude that:
Eurocode 0: Basis of Structural Design. EN 1990.
– for MCS with n1 and n2 the results are CEN. 2002b. European Committee for Standardization.
consistent; Eurocode 1: Actions on Structures. EN 1991.
– FORM results give very acceptable approxima- CEN. 2007. European Committee for Standardization.
Eurocode 7: Geotechnical Design. EN 1997.
tions to MCS results;
Cornell, A.C. 1969. A probability-based structural code.
– the results also show a clear relationship (expo- Journal American Concrete Institute 66(12): 974–985.
nential type) between the reliability and the Hasofer, A.M. & Lind, N.C. 1974. Exact and invariant
load; second-moment code format. Journal of the Engineer-
– if one considers the hypothetical action (800 kN) ing Mechanics Division, ASCE 100(EM1): 111–121.
it is possible that it does not meet the require- Henriques, A.A., Calheiros, F. & Figueiras, J.A. 1999.
ments, that usually recommend a reliability Probabilistic modelling of nonlinear behaviour of
index between 2.5 and 4.0; concrete structures. In Schueller, & Kafka (Ed.), Proc.
– and it is also possible to see that the uncertainty ESREL, ‘99 (pp. 495–500) Munich, Germany: A. A.
Balkema.
that has more influence in the reliability is the
Holicky, M., Markova, J. & Gulvanessian, H. 2007. Code
model error, being the other uncertainties con- calibration allowing for reliability differentiation and
siderably less important. production quality. Proc.10th International Confer-
ence on Applications of Statistics in Civil Engineering,
Tokyo, Japan. Taylor & Francis.
Honjo, Y. & Setiawan, B. 2007. On conditional estima-
4 CONCLUSIONS tion accuracy of spatial average of soil proprieties
and sample size. Proc.10th International Conference on
This paper presented a simple methodology for a Applications of Statistics in Civil Engineering, Tokyo,
reliability-based safety evaluation of an axial pile Japan. Taylor & Francis.
foundation, based on FORM and MCS. These are Honjo, Y., Suzuki, M., Shirato, M. & Fukui, J. 2002.
believed to be easy methodologies and following Determination of partial factors for a vertically loaded
the trend of reliability-based assessments. They pile based on reliability analysis. Soils and Foundations
can help support the design of pile foundations 42(5): 91–109.
Honjo, Y., Jliati, M.N. & Ishino, J. 2007. Effects of spa-
and also try to eliminate the possible confusions
tial variability and statistical estimation error in pre-
and difficulties that traditional reliability method- diction of settlement of shallow foundation. Proc.10th
ologies used in structures can cause to geotechni- International Conference on Applications of Statistics
cal designers in practice. This paper also presented in Civil Engineering, Tokyo, Japan. Taylor & Francis.
an application example, for which calculations Honjo, Y., Hara, T. & Kieu Le, T.C. 2010. Level III
were repeated for different loads. One of the main reliability based design of examples set by ETC10.
results is the relationship between probability of Proc. 2nd International Workshop on the Evaluation of
failure and the load value and the high influence Eurocode 7, Pavia, Italy.
of the model error in these results. Furthermore, ISO. 1998. International Organization for
Standardization. ISO 2394: General principles on reli-
the actions uncertainties were demonstrated to
ability for structures.
have not an important role in the reliability. This Jones, A.L., Kramer, S.L. & Arduino, P. 2002.
type of reliability-based assessments allows a more Estimation of uncertainty in geotechnical properties
rational way to deal with uncertainties of a prob- for performance-based earthquake engineering. Pacific
lem, instead of just introducing safety factors into Earthquake Engineering Research Center, University
the calculation. of Washington, USA.
468
469
C.Y. Ou
Department of Construction Engineering, National Taiwan University of Science and Technology, Taipei, Taiwan
ABSTRACT: Cross walls are commonly used for constraining the excavation-induced deformations as
well as the subsequent building damages in Taiwan. Based on the real case histories, a regression model
that is capable of predicting maximum wall displacements for excavations with or without cross walls is
proposed. With the consideration of uncertainties in the proposed model, a simplified Reliability-Based
Design (RBD) method in the form of partial factors for the Serviceability Limit State (SLS) design of
excavations is developed. The partial factors are calibrated by rigorous reliability theory and are verified
to produce design outcomes that meet the target reliability index. The use of cross walls which is extremely
effective in achieving the target reliability level is highlighted.
471
+ a ⎡⎣ l ( H )⎤⎦ + a ⎡ ln (L )⎤ + a ⎡⎣ ln (K ) ⎤⎦
2 2 2
′
where L′, B, and He (in meters) are respectively 6 e ⎣ ⎦ 7 8
472
473
δ m,d ( mm )
*
⎛ a0 + a1B + a2 l ( H e ) ⎞
⎜
( )
l ⎡ γ su × su σ v ⎤ + a4 ln L ⎟
⎜ + a3ln ⎣
′
m⎦
′
( )
⎟
⎜ ⎟
⎜ + 5l ( ) + a6 ⎡⎣ l ( )⎤⎦ ⎟
2
= exp ⎜ 2⎟ × γ
( )
2
⎜ + a7 ⎡⎣ l L ⎤⎦ + a8 ⎡⎣ l ( )⎤⎦ ⎟ δ
′
⎜ ⎟
⎜ ⎣ ( m
)
l ⎡ γ su × su σ v′ ⎤ ln L′
⎜ + 9 ln
⎦ ( ) ⎟
⎟
⎜⎝ + 10 ln ( )
l L ln (
′
) ⎟⎠
(11)
(
γ δ = exp ε −η ) = exp ⎡⎣σ ε × Φ 1 (1 − η)⎤⎦ (12)
γ su = exp ⎡(e1 )η ⎤ = p ⎡⎣σ1 × Φ −1 ( η) ⎤⎦
⎣ ⎦
474
475
476
A. Puissant
Laboratoire Image, Ville, Environnement, CNRS and University of Strasbourg, Strasbourg, France
J.-P. Malet
Laboratoire Image, Ville, Environnement, CNRS and University of Strasbourg, Strasbourg, France
Institut de Physique du Globe de Strasbourg, CNRS and University of Strasbourg, Strasbourg, France
ABSTRACT: In this study, landslide risk to buildings is assessed through an index-based, semi-
quantitative approach by integrating information obtained from multiple sources: landslide inventory and
susceptibility maps; information concerning the exposed buildings and their vulnerability. Once the risk
value of the individual buildings have been determined, a methodology based on the geostatistical tech-
nique of Indicator Kriging (IK) is proposed in order to cluster the exposed buildings into zones of equal
risk. Each risk zone consists of the areas where the probability of exceedance of the risk zone value—that
is produced through geostatistical analysis—is greater than a given threshold. The final risk zone map is
produced by the synthesis of each risk zone. The proposed methodology is applied to the Barcelonnette
Basin, where areas of low, moderate and high landslide risk zones are determined. The risk zone map
produced objectively can provide useful guidelines for risk management design.
479
480
4 IMPLEMENTATION TO THE
3.2 Indicator Kriging for locating risk zones BARCELONNETTE AREA
Indicator Kriging can be used to cluster the
exposed elements in groups of equal risk. The 4.1 The study area
whole procedure that can be totally implemented The test site (Fig. 1) is located on the north-facing
in a GIS environment is summarized in the fol- hillslope of Barcelonnette Basin extending over an
lowing steps. area of about 100 km2 and it is representative of
Once the risk value, either quantitative or quali- the climatic, lithological, geomorphological and
tative, of the buildings has been determined, the landcover conditions common to several regions
risk R(x) at location x (x∈D, D the study area) of the South French Alps. Barcelonnette Basin is
can be considered as a regionalized variable whose a highly landslide prone area characterized by a
value is equal to the building’s risk value inside the large variety of slope movements (e.g., Maquaire
building’s area and zero otherwise. et al. 2003, Malet et al. 2005, Thiery et al. 2007),
The number of the desired risk zones nk form-
ing the risk zone map has to be selected, along
with the risk value Ro,k of each zone. For the kth
risk zone, the indicator function I(x; Ro,k) and the
variogram γk that corresponds to its spatial varia-
tion is determined. Then, by applying the kriging
equations at different locations xo of the field, an
estimate of the indicator function I(xo; Ro,k) at the
whole study area can be calculated. At this point
it should be noted that in order to provide moni-
tonicity, a single variogram—that of the median
of thresholds—can be used for the whole analysis
(Cressie 1993) instead of estimating different var-
iograms for each risk zone. Also, the kriging esti-
mations equal the input data values. To produce
more smooth risk zones that can include exposed
buildings of different risk value this property of
kriging should be relaxed. Thus, a nugget effect
ne, i.e. an offset at the origin, fully attributed
to measurement errors should be added to the
variogram.
According to the geostatistical theory, the indi-
cator function I(xo; Ro,k) is an estimate of the prob-
ability of the R(x) not to exceed the threshold Ro at
that point. However, the risk zones should contain
all the areas where R(x) is equal or greater than Figure 1. (a) Typical landscape of the north-facing
Ro, so, it is the complementary event that is actu- hillslope of the Barcelonnette Basin and (b) shaded relief
ally of interest. Thus, a map with the probability map and distribution of landslides (Thiery et al. 2007).
481
482
To validate the selected parameters, a 10% of the occur. So, an emergency plan should include extra
input data was randomly selected and excluded caution for this area in case of a landslide event.
from the analysis. The selected exponential vari-
ogram was used with the remaining data to predict
the indicator function at the excluded locations. 5 CONCLUSIONS
The very small value of the mean absolute error
(0.03) verifies the efficiency of the adopted vari- The present paper provides an index-based,
ogram and, in general, the legitimacy of the use semi-quantitative approach for spatial analysis
of geostatistics for the determination of the risk of landslide risk at medium scales by combining
zones. Also, to allow for a more smooth risk zone, information on hazard analysis with the character-
a nugget effect ne equal to 50% of the sill parameter istics and function of the elements at risk as well as
is artificially added to the selected variogram. The their vulnerability. The proposed approach is quite
same variogram is used for all three risk zones. general in concept and depends primarily on data
To finally determine the risk zones, a value of Qb availability and the expertise of the user.
equal to 0.2 was selected. Once the risk value of each individual element
The risk zones are illustrated in Figure 3. No at risk (building in our case) is assessed, this study
high risk zone exists because the buildings that presents a geostatistical framework based on indi-
exhibit high risk are few and spatially scattered. cator kriging equations, which offers an automatic
Six moderate risk zones—the larger one at the procedure to objectively cluster the elements at risk
North-West side of the Barcelonnette Basin cov- in zones of equal risk. Although the procedure has
ers an area of about 445,000 m2—and multiple been applied to landslide hazard, it can be eas-
low risk zones are formed. The risk zone map ily extended to other hazards as well. The whole
highlights the prone areas on the basin. Indeed, in procedure can by fully implemented in an ArcGIS
the Barcelonnette Basin, the South-West moder- environment.
ate risk zone encloses a vast area where the most The resulting risk zone map offers an inte-
severe and densely-spaced damages are expected to grated risk representation necessary for efficient
483
484
ABSTRACT: Tunnel construction can cause deformations of the surrounding ground, which endanger
buildings and other structures located in the vicinity of the tunnel. The prediction of these deformations
and the damages to the buildings is difficult due to limited knowledge of the geotechnical conditions and
uncertainty in the response of the structures to the settlements. This paper presents a probabilistic model
for prediction of the damage to buildings due to tunneling, which combines the volume loss method with
the equivalent beam model. It furthermore proposes a probability-based method for determination of the
limiting value of settlement that is used for control purposes during the tunnel construction. Updating of
the limiting settlement with measurements gathered during the construction is described. The proposed
methodology is applied to a masonry building affected by the construction of the L9 metro line tunnel
in Barcelona.
485
⎛ y2 ⎞ s( y ) ⋅ y
s( y ) = Smax exp ⎜ − 2 ⎟ (1) u( y ) = (5)
⎝ 2i cos θ ⎠ z0
du( y )
where i is the location of the inflection point (hor- ε h ( y) = (6)
izontal distance from tunnel axis), θ is the angle dy
between the modeled plane and the perpendicular
plane and Smax is the maximum settlement in the The model errors Eε br d Eε ddr are considered
center of the Gaussian curve, i.e. above the tunnel as multiplicative RVs with mean value equal to 1.
axis. Smax can be calculated as: They result from the assumption of linear elastic-
ity, the position of the neutral axis and the omis-
sion of the presence of openings.
VL d 2
Smax f1(VL , K ) = (2) Maximum bending (εbmax) and shear (εdmax)
3.192 ⋅ K z0 strains in the equivalent beam are calculated as:
486
Limiting tensile
Category of damage Normal degree of severity Typical damage strain (εlim)(%)
height (t = H/2). In case of hogging deflection, the Reliability Method (FORM)) were used, separate
LSFs for ε ε br ε ε dr ε ε br , ε ε dr should be defined
sag sag hog hog
neutral axis is assumed to be at the top fiber (t = H).
Strains are calculated in the most critical fiber from and the failure event should be described as a series
the position of the neutral axis, so that a = t in both system.
cases. The damage on the buildings is determined
depending on the maximum strain εmax:
3 DETERMINATION OF LIMITING
SETTLEMENT
ε max = max ⎡⎣ε εsabrg ε εsag
dr
ε εhog
bbr
, ε εhog ⎤
dr ⎦
(9)
The measured maximal settlement above the crown
of the tunnel Sm equals:
where ε εsag
br
ε εhog
are the maximum bending strains
br
in sagging and hogging respectively,
p both obtained Sm Smax + E f Em = Smax EE (12)
using Eq. (3), and ε ε dr ε ε dr are the maximum shear
sag hog
487
This probability can be evaluated using a MC The proposed method is applied to a case study
simulation for different values of sm. The limiting of the L9 metro line construction in Barcelona.
settlement value slim ensuring Eq. (13) is then found The damage produced by the tunnel construction
iteratively. to a complex of masonry buildings from the late
488
489
Figure 2. Scatter plot of volume loss VL and max. Figure 4. Scatter plot of settlement Smax and max.
strain εmax. strain εmax.
Figure 3. Scatter plot of parameter K and max. strain Figure 5. Conditional probability of failure for different
εmax. values of measured settlement, Pr ( ε 0. % | sm ).
llim
490
5 CONCLUDING REMARKS
REFERENCES
The paper presented a computationally efficient
model for probabilistic prediction of building dam- Attewell, P.B., Yeates, J. and Selby, A.R. (1986) Soil move-
age due to tunnelling that is applicable in engineer- ments induced by tunneling and their effects on pipelines
ing practice (Sec. 2). Further, a novel method for and structures. Blackie Academic & Professional.
determining the limiting settlement was presented Boscardin, M.D. and Cording, E.J. (1989) Building
response to excavation-induced settlement. Jour of
(Sec. 3), which is a more systematic and traceable Geo. Eng. 115(1–21).
reliability-based approach with an explicit ration- Burland, J.B. (2008) The assessment of the risk of dam-
ale than the deterministic methodology typically age to buildings due to tunnelling and excavations. J.T.
used in practice. Additionally, the reliability-based Payma Cotas: Mov. de edificios inducidos por excava-
approach allows to incorporate measurements made ciones. p. 3.
during the construction. The proposed procedure Burland, J.B., Broms, B. and De Mello, V.F.B. (1977)
was demonstrated on a case study of a tunnel con- Behaviour of foundations and structures. Proc. 9th
struction in Barcelona (Sec. 4). First, the influence Int. Conf. on Soil Mech. and Found. Eng., 2, 495–546.
of the uncertainty in the model parameters (volume Burland, J.B. and Wroth, C.P. (1974) Settlement of build-
ings and associated damage. London: Pentech Press.
loss, shape parameter of the settlement through, Camós, C., Molins, C. and Arnau, O. (2012) A case
Young and shear moduli of the building material) study of damage on masonry buildings produced by
on the estimated damage was presented. Second, tunneling induced settlements. Int. Jour. of Architec-
the value of the limiting settlement was determined tural Heritage.
with the approximate approach as 22 mm. Third, Papaioannou, I. and Straub, D. (2012) Reliability updat-
the limiting settlement was determined more pre- ing in geotechnical engineering including spatial vari-
cisely using an advance realiability-based approach ability of soil. Computers and Geotechnics 42: 44–51.
as 23 mm. Both approaches provide similar values Peck, R.B. (1969) Deep excavations and tunneling in soft
and the simpler method appears to be satisfactory ground. SOA Report, 7th Int. Conf. SM&FE.
Straub, D. (2011) Reliability updating with equality infor-
for practical applications. Both of these values are mation. Probabilistic Engineering Mech. 26: 254–258.
more strict than the value that was used in the real TYPSA. (2003) Análisis de los movimientos del terreno
case, where a settlement of up to 24 mm was con- producidos por la excavación mecánica del túnel de
sidered to be safe. The reason for this difference is la L9 en la zona de Santa Coloma de Gramanet. Pro-
the fact that in the real case, the uncertainties in the jecte de construcció de la Línia 9 de metro de Barcelona.
ground parameters and building parameters were Tram 4t Bon Pastor—Can Zam. Infraestructura i Esta-
not considered and some unfavorable values of cions. Annex 8.
491
ABSTRACT: Many regions of the world are exposed to and affected by several natural hazards.
Implementation of effective risk management strategies in these areas requires that all relevant threats
are assessed and considered. Compared to single-risk analysis, the examination of multiple risks poses a
range of additional challenges due to the different characteristics of hazards. This paper summarizes pre-
vious research on multi-risk assessment and proposes a new three-level framework for multi-risk assess-
ment that could account for the possible interactions among the threats. The framework is developed as
part of the EU FP7 Collaborative Research Project MATRIX. The first level is a simple flow chart that
guides the user in whether a multi-hazard, multi-risk approach is required for the problem at hand. The
second level is a simplified, semi-quantitative approach to explore if a detailed assessment is needed. The
third level is a detailed quantitative multi-risk analysis based on Bayesian networks. The key components
of this framework, such as assessment of cascading hazards, time-dependent vulnerability estimation,
and the choice of the required level of sophistication are addressed in the paper. The multi-risk assess-
ment procedure outlined in the paper integrates the results of risk posed by each threat, cascade effect,
and appropriate consideration of uncertainties, to provide a rational estimate of multiple risks. Simple
examples that demonstrate the application of the method are presented in the paper.
493
494
495
496
497
REFERENCES
498
S. Moriguchi
Tohoku University, Sendai, Miyagi Prefecture, Japan
N. Asano
Chubu University, Kasugai, Aichi Prefecture, Japan
ABSTRACT: This paper presents a framework of risk management of rockfall using real survey data.
Hida area located in Gifu prefecture in Japan is selected as a target area. In the evaluation process of
rockfall provability, relative probability is calculated based on results of the survey data using the logistic
regression analysis. Then, absolute probability is calculated using history data of occurrence of rockfall
in the target area. In the evaluation process of the economic loss, road closure induced by rockfall is con-
sidered. Direct and indirect losses are taken into consideration. A trial calculation with different strate-
gies were conducted based on calculated risk, and then it was summarized the obtained results are useful
information for rockfall risk management.
Rockfall is one of the serious slope disasters. Once 2.1 Target area
rockfall take place, it has potential to cause dam-
Hida area, northern part of Gifu prefecture in
age on infrastructures and loss of human lives.
Japan (Fig. 1), is selected as a target area. In the
Japanese governments, researchers, and engineers
area, mountain areas are widely distributed, and
have been making effort to minimize the damages.
There are however still huge numbers of danger-
ous slopes that has high risk of rockfall in Japan.
Furthermore it is becoming difficult to construct
new countermeasures due to the budget cut of
public works, and maintenance cost of existing
countermeasures is becoming larger year by
year. It is therefore quite important to develop a
method which realizes a strategic decision mak-
ing of the maintenance and new construction of
countermeasures.
In order to overcome the problem mentioned
above, this study presents a framework of risk
management of rockfall. Risk of rockfall is cal-
culated based on real survey data. In this paper,
first, setting conditions, such as definition of
risk, study area, and assumptions are explained,
then calculation procedures of risk of rockfall
are explained. This paper finally shows a result
of a trial calculation is shown. Based on the
result, effectiveness of the proposed framework
is discussed. Figure 1. Study area.
499
500
Gero MR 186 51 5 7
OB 399 68 12
NM 204 10 0
Takayama MR 275 109 7 17
OB 517 380 7
NM 133 36 0
Furukawa MR 277 36 11 16
OB 126 20 5
NM 186 11 3
Figure 3. Histogram of rockfall probability (Gero). Figure 5. Histogram of rockfall probability (Gero).
501
502
503
effectively, but the strategy of plan B is not best such as fairness, are not considered in proposed
solution of risk management. The reason is that framework. In addition, as mentioned in this
the strategy depends on only economic efficiency, paper, some big assumptions are still included. It
and other viewpoints, such as fairness, are not is therefore necessary to develop the framework in
considered. It is therefore important to introduce the future works.
another indexes in future study.
REFERENCES
6 CONCLUSION Cabinet office (Japanese Government), A report of eco-
nomic analysis of damage of traffic accident, 2007. (in
This paper presents a framework of risk manage- Japanese)
ment of rockfall risk using real survey data. Road Hashimoto, T., Kurihara, M., Inoue, K., Iwasaki, Y. &
closure is defined as risk event, and the provability Fujimoto, A. 2002. Relationship between the carry-
of occurrence and economic losses are calculated. ing time and the survival rate in emergency medical
Finally, the trial calculation with different strategies service. Journal of Japanese Society for Emergency
of taking countermeasure was shown. According Medicine 5: 285–292. (in Japanese)
Public Works Research Institute (PWRI), A manual of
to the result of the trial calculation, it could be risk analysis and risk management (draft version),
summarized the obtained results are useful infor- Records of PWRI, No. 3926, 2004. (in Japanese)
mation for risk management of rockfall. The Ministry of Land, Infrastructure, Transport and
Although effectiveness of the proposed frame- Tourism, A manual of cost-benefit analysis, 2008. (in
work was validated, some important viewpoints, Japanese)
504
ABSTRACT: Rijkswaterstaat, the largest public client organization for infrastructural projects in
The Netherlands, is responsible for the design, construction, management, and maintenance of the main
civil infrastructure in The Netherlands. Their key success factor is realizing infrastructure projects, not
only within time, budget, safety and quality standards, but also with a minimum of hindrance for the
public during construction. For these reasons, Rijkswaterstaat initiated the Dutch Geo-Impuls programme
for reducing geotechnical failure. Consequently, this public organization embraced Geotechnical Risk
Management (GeoRM) as the leading geotechnical working method and is serious about implementing
GeoRM in its own organization. This paper presents the selected approach for the GeoRM implementation
process, which includes embedding GeoRM in existing formal working procedures and empowering
geotechnical engineers and managers to adopt and apply GeoRM in their day to day activities. The main
hurdles encountered and solutions applied are discussed. Finally, the paper draws the main conclusions
for implementing GeoRM in public client organizations in the construction industry.
505
506
507
508
509
REFERENCES
510
511
512
Done
1 Generate and protect value Make all geotechnical risks in each project phase explicit,
including risk effects and risk remediation measures
2 Participate in decision making Make a geotechnical risk file from the start of the project
in all project phases and use it for decision making
3 Make geotechnical uncertainty Include geotechnical sensitivity analyses with margins
explicit in project reports
4 Work systematically, Include GeoRM explicitly in the project planning and
structured, and in time reserve resources for it
5 Use all available information Work from a general level to a detailed level, from using
geological maps to geotechnical monitoring
6 Work transparently together Indicate and communicate any dependencies of geotechnics
with all stakeholders with other disciplines in the project
7 Include the role of the human Make differences in organizational culture of all involved
factor project parties visible and feasible
8 Use experiences and lessons Use all available and relevant project evaluations, risk
for continuous improvement checklists, and experiences from professionals
513
514
No. Action
No. Action
1 Give specific attention to geotechnical risks during project risk management sessions.
2 Organize internal and interdisciplinary GeoRM workshops for a number of projects.
3 Demonstrate to the internal organization how to apply interdisciplinary GeoRM, including its benefits.
4 Make geotechnical risks, including its inter-disciplinary remediation actions, explicit in the project risk registers.
5 Provide internal checklists for geotechnical risks that demonstrate their inter-disciplinary character.
515
No. Action
1 Apply internal lobby activities in the organization for getting GeoRM embedded in existing processes.
2 Include geotechnical engineering with GeoRM in primary process descriptions.
3 Provide a clear internal protocol for applying GeoRM.
4 Include the GeoRM process in the formal quality management protocol.
5 Formally embed GeoRM in the geotechnical department.
6 Consider together with the geotechnical consultancy team how to formalize GeoRM in projects.
7 Formally apply the relevant GeoRM tools in projects.
8 Design and apply formal GeoRM reporting structures.
9 Organize formal sessions about the GeoRM embedment in the organization.
10 Develop a formal scan for deriving the top 3 geotechnical risks in each project.
No. Action
1 Make non-technical decision makers enthusiastic by presenting the specific GeoRM benefits for them.
2 Identify geotechnical risks with the entire project environment and stakeholders at the start of the project.
3 Integrate geotechnical risks systematically and explicitly in contract between project parties.
4 Show to the project environment and stakeholders the need for and specific benefits of GeoRM.
Waterway GeoRiskScan 1, 4
Third Party Expertise
Motorway Risk-based soil investigation 1, 3
Bridge construction Geotechnical baseline report 1, 2
Tunnelling GeoRM Communication protocol 1, 4
516
5 CONCLUSIONS
REFERENCES
The three Dutch and industry-wide workshops, as
Bles, T.J., van Staveren, M.Th., Litjens, P.P.T. &
presented in this paper, demonstrate that the status Cools, P.M.C.B.M. 2009. Geo risk scan: getting
of GeoRM implementation, as well as the ways to grips on geotechnical risks. In Y. Honjo et al. (eds.),
increase the degree of GeoRM implementation, is Geotechnical Risk and Safety; Proc. 2nd Int. Symp.,
rather organization specific. Between and within Gifu, 11–12 June 2009. London: Taylor and Francis.
organizations, quite different degrees of GeoRM Cools, P.M.C.B.M. 2011. The Geo-impuls programme
implementation emerged by applying the quick reducing geotechnical failure in The Netherlands. In
scan on the degree GeoRM key conditions pres- N. Vogt et al. (eds.), Geotechnical Safety and Risk; Proc.
ence. Obviously, similar types of workshop can 3rd Int. Symp., Munich, 2–3 June 2011. Karlsruhe:
be easily organized in other countries with simi- Bundesanstalt fur Wasserbau.
ISO 2009. International organization for standardization
lar programs for reducing geotechnical failure, or 31000:2009 risk management—principles and guide-
within individual companies for strengthening their lines. Geneva: ISO.
GeoRM implementation capabilities and effects. van Staveren, M.Th. 2006. Uncertainty and ground
The results and examples in this paper demon- conditions: a risk management approach. Elsevier,
strate that a lot of pragmatic actions for embed- Oxford.
ding GeoRM are readily available. Therefore, it is van Staveren, M.Th. 2009a. Extending to geotechnical
simply up the responsible engineers and managers risk management. Georisk 3(3): 174–183.
to apply the quick scan for assessing the degree of van Staveren, M.Th. 2009b. Risk, innovation & change:
presence of GeoRM key conditions in their organ- design propositions for implementing risk manage-
ment in organizations. Köln: Lambert Academic
izations, and to select and apply suitable actions Publishing.
for strengthening any yet underdeveloped key van Staveren, M.Th., Litjens, P.P.T. & Heerema, J.J.
conditions. By doing so, considerable steps can be 2013. Implementing geo risk management in a client
made with implementing GeoRM in client organi- organization. In Geotechnical Risk and Safety; Proc. 4th
zations, engineering firms, knowledge institutes, Int. Symp., Hong Kong, 2–4 December 2013 (in press).
517
M.-H. Wu
Department of Civil and Environmental Engineering, National University of Kaohsiung, Kaohsiung, Taiwan
J.P. Wang
Department of Civil and Environmental Engineering, Hong Kong University of Science and Technology,
Hong Kong, China
ABSTRACT: In recent years, global warming has resulted in serious extreme climate events. Natural
disasters occur more frequently world widely, and the corresponding impacts are more and more
serious, for example of the tragedy of Hsiaolin Village, Kaohsiung due to the impact of Typhoon
Morakot in 2009. In this research, the study areas are primarily those who were suffering from typhoon
storm-related disasters and are characterized with high disaster potential. A series of field reconnais-
sance works were conducted in flood and slopeland disaster-prone areas in Kaohsiung city in order
to analyze the causative factors of the hazards. There are ninety eight communities and ninety nine
villages under investigation for flood and slopeland disaster potential analyses. With the environmen-
tal index buildup and assessment, the districts of Shanlin, Namaxia, Meinong, Maolin, Taoyuan and
Liugui are categorized as high slopeland disaster-prone districts, and the causative factors are exist-
ence of serious landslides and poor vegetated bare slopes. The districts of Yongan, Mituo, Gangshan,
Qiaotou, Nanzi, Dashe and Renwu are classified as high flood disaster-prone districts, and the causes
are abundant rainfall and poor drainage system. The results obtained from this study will contribute
to keep abreast of high disaster-prone areas and to promote the corresponding disaster prevention and
preparedness tasks.
519
520
521
522
score
Total
the study also reported that all of the mountain-
8
9
8
8
ous clusters have no substitute access roads, which
potential
means once the main roads are interrupted, the
Flood
clusters would be isolated.
zone
The statistical analysis for geographical features
1
1
1
1
flood prone of slope hazard prone areas showed that 81% of the
sites have more than one sensitive geological char-
Approved
facility
District
Poor drainage,
environmental
District score Vulnerability
pumping
dalay
Yong An 4 High
Mi Tuo 4
1
1
1
1
Gang Shan 4
overflow,
silting or
blocking
Hu Nei 4
Stream
Yan Chao 4
Lu Jhu 4
0
0
0
0
A Lian 4
Flooding
Da She 4
history
Ren Wu 4
Cian Jin 4
1
1
1
1
Cian Jhen 4
San Min 4
subsidence
or ground
Lowland
Ling Ya 4
Mei Nong 4
Village flooding indices grading form.
Zih Guan 4
1
1
1
1
reservoir spillway
Da Shu 3.8
Cie Ding 3.7
Niao Song 3.5
Tian Liao 3.3
Da Liao 3.3
1
1
1
1
Jiugang
Wei Sin
Gu Shan 3 Medium
Village
Ci Jin 3
Tadi
Liou Guei 3
Sin Sing 3
Table 3.
District
523
score
Total
6.2
7.2
assessment.
10
District
subsitute
environmental
access
road
District score Rank
No
1.0
1.0
1.0
debris flow Na Ma Sia 4 High
Mei Nong 4
dip slope
Potential
Mao Lin 4
creek,
0.0
1.0
1.0
Building foundation
subsidence, retaining
crack
0.0
1.0
0.0
or valley
0.2
1.0
0.2
5 CONCLUSIONS
Guanglin
524
525
K. Abbeche
Department of Civil Engineering, Laboratory LARHYA, University of Batna, Batna, Algeria
ABSTRACT: Collapsible soils types are found in natural unsaturated soils, which are characterized by
a metastable structure, and undergo an abrupt collapse when they are flooded (with or without loading),
which causes important damages. Based on the structural composition of these soils, samples made up of
various proportions of sand and fine particles were reconstituted. The first phase of work consists in the
experimental determination of the geotechnical characteristics of these samples. It is thereafter proposed
a test program based primarily on the use of the cone penetrometer and the ultrasonic apparatus as the
approaches for the prediction of the factors which influence collapse. The results obtained clearly show
the influence of certain parameters such as; water content and the energy of compaction on the collapse
potential and the ultrasonic speed.
529
Materials Characteristics
530
3 USED MATERIALS
531
4 TEST PROGRAM
532
533
The decrease of collapse is more obvious that 6.2.2 Interpretation of the penetrations
the moisture content increases. In the same condi- versus ωopt/ω0
tions of compactness and moisture content of the Collapsible soils are characterized by the condition
soil containing the greatest percentage of kaolin ωopt/ω0 > 1, (Holtz & Hilf 1961). Analysis of
exhibit greatest collapse potential. the penetration curves versus the ratio ωopt/ω0,
These results agree with those of Lawton et al. Figure 11 confirm the existence of two distinct
(1989) and Ayadat et al. (1998) and confirm the behaviors and separated by the line ωopt/ω0 = 1.
observations of Abbeche et al. (2007). One can On the left of this line, the penetration knows a
conclude that the reconstructed soils at the labora- gradual decrease for then growing in a roughly
tory hold a similar behavior to those met in situ, regular way as one moves away from the limit sepa-
therefore suitable for the series of tests suggested. rating the collapsible soils (ωopt/ω0 >1) of the non
534
535
536
537
ABSTRACT: This study represents an approach for investigating the effect of riverbed scouring
on piles of a reinforced concrete bridge in the event of an earthquake. A finite element model of box
girder bridge with 9 spans was developed for analysis of this phenomenon. Loading was generated
due to earthquake in along the traffic and perpendicular to the traffic direction. Scouring depth of
the pier and rived bed scour percentage with main river width are taken as parameter of the study.
Flexural stress in pile increases by 137% and axial stress was found increasing by 11% about original
no scoured condition for earthquake along the traffic direction while the flexural stress increases by
87% and axial load increases by 5% about original no scoured condition for earthquake in parallel
direction of traffic.
539
540
4.1 Assumptions
The investigation was carried out on the basis of
some assumptions to avoid complications. These
assumptions are as follows:
I. Materials involved in the study were elastic.
II. Only lateral force that was encountered by the
structure is the Earthquake.
III. In analysis for earthquake only dead load of
the structure was considered. No live load was
taken into consideration.
IV. Linearly increasing soil subgrade modulus
was considered.
Figure 4. A typical cross section of deck slab (all
V. Soil scouring profile around the bridge piers
dimensions are in mm).
was assumed to be parabolic.
VI. All the calculations were made on the basis of
pure loading. No scaling of load was done for
single point response spectrum analysis. Table 1. Material property.
541
EI
Stiffness factor, R = 4 (1)
kB
EI
Stiffness factor, T = 5 (2)
nh
Soil modulus
542
6 STUDY PARAMETERS
7 STUDY OBSERVATION
543
for 60% of the riverbed scouring and the value is scouring and the change in moment and axial load
almost 1850 kN-m. With no scouring condition in the piles. As the riverbed scouring progresses,
the moment is around 820 kN-m and the percent- the fixing length of the pile also moves down-
age increment of moment is 137%. Although the ward making the bridge more susceptible to
percentage increment of moment is more in earth- lateral loading. For that reason, in the event of an
quake in perpendicular direction of traffic but earthquake with scoured riverbed condition, the
earthquake in parallel direction of traffic is more moment in the pile increases and the axial load
vulnerable because of the higher moment value in decreases.
the latter case.
8 CONCLUSIONS
7.2 Axial load variation in corner pile
Figures 10 and 11 shows axial load variation, in From above results and discussions it can be con-
corner pile of maximum scoured pier, for earthquake cluded that bridge pier scouring incorporated with
in parallel direction of traffic and for earthquake earthquake will have adverse effect in the overall
in perpendicular direction of traffic respectively. integrity and health of the structure. So proper
It has been found that for 10% and 40% of river- protective measures should be taken for the soils
bed scouring, axial load in pile decreases for earth- around the bridge piers to prevent scouring.
quake both in perpendicular and parallel direction
of traffic. For 60% of riverbed scouring, axial load
in pile increases by a small margin. The increment REFERENCES
percentage is about 5% for earthquake in parallel
Pagán-Ortiz, J.E. 2003. Scour program update. Proc.,
direction traffic and 11% for earthquake in perpen- Midwestern Hydraulic Engineering Conf., U.S.
dicular direction of traffic about the axial load in Department of Transportation, Federal Highway
original no scouring condition. Administration, Washington, D.C.
From the above numerical results, it can be Tomlinson, M.J. 1994. Pile Design and Construction
inferred that there is a relation between the riverbed Practice. 4th Edition, E & FN Spon.
544
ABSTRACT: The excavation design for the new ship lock in Minden resulted challenging due to the
immediate vicinity of the 14 m deep excavation pit to existing structures and surface water bodies. An
alternative to a fully braced variant consisted of anchored support walls and slopes to the opposite side.
The hydrogeological setting consists of fill materials and Quaternary sands above bed-rock. Borehole pulse
and pumping tests revealed preferential flow along a system of fissures in the claystone in which river stage
fluctuations can propagate very fast. The excavation design in a setting with a fair hydraulic connection
to the river demanded an effective groundwater containment system and a robust dewatering. In order
to minimize the level of risk an active-passive dewatering scheme was developed. The excavation design
proved to be robust and yielded significant cost savings compared to a fully braced support walls variant.
2 GROUND CHARACTERIZATION
545
Figure 2. Cross section of the excavation pit for the new Figure 3. Arrangement of the containment system
ship lock in Minden. around the excavation pit.
546
6 DEWATERING SCHEME
Due to the limited anchor forces of 1000 kN accept- through the bored pile wall. Figures 7 and 8 show
able groundwater heads (according to the static schematically the active and passive dewatering
loading design of the bored pile wall) had to be approach respectively.
met. There are essentially two approaches to con- The passive operation has a considerable advan-
trol actions from groundwater on excavation sup- tage in the case of failure of the pump system
port walls. On the one hand, “actively” by means (e.g. by lightning strike etc.) as no inadmissible
of groundwater lowering through wells outside the head differences can build up behind the pile wall.
excavation pit, and on the other hand, “passively” In the worst case the pit is inundated in contrast
discharging the groundwater into the excavation pit to the active operation mode where a pumping
taking advantage of the groundwater head differ- breakdown can result in a rapid rise of groundwa-
ence outside and inside the pit. The groundwater ter heads and loadings on the excavation support
discharge and so the lowering can be considerably walls.
enhanced if groundwater inflow is short-circuited Further, a passive dewatering provides opera-
tional benefits since no redundant pump system
and permanent monitoring is needed as is the
case during active groundwater lowering to mini-
mize the risk of inadmissible static heads/loading
behind the support wall due to unanticipated fail-
ure of the pumps.
On the other hand a passive dewatering scheme
right from the start of the excavation on was not
appropriate. This would require groundwater
discharge through the bored pile wall at several
excavation levels. Therefore an active dewatering
based on vertical extraction wells (see Fig. 6) was
established during the excavation. After reach-
ing the final excavation level, the dewatering was
shifted from active to passive operation. To short-
circuit groundwater through the pile wall inclined
boreholes were drilled from the excavation base
to the vertical extraction wells and supplied with
a casing. Taking advantage of the head difference
outside and inside the pit, groundwater flows is
Figure 4. Schematic illustration of the dewatering sys- captured outside the pit by the vertical wells, dis-
tem during active operation at the final excavation level. charged through the wall in the obliquely-installed
547
548
Figure 11. Initial withdrawal rates [m3/h] at the start of Figure 12. Well drawdown [NN+m] at the start of the
active dewatering. Well drawdown at NN+27 m. passive operation.
549
550
ABSTRACT: In Hong Kong, it is a statutory requirement to provide a monitoring plan where the con-
struction works may affect any nearby building, structure, land, street or services. The traditional use of
the triple-A response values has many drawbacks when the projects are large and complicated. A modified
instrumentation monitoring plan is therefore proposed. Instruments are grouped to serve two purposes.
The first grouping is for the protection of the external sensitive receivers. They are categorized into types
for different response actions as different sensitive receivers can tolerate different extent of settlement
and angular distortion. The second grouping is for the performance review of the excavation and lateral
support systems. Another set of response values are determined from the predicted movements of the sys-
tems at various critical stages. Failure mechanisms are identified from the geological models and PLAXIS
analysis. Mitigation and contingency measures are considered in advance.
551
552
553
4.2 For the protection of the external sensitive and contingency measures are provided in advance
receivers by the Contractor and approved by the Engineer
and BD/GEO. During excavation, the design
For the external sensitive receivers, settlements
assumptions and design parameters are validated.
and angular distortions are the basic criteria
Monitoring instruments along a design section are
for the AAA values, which are the same as the
grouped for easy interpretations. Should the moni-
traditional AAA system and that recommended
tored ELS system movements and ground settle-
in the PNAP APP-137 (February 2012 revision).
ments deviate from the predicted trends, mitigation
However, different types of sensitive receivers can
measures are implemented. Should the monitored
tolerate different extent of settlement and angular
ELS system movements and ground settlements
distortion, and require different response actions.
reach the tolerable limits, contingency measures are
They can be categorized as in Table 1.
implemented and the Work is ceased immediately.
More thorough investigation report with further
4.3 For the performance review action plans are submitted to BD/GEO for approval
of the ELS systems before resumption of the Works.
In most situations, the estimated ground settle-
ment caused by the movement of the ELS system
5 MANAGEMENT OF THE
may not bear any direct relationship with the tra-
INSTRUMENTATION MONITORING
ditional AAA values for the ground settlements.
If, for example, the adjacent ground is bare or the
5.1 Management team
utilities are supported, the tolerable limit of ground
settlement can be larger than the traditional third In order to safeguard the unawareness of probable
A level, and the excavation can still be carried out catastrophic failure, it is the contract requirement
in a safe manner. Despite this, investigation reports to establish a joint monitoring team. The team
with mitigation measures are still required to be comprises the instrumentation team, geotechni-
submitted to BD/GEO and a new third A level has cal engineers, and the designer if necessary, of the
to be approved before resumption of the Works. Contractor, and the geotechnical engineers, con-
In the modified instrumentation monitoring struction engineers and detailed design consultant
plan, the concept of minimum expected value, from the Client, and an Independent Monitoring
most probable value, and maximum expected Consultant (IMC).
value in the MTR Practice Note PIM/PN/11-22
is extended such that a new set of AAA values for
5.2 Management process
design review are set from the predicted movements
at several critical stages. This is similar to the tol- During construction, the Contractor would
erable ground settlement limits given in the PNAP measure the instrument readings and upload
APP-137 (February 2012 revision). All probable them to a unified web-based instrumentation
failure mechanisms are identified, and mitigation database daily at a specified format agreed with
554
555
556
REFERENCES
557
Supandi
Mining Engineering Department, STTNAS University, Yogyakarta, Indonesia
H. Hidayat
Geotechnical Department, PT Borneo Indobara, South Kalimantan, Indonesia
ABSTRACT: Bedding geometry is important in supporting the slope stability, particularly on low-wall
part of coal mining. Elements that are part of geometry are bedding slope, slope height and bedding ratio.
The key to this analysis is inserting geology conditions mainly bedding geometry. Bedding ratio value is
calculated based on modeling geology in which more detailed modeling geology factor can be identified
better analysis in slope stability. Slope stability analysis is assembled based on correlation between dip
of bedding, bedding height and factor of safety. The result from this study is increasing bedding ratio
value will be followed by increasing the factor of safety and bedding height. Changes in dip of bedding
slope will impact on the slope stability and respond changes on factor of safety and bedding height. The
final of the study is providing geo-mechanic parameters purposed for mine planning design based on site
characteristic geology.
559
560
561
562
H. Wang & Z. Wu
Ranken Railway Construction Group Co., Ltd., Chengdu, China
D. Wang
Tianjin Institute of Urban Construction, Tianjin Key Laboratory of Soft Soil and Engineering Environment,
Tianjin, China
ABSTRACT: More risks associated with subway transfer station construction than typical
excavations due to larger depth and more interactions. Besides elaborate design and careful arrangements
for constructions, monitoring is a critical process to control, predict and mitigate risks emerging possibly
during construction, as well provides updated information to revise the design in a timely manner. This
paper introduces a famous subway construction project-Jiaomen West Station, a transfer station for the
existing subway Line 4 (M4) and the proposed new subway Line 10 (M10) in Beijing. Based on the project
information, this paper identifies the risks associated to the project, introduces the methods to reduce
and mitigate the risks; monitoring data are presented and discussed in association with the mitigation
methods. Information presented in this paper and lessons learned from this case study are useful to the
similar projects.
563
564
Table 1. The geological conditions of Jiaomen West 4. Construction management and supervision
transfer station. Poor construction management may destroy a
good project. Not only it is likely to reduce the
Layer Descriptions project profits but also may cause project acci-
Backfill, clay mixed silt, bricks, cobbles, loose dents or failure. Good supervision, however,
2
and inhomogeneous may prevent this happening.
Silty clay, yellow-brown, soft to stiff, wet, 5. Construction monitoring
2
highly compressible Underground work has high potential to cause
Cobble and gravels mixed with fine sands, settlement of ground, displacement of vertical
medium dense, wet, size 2–5 cm, max. size strut, deflection of bracing system, damage to
10 cm, the utility pipes, etc. Rigorously monitoring is
1 Clayey silt and silty clay, yellow, medium stiff, very helpful to indentify and control such risk,
medium dense, medium compressible finally may make them avoidable. However, poor
2 Fine to medium sands, yellow, dense, wet, low monitoring would not achieve the expected role.
compressible
Cobbles mixed with fine and medium sands,
mixed color, dense, wet, low compressible.
Size 2–6 cm, max. size 10 cm, 3 MITIGATION MEASURES:
1 Clayey silt and silty clay, yellow, medium stiff, CONSTRUCTION AND MONITORING
wet, low to medium compressible
PLANS
2 Fine to medium sands, yellow, dense,
saturated, low compressible
Cobbles mixed with fine and medium sands, 3.1 Construction plans
mixed color, dense, wet to saturated, low Bench-cutting method was proposed to undercut the
compressible, size 4–8 cm, max. size 10 cm soils below the bottom slab of Jiaomen West Station
for M4. The maximum progress length for each
tion bracing system design, strut selection and excavation cycle was limited to 4 m. In order to sta-
design as well as dewatering, particularly, the bilize and facilitate the soil excavation, the site was
work mentioned above are more associated with dewatered with stages described below; a shaft in east
geological and geotechnical conditions as well side of the station was built to provide construction
as construction workmanship. space and mobilize personnel and equipment to the
3. Construction experience and workmanship designed ground level; surroundings soils were then
These involve in every process of the project stabilized by grouting cement-silicate mortar into the
construction. Particularly, the construction ground along the long PVC sleeve valve pipes and
workmanship for the engineering work men- pre-support systems, utilizing short perforated pipes
tioned in item 2 is more important, even critical, to form umbrella vault, were installed. The details of
to the entire project. the construction plan are described below:
565
Figure 3. Stabilizing the soils surrounding the twin tun- Figure 5. Umbrella vault system and soil stabilization
nels for M10. (longitudinal section view).
566
design. Table 2 summarizes the items to monitor Table 3. Criteria for alarm warning issuance.
and the control criterion utilized for monitoring in
construction. Due to the length limitation of the Alarm
paper, the plan of the monitoring points is shown warning level Criterion for alarm warning issuance
in Figure 1 partially. In the plan, DB means the Yellow Both absolute value and the rate of the
points monitoring the settlement of ground, JCJ monitored data exceeding 70% of
means the points for monitoring the building set- the control criterion, or one of them
tlement and GXC means the points for monitoring exceeding 85% of the control criterion.
the settlement of the utility pipes. Orange Both absolute value and the rate of the
Project management classified the monitor- monitored data exceeding 85% of
ing data to three levels to issue alarm warning, the control criterion, or one of them
1) yellow; 2) orange; and 3) red. The classification exceeding 100% of the control criterion.
is based on the absolute value and the rate of the Red Both absolute value and the rate of the
monitored data. Table 3 summarizes the criterion monitored data exceeding 100% of
the control criterion, or the rate of the
of the three levels for warning issuance. monitored data varies dramatically.
567
subsidence accordingly, the maximum settlement construction arrangement are important to a safe
occurred right behind the cutting face at point and successful project; 2) monitoring is a critical
DB-05-3, which is above the right twin tunnel. process to construction; 3) covered undercutting
Monitored data also shows the subsidence of has more risk potential that may lead to more
ground around the east and west exaction pit. ground subsidence than open excavations; and
However, the rate of the settlement around excava- 4) with respect to covered undercutting, construc-
tion is small and the maximum settlement occurred tion quality and overburden depth may play a role
at DB-13-03 (east part of the excavation pit of the for ground subsidence. Lessons learned from this
east hall), but not exceeding the control criterion. project are: 1) every underground project is unique,
Analysis indicated that the backfill behind the even in the same project, covered undercutting for
bracing system is not compacted well and is more transfer passage and twin tunnels may lead to dif-
influenced by the excavation, the consolidation ferent settlement; 2) attention should be paid to
of the backfill also contributed to the settlement more details no matter how small the project is.
observed. As for the movement of track and track Actually, auxiliary parts of the project may deserve
bed structures, it was caused by the subsidence of more attention though people usually do to big
the soils surrounding the twin tunnels. Such move- and main part.
ment is probably inevitable as long as the ground
loss occurred, and probably, the control criterion REFERENCES
for such movement will govern the stabilization
design for the covered undercutting. Deng, X.P., Li, Q.M. & Zhou, Z.P. 2010. Statistical analy-
One of the interests is the monitoring data for sis of the accidents of subway construction in China.
transfer passage (which was not shown in Fig. 1 Chinese J. Statistics and Decision 9: 87–89.
due to length limitation of this paper) exceeds the Gatti, M.C. & Cassani, G. 2007. Ground loss control in
criterion. The dimension of the cross-section, the EPB TBM tunnel excavations. In J. Bartak, I. Hrdina,
treatment and excavation methods for the transfer G. Romancov & J. Zlaml (eds), Underground Space-
passage are similar with those applied in twin tun- the 4th Dimension of Metropolis; Proceedings of the
nels, except the overburden depth is almost half of World Tunnel Congress 2007 and 33rd ITA/AITES
Annual General Assembly, Prague, May 2007. London:
the covered undercutting for the twin tunnel. Thus, Taylor and Francis Group.
the reasonable explanations could be 1) shallower Huang, M.T. 2011. Dewatering practice for a deep sub-
overburden depth lead to more obvious ground way station excavation. Chinese J. Modern Urban
subsidence observations and 2) poorer construc- Transit 3: 58–61.
tion quality lead to more subsidence observed. Hulme, T.W., Shirlaw, J.N. & Hwang, R.N. 1990.
Settlement during the underground construction of
the Singapore MTR. In Proc. 10th South East Asian
5 CONCLUSIONS Geotechnical Conference, Taipei, 1990.
Li, T., Xu, L.P., Chen, H. & Gu, G.R. 2006. Quantitative
Analysis of risks of foundation excavations. Chinese J.
This paper identifies the risks associated with a Geotech. Eng. 11(28): 1917–1920.
subway transfer station construction in Beijing. Sejonaha, J., Jaruskova, D., Spackova, O. & Novotna, E.
Mitigation plans are introduced and the monitor- 2009. Risk quantifications for tunnel excavation
ing data are presented. Conclusions are made from process. World Academy of Sciences, Engineering, and
the project that 1) elaborate design and careful Technology 58: 393–401.
568
N.R. Wightman
Sinclair Knight Merz (HK) Limited
Formerly SMEC Asia Ltd.
A.D. Mackay
Nishimatsu Construction Company Limited
Formerly SMEC Asia Ltd.
ABSTRACT: To increase power output from the Magat Pumped Storage Power Station, Luzon Island in
the Philipines, additional power generating equipment will be installed immediately adjacent to the existing
pump storage system. A concept design for the open cut excavation needed to install the equipment involved
rock excavation to an average of 60 m below the existing ground level with a 22 m span. To allow equip-
ment placement and future maintenance access into the excavation a 325 tonne load capacity gantry crane
spanning the excavation was proposed. The crane will run on rails placed onto rock with sufficient bearing
capacity to support the dead and live loads exerted during operation. Based on a site inspection, literature
search and preliminary ground investigation the rock mass bearing capacity will depend upon the presence
and characteristics of the rock joints. These included potentially persistent joints, continuing the length
of the excavation, with calcite coating and an adverse orientation with respect to the excavation, trending
parallel to the sides of the excavation and dipping about 60 degrees into the excavation. The calcite poten-
tially reduces the joint surface shear strength to a friction angle of about 20 degrees. To increase the shear
resistance it is intended to install six high strength proprietary rock anchors (each with 16 strand multi-core
tensioned cables and capable of imposing loads over 2,200 kN) perpendicular to the joint orientation. This
paper summarizes the anticipated ground conditions, the potential bearing capacity of the ground support-
ing the crane and the ground anchor installation proposed to improve the rock mass strength.
569
570
excavation was therefore positioned with sufficient Figure 8. Areas of major jointing in the vicinity of the
site.
space to allow the installation of the cable anchors.
Figures 6 and 7 present the proposed open cut
excavation in relation to the existing dam.
571
572
573
Upper bound finite element method for seismic stability limit analysis
of rockfill dams
ABSTRACT: In this study, an upper bound limit analysis Finite Element Method (FEM) is developed
to study the seismic stability of rockfill dams. Considering the large value of the internal friction angel
and the non-linear shear strength parameters of rockfill materials, a static form, which is corresponding
dual Second-Order Cone Programming (SOCP) problem of the upper bound limit analysis, is formulated
with constraints based on the yield criterion, flow rule, boundary conditions and the energy-work balance
equation. The upper bound solution of critical seismic coefficient is then obtained by a state-of-the-art
interior-point algorithm iteratively. Compared with the limit equilibrium method, the proposed method is
applied to the seismic stability problem of a typical rockfill dam. The results demonstrate the correctness
and practical value of the proposed method. From the displacement field obtained by solving the SOCP,
we can predict the failure surface of the dam slope easily.
575
p∑σ : ε
d p ( ) sup (4)
σ ∈F
E { d +∞ } (5)
576
where
where D is the dimension of the tensors and δ is
Kronecker’s δ . For plane strain conditions of the ci ℜ3 c iT = ⎡⎣ ki 0 0 ⎤⎦
Mohr-Coulomb criterion, D = 2, a = sin ϕ , and
k c cos ϕ , where c is the cohesion and ϕ is the z i ∈ℜ3 z iT = ⎡⎣ λi ( e ired )T ⎤⎦
ℜ3×3 A i = diag [ai 1 1]
internal friction angle.
Ai
T
ℜ3× NZ Bi = ⎡⎣ B m,i
Bi ∈ℜ B d,i ⎤⎦ (15)
3 FEM FORMULATION FOR UPPER
BOUND LIMIT ANALYSIS NZ is the total number degrees of freedom and
z i ∈ K i is the second-order cone constraint. Sets K i
Considering a plane strain structure divided into are quadratic cones of the following form:
NE finite elements, the optimisation problem
Equation 6 can be transformed as follows: K { ℜd : ≤ x1, x1 ≥ } (16)
2 :d
NE
min ∑ ∫A i
λdAi − q T0 u Equation 14 is a standard SOCP problem, which
can be transformed into the corresponding dual
i =1
ai λi = B m
m,i u (i = 1, , NE) max kc
e ired B d,i u (i = 1, …, NE) s.t. ( ym,i , sired ) ∈ K i (i = 1, …, NP)
q1T u =1 (11) ym,i + ai m ,i = ki* (i = 1, …, NP)
NP NP
where ∑ G m,
m σ m ,,ii ∑ Gd,i,i sired β q1 = q 0 (17)
i =1 i =1
θi = B ,i u (12)
where
T
⎡2exx 2exy ⎤⎦ = B dd,i u (13) G m, B Tm,i
⎣
G d, B dT,i
The matrices B m,i and B d,i can easily be obtained
from the typical strain-displacement relations. In ki* ki / ηi
geotechnical engineering, the displacement bound- ηi = Aee,,i / 3 (18)
ary is typically that u = 0 on S. For conciseness, we
assume that the above-mentioned conditions have It has been shown that the dual formulation
already been satisfied in Equation 11 and in what Equation 17 can be considered as a static form
follows. of upper bound analysis by Makrodimopoulos &
As mentioned above, the flow rule only needs to Martin (2007a). The main advantages of solving
be enforced at the three vertices of each six-node the dual problem of the upper bound limit analysis
triangular element to ensure that it holds through- are not only the dual form is robust and stability
out the element. For this reason, we take the ver- but also the stress distribution of structure can be
tices of each triangle as flow rule points. Thus, for obtained, which is used to describe the nonlinear
an NE-element structure, the total number of the failure criterion of rockfill materials.
flow rule points must be NP, where NP = 3NE. At present, large-scale SOCP problems can be
Now, Equation 11 can be formulated as follows, solved effectively using the interior-point method.
and more details can be found in the work by Based on this method, various state-of-the-art
Makrodimopoulos & Martin (2007a): algorithms are developed by researchers in
577
578
ACKNOWLEDGMENT
kc
REFERENCES
Project 1 Project 2 Project 3
Bishop, A.W. 1955. The use of the slip circle in the
Coarse mesh 1.139 0.948 1.369 stability analysis of slopes. Géotechnique 5(1): 7–17.
Fine mesh 1.128 0.945 1.349 Chen, L. & Chen, Z. 2007. Effect of nonlinear strength
Bishop’s method 1.09 0.92 1.32 of rockfill on slope stability of high earth-rock dam.
Rock and Soil Mechanics 28(9): 1807–1810.
579
580
an informa business