Settle3D Liquefaction Theory Manual

© Rocscience Inc., 2016

Table of Contents

1 Introduction ............................................................................................................ 4
2 Theory ..................................................................................................................... 4
3 Cyclic Stress Ratio (CSR) ....................................................................................... 5
4 Stress Reduction Factor, rd .................................................................................... 5
5 Magnitude Scaling Factor, MSF ............................................................................ 8
6 Standard Penetration Test (SPT) Based Calculations........................................ 10
6.1 Pre-Defined Triggering Methods ................................................................... 10
6.2 SPT-N Value Correction Factors ................................................................... 11
6.2.1 Overburden Correction Factor, CN .......................................................... 12
6.2.2 Hammer Energy Efficiency Correction Factor, CE ................................. 14
6.2.3 Borehole Diameter Correction Factor, CB ............................................... 14
6.2.4 Rod Length Correction Factor, CR ........................................................... 15
6.2.5 Sampler Correction Factor, CS ................................................................ 16
6.3 Cyclic Resistance Ratio (CRR) ....................................................................... 16
6.3.1 Seed et al. (1984) ...................................................................................... 17
6.3.2 NCEER (1997) .......................................................................................... 17
6.3.3 Idriss and Boulanger (2004) .................................................................... 19
6.3.4 Cetin et al. (2004) – Deterministic .......................................................... 19
6.3.5 Japanese Bridge Code (JRA 1990) .......................................................... 19
6.3.6 Cetin et al. (2004) – Probabilistic ............................................................ 20
6.3.7 Liao et al. (1988) – Probabilistic .............................................................. 20
6.3.8 Youd and Noble (2001) – Probabilistic .................................................... 21
6.4 Relative Density, DR ....................................................................................... 22
6.5 Fines Content Correction ............................................................................... 23
6.6 Overburden Correction Factor, K ................................................................ 24
6.7 Shear Stress Correction Factor, K ............................................................... 26
7 Cone Penetration Test (CPT) Based Calculations .............................................. 27
7.1 Robertson and Wride (1997) ........................................................................... 27
7.1.1 Calculating IC ........................................................................................... 27
7.2 Modified Robertson and Wride (1998) ........................................................... 29
7.2.1 Calculating IC ........................................................................................... 29
 7.3 Idriss and Boulanger (2004) ........................................................................... 30
7.3.1 Calculating qc1N ........................................................................................ 30
7.4 Idriss and Boulanger (2014) ........................................................................... 31
7.5 Moss et al. (2006) – Deterministic ................................................................. 32
7.6 Moss et al. (2006) – Probabilistic ................................................................... 33
8 Shear Wave Velocity (Vs) Based Calculations ..................................................... 33
8.1 Andrus (2004) ................................................................................................. 33
8.2 NCEER (1997) ................................................................................................ 34
8.3 Juang et al. (2001) Probabilistic .................................................................... 34
9 Post-Liquefaction Lateral Displacement ............................................................. 36

Methods ........................................................................................ 36
9.1 Ground Profile ................................................................................................ 36
9.2 SPT
9.2.1 Zhang, Robertson, and Brachman (2004)................................................ 36
9.2.2 Tokimatsu and Yoshimi (1983)................................................................ 38
9.2.3 Shamoto et al. (1998) ............................................................................... 39
9.2.4 Wu et al. (2003) ........................................................................................ 42
9.2.5 Cetin et al. (2009) ..................................................................................... 42
9.3 CPT Methods ........................................................................................ 43
9.3.1 Zhang, Robertson, and Brachman (2004)................................................ 43
9.3.2 Yoshimine et al. (2006) ............................................................................ 44
9.4 VST Methods ........................................................................................ 45
10 Post-Liquefaction Reconsolidation Settlement ................................................. 46
10.1 SPT Methods........................................................................................... 46
10.1.1 Ishihara and Yoshimine (1992) ............................................................ 46
10.1.2 Tokimatsu and Seed (1984) .................................................................. 48
10.1.3 Shamoto (1984) ..................................................................................... 49
10.1.4 Wu et al. (2003) ..................................................................................... 52
10.1.5 Cetin et al. (2009).................................................................................. 52
10.2 CPT Methods .......................................................................................... 53
10.3 VST Methods .......................................................................................... 54
11 References .......................................................................................................... 55
1 Introduction

Settle3D offers different methods of calculating the factor of safety associated with
liquefaction resistance, probability of liquefaction, and the input parameters
required for those calculations. This manual also describes the calculating of lateral
spreading displacement as well as the vertical settlement due to liquefaction.

2 Theory

The use of in situ “index” testing is the dominant approach for assessment of the
likelihood of “triggering” or initiation of liquefaction. The methods available in
Settle3D are:

- Standard Penetration Test (SPT)

- Cone Penetration Test (CPT)
- Shear Wave Velocity (VST)

The potential for liquefaction can be evaluated by comparing the earthquake

loading (CSR) with the liquefaction resistance (CRR), expressed as a factor of safety
against liquefaction:

=
.

where

. = cyclic resistance ratio for an earthquake with magnitude 7.5

CSR = cyclic stress ratio
= magnitude scaling factor
= overburden stress correction factor
= ground slope correction factor
3 Cyclic Stress Ratio (CSR)

The cyclic stress ratio, CSR, as proposed by Seed and Idriss (1971), is defined as the
average cyclic shear stress, , developed on the horizontal surface of soil layers
due to vertically propagating shear waves normalized by the initial vertical effective
stress, , to incorporate the increase in shear strength due to increase in effective
stress. By appropriately weighting the individual stress cycles based on laboratory
test data, it has been found that a reasonable amplitude to use for the “average” or
equivalent uniform stress, , is about 65% of the maximum shear stress.

τ
= = 0.65 $ $ %&
! "
#
2

where

! "
#
= maximum horizontal ground surface acceleration (g)
= gravitational acceleration
= total overburden pressure at depth z

%&
= effective overburden pressure at depth z
= stress reduction factor

This equation is used to calculate CSR for all three analysis types.

4 Stress Reduction Factor, rd

The stress reduction factor, rd, is used to determine the maximum shear stress at
different depths in the soil. Values generally range 1 at the ground surface to lower
values at larger depths.

The SPT, CPT, and VST methods use the same rd formulations. The following are
provided in Settle3D:

- NCEER (1997)
- Idriss (1999)
- Kayen (1992)
- Cetin et al. (2004)
- Liao and Whitman (1986b)

NCEER (1997)

%& = 1.0 − 0.00765* +,% * ≤ 9.15/ 3

Idriss (1999)

ln2%& 3 = 42*3 + 62*3 7

*
42*3 = −1.012 − 1.126 sin ; + 5.133=
11.73
*
62*3 = 0.106 + 0.118 sin ; + 5.142=
11.28

* depth in meters ≤ 34/

where
=
7 = earthquake magnitude

For depths greater than 34m, rd = 0.5.

