4 - Use of DeepFEM - Finite Element Analysis in DeepEX
4 - Use of DeepFEM - Finite Element Analysis in DeepEX
Issued: 27-Sep-2019
www.deepexcavation.com
www.deepex.com
DEEP EXCAVATION 1
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
Table of Contents
1. Introduction ........................................................................................................................................ 3
2. Applying an Analysis method to a DeepEX model .............................................................................. 4
2.1. Conventional Limit Equilibrium Analysis Method (LEM). ............................................................ 4
2.2. Non-Linear (Beam on elastoplastic foundations) Method (NL). .................................................. 5
2.3. Limit Equilibrium and Non-Linear Analysis Combination Method (LEM+NL). ............................. 5
2.4. Finite Element Analysis Method (FEM)........................................................................................ 6
2.5. Limit Equilibrium and Finite Element Combination Method (LEM+FEM). .................................. 6
3. Defining the Mesh Analysis Options ................................................................................................... 7
4. Defining Convergence Options ........................................................................................................... 9
5. Defining Soil Properties for FEM Analysis ......................................................................................... 10
5.1. General Soil Properties .............................................................................................................. 11
5.2. Soil Models and Behavior Properties ......................................................................................... 12
6. Defining Wall Friction to the Model.................................................................................................. 13
7. Common Reasons why a FEM Analysis might not Converge ............................................................ 14
8. Performing FEM Analysis .................................................................................................................. 17
9. Finite Element Analysis results in DeepEX ........................................................................................ 17
10. Example 1: Braced Excavation - LEM, Non-Linear and FEM Analysis Comparison ......................... 21
11. Example 2: Slope Surface close to a Building: Slope Stability and FEM Settlements ..................... 26
12. Example 3: Stepped Anchored Walls – FEM Analysis ..................................................................... 31
DEEP EXCAVATION 2
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
1. Introduction
The Finite Element Method (FEM) is a numerical method for solving engineering problems.
Over the last years, it is adopted by several structural and geotechnical design software
programs and it offers an advanced solution to complex soil retaining problems. FEM uses
mesh generation techniques for dividing a complex problem into smaller parts, called finite
elements, that can be described by simpler equations. Usually, a Finite Element software is
used in order to assemble the finite elements and simulate the entire model.
DeepEX software uses our Finite Element Analysis method engine, called DeepFEM. This
option is implemented in all software versions since 2019, and it can be activated as an
additional, optional module within the software, allowing us to model full soil-structure
interaction considering all construction stage effects. The following options are available:
• Elastic and inelastic soil constitutive laws with associative or non-associative Failure
Criteria (e.g. Mohr-Coulomb, Drucker-Prager, Soil Hardening).
• Drained and Undrained modelling of clays with excess pore pressures.
• External loads: Prescribed, surface Strip and Linear Surcharges, 3D Loads (Buildings,
Footings, Area Loads).
• Automated, high-quality mesh generation with quadratic triangular finite elements.
• Inelastic zero thickness interface elements automatically generated on soil to
structural element interfaces.
• Steady state seepage analysis where the user can define boundary conditions, drains
and wells.
• Expert decision making on the estimation of the FEM analysis parameters.
• Calculation of Surface Settlements.
• Deadman Wall Systems.
• Design of Anchored walls and Braced excavations.
• Analysis of Multiple walls.
With the Finite Element Analysis method activated, DeepEX offers a complete set of all
analysis methods in deep excavation. The purpose of this document is to explain how the FEM
analysis method can be defined and used in DeepEX.
DEEP EXCAVATION 3
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
The following sections present the available analysis methods in DeepEX software.
In Limit Equilibrium Analysis, the retaining wall is analyzed to provide moment and force
equilibrium, when possible. Support reactions are also calculated, typically by using the
tributary area method.
DEEP EXCAVATION 4
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
On the reloading part, every soil model has a linear reloading elasticity parameter. Such a
parameter should typically range from 2 to 4 times the loading elasticity value (with average
3). In excavations, the reloading elasticity parameter typically describes the remaining soil
below the excavation while the loading elasticity is mostly applicable for soil on the retained
side. In a non-linear analysis the excavation models reduced to a plane problem, in which a
unit wide slice of the wall is analyzed, as outlined in the Figure below. Therefore, DeepEX is
not suitable to model excavation geometries in which three-dimensional effects may play an
important role. In the modelling of the soil-wall interaction, the very simple yet popular
Winkler approach is adopted. The retaining wall is modelled by means of beam elements with
transversal bending stiffness EI; the soil is modelled by means of a double array of
independent elastoplastic springs; at each wall grid point, two opposite springs converge at
most.
In stepped excavations and deadman wall systems, a part of the passive pressures of the back
wall is transferred to the front wall as an active impact load when the LEM or the LEM+NL
methods are used. The magnitude of this additional load is affected by several parameters,
like the depth of the back wall (thus the passive pressures), the distance between the two
walls etc.
DEEP EXCAVATION 5
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
DeepFEM can be used within the DeepEX software interactive interface, to analyze composite
models like braced excavations with struts and rakers. Anchored walls. Deadman wall systems
and more. It can calculate all analysis results – soil and water pressures, wall moment shear
and displacement diagrams, support reactions, structural and geotechnical ratios, surface
settlements and more. The results can be presented in tables and graphically on the model
area for each stage.
NOTE: The additional FEM Analysis module is required for this option.
NOTE: The additional FEM Analysis module is required for this option.
DEEP EXCAVATION 6
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
In the dialog that appears (see Image above), the following options are available:
DEEP EXCAVATION 7
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
• Level of Mesh
This option defines the level of refinement of the generated mesh (size of finite elements on
the model area). The level of mesh can be Coarse, Medium or Fine. This option affects the
finite element convergence and accuracy of results. A Fine Mesh would produce more
accurate results, but the analysis takes longer.
With this option selected, further mesh refinement is imposed on the soil layer boundaries
and soil body surface of each stage. The length variable DX corresponds to the maximum finite
element dimension allowed on the aforementioned locations.
With this option the user can select the preferred algorithm for the solution of the algebraic
non-linear equation system formed at each stage (or sub-step of stage) by the fem method.
When exponential constitutive laws are used, the option of Newton with line-search is
selected by default. The available options are:
This option defines the Hyperbolic Hardening model that will be used in FEM. For more
information, see Section 5 below. The following options are available.
Approximate Solution – This option makes the analysis faster and easier to converge
Exact Theoretical – This option makes the analysis more accurate but slower and more
difficult to converge
DEEP EXCAVATION 8
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
DEEP EXCAVATION 9
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
In the dialog that appears (see Image above), the following options are available:
This option enables the Crisfield’s Formulation of the Arc-length method for the solution of
the nonlinear equation system. The use of the arc-length method could provide better
convergence in problems with snapping behavior (work softening behavior of the constitutive
laws).
By selecting this option, we can define a custom convergence to our model analysis. If not
selected, the software will use the default tolerance, which is 3%. Any value more than 3%
should be used with caution.
This option, when selected, imposes a tension cut-off on all soil constitutive laws. (i.e.
principal stresses are always in compression).
• Define the Initial substep increment and the substep relative tolerance
This option enables the use of sub-stepping in the solution of each stage. More specifically,
the unbalanced forces calculated upon the initiation of each stage, are proportionally
eliminated in multiple steps. This option aids on the easier convergence of a stage, as each
sub-step converges on the equilibrium path and scenarios of divergence from the path are
avoided. The user can select the initial step increment, and when necessary the software will
further decrease it. Additionally, the user can select the normalized unbalanced force
tolerance for each sub-step. The When exponential constitutive laws are used, this option is
selected by default.
DEEP EXCAVATION 10
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
NOTE: For all soil layers (included sands) it is recommended to use a minimum cohesion value
of 5 psf (or 0.25 kPa). This is extremely beneficial to the FEM model convergence. If the model
is not converging, we suggest to try to increase the sands cohesion value.
Image: Access Edit Soils dialog and Define General Soil Properties
DEEP EXCAVATION 11
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
Elastoplastic Model: Single stiffness model with linear elasticity in combination with a non-
associative Mohr-Coulomb plasticity rule.
Small Strain Hardening: This model is currently under development for DeepFEM. When this
option is selected to any soil, DeepEX uses the Exponential model for the FEM analysis.
Subgrade Modulus: This option does not exist in FEM analysis (it works in the Non-Linear
method with the soil springs). When this model is selected, DeepEX uses the Elastoplastic
Model in FEM analysis.
DEEP EXCAVATION 12
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
NOTE: A use of a high wall friction percentage is beneficial to the FEM convergence. Typically,
we should use 30 to 70% of the soil friction.
DEEP EXCAVATION 13
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
• Project Geometry
In DeepEX, we should avoid totally vertical soil cuts. When defining the soil surface points, we
have to make sure that we use slopes. In general, an unsupported slope surface angle cannot
exceed the soil friction angle. If the cut is totally vertical, or quite bigger than the soil friction
angle, the model might collapse before the analysis goes through.
• Construction Stages
In DeepEX, we recommend the use of multiple construction stages for all analysis methods.
This allows us to review the results in all construction stages and take important decisions
about how to effectively optimize our models (increase/decrease structural sections, move
support locations etc.).
Apart from model optimization reasons, staging is important for the convergence of our
advanced analysis methods (Non-Linear and Finite Element Analysis). In these methods, an
initial stage with no excavation on either wall side should be used (the soil surface points on
the 2 sides of each wall should have the exact same elevation coordinate Z). After the initial
stage, we recommend the use of strict staging, where only a basic action should be executed
(excavation, support installation, backfill operation etc.). Dewatering and external load
activation can be performed in any construction stage, along with the basic action.
DEEP EXCAVATION 14
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
• Wall Friction
As mentioned in Section 6 above, the use of high values of wall friction (30 to 70% of the soil
friction) are recommended (and in most cases required) for the convergence of the FEM
analysis.
DEEP EXCAVATION 15
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
If a model is not well embedded according to the Limit Equilibrium Analysis, it is possible to
not converge in Non-Linear and Finite Element Analysis. For this reason, it is recommended
to analyze the model with Limit Equilibrium Method, optimize it (soil properties, wall depth
etc.), and then assign and run a more advanced method (NL or FEM).
While performing LEM analysis, we have to check all construction stages and make sure that
the moment and shear diagrams close below the excavation. We also need to check the basal
stability and the wall embedment safety factors. Having adequate wall embedment FS is
beneficial to the convergence of the NL and FEM analysis.
DEEP EXCAVATION 16
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
DEEP EXCAVATION 17
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
Image: Wall Moment and Shear Diagrams – Support Reaction (FEM Analysis)
DEEP EXCAVATION 18
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
DEEP EXCAVATION 19
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
DEEP EXCAVATION 20
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
10. Example 1: Braced Excavation - LEM, Non-Linear and FEM Analysis Comparison
A. Introduction – Model Description – Soil Properties – Structural Sections
In this example, we will examine a case of a braced excavation system (sheet piles braced
with struts), supporting a 30 ft excavation. A 0.6 ksf external load (Strip surcharge, developed
for 30 ft, starting 2 ft behind the wall) will be used on the top, on one side of the excavation.
The water table is at El: -15 ft.
The model will be generated and analyzed with DeepEX – Shoring Design Software, using all
available analysis methods (Limit Equilibrium, Non-Linear and Finite Element Analysis).
The following tables present the soil properties, the wall section properties and the support
properties respectively.
Table: Soil Properties and Stratigraphy
Layer Dry Total Friction Cohesion Modulus
Top Soil Unit Unit Angle c’ of
Soil Model Exponent
Elevation Type Weight Weight Elasticity
(ft) (pcf) (pcf) (deg) (psf) (ksf)
Wall Depth 50 ft
Wall Section AZ 19
DEEP EXCAVATION 21
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
DEEP EXCAVATION 22
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
Image: Water Pressures method and Defined Wall Friction for each wall
Design Standards
In this model we have used steel sheets and steel struts. In the Design tab of the software,
we will select to use AISC 360-16 Allowable method for the calculation of the structural
capacities of all steel members.
DEEP EXCAVATION 23
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
DEEP EXCAVATION 24
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
Image: Analysis and Checking Summary Table – Critical Results for each Design Section
DEEP EXCAVATION 25
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
11. Example 2: Slope Surface close to a Building: Slope Stability and FEM Settlements
A. Introduction – Model Description – Soil Properties – Structural Sections
Many times, in urban areas, we have to excavate a steep slope close to existing buildings. In
this case, it is recommended to examine both the slope stability, and the surface settlements
that might appear below the building. This design example presents a case with a 45deg slope
surface. An existing building is located on the top of the bench.
Two soil layers are considered in this model. This example focuses on the slope stability
analysis of the slope surface. The following table presents the soil properties and the model
stratigraphy.
For simplicity in this example, a circular slope stability analysis with the Morgenstern-Price
(General Limit Equilibrium) method is presented. More complex failure surfaces, possibly
using an automatic search approach, should also be examined in real projects.
In the same model, we will perform Finite Element Analysis, so that we can examine the soil
stresses, displacements and settlements below the building.
On the model area of DeepEX, we can right click on the left side of the wall and select to add
a left bench surface. In the dialog that appears, we can define the bench slope angle (45 deg),
the surface elevation on the bench top (20 ft), and where the bench distance from the wall (1
ft).
DEEP EXCAVATION 26
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
Add a Building
A building can be added as an external load. In the Building Wizard dialog, we can specify the
building location, geometry (length, height, number of floors, number of columns pes axis),
as well as, the dimensions, thicknesses and materials of all building floors, columns, walls and
footings.
DEEP EXCAVATION 27
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
In the Analysis tab of the software, we select to use the Finite Element Analysis method for
this project. In the Mesh Options dialog, we will select to use the exact theoretical solution
for the soil hardening model. In this example we will not use a user-defined convergence
tolerance. Finally, we will not use a wall friction in this example.
DEEP EXCAVATION 28
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
In the Slope tab of DeepEX, we can select to perform slope stability analysis. In this example,
we will use a circular slope surface, so we will define a rectangle of possible slope centers
close to the slope, and we will also define the surface radius search limits. The slope stability
safety factor will be calculated with the General Limit Equilibrium Method (Morgenstern-
Price).
DEEP EXCAVATION 29
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
DeepEX can analyze and present results for all construction stages and created design
sections, allowing the designer to review the full project history and take important decisions
for optimizing a solution while ensuring adequate safety factors. All results can be presented
in Tables and graphically, on the model area. The following images illustrate some of the
calculated results, both from the slope stability analysis (critical slope surface and safety
factor), and from the finite element analysis (settlements below the building).
Image: Slope Stability Analysis Results (Critical Slope Surface – Slope Stability FS)
DEEP EXCAVATION 30
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
In this example, we will examine a case of a stepped excavation system (anchored sheet piles),
supporting a 60 ft excavation. The water table is at El: -20 ft.
The model will be generated and analyzed with DeepEX – Shoring Design Software, using the
Finite Element Analysis method. The following tables present the soil properties, the wall
section properties and the support properties respectively.
Table: Soil Properties and Stratigraphy
Layer Dry Total Friction Cohesion Modulus
Top Soil Unit Unit Angle c’ of
Soil Model Exponent
Elevation Type Weight Weight Elasticity
(ft) (pcf) (pcf) (deg) (psf) (ksf)
Wall Depth 40 ft 45 ft
Wall Section AZ 26 AZ 26
Level 2 Wall 1 -18 ft Tieback 5 Strands 0.5’’ @ 10’ Strands 270 ksi
Level 3 Wall 2 -38 Tieback 6 Strands 0.6’’ @ 10’ Strands 270 ksi
Level 4 Wall 2 -48 Tieback 6 Strands 0.6’’ @ 10’ Strands 270 ksi
DEEP EXCAVATION 31
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
DEEP EXCAVATION 32
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
Image: Water Pressures method and Defined Wall Friction for each wall
Design Standards
In this model we have used steel sheets and steel struts. In the Design tab of the software,
we will select to use AISC 360-16 Allowable method for the calculation of the structural
capacities of all steel members.
DEEP EXCAVATION 33
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
In the Analysis tab of the software, we select to use the Finite Element Analysis method for
this project. In the Mesh Options dialog, we will select to use the exact theoretical solution
for the soil hardening model. In this example we will not use a user-defined convergence
tolerance. Finally, we will not use a wall friction in this example.
DEEP EXCAVATION 34
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
DeepEX can analyze and present results for all construction stages and created design
sections, allowing the designer to review the full project history and take important decisions
for optimizing a solution while ensuring adequate safety factors. All results can be presented
in Tables and graphically, on the model area. The following images illustrate some of the
calculated results from the finite element analysis (surface settlements, support reactions,
wall moment and displacement diagrams).
DEEP EXCAVATION 35
DEEPFEM – FINITE ELEMENT ANALYSIS IN DEEPEX
DEEP EXCAVATION 36