Geotechnical Risk Management and Reliability Based Design-Lessons Learned
Geotechnical Risk Management and Reliability Based Design-Lessons Learned
Geotechnical Risk Management and Reliability Based Design-Lessons Learned
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 99
(working) stress design and gain acceptance of it by geotechnical practitioners. The intent was to
keep things initially simple in order for practitioners to better understand fundamental concepts
and theoretical principles. As practitioners become more accustomed with LRFD and reliability
based design concepts, future code editions would introduce sophistication.
The degrees of sophistication are evident when the CHBDC (2014) is compared with the
previous editions in 2000 and 2006. Examples of sophistication include the introduction of
geotechnical resistance factors based on Level of Understanding of site and ground conditions
and Consequence Factor to account for the consequence of failure of a structure (Fenton et al.
2016). The rationale in support of these two significant and fundamental changes is that a higher
geotechnical resistance factor should apply when the geotechnical engineer has an improved
level of understanding of ground conditions due to more comprehensive site investigation and
analyses, or when the consequence of failure is lower than the typical case. Previously, a single
geotechnical resistance factor for ultimate limit states applied regardless of level of
understanding and consequence of failure (e.g., for bearing resistance of shallow foundations, the
geotechnical resistance factor was 0.5 and for pile axial compression resistance, the geotechnical
resistance factor ranged from 0.4 (static analysis) to 0.6 (static load test). CHBDC (2014) also
introduces geotechnical resistance factors other than 1.0 for serviceability limit states. In earlier
versions of CHBDC and NBCC, geotechnical resistance factors for serviceability limit states
were taken as 1.0.
For details of the changes and revised geotechnical resistance factors, the reader is referred to
CHBDC (2014) and Fenton et al. (2016). For example, for bearing resistance of shallow
foundations, the geotechnical resistance factor now varies from 0.45 (low understanding) to 0.65
(high understanding). For pile axial compression resistance, the geotechnical resistance factor
ranges from 0.35 (low understanding - static analysis) to 0.7 (high understanding - static load
test).
The next edition of the NBCC and CFEM will similarly be updated to reflect the provisions in
CHBDC (2014). This will be done so that consistency (harmony) is obtained amongst Canadian
codes and the CFEM, which is frequently referenced by the codes.
This incremental approach appears to have worked reasonably well, though implementation
issues such as those presented and described in this paper exist. In any event, for the reasons
presented in Becker (1996), the author believes that the factored overall resistance LRFD
approach has been much better received and accepted by geotechnical practitioners than the
factored strength approach. The author is of the opinion that with sufficient time and experience,
practitioners will feel increasingly comfortable with LRFD and realize its benefits over that of
allowable (working) stress design. It is noted that structural engineers also went through similar
challenges, albeit almost 40 years ago, when structural engineering design switched from
allowable (working) stress design to limit states design (LRFD).
LRFD applied to retaining walls has also received significant attention over the past 10 to
20 years. Publications such as Bathurst et al. (2012) and others summarize this work.
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 100
LRFD FRAMEWORK
The general LRFD design equation is:
ФR ≥ Σα S [1]
where ФR is the factored geotechnical resistance, Ф is the geotechnical resistance factor (values
of less than 1.0), R is the nominal (characteristic) geotechnical resistance, Σα S is the
Downloaded from ascelibrary.org by University of California, San Diego on 06/02/17. Copyright ASCE. For personal use only; all rights reserved.
summation of factored load effects for a given load combination, and α is the load factor
(usually greater than 1.0) for nominal load effects S (e.g., dead load due to weight of structure
or live load due to wind), and i represents various types of loads such as dead load, live load and
wind load.
The design equation can be visualized by inspecting the interaction of the probability distribution
curves for resistance and load effects, as shown schematically in Figure 1. The design intent is to
achieve a specified reliability index (or probability of failure) that is related to the size of the
shaded area shown in Figure 1, which corresponds to a failure condition (i.e., resistance is less
than applied loading). It should be noted that the resistance and load effects are assumed to be
independent variables, which is approximately true for the case of static loading condition
associated with foundations. The characteristic values for design are related to the mean values
through the factors kr (the ratio of mean value to characteristic value for geotechnical resistance)
and the factor ks (the ratio of mean value to specified (characteristic) value for load effects).
Typically, kr values are equal to or greater than 1.0 and ks values are less than 1.0. The terms kr
and ks are also referred to as bias factors. The bias factor is 1.0 if the mean value is used as the
characteristic value or when the predicted mean resistance is the same as the measured mean
average resistance.
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 101
The values of Ψ range from 0.9 to 1.15 for high and low consequence structures, respectively.
For typical consequence structures, Ψ = 1.0. The rationale for Consequence Factor and its values
are provided in CHBDC (2014) and Fenton et al. (2016).
IMPLEMENTATION ISSUES AND CHALLENGES
Since LRFD-based codes became mandatory in Canada, the following key implementation issues
and challenges have arisen:
Downloaded from ascelibrary.org by University of California, San Diego on 06/02/17. Copyright ASCE. For personal use only; all rights reserved.
• Some owners have observed that costs of foundations and geotechnical systems have
increased relative to those produced previously by allowable (working) stress design.
• Geotechnical practitioners do not appear to fully understand the principles of LRFD
design and reliability based design. In particular, there appears to be a lack of
understanding of the characteristic value, resulting in inappropriate selection of
geotechnical characteristic values for design. Geotechnical practitioners seem to lack
confidence in the LRFD approach and, as a result, they generally tend to select overly
conservative characteristic values that result in conservative designs and increased costs.
• There appears to be a lack of effective communication and interaction between structural
and geotechnical engineers. This has contributed to the geotechnical practitioner’s angst
when providing what seems to them to be very high factored bearing resistance values
for ultimate limit states, which are of the order of three times higher than typical
allowable bearing pressure values. The tendency of practitioners is to provide values
lower than what is truly appropriate for ultimate limit states. A reason for this is that
allowable (working) stress design does not explicitly separate ultimate and serviceability
limit states design recommendations, though this distinct consideration is inherent in
limit states design. In the author’s opinion, a large portion of geotechnical practitioners’
angst is that they do not know what the structural engineer will do with the apparently
very high factored bearing resistance values. However, this angst should be easily
resolved because in many cases serviceability limit states considerations often control
design. Allowable bearing pressure design values are frequently very similar to
serviceability limit states design values.
• Both structural and geotechnical engineers have been reluctant to adopt and embrace
LRFD for foundations and geotechnical systems. The author has experienced structural
engineers requesting allowable (working) stress design parameters instead of LRFD
recommendations. This is not conducive to the encouragement of geotechnical
practitioners to adopt and embrace LRFD design.
• It appears that most owners and their consultants do not realize that Canadian codes
allow the development of project and site specific load and geotechnical resistance
factors, based on an appropriate reliability index (probability of failure) for the project.
However, the size and scale of project needs to be sufficiently large or special to warrant
the investment made to develop project and site specific values.
The following sections describe case records (projects) of the author to demonstrate some of the
above implementation issues and challenges. However, before doing so it is useful and insightful
to examine briefly statistical methods and their use in reliability based design.
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 102
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 103
These equations are used in the next section to estimate approximate mean value when maximum
value is known.
FROST UPLIFT ON LIGHTLY LOADED PILES
This implementation issue discusses the standard practice of geotechnical engineers tending to
report the maximum frost depth that would be anticipated within the design life of a project. This
practice is not consistent with the intent of the characteristic value being the mean value or close
Downloaded from ascelibrary.org by University of California, San Diego on 06/02/17. Copyright ASCE. For personal use only; all rights reserved.
to the mean value of the expected frequency distribution for frost depth. The geotechnical
resistance factor for uplift of 0.3 given in the NBCC and CHBDC is based on the characteristic
value being the mean value or close to the mean and the specified load factors for type of load
and load combinations (e.g., in the range of 1.25 to 1.7).
LRFD design of lightly loaded piles to resist uplift due to frost action is an example of
misalignment and misunderstanding of the appropriate selection of characteristic value for frost
depth in geotechnical reports. For sites with reasonably deep frost depth, the design of lightly
loaded piles, such as piles to support pipe racks in petroleum process facilities in Alberta, is often
controlled by resistance against frost uplift. Traditional allowable (working) stress design has
worked well and led to an empirical design guidance that, in most cases, the length of a pile
(with minimal externally applied axial compression load) to resist frost uplift is approximately
three times the anticipated maximum frost depth. In Alberta, the required length of pile was often
in the range of 5 m to 10 m, depending on maximum frost depth.
Following the implementation of LRFD-based codes in Canada, the author has received many
phone calls and emails from structural engineers saying that when they followed the code they
calculated excessively long piles to resist frost uplift. In some cases piles as long as 20 m or
longer were calculated. Clearly something is not right.
The reason for the outcome of very long piles is that the design values given for frost depth in
geotechnical reports are usually based on maximum (extreme) values that lie in the tails of their
probability density functions (Figure 2). The estimated maximum frost depth is an important
parameter for the design of buried infrastructure, such as unheated water supply lines; however,
the maximum frost depth is not directly appropriate for pile design to resist frost action uplift.
The basis of the derived geotechnical resistance factor of 0.3 for uplift action on a pile in the
CHBDC and NBCC is that characteristic values for adfeeze (skin friction of frozen soil) and
frost depth be based on its mean value or close to the mean value, not the estimated maximum
value.
In order to resolve the situation for the excessively long piles, the author’s immediate response
was that, for the reported values of maximum frost depth, a geotechnical resistance factor and
load factor of 1.0 be used when checking for resistance against uplift to frost action. The reason
for this is that resistance and load factors of 1.0 seem reasonable when extreme (maximum)
values are used for the characteristic value for frost depth. The value of 0.3 is appropriate when
uplift is caused by external loads such as wind and other loads leading to eccentric loading, used
in combination with code-specified load factors.
Frost depth is a function of many factors including soil type, water content, degree-freezing days,
snow cover, albedo and others. A number of relationships to estimate frost depth has been
developed and reported in the technical literature (e.g. CFEM 2006). Uplift action on a pile is a
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 104
function of the frost depth and the adfreeze acting of the circumference of the pile shaft within
the frost depth (Figure 3) and is given as:
Fu = fadfdC [5]
where Fu = uplift force due to frost action (kN), fad = adfreeze (kPa), fd = anticipated frost depth
(m), and C = circumference of the pile shaft within the frost depth (m).
Downloaded from ascelibrary.org by University of California, San Diego on 06/02/17. Copyright ASCE. For personal use only; all rights reserved.
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 105
From Equation [7], the following equation is developed for the case of maximum frost depth
fdmax:
ФL ≥ (fad/fs)fdmax [8]
Frost action uplift is a load and the structural engineer would apply a load factor (e.g., = 1.25)
to the uplift force due to frost action (Fu). If Ф = 0.3, = 1.25 and (fad/fs) ranging from 1.4 to 2.5
are put into Equation [8], L ranges from ≥ 5.8fdmax to ≥ 10.4fdmax, which is consistent with
Downloaded from ascelibrary.org by University of California, San Diego on 06/02/17. Copyright ASCE. For personal use only; all rights reserved.
structural engineers calculating pile lengths that are much longer than pile lengths based on
empirical guidance (i.e., L ≥ 2fdmax). For example, if Ф = 0.3 and α = 1.25 with maximum frost
depth of 2.5 m is used in the calculation, L becomes 20.8 m and a 23 m long pile is calculated to
be required, which is excessively long and not appropriate.
If Ф = = 1.0 and (fad/fs) = 2, Equation [8] becomes L ≥ 2fdmax and is consistent with empirical
guidance experience. This seems reasonable because for extreme loading cases, resistance and
load factors of 1.0 generally apply. In order to obtain L ≥ 2fdmax when Ф = 0.3, needs to be
0.3. This is a very low load factor that would probably cause angst for structural engineers.
The above approach is examined below in more detail using reliability design concepts to assess
appropriate values for geotechnical resistance and load factors.
The First Order Second Moment (FOSM) method (Becker 1996, FHWA 2001), is a simple
approximate equation that is considered to capture the essence of the problem and serves as an
approximate basis to theoretically interrogate the suitable value of geotechnical resistance factor
and load factor for frost uplift loading condition. The FOSM relationships for load and
geotechnical resistance factor is as follows (Becker 1996):
Ф = kre-θβVr [9]
= kseθβVs [10]
where is a separation factor of 0.75, kr and ks are bias factors (typically in range of 1.0 to 1.1
for kr and 0.9 to 1.0 for ks), β is the target reliability index (typically 2.5 to 3.0), and Vr and Vs
are the coefficient of variation of the resistance and load distributions, respectively.
Although the author does not have data to support the opinion, it may be reasonable to consider
that coefficient of variation for extreme loads would be small. If it is assumed that the coefficient
of variation is small (say 0.05) the value of from the FOSM Equation [10] becomes
approximately 1.0. This provides additional support that = 1.0 is appropriate for design based
on maximum frost depth.
As an example of an initial step for examining and assessing an appropriate and consistent
methodology for designing lightly loaded piles to resist uplift forces due to frost action, the
climatic data for Calgary from 1970 to 2010 was used to predict annual values of frost depth for
a typical clayey/silty soil in Calgary. The frost depths would be greater for sandy and gravelly
soils. The freezing index in terms of Celsius degree-days is a key factor in predicting frost depth.
The freezing index varied from about 250 Celsius degree-days to about 1,450 Celsius degree-
days, with a mean value of about 800 Celsius degree-days (Figure 4). The corresponding
predicted frost depth varied from about 0.7 m to 2.3 m as shown in Figure 5. The mean value is
approximately 1.4 m, with a standard deviation of 0.37 m and a coefficient of variation of 0.26.
Substituting these values into Equations [3] and [4] results in a maximum frost depth of 2.5 m.
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 106
The results of this analysis agrees with the state-of-practice in Calgary where a maximum frost
depth of about 2.5 m is often stated in geotechnical reports.
Downloaded from ascelibrary.org by University of California, San Diego on 06/02/17. Copyright ASCE. For personal use only; all rights reserved.
Figure 4. Histogram and normal distribution – annual freezing index, Calgary, Alberta.
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 107
Given the above examination, the following design options become available so that excessively
long piles would not be calculated to resist frost action uplift:
• Use Ф = 1.0 and α = 1.0 with maximum frost depth, fdmax.
• Use Ф = 0.3 (possibly 0.4) and α = 1.25 with mean frost depth, fdmean.
The CHBDC and NBCC technical committees are in the process of rigorously examining and
assessing the most appropriate and consistent methodology for designing lightly loaded piles to
Downloaded from ascelibrary.org by University of California, San Diego on 06/02/17. Copyright ASCE. For personal use only; all rights reserved.
resist uplift forces due to frost action. The author is of the opinion that it makes sense that the
characteristic value for frost depth be the mean or close to the mean value to be consistent with
current LRFD methods. If this is the approach is to be adopted in the CHBDC, NBCC and
CFEM, geotechnical engineers will need to also report the mean value of frost depth instead of
just the maximum frost depth. In the interim, geotechnical reports could state that for the
maximum frost depth reported, a load and resistance factor of 1.0 should be used in the
calculation of required pile length to resist uplift forces due to frost action.
CHARACTERISTIC VALUE
In the author’s experience, the inappropriate selection of characteristic value by geotechnical
practitioners is the greatest implementation issue in connection with LRFD-based codes. The
following discussion is intended to assist the reader in better understanding how to select
characteristic values for design.
Although the values of the geotechnical resistance factors are important, the real essence of
design and the key question is: What is the appropriate characteristic value? The definition and
basis for the appropriate selection of the characteristic value has not received the attention it
rightly deserves. The technical literature and design codes in general give little or no guidance
for the appropriate selection of geotechnical characteristic values. In the author’s opinion, the
assessment and selection of the appropriate characteristic value is a most important part of the
geotechnical design process.
A description and discussion of the characteristic value is provided by Becker (1996) and CFEM
(2006). The characteristic value is defined as the geotechnical engineer’s best assessment of the
most likely representative (unfactored) value that controls a specific limit state. The value needs
to account for all factors that have influence on the parameter or property within the volume of
ground (zone of influence) under consideration.
Selection of Appropriate Characteristic Value
Let’s interrogate the above definition to gain a better sense of what the words mean. From the
above definition, the characteristic value for undrained shear strength may not directly be the
mean value or a value close (cautious estimate) to the mean of the actual field or laboratory
measurements. Measured values often require adjustments to take into account factors that have
influence on the measured parameter. For example, field vane strength for high plasticity clay
should be adjusted using Bjerrum or Aas correction factors (CFEM 2006). In clay deposits
consisting of regular interlayers of stiffer and weaker soil, such as varved clay, the field vane
strength may measure a strength that is too high because the vertical failure surface imposed by
the vane passes through the stiffer and weaker layers (Becker et al. 1988). A slope stability
failure surface often has a relatively large horizontal plane component and would tend to pass
through the weaker layers in the deposit. The direct use of the field vane strength along
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 108
horizontal failure surfaces would not be representative and as such is not the characteristic value
for horizontal failure surfaces. In the author’s experience, a suitable characteristic value for
horizontal weak layers in varved clay deposits may be the strength corresponding to 0.75 times
the measured vane strength.
Another example is the direct use of strength measured using small sized tested specimens taken
from a fissured clay or fractured bedrock. The measured value is often the intact strength, which
Downloaded from ascelibrary.org by University of California, San Diego on 06/02/17. Copyright ASCE. For personal use only; all rights reserved.
is much higher than the overall operational or mobilized strength of the ground mass that is
largely controlled by the fissures and fractures (Lo 1970).
Factors such as fabric anisotropy and stress anisotropy as described and discussed in the
technical literature (e.g., Bjerrum 1972, Becker et al. 1984) need to be considered when
assessing appropriate values for the characteristic value. All factors that influence the
geotechnical parameter or property must be taken into account. When assessing characteristic
values for strength, the key consideration by the geotechnical engineer is suitably assessing the
mobilized strength for the given limit state under consideration.
It is important to realize that the above considerations also apply when selecting suitable design
values for use in traditional allowable (working) stress design.
What is meant by “within the volume of ground (zone of influence) under consideration”?
Figure 6 illustrates the intended meaning of zone of influence. For a given site, there is more than
one characteristic value for strength and displacement (deformation) parameters. The
characteristic value is not simply the mean (or cautious estimate of the mean) based on all the
test results from the site investigation program. For example, for the design of spread footings,
the test results within the stress bulb (zone of influence) would be considered in assessing the
characteristic value. For deep foundation design it would be along and deeper than the
anticipated pile length. If the ground is layered there would be characteristic values for each of
the individual layers. For embankments, stress anisotropy considerations and parameters
associated with the anticipated field stress range for the loading condition need to be taken into
account. This similarly applies to excavations and cuts, except that parameters associated with
unloading condition and anticipated stress range would be considered when selecting appropriate
values for characteristic values. The primary reason for these considerations is that ground
behaviour exhibits significant non-linearity.
The geotechnical engineer also needs to be cognizant of the interrelationship between resistance
and load factors and characteristic value when selecting characteristic geotechnical parameters
for design purposes. Currently, it appears that practicing geotechnical engineers may not be
aware of the importance of selecting the characteristic value in a specific manner as described in
the following paragraph.
Load and resistance factors and characteristic values are interrelated. Load and resistance factors
have been derived (calibrated) based on characteristic values that have been defined in a specific
manner. For consistent and rational design in practice, the selection of a characteristic value for
geotechnical resistance needs to be made in the same manner as that used to derive the resistance
factor. The geotechnical resistance factors were derived using characteristic values that are close
to the mean value; therefore, a characteristic value of a geotechnical property that is similarly
close to the mean should be used in the calculation of geotechnical resistance. A value closer to
the lower bound value should not be used because the uncertainty associated with a given design
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 109
parameter has already been incorporated into the numerical values of the resistance factors.
Should a value closer to the lower bound strength be chosen, an additional level of conservatism
enters into the design. Although an additional level of reliability and safety may be viewed as
desirable, it could have significant cost implications as discussed later in this paper.
Downloaded from ascelibrary.org by University of California, San Diego on 06/02/17. Copyright ASCE. For personal use only; all rights reserved.
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 110
As an example, a new highway was built with its alignment extending several tens of kilometres
across undeveloped terrain involving extensive deposits of glacio-lacustrine and glacio-fluvial
sediments. The site investigation program consisted of sampled boreholes with frequent field
vane tests and was complemented by a reasonably detailed laboratory testing program. The site
investigation program was considered consistent with good state-of-practice. The highway
embankment performed well except for a small reach where a series of embankment foundation
failures took place. In the subsequent forensic investigation it was found that the data along the
alignment was compiled and plotted collectively, including plots of field vane test results. The
plasticity of the clay deposits was such that adjustment to measured vane strength (e.g. Bjerrum
or Aas correction factors) was not required.
Although the details are not certain, it appears that all test results were plotted using the same
symbol for each borehole rather than specific and unique symbols to denote each borehole along
the alignment. The range in the collective strength vs. depth profiles did not raise flags of
significant variability. Figure 7a shows the design strength profile.
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 111
The profile corresponds to a trend line slightly less than the mean and as such could be viewed as
being consistent with a ‘cautious estimate’ of the mean. The author has provided the
corresponding probability density function plotted at a given elevation within the strength vs.
depth profile. The forensic investigation into the embankment failures showed that the failed
section of the alignment contained three boreholes in which the measured field vane strength
comprised the lower bound of the plotted strength profile as shown by the symbol “x” in Figure
7b. The figure also shows that for the project, there should have been at least two separate
Downloaded from ascelibrary.org by University of California, San Diego on 06/02/17. Copyright ASCE. For personal use only; all rights reserved.
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 112
strength parameters should be grouped into separate loading and unloading probability density
functions – not lumped all together.
Lesson Learned: Parochial knowledge and long standing local practice (i.e. “we always do it
this way”) can stand in the way of good engineering and improvement in terms of economic
benefit and efficiency. The application of statistics and probability concepts is not to replace
good understanding of fundamental soil behaviour and theoretical considerations. Statistics and
Downloaded from ascelibrary.org by University of California, San Diego on 06/02/17. Copyright ASCE. For personal use only; all rights reserved.
probability are useful tools if properly applied. They must not become a substitute for
understanding fundamental behaviour, loading conditions and failure mechanisms, geology and
geological depositional processes.
RELIABILITY BASED DESIGN CONCEPTS APPLIED TO SLOPE STABILITY
The CHBDC (2014) specifies geotechnical resistance factors (Ф) for slope stability. For the
temporary condition, Ф ranges from 0.7 to 0.8, depending on degree of understanding of
geotechnical conditions at a given site. For the permanent condition, the corresponding values
for Ф range from 0.6 to 0.7 (Fenton et al. 2016). A potential implementation issue is that the
result from geotechnical limit equilibrium analysis is expressed as a global factor of safety, FS.
Thus, it is necessary to relate FS and geotechnical resistance factor. Becker (1996) developed,
through calibration by fitting with allowable (working) stress design, the following approximate
relationship between global factor of safety (FS) and load and geotechnical resistance factors:
FS = α/Ф [11]
In geotechnical limit equilibrium analysis, α = 1.0 (i.e., no load factor applied to unit weight) and
thus FS = 1/Ф. For the range of Ф from 0.6 to 0.8, the corresponding calculated value for FS
ranges from 1.67 to 1.25.
Now let’s examine the above relationship using reliability based design concepts. Meyerhof
(1970 and 1995) compared FS, probability of failure (pf), reliability index (β) and coefficient of
variation (Vr). For earthworks stability, the corresponding value of pf is in the range of 0.005 to
0.001 and β is in the range of 2.3 to 3.0 when FS = 1.5 and Vr = 0.15.
Table 1 summarizes the relationship between FS and Vr for β in the range of 2.3 to 3.0. The
range in FS in Table 1 is generally consistent with the range of 1.25 to 1.67 using the relationship
that FS is approximately 1/Ф (when α = 1.0), obtained by calibration by fitting with allowable
(working) stress design. This agreement provides a degree of confidence that FS = α/Ф is a
reasonable engineering approximation.
Table 1. Summary of FS and Vr values for reliability index in the range of 2.3 to 3.
Vr FS
0.10 1.3
0.15 1.5
0.20 1.8
Using the FOSM Equation [9] and substituting the values of kr ranging from 1.0 to 1.1,
= 0.75, β ranging from 2.3 to 3.0 and Vr ranging from 0.1 to 0.2, results in Ф ranging from
0.64 to 0.93, which is generally consistent with the range in Ф values cited in CHBDC (2014).
However, if the value of Vr = 0.3 is used, Ф ranging from 0.51 to 0.66 is obtained, and the
corresponding values for FS range from 2.0 to 1.5. Therefore, the geotechnical engineer should
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 113
be cautious when there are data to establish the Vr of the strength parameters and the value is
approaching 0.3 or higher. In this case the use of a lower Ф or higher FS may be warranted when
the characteristic value used for analysis and design corresponds to a bias factor in the range of
1.0 to 1.1.
The CHBDC (2014) introduces Consequence Factor, Ψ as shown in Equation [2]. The equation
for FS then becomes 1/( Ψ Ф). So for a structure with high consequence ( Ψ = 0.9) and typical
Downloaded from ascelibrary.org by University of California, San Diego on 06/02/17. Copyright ASCE. For personal use only; all rights reserved.
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 114
Figure 8a. Histogram and normal distribution – peak effective friction angle (c’ = 0 kPa).
Figure 8b. Histogram and normal distribution – residual effective friction angle
(c’ = 0 kPa).
From Equation [9], Vr is the controlling parameter for fixed values for β, kr and θ. Table 3
summarizes the results for the cases of kr = 1.0, = 0.75 and β = 2.3 and 3.0, and Vr ranges by a
factor of 2 (i.e., 0.30 to 0.15 and 0.2 to 0.1). When Vr is changed from 0.2 to 0.1 there is an
approximate change in the value of FS of 0.25 or about 20%. Similarly, when Vr changes from
0.30 to 0.15, there is an approximate change in the value of FS of 0.45 or about 35%. Therefore,
if the basis of design is 1.3 for the most likely case or peak strength parameters, the use of 1.1 for
residual strength or reasonable worst case parameters appears to be appropriate on the basis that
the Vr may be a factor of about 2 lower than for peak strength parameters.
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 115
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 116
In the author’s experience and opinion, a primary cause of increased foundation costs is the
inappropriate selection of characteristic value as discussed in earlier sections of this paper. In
summary, most practitioners tend to stay with the procedure and basis of selecting design values
such as strength as they did when using allowable (working) stress design. In many cases they
would tend to select characteristic strength values (profiles) closer to lower bounds than the
mean value or a cautious estimate of the mean. When a lower bound strength is used, an
additional level of conservatism enters into the design, which tends to increase costs. As
Downloaded from ascelibrary.org by University of California, San Diego on 06/02/17. Copyright ASCE. For personal use only; all rights reserved.
discussed previously in this paper, the geotechnical resistance factors listed in the NBCC and
CHBDC are based on characteristic values that are close to the mean.
The owner’s experience is also that most practitioners say that they must follow the code and use
the geotechnical resistance factors provided in the code. The NBCC does not include the
geotechnical resistance factors within the code (legally binding) portion of the document. The
factors are provided in the Commentary to the code (not legally binding) with reference to the
Canadian Foundation Engineering Manual (CFEM 2006). The language (provisions) in the
NBCC and CHBDC allows the engineer the freedom to be innovative. The codes permit the
engineer to develop project and site specific values of reliability index and load and resistance
factors, based on sound engineering principles and fundamental theoretical considerations.
However, the development of project and site specific criteria requires a reasonably significant
effort in terms of resources (schedule and cost). In many cases, there is not enough economic
incentive to do this because there is insufficient cost savings to justify the effort. For large
projects, however, substantial cost savings can be realized, thereby justifying the investment.
While staying within the intent of codes, there is an opportunity to increase the values of
geotechnical resistance factors provided that: (i) site-specific geotechnical uncertainties are
characterized; (ii) reliability of the geotechnical prediction method is quantified; and (iii) the
target reliability index is rationally selected and meets the intent of the codes. The selection of a
suitable project specific target reliability index is ultimately the responsibility of the owner,
though the consultant’s role and contributions are important.
The following case record shows the key elements of how project and site specific values can be
developed within the framework of reliability based design and within the intent of the NBCC.
The details of the development of project and site specific reliability index and geotechnical
resistance factors are provided in Thomson et al. (2016) and (2017). The following discussion
provides a summary of that work.
Bitumen Processing Facility
A proposed bitumen processing facility in Northern Alberta includes many low to no occupancy
structures that support pipe racks and processing units, and provide protection to pumps and
other mechanical/electrical components. Several thousand steel piles will be driven during
construction of the facility. The owner put the challenge to the geotechnical consultant to provide
foundation design recommendations that comply with code requirements and result in
constructed foundations that are cost effective to the fullest extent practicable.
Values for reliability index and associated geotechnical resistance factors for design of structures
are given in applicable codes (e.g., NBCC). However, the use of the building code may not be
strictly applicable to structures with limited to no occupancy that do not satisfy the definition of
building or do not have high failure consequence should they not perform as expected. The
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 117
reliability for such structures can be lower than that normally used in building codes. The direct
use of geotechnical resistance factors identified in a building code in such circumstances would
tend to produce a foundation design that is too conservative and expensive.
The ground conditions at the site comprised sandy silty clay till containing dense to very dense
sand and silty sand interlayers. The till was cohesive with intermediate plasticity, varied from
firm to hard and contained varying proportions of sand and gravel. The geotechnical
Downloaded from ascelibrary.org by University of California, San Diego on 06/02/17. Copyright ASCE. For personal use only; all rights reserved.
investigation included 23 sampled boreholes, 41 cone penetration tests with pore pressure
measurement (CPTs), and a laboratory index, shear strength and stress-deformation testing
program. Groundwater monitoring instrumentation was also installed.
In addition to the reasonably comprehensive site investigation program, a static compression pile
load test program consisting of five tests was carried out in conformance with ASTM D1143 as
part of design to develop site-specific geotechnical resistance factors. The test and reaction piles
had diameter of 324 mm and 406 mm with wall thicknesses of 9.5 mm and 12.7 mm,
respectively. The piles were driven open-ended using a hydraulic hammer. Pile embedment
lengths were pre-selected to vary between 13 m and 18 m. High strain dynamic testing was
completed at end of initial drive and beginning of restrike on all test and reaction piles using a
Pile Driving Analyzer (PDA) in conformance with ASTM D4945, followed by Case Pile Wave
Analysis Program (CAPWAP) assessment of selected hammer blows.
Inherent (spatial) uncertainty was quantified by predicting the geotechnical axial compression
resistance of piles with various diameters and lengths at each individual boring and CPT
investigation location within the site. The geotechnical resistance was predicted with methods
commonly used in Canadian practice (e.g., unit shaft resistance estimated as a proportion of soil
undrained shear (CFEM 2006)). A typical frequency distribution of predicted geotechnical
resistance is shown in Figure 9 for a 406 mm diameter pile with a 15 m embedment length.
Comparison between predicted geotechnical axial compression resistance and measured axial
compression resistance showed a bias factor (measured/predicted) ranging from 0.80 to 1.64,
with a mean value of 1.16 and coefficient of variation of 0.18.
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 118
Assuming that the sources of uncertainty are statistically independent, the overall coefficient of
variation (V ) was calculated as (FHWA 2001):
V = V +V +V [12]
where Vinh, Vmeas and Vpred are the coefficients of variation for inherent (geological/spatial),
measurement and prediction model uncertainty, respectively. The results of the work showed that
Downloaded from ascelibrary.org by University of California, San Diego on 06/02/17. Copyright ASCE. For personal use only; all rights reserved.
V was 0.16 (see Table 4), V was 0 (assumed) and Vpred was 0.18. The resulting value of
Vr was 0.24.
Table 4. Inherent (spatial) uncertainty parameters (from Thomson et al. 2016).
Predicted Mean Standard
Pile Diameter Embedment Length
Resistance Deviation
(mm) (m)
(kN) (kN)
13 2,190 378 0.17
406 15 2,520 448 0.18
17 2,722 406 0.15
13 1,689 261 0.16
324 15 1,948 311 0.16
17 2,122 295 0.14
After geotechnical uncertainty had been quantified, appropriate target reliability indices were
developed by the consultant in consultation with the owner, with consideration of project risks
and redundancy in pile groups. The owner selected values of reliability index, β, ranging from
2.0 to 3.5, depending on the importance category of the structure and pile redundancy
considerations.
For the calculation of geotechnical resistance factors, the FOSM Equation [9] was used for the
selected β values, with kr (bias factor) = 1.16, Vr = 0.24, and = 0.75. The resulting geotechnical
resistance factors vary from 0.62 to 0.81 as shown in Table 5.
NBCC 2010 suggests a geotechnical resistance factor of 0.6 for analysis of deep foundations
using static load test results. The site specific geotechnical resistance factors in Table 5 are
greater than the values suggested in NBCC 2010 as a result of the site specific assessment of
geotechnical uncertainty, combined with project specific decisions on appropriate target
reliability indices. The geotechnical resistance factors in Table 5 are specific to the site and are
also tied to the geotechnical prediction method that was used, which is incorporated into
resistance bias factor using the high strain dynamic and static pile load test results. Different
prediction methods, or different methods of selecting the characteristic values, would likely
result in different geotechnical resistance factors being derived.
It was estimated by the owner that applying the results of this study would decrease the
foundation construction costs for this project by 40% and save millions of dollars.
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 119
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 120
experience.
ACKNOWLEDGEMENTS
The author acknowledges the support of Golder Associates Ltd. and the contributions of his
colleagues who worked on the projects and case records mentioned and described in this paper,
in particular Drs. Peter Thomson, Paul Dittrich and Graeme Skinner who reviewed the paper and
provided valuable and insightful comments and suggestions. Special thanks are also extended to
Ms. Sarah Bungay, Dr. Masoumeh Saiyar and Mr. Bob McDonald who assisted in the
preparation of this paper. The thoughtful review comments of Dr. Gordon Fenton are also greatly
appreciated and valued.
REFERENCES
Bathurst, R.J., Huang, B. and Allen, T.M. 2012. LRFD calibration of the ultimate pullout limit
state for geogrid reinforced soil retaining walls. ASCE International Journal of
Geomechanics, Special Issue on Geosynthetics 12(4): 399-413.
Becker, D.E. 1996a. The Eighteenth Canadian Geotechnical Colloquium: Limit States Design for
Foundations I. An overview of the foundation design process, Canadian Geotechnical
Journal, Vol. 33, No. 6, 956-983.
Becker, D.E. 1996b. The Eighteenth Canadian Geotechnical Colloquium: Limit States Design
for Foundations II. Development for the National Building Code of Canada, Canadian
Geotechnical Journal, Vol. 33, No. 6, 984-1007.
Becker, D.E. 2006. Limit States Design Based Codes for Geotechnical Aspects of Foundations in
Canada. International Symposium on New Generation Design Codes for Geotechnical
Engineering Practice, Taipei, Taiwan. November 2 and 3, 2006.
Becker, D.E., Crooks, J.H. and Been, K. 1988. Interpretation of the Field Vane Test in Terms of
In-Situ and Yield Stresses in Vane Shear Strength Testing in Soils, Field and Laboratory
Studies, Richards, A.F, Editor, ASTM. pp. 71-87.
Becker, D.E., Crooks, J.H.A., Jefferies, M.G. and McKenzie, K.J. 1984. Yield behaviour and
consolidation, Part 2: strength gain. Proceedings ASCE Geotechnical Engineering Division
Symposium on Sedimentation and consolidation models: predictions and validation, Ed.
Young & Townsend, pp. 382-398.
Bjerrum, L. 1972. Embankment on soft ground. Proceedings of the ASCE Specialty Conference
on Performance of Earth and Earth-Supported Structures. Purdue University, USA, 2: 1.54.
Canadian Geotechnical Society. 2006. Canadian Foundation Engineering Manual (CFEM) 4th
edition. Edited by D. Becker and I. Moore. BiTech Publishers, Vancouver, BC.
© ASCE
Geo-Risk 2017
Geo-Risk 2017 GSP 282 121
CHBDC – Canadian Standards Association. 2000. Canadian Highway Bridge Design Code – A
National Standard of Canada. CAN/CSA Standard S6-00.
CHBDC – Canadian Standards Association. 2014. Canadian Highway Bridge Design Code – A
National Standard of Canada. CAN/CSA Standard S6-14.
Federal Highway Administration (FHWA). 2001. Load and Resistance Factor Design (LRFD)
for Highway Bridge Substructures: Reference Manual and Participant Workbook.
Downloaded from ascelibrary.org by University of California, San Diego on 06/02/17. Copyright ASCE. For personal use only; all rights reserved.
© ASCE
Geo-Risk 2017