Kayen (1992)

%& = 1 − 0.012*
5

where

* = depth in meters

Cetin et al. (2004)

−23.013 − 2.949 ! " + 0.999 7 + 0.0525BC,DE!

∗
H1 + V
16.258 + 0.201I J.KLD@M&NJ.J O PQ,RST N . OUG
∗
%& @*, 7, ! " , BC,DE! G
∗
= ± XYZ
−23.013 − 2.949 ! " + 0.999 7 + 0.0525BC,DE!
∗
H1 + V
16.258 + 0.201I [email protected] O PQ,RST N . OUG
∗

for z<20 m (65ft)

%& @*, 7, ! " , BC,DE! G

∗

−23.013 − 2.949 ! " + 0.999 7 + 0.0525BC,DE!

∗
H1 + V
16.258 + 0.201I [email protected] O PQ,RST N . OUG
∗
= − 0.00462* − 203
−23.013 − 2.949 ! + 0.999 7 + 0.0525BC,DE!
∗
H1 + V
"
16.258 + 0.201I [email protected] O PQ,RST N . OUG
∗

± X YZ

for z≥20m (65ft)

 XYZ 2*3 = * J.O × 0.0198 for z<12m (40ft)
JJ

XYZ 2*3 = 12J.O × 0.0198 for z≥12m (40ft)

JJ
6

where

σ]^
*
_
= standard deviation (assumed to be zero)
= depth in meters
! "
BC,DE!
∗
= gravitational acceleration
= site shear wave velocity over the top 12m

- If the site stiffness estimation is difficult, take BC,DE!

∗
` 150-200m/s.
Notes:

- For very soft sites with BC,DE! less than 120m/s, use a limiting stiffness of
∗

- For very stiff sites, BC,DE!

∗
120m/s in calculations.
with stiffness greater than 250m/s, use 250m/s as
the limiting value in calculations.

Liao and Whitman (1986b)

%& = 1.0 − 0.00765* for z≤9.15m 7

%& = 1.174 − 0.0267* for 9.15 m<z≤ 23 m

where

* = depth below ground surface in meters

5 Magnitude Scaling Factor, MSF

If the magnitude of the earthquake is not 7.5, then the CRR values need to be
corrected for earthquake magnitude. The following corrections are available:

- Tokimatsu and Seed (1987)

- Idriss (1999)
- Idriss and Boulanger (2014) – SPT and CPT only
- Andrus and Stokoe (1997)
- Youd and Noble (1997) – SPT only
- Cetin (2004)
- Idriss (NCEER)

Tokimatsu and Seed (1987)

= 2.5 − 0.2 8

Idriss (1999)

= 6.9 exp ;− L = − 0.058 ≤ 1.8

d
9

This method can also be found in Idriss and Boulanger 2004 and 2008.

Idriss and Boulanger (2014)

=1+2 − 13 8.64 exp − $ − 1.325$

! "
4
10

efDghC K
= 1.09 + ; = ≤ 2.2
! "
180

2jD 3UJhC
E
= 1.09 + i k ≤ 2.2
! "
31.5

Andrus and Stokoe (1997)

MK.K
= $
7
7.5
11
Youd and Noble (1997)

The summary of the 1996/1998 NCEER Workshop proceedings by Youd and Idriss
(2001) outlines various methods for calculating the MSF and provide
recommendations for engineering practice.
The following MSF values are for calculated probabilities of liquefaction, the
equation for which is also shown.

op
l,#mn2op 3 = ln $ = −7.0351 + 2.1738 − 0.26782jD 3UJhC + 3.0265 ln2 3
1 − op 7

10K.OD
+,% op < 20% = L. K
+,% 7 <7

10K. L
+,% op < 32% = L.KK
+,% 7 <7

10L.ED
+,% op < 50% = L.OD
+,% 7 < 7.75

10E.EL
+,% 7 ≥ 7.5 = E. U
7

12

Cetin (2004)

7.5 E.ED
= $
7
13

Idriss (from NCEER report)

10E.EL
= E. U
14
6 Standard Penetration Test (SPT) Based Calculations

This section summarizes the methods available for calculating liquefaction

resistance based on SPT data. The following are presented:

- SPT N-Value Correction Factors

- Cyclic Resistance Ratio (CRR)
- Relative Density (DR)
- Fines Content Correction
- Overburden Correction Factor
- Shear Stress Reduction Factor

SPT-based calculations can be carried out two ways in Settle3D:

1. Pre-defined Triggering Methods – Users choose one of four pre-defined

methods for calculating liquefaction. When one of the pre-defined options are
chosen, the correction factors and triggering method are automatically
selected according to the method and cannot be modified.
2. Customized Triggering Methods – Users can select any combination of
correction factors and triggering methods.

6.1 Pre-Defined Triggering Methods

The following pre-defined triggering methods are available in Settle3D:

1. Youd et al. (2001)

2. Idriss and Boulanger (2008)
3. Cetin et al. – Deterministic (2004)
4. Cetin et al. – Probabilistic (2004)

The table below outlines the options that are automatically selected when each pre-
defined triggering method is used.

Triggering Methods
Youd et al. Idriss and Cetin et al. Cetin et al.
(2001) Boulanger (2004) – (2004) –
(2008) Deterministic Probabilistic
Triggering Method NCEER Idriss and Cetin et al. Cetin et al.
(1997) Boulanger (2004) (2004)
(2004) Deterministic Probabilistic
Depth Depth Liao & Idriss and Liao & Liao &
Correction Correction Whitman Boulanger Whitman Whitman
(1986) (2004) (1986) (1986)
 Sampling Standard Standard Standard Standard
Method
Advanced MSF Idriss Idriss and None None
Settings (1999) Boulanger
(2008)
Stress Idriss Idriss Cetin et al. Cetin et al.
Reduction (1999) (1999) (2004) (2004)
Factor
Relative Skempton Idriss and Skempton Skempton
Density (1986) Boulanger (1986) (1986)
(2003)
Fines Youd et al. Idriss and Cetin et al. Cetin et al.
Content (2001) Boulanger (2004) (2004)
Correction (2008)
K sigma Hynes and Idriss and Cetin et al. Cetin et al.
Olsen Boulanger (2004) (2004)
(1999) (2008)
K alpha None None None None

6.2 SPT-N Value Correction Factors

Before the CRR can be calculated, the N values obtained from the SPT must be

borehole diameter, and hammer energy efficiency, resulting in a2jD 3UJ value. The
corrected for the following factors: overburden, rod length, non-standard sampler,

equation below illustrates the correction.

jUJ = j t u v w

15

2jD 3UJ = jUJ g 16

Table 1: Summary of Correction Factors for Field SPT-N Values

Factor Equipment Term Correction

Variable
Overburden
CN Section 6.2.1
Pressure
Donut hammer 0.5-1.0
Safety hammer 0.7-1.2
Energy Ratio CE
Automatic
0.8-1.3
hammer
65 mm -115 mm 1.0
Borehole Diameter 150 mm CB 1.05
200 mm 1.12
<3 m 0.75
3m–4m 0.80
Rod Length 4m-6m CR 0.85
6 m -10 m 0.95
10 m – 30 m 1.00
Standard Sampler 1.0
Sampling Method Sampler without CS
1.0-1.3
Liner

6.2.1 Overburden Correction Factor, CN

The overburden correction factor adjusts N values to the N1 value that would be
measured at the same depth if the effective overburden stress was 1 atm.

The following formulations are available:

- Liao and Whitman (1986a)

- Bazaraa (1967)
- Idriss and Boulanger (2004)
- Peck (1974)
- Kayen et al. (1992)

Liao and Whitman (1986a)

o
J.

g =i k
J
17
Bazaraa (1967)

4
= +,% ≤ 1.5
g
1+2 J
J

4
18
= +,% > 1.5
g
3.25 + 0.5 J
J

m` my z`+

g ≤ 2.0

Idriss and Boulanger (2004)

o J. OLMJ.J UO|2gR 3}~

g = $ ≤ 1.7
{

19

2jD 3UJ ≤ 46

Peck, Hansen and Thorburn (1974)

2000
g = 0.77 log i k
J
20

J m` my zo ≤ 282 z•

