GEO-SLOPE International Ltd, Calgary, Alberta, Canada www.geo-slope.

com

An Introductory QUAKE/W Example

1 Introduction
This is a fairly simple introductory example. The prime purpose is to demonstrate to new QUAKE/W
users how to get started, to introduce the usual type of data required to do a QUAKE/W analysis, and to
illustrate what can be done with QUAKE/W.
The problem is quite simple, so that it can be analyzed with the GeoStudio Student Edition.

2 Problem configuration
Figure 1 shows the problem configuration. Basically it is an earth dam founded on an 8-metre stratum of
fairly loose, silty soil. The embankment is 5-m high with 2h:1v side slopes. The dam retains a reservoir
with a full supply level (FSL) at El. 12 m. Other relevant dimensions are shown in the following sketch.
14

12 2
1
10 Under-drain

8
Elevation

-2
-2 0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44

Distance

Figure 1 Tutorial problem configuration

The dam has a granular under drain to control the seepage through the dam. The drain itself is actually not
included in the analysis, but is considered when defining the phreatic surface.
On the downstream side, the water table is at the ground surface.
For the analysis, the geometry can be represented by two GeoStudio regions; one region for the
foundation and the second region for the embankment.

3 Objective
The objective of the analysis is to investigate the response of the ground and the structure when it is
subjected to an earthquake, and to examine the possibility that there may be some generation of excess
pore-pressures, which in turn could lead to some liquefaction.

4 Earthquake
The structure will be subjected to earthquake shaking according to the time-history record shown in
Figure 2. The peak acceleration is 0.12 g and the duration is 10 seconds.

QUAKE/W Example File: QUAKE Tutorial.docx (pdf) (gsz) Page 1 of 10

 GEO-SLOPE International Ltd, Calgary, Alberta, Canada www.geo-slope.com

0.15

0.10
Acceleration ( g )

0.05

0.00

-0.05

-0.10
0 2 4 6 8 10

Time (sec)

Figure 2 Earthquake time-history record

5 Initial insitu stresses

The first step in any QUAKE/W analysis is to establish the insitu stress state conditions that exist before
the earthquake occurs. Stated another way, it is necessary establishing the initial static stresses before
subjecting the structure to the dynamic action.
In QUAKE/W the initial insitu static stress can be computed with the Initial Static analysis type.
The two most important soil properties required for the Initial Static analysis is the total unit weight of the
materials and Poisson’s ratio ν.
14

12 2
1
10 Under-drain

8
Elevation

-2
-2 0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44

Distance

Figure 3 Boundary conditions for the Initial Static analysis

QUAKE/W Example File: QUAKE Tutorial.docx (pdf) (gsz) Page 2 of 10

 GEO-SLOPE International Ltd, Calgary, Alberta, Canada www.geo-slope.com

Poisson’s ratio is important because it influences Ko - coefficient of earth pressure at rest. For a 2-D
analysis like we are doing here, Ko is equal to ν / (1-ν). Poisson’s ratio is specified as 0.334 (1/3) for this
problem which represents a Ko condition equal to 0.5.

5.1 Boundary conditions for insitu analysis

Along the vertical ends of the problem, we must specify the equivalent of rollers. This is done by
specifying a zero x-displacement condition along these edges. The ground is free to move in the vertical
direction, but is fixed in the horizontal direction.
The problem is fixed along the base. This is accomplished by specifying the displacement to be zero in
both the horizontal (x) and vertical (y) directions.
To compute the correct total and effective stresses, it is necessary to apply the weight of the reservoir
water as a boundary condition. The arrows perpendicular to the ground surface in Figure 3 graphically
represent the water pressure. The water weight or pressure is specified with a Fluid Pressure type of
boundary. With this type of boundary condition, it is necessary to specify the FSL elevation and the unit
weight of the fluid (water in this case).

5.2 Initial pore-pressure conditions

The initial pore-pressures are defined by and computed from a specified watertable. The watertable goes
along the reservoir water surface as it naturally should, and then follows a straight line in the embankment
from the FSL elevation to the upstream tip of the under drain, and then follows a horizontal line to the
right end of the problem at the downstream ground surface elevation.

5.3 Finite element mesh

