363 - Reliability Based Calibration of Foundation Strength Factor Using Full-Scale Test Data - A Guide For Design Engineers
363 - Reliability Based Calibration of Foundation Strength Factor Using Full-Scale Test Data - A Guide For Design Engineers
Working Group
B2.07
December 2008
WG B2.07
WG B2.07 members
Regular Members
Mike Vanner (UK), Anthony Digioia (USA), Branko Zadnik (Slovenia), Jean-Pierrre Chouteau
(Canada), Asim Haldar (Canada), Eugene Dembicki (Poland), Aureo Ruffier (Brazil), Sylvain
Muyard (France), Mike Buckley (Ireland), Maurizio Leva (Italy), Michael Pietschke (Germany), John
P. Sivertsen (Norway), Boud Boumecid (UK), Toshiyuki Ozaki(Japan)
Corresponding Members
Kari Sahla (Finland), Graeme Paterson (Australia), N Sabri (Switzerland), Willem Combrinck (South
Africa), Patrick Eychenne (France)
Copyright©2008
“Ownership of a CIGRE publication, whether in paper form or on electronic support only infers right
of use for personal purposes. Are prohibited, except if explicitly agreed by CIGRE, total or partial
reproduction of the publication for use other than personal and transfer to a third party; hence
circulation on any intranet or other company network is forbidden”.
Disclaimer notice
“CIGRE gives no warranty or assurance about the contents of this publication, nor does it accept any
responsibility, as to the accuracy or exhaustiveness of the information. All implied warranties and
conditions are excluded to the maximum extent permitted by law”.
3.7.2 Site specific test data along the line route..................................... 3-14
3.8 Summary ................................................................................................3-15
4.0 CALIBRATION OF FOUNDATION STRENGTH FACTOR ( m - MODEL
OF BUCKLEY)..................................................................................................... 4-1
4.1 General .....................................................................................................4-1
4.2 m - Model – No Correlation (Repeated Test Data)................................4-1
4.3 m - Model-With Correlation (Scattered Test Data)................................4-3
4.3.1 Discussion ....................................................................................... 4-5
4.4 Example Problem- ESB Data (Buckley, 1994).......................................4-6
4.4.1 Background ..................................................................................... 4-6
4.4.2 Defining foundation capacities ....................................................... 4-7
9.5 Correlation check (Test No. 1- equation 4.1 and replace y by R)........9-5
9.6 Regression Analysis (variance of residuals is constant) ....................9-5
9.7 “Zero Intercept” Check – [Test No. 2 - equation 5.11] .........................9-7
9.8 Variance of residuals-Test (Test No. 3- Figures 5.4a and 5.4b) ..........9-9
9.9 Predicted Test Capacity (sample size effect included) ........................9-9
9.10 Comparisons of foundation strength factor values (Case C and Case
D) .................................................................................................................9-13
9.11 CIGRE (2002)..........................................................................................9-14
SOMMAIRE
Key Words:
Mots clé:
The author would like to acknowledge the contributions of the following individuals during
the preparation of this report. The author would like to thank Electricity Supply Board
(ESB), Ireland for providing the uplift test data on undercut pad and chimney foundation.
The author would also like to thank Dr. Jan Rogier from Belgium and Dr. Ruy Menezes
from Brazil for reviewing the report and providing valuable comments. Finally, the author
sincerely appreciates the long term support provided by the Senior Leadership Team of
Newfoundland and Labrador Hydro.
1.0 Introduction
Overhead line support structures are subjected to the effects of meteorological loads as
well as to the longitudinal loads due to broken conductor and unbalanced ice shedding.
Forces transmitted at the foundation level are often uplift (pullout) forces with small
transverse shear forces for tangent lattice towers or large overturning moments with
corresponding shear forces for single pole, angle or dead end structures. Foundation
types are normally spread footings and steel grillages for lattice towers and drilled shafts
for single poles, angle or dead end structures. Figure 1 depicts a typical undercut pad and
chimney foundation under uplift load.
Backfill
G.L.
Native
W soil Soil
Frustum
Angle, Assumed
Failure
Surface
Since the meteorological loads are probabilistic in nature, the transmitted forces acting on
the foundation are also probabilistic and are subjected to a high level of uncertainties. In
addition, spatial variation of soil strength parameters also contributes a high level of
uncertainty. Therefore, the foundation design must consider all these uncertainties
quantitatively.
2.0 Objective
The primary objective of this study is to present method(s) of undertaking the calibration of
the theoretical geotechnical design model and hence determine the characteristic capacity
of the foundation, based on full-scale test data. To develop this calibration methodology, a
framework to determine the prediction interval of the foundation test capacity based on a
limited number of full-scale test data is developed using a linear regression model. The
prediction interval provides a confidence on a new test capacity value predicted given a
single geotechnical design capacity value, Rn . The lower bound of the prediction interval is
directly linked to the foundation strength factor, φF which can be used in reliability based
design equation for a single foundation component. The foundation strength factor with an
exclusion limit is defined as the ratio of the lower limit of the prediction interval (predicted
test capacity) with a certain confidence level (analogous to characteristic capacity) and the
3.0 Methodology
To develop a reliability based calibration methodology for foundations, the four steps
identified here are followed:
• Data analysis and calibration based on m - model (Buckley, 1994) and linear
regression model
IEC 60826 provides a semi probabilistic design equation in terms of the limit load acting on
a component and the strength of the component. The design equation (1) relates the
characteristic capacity, RC to the appropriate load effects, QT having a specified return
period of T-years with a global strength factor, φ R .
φ R RC ≥ QT (1)
Area, α
L
s R = Standard Deviation
Rc Rn
Figure 2 depicts the characteristic capacity based on e% exclusion limit for a normal
distribution. For a 10% exclusion limit the characteristic capacity, RC ,10% is
where R and vR are the mean and the coefficient of variation (ratio of the standard
deviation, s R and the mean value, R ) of the test capacity values for a finite sample size,
n. The term kα is a measure of how far the point “L” is from the mean to obtain the 10%
exclusion limit ( kα =1.28). For a finite sample size, n < 30, kα should be corrected based
on the Student t-distribution. The exclusion limit 10% implies that the characteristic
capacity is for the 90% confidence level and that there is a 10% probability that the
capacity could fall below this value.
Within the framework of IEC 60826 (2003), the characteristic capacity RC should be
related to the geotechnical design capacity Rn and this could be done through calibration.
If the full-scale test data is available, the appropriate foundation strength factor φ F
RC
φF = (3)
Rn
can be used to relate the geotechnical design capacity Rn to the characteristic capacity
RC with a defined exclusion limit (5% for ASCE and 10% for IEC 60826). The
geotechnical design capacity Rn in figure 2 is derived from analytical and/or semi
empirical relationship or numerical models. Depending on the calibration assumptions, the
relationship could be linear or non linear.
Once the geotechnical design capacity Rn has been determined, the foundation designer
is given guidance on how to undertake a full-scale foundation test programme. Several
steps need to be considered before a full-scale foundation test programme can be
commenced. These include: (a) selection of the test site(s), (b) selecting the applicable
foundation geotechnical design model(s) to determine the foundation’s geotechnical
400
300
200
100
0
0 5 10 15 20 25 30 35
Displacement (mm)
3.4 Data analysis and calibration based on m - model and linear regression
model
3.4.1 General
∑ ( xi − x )( yi − y )
r= (4)
(n − 1) s x s y
where x and y are the mean values and s x and s y are the standard deviations of the
geotechnical design capacity values xi and the test capacity values yi respectively. A
low correlation value ( r 2 ) would imply that there is no linear relationship between the
Linear association
Correlation Value, r 2
between x and y
weak
Low , r 2 ≤ 0.3
moderate
Moderate (0.3 ≤ r 2 ≤ 0.8)
strong
High, r 2 > 0.8
Test Capacity, y
Data Banded
Around The
Line With
High sxy
Data Banded
Around The
Line With
Low s xy
3.4.2 m - model - no correlation (repeated test data at one specific test site)
Since xi = x , for the repeated test data at one specific site, the correlation coefficient r=0.0
(equation 4) indicates that the data is unique to a specific test site. In this case, the
3.4.3 m - model - with correlation (scattered test data along a line route)
Once the user has determined the test and geotechnical design capacities for the various
test sites (s) and foundation sizes, the user can choose a simple calibration model based
Test Capacity
on m = ratio(s) following (Buckley, 1994). Using this
Geotechnical Design Capacity
approach, it can be shown that the foundation strength factor, φ F for a large sample size
(n >30) is
RC
φF = = m (1- kα * vm ) (5)
Rn
where
However, to use this model, the user needs to ensure that the correlation of the data is
high (Table 1) and the coefficient of variation of the model error, vm is low. Provided these
two conditions are met, it can be shown that the m -model is a subset of the general linear
regression model (figure 5). The specific conditions are: (1) zero intercept and (2) non
constant variance of residuals, where the residual error values vary linearly with
geotechnical design capacity values (x-values). The residual error value ε defines the
difference between the actual test capacity and the estimated mean test capacity value
from the linear regression model. For other situations where this is not the case, the
foundation designer should undertake a more detailed analysis based on a linear
regression model; either to ensure that the intercept effect can be neglected or
alternatively to improve the theoretical design model. If the intercept effect cannot be
neglected, then an alternate methodology based on a linear regression model is also
presented illustrating how to include this in the analysis and to estimate the predicted
capacity with the required confidence interval.
In the linear regression model, the designer determines the slope, intercept and the
variance of residuals. Figure 5 depicts this information. Figure 5 also presents the
confidence interval on the estimated mean test capacity value, as well as, the prediction
interval on a new test capacity value given a geotechnical design capacity value. This
figure is very similar to figure 2, except that in figure 2 the characteristic capacity is
computed for the test data only; while, in figure 5 this is determined from the linear
relationship between the test capacity and geotechnical design capacity values.
130
Prediction Interval
Estimated Mean Test Capacity,
120 y^ at x o
Residual
Error, ε
*
100
A
A simple statistical test is presented to check for zero intercept condition ("a" on figure 5)
with a slope, m (b on figure 5) and a graphical method to test for variance trend of
residuals (figure 6). For a non constant variance of residuals, an appropriate weight
1
( wi = ) can be used to make a correction to the original dataset (Ang and Tang, 1975).
xi 2
The non constant variance of residuals is defined as the residual error value, ε in figure 6b
increasing with the geotechnical design capacity value. Moreover, the effect of constant or
non constant variance of residuals is presented and its effect is emphasized in analyzing
the data.
C onstant V ariance
of R esiduals
ε) Residual Errors (
ε = 0.0
ε)
Shaped
Residual Errors (
The foundation strength factor, φ F under the linear regression model is defined as the
ratio between the lower limit of the prediction interval, y p (figure 5, predicted capacity)
and the geotechnical design capacity value, x0 used to determine this capacity.
yp
φF = (6)
x0
Since there are two options (constant versus non constant variance of residuals) under
each of these conditions (zero or non zero intercept condition), the user will encounter four
possible cases for analyzing the data. Table 2 can be used as a guide for calibrating the
foundation strength factor value for a specific foundation type. Each case in Table 2 will
provide a lower limit on the predicted capacity. The equation for each case is given in the
technical brochure. The test for zero intercept is definitive. However, the test for constant
versus non constant variance of residuals is not quantitative and is rather subjective;
because the designer needs to plot the residual error values against the geotechnical
design capacity values and derive a trend line similar to figures 6. A flow diagram is
developed to assist the foundation engineer in choosing a particular case.
Variance of Residuals
Intercept
Constant Non constant
Non Zero Case A Case B
Zero Case C Case D
yi
Compute mi =
xi
High or Moderate
Correlation
Slope &
Intercept Test 2: Slope only
Test
Test 2:
2:
Zero Intercept
Eq.
Eq. 4.14
Check
Test 3: Test 3:
Chk Figs. 6(a) Chk Figs. 6(a)
& 6(b) & 6(b)
Calculate Foundation
Strength Factors to Provide
Conservative Values -
Equation 6
Test 1: Test for Correlation END
__
Test 2: Simultaneous Test for Intercept = 0.0 and Slope = m
Test 3: Visual Comparison to Test Constant Variance of
Residuals vs. Non Constant Variance of Residuals
*Constant COV - Special Case
(Residual values banded around zero line or "fan
- CIGRE (2002)
shaped" - use judgment)
4.1 m -model
Electricity Supply Board (ESB), Ireland provided data for pad and chimney foundations
tested under various soils conditions. The nine uplift foundation tests carried out in
cohesionless soils are initially analyzed using the m -model. The mean values of the test
and the geotechnical design capacities are 378 kN and 295 kN respectively, while the
coefficients of variation are approximately 35%. The m and vm values are 1.30 and 0.156
respectively.
Figure 8 presents the data for regression analysis. The correlation value r 2 is 0.81 which
is high (Table 1). The 300 frustum cone model in general, appears to be conservative
because of large bias, i. e. the model in general underestimates the test capacity value by
30%. A linear regression model (CASE A) is run first following the flow diagram in figure 7
and Table 2. The statistical test for intercept indicated that the zero intercept assumption
and the slope, m =1.30 ( m -model) can be accepted. Therefore, the designer has to
choose between Case C and Case D respectively (figure 7).
600
Test Capacity (kN)
400
Least Square Fit
200
0
0.0 200.0 400.0 600.0
Geotechnical Design Capacity (kN)
Next, figure 9 presents the residual plot. Two lines are drawn to see a trend line that the
variance of residuals is non constant (similar to figure 6b) and it is not fully definitive that
this is the case. Therefore, the designer should try both cases C and D for determining the
foundation strength factor values. Figure 10 depicts the comparison of foundation strength
-100
-150
Geotechnical Design Capacity (kN)
1.200
Foundation Strength Factor
1.000
0.800
0.600
Data range- 150kN Foundation strength
0.400
to 499kN factor(Constant),φF = 1.00
Case D
Foundation strength
0.200
factor(Variable),φF -Case C
0.000
0 100 200 300 400 500 600 700 800 900
Geotechnical Design Capacity (kN)
This study shows clearly that there are a number of steps one needs to consider before a
foundation strength factor, for a particular foundation type, can be estimated based on a
limited number of full-scale test capacity data and the theoretical foundation model
predicted geotechnical design capacity. The foundation engineer also needs to understand
that each dataset is different because of various test conditions, as well as model
parameters imposed on them. The WG B2-07 recommends that a clear and consistent
methodology based on accepted statistical theories, supported by valid assumptions, will
provide a good basis for testing any database and its practical use for reliability based
design of foundations. Simple tests, as proposed in this document include: correlation
test, simultaneous test for zero intercept and slope, trend line for variance of residuals
variation etc.
Based on these tests, the foundation designer can clearly decide whether to use the m -
model (CIGRE, 2002) or to use the linear regression model (alternative method),
applicable to those databases where the results clearly show that the assumptions made
under the m -model may not be valid. A flow diagram is developed to assist the foundation
engineer in choosing a particular option.
• The Buckley method relies on two very specific conditions which imply that the
intercept should be zero (figure 5) and the residual error values should vary with
the geotechnical design capacity values in an increasing manner (figure 6b).
• To apply this CIGRE method to a wide range of reliability based foundation design
calibrations with various degrees of correlation, the designer must be aware that
unless these specific conditions are met, the prediction interval with a high degree
of confidence (90% for IEC 60826 and 95% for ASCE) will produce results that
may not be valid.
Finally, it is shown that the general linear regression model with the unique test for “zero
intercept” is quite applicable to the analysis of a wide range of foundation test data with
6.0 References
Ang, A. H. S and Tang, W.H. 1975 Probability Concepts in Engineering Planning and
Design, Volume 1, John Wiley & Sons
E*
Figure 2.4 - Sample size versus ( ) for various confidence levels...............................2-8
sR
Figure 3.1 - Typical line route showing the test sites (Topographic Map).........................3-2
Figure 3.2 - Typical undercut pad and chimney foundation ..............................................3-3
Figure 3.3 - Failure in native soils .....................................................................................3-4
Figure 3.4 - Failure in backfill soils....................................................................................3-5
Figure 3.5 - Shear model (IEEE, 1985).............................................................................3-6
Figure 3.6 - Typical subsurface profile of a specific test site.............................................3-7
Figure 3.7 - Subsurface soil profiles with assumed failure surface ...................................3-7
Figure 3.8 - Grillage foundation under uplift load (Prasad & Haldar, 2001) ....................3-10
Figure 3.9 - Typical load displacement plot for a pad and chimney foundation ..............3-11
Figure 3.10 - Repeated test data for an undercut pad and chimney foundation-project
specific ....................................................................................................................3-14
Figure 3.11 - Uplift test -scattered test data along the line route (damage limit capacity at
10mm displacement) ...............................................................................................3-15
Table 5.7 –Foundation Strength Factor Values (Linear Regression Model) ...... 5-20
Table 5.8 - Typical Dimensions (m) of a Pad and Chimney Foundation (Figure
3.2) ............................................................................................................. 5-21
Table 9.1: Methodology to Follow the Flow Diagram ........................................... 9-3
Foundation - a structure below the ground to transfer support loads to the soils.
Null hypothesis -a proposition which is believed to be true; this could be that the line
passes through the origin while the alternative hypothesis could be the line does not pass
through the origin.
Significance level - a probability measure to reject a null hypothesis even though the
hypothesis is in fact true.
Exclusion limit - a variable defined in terms of percentage (%) to define the low side of
the mean value from a probability distribution.
Confidence interval - provides a probability measure of how close the population mean is
to the sample mean given a finite sample size; population mean is normally unknown and
is an unobservable quantity.
Prediction interval - provides a probability measure that a new test value drawn randomly
from a distribution will lie within a certain interval given a finite sample size.
Reliability index - a measure of structural reliability taking into account the variations in
load and strength.
Undercut foundation – a type of foundation where the soil above the base of the
foundation is undercut, such the foundation is in direct contact with the soil; very common
in European countries.
Frustum angle - an angle used to define the failure surface in the frustum model for
determining the uplift capacity.
Limit load - climatic load corresponding to a return period, T-years without imposing any
load factor.
Damage limit (working limit) - the load at which some permanent (irreversible)
displacement is observed and is similar to serviceability limit state in building codes.
Ultimate capacity - the load at which the foundation has undergone a significant
movement and a complete failure will happen if this load is exceeded.
Residual - difference between the actual test capacity value and the estimated mean test
capacity value from the linear regression line.
Constant variance of residuals - square of the residuals provide a constant value across
all values along x-axis (geotechnical design capacity).
Non constant variance of residuals - square of the residuals vary with values along x-
axis (geotechnical design capacity).
Foundation Strength factor - provides a measure of the relationship (linear or non linear)
between the characteristic capacity and the geotechnical design capacity using full-scale
test data.
Weighted residual - a term used in regression analysis to correct the data when the
variance of residuals is non constant.
1.1 General
Foundations for overhead line supports are frequently subjected to the effects of
various types of meteorological loads, i.e. due to wind, ice and combined wind and
ice loads on conductors and the support. Apart from these loads, supports are also
subjected to broken conductor loads as well as longitudinal unbalanced loads due to
ice shedding. The support foundation design must take into account the effects of all
these forces.
V1
T1
L1
T2 T2
L2 L2
V2 V2
T3 T3
L3 L3
V3 V3
T4 T4
L4 L4
V4 V4
T5
T6
T7
Forces transmitted at the foundation level are often uplift (pullout) forces with small
transverse shear forces for tangent lattice towers or large overturning moments with
corresponding shear forces for single pole, angle or dead end structures. Foundation
types are normally spread footings and steel grillages for lattice towers and drilled
shafts for single poles, angle or dead end structures. Figure 1.1 depicts a typical
double circuit tower with applied forces at the conductor attachment points. Figure
1.2 presents the resultant loads on a typical foundation, while figure 1.3 presents the
various resistive forces acting on this foundation.
Since the variable forces acting on the conductor and tower are probabilistic in
nature, the transmitted forces acting on the foundation are also probabilistic and are
subjected to high levels of uncertainties. In addition, soil conditions and
characteristics such as density, shear resistance, water level, vary also significantly
from one site to the other and also contribute a high level of uncertainty. To resist the
tower forces, the foundation must be designed based on the geotechnical design
capacity which is normally obtained from a theoretical design model or an empirical
method.
L2V
L2H
L2
Applied
forces
Applied
forces
Foundation = Resisting
weight forces
1.2 Objective
• Section 1 presents the introduction and outlines the objective of the study.
• Section 6 discusses the general criteria for extrapolation and its limitation
within the context of a linear regression model. A design example for
foundation characteristic capacity (predicted) is presented to show the effects
of interpolation versus extrapolation.
2.1 General
n
R = (∑ Ri ) / n (2.1)
1
1 n
sR 2 = ∑ ( Ri − R )2
n −1 1
(2.2)
s
vR = R (2.3)
R
Normally the test capacity value (R), is a random variable and is assumed to follow a
distribution (normal, log normal etc.). The population mean resistance ( μ ), is defined
as the estimated mean test capacity value of this random variable with a standard
deviation ( σ ), which are normally unknown. The variance of residuals is defined as
the square of the standard deviation and is a measure of the dispersion around the
estimated mean test capacity value.
The mean of a sample drawn randomly from the population is a point estimate of the
population mean ( μ ), and will vary for different samples. Given that the sample
mean is unlikely to be equal to population mean, how reliable R is the measure of μ .
The confidence interval provides this measure. In theory, as the number of samples
approaches infinity ( n → ∞ ), there will be no difference between the sample
statistics and population statistics. However, a large sample size is not practical for
Figure 2.1 depicts a typical normal probability density function of R, as f (R) and is
described by the following equation:
1 (R−R ) 2
− [ ]
1 2 sR
f ( R) = e (2.4)
sR 2π
Figure 2.1 also presents the concept of two sided confidence interval (CI), confidence
level and the level of significance. The significance level is based on the area under
the shaded region. The total area under the curve is 1.0.
Probability Density
f (R)
Function of R-
α = Level of significance Equation (2.4)
Area = Confidence level = 1 - α
Figure 2.1 - Two sided confidence interval and significance level (Rao, 1982)
The confidence interval provides an estimate where the true value may lie within the
upper and lower limits. For this case the confidence level is the area under the
probability density curve between the upper confidence (U) and the lower confidence
(L) limits marked on figure 2.1. Once the confidence level (p), is stated, the
probability of true mean (population mean) lying in the confidence interval is p.
Therefore, there is a probability of (1-p) % that the mean value will not lie within the
confidence interval. This quantity is also defined as the significance level, α (Rao,
1982) and is given by equation (2.5)
α = 1-p (2.5)
Table 2.1 - kα / 2 & kα - Values for various confidence levels (Benjamin &
Cornell, 1970)
s
CI= R ± kα / 2 R (2.6)
n
The population standard deviation ( σ ), is not known with certainty when the sample
size is less than 30 (n < 30). In this case, the two sided confidence interval (CI) is
obtained from equation (2.7)
sR
CI= R ± tα / 2,ν (2.7)
n
Table 2.2 provides some typical tα / 2,ν values for two confidence levels, p= 90%
and 95% respectively. This is normally the case when the sample size, n is less than
30.
Table 2.2 - tα / 2,ν & tα ,ν Values for figures 2.1 and 2.2 (Benjamin &
Cornell, 1970)
tα / 2,ν tα ,ν
Degrees of
Freedom, ν 90% 95% 90% 95%
( α / 2 = 0.05) ( α / 2 = 0.025) ( α = 0.1) ( α = 0.05)
1 6.32 12.7 3.08 6.32
2 2.92 4.30 1.89 2.92
3 2.35 3.18 1.64 2.35
4 2.13 2.78 1.53 2.13
8 1.86 2.31 1.40 1.86
10 1.81 2.23 1.37 1.81
11 1.80 2.20 1.36 1.80
30 1.70 2.04 1.31 1.70
∞ 1.64 1.96 1.28 1.64
Considering the effect of sample size n = 5 , equation (2.7) is used with tα / 2,ν =
2.13 for 90% confidence level (table 2.2). The revised calculation shows that the
38.5
lower and upper values are 224 kN and 298 kN respectively. ( 261 ± 2.13 x =
5
261 ± 37 kN); however, it is noted that the error margin ( ± E*) has now increased by
33% (37 kN versus 28 kN) because of the sample size effect. For the case when the
sample size is small, the sample size and the confidence interval width, ± E* =
sR
± tα / 2,ν are interrelated because tα / 2,ν depends on ν which is also
n
dependent on n (table 2.2). This will be discussed further in section 2.5 in
determining the sample size for foundation tests.
Similar to confidence interval (CI), the foundation designer can also consider an
interval such that there is a p% probability that a new test capacity drawn randomly
from the population will occur within the interval. This is called the prediction interval.
The confidence interval is for a population parameter (unobservable quantity and
unknown), while the prediction interval is for a single test capacity value (future
1
PI= R ± kα / 2 sR [1 + ]0.5 (2.8)
n
The prediction interval (PI) for a new test capacity value, Rn +1 for a small sample
size (n less than 30) is
1
PI= R ± tα / 2,ν sR [1 + ]0.5 (2.9)
n
The prediction interval (PI) is wider than the confidence interval (CI in equation 2.6)
because Rn +1 is a variable quantity while the sample mean, ( R ), is a constant
value. For large sample size n, equation (2.8) can be approximated as
PI= R ± kα / 2 sR (2.10)
However, in certain applications, the prediction interval (PI) can be used to derive a
guaranteed minimum test capacity with a specified confidence level. If the designer is
only interested in the minimum value (lower limit of the prediction interval) then this
becomes R − kα sR (one sided). Equation (2.10) can be used to provide this capacity
considering the lower limit only and figure 2.2 presents this case as a one sided
interval and kα / 2 is replaced by kα (table 2.1).
For example, if we have foundation test data of a large sample size with mean, ( R ),
s
and a coefficient of variation, ( vR = R ), then the characteristic foundation capacity
R
with an exclusion limit, e % (similar to significance level, α ) can be estimated from
equation (2.11)
where kα is a measure from the mean value to the point L, on the x-axis where a
specified exclusion limit is sought. The exclusion limit e=10% implies that the mean
characteristic foundation capacity is guaranteed for 90% confidence level and there
is a 10% probability that the mean capacity could fall below this minimum value. This
implies that the area under the curve from left of point “L” is 0.10 and the area under
the curve to the right of point “L” is 0.90. The geotechnical design capacity, Rn in
figure 2.3 is calculated based on analytical and/or semi empirical relationship or
numerical models (Buckley, 1994, Digioia et al, 1993).
f (R)
α = Level of significance
sR
x
0 L
(Lower Confidence
Limit)
k α sR from the One Sided Confidence Interval
mean to define
the interval
width
Figure 2.2 - One sided confidence interval, p and significance level α (Rao, 1982)
Area, α
L
s R = Standard Deviation
Rc Rn
The characteristic capacity with a 10% exclusion limit for the five tests presented
earlier is obtained from equation (2.11).
This capacity is based on a kα value of 1.28 when the sample size effect is not
included. Since the number of test is less than 30 (n < 30), tα ,ν value is used to
correct for the small sample size effect. Using tα ,ν value = 1.53 from Table 2.2 (for
α = 0.10 and ν = n-1 = 4), the revised estimate of the characteristic capacity is:
Quite often, the question is asked “what sample size do I need to estimate the
sample mean test capacity value which will be accurate to within ± E* of the
population mean value”. ± E* defines the error margin in Figure 2.1.
In order to match the accuracy between the sample mean and the population mean
by, ± E* the required sample size is:
k s
n = ( α / 2 R )2 (2.12)
*
E
However, if the sample size, n becomes less than 30 then a correction needs to be
made based on Student t-distribution. This will be discussed in the next section.
Therefore, to estimate the sample size, n we need to know the confidence level
(Figure 2.1), error margin and the standard deviation. Often we need to estimate the
standard deviation from prior knowledge or from the sample data.
For example, the engineer wants to estimate how many tests need to be done in
order to get an accuracy of E* = ± 28 kN with a 90% confidence level (p=0.9) when
k s 1.64 x38.5 2
n = ( α / 2 R )2 = ( ) = 4.98 ≅ 5. This was the sample size given in
*
E 28
Section 2.2.
600
550
500
kα =1.28, 90%
450
confidence Level
400
Sample Size,n
350 90%
300 kα =2.33, 99% 95%
250 confidence Level 99%
200
150
100
50
0
0 0.1 0.2 0.3 0.4 0.5
(Error/Standard Deviation)
E*
Figure 2.4 - Sample size versus ( ) for various confidence levels
sR
However, if the engineer decreases the accuracy of the prediction by increasing the
E* value to E * = ± 40kN, the estimated sample size becomes:
k s 1.645 x38 2
n = ( α / 2 R )2 = ( ) = 2.44 ≅ 3
*
E 40
Therefore, the sample size, n starts decreasing as we increase the error margin.
tα / 2,ν s R 2
n=( ) (2.13)
E*
The problem arises in solving this equation (2.13) for n because tα / 2,ν is not known
a-priori; since tα / 2,ν is dependent on sample size through ν because ν = n-1.
However, a trial and error procedure can be used in solving equation (2.13).
Let us assume that tα / 2,ν initially is equal to kα / 2 which is 1.645 for 90%
confidence level. Equation (2.13) provides a value of n=5 which gives tα / 2,ν = 2.13
from Table 2.2. Using this value in equation (2.12), we find n=8.57 ≅ 9. Therefore,
the revised tα / 2,ν = 1.86 and we obtain n=8.54 in the second iteration. Therefore
we need at least n=9 tests to maintain the same level of error margin of ± 28 kN.
Figure 2.4 presents this when the standard deviation is known (n > 30) and the
sample size estimate is based on equation (2.12). The foundation engineer can also
estimate sample size for (n < 30) following the iterative scheme presented before.
Current practices for design of overhead line support foundations are primarily based
on deterministic assumptions where the load effects on the foundation and the
foundation capacity are both defined in terms of single values. The load effects and
the capacity could be either defined in terms of the geotechnical design capacity or
the ultimate capacity with appropriate factors of safety. In simple form, this is
expressed as:
Rn
≥ ( D n +WI n ) (2.14)
FS
or
ψ *Rn ≥ ( Dn + WI n ) (2.15)
where:
IEC 60826 (2003) provides a basic framework where the semi probabilistic design
equation is given in terms of the climatic load acting on a component and the
strength of the component. The basic design equation relates the characteristic
strength (capacity) value to the appropriate load effects having a specified return
period of T years. In simple form, this can be expressed as:
RC ≥ QT (2.16)
The design equation is based on the assumption that the distribution of load is
Gumbel type 1 and the strength follows a normal distribution. The design equation is
valid with a constant annual failure probability between (1/T to 1/2T) provided the
coefficient of variation of strength is 0.05 ≤ vR ≤ 0.20 and of the load effect is
0.20 ≤ vQ ≤ 0.50 respectively.
Since IEC 60826 recommends the use of annual failure probability of (1/T to 1/2T), it
is important that the coefficient of variation of strength be kept under 0.3 to meet the
above criteria. For foundation, this may not be always possible. The capacity is
normally based on the available test data with a typical exclusion limit of 10% with a
50 year return period normally assumed for the load. IEC 60826 also recommends
larger return periods for higher reliability levels.
φR RC ≥ QT (2.17)
Within the framework of IEC 60826, the characteristic capacity should be related to
the model predicted geotechnical design capacity and this could be done through
calibration. If the full-scale test data are available, the appropriate foundation strength
factor φF can be used to relate the geotechnical design capacity to the full-scale test
capacity with a defined exclusion limit (5% for ASCE, 10% for IEC 60826).
φF Rn = QT (2.18)
where φF is a function of the statistical properties of the test and the geotechnical
design capacities. The details of computing the calibration factor, φF based on m -
model and linear regression model will be discussed in sections 4.0 and 5.0
respectively.
Although the ASD approach has served the utility industry reasonably well,
increasingly utilities’ design practices are moving towards a reliability based design
methodology, in order to quantify the foundation reliability in a predictable manner. A
simple reliability based design format is the use of load and resistance factor design
(LRFD) approach where the geotechnical design capacity and the load effects are
factored appropriately to account for a specified target reliability index β . The
reliability index β provides a measure of reliability and is directly related to the
failure probability which can be determined from normal tables (Melchers, 1999).
ϕ Rn ≥ γ D D n +γ WI WI n (2.19)
where
ϕ = the capacity reduction factor or resistance factor and is always less than
1.0,
γ D = load factor for dead load Dn ,
γ WI = load factor for the combined wind and ice loads WI n .
If ice is not present, γ WI can be replaced as γ W , to represent the load factor for wind
load effect, Wn . These factors are calibrated with respect to existing design
practices, which are often deterministic in nature or against a design format
prescribed, by the standard or code. For a simple log normal distribution of load
effects and strength values, these factors can be computed analytically (Melchers,
1999).
2.7 Summary
This section provided a basis for determining (1) the confidence interval on mean
value and (2) the prediction interval on a future test value. The lower limit of the
prediction interval can be used to determine the guaranteed minimum foundation
capacity, with a specified confidence level; while, the characteristic capacity is
defined as the low side of the mean foundation capacity with a specified exclusion
limit. For a large sample size, both will provide almost an identical value. The
characteristic capacity is applicable to single variable data (test capacity only); while,
the concept of prediction interval based on a linear regression model can also be
applied to the problem where more than one variable is involved, (such as the test
3.1 General
This section of the report presents a brief review of various design models in
determining the geotechnical design capacity of a pad and chimney foundation under
uplift condition. The soil strength parameters are identified for capacity calculation
followed by the field data collection programme. A short description of the field test
programme is presented which outlines the various steps that need to be followed in
collecting the test data as well as the data interpretation and analysis using the
methodology described in IEC 61773.
The transmission line engineer can do these tests in three different ways. The
engineer can carry out a number of tests with a foundation of fixed dimensions, in
homogeneous soils, with fixed installation technique to ensure that all tests will
produce similar results in soils of similar categories. This would provide a good
coefficient of variation of a specific foundation type and can be compared with the
coefficient of variations of other line components such as the support, conductor, and
insulator to coordinate the strength as per IEC 60826.
Alternatively, the engineer can vary the dimensions and carry out these tests in
similar soils to assess the effect of the geometry and can estimate the coefficient of
variation of the capacity.
Finally, the engineer can test a typical foundation of fixed dimensions, in various soil
types with repeated tests, thereby obtaining a wide range of coefficient of variations,
which will be more representative of a typical transmission line route. Of course, the
last two approaches will cost the utility more compared to the test carried out in one
soil type.
If we assume that the engineer has designed a typical foundation and wants to
determine the actual or characteristic capacity with certain confidence limits, then it is
necessary to test the foundation for a wide range of soil conditions to provide a good
level of confidence in the predicted capacity. In testing a specific foundation type, the
engineer needs to ask a few questions and these are listed as follows.
Assuming that the specific foundation type is an undercut pad and chimney
foundation, and the foundation is to be tested under uplift load, the first thing the
designer needs to ask is “what design model is appropriate for predicting the
geotechnical design capacity of this type of foundation”. Figure 3.2 depicts a typical
undercut pad and chimney foundation.
An undercut foundation is where the base of the foundation is slightly wider than the
width at the neck level. The foundation is excavated by machine to the width ‘C’ in
figure 3.2. The bottom of the excavation is widened to the width ‘B’ in figure 3.2 using
a handheld shovel. The undercut is usually 0.1 m on each side of the foundation.
Table 3.1 presents the dimension of a typical pad and chimney foundation.
F G.L.
C
Soil
E
A
D
B Undercut
Type A B C D E F
Square 2.2 2.0 1.8 1.0 0.8 0.1
Based
The failure mode will often depend on the type of the installation used. For example,
if the foundation excavation is done in such a way that the excavated hole is sized
sufficiently to place the foundation vertically or cast in place vertically with the
surrounding space filled by compacted backfill, then the failure zone would most
likely be covered by both native as well as the backfill soils. Figure 3.3 depicts a
situation where the failure zone is in the native soils. However if the excavation is
“open pit type” and is sufficiently large then it is possible to have the failure only
within the backfill zone (Figure 3.4). In this case, the backfill soil properties would
determine the geotechnical design capacity. This is more appropriate for sites with
boulder, cobbles etc. where often the excavation cannot be done narrowly or where
the walls of the excavation cannot stand up sufficiently safely.
In this particular model, it is assumed that an inverted frustum starts at the base of
the undercut foundation and slopes at an angle, ψ to the vertical to the surface (refer
to figures 3.3 and 3.4) and the weight of the soils within this frustum and the
foundation weight itself resist the uplift force. Mathematically, this is described as
follows:
U = VS ρ S + VC ρC (3.1)
where VS and VC are the volume of the soils and the concrete within the frustum
while ρ S and ρC are the unit weight of the soils and the concrete respectively. U is
the uplift force.
Backfill
G.L.
Native
W soil Soil
Frustum
Angle, Assumed
Failure
Surface
Backfill G.L.
Wsoil
Native
Soil
Frustum Assumed
WF Failure Surface
Angle,
Within Backfill
In this model, it is assumed that the uplift force is resisted by shear force acting along
the vertical shear surface circumscribing the base of the foundation, the weight of the
concrete foundation, WC and the weight of the soil above the concrete base within
the frustum, W S . The uplift force is determined as:
U = Sτ + WS +WC (3.2)
Although a constant shear strength value of 19.6 kN/m2 was used by Killer (1953),
Buckley (1994) recommended varying shear strength values with depth to ensure
that the predicted ultimate capacity is more realistic and comparable with the test
capacity. Mors (1965) presented an empirical expression for the frictional resistance
(shear capacity):
F0 = SτAη (3.3)
FO WC WS FO
A
The strength parameters required for an undercut pad and chimney foundation under
uplift load are :( 1) soil density, ( ρ S ), (2) concrete unit weight, ( ρC ), and (3) shear
strength, ( τ ), along the depth depending on the model used.
Once the number of sites along the line route (Figure 3.1) has been chosen, soil
densities can be obtained by using a sand cone test in the test pit(s) near the
installation site. Alternatively, a standard penetration test (SPT) can be carried out in
undisturbed soils, to measure the penetration resistance in terms of blow counts (N-
value); which can be related to soil characteristics and strength parameters. The
sand cone test can be done at various depths for easy profiling; while, the SPT test
provides a boring log which can be further interpreted for soil density and angle of
internal friction (Brink et al, 1982). The shear strength values, τ can be estimated
once the soil classification and the angle of internal friction, φ has been determined.
For more detail evaluation, sample needs to be tested in laboratory under triaxial
load condition. If the excavated hole is entirely backfilled, the soil properties of the
backfill should be obtained from historical data or specific test in the backfill soil.
Backfill
N=30
G.L.
Native
Soil
Layer 1, N=20
Native
Soil
Layer 2, N=30
Based on the above information, the following properties (IEEE, 1985) are chosen in
figure 3.7 for the native soils and the backfill soils.
Backfill
(N=30)
ρS = 17.2 kN/m 3 G.L.
Native
(N=20)
Soil
Assumed Layer 1,
Failure ρS = 16.7 kN/m3
Surface
Native
Frustum Angle Soil (N=30)
= 25°
Layer 2,
ρS = 17.2 kN/m3
VS = VS , F - VC , B (3.5)
where
A 2
VS , F = [G + B 2 + ( B * G )] (3.6)
3
G= B+ 2 (A*tan ψ ) (3.7)
and
D 2
VC , B = [ B + C 2 + ( B * C )] + (( A − D) * ( E 2 )) (3.8)
3
where VC , B represents the volume of concrete below the ground line and VS , F is
the volume of the frustum. The total concrete volume is given by equation (3.9).
D 2
VC = [ B + C 2 + ( B * C )] + (( A − D) * ( E 2 )) + ( E 2 * F ) (3.9)
3
The above formulae are approximate and are based on a square base foundation. If
the base is rectangular, the appropriate adjustment should be made to equations
(3.6) to (3.9). The calculated geotechnical design capacity for the foundation type
shown in figure 3.2 is 375 kN. A sensitivity study was done by varying the failure
angle and the result is presented in Table 3.2.
Since this report describes the uplift test for a pad and chimney type foundation we
will describe the steps required to carry out this type of test only. Other foundation
tests such as drilled shaft under moment may require additional information and can
be obtained from the test report (EPRI 2197, 1982). However, the overall steps would
be very similar to those described here.
Section 2.5 has provided the required background information on how to estimate the
sample size for cases where the standard deviation of foundation test capacity data
is either known or unknown. The user can use this information as a guide to
estimate the sample size for specified confidence level on the estimated mean test
capacity value of foundation test capacity data (figure 2.4). For sample size less than
30, a correction needs to be made based on Student–t distribution.
As mentioned earlier, the user needs to know whether the tests will be done at one
specific site or at several sites along the transmission line route (Figure 3.1). The
key parameter for selecting a test site could be the access as well as its typical soil
characteristics and distributions with respect to line routing.
Depending on the testing program, the subsurface investigation could vary between
simple to very complex. For example, the uplift test program requires the
determination of a number of soil strength parameters. These are: soil types
(cohesive versus. cohesionless), relative densities, unit weight and the c − φ
properties. This information can be obtained from a sand cone test, SPT test or any
other suitable test program (laboratory or in-situ).
The foundation is excavated by machine to the width ‘C’ in Figure 3.2. The bottom of
the excavation is widened to the width ‘B’ using a handheld shovel. Once the
foundation is in place, the surrounding hole is backfilled with well compacted granular
materials. The strength parameters of the compacted backfill soil at each layer
(300mm to 450mm) can be determined from a nuclear densiometer test, sand cone
test, etc. The detailed installation procedure for various foundation types can be
found in CIGRE (2006) technical brochure.
Figure 3.8 presents a loading test arrangement of a grillage foundation under inclined
uplift load. The test was arranged to determine the anchor capacity of a grillage
Figure 3.8 - Grillage foundation under uplift load (Prasad and Haldar, 2001)
The loading sequence can be simply one load cycle where the load is applied until
the foundation fails or it could be cyclic where it follows a loading-unloading
sequence. The rate of loading can also be varied. In this study, a cyclic loading was
used.
For typical uplift test, the user needs to record the load and displacement of the
foundation. The load can be measured either using a dynamometer or electronically
through a load cell. Similarly, the displacement can be measured using a survey
transit, dial gages or using LVDT. The engineer can use any one of these or all
depending on the accuracy and the verification desired in the measurement process.
This will include taking out any offsets or any other transformations required in the
load and displacement measurements.
Once the data has been reduced, the final test result can be plotted in terms of load
versus displacement as shown in Figure 3.9. In the next section, the user will be
presented with the methodology to determine the ultimate capacity of the foundation.
The capacity is performance based and the method follows IEC 61773.
400
300
200
100
0
0 5 10 15 20 25 30 35
Displacement (mm)
Figure 3.9 - Typical load displacement plot for a pad and chimney foundation
For a given load displacement data set, the utility engineer can determine the elastic
limit capacity and the ultimate limit capacity following the procedures described in the
IEC 61773 entitled “Testing of Structure Foundations for Overhead Lines”. However,
the definition of working limit (damage limit) capacity is more subjective.
Buckley (1994) carried out a number of foundation tests in various soil conditions in
Ireland. The elastic limit load was defined as the load when the foundation movement
upon unloading was insignificant. The working limit load was defined as the load
associated with a 10 mm foundation movement. Details on the selection of 10 mm
displacement are given in section 4.4.2. This is similar to “Damage load” defined in
IEC 60826. The ultimate load was defined as the final pull out of the foundation.
Since one can get a large movement under ultimate condition, it was reported that
this load is not practical to use because the tower may not be able to sustain such a
large movement. Buckley (1994) has recommended a 10 mm displacement as the
determination of the working limit capacity for spread footing. For directly embedded
single pole structure, the capacity associated with a 20 rotation can be used to define
the damage limit (IEC 60826, 2003).
The same is also true if the user chooses to use different displacement criteria for the
working limit capacity or damage limit as per IEC 60826 (2003). IEC Standard
provides various methods to determine the limit capacities. Table 3.3 presents the
comparison of ultimate capacities based on various methods described in IEC 61773.
Slope
Tangent
90% Tangent
Methods Intersection Parabolic Hyperbolic
Criteria (ST)
(TI)
4 10
mm mm
Estimated
Ultimate
323 542 620 625 362 450
Capacity
(kN)
In the above sections, a specific example (figure 3.9) is presented to determine the
various test capacities of an undercut pad and chimney foundation tested by ESB in
Ireland. However, depending on the method used, the estimated ultimate capacity of
this foundation can vary considerably; it could range from 323kN to 625kN. The
mean value is 487kN with a coefficient of variation of 26.5%.
The calculated geotechnical design capacity can also vary significantly depending on
the specific failure model used. As stated earlier, the overall variation of the
foundation capacity is a direct function of the variation of soil strength parameters.
The uncertainty in the model predicted geotechnical design capacity, Rn can be
accounted theoretically only through the joint density function of all these parameters
in a probabilistic sense. However, in practice this is not possible because the joint
density function of all these parameters is often complex and is unknown. Obviously,
Test Capacity
the designer can obtain different ratios of m =
Geotechnical Design Capacity
depending on what methods are used in determining these limit capacities.
Let us assume that the design engineer wants to test a pad and chimney foundation
(figure 3.2) to determine the characteristic capacity with a defined exclusion limit. The
engineer has a number of choices. The engineer can test the foundation either in one
specific site (Site A in Figure 3.1), or can arrange a testing programme along the line
route at different site locations. The test can be repeated at each site to obtain a
dispersion of the test capacity (Sites A, B and C in Figure 3.1), which can then be
used in determining the confidence level. Also, by testing at various locations along
the line route, the engineer will obtain the limit capacities, which are more
representative of the entire line route.
3.7.1 Site specific repeat test data for single geotechnical design capacity
( Rn )
The geotechnical characteristics of SITE A are obtained by SPT tests and Table 3.4
presents the data. The dimensions of this test foundation are presented in Table 3.1.
Using equations (3.1) and (3.5 to 3.8), the calculated geotechnical design capacity
based on a ψ =300 frustum angle (Section 3.2.1) is 439kN. To calculate this
capacity, the minimum value of the soil density used is 16.7kN/m3 (figure 3.7). The
concrete unit weight, ρC is assumed to be 22.6kN/m3. The engineer performed three
repeated tests at this site (Site A in figure 3.1) and obtained working limit capacities
at 10 mm displacement as 500kN, 540kN and 600kN respectively. These values can
be obtained from actual tests similar to the plot presented in Figure 3.9. Figure 3.10
depicts the plot, which shows the scatter in the test capacity with respect to one fixed
geotechnical design capacity value.
Relative Angle of
N-Values Internal Unit Weight
Density Density, Dr in
(SPT) Friction, φ (kN/m3)
Percent
Since the test capacities are representative for one single geotechnical design
capacity value, there is no correlation between the test capacity and the geotechnical
design capacity. A correlation defines the degree of measure of the linear
relationship between the test capacity y and the geotechnical design capacity x. In
this case, it relates the linearity between the test capacity and the geotechnical
design capacity. In statistical data analysis, the degree of predictability will often
depend on the quality of the dataset (mutual dependence between test capacity and
geotechnical design capacity) and this is often measured by the statistical correlation
coefficient r. This will be discussed in detail later in sections 4.0 and 5.0.
600
580
Test Capacity (kN)
560
540
520
500
480
0 100 200 300 400 500
Figure 3.10 - Repeated test data for an undercut pad and chimney foundation-
project specific
Alternatively, the engineer decides to test this foundation type at three different
locations along the line route considering that the capacities will be more
representative of the soil conditions along the line route depicted, in Figure 3.1.
These sites are marked in figure 3.1 as Site A, Site B and Site C; respectively. Based
on SPT test data, Table 3.5 provides the soil characteristics and the estimated soil
densities.
The model foundation shown in Figure 3.2 is again used to calculate the geotechnical
design capacity for the three different soil types. Using the frustum angle model, the
calculated geotechnical design capacities are 379kN, 439kN and 498kN respectively
for loose, medium and dense soil conditions. The damage capacities at 10 mm
displacement are assumed to be 450kN, 540kN and 680kN respectively. Figure 3.11
Test Capacity
mi = (3.10)
Geotechnical Design Capacity
600
Test Capacity (kN)
500
400
300
450 degree
200 line
100
0
0 100 200 300 400 500 600
Figure 3.11 - Uplift test -scattered test data along the line route (damage limit
capacity at 10mm displacement
3.8 Summary
This section presented a review of the various steps involved in the assessment of
foundation capacities from test data as well as from theoretical model predicted
design capacities. The latter is defined as geotechnical design capacity value. This
includes: (a) determination of the geotechnical design capacity based on specific
model(s), (b) field test programme to determine the test capacity and associated
strength parameters for the model, (c) sample size requirement for the test
programme, and (d) testing programme and interpretation of data according to IEC
4.1 General
This section of the report presents a framework to calibrate the foundation strength
factor based on a limited number of full-scale test data. The methodology is based
on a simple mathematical model where the statistical properties of mi values (Table
3.6) are used following Buckley (1994). Later, various assumptions of this model are
examined and a linear regression model is proposed for general application. The
correlation effect is discussed and its importance in the data analysis is emphasized.
The uncertainty in the prediction is primarily introduced through two factors: (i) model
error and (ii) data error due to finite sample size. However, to develop a general
methodology for the calibration of the foundation test data against the model
predicted geotechnical design capacity (geotechnical foundation design model), the
foundation designer needs to assess the quality of the dataset (correlation) and its
influence on the predicted test capacity. The predicted test capacity is analogous to
the characteristic capacity described in section 2.4.
In statistical data analysis, the degree of predictability will often depend on the quality
of the dataset (mutual dependence between test capacity and geotechnical design
capacity) and this is often measured by the statistical correlation coefficient, r.
∑ ( xi − x )( yi − y )
r= (4.1)
(n − 1) s x s y
where s x , s y are the standard deviations of the geotechnical design capacity values
xi and the test capacity values yi ; x and y are the mean values of the geotechnical
design and test capacities respectively. The last two symbols are also analogous to
Rn (mean geotechnical design capacity value) and R (mean test capacity value).
The variables x (or y ) and s x (or s y ) are computed from equation (2.1) and
(2.2) replacing Ri and R by xi (or yi ) and x (or y ) respectively. For convenience,
x and y will represent geotechnical design capacity and test capacity values. For
example, figures 4.1 (a) and (b) present two datasets with very different correlation
properties.
Two scenarios are considered. In the first scenario, it is assumed that the tests have
been carried out only at specific site with one type of foundation (repeated test data)
while in the second scenario the data is available for multiple test sites along the
transmission line route. This is identified as the scattered test data (refer to section
3.7).
The m - model of Buckley is based on the ratio of the test capacity and the
geotechnical design capacity.
n
y
∑ ( xi )
i
m= 1 (4.3)
n
The standard deviation and the coefficient of variation of the mi values are
computed from equations (4.4) and (4.5)
1 n
2
sm = ∑
n −1 1
(mi − m ) 2 (4.4)
s
vm = m (4.5)
m
Since xi = x , for the repeated test data presented in Section 3.7.1, the correlation
coefficient is zero ( r = 0.0 ) from equation (4.1) indicating that the data is unique.
Therefore, the exclusion limit for the repeated test data is computed from the test
data only. In this case, the three data points are 500kN, 540kN and 600kN
respectively (figure 3.10). The geotechnical design capacity for all three tests is
439kN (reference to section 3.7.1).
Using equations (2.1-2.3), the mean test capacity R = 556 kN and the coefficient of
sR 115.9
variation, vR = = = 0.20
R 556
Therefore, the characteristic capacity with a 10% exclusion limit (90% confidence
level, one sided), is obtained from equation (2.11):
RC ,10% = R (1- k α vR )
Test Capacity, y
Data Banded
Around The
Line With
High s xy
where tα ,ν = 1.90 from Table 2.2. Therefore, φF = 0.785 considering sample size
effect; correspondingly the foundation strength factor is decreased by 20% when the
sample size effect is considered.
Since the data for the repeated tests (reference to figure 3.10) are only valid for one
specific test site and for one size of foundation, it has no correlation with respect to
other test site locations along the line route unless the foundation designer can be
sure of another location where the subsurface information is similar to the test
location and the foundation has a similar aspect ratio (depth/width). For undercut
foundation with rectangular base, an equivalent diameter D0 may be considered as
alternative for the width.
In the second case when the data is scattered and the tests are undertaken at three
different locations along the transmission line route (sites A, B and C in figure 3.1),
the foundation designer can expect that the data from various test sites will have a
degree of correlation. However, the degree of correlation will depend on how well the
geotechnical design model can predict the foundation capacity at the various test
Based on figure 3.11 and the mi data presented in table 3.6, it is obvious that the
geotechnical design capacity based on a 300 frustum angle model is under-predicting
the test capacity (working limit at 10 mm), since all three geotechnical design
capacity values are below the 450 line passing through the origin. The mean, m
(slope) and the coefficient of variation, vm of these tests are obtained from equations
(4.2 – 4.5) and are 1.26 and 0.085 respectively. Figure 4.2 presents the average
slope, m . This means that on average, the model predicted geotechnical design
capacity underestimates the test capacity by 26%.
Rc
φF = = m (1- kα * vm ) (4.7)
Rn
The characteristic capacities for 5% (ASCE) and 10% (IEC 60826) exclusion limits
can be obtained initially first from equation (4.7) neglecting the sample size (CIGRE
2002). The kα values for 5% and 10% exclusion limits (95% and 90% confidence
levels) are obtained from Table 2.2 and these values are 1.645 and 1.282
respectively.
and
For finite sample size, (n less than 30, Section 2.1), kα should be replaced by tα ,ν
in equation (4.7) and
Rc
φF = = m (1- tα ,ν * vm ) (4.8)
Rn
The tα ,ν value is obtained from Table 2.2. Table 4.1 summarizes the foundation
strength factor values with the plot shown in figure 4.2. Correspondingly, the
foundation strength factor values are reduced by 6% and 15% when the sample size
effect is considered (ν = n − 1 = 3 − 1 = 2) .
and
1.26
600 kα =1.28, sample
size neglected
500
400 90% lower confidence
300 level-10% exclusion
200 limit, tα,ν =1.90,
100 sample size included
0
0 200 400 600
Geotechnical Design Capacity(kN)
4.3.1 Discussion
The simplicity of the above model (equation 4.3) is that once the relationship
between the characteristic capacity, RC and the geotechnical design capacity
value, Rn is established through the foundation strength factor, φF , the geotechnical
design capacity is easily determined from the characteristic capacity with a defined
exclusion limit (e %). However, the method is based on the following five key
assumptions.
• The relationship between the foundation test capacity versus model predicted
geotechnical design capacity is derived from a linear regression model.
• The intercept is assumed to be zero, i.e., the regression line passes through
the origin.
• The residual error, ε i (the difference between the estimated mean test
capacity and the actual test capacity), is non constant, i.e. it increases as the
• The sample size, n is large and the data follow a normal distribution.
• Correlation of the data is high and the coefficient of variation, vm is small (<
0.10).
A low correlation value, r2 would imply that there is no linear relationship between the
foundation test capacity and the geotechnical design capacity. Table 4.2 presents
some suggested correlation values and the degree of linear association between the
test capacity, y and the model predicted geotechnical design capacity, x.
Linear association
Correlation Value, r 2
between x and y
weak
Low, r 2 ≤ 0.3
moderate
Moderate (0.3 ≤ r 2 ≤ 0.8)
strong
High, r 2 > 0.8
Each of the above assumptions requires to be clearly understood. There are cases
where these assumptions may not apply. For example, the assumption that the line
passes through the origin needs to be validated. The assumption that the residual
errors vary in a non-constant manner also needs to be validated. Without these
validations, the lower limit on the prediction interval of the test capacity (analogous to
characteristic capacity with exclusion limit, section 2.4) for a single geotechnical
design capacity value, Rn and the extrapolation of the data across any geotechnical
design capacity values (data along x-axis) could be problematic and could lead to
erroneous conclusions. This will be discussed in detail in section 5 under “linear
regression model” and in section 6 under “interpolation versus extrapolation”.
4.4.1 Background
Electricity Supply Board (ESB), Ireland carried out a total of 68 full-scale foundation
tests during the period 1984 to 1989. Fifty four (54) of these foundations were of the
pad and chimney type and cast, where possible against undisturbed soil and
undercut where feasible. These tests were carried out on foundations newly installed
primarily for research purposes; while the remaining fourteen (14) old foundations
were also pad and chimney type and made available for testing because of the
removal of an existing 40 year old line. These were installed in the disturbed soils (c-
φ ). All of the fifty four (54) foundations tested were installed in various types of soils
and using different installation techniques. Buckley (1994) has reported in detail the
foundation test results and the next paragraph is a direct quote from his recent
communication (Buckley, 2007).
Except for the existing 14 foundations, all the research foundations were installed in
cohesive and/or cohesionless soils. Nine (9) of these installed foundations were in
granular soils where the foundations were above the ground water table; while 22 of
the research foundations were in soils partly below the ground water table and the
soils were overconsolidated glacial tills (c- φ soils). The soil conditions at the existing
fourteen foundation sites varied from cohesive to silty soils (c- φ ). The majority of
these foundations was pad and chimney type (figure 3.2) and is classified as
undercut foundations.
All tests were completed to determine the ultimate capacity of the foundations under
uplift condition. Uplift loads were applied using high capacity cranes, while the
horizontal (shear) loads were applied via a tractor mounted winch. The typical testing
period was 20 minutes and in some instances loads were cycled. The tests were
carried out until the foundation failed, although in some instances the ultimate
foundation capacity could not be reached due to limitations of the test equipment.
The ESB study defined the foundation failure capacity in terms of the foundation
movement. These were: (i) elastic limit load, (ii) working limit load and (iii) ultimate
limit load. The elastic limit load was defined as the load when the foundation
movement upon unloading was insignificant. The working limit load was defined as
the load associated with a 10 mm foundation movement (similar to “Damage load”
defined in IEC 60826) and was based on the following criteria:
• The load was normally 60% -70% of the ultimate load determined based on
the other criteria such as the tangent intersection, slope tangent etc (IEC
61773).
• A number of these foundations initially tested with the lower capacity crane
equipment (50 tonnes) experienced movement more than 10 mm but did not
fail. These were later tested again to failure with the larger equipment (200
tonnes) and showed the same results under the repeat tests.
• Earlier work by (Mors, 1964) indicated that the tower can sustain a 10 mm
differential leg movement without any damage to the tower.
The ultimate load was defined as the final pull-out capacity of the foundation. Since
usually the foundation displacement under ultimate condition is large, it was reported
that is not practical to use this load because the tower would not be able to sustain
such a large displacement. Therefore, in all subsequent data analyses, a working
limit load at 10 mm is used as the ultimate capacity.
The nine data points for the “Category 2” foundation test are analyzed here using the
m -model (equations 4.7 and 4.8). Table 4.2 provides the data, i.e. the geotechnical
design capacity based on a 300 frustum angle model (section 3.1.1) and the
measured test capacity at 10 mm displacement (working limit capacity) for each test.
This data is taken directly from ESB report. The unit weights used are 22.56kN/m3
for concrete and 15.7kN/m3 for soil respectively. Table 4.3 also presents two
dimensions for each foundation, the depth and the equivalent diameter at the base.
The equivalent diameter was used since all of the tested foundations had a
rectangular base (Buckley, 2001).
Foundation Type ML1A ML1B ML1C ML1D ML2A ML2B ML4B ML8A ML8B
Depth-A (m) 1.5 1.7 1.9 2.1 2.3 2.3 2.7 2.2 2.3
Equivalent Base Diameter-D0 (m) 1.53 1.34 1.63 1.49 1.54 1.52 1.65 1.66 1.59
A 0.98 1.27 1.65 1.41 1.49 1.51 1.63 1.32 1.45
Aspect Ratio- ( )
D0
Model Predicted –Geotechnical 154** 163 249 270 332 328 489 332 345
design capacity (kN)
Test Capacity at 10 mm 204* 248 364 416 389 306 680 383 418
displacement (kN)
Test Capacity 1.32 1.52 1.46 1.54 1.17 0.93 1.39 1.15 1.21
mi =
Geotechnical Design Capacity
* Numbers have been rounded to nearest digit
Table 4.4 presents the data statistics. The mean values of the test and the
geotechnical design capacity are 378kN and 295kN respectively, while the coefficient
of variation values are approximately 35%. The m and vm values are 1.30 and
Average Coefficient of
n x sx y sy m vm Aspect Variation of
Ratio Aspect Ratio
Table 4.5 presents the foundation strength factor values for 5% and 10% exclusion
limits (prediction with 95% and 90% confidence levels respectively).
Table 4.5 – Foundation Strength Factor Values for 9-point dataset (ESB
Data)
4.5 Summary
This section of the report presented a simple mathematical model to compute the
foundation strength factor following Buckley (1994). The mathematical model
requires the statistical properties of the m-values defined as the ratio of the
foundation test capacity and the geotechnical design capacity. The various
underlined assumptions of the model are presented and a need for a more robust
linear regression model has been identified when the specific assumptions are not
met. This is presented in Section 5 under the linear regression model (“alternative
method”). An example design problem for a pad and chimney foundation is used in
section 5 to demonstrate the applications of simple model as well as the ‘alternative
model’ based on a linear regression analysis.
5.1 General
The prediction interval is very similar to the confidence interval estimate for the
mean test capacity value (reference to section 2.3); except that the prediction interval
is used to predict a foundation test capacity (new value) for a given single foundation
geotechnical design capacity value. The lower limit on the prediction interval in the
linear regression model is analogous to the exclusion limit on characteristic capacity
described in section 2.4.2. This linear regression model can be applied to a wide
variety of foundation test data, where the correlation value may not be high (table
4.1) and the intercept from the linear regression equation may have an effect.
5.2.1 General
In the linear regression model, the test capacities (y-values) are plotted against the
foundation model predicted geotechnical design capacities (x-values) and a straight
line is fitted through the data to ensure that the residual error, ε i , is minimum. Figure
5.1 depicts a least square fit of a typical test dataset. The residual error, ε i given in
equation (5.1), is defined as the difference between the estimated mean test capacity
)
value, yi (given in equation 5.5), and the actual test capacity value, yi and is shown
in Figure 5.1.
)
ε i = yi - yi (5.1)
This residual error, ε i , should be minimized to ensure that there is a ‘best fit’ line to
the dataset. The regression model presented is for the case where each test
b=
∑ wi ( xi − x ) yi (5.2)
∑ wi ( xi − x )2
a = y − bx (5.3)
where wi is the appropriate weight applied to each data point (Ang and Tang, 1975).
Normally, the variance of residuals in the linear regression model is assumed to be
constant and wi =1.0, unless the residual varies in a non constant manner; this point
will be discussed later. The sum is carried out based on a sample size of n.
The slope, b, can also be computed from the following relationship if we know the
correlation coefficient, r from Equation 4.1:
sy
b=r (5.4)
sx
The estimated mean test capacity value along the regression line corresponding to
)
each xi , is yi , and can be calculated from equation (5.5) as:
)
yi = a + b( xi ) (5.5)
110
Error, ε
90
Estimated, y
80
70
Intercept, a Slope, b
15 20 25 30 35 40 45 50 55 60
n
)
∑ wi ( yi − yi )2
s xy 2 = 1 (5.6)
n−2
)
where yi is the estimated mean test capacity value of yi , and is based on Equation
(5.5) and wi =1.0. For non constant variance of residuals where the residuals vary
with xi values, wi is dependent on the predictor values along the x-axis (Ang and
Tang, 1975). Depending on the weight, the unit of s xy could be dimensional or non-
dimensional.
The variance of residuals can be also estimated directly from the data without doing
a regression analysis and is given by equation (5.7):
n −1 2
s xy 2 = [ ]s y (1 − r 2 ) (5.7)
n−2
For large n,
s xy 2 → s y 2 (1 − r 2 ) (5.8)
) 1 ( x0 − x ) 2
CI = y ± kα / 2 * s xy [ + ]0.5 (5.9)
n n
∑ wi ∑ wi ( xi − x ) 2
1 1
In this equation, wi = 1.0 because the variance of residuals is constant. For finite
sample size, kα / 2 should be replaced by tα / 2,ν where ν = n-2 degrees of
freedom. The equation is very similar to equation (2.6) except that the estimated
)
mean test capacity y is a conditional mean given a geotechnical design capacity
value x0 . The term s xy is the standard deviation of the residual error of the combined
data (square root of equation 5.6), which includes both the actual test capacity and
For repeated test data at a specific site, x = x0 (sections 3.7.1 and 4.1), r = 0.0
(equation 4.1), b= 0.0 (equation 5.4), a= R (mean value of the test data itself,
equation 5.3), s xy = sR (square root of equation 5.8) and the equation (2.6) is
obtained from equation (5.9). Therefore, the confidence level on the mean value of
the test capacity data referred in section 2.3 (equation 2.6, for test capacity data
only) is a subset of equation (5.9) under these specific conditions.
2
The prediction interval is not only dependent on the variance of residuals s xy , but
)
also on the uncertainty in the estimated mean test capacity value, y given a
geotechnical design capacity value, x0 (equation 5.5). The lower limit of the
prediction interval of the predicted test capacity y p , based on a single geotechnical
design capacity x0 , (same as Rn ) is obtained from equation (5.10)
) 1 ( x0 − x ) 2
y p = y - kα s xy [ 1 + + ] 0.5 (5.10)
n n
∑ wi ∑ wi ( xi − x ) 2
1 1
The equation (5.10) is very similar to equation (2.8) except that the prediction is
based on the lower limit of the interval given a geotechnical design capacity value x0 .
This is depicted in figure 5.2 where it is shown that for a specified geotechnical
design capacity value x0 , the lower limit of the prediction interval point B, is similar to
point L, in figure 2.3 (the characteristic capacity RC ). The only difference is that the
computation of the lower limit of the prediction interval in figure 5.2 is based on a
given conditional geotechnical design capacity value x0 ; while in the latter case, it is
done on the test data directly (figure 2.3).
For repeated test data at a specific site x = x0 (sections 3.7.1 and 4.1),
r = 0.0 (equation 4.1), b= 0.0 (equation 5.4), a= R (mean test capacity value of the
test data itself, equation 5.3), s xy = sR (equation 5.8) and the lower limit of equation
140
130
Prediction Interval
Estim ated M ean Test Capacity,
120 y^ at x o
Test Capacity, y (kN)
110
Residual
Error, ε
*
100
A
Figure 5.2 - Regression plot with confidence interval on mean test capacity
values and prediction interval for new test capacity values
Equation (5.8) indicates that the variance of residuals is dependent on the variance
2
of the test data s y (test capacity data along y-axis) and the correlation value r 2
between the test data and geotechnical design capacity. The correlation value can
be obtained from equation (4.1). If a dataset has a high degree of correlation ( r 2 →
1.0), then every test data point is almost equal to the model predicted geotechnical
design capacity (refer to figure 4.1b). In this case, the slope of the line b → 1.0
and y ≅ x and s y ≅ s x . The intercept will be close to zero (equation 5.3). In this
case the user can assume that the line passes through the origin. The correlation
value will reflect the realistic characteristic of a design model if the data covers
geotechnical design capacity values both in the low (near x =0.0) and high ranges
within a sample size.
To validate the linearity issue, the user needs to plot the residual errors ε i (equation
)
5.1) against the estimated mean test capacity values yi obtained from equation
(5.5). A typical graph for the residual error values plotted against the estimated mean
test capacity values ( ŷi values) is shown in figure 5.3. The residual errors should
Linearity Test
300
200
Residual Errors (kN)
100
0
-100 0.0 200.0 400.0 600.0 800.0
-200
-300
-400
Estimated Mean Test Capacity (kN)
As the correlation value tends to decrease ( r 2 ≤ 0.8, table 4.2) with increasing
coefficient of variation of the model predicted error ( vm ≥ 0.1), consideration should
be given to include the intercept effect in the prediction interval calculation. For
foundation calibrations where the dataset has moderate correlation ( r 2 ≤ 0.8 ), the
)
intercept has a significant influence in the estimated mean test capacity value y i . In
this case, by forcing the line to pass through the origin arbitrarily, the foundation
designer will make the prediction optimistic (less conservative); while in reality the
designer will not be able to validate the predicted foundation test capacity with a very
high confidence level.
For example, if a database has a poor correlation ( r → 0.0 ) and the second
assumption (zero intercept) is used arbitrarily, this will make the results acceptable
when in reality the test capacity values along the y-axis have no linear relationship
with corresponding geotechnical design capacity values along x –axis, due to a lack
of correlation (Table 4.1). If the foundation design engineer considers the feasibility
) )
in reducing the model from yi = a + b( xi ) to yi = m ( xi ), where the line passes
through the origin, then in principle a null hypothesis ( H 0 ) should be tested
simultaneously at a specified significance level, α (Benjamin and Cornell,1970). A
null hypothesis implies that the proposition
H 0 : a=0.0 and b= m
A specific test can be used to reject simultaneously a null hypothesis that the
intercept, a = a* (arbitrary value) and the slope b = m ( H 0 : a= a* and b= m ) using
equation (5.11).
*
If the calculated F value is greater than the allowable value Fα ,2, n − 2 obtained from
F-distribution, with a significance level, ( α ), and ν = n − 2 degrees of freedom then
the hypothesis is rejected (Bowker and Lieberman, 1972, Sneddon, 2003). In this
case, the intercept should be included in the analysis and restrictions imposed on
how far the data can be extrapolated. This is particularly important for the
database where the correlation is moderate ( 0.3 ≤ r 2 ≤ 0.8, table 4.2) and the data
has large coefficients of variation, i. e. all values of vm , vx , v y are considerably
greater than 0.1. This situation can also arise when the database does not contain
enough points near x=0.0 (low end). It is important that the foundation design
engineer collects test data which at least will cover the low as well as high ranges to
ensure that the correlation value truly represents the main characteristics of the
design model.
The foundation designer should validate that equation (4.7) is satisfied for a* = 0.0
and b = m for the m model presented by Buckley (1994). This will be discussed
later in the example problem presented in section 5.6.
In the linear regression model, it is normally assumed that the residual errors across
all geotechnical design capacity values remain fairly constant as depicted in figure
5.4a. However for many datasets, it is possible for the residual errors to vary linearly
or non linearly with the geotechnical design capacity values along the x-axis. In this
case, it is known as non constant variance of residuals. Figure 5.4(b) presents this
Constant Variance
of Residuals
ε)
Residual Errors ( Y
ε = 0.0
Non Constant
Variance of
Residuals Fan
ε)
Shaped
Residual Errors (
For non constant variance of residuals, the assumption is that the variance of
residuals is proportional to the square of the x-values and an appropriate weight wi =
1
is applied to the original dataset and the regression analysis is carried out on the
xi 2
modified data (Ang & Tang, 1975). Figure 5.5 presents a typical residual plot against
the geotechnical design capacity values for a large dataset.
200
Residual Errors(kN)
100
0
0 200 400 600 800
-100
-200
-300
-400
Geotechnical Design Capacity (kN)
Figure 5.5 - Non constant variance of residuals trend for a larger dataset
5.3.1 General
The basic linear regression model can be used to provide the confidence interval (CI,
equation 5.9) on the estimated mean test capacity ŷ (equation 5.5 for xi = x0 ) , as
well, as the lower limit of the prediction interval of a new test capacity y p (equation
5.10), given a single geotechnical design capacity value x0 (figure 5.2)
The foundation strength factor φ F is defined as the ratio between the lower limit of
the predicted test capacity y p and the geotechnical design capacity value x0 used to
determine this capacity (figure 5.2, Haldar, 2001, Ang, 2002). If the predicted test
capacity with a specified confidence level is y p and the geotechnical design
capacity value is x0 , then the foundation strength factor is:
yp
φF = (5.12)
x0
So far we have discussed two scenarios; these are: (1) the line passes through the
origin and (2) the line does not pass through the origin and therefore has an
intercept. Since there are a further two options (constant versus non constant
variance of residuals), under each of these conditions (zero or non zero intercept
condition), the foundation designer will end up with four possible cases for analyzing
the data. This is described in Table 5.2. Each case in Table 5.2 will provide a lower
The foundation designer can still obtain the foundation strength factor value using
equation (5.12) provided the predicted test capacity is known from any one of the
cases in Table 5.2. Equation (5.11) provides a definite test whether to include the
intercept in the analysis or not. However, the test for constant variance of residuals
versus non constant is not quantitative and is rather subjective because the designer
needs to plot the residuals against the geotechnical design capacity values and to
derive a trend line similar to figures 5.4a and 5.4b respectively. If it is not clear that
the residuals vary either in a constant or non constant manner then the designer will
end up with Cases A and B as one group or Cases C and D as the other group
respectively. This will provide the upper and lower limits on the foundation strength
factor values or in the predicted test capacity values. An example design problem
based on ESB data is shown later to indicate how the designer can select a
particular option.
The following section summarizes formulae for the four cases that the user may face
in analyzing the data. The formulae are presented to predict the lower limit of the
foundation test capacity y p for a given geotechnical design capacity value x0 for
conditions with or without zero intercept and with or without constant variance of
residuals. All these formulae also consider the correlation issue and the coefficients
of variation of the data. The various parameters such as slope b, variance of
2
residuals s xy , estimated mean test capacity ŷi , confidence interval (two sided) on
mean test capacity values CI and the lower limit of the predicted test capacity for a
single geotechnical design capacity value x0 with a given confidence level y p are
presented for each case. Often, the predicted test capacity value will depend on the
2
intercept a, slope b, variance of residuals s xy , the correlation coefficient r and the
weight wi used to estimate the variance of residuals. Since the intercept can be
either zero or can have a finite value and the variance of residuals can vary either in
a constant or a non constant manner (figures 5.4a and 5.4b), the designer will likely
encounter at least two of the four cases as presented in Table 5.2. A flow diagram is
also presented later to provide a guidance to the design engineer in the data
analysis process.
This case is applicable to problems where the intercept requires to be included and
the residuals plotted against the geotechnical design capacity values produce a band
around the ε i = 0 line (Figure 5.4a, constant variance of residuals). The slope b,
intercept a, estimated mean test capacity value ŷi , variance of residuals
5.3.3 CASE B - Intercept is non zero and the variance of residuals is non
constant
This case is applicable to problems where the intercept needs to be included and the
residual does follow a non constant variation according to Figure 5.4b. In this case,
there is a need to include both intercept and slope in estimating the prediction
interval. The slope b, intercept a, estimated mean test capacity value ŷi , variance of
2
residuals s xy , confidence interval (CI) on the estimated mean test capacity value
ŷ and the lower limit on the prediction interval (PI) of the predicted test capacity y p
given a single geotechnical design capacity x0 can be obtained from equations
(5.2),(5.3), (5.5),(5.6), (5.13) and (5.14) respectively (figure 5.2) with the appropriate
1
weight, wi = .
xi 2
1 ( x0 − x ) 2
Confidence interval CI = ŷ ± kα / 2 * s xy [ + ] 0 .5 (5.13)
n n
∑ wi ∑ wi ( xi − x ) 2
1 1
2
x 1 ( x0 − x ) 2
Predicted test capacity y p = ŷ - kα * s xy [ 0 + + ] 0.5 (5.14)
1 n n
∑ wi ∑ wi ( xi − x ) 2
1 1
This case is applicable to problems where the intercept can be neglected based on
equation (5.11) and the residual does not increase with geotechnical design capacity
values (Fig. 5.4a). In this case, we only need to include the slope in estimating the
prediction interval. The weight wi = 1.0 , because the variance of residuals is
2
constant. The slope b, variance of residuals s xy , estimated mean test capacity
value ŷi , confidence interval (CI), on estimated mean test capacity value ŷ and the
lower limit on the prediction interval (PI) of the predicted test capacity y p , given a
single geotechnical design capacity value x0 with a specified confidence interval can
be obtained from equations (5.15), (5.16), (5.17), (5.18) and (5.19) respectively
(figure 5.2 with the intercept being zero).
b=
∑ wi ( xi ) yi (5.15)
∑ wi ( xi )2
( x0 ) 2
Confidence interval CI = ŷ ± kα / 2 * s xy [ ] 0.5 (5.18)
n
∑ ( xi ) 2
wi
1
( x0 ) 2
Predicted test capacity y p = ŷ - kα * s xy [ 1 + ] 0.5 (5.19)
n
∑ ( xi )2 wi
1
5.3.5 CASE D - Slope only with non constant variance of residuals (Similar
to m -model)
∑ wi ( yi − yˆi )
2
2
s xy = = sm 2 (5.21)
n −1
1 0 .5
Confidence interval CI= ŷ ± kα / 2 * sm x0 [ ]
n
v
= m (1 ± kα / 2 m ) x0 (5.23)
n
1 0 .5
Predicted test capacity y p = ŷ - kα * sm x0 [ 1 + ]
n
Equation (5.24) is similar to equation (4.3) presented by Buckley (1994) but includes
1 0.5
a factor [1+ ] which was omitted in the earlier study. For small sample size, this
n
may create differences in the foundation strength factor values. For large sample
size, the term under the second bracket in equation (5.24) is close to 1.0. This is
similar to the difference in the predicted test capacity value and the characteristic
capacity presented in section 2.4. The foundation strength factor obtained from
equation (5.24) will be constant across all model predicted geotechnical design
capacity values (Sneddon, 2003). If this is not the case, then the foundation designer
may find that the foundation strength factor will be a function of the geotechnical
design capacity since the predicted test capacity is also a function of the
geotechnical design capacity (Cases A, B and C in table 5.2).
In all the above equations, kα values should be obtained from Table 2.1 for one
sided confidence level or use tα ,ν when the sample size is finite (n < 30, Table 2.2).
For cases A and B, ν = n- 2 while for cases C and D, ν = n-1.
Equation (5.23) is similar to equation (2.6) where the confidence level is sought on
mean test capacity data only. Equation (5.23) provides the confidence interval on the
estimated mean test capacity (conditional mean value of ŷ given a geotechnical
design capacity value x0 (or Rn ), ŷ I x0 as long as they are related through a
regression line with a slope b= m and a coefficient of variation vm . Equation (5.24)
similarly provides the lower limit on the prediction for a single geotechnical design
capacity value x0 This also corresponds to equation (2.8), where a minimum
strength (characteristic capacity with e% exclusion limit) is obtained from the test
data directly.
The foundation designer needs to consider two questions before the equation (5.24)
(also applicable for the simple model, equation 4.7), can be applied generally. These
are:
• What is the upper limit of vm beyond which the method is not valid?
• How far do we go along the x-axis to extrapolate the data in providing a
prediction interval for the test capacity?
To address the first question, the coefficient of variation cannot be greater than
1
vm ≤ [ ] (5.25)
1
tα ,ν [1 + ]0.5
n
This will ensure that the foundation strength factor is positive and provides a
reasonable value. Of course it can be argued that even with a very high m value, it
is always possible to achieve a positive foundation strength factor value when the
term inside the bracket (equation 5.24) is very low. However, in such a case the
foundation designer should not use the geotechnical design model, because it is
underestimating significantly the test data and correlation will be low (table 4.2).
Table 5.3 presents some typical values of vm to ensure that the foundation strength
factor does not become negative. For other intermediate value, use equation (5.25).
If the test data and the model predicted geotechnical design capacity produce a vm
value greater than the allowable limit shown in Table 5.3, the foundation strength
factor φF could become negative. The foundation design engineer needs to be
careful when the data is extrapolated outside the range. Even when the vm values
are close to the values shown in Table 5.3, the foundation strength factor value will
be low and probably unrealistic indicating that the error is introduced because of the
large coefficient of variation. Similar criteria can also be used for other cases from
the above equations or plotting the data graphically. Section 6 will provide some
guidelines on interpolation and extrapolation.
In this section, the 9-point ESB test data (table 4.3) is re-analyzed using the linear
regression model where two major assumptions that the intercept is zero and the
variance of residuals is non constant are examined. Figure 5.7 presents the
regression plot of the data and the correlation value.
Table 5.4 presents the data statistics. The model, in general, appears to be
conservative because of large bias, i. e. the geotechnical design model on average
underestimates the test capacity value by 30%. The correlation value, r 2 is 0.81 is
high (equation 4.1 - Test 1 in flow diagram, table 4.2). A linear regression model
(CASE A) is run first following the flow diagram in figure 5.6. The standard deviation
yi
Compute mi =
xi
High or Moderate
Correlation
Regression Analysis (Figure 5.2)
Slope (5.2), Intercept (5.3), Estimated Test Capacity
Value (5.5), and Variance of Residuals (5.6)
Slope &
Intercept Test 2: Slope only
Eq. 4.14
5.11
Test 3: Test 3:
Chk Figs. 5.4(a) Chk Figs. 5.4(a)
& 5.4(b) & 5.4(b)
Calculate Foundation
Strength Factors to Provide
Conservative Values -
Equation 5.12
Test 1: Test for Correlation END
__
Test 2: Simultaneous Test for Intercept = 0.0 and Slope = m
Test 3: Visual Comparison to Test Constant Variance of
Residuals vs. Non Constant Variance of Residuals
*Constant COV - Special Case
(Residual values banded around zero line or "fan
- CIGRE (2002)
shaped" - use judgment)
600
Test Capacity (kN)
0
0.0 200.0 400.0 600.0
Geotechnical Design Capacity (kN)
s xy
n r2 x sx y sy a b ∑ xi 2 F*-value m vm
9 0.81 296 102 378 135.4 25.8 1.19 63 0.87E+05 0.164 1.3 0.154
< 3.26
Next, the residual values for the nine data points are plotted against the geotechnical
design capacity values and Figure 5.8 presents this. Two lines are drawn to see a
trend line where the variance of residuals is non constant (similar to figure 5.4 b) and
it is not fully definitive that this is the case. Therefore, according to the flow diagram
in figure 5.6, the designer should try both Cases C and D for determining the
foundation strength factor. Obviously under Case C the foundation strength factor
will not be constant and rather will be a function of the geotechnical design capacity;
because the predicted test capacity value is also dependent on the geotechnical
design capacity value (equation 5.19)
-100
-150
Geotechnical Design Capacity (kN)
The lower limit of the prediction interval with (90%) confidence level is presented for
both cases, in tables 5.5 and 5.6 respectively. Figures 5.9 and 5.10 present the plots.
Figure 5.11 depicts the comparison of foundation strength factor, φF based on
CASE C and CASE D scenarios; it is interesting to note that for CASE D where the
foundation strength factor is constant across all geotechnical design capacity values,
the interpolation may not be a problem. However, for CASE C scenario, where the
residual is interpreted as constant, (figure 5.4a), the foundation strength factor
becomes strongly dependent on the geotechnical design capacity value (x-value)
chosen as the predictor values and further away from the mean value (297 kN).
Therefore, in this situation, the foundation designer may want to consider the lower
value of the two foundation strength factors to ensure that the foundation design is
700
Residuals=3515
600
500 Original 90% Lower
400 Data Points Prediction Line -10%
300 Exclusion Limit
200
100
0
0.0 100.0 200.0 300.0 400.0 500.0 600.0 700.0 800.0 900.0 1000.0
Figure 5.9 - Predicted test capacity with 90% confidence level (Case C, constant
variance of residuals)
1
Table 5.6 (Case D- m =1.30, s xy = sm =0.20, tα ,ν =1.40, wi = )
xi 2
If the foundation designer wants to correct for the small sample size, then tα ,ν value
should be based on Student t-distribution (reference to Table 2.2) and the foundation
strength factor becomes
This value is shown in Table 5.7. The value is slightly different from the one
1 0.5
presented in Table 4.5 because the term [1+ ] in equation (5.24) was not
n
included in the earlier analysis.
700
600
500
400 90% Lower prediction
300 line-10% exclusion limit,
Sample size included
200
100
0
0.0 100.0 200.0 300.0 400.0 500.0 600.0 700.0 800.0 900.0 1000.0
Figure 5.10 - Predicted test capacity with 90% confidence level (CASE D, non
constant variance of residuals)
1.000
0.800
0.600
Constant Foundation
0.400 Data Strength Factor,φ F =
R 1.004
Variable Foundation
0.200
Strength Factor Values,φ F
0.000
0 100 200 300 400 500 600 700 800 900
A base case model for an undercut foundation (refer to figure 3.2 and table 5.8) was
used to obtain the uplift capacity based on a 300 cone model. The unit weights for
soil and the concrete were taken as 15.7kN/m3 and 22.6kN/m3 respectively. The
calculated capacity was 508kN. Next, a sensitivity study was undertaken by
changing the depth while keeping the base width constant. Figure 5.12 presents the
data and it is apparent that the tripling of the depth has a seven fold increase in the
foundation uplift capacity. The equation for computing the depth for a given
geotechnical design capacity is also given on the graph. The correlation value is
0.9999 (figure 5.12).
600
400
y = 73.586x2 + 8.2967x + 30.364
200
R2 = 0.9999
0
0 1 2 3 4
Change in Depth (A in m)
Type A B C D E F
Square 2.5 2.0 1.9 0.2 0.8 0.1
Based
The design equations (2.18) and (5.12) provide a relationship between the predicted
capacity y p (similar to RC ), and the geotechnical design capacity x0 , (same as Rn )
RC = φF Rn
For example, if the foundation design load due to extreme wind is QT = 500 kN, and
the foundation designer wants to provide a 90% confidence level based on nine ESB
tests, then the geotechnical design capacity required is (Table 5.5 , φF =1.00 for
90% confidence, 10% exclusion limit)
QT 500
Rn = = = 500kN
φFC
1.00
Once the geotechnical design capacity is known, the required foundation depth can
be estimated using figure 5.12. For QT = 500 kN, Case C provides a geotechnical
design capacity of 465 kN. This can be computed from equation (5.19) by
substituting y p = QT and solving for x0 . The foundation designer can also compute
the geotechnical design capacity from figure 5.12 for a predicted test capacity, y p =
500 kN following the lower limit of the prediction interval line (green line). The results
show a 7% difference in capacity between Case C and Case D respectively. The
corresponding design depth can be estimated either from figure 5.12 or from the
equation given on figure 5.12 for these two geotechnical design capacity values.
These are: design depths 2.37 m versus 2.47 m (5.0% difference).
5.8 Summary
This section of the report presented a framework to determine the prediction interval
of the foundation test capacity based on a limited number of full-scale test data. The
lower limit of the prediction interval is analogous to characteristic capacity with a
specified exclusion limit (section 2.4) and can be linked to the foundation strength
factor (equation 5.12). The prediction interval is based on the calibration of data on
full-scale test capacity against model predicted geotechnical design capacity. The
6.1 General
The foundation strength factor, φF , obtained from a linear regression model has a
valid range and is normally dependent on the predictor data range during the
calibration process (i.e. geotechnical design capacity values along the x-axis x0 or
Rn values). Any interpolation and/or extrapolation should be done properly to ensure
that the confidence level on the estimated mean test capacity value, ŷ from equation
5.5), and the prediction interval used to determine a new foundation test capacity y p ,
for a single geotechnical design capacity value x0 , are obtained accurately.
6.2 Interpolation
The relationship between the predicted test capacity y p for a given geotechnical
design capacity, x0 and the predictor (geotechnical design capacity, x0 ) used in the
calibration model applies strictly when the predictors have similar values that were
used in the calibration model. If this is the case, then it is called interpolation.
800
750 Original Data Points
700
Regression Line,
Test Capacity (kN)
650
600 slope=1.30
550
500 Interpolation Point, x=400,
450
400 Estimated Mean Test
350 Capacity=512, Equation 5.22
300
250
200 Predicted Capacity, 90%
150
100 confidence level, 401 -Equation
50 5.24
0
0 100 200 300 400 500 600 700 800 900 1000
6.3 Extrapolation
The extrapolation process is making a prediction outside the range of values of the
predictor used in the sample to generate the model. Therefore, the further the value
of the predictor is moved from the mean value, the larger the error in the prediction
(reference to section 5.2.2). Besides, there is no guarantee that the linear
relationship holds outside the data range (Ang and Tang, 1975).
Figure 6.2 presents a case where the extrapolation is done for a geotechnical design
capacity of 800 kN along the regression line (CASE D model). The estimated mean
test capacity is 1040 kN (equation 5.22). On the same graph, the predicted test
capacity (90% confidence level-10% exclusion limit) is also depicted which has a
value of 849 kN (equation 5.24, tα ,ν = 1.90 ) where the slope is only considered (the
line passes through the origin) and the variance of residuals is assumed to be non
constant.
Data
800 points
700
600 Predicted capacity for a
500 geotechnical design
400 capacity, 800 kN with 90%
Regression
300 line,slope=1.30 confidence -Equation 5.24
200
100
0
0.0 100.0 200.0 300.0 400.0 500.0 600.0 700.0 800.0 900.0 1000.0
Geotechnical Design Capacity (kN)
Figure 5.11 presents the foundation strength factor plot for various geotechnical
design capacity values. The geotechnical design capacity values for ESB data varies
However, if the foundation designer believes that the variance of residuals is fully non
constant based on figure 5.8, then a constant foundation strength factor value can be
used for the entire data range (150 kN to 500 kN). Of course to extrapolate outside
the data range, the assumption still is that the linearity is still valid in this range. As
mentioned earlier, this is only valid when the line passes through the origin (no
intercept) and the variance of residuals is strictly non constant. If this is not the case,
the foundation strength factor can be dependent on the geotechnical design capacity
values for certain data range.
6.5.1 Interpolation
The lower and upper ranges of the design value are 400kN to 418kN
6.5.2 Extrapolation
The lower and upper ranges of the geotechnical design capacity value are 800kN to
908kN respectively.
The user needs to be careful when extrapolating the data outside the range,
particularly when the variance of residuals is constant (Case A and Case C). Also,
there are no published guidelines how far one can go to extrapolate the predictor
value. Correspondingly, it is recommended here that for constant variance of
residuals case, we do not go normally beyond the xL - s x or xU + s x where xL and
xU are the lower and upper value of the predictor data range and s x is the standard
deviation of the geotechnical design capacity. The foundation designer therefore
needs to use judgment in extrapolating the data outside the range assuming that the
linearity relationship holds. As the engineer gets more data from future tests, this
data can then be added to the database and the analysis updated. The designer
should collect data which covers the low and high ranges, particularly data near the
low end (x=0.0); which could allow the intercept effect to be minimum provided the
correlation is high. On the other had if the designer collects data only in the high end
and the correlation is moderate, then most likely the intercept effect would be
This section presented the limitation of the data extrapolation and some guidelines
how to do this extrapolation. An example problem is shown using the foundation
strength factor obtained in Section 5.0 (Tables 5.5 and 5.6).
6.6 Summary
7.1 Summary
This report presents the various steps necessary to determine the lower limit of the
prediction interval of the foundation test capacity (similar to characteristic capacity)
and to calibrate the foundation strength factor based on a limited number of full-scale
foundation test data. The calibration is undertaken within the framework of a linear
regression model. Initially, a brief review is presented with respect to confidence
interval on the sample mean of the test capacity values and the prediction interval of
a new test capacity value given a sample data of finite size. The confidence interval
is used in the determination of required sample size for a foundation test programme
while the lower limit of the prediction interval is linked to the characteristic capacity
with a specified exclusion limit. This is followed by a brief review outlining the
deterministic as well as the semi-probabilistic design methodologies such as ASD
(Allowable stress design), IEC 60826 (2003) and ASCE (LRFD).
Once the foundation’s test capacity and the geotechnical design capacity have been
determined for specific site(s), the foundation designer can choose a simple
Test Capacity
calibration model based on m = ratio(s). However, to
Geotechnical Design Capacity
use this simple model ( m -model) the foundation designer needs to ensure that the
correlation value of the data is high (table 4.2) and the coefficient of variation of the
m-values, v m (model error) is low.
For other situations where this is not the case, the foundation designer may need to
carry out a more detailed analysis (full linear regression model) to ensure that the
intercept (ordinate at origin) effect can be neglected. This can be undertaken using a
simple statistical test. If the intercept effect cannot be neglected, then the linear
regression model (an alternative methodology) is also presented on how to include
Finally a design example is presented based on ESB test data on undercut pad and
chimney foundation. A discussion on the limitation of data extrapolation is also
presented.
7.2 Conclusions
• The CIGRE Method relies on two very specific conditions which imply that the
intercept should be zero and the residual (the difference between the actual
test capacity and the estimated mean test capacity) should vary with the
geotechnical design capacity values in an increasing manner.
• In addition, it is shown that the CIGRE method does not consider the error
due to small sample size in calculating the prediction interval. CIGRE method
is similar to Case D (equation 5.24) but does contain an additional
1
term [1 + ]0.5 . Therefore the CIGRE method is modified to include the effect
n
of small sample size. The difference in the foundation strength factor value
could be 2 % to 6% depending on the sample size, n and v m values.
• However; if the correlation value, r 2 is greater or less than equal to 0.8 (table
4.2) and the coefficient of variation v m starts to increase (> 0.10), the
designer needs to carry out a full regression analysis to choose a particular
option and/or options from the flow diagram (figure 5.6) to ensure that the
analysis result provides not only a predicted capacity which produces a
conservative design, but also that the design is optimal (economical).
• Finally, it is shown that the general regression model with the unique test for
“zero intercept” is quite applicable to the analysis of a wide range of
foundation test data with various degree(s) of correlation, coefficient of
variations and m values. The methodology presented here is quite robust
and provides a consistent basis for analyzing a wide range of test data. The
foundation designer is able to review the details of these calculations in order
to understand the intricate relationships of various parameters that contribute
to the estimation of the foundation strength factor, φ F values. The foundation
designer can make an informed decision to choose a particular option or
options in determining the predicted test capacity of the foundation with a
specified exclusion limit.
Benjamin, J. and Cornell, Allin 1970 Probability, Statistics and Decisions for
Civil Engineers, McGraw Hill
Brink, A.B.A, Patridge, T.C. and Williams, A.A.B. 1982 Soil Survey for
Engineering, Clarendom Press, Oxford, 378 p
Davidson, H.L. 1982 Laterally Loaded Drilled Pier Research, Volumes 1 and 2,
EPRI EL 2197, RP 1280-1, prepared by GAI Consultant, USA
Digioia Jr., A.M., 2000 Reliability Based Design and Assessment of Foundations
for Transmission Line Structures, IEEE T & D World Expo, April, 26-28, Cincinnati,
Ohio
Prasad, Yenumula and Haldar, Asim 2001 Full-scale Tests On Grillage Anchor
Foundations for High Voltage Lines, Proc. Canadian Geotechnical Conference,
Calgary
This section will provide the user a step-by-step procedure for analyzing the data
following the flow diagram in figure 5.6. Each step is shown clearly so that the user
can understand the calibration process in determining the foundation strength factor
for a specific foundation dataset. The user is still required to use some judgment in
selecting a particular option or a combination of options to ensure that the foundation
strength factor is not only conservative but also optimum. Details of the calculations
are presented here. All calculations are based on 90% confidence level (10%
exclusion limit). Figure 9.1 presents the three point scattered test data along a
transmission line route (section 3.6.2, reproduced figure 3.11). Figure 9.2 presents
the flow diagram (reference figure 5.6) which will be followed in computing the
foundation strength factor for this dataset. Table 9.1 presents the various parameters
that need to be computed to determine the foundation strength factor. A summary of
the methodology is presented in two steps.
600
Test Capacity (kN)
500
400
300
450 degree
200 line
100
0
0 100 200 300 400 500 600
In the first step, a routine linear regression analysis is carried out by minimizing the
residual errors across the all data points. Underlying assumption of the linear
regression analysis is that the regression considers an intercept and the variance of
residuals is constant. To ensure that the intercept can be neglected, a null hypothesis
is tested and a decision is made either to include or exclude the intercept. In the
second step, residual error values are plotted to test the assumption of constant
variance of residuals. If this is not validated, then an appropriate weight is applied to
each data point and the regression analysis is carried out on the transformed data,
with or without the intercept. Table 9.2 presents the necessary equation numbers for
four cases referred in Table 5.2.
yi
Compute mi =
xi
High or Moderate
Correlation
Regression Analysis (Figure 5.2)
Slope (5.2), Intercept (5.3), Estimated Test Capacity
Value (5.5), and Variance of Residuals (5.6)
Slope &
Intercept Test 2: Slope only
Eq. 4.14
5.11
Test 3: Test 3:
Chk Figs. 5.4(a) Chk Figs. 5.4(a)
& 5.4(b) & 5.4(b)
Calculate Foundation
Strength Factors to Provide
Conservative Values -
Equation 5.12
Test 1: Test for Correlation END
__
Test 2: Simultaneous Test for Intercept = 0.0 and Slope = m
Test 3: Visual Comparison to Test Constant Variance of
Residuals vs. Non Constant Variance of Residuals
*Constant COV - Special Case
(Residual values banded around zero line or "fan
- CIGRE (2002)
shaped" - use judgment)
n
450 540 680
m = (∑ mi ) / n = [ + + ] /3 = [1.19 +1.23+1.37]/3 = 1.261
1
379 439 498
=0.009
s 0.093
v m = rm = = 0.0737
m 1.261
n
R = (∑ Ri ) / n = (450+540+680)/3 = 556.7 kN
1
1 n 1
sR 2 = ∑
n −1 1
( Ri − R ) 2 =
3 −1
[(450 − 556) 2 + (540 − 556) 2 + (680 − 556) 2 ]
= (115.9)2 kN
s 115.9
vR = R = = 0.208
R 556.7
Similarly the mean, standard deviation and the coefficient of variation of x-Values
(Geotechnical design capacity) are x = 438.7 kN, s x =59.5 kN and v x = 0.136. The
correlation coefficient r , value is calculated directly from equation (4.1).
n
1 /(n − 1)[∑ ( xi − x )( yi − y )
r= 1
sxs y
1/ 2[(379− 438.7)(450− 556.7) + (439− 438.7)(540− 556.7) + (498− 438.7)(680− 556.7)]
=
59.5x115.9
= 0.992
Therefore, the correlation value r 2 = 0.983 is high and the linear association
between the test capacity y and the geotechnical design capacity x is strong (table
4.2). Therefore, the dataset passed the Test 1 according to the Flow diagram (Figure
9.2).
a = y − bx
Geotech
Damage Square of
nical )
Test design
Limit yi = a + b( xi ) Residual Residuals
Number Capacity = ( yˆi − yi ) xi 2
capacity = ( yˆ i − yi ) 2
yi
xi
1 379 450 441.4 8.58 73.65 143641
2 439 540 557.3 -17.31 299.65 192721
3 498 680 671.3 8.72 76.17 249004
Sum 1.7e-13 449.49 584366
= ∑
n
)
∑ wi ( yi − yi )2 [(450 − 441) 2 + (540 − 557) 2 + (680 − 671) 2 ]
s xy 2 = 1 = = 449.49;
n−2 (3 − 2)
Therefore, s xy = 21.1;
600
Test Capacity (kN)
400
200
0
0.0 200.0 400.0 600.0
Geotechnical Design Capacity (kN)
Regression line
600 Original Data
with slope
Test Capacity (kN)
Points
only,1.26 (m-bar)
400
Regression Line With
200 Intercept, a=-290 and
slope, b=1.93
0
0.0 200.0 400.0 600.0
-200
-400
Geotechnical Design Capacity (kN)
Figure 9.4 - Linear regression plot (with intercept and without intercept)
where
First Term-
n (a-0.0)^2
Second Term-
2n x (a-0.0)(b- m ) ∑
Third Term- ( xi 2 )(b − m ) 2 F* = ∑ Sum
3(-290.64)^2= 2*3*438(-290-0.0)(1.93-1.26) = - 584366*(1.93-1.26) 2 3378
253415 483114 = 233078
F* =
∑ Sum =
3378
= 3.78 < Fα ,2, n − 2 = 49.50 (Table 5.1)
2 * s xy 2 2 * 212
Therefore, we accept the Hypothesis that the intercept is zero and the slope is equal
to m = 1.261 ( m - model, Section 4); All results are summarized in Table 9.5.
Table 9.7- Data statistics for 3 – point scattered test data problem
n R2 x sx y sy a b s xy ∑ xi 2 F-value m vm
3 0.983 438 59 556 116 -290 1.93 21 584366 3.78 1.26 0.0737
< 49.50
Note:
Similar calculations can be done for ESB data and Table 5.4 can also be
produced following the above steps.
-10
-15
-20
Geotechnical Design Capacity Values (kN)
Geotechnical
Test
Test design
Capacity, wi ( xi ) yi wi ( xi ) 2 yˆi = bxi ( yˆ i − yi ) 2
Number capacity, xi
yi (kN)
(kN)
1 379 450 170550 143641 484 1155
2 439 540 230760 192721 560 425
3 498 680 338640 248004 638 1940
Sum 746350 584366 3520
b=
∑ wi ( xi ) yi =
746350
=1.277
∑ wi ( xi )2 584366
s xy 2 =
∑ wi ( yi − yˆi )2 =
3520
=1760
n −1 2
Therefore, s xy = 41.9
The predicted test capacity is computed from equation (5.19) for kα = tα ,ν (small
sample size n=3)
( x0 ) 2
yp = ŷ - tα ,ν * s xy [1+ ] 0.5
n
∑ ( xi )2 wi
1
For a sample size n=3, tα ,ν = 1.90 and is obtained from Table 2.2 (90% value for one
sided interval). Therefore, the second term in the above equation for x0 =200 kN
becomes
( x0 ) 2 200 2 0.5
P= tα ,ν * s xy [ ] 0.5 = 1.90 * 42 *[1+ ] = 82.39
n 584366
∑ ( xi )2 wi
1
Table 9.9 presents the data for four geotechnical design capacity values.
Geotechni Estimated ( x0 ) 2 y p = ŷ -P
cal design Mean Test P= tα ,ν * s xy [ 1+ ] 0 .5
n (equation 5.19)
capacity,
x0
Capacity
Value
∑ ( xi )2 wi
1
yˆ = bx0 (second term in equation 5.19)
(first term in
equation 5.19)
200 255 82 173
400 511 90 421
600 766 101 665
800 1022 115 907
800
720 Regression
640
Line Predicted Capacity -
560
480 90% Confidence
400 level
320
240
160
80
0
0.0 100.0 200.0 300.0 400.0 500.0 600.0 700.0 800.0 900.0 1000.0
1
CASE D ( wi = ) – (equations 5.20 -5.24)
xi 2
b=
∑ wi ( xi ) yi = m=
3.78
= 1.261;
∑ wi ( xi )2 3
The “hat” on y-values represents the estimated mean test capacity and is obtained
from equation (5.21) as yˆ = m x0 where m is the slope. Therefore the coefficient of
s m 0.093
variation v m = = =0.0737. For a geotechnical design capacity value x0 =
m 1.26
200kN, the estimated mean test capacity value is obtained from equation (5.22) in
which m =1.26 and the estimated mean test capacity value is 252 kN. Equation
(5.24) provides
1 0 .5
yp = m x0 - tα ,ν sm x0 [ 1 + ]
n
1
Therefore, for x0 = 400kN, y p = {1.26 * 400- 1.90 * 0.0735* [1+ ]0.5 * 400} =
3
422.78kN
700
600
500
400 Predicted Capacity
(90% Confidence
300
Level)
200
100
0
0.0 100.0 200.0 300.0 400.0 500.0 600.0 700.0 800.0 900.0 1000.0
Geotechnical Design Capacity (kN)
tα,ν=1.90from
1.2
Foundation Strength
Constant φF
"Case C"
0.6
0.4 "Case D"
0.2 Variable φF
0
0 200 400 600 800 1000
Geotechnical Design Capacity (kN)
700
600
Test Capacity (kN)
500
400
φF =1.084
300 0
45
200
degree
100 line
0
0 100 200 300 400 500 600
Geotechnical Design Capacity (kN)
Figure 9.9 – Foundation Strength factor (Buckley, 1994 & CIGRE, 2002)