Kayen et al. (1992)

2.2
g = ≤ 1.7
1.2 + o {

21
6.2.2 Hammer Energy Efficiency Correction Factor, CE

The energy efficiency correction factor is calculated using the measured energy ratio
as follows.

‚ !
=
w
60
22

It varies from 0.5-1.3. The ranges are taken from Skempton (1986).

Hammer Type CE
Donut hammer 0.5-1.0
Safety hammer 0.7-1.2
Automatic hammer 0.8-1.3

More specifically,

Hammer Type CE
Automatic Trip 0.9-1.6
Europe Donut Free fall 1.0
China Donut Free Fall 1.0
China Donut Rope& Pulley 0.83
Japan Donut Free Fall 1.3
Japan Donut Rope& Pulley 1.12
United States Safety Rope& pulley 0.89
United States Donut Rope& pulley 0.72
United States Automatic Trip Rope&
1.25
pulley

6.2.3 Borehole Diameter Correction Factor, CB

The following table, from Skempton (1986) summarizes the borehole diameter
correction factors for various borehole diameters.

Borehole Diameter (mm) CB

65-115 1.0
150 1.05
200 1.15
6.2.4 Rod Length Correction Factor, CR

The rod length correction factor accounts for how energy transferred to the
sampling rods is affected by the rod length.

6.2.4.1 Youd et al. (2001)

The following table from Youd et al (2001) summarizes the rod correction factor for
various rod lengths. The rod length above the ground is added to the depth to
obtain the total rod length before choosing the appropriate correction factor.

Rod Length (m) CR

<3 0.75
3-4 0.80
4-6 0.85
6-10 0.95
10-30 1.00

6.2.4.2 Cetin et al. (2004)

The figure below illustrates the recommended CR values (rod length from point of
hammer impact to tip of sampler). Note that Cetin assumes a length of 1.2m for rod
protrusion, and this is added to the depth before the correction factor is calculated.

Figure 1: Recommended Cr Values

6.2.5 Sampler Correction Factor, CS

The sampler correction factor is applied in cases when the split spoon sampler has
room for liner rings, but those rings were not used.

For the standard sampler, with a liner, the correction is 1.0.

For samplers without liners, the correction factor Cs ranges from 1.0-1.3 (NCEER,
1997). The following CS values are implemented.

u = 1.1 jD,UJ ≤ 10
CS Condition Reference

jD,UJ
(Cetin et al, 2004)
= 1+ 10 ≤ jD,UJ ≤ 30
u
100
(Cetin et al, 2004)
u = 1.3 jD,UJ ≥ 30 (Cetin et al, 2004)

6.3 Cyclic Resistance Ratio (CRR)

The cyclic resistance ratio is the other term required to calculate the factor of safety
against liquefaction. The cyclic resistance ratio represents the maximum CSR at
which a given soil can resist liquefaction.

The equation for CRR, corrected for magnitude, is

= .
23

The following methods of calculating CRR are available:

- Seed et al. (1984)

- NCEER (1997)
- Idriss and Boulanger (2004)
- Cetin et al. (2004) Deterministic
- Japanese Bridge Code (JRA 1990)
- Cetin et al. (2004) Probabilistic
- Liao et al. (1988) Probabilistic
- Youd and Noble (2001) Probabilistic
6.3.1 Seed et al. (1984)

Figure 2: Liquefaction boundary curves -

Correlation of (N1)60 values and CRR (M=7.5) (Seed et al. (1984)

6.3.2 NCEER (1997)

The curves recommended by Youd and Idriss (2001) / NCEER (1997) are based on
the Seed et al. (1984) curves.
 Figure 3: Simplified Base Curve Recommended for Calculation of CRR from SPT
data, with Empirical Liquefaction Data (modified from Seed et al., (1985)

The equation implemented in Settle3D is:

1 2jD 3UJhC 50 1
= + + −
.
34 − 2jD 3UJhC 135 ƒ102jD 3UJhC + 45„ E 200
24
6.3.3 Idriss and Boulanger (2004)

Idriss and Boulanger (2004) recommend the following equation:

2jD 3UJhC 2jD 3UJhC 2jD 3UJhC 2jD 3UJhC

E K L
= exp † +i k −i k +i k − 2.8‡
d… . , …D
14.1 126 23.6 25.4
25

6.3.4 Cetin et al. (2004) – Deterministic

The following equation is used to calculate CRR for a given probability of

liquefaction. The correction for fines content is built into the equation.

2jD 3UJ 21 + 0.004 3 − 29.53 ‰y2 73 − 3.70 ‰y $ + 0.05 + 16.85 + 2.70ŠMD 2op 3
o
22jD 3UJ , 7, , , op 3 = exp ˆ ‹
13.32

26

6.3.5 Japanese Bridge Code (JRA 1990)

This method is based on both the equivalent clean sand N value as well as the
particle size distribution.

Note that in the equation below is in z#/•/E.

2jD 3UJhC 0.35

= 0.0882Ž + 0.255 log $+ K +,% 0.05// ≤• < 0.6//
d… . , …D
+ 0.7 •J J

2jD 3UJhC
= 0.0882Ž − 0.05 +,% 0.6// ≤ • < 2//
d… . , …D
+ 0.7 J

K = 0 +,% < 40%

K = 0.004 − 0.16 +,% ≥ 40%

27
6.3.6 Cetin et al. (2004) – Probabilistic

Similar to the deterministic method, the Cetin et al. (2004) Probabilistic method has
the fines content correction built into the PL formulation.

2jD 3UJ 21 + 0.004 3 − 13.32 ln@ •‘ G − 29.53 ln2 73 − 3.70 ln $ + 0.05 + 16.85
o
op @2jD 3UJ , •‘ , 7, , G = Φ “− ”
2.70

28

6.3.7 Liao et al. (1988) – Probabilistic

Figure 4: Probabilistic SPT-based liquefaction triggering (Liao et al. 1988)

6.3.8 Youd and Noble (2001) – Probabilistic

The Youd and Noble (2001) formulation is outlined below.

op
l,#mn2op 3 = ln $ = −7.0351 + 2.1738 − 0.26782jD 3UJhC + 3.0265 ln2 3
1 − op 7

29

Figure 5: Probabilistic SPT-based liquefaction triggering (Youd and Noble, 1997)

6.4 Relative Density, DR

The relative density, DR, of a soil is used in the calculation of the overburden
correction factor, CN. The following methods are available:

- Skempton (1986)
- Ishihara (1979)
- Tatsuoka et al. (1980)
- Idriss and Boulanger (2003)
- Ishihara, Yasuda, and Yokota (1981)

Skempton (1986)

jD,UJ = 41 ∗ •t E 30

Ishihara (1979)

•t = 0.9 ∗ @jD,UJ + 14 + 6.51 logDJ G 31

Tatsuoka et al. (1980)

•t = 0.9 ∗ @jD,UJ + 14 + 6.51 logDJ G 32

Idriss and Boulanger (2003)

jD,UJ
•t = Ž
46
33

Ishihara, Yasuda, and Yokota (1981)

0.5
•t = 0.0676|jD,UJ + 0.085 ‰,#DJ $
•J
34
6.5 Fines Content Correction

The following fines content correction methods are available:

- Idriss and Boulanger (2008)

- Youd et al. (2001)
- Cetin et al. (2004)

Idriss and Boulanger (2008)

2jD 3UJhC = 2jD 3UJ + Δ2jD 3UJ

9.7 15.7 E
Δ2jD 3UJ = exp i1.63 + − $ k
+ 0.01 + 0.01
35

Youd et al. (2001)

2jD 3UJhC = 4 + 62jD 3UJ

4 = 0 +,% ≤ 5%

190
4 = exp H1.76 − E
$V +,% 5% < < 35%

4 = 5.0 +,% ≥ 35%

6 = 1.0 +,% ≤ 5%

D.
6 = –0.99 + i k— +,% 5% < < 35%
1000

6 = 1.2 +,% ≥ 35%
36

Cetin et al. (2004)

2jD 3UJhC = 2jD 3UJ ˜™gwu

= 21 + 0.004 3 + 0.05 i k +,% 5% ≤ ≤ 35%

˜™gwu
jD,UJ
37
6.6 Overburden Correction Factor, K

In addition to magnitude, the CRR can be corrected for overburden. The CRR of
sand depends on the effective overburden stress; liquefaction resistance increases
with increasing confining stress.

There are three options available for SPT:

- Hynes and Olsen (1999) (NCEER)
- Idriss and Boulanger (2008)
- Cetin et al. (2004)

Hynes and Olsen (1999)

2šMD3
=i k
{
o

+ = 0.7 − 0.8 +,% 40% < %I‰ nm›I œIy`mn• < 60%
+ = 0.6 − 0.7 +,% 60% < %I‰ nm›I œIy`mn• < 80% 38

The parameter f is a function of site conditions, and the estimates below are
recommended conservative values for clean and silty sands and gravels.

Figure 6: Recommended curves for estimating for engineering practice (from

NCEER 1996 workshop)
Idriss and Boulanger (2008)

This method is essentially the same as the one found in Idriss and Boulanger
(2004), except that the limit for K is higher.

= 1− ln i k ≤ 1.1
{
o

1
=
@18.9 − 2.55|2jD 3UJ G
39

Cetin et al. (2004)

The following figure illustrates the recommended values.

Figure 7: values, shown with NCEER recommendations (for n=0.7 and DR<60%)
for comparison
6.7 Shear Stress Correction Factor, K

is the static shear stress correction factor, used to correct CRR values for the
effects of static shear stresses. The only option available in Settle3D for this factor
is from Idriss and Boulanger (2003).

Ÿt
= + ž exp $
•

= 1267 + 6364 E − 634 exp243 − 632 exp2−43

ž = expƒ−1.11 + 12.34 E + 1.31 ln24 + 0.00013„

• = 0.138 + 0.1264 + 2.524 K

1
Ÿt = − •t
100•
− ln ; o =

4 ≤ 0.35

−0.6 ≤ Ÿt ≤ 0
40

where

•t
•′
= relative density

= empirical constant which determines the value of •′ at which dilatancy is

= mean effective normal stress

suppressed and depends on the grain type (Q~10 for quarz and feldspar, 8 for

o
limestone, 7 for anthracite, and 5.5 for chalk; Settle3D uses 8)

4
= atmospheric pressure
= tan of slope angle.
7 Cone Penetration Test (CPT) Based Calculations

The following methods are available in Settle3D for determining triggering of

liquefaction:
- Robertson and Wride (1997)
- Modified Robertson and Wride (1998)
- Boulanger and Idriss (2004)
- Boulanger and Idriss (2014)
- Moss et al. (2006) – Deterministic
- Moss et al. (2006) – Probabilistic

As mentioned in previous section, the magnitude scaling factor (MSF) and stress
reduction factor (rd) equations are the same as for SPT. These equations can be
found in sections 4 and 5.

7.1 Robertson and Wride (1997)

The following methods are employed in the Robertson and Wride (1997) triggering
method:
1. Calculate Ic using the procedure outlined in the NCEER summary report.
2. Calculate qc1N using the n value from the Ic calculation.

Depths with ehDghC ≥ 160 are considered not liquefiable.

3. Calculate qc1Ncs, with Kc calculated based on the NCEER recommendation.

h = 1.0 +,% ¢h ≤ 1.64

h = −0.403¢hL + 5.581¢hK − 21.63¢hE + 33.75¢h − 17.88 +,% ¢h > 1.64

41

ehDghC = h ehDg 42

4. Calculate CRR based on Robertson and Wride (1997).

ehDghC
= 0.833 £ ¤ + 0.05 m+ ehDghC < 50
.
1000

ehDghC K
= 93 £ ¤ + 0.08 m+ 50 ≤ ehDghC < 160
.
1000
43

7.1.1 Calculating IC

The soil behavior type index, Ic, is calculated using the following equation:
¢h = ƒ23.47 − log2 33E + 21.22 + log2 33E „J. 44

where

+C
= H V ∗ 100%
eh − {
45

eh − o ¥
=H V– $ —
{
o {
46

The recommended procedure for calculating the soil behavior type index is iterative,
as outlined in the NCEER summary report (Robertson and Wride, 1997).

1. Assume n=1.0 and calculate Q using the following equation.

eh − o D.J
eh −
=H $– $ —=H V
{ {
o { {
47

Calculate Ic using the equation in the previous section.

2. If Ic > 2.6 (or the user-defined Ic,max), the soil is clayey and not susceptible to
liquefaction.

3. If Ic < 2.6 (or the user-defined Ic,max), recalculate Q using n = 0.5, and
recalculate Ic.

4. If Ic < 2.6 (or the user-defined Ic,max), the soil is non-plastic and granular. No
further calculation is required.

5. If Ic > 2.6 (or the user-defined Ic,max), the soil is probably silty. Calculate qc1N
using the equations below, with n = 0.7 in the equation for CQ.

eh
ehDg = $ ≤ 254
«
o
48

o
¥

« =i k ≤ 1.7
J
49

6. Calculate Ic using the qc1N value calculated in (5).

7.2 Modified Robertson and Wride (1998)

The following methods are employed in Modified Robertson and Wride (1998):
1. Calculate Ic using the procedure outlined in Robertson and Wride (1998).
2. Calculate qc1N using the n value from the Ic calculation.

Depths with ehDghC ≥ 160 are considered not liquefiable.

3. Calculate qc1Ncs, with Kc calculated based on Robertson and Wride (1998).

h = 0 +,% ≤ 5%

h = 0.02672 − 53 +,% 5 < < 35%

h = 0.8 +,% ≥ 35% 50

ΔehDg = ehDg
h
1− h
51

ehDghC = ehDg + ΔehDg 52

4. Calculate CRR based on Robertson and Wride (1997).

ehDghC
= 0.833 £ ¤ + 0.05 m+ ehDghC < 50
.
1000
53

ehDghC K
= 93 £ ¤ + 0.08 m+ 50 ≤ ehDghC < 160
.
1000
54

7.2.1 Calculating IC

The recommended procedure for calculating the soil behavior type index is iterative,
as outlined in Robertson and Wride (1998).

1. Assume n=1.0 and calculate Q and Ic as outlined in Section 7.1.1. If ¢h > 2.6
then the point is considered not liquefiable.

2. If ¢h ≤ 2.6, calculate qc1N using n=0.5, and recalculate Ic using qc1N.

3. If the recalculated ¢h ≤ 2.6, the value of Ic calculated with n=0.5 is used. If Ic

iterates around 2.6 depending on n, then use n=0.75 to calculate qc1N and Ic.
7.3 Idriss and Boulanger (2004)

The following methods are employed in the Idriss and Boulanger (2004) triggering
method:
1. Calculate qc1N according to Idriss and Boulanger (2004) iterative procedure.
2. Calculate Kc, based on Idriss and Boulanger (2004).

1
= ≤ 0.3; ehDg ≤ 211
37.3 − 8.272ehDg 3J.EUL
55

= 1− ln i k ≤ 1.0
{
o
56

3. Calculate qc1Ncs, based on Idriss and Boulanger (2004).

ehDg 9.7 15.7 E

ΔehDg = ;5.4 + = exp i1.63 + − $ k
16 + 0.01 + 0.01
57

ehDghC = ehDg + ΔehDg 58

4. Calculate CRR based on Idriss and Boulanger (2004).

ehDghC ehDghC E ehDghC K ehDghC L

= exp +; = −; = +; = − 3$
d… . , -® …D
540 67 80 114
¯

59

7.3.1 Calculating qc1N

The following iterative procedure is used to calculate qc1N:

1. Calculate qc1N using n=1.0.

2. Recalculate qc1N using the following equation for n:

y = 1.338 − 0.2492ehDg 3J.EUL 60

A total of 100 iterations are performed, after which the last calculated value of qc1Nis
used.
7.4 Idriss and Boulanger (2014)

The following methods are employed in the Idriss and Boulanger (2014) triggering
method:
1. Calculate qcN = qt/Pa.

e¹ = eh + ºE 21 − 3 61

where a is the cone area ratio.

e¹
ehg =
o
62

2. Calculate qc1Ncs according to Idriss and Boulanger (2008). This is an iterative

procedure, as outlined below.

o !
g = $ ≤ 1.7

/ = 1.338 − 0.2492ehDghC 3J.EUL 63

ehDg = g ehg 64

ehDg 9.7 15.7 E

ΔehDg = ;11.9 + = exp i1.63 − − $ k
14.6 +2 +2
65

ehDghC = ehDg + ΔehDg 66

3. Calculate according to Idriss and Boulanger (2014).

1
= ≤ 0.3
37.3 − 8.272ehDghC 3J.EUL
67

= 1− ln i k ≤ 1.1
o
68
 4. Calculate CRR based on Idriss and Boulanger (2014).

ehDghC ehDghC E ehDghC K ehDghC L

= exp +; = −; = +; = − 2.80$
d… . , -̄ …D ¹!
113 1000 140 137
69

7.5 Moss et al. (2006) – Deterministic

The following methods are employed in the deterministic Moss et al. (2006)
triggering method:
1. Calculate qc1, with c calculated according to the method outlined in Moss et
al. (2006).

šS
• = +D $
š
+K
70

+C
= $ ∙ 100
š
eh −

+D = ¼D ∙ eh S
"

+E = −@•D ∙ eh S + •K G
½

+K = ž`2log210 + eh 33 ¾R

where

¼D = 0.78; ¼E = −0.33; •D = −0.32; •E = −0.35; •K = 0.49; *D = 1.21

ehD = ‘ ∙ eh

o h
‘ = $ ≤ 1.7

2. Calculate CRR according to Moss et al. (2006), based on a 50% probability of

liquefaction.

= exp{Àeh,D
D.JL
+ eh,D @0.110 ∙ š G + @0.001 ∙ š G + •@1 + 0.850 ∙ š G − 0.848
∙ ln2 7 3 − 0.002 ∙ ln2 3 − 20.923 + 1.632 ∙ ΦMD 2op 3Á /7.177}
71
7.6 Moss et al. (2006) – Probabilistic

The following methods are employed in the probabilistic Moss et al. (2006)
triggering method:
1. Calculate qc1, with c calculated according to the method outlined in Moss et
al. (2006). The calculations are outlined in the section above.
2. Calculate PL according to Moss et al. (2006), based on the user-defined Factor
of Safety, or calculate CRR based on the user-defined probability of
liquefaction. The CRR calculation method is outlined above.

eh,D
D.JL
+ eh,D @0.110 ∙ š G + @0.001 ∙ š G + •@1 + 0.850 ∙ š G − 7.177 ∙ ln2 3
Å −0.848 ∙ ln2 7 3 − 0.002 ∙ ln2 3 − 20.923 È
op = Φ{− Ä Ç
1.632
à Æ
72

8 Shear Wave Velocity (Vs) Based Calculations

The magnitude scaling factor (MSF) and stress reduction factor (rd) equations are
the same as for CPT and SPT. These equations can be found in sections 4 and 5.

The following methods are available in Settle3D for determining triggering of

liquefaction based on shear wave input:
- Andrus (2004)
- NCEER (1997)
- Juang et al. (2001) Probabilistic

Before triggering, the input vs value is normalized to vs1 as follows:

o J.E
BCD = BC $
73

8.1 Andrus (2004)

The following methods are employed in the Andrus (2004) triggering method:
1. Calculate Vs1cs using the formulation for Kfc from Juang et al.
2. Calculate CRR according to Andrus (2004).
 BCDhC E 1 1
= 0.022 H V + 2.8 H − V
.
100 215 − BCDhC 215
74

You can also account for an overburden correction factor. The Idriss and Boulanger
(2004) equation is as follows:

= 1− ln i k ≤ 1.1
o

1
= ≤ 0.3
18.9 − 3.12BCDhC /1003D.É U

75

8.2 NCEER (1997)

The following methods are employed in the NCEER method:

1. Calculate CRR according to NCEER recommendations.

BCD E ž ž
= $ + −
100 BCDh − BCD BCDh
76

where = 0.03 and ž = 0.9, and

Vs1cs = 220 for FC < 5%

Vs1cs = 210 for FC < 35%
Vs1cs = 200 for all other FC values

2. Calculate Vs1cs according to Juang et al, and calculate if desired.

8.3 Juang et al. (2001) Probabilistic

The Juang et al. (2001) method is outlined below:

1. Calculate Vs1cs.

BCDhC = šh BCD
77

where
 šh = 1, for ≤ 5%
šh = 1 + Ê2 − 53 for 5 < < 35%
šh = 1 + 30Ê for ≥ 35%

BCD BCD E
Ê = 0.009 − 0.0109 $ + 0.0038 $
100 100

2. Calculate PL based on the user-defined Factor of Safety, or calculate CRR

based on the user-defined probability of liquefaction.

op
ln H V = 14.8967 − 0.0611BCDhC + 2.6418 ln2 3
1 − op
78

If the PL is calculated, no further calculations are performed. If FS is being

calculated based on the CRR, then can be calculated.
9 Post-Liquefaction Lateral Displacement

The post-liquefaction lateral spreading is calculated by integrating the maximum

shear strain values over depth.

l•¢ = ËJ TÍÎ Ì! ∙ œ*
¾
" 79

9.1 Ground Profile

Zhang et al. (2004) proposed a method for estimating liquefaction-induced lateral

displacements based on the ground slope and/or free face height and distance to a
free face.

For a gently sloping ground without a free face:

l• = 2 + 0.23 ∙ l•¢
+,% 0.2% < < 3.5%
80

l MJ.O
l• = 6 ∙ $ ∙ l•¢
Ï
+,% 4 < l/Ï < 40 81

9.2 SPT Methods

The following methods are available for calculating the maximum shear strain,
when SPT data is used:
- Zhang, Robertson, and Brachman (2004)
- Tokimatsu and Yoshimi (1983)
- Shamato et al. (1998)
- Wu et al. (1993)
- Cetin et al. (2009)

9.2.1 Zhang, Robertson, and Brachman (2004)

In this method, the relative density (Dr) is first calculated based on the method
selected by the user.

The curves shown in the figure below are interpolated to determine the correct
maximum shear strain.
Figure 8: Relationship between maximum cyclic shear strain and factor of safety for
different relative densities
9.2.2 Tokimatsu and Yoshimi (1983)

The curves shown below are interpolated to determine the correct maximum shear
strain.

Figure 9: Shear strain induced by earthquake shaking

9.2.3 Shamoto et al. (1998)

For this method, one of three graphs is used to interpolate the maximum shear
strain.

For FC < 10%, the graph below is used.

Figure 10: Relationship between normalized SPT-N value and shear strain potential
for clean sands
For FC < 20%, the graph below is used.

Figure 11: Relationship between normalized SPT-N value and shear strain potential
for the case of FC=10%
For FC > 20%, the graph below is used.

Figure 12: Relationship between normalized SPT-N value and shear strain potential
for the case of FC=20%
9.2.4 Wu et al. (2003)

The graphs below are interpolated to find the maximum shear strain.

Figure 13: Estimation of cyclically induced deviatoric strains

9.2.5 Cetin et al. (2009)

The steps for calculating the maximum shear strain according to Cetin et al. (2009)
are outlined below.
1. Calculate according to Hynes and Olsen (1999). The formula can be found
in Section 6.6.
2. Calculate the relative density, Dr, according to the method selected by the
user.
 3. Calculate Kmc.

!h = −3 × 10M ∙ •tE + 0.0048•t + 0.7222 82

4. Calculate CSRss201D1.

CC,EJ,D,Ð,D = ∙ ∙ !h
83

5. Calculate Ì! " .

−0.025jDUJhC + ln@ CC,EJ,D,Ð,D G + 2.613

Ì! =
"
0.004jDUJhC + 0.001
84

where

5 ≤ jDUJhC ≤ 40; and 0.05 ≤ CC,EJ,D,Ð,D ≤ 0.6

9.3 CPT Methods

The following methods are available for calculating the maximum shear strain,
when CPT data is used:
- Zhang, Robertson, and Brachman (2004)
- Yoshimine (2006)

9.3.1 Zhang, Robertson, and Brachman (2004)

The relative density is first calculated according to Tatsuoka et al. (1990).

•Ñ = −85 + 76 log2ehDg 3 85

where ehDg ≤ 200.

The graph below is then used to determine Ì! " .

Figure 14: Relationship between maximum cyclic shear strain and factor of safety for
different relative densities

9.3.2 Yoshimine et al. (2006)

The Yoshimine et al. (2006) method is based on and a limiting shear strain.

= −11.74 + 8.342ehDghC 3J.EUL − 1.3712ehDghC 3J. EO

86

where efD¥fu ≥ 69.

ÌÒÓ! = 1.85922.163 − 0.4782ehDghC 3J.EUL 3K ≥ 0 87

The maximum shear strain is calculated as follows.

Ì! " = ÌÒÓ! m+ <

1−
Ì! = min †ÌÒÓ! , 0.03522 − 3 $‡
"
−

88
9.4 VST Methods

The and ÌÒÓ! expressions from Yoshimine et al. (2006) and Idriss and Boulanger
(2008) were adapted for shear wave velocity by Yi (2010).

BCDhC D.É U
BCDhC K.É E
= 0.032 + 0.836 $ − 0.190 $
100 100
89

where BCDhC ≥ 150 m/s.

BCDhC M . K
ÌÒÓ! = min –0.5, 7.05 $ —≥0
100
90

Ì! " is calculated as follows:

Ì! " = ÌÒÓ! m+ <

1−
Ì! = min †ÌÒÓ! , 0.03522 − 3 $‡
"
−

91
10 Post-Liquefaction Reconsolidation Settlement

The post-liquefaction settlement is calculated by integrating the volumetric strain

values over depth.
¾TÍÎ
=Ô Õ ∙ œ*
J
92

10.1 SPT Methods

The following methods are available for calculating Õ when SPT data is used:
- Ishihara and Yoshimine (1992)
- Tokimatsu and Seed (1984)
- Shamato (1984)
- Wu et al. (2003)
- Cetin et al. (2009)

10.1.1 Ishihara and Yoshimine (1992)

The following formulation is used to calculate the volumetric strain:

Õ = 1.5 ∙ exp2−2.5•t 3 ∙ min20.08, Ì! " 3

93

where DR is calculated according to the method specified by the user, and Ì! " is
calculated according to Zhang, Robertson, and Brachman (2004).
Figure 15: Ishihara and Yoshimine (1992) method for predicting volumetric and
shear strain
10.1.2 Tokimatsu and Seed (1984)

The figure below is used to interpolate a value of Õ .

Figure 16: Relationship between CSR, N160, and volumetric strain

10.1.3 Shamoto (1984)

One of three graphs is used to find Õ .

For FC<10%:

Figure 17: Relationship between normalized SPT-N, dynamic shear stress ratio, and
volumetric strain for clean sands
For FC<20%:

Figure 18: Relationship between normalized SPT-N, dynamic shear stress ratio, and
volumetric strain for FC=10%
For other fine content values:

Figure 19: Relationship between normalized SPT-N, dynamic shear stress ratio, and
volumetric strain for FC=20%
10.1.4 Wu et al. (2003)

The following graph is used to find Õ .

Figure 20: Correlations between CSR, N160cs, and reconsolidation volumetric strain
(Wu et al., 2003)

10.1.5 Cetin et al. (2009)

The Cetin et al. (2009) method incorporates a depth factor. With the depth factor,
the contribution of layers to settlement at the surface decreases as the depth of the
layer increases, and beyond a certain depth (zcr) the settlement of an individual
layer cannot be traced at the ground level. It was determined that the threshold
depth is 18m.

The steps for calculating the maximum shear strain according to Cetin et al. (2009)
are outlined below:
1. Calculate according to Hynes and Olsen (1999).
2. Calculate relative density, Dr, according to the method selected by the user.
 3. Calculate Kmc, and CSRss,20,1,D,1.

!h = −3 × 10M ∙ •tE + 0.0048•t + 0.7222

CC,EJ,D,Ð,D = ∙ ∙ !h
94

4. Calculate the critical depth factor, DF.

*
• =1−
*hÑÓ¹Óh Ò
95

where *hÑÓ¹Óh Ò = 18 /.

5. Calculate Õ , corrected for depth.

780.416 ln@ CC,EJ,D,Ð,D G − jDUJhC + 2442.465

Õ = 1.879 ln i k + 5.583
J
636.613jDUJhC + 306.732
96

Õ =• ∙Õ J 97

where the following limits apply:

5 ≤ jDUJhC ≤ 40; 0.05 ≤ CC,EJ,D,Ð,D ≤ 0.60; 0% ≤ Õ ≤ 5%

Note that it is left to the user to determine the normalized settlement.

10.2 CPT Methods

When CPT input data is used, the strain is calculated according to Yoshimine et al.
(2006).

Õ = 1.5 ∙ exp22.551 − 1.1472ehDghC 3J.EUL 3 ∙ min20.08, Ì! " 3 98

where Ì! " is calculated using the Yoshimine et al. (2006) formulation.

10.3 VST Methods

Yi (2010) adapted Ishihara and Yoshimine (1992) for Vs data, and the following
formulation for reconsolidation strain is used.

BCDhC D.É U
Õ = 1.5 ∙ exp i−0.449 $ k ∙ min20.08, Ì! " 3
100
99
11 References

Andrus, R., Stokoe, K. H. (1997), "Liquefaction Resistance Based on Shear

Wave Velocity", Proceedings of NCEER Workshop on Evaluation of
Liquefaction Resistance of Soils.

Boulanger, R. W., (2003a). "Relating Kα to relative state parameter index." J.

Geotechnical and Geoenvironmental Eng., ASCE 129(8), 770–73.

Boulanger, R. W., and Idriss, I. M. (2004). "State normalization of penetration

resistances and the effect of overburden stress on liquefaction resistance." Proc.,
11th Intl. Conf. on Soil Dynamics and Earthquake Engineering, and 3rd Intl. Conf.
on Earthquake Geotechnical Engineering, Doolin et al., eds, Stallion Press, Vol. 2,
pp. 484-491.

Cetin, K. O., and Bilge, H. T. (2012) “Performance-based assessment of magnitude

(duration) scaling factors.” J. Geotech. Geoenviron. Eng.,138(3), 324–334.

Cetin, K. O., Bilge, H. T.,Wu, J., Kammerer, A. M., and Seed, R. B. (2009).
“Probabilistic models for cyclic straining of saturated clean sands.” J. Geotech.
Geoenviron. Eng., 135(3), 371–386.

Cetin, K. O., Bilge, H. T.,Wu, J., Kammerer, A. M., and Seed, R. B. (2009).
“Probabilistic Model for the Assessment of Cyclically Induced Reconsolidation
(Volumetric) Settlements.” J. Geotech. Geoenviron. Eng., 135(3), 387–398.

Cetin K.O., Seed R.B., Der Kiureghian A., Tokimatsu K., Harder L.F. Jr, Kayen
R.E., MossR.E.S. (2004), “SPT-based probabilistic and deterministic assessment of
seismic soil liquefaction potential”, Journal of Geotechnical and Geoenvironmental
Engineering, ASCE, 130(12), 1314-1340.

Hynes, M. E., and Olsen, R. S. (1999), “Influence of confining stress on liquefaction

resistance.” Proc., Int. Workshop on Phys. And Mech. Of Soil Liquefaction,
Balkema, Rotterdam, The Netherlands, 145-152.

Idriss I. M., 1999, "An update to the Seed-Idriss simplified procedure for evaluating
liquefaction potential in Proceedings, TRB Workshop on New Approaches to
Liquefaction, Publication No. FHWA-RD-99-165, Federal Highway Administration,
January.

I. M. Idriss and R. W. Boulanger, “"Estimating Kα for use in evaluating cyclic

resistance of sloping ground." Proc. 8th US Japan Workshop on Earthquake
Resistant Design of Lifeline Facilities and Countermeasures against Liquefaction,
Hamada, O’Rourke, and Bardet, eds., Report MCEER-03-0003, MCEER, SUNY
Buffalo, N.Y., 2003, 449-468.

Idriss IM, Boulanger RW., Semi-empirical procedures for evaluating liquefaction

potential during earthquakes. Proc., 11th International conference on soil dynamics
and earthquake engineering, and 3rd International conference on earthquake
geotechnical engineering, vol. 1. Stallion Press; 2004. p. 32–56.

Idriss, I. M., and Boulanger, R. W. (2008). Soil liquefaction during earthquakes.

Monograph MNO-12, Earthquake Engineering Research Institute, Oakland, CA,
261 pp.

Ishihara, K. (1977), "Simple Method of Analysis for Liquefaction of Sand Deposits

During Earthquakes", Soils and Foundations, Vol. 17, No. 3, September 1977, pp. 1-
17.

Ishihara, K., Shimuzu, K., and Yamada, Y. (1981), “Pore Water Pressures Measured
in Sand Deposits During an Earthquake”, Soils and Foundations, Vol. 21, No. 4, pp.
85-100.

Ishihara, K., and Yoshimine, M. _1992_. “Evaluation of settlements in sand deposits

following liquefaction during earthquakes.” Soils Found., 32(1), 173–188.
JRA (1990) , Specification for Highway Bridges: Part V- Seismic Design. Japan
Road Association, Tokyo.

Juang, C. H., Fang, S. Y., Khor, E. H. (2006) “First-Order Reliability Method for
Probabilistic Liquefaction Triggering Analysis Using CPT”, J. Geotech. Geoenviron.
Eng., 132(3), 337-350.

Kayen, R. E, Mitchell, J. K., Seed, R. B.’ Lodge, A., Nishio, S., and Coutinho, R.
(1992), "Evaluation of SPT-, CPT-, and shear wave-based methods for liquefaction
potential assessment using Loma Prieta data", Proc., 4th Japan-U.S. Workshop on
Earthquake-Resistant Des. Of Lifeline Fac. And Counterneasures for Soil
Liquefaction, Vol. 1, 177-204.

Liao, S. S. C., Veneziano, D., Whitman, R.V. (1988), "Regression Models for
Evaluating Liquefaction Probability", Journal of Geotechnical Engineering,
ASCE, Vol. 114, No. 4, pp. 389-409.

Liao, S.S.C. and Whitman, R.V. (1986a). "Overburden Correction Factors for SPT in
Sand" Journal of Geotechnical Engineering, Vol. 112, No. 3, p. 373 - 377.
Liao, S.S.C. and Whitman, R.V. (1986b). "Catalogue of A Liquefaction and Non-
Liquefaction Occurrences During Earthquakes" Research Report, Dept. of Civil
Engineering, M.I.T., Cambridge, MA.

Meyerhof, G. G., 1957. Discussion on research on determining the density of sands

by spoon penetration testing, in Proceedings, 4th International Conference on Soil
Mechanics and Foundation Engineering, London, Vol. 3, p.110.

Moss, R. S. E, Seed, R. B., KAyen, R. E., Stewart, J. P., Der Kiureghian A., Cetin, K.
O. (2006) “CPT-Based Probabilistic and Deterministic Assessment of In Situ
Seismic Soil Liquefaction Potential”, J. Geotech. Geoenviron. Eng., 132(8), 1032-
1051.

NCEER, 1997, "Proceedings of the NCEER Workshop on Evaluation of Liquefaction

Resistance of Soils", Edited by Youd, T. L., Idriss, I. M., Technical Report
No. NCEER-97-0022, December 31, 1997.

Peck, R B Hanson, W E & Thornburn, T H (1974) Foundation engineering Pub:

John Wiley, New York

Robertson, P. K., Wride (Fear), C. E.,(1998) “Evaluating cyclic liquefaction potential

using the cone penetration test”, Can. Geotech. J. 35: 442–459 (1998).

Seed, H. B., Idriss, I. M. (1971), “Simplified Procedure for Evaluating Soil

Liquefaction Potential”, Journal of the Soil Mechanics and Foundations Division,
ASCE, Vol. 97, No SM9, Proc. Paper 8371, September 1971, pp. 1249-1273.

Seed, H. B., Idriss, I. M. (1982), "Ground Motions and Soil Liquefaction

During Earthquakes", Earthquake Engineering Research Institute Monograph
Series.

Seed, H. B., Idriss, I. M., Arango, I. (1983), "Evaluation of Liquefaction

Potential Using Field Performance Data", Journal of Geotechnical Engineering,
ASCE, Vol. 109, No. 3, pp. 458-482.

Seed, H. B., Tokimatsu, K., Harder, L. F., Chung, R. M. (1984), "The Influence of
SPT Procedures in Soil Liquefaction Resistance Evaluations", Earthquake
Engineering Research Center Report No. UCB/EERC-84/15, University of
California at Berkeley, October, 1984.

Shamoto, Y., Zhang, J., and Tokimatsu, K. (1998). “New charts for predicting large
residual post-liquefaction ground deformation.” Soil Dyn. Earthquake Eng., 17_7–
8_, 427–438.
Skempton, A.W. 1986. Standard penetration test procedures and the effects in
sands of overburden pressure, relative density, particle size, ageing and
overconsolidation. Geotechnique 36(3): 425-447.

Tokimatsu, K., and Seed, H. B. _1984_. “Simplified procedures of the evaluation of

settlements in clean sands.” Rep. No. UCB/GT-84/16, Univ. of California, Berkeley,
Calif.

Tokimatsu, K., and Seed, H. B., 1987. Evaluation of settlements in sands due to
earthquake shaking, J. Geotechnical Eng., ASCE 113 (GT8), 861-78.

Wu, J., Seed, R. B., and Pestana, J. M. (2003). “Liquefaction triggering and post
liquefaction deformations of Monterey 030 sand under unidirectional cyclic simple
shear loading.” Geotechnical Engineering Research Rep. No. UCB/GE-2003/01,
Univ. of California, Berkeley,Calif.

Youd, T. L., Hansen, C. M., and Bartlett, S. F., 2002. Revised Multilinear regression
equations for prediction of lateral spread displacement, J. Geotechnical and
Geoenvironmental Eng. 128(12),1007-017.

Youd, T. L., Idriss, I. M., Andrus, R. D., Arango, I., Castro, G., Christian, J. T.,
Dobry, E., Finn, W. D. L., Harder Jr., L. F., Hynes, M. E., Ishihara, K., Koester, J.
169 P., Liao, S. S. C., Marcusson III, W. F., Martin, G. R., Mtchell, J. K., Moriwaki,
Y., Power, M. S., Robertson, P. K., Seed, R. B., and Stokoe II, K. H. (2001).
“Liquefaction resistance of soils: Summary report from the 1966 NCEER and 1998
NCEER/NSF workshops on evaluation of liquefaction resistance of soils” J.
Geotechnical and Geoenvironmental Eng., 124(10), 817-833.

Youd, T. L., Noble, S. K. (1997), "Liquefaction Criteria Based on Statistical

and Probabilistic Analyses", Proceedings of the NCEER Workshop on
Evaluation of Liquefaction Resistance of Soils, December 31, 1997, pp. 201-205.

G. Zhang; P. K. Robertson, M.ASCE; and R. W. I. Brachman (2004).“Estimating

Liquefaction-Induced Lateral Displacements Using the Standard Penetration Test
or Cone Penetration Test, J. Geotechnical and Geoenvironmental Eng. 130(8), 861-
871.