Once the problem has been completely defined with respect to the geometry, the boundary conditions and
the material properties, it is necessary to check that the finite element mesh is appropriate for the analysis.
The mesh can be viewed and modified with the Draw Mesh Properties command. In this case, it was
specified that the global element size should be about one 1 m. With this one and only user-specified
criterion, the default mesh is as presented in Figure 4, which is deemed adequate for this particular
analysis.

Figure 4 Finite element mesh

QUAKE/W Example File: QUAKE Tutorial.docx (pdf) (gsz) Page 3 of 10

 GEO-SLOPE International Ltd, Calgary, Alberta, Canada www.geo-slope.com

5.4 Insitu results

There are many different ways to examine the resulting insitu stresses. One typical way is to look at the
total and effective vertical stress contours, as in Figure 5 and Figure 6. Notice how the total stress under
the reservoir where the ground is horizontal is about 40 kPa. The corresponding effective stress contour is
near zero, as it correctly should be.

20

60

100

140

180

Figure 5 Total vertical stress contours

10 10

30

50 70

90

Figure 6 Effective vertical stress contours

The unit weight of water has been specified as 10 kN/m3. This is purely for convenience. With a round number like
this, it easy to make mental calculation for checking that the results are meaningful and as intended.

QUAKE/W Example File: QUAKE Tutorial.docx (pdf) (gsz) Page 4 of 10

 GEO-SLOPE International Ltd, Calgary, Alberta, Canada www.geo-slope.com

6 Dynamic Analysis
Once the insitu static stresses have been established, the next step is to do the dynamic or shaking
analysis.
With the KeyIn Analyses command, a new QUAKE/W analysis is created where the “Parent” is the
previous Initial Static analysis. Both the initial stress conditions and the pore-pressure conditions come
from the Parent (previous) analysis.
The Equivalent Linear Dynamic analysis type is used here.

6.1 Boundary conditions

The boundary conditions at the vertical ends of the problem have to be changed for the dynamic analysis.
Now the vertical movement is fixed, but the ground is allowed to move laterally. Notice the different
boundary condition symbols along the ends in Figure 7. These conditions allow the ground to sway from
side to side when the horizontal earthquake accelerations are applied.

Figure 7 Boundary conditions for the dynamic analysis

6.2 History node

In QUAKE/W, selected points can be flagged where the results will be saved for each and every time step
while integrating through the earthquake record. They are known as History Points. Two History Points
have been defined here. They are the two small red squares shown in Figure 8.

Figure 8 Locations of History Points

QUAKE/W Example File: QUAKE Tutorial.docx (pdf) (gsz) Page 5 of 10

 GEO-SLOPE International Ltd, Calgary, Alberta, Canada www.geo-slope.com

6.3 Stiffness and damping properties

The constitutive behavior of the soil will be treated as being linear elastic.
The shear modulus Gmax is 5000 kPa for both the foundation and embankment soil.
The damping ratio is a constant 0.1 (10%).

6.4 Pore-pressure functions

To make compute excess pore-pressures that may arise due to the shaking, it is necessary as a minimum
to define two functions. They are called the Pore-Pressure Ratio function and the Cyclic Number
function. For this introductory example, we can use the sample functions shown in Figure 9 and Figure 10
(these sample functions are included with the software). The Cyclic Number function in Figure 10 is the
sample function for loose sand.

PWP dunction
1.0

0.8
PWP Ratio

0.6

0.4

0.2

0.0
0.0 0.2 0.4 0.6 0.8 1.0

Cyclic Number Ratio N/NL

Figure 9 Pore-pressure ratio function

Cyclic # function
1000

100
Cyclic Number

10

1
0.1 0.2 0.3 0.4

Shear Stress Ratio

Figure 10 Cyclic number function

QUAKE/W Example File: QUAKE Tutorial.docx (pdf) (gsz) Page 6 of 10

 GEO-SLOPE International Ltd, Calgary, Alberta, Canada www.geo-slope.com

The embankment material is deemed to be non-liquefiable soil, and so no pore-pressure functions are
specified for this material.

6.5 Time integration

QUAKE/W is formulated on the basis of a time integration scheme. This means that QUAKE/W steps
through the earthquake record at a specified time interval and does a finite element analysis for each time
step. Usually, the time steps are equal to the time interval of the data points in the earthquake record.
In this analysis, the earthquake record data points have a time interval of 0.02 seconds, making a total of
500 data points for the 10-second duration. This means there will be 500 time steps, or 500 finite element
analyses.
By default, QUAKE/W stores data at the 20 highest peaks in the earthquake record. In addition, it has
been specified that the data should be saved every 50th time step.

