Figure 2.1 Top View of The Excavation Pit: 5.0 M 5.0 M 5.0 M 5.0 M
Figure 2.1 Top View of The Excavation Pit: 5.0 M 5.0 M 5.0 M 5.0 M
Figure 2.1 Top View of The Excavation Pit: 5.0 M 5.0 M 5.0 M 5.0 M
2 EXCAVATION IN SAND
This tutorial describes the construction of an excavation pit in soft clay and sand layers.
The pit is a relatively small excavation of 12 by 20 m, excavated to a depth of 6.5 m
below the surface. Struts, walings and ground anchors are used to prevent the pit from
collapsing. After the full excavation, an additional surface load is added on one side of
the pit.
4.0 m
50.0 m
(30 20) (50 20)
(34 19) (41 19)
80.0 m
The proposed geometry for this exercise is 80 m wide and 50 m long, as shown in Figure
2.1. The excavation pit is placed in the center of the geometry. Figure 2.2 shows a cross
section of the excavation pit with the soil layers. The clay layer is considered to be
impermeable.
Objectives:
• Using the Hardening Soil model
• Modelling of ground anchors
• Using interface features
• Defining over-consolidation ratio (OCR)
• Prestressing a ground anchor
• Changing water conditions
• Selection of stress points to generate stress/strain curves
• Viewing plastic points
z=0 Fill
z = -1
Sheet pile walls
z = -4 Sand
Sand
z = -20
Figure 2.2 Cross section of the excavation pit with the soil layers
2.1 GEOMETRY
Project properties
• Start a new project.
• Enter an appropriate title for the project.
• Define the limits for the soil contour as xmin = 0, xmax = 80, ymin = 0 and ymax = 50.
• In the Interfaces tabsheet, select Manual in the Strength box and enter a value of
0.65 for the parameter Rinter . This parameter relates the strength of the interfaces to
the strength of the soil, according to the equations:
ci = Rinter csoil and tanϕi = Rinter tanϕi ≤ tanϕsoil
Hence, using the entered Rinter -value gives a reduced interface friction and interface
cohesion (adhesion) compared to the friction angle and the cohesion in the adjacent
soil.
Hint: When the Rigid option is selected in the Strength drop-down, the interface
has the same strength properties as the soil (Rinter = 1.0).
» Note that a value of Rinter < 1.0, reduces the strength as well as the stiffness
of the interface (Section 6.1.6 of the Reference Manual).
• In the Initial tabsheet, define the OCR -value according to Table 2.1.
• Click OK to close the window.
• In the same way, define the material properties of the "Sand" and "Soft Clay"
materials as given by Table 2.1.
• After closing the Material sets window, click the OK button to close the Modify soil
layers window.
• In the Soil mode right click on the upper soil layer. In the appearing right hand
mouse button menu, select the Fill option in the Set material menu.
• In the same way assign the Soft Clay material to the soil layer between y = -9.5 m
and y = -11.0 m.
• Assign the Sand material to the remaining two soil layers.
• Proceed to the Structures mode to define the structural elements.
Hint: The Tension cut-off option is activated by default at a value of 0 kN/m2 . This
option is found in the Advanced options on the Parameters tabsheet of the
Soil window. Here the Tension cut-off value can be changed or the option
can be deactivated entirely.
Table 2.4 Material properties for the embedded beams (grout body)
Parameter Name Grout Unit
Young's modulus E 3 · 107 kN/m2
Unit weight γ 24 kN/m3
Beam type − Predefined −
Predefined beam type − Massive circular beam −
Diameter Diameter 0.14 m
Axial skin resistance Type Linear −
Skin resistance at the top of the Tskin,start,max 200 kN/m
embedded beam
Skin resistance at the bottom of the Tskin,end,max 0.0 kN/m
embedded beam
Base resistance Fmax 0.0 kN
Hint: The colour indicating the material set assigned to the entities in project can
be changed by clicking on the Colour box of the selected material set and
selecting a colour from the Colour part of the window.
The remaining grouted anchors will be created by copying the defined grouted anchor.
Click on the Select button and click on all the elements composing both of the
ground anchors keeping the Ctrl key pressed.
Use the Create array function to copy both ground anchors (2 embedded beams + 2
node-to-node anchors) into a total of 4 complete ground anchors located at y = 24
and y = 28 by selection the 1D, in y direction option in the Shape drop-down menu
and define the Distance between columns as 4 m.
Multi-select all parts of the ground anchors (8 entities in total). While all parts are
selected and the Ctrl key is pressed, click the right mouse button and select the
Group from the appearing menu.
In the Model explorer tree, expand the Groups subtree by clicking on the (+) in front
of the groups.
• Click the Group_1 and rename it to "GroundAnchors".
Hint: The name of the entities in the project should not contain any space or
special character except "_" .
To define the sheet pile walls and the corresponding interfaces, follow these steps:
Select all four vertical surfaces created as the volume was decomposed. Keeping
the Ctrl key pressed, click the right mouse button and select the Create plate option
from the appearing menu.
Create a data set for the sheet pile walls (plates) according to Table 2.5. Assign the
data sets to the four walls.
• As all the surfaces are selected, assign both positive and negative interfaces to
them using the options in the right mouse button menu.
Hint: The term 'positive' or 'negative' for interfaces has no physical meaning. It
only enables distinguishing between interfaces at each side of a surface.
Non-isotropic (different stiffnesses in two directions) sheet pile walls are defined. The
local axis should point in the correct direction (which defines which is the 'stiff' or the 'soft'
direction). As the vertical direction is generally the stiffest direction in sheet pile walls,
Hint: The first local axis is indicated by a red arrow, the second local axis is
indicated by a green arrow and the third axis is indicated by a blue arrow.
More information related to the local axes of plates is given in the Reference
Manual.
Hint: The Enhanced mesh refinements are automatically used in mesh generation.
More information is available in Section 7.1.3 of Reference Manual.
The calculation consists of 6 phases. The initial phase consists of the generation of the
initial stresses using the K0 procedure. The next phase consists of the installation of the
sheet piles and a first excavation. Then the walings and struts will be installed. In phase
3, the ground anchors will be activated and prestressed. Further excavation will be
performed in the phase after that. The last phase will be the application of the additional
load next to the pit.
• Click on the Staged construction tab to proceed with definition of the calculation
phases.
• The initial phase has already been introduced. Keep its calculation type as K0
procedure. Make sure all the soil volumes are active and all the structural elements
are inactive.
Add a new phase (Phase_1). The default values of the parameters will be used for
Add another phase (Phase_4). The default values of the parameters will be used for
this calculation phase.
Select the soil volume to be excavated in this phase (between z = −1 and
z = −6.5).
In the Selection explorer under WaterConditions feature, click on the Conditions and
select the Dry option from the drop-down menu.
Display the contour lines for steady pore pressure distribution. Make sure that the
Legend option is checked in View menu. The steady state pore pressure distribution
is displayed in Figure 2.5. Scroll the wheel button of the mouse to zoom in or out to
get a better view.
• Change the legend settings to:
• Scaling: manual
• Maximum value: 0
• Number of intervals: 18
Figure 2.5 Preview of the steady state pore pressures in Phase_4 in a cross section
Hint: The visualization settings can be changed from the menu. For more
information refer Section 8.5.2 of Reference Manual .
• Click the Search closest button. The number of the closest node and stress point
will be displayed.
• Click the checkbox in front of the stress point to be selected. The selected stress
point will be shown in the list.
• Select also stress points near the coordinates (37.5 19 -5), (37.5 19 -6) and (37.5 19
-7) and close the Select points window.
• Click the Update button to close the Output program.
Start the calculation process.
Save the project when the calculation is finished.
Hint: Instead of selecting nodes or stress points for curves before starting the
calculation, points can also be selected after the calculation when viewing
the output results. However, the curves will be less accurate since only the
results of the saved calculation steps will be considered.
» To plot curves of structural forces, nodes can only be selected after the
calculation.
» Nodes or stress points can be selected by just clicking them. When moving
the mouse, the exact coordinates of the position are given in the cursor
location indicator bar at the bottom of the window.
After the calculations, the results of the excavation can be viewed by selecting a
calculation phase from the Phases tree and pressing the View calculation results button.
Select the final calculation phase (Phase_5) and click the View calculation results
button. The Output program will open and will show the deformed mesh at the end
of the last phase.
• The stresses, deformations and three-dimensional geometry can be viewed by
selecting the desired output from the corresponding menus. For example, choose
Plastic points from the Stresses menu to investigate the plastic points in the model.
• In the Plastic points window, Figure 2.6, select all the options except the Elastic
points and the Show only inaccurate points options. Figure 2.7 shows the plastic
points generated in the model at the end of the final calculation phase.
The graph will now show the major principal strain against the major principal stress.
Both values are zero at the beginning of the initial conditions. After generation of the
initial conditions, the principal strain is still zero whereas the principal stress is not zero
anymore. To plot the curves of all selected stress points in one graph, follow these steps:
• Select Add curve → From current project from right mouse button menu.
• Generate curves for the three remaining stress nodes in the same way.
The graph will now show the stress-strain curves of all four stress points (Figure 2.9). To
see information about the markers, make sure the Value indication option is selected
from the View menu and hold the mouse on a marker for a while. Information about the
coordinates in the graph, the number of the point in the graph, the number of the phase
and the number of the step is given. Especially the lower stress points show a
considerable increase in the stress when the load is applied in the last phase.
To create a stress path plot for stress node (37.5 19 -1.5) follow these steps:
Hint: To re-enter the Curve generation window (in the case of a mistake, a desired
regeneration or a modification), the Curve settings option from the Format
menu can be selected. As a result the Curves settings window appears, on
which the Regenerate button should be clicked.
» The Chart settings option in the Format menu may be used to modify the
settings of the chart.
Figure 2.10 Vertical effective stress (σ 'zz ) versus horizontal effective stress (σ 'yy ) at stress node
located near (37.5 19 -1.5)