7 Results of the dynamic analysis

As noted earlier for the Static analysis, there are many ways of presenting and viewing the results, and
this is particularly true the dynamic part of the analysis. Only a sampling of what can be done and is
available will be presented here.

7.1 Response at dam crest

We had defined a History Point at the crest of the dam and at the base of the problem. With the
CONTOUR Graph command, we can create graphs specific to the History Points. It is possible, for
example, to create graphs of x-accelerations versus time for the crest and the base and plot them together
at the same time as in Figure 11. The base values are a replication of the input; that is, the input
earthquake record. Recall that the input record has a peak of 0.12 g. As is evident in Figure 11, the peak
horizontal accelerations at the crest reach values greater than 0.3 g, indicating there is some amplification
in the motion at the crest.

0.4

0.3

0.2
X-Acceleration (g)

0.1

0.0

-0.1

-0.2

-0.3

-0.4
0 2 4 6 8 10

Time (sec)

Figure 11 Horizontal accelerations at the crest compared with the base input

QUAKE/W Example File: QUAKE Tutorial.docx (pdf) (gsz) Page 7 of 10

 GEO-SLOPE International Ltd, Calgary, Alberta, Canada www.geo-slope.com

7.2 Lateral movement

Figure 12 and Figure 13 show lateral displacement profiles at the center-line of the dam. Figure 12 shows
the total or absolute displacements. Absolute displacements include the solid body motion at the base of
the problem, together with any additional movements relative to the base. Figure 13 shows the
displacements relative to the fixed base.

Lateral displacement
14

12

10
Y (m)

0
-0.02 0.00 0.02 0.04 0.06 0.08

X-Displacement (m)

Figure 12 Absolute lateral displacement profiles at the center line of the dam

Lateral displacement
14

12

10
Y (m)

0
-0.03 -0.02 -0.01 0.00 0.01 0.02 0.03

Relative X-Displacement (m)

Figure 13 Relative lateral displacement profiles at the center line of the dam

In a QUAKE/W analysis, it is very important to be mindful of the fact that only relative displacements create
dynamic shear stresses and only dynamic shear stresses lead to the generation of excess pore-pressures.

QUAKE/W Example File: QUAKE Tutorial.docx (pdf) (gsz) Page 8 of 10

 GEO-SLOPE International Ltd, Calgary, Alberta, Canada www.geo-slope.com

Figure 14 and Figure 15 show the relative displacements as a deformed mesh at two different times.
These two diagrams illustrate the swaying of the ground during the earthquake.

Figure 14 Relative displacements at 4.0 seconds into the shaking (100x exaggeration)

Figure 15 Relative displacements at 4.52 seconds into the shaking (100x exaggeration)

Click this link to watch a movie of the motion

7.3 Liquefaction
In an Equivalent Linear analysis, a key parameter obtained is the Cyclic Stress Ratio (CSR). This number
is used together with the Cyclic Number function discussed earlier (Figure 10) to indicate the possibility
of liquefaction. Without going into much detail, generally the higher the CSR, the higher the possibility of
liquefaction. In this example, CSR’s greater than about 0.2 will indicate liquefaction is possible, as
demonstrated in Figure 17.

0.3
0.3
0.2 0.15 0.25
0.
2

Figure 16 Cyclic Stress Ratio contours

QUAKE/W Example File: QUAKE Tutorial.docx (pdf) (gsz) Page 9 of 10

 GEO-SLOPE International Ltd, Calgary, Alberta, Canada www.geo-slope.com

0.3
0.3
0.2 0.15 0.25
0.
2

Figure 17 Zones of potential liquefaction

Another interesting piece of information is the generated excess pore-pressures, as shown in Figure 18.
Of interest are the high excess pore-pressures on the downstream side, where there is pocket of excess
pore-pressures that exceed 30 kPa.

10
20
10 25
15
30
20

25

Figure 18 Contours of excess pore-pressure

8 Concluding remarks
This introductory example is rather simplistic and not all that realistic in some respects, but it nonetheless
vividly illustrates the powerful features and capabilities available in QUAKE/W.

QUAKE/W Example File: QUAKE Tutorial.docx (pdf) (gsz) Page 10 of 10