3D2013 2 Reference
3D2013 2 Reference
Reference Manual
2013
Build 6448
TABLE OF CONTENTS
TABLE OF CONTENTS
1 Introduction 11
2 General information 13
2.1 Units and sign conventions 13
2.2 File handling 15
2.3 Help facilities 15
10 Curves 299
10.1 Selecting points for curves 299
10.1.1 Mesh point selection 299
10.1.2 Pre-calculation points 300
10.1.3 Post-calculation points 300
10.2 Generating curves 301
10.2.1 Load-displacement curves 302
10.2.2 Force-displacement curves 303
10.2.3 Displacement-time or force-time curves 304
11 References 315
Index
Appendix A - The format required for the files to import CPT data 321
1 INTRODUCTION
2 GENERAL INFORMATION
Before describing the specific features in the different parts of the PLAXIS 3D user
interface, the information given in this chapter applies to all parts of the program.
It is important in any analysis to adopt a consistent system of units. At the start of the
input of a geometry, a suitable set of basic units should be selected. The basic units
comprise a unit for length, force and time. These basic units are defined in the Model
tabsheet of the Project properties window in the Input program. The default units are
meters [m] for length, kiloNewton [kN] for force and day [day] for time. Table 2.1 gives an
overview of all available units, the [default] settings and conversion factors to the default
units. All subsequent input data should conform to the selected system of units and the
output data should be interpreted in terms of the same system. From the basic set of
units, as defined by the user, the appropriate unit for the input of a particular parameter is
generally listed directly behind the edit box or, when using input tables, above the input
column. In all of the examples given in the PLAXIS manuals, the standard units are used.
Table 2.1 Available units and their conversion factor to the default units
Length Conversion Force Conversion Time Conversion
mm = 0.001 m N = 0.001 kN s (sec) = 1/86400 day
cm = 0.01 m [kN] = 1 kN min = 1/1440 day
[m] =1m MN = 1000 kN h = 1/24 day
km = 1000 m lbf (pounds force) = 0.0044482 kN [day] = 1 day
in (inch) = 0.0254 m kip = 4.4482 kN
ft (feet) = 0.3048 m
yd (yard) = 0.9144 m
For convenience, the units of commonly used quantities in two different sets of units are
listed below:
Int. system Imperial system
Basic units: Length [m] [in] or [ft]
Force [N] [lbf] or [kip]
Time [day] [day]
Sign convention
Stresses computed in PLAXIS 3D are based on the Cartesian coordinate system shown
in Figure 2.1. In all of the output data, compressive stresses and forces, including pore
pressures, are taken to be negative, whereas tensile stresses and forces are taken to be
positive. Figure 2.1 shows the positive stress directions.
σzz
z
σzy
σzx σyz
σxz
σyy
y σxy σyx
x σxx
All file handling in PLAXIS is done using a modified version of the general Windows® file
requester (Figure 2.2).
With the file requester, it is possible to search for files in any admissible folder of the
computer (and network) environment. The main file used to store information for a
PLAXIS project has a structured format and is named <project>.P3D, where <project> is
the project title. Besides this file, additional data is stored in multiple files in the sub-folder
<project>.P3DAT. It is generally not necessary to enter such a folder because it is not
possible to read individual files in this folder.
To inform the user about the various program options and features, PLAXIS 3D provides
a link in the Help menu to a digital version of the Manuals. A link is provided to the
reference document of the commands used in the program as well. Moreover, the Help
menu may be used to generate a file with software license information as stored in the
security lock (to be used for license updates and extensions). A more detailed description
of the Help menu of the Input and Output program is given in Section 3.3.11 and Section
8.2.11 respectively.
Many features are available as buttons in a toolbar. When the mouse pointer is
positioned on a button for more than a second, a short description ('hint') appears,
indicating the function of the button. Feedback from the program on the user actions is
given in the command line area. For some input parameters side panels appear to help
the user decide which value to select.
To carry out a finite element analysis using the PLAXIS 3D program, the user has to
create a three dimensional geometry model composed of points, lines, surfaces,
volumes, and other components and specify the material properties and boundary
conditions. This is done in the first two tabsheets (Geometry modes) of the Input
program. The mesh generation and the definition of the calculation phases is done in the
last three tabsheets (Calculation modes) of the Input program.
This icon represents the Input program. At the start of the Input program the Quick
select window appears in which a choice must be made between the selection of
an existing project and the creation of a new project (Figure 3.1).
Project
The title, folder and the file name of the project are available in the Project group box
available in the Project tabsheet.
Title The defined title appears as a default name for the file of the
Comments
The Comments box in the Project tabsheet gives the possibility to add some extra
comments about the project.
Company logo
Clicking the Company logo box in the Project tabsheet will open a file requester from
which the desired file can be selected. This logo can be included in the output plots.
General options
The general options of the project are available in the Project tabsheet of the Project
properties window.
Units
Units for length, force and time to be used in the analysis need to be specified. These
basic units are entered in the Model tabsheet of the Project properties window (Figure
3.3).
The default units, as suggested by the program, are m (meter) for length, kN (kiloNewton)
for force and day for time. The corresponding units for stress and unit weights are listed
in the box below the basic units.
All input values should be given in a consistent set of units (Section 2.1). The appropriate
unit of a certain input value is usually given directly behind the edit box, based on the
basic set of units.
Gravity
By default, the earth gravity acceleration, g , is set to 9.810 m/s2 , assuming the default
basic length unit is [m]. The direction of gravity coincides with the negative z -axis.
Gravity is implicitly included in the unit weights given by the user (Section 6.1).
Contour
At the start of a new project the user needs to specify the contours of the geometry
model. The initial setting of the xmin , xmax , ymin and ymax parameters set the outer
horizontal boundaries of the geometry model. The initial view of the draw area is such
that the model dimensions are fully visible within the draw area. The aforementioned
parameters are entered in the Model tabsheet of the Project properties window. The
contour of the geometry can also be modified in the draw area when the Soil mode is
selected.
available folders and to select the desired PLAXIS project file (*.P3D). After the selection
of an existing project, the corresponding geometry is presented in the main window.
An existing PLAXIS 3D project can also be read by selecting the Open project option in
the File menu. In the file requester, the type of the file is, by default, set to 'PLAXIS 3D
files (*.P3D)'.
The options for compression and volume size are available in the Archive
options window (Figure 3.5), displayed by clicking the button in the Purpose box.
The Content box displays the options for the information to be included in the archive is
shown. The options available are:
The general layout of the Input program for a new project is shown in Figure 3.6. The
main window o the Input program contains the following items:
Title bar
The name of the program and the title of the project is displayed in the title bar. Unsaved
modifications in the project are indicated by a '∗' in the project name.
Menu bar
The menu bar contains drop-down menus covering the options available in the Input
program.
General toolbar
The general toolbar contains buttons for general actions such as disk operations, printing
and model display arrangement. It also contains a button to start the Output program.
Hint: If the mouse is moved over a button in a toolbar, a hint about the function of
this button is displayed.
Mode tabs
The mode tabs are used to separate different modelling modes. The following tabs are
available:
Soil The soil stratigraphy is defined.
Structures The structural components of the model are defined.
Mesh The mesh properties are defined and the geometry model is
discretized and transformed to a finite element model.
Water levels Water levels are defined.
Staged construction The calculation settings are defined and the project is calculated.
For a more detailed description of the modes in the Input program see Section 3.4.
Hint: The first two modes (Soil and Structures) are referred to as the Geometry
modes, whereas the latter three modes (Mesh, Water levels and Staged
construction) are referred to as the Calculation modes.
Draw area
The draw area is the drawing sheet on which the geometry model is created and
modified. The geometry model can be created by means of the mouse and using the
buttons available in the side toolbar. The buttons in the side toolbar vary according to the
active mode. A more detailed description on the draw area is given in Section 3.5.
Side toolbar
The buttons of the features available in the active mode are located in the side toolbar. A
hint on the function of the button is displayed when the pointer is located on it for a few
seconds. The side toolbar also contains options to select one or more objects in the draw
area.
Status bar
The status bar displays information about the location of the mouse cursor in the draw
area, the model object it is located on and a hint on snapping when a drawing tool is
selected.
Command line
PLAXIS 3D enables to carry out actions using keyboard input by typing the
corresponding commands in the command line. On the other hand, all the actions carried
out using the mouse are translated into commands. All the commands executed in the
project are displayed when the Model history tab is clicked, whereas only the commands
executed in the active session together with the program's feedback are displayed when
the Session tab is clicked.
Information about the available commands is provided when the Command reference
option of the Help menu is selected.
Explorers
Information about the model and the project is given in the explorers. The availability and
the functionality of the explorers depend on the active mode.
The explorers available in the input program are:
Selection explorer It displays the properties of the selected geometry objects and
the features assigned to them and gives the possibility to
hide/show, activate/deactivate these objects (Section 3.7.2).
Model explorer It displays the properties of all the geometry objects and the
features assigned to them in the model and gives the possibility
to hide/show, activate/deactivate these objects (Section 3.7.1).
Phases explorer It shows the list of the calculation phases and enables adding,
Hint: Selection explorer and Model explorer are referred to as Object explorers.
The menu bar of the Input program contains drop-down menus covering most options for
handling files, transferring data, viewing graphs, creating a geometry model, generating
finite element meshes and entering data in general.
The availability of the menus depends on the active mode (Section 3.4). The menus
available in the Input program are:
Select multiple objects To select entities in the model located in the region defined using
the tool Section 3.5.2.
Move object To relocate a user water level in the model Section 5.1.1.
Create water level To create a water level in the model.
Preview phase To preview the phase selected in the Phases explorer.
The Water levels menu is available only in the Water levels mode.
The modelling process is completed in five modes. The mode tabs are shown in the
mode bar. The modes are separated into Geometry and Calculation modes.
The draw area displays the physical model and is automatically updated after any change
in the geometry. The buttons in the toolbar located at the left side of the draw area (side
toolbar) can be used to create and modify the model. The availability of the tools depends
on the active mode. These features are described in the following chapters depending on
their function.
Hint: Right clicking can be used to select entities in the model. After right clicking
an entity or a multiple selections of entities either in the draw area or
explorers, a menu appears displaying the name of the selection and different
options depending on the selection. Clicking the name of the object in the
menu copies it to clipboard. The copied name can be used in the command
line (Section 3.6).
View options
The display of the supporting tools (symbol, grid, global axes, etc.) in the view and their
presence in the draw area is defined in the View tabsheet of the Visualization settings
window.
Symbol size: The size of the symbols is defined in the Symbol size dialog box. The
default values are as shown in Figure 3.7. These values can be used as a reference to
get the preferred display.
Grid: The grid spacing and the number of intervals are available in the Grid dialog box.
To display the grid and ruler select the corresponding option in the Options menu. Points
indicating the grid are shown in the rulers which are displayed at the mouse cursor when
a drawing tool is selected.
Global axes: The availability, the display and the origin of the global axis in the draw
area are defined in the Global axes box.
3D view: The view of the model can be arranged according to the selected option in the
Visualization settings window. The options available are:
Most applicable The perspective view of the model is displayed only when the
corresponding option is selected as default orientation (Section
3.5.4). For the other default view options the objective view of
the model is displayed.
Always perspective The perspective view of the model is displayed regardless of the
options selected as default orientation (Section 3.5.4). The
perspective view mimics real life such as closer objects appear
larger than far away objects, and parallel lines vanish toward the
horizon.
Always objective The model will always be displayed such that two same-size
objects will appear the same size no matter how far apart they
are from each other or from the camera.
Default visibilities
The default visibility of the objects in the model can be defined in the Visibility tabsheet
(Figure 3.8).
The available options are:
The objects are displayed in full opacity. For individual objects in the model, right
click the object then click Show.
The objects are displayed in 70 percent opacity (30 percent transparent). For
individual objects in the model, right click the object, point to Custom and then click
Show 70%.
The objects are displayed in 30 percent opacity (70 percent transparent). For
individual objects in the model, right click the object, point to Custom and then click
Show 30%.
The objects are displayed as wireframes. For individual objects in the model, right
click the object, point to Custom and then click Show as wireframe.
The objects are made completely invisible (hidden). For individual objects in the
The Movement limitation window pops up as one of the default views, except the first
(Perspective view), is selected. The movement limitation simplifies defining three
dimensional geometry. The movement limitation is cancelled by clicking on the
Perspective view button. Definition of geometry is described in detail in Section 3.5.5.
Non-default views
The model in the draw area can be panned or rotated using the left and the middle
mouse buttons. The options are available either as buttons in the general toolbar or as
options in the View menu.
Click the Pan camera button to pan holding the left mouse button pressed. To rotate
the view, keep the <Ctrl> key pressed while dragging the view.
Click the Rotate camera button to rotate holding the left mouse button pressed. To
pan the view, keep the <Ctrl> key pressed while dragging the view.
3.5.5 DRAWING IN 3D
Defining 3D geometry using mouse input is a difficult task considering the fact that the
location in the direction perpendicular to the view cannot be precisely defined. This part
of the manual introduces the tools used in PLAXIS 3D to enable drawing in 3D. The
drawing procedure is illustrated with an example.
The drawing tools available in the side toolbar of the Geometry modes enable defining
the geometry of the model. When one of the drawing tools is selected, the location of the
mouse cursor in the draw area is displayed in the status bar.
• Click at the commandline. Note that the first point in the current line segment is
already defined. Type 0 0 0 to define the last point and press <Enter>.
• Right-click to finalize the process. The resulting geometry is shown in Figure 3.11.
Hint: The Movement limitation window can be closed by selecting the Perspective
view option of the default views.
At the bottom of the input application, a command line panel is available. The command
line provides an alternative way of executing operations, by entering text commands. In
fact, all operations in the input application generate a text command, that is sent to an
interpreter, before being processed.
A feedback panel is situated above the command line panel. Here, the executed
commands are shown, together with the feedback of the command. Successful
commands return green colored feedback, consisting of details on newly generated
geometry or elements. Erroneous commands return red colored feedback, containing an
error report.
A text command consists of a command name, followed by a number of parameters.
Together, they constitute the signature of the command. The number and type of
parameters depends on the particular command, and can be zero or larger.
Commands can be divided into two types:
Global commands Operate on global objects (like the project or the geometry) that
are not necessarily named in the command, and may or may not
have parameters. An example is the undo command.
Target commands Operate on specifically named (target) objects, such as
geometric entities and materials. An example of such command
is the set command.
Command names consist of a simple string. Many commands have a short name in
addition to the standard, explicit long name (for example point and pt). Any command
name can be prefixed with one underscore in order to distinguish it from any objects
bearing the same name (the commands undo and _undo are identical). The application
will always generate the commands with an underscore prefix. It is safe to omit the prefix
in hand written code. There are also commands that start with two underscores. These
are mostly intended for debugging, troubleshooting or advanced automation purposes.
The underscores are in those cases not optional.
Command parameters can be of two different types:
References to objects Represented by strings containing the name of an object. The
names of such objects may not start with an underscore (_).
Value Can be of many data types, such as strings, integers, floating
point numbers and enumerations. String value parameters must
be between single (') or double (") quotes, and may start and end
with either one, two or three of these quotes. Examples of valid
string value parameters are: "hello", "hello 'world'!", "Young's
"modulus"". Floating point numbers must always use the dot
decimal notation, even on systems that normally use comma as
decimal separator. Enumerations may be represented either by a
string value, or by the corresponding integer index of that value.
In most cases where a value is expected, a property of that type
is also allowed - this is indicated by a trailing single quote. For
example when trying to set the x coordinate of a point, this can
be done as set point_1.x 5.2 or as set point_1.x point_2.x.
In some cases it is allowed to place parameters between parentheses, which is indicated
by a trailing single quote in the parameter signature of a command. For example a point
can be made with the command point 1 2 3 or as point (1 2 3). The parentheses are
optional and are mostly useful for legibility, but may in some cases be necessary in order
to distinguish between different types of calls to a given method.
In PLAXIS 3D, a global command is included that shows information on the commands
that can be executed: cms. This command lists the signatures of all commands available
in PLAXIS 3D.
As an example, the signature for the delete command is shown here:
delete (del)
Material'
In the first line, the command name and its abbreviation are shown. Each of the following
lines shows a different set of parameters that can be used with this command. In this
example there are five different signatures for the delete command. The delete command
can be either used to delete a single material (referenced by its name), a set of features,
a set of geometric elements, a set of boreholes, or a single soil layer. Note that the
signatures do not allow a feature, a geometric element and a borehole to be deleted in
cylinder (cyl)
Number' Number'
Number' Number' Integer' <Coords: Number' Number' Number'>' <Coords: Number' Number' Number'>'
If we take a closer look at the bottommost line, the signature consists of two numbers, an
integer value and two sets of coordinates. The first two parameters represent the radius
and height parameters, respectively. Then an integer value representing the accuracy of
the cylinder is expected, followed by the position of the cylinder, and a vector
representing the orientation of the cylinder.
An example of the use of the cylinder command is:
• cylinder 3 12 8 (0 0 0) (0 0 1)
• cylinder 3 12 8 (point_1.x point_1.y point_1.z) (0 0 1)
An extended overview and description of all available commands is accessible via the
Help menu in the user interface of the PLAXIS 3D Input program.
>line (1 2 0)(5 1 0)(5 3 0)(4 7 0) # Create points with coordinates specified in the brackets and lines connecting them
>tabulate Points "x y z" # List the created points and their coordinates
Object x y z
Point_1 1 2 0
Point_2 5 1 0
Point_3 5 3 0
Point_4 4 7 0
Note that indexing is done according to the order of creation of the geometric entities.
The indexing of features such as soil, structural elements, loads or prescribed
displacements, is NOT done according to the order of their assignment to geometry but
according to the order of creation of the geometric entities they are assigned to.
3.7 EXPLORERS
Explorers display information about the physical model, the calculation phases and their
configuration.
Model explorer: Information related to all the objects in the physical model is
given.
Selection explorer: Information related to the object (or group of objects) selected in
the draw area is given. From a group selection, only the
information valid for all the objects in the group is displayed.
Phases explorer: A list of the calculation phases defined in the project is displayed.
The Phases explorer is available only in Calculation modes.
However, as phases are defined only in the Phase definition
modes, it is greyed out in the Mesh mode.
Hint: Model explorer and Selection explorer are referred to as Object explorers.
Attributes library contains the set of water levels that have been created from the
information in the boreholes, the soil clusters and the user-defined water levels.
Dynamic multipliers All the displacement and load dynamic multipliers defined in the
project are listed under Dynamic multipliers.
Flow functions All the flow functions defined in the project are listed under Flow
functions.
Materials All the material sets assigned to the entities in the model are
listed under Materials. The identification and the colour
representing the material data set are displayed.
Water levels All the water levels created in the model are listed.
Borehole water levels The water levels generated according to
the water conditions defined in the
boreholes are listed.
User water levels The water levels created in the Water
levels mode are listed.
Note that the Water levels subtree is available only in the
Calculation (green) modes.
Since these attributes are Global, any change in such an attribute will affect the model as
a whole. If it is the intention to change an attribute only in a particular calculation phase, it
is advised to duplicate the attribute using the corresponding option in the right-hand
mouse menu, and make the change in the copy.
Geometry: All the geometric entities created in the model are listed under Geometry.
Note that when a geometric entity is created, its sub-entities are automatically created by
the program. When a line is created, the first (starting) point and the second (last) point
are automatically created. New geometric entities are also automatically created when a
geometric entity is decomposed.
Boreholes: Under the Borehole heading, all the boreholes created in the model are
listed. For each borehole, the x and y coordinates of each borehole and the Head level
are given.
Information about boreholes is available in the Model explorer only in the Geometry
modes. The values of coordinates and Head can be modified in the corresponding cells.
Connections: All the connections created explicitly by the user in the Staged
construction mode are listed. This option is available in the Model explorer only in the
Phase definition modes.
Beams: All the beams created in the model are listed. Information about the assigned
material is given for each beam.
Embedded piles: All the embedded piles created in the model are listed. Information
about the assigned material and the connection type is given for each embedded pile.
Soils: All the soil created in the model are listed. Information about the assigned
material, volumetric strain and water conditions is given for each soil.
Plates: All the plates created in the model are listed. Information about the assigned
material is given for each plate.
Geogrids: All the geogrids created in the model are listed. Information about the
assigned material is given for each geogrid.
Surface loads: All the surface loads created in the model are listed. The distribution
and the values of the load components can be defined.
Contractions: All the contractions created in the model are listed. The distribution and
the values of the contractions can be defined.
Surface prescribed displacement: All the surface prescribed displacements created in
the model are listed. The distribution and the values of the prescribed displacement
components can be defined.
Line loads: All the line loads created in the model are listed. The distribution and the
values of the load components can be defined.
Line prescribed displacement: All the line prescribed displacements created in the
model are listed. The distribution and the values of the prescribed displacement
components can be defined.
Point loads: All the point loads created in the model are listed. The values of the load
components can be defined.
Point prescribed displacement: All the point prescribed displacements created in the
model are listed. The values of the prescribed displacement components can be defined.
Interfaces: All the positive and negative interfaces created in the project are listed.
Either the material of the adjacent soil or a new material can be assigned. Permeability
condition can be specified.
Node-to-node anchors: All the node-to-node anchors created in the model are listed.
Information about the assigned material is given for each node-to-node anchor.
Fixed-end anchors: All the fixed-end anchors created in the model are listed.
Information about the assigned material is given for each fixed-end anchor. The
components and the equivalent length can be defined for each fixed-end anchor.
Surface drains: All the surface drains created in the model are listed. The head can be
defined for each surface drain.
Surface groundwater flow BC's: All the surface groundwater flow boundary conditions
created in the model are listed. Information about the assigned behaviour (Seepage,
Closed or Head), is given for each boundary condition. The assigned behaviour can be
changed by selecting the option in the Behaviour drop-down menu.
Wells: All the wells created in the model are listed. The behaviour (Extraction,
Infiltration), the discharge (|qwell |) and the minimum head (hmin ) can be defined for each
well.
Line drains: All the line drains created in the model are listed. The head can be defined
for each line drain.
Model conditions: The general boundary conditions of the model as a whole for each
calculation phase can be specified in the Model conditions subtree in Model explorer.
The conditions at the extremities of the model are considered. This makes it possible to
make a quick selection of general boundary conditions that apply to the model as a
whole. The Model conditions subtree is available in the Calculation modes in the Input
program. Modification of the information given in the Model conditions subtree is possible
only in the Water levels and Staged construction modes. Note that any change in the
model boundaries is applied only to the phase selected in Phases explorer.
The options available are:
Deformations A set of general fixities is automatically applied to the boundaries
of the geometry model for the selected calculation phase. More
information on default fixities is given in Section 7.9.9. The
default general fixities are removed when the corresponding
option is set to False. In that case it is required to set the
appropriate boundary conditions manually.
Dynamics The model conditions for a dynamic analysis can be defined at
the extreme boundaries of the model. The options available are
None and Viscous. Besides the boundary conditions the
relaxation coefficient (C1 and C2 ) can be defined for each phase.
For more information on boundary conditions for Dynamic
calculations see Section 7.9.9.
GroundwaterFlow The model conditions for a groundwater flow calculation, a
consolidation analysis or a fully coupled flow-deformation
analysis can be defined at the extreme boundaries of the model.
The options available are Open and Closed.
Precipitation This option can be used to specify a general vertical recharge or
infiltration (q) due to weather conditions to all boundaries that
represent the ground surface. Other types of flow boundary
conditions can be defined using the corresponding feature in the
Structures mode (Section 5.7).
Water The global water level for a selected calculation phase can be
assigned in the Water subtree under Model conditions in Model
explorer. More information of the usage of the global water level
is given in Section 7.5.2.
The soil stratigraphy is defined in the soil mode using the Borehole feature of the
program. When a new project is created, the soil contour defined in the Project properties
window is displayed in the draw area.
The model contour can be defined using the options available in the side toolbar.
To modify the model contour click the Adjust soil contour button. The contour is
defined in the draw area. The buttons in the expanded menu can be used to modify the
model contour.
Click the Move contour points/lines button and drag any point or line in the contour
to a new location.
Click the Insert contour points button, click the contour where a new point is to be
defined and drag it to the desired position of the new point. The soil contour will be
modified accordingly.
Click the Delete contour points button and click the contour point that is to be
deleted.
Besides the expanded menu, the Surface points window pops up displaying the
coordinates of the points defining the model contour. Points defining the contour can be
relocated by entering the coordinates of the new location in the table. The options to
move, insert and delete contour points are available in the menu appearing as the
contour point number is right clicked in the table. A more detailed description on Surface
points window is given in the Section 5.2.4.
Boreholes are locations in the draw area at which the information on the position of soil
layers and the water table is given. If multiple boreholes are defined, PLAXIS 3D will
automatically interpolate between boreholes, and derive the position of the soil layers
from the borehole information. Each defined soil layer is used throughout the whole
model contour. In other words, all soil layers appear in all boreholes. The top and the
bottom boundaries of the layers may vary through boreholes, making it possible to define
non-horizontal soil layers of non-uniform thickness as well as layers that locally have a
zero thickness (Figure 4.5).
Hint: A warning message appears in the Modify soil layers window when multiple
boreholes are located at the same position.
• To select the source borehole click the corresponding check box. The thickness of
the layers in the new borehole is equal to that in the source borehole.
• Select the appropriate option in the Shape drop-down menu.
• Define the number of columns (and/or rows) and the distance between them. A
detailed description of the Array feature is given in Section 5.1.4.
The first borehole, created at the start of a new project, contains no soil layers. Soil layers
in boreholes are defined in the Modify soil layers window.
The Modify soil layers window contains the following items (Figure 4.2):
Soil column: Graphs of all boreholes with indication of their locations, water
head, layer boundary levels and soil material.
Upper buttons: Buttons to add, insert or delete soil layers.
Soil layers tabsheet: Table showing the boundaries and the materials of the soil
layers.
Water tabsheet: Table showing the water conditions in every layer and the top
and bottom pressure values per borehole in each layer boundary.
Initial conditions tabsheet:
Table showing the name, material model and the initial stress
condition parameters, OCR , POP , K 0x and K 0y for all soil
layers.
Surfaces tabsheet: This tabsheet enables the import of predefined top and bottom
soil surfaces.
Field data tabsheet: This tabsheet enables the import of CPT measurements. More
information on the usage of CPT logs to generate the soil
stratigraphy is given in Section 4.3.5.
Bottom buttons: Buttons to add, select, visualize or sort boreholes, to open the
material database and to accept (OK ) the soil layers. Pressing
Ok closes the Modify soil layers window.
If a certain layer does not exist in the current borehole but it exists in other boreholes, it
should be eliminated only in the soil column of the current borehole by setting its bottom
level equal to its top level.
The program generates water levels according to the information provided in the Water
tabsheet. Water levels can represent external water levels as well as phreatic levels in the
soil. The water level created according to the Head levels in the boreholes is
automatically assigned as Global water level in the calculation phases. A more detailed
description is given in Section 7.8.
By default, positive pore water stress (suction) in the unsaturated zone above the phreatic
level is ignored. However, PLAXIS can deal with suction in calculations. This requires the
selection of an appropriate soil water retention curve in the material data set (Section
6.1.3). It also requires that suction is allowed during the calculations (Section 7.7.4).
Figure 4.7 Initial conditions tabsheet of the Modify soil layers window
4.3.4 SURFACES
By default, the ground level surface and the bottom layer of the soil are modelled using
the information provided by boreholes, as indicated by the From boreholes option in the
Surfaces tabsheet (Figure 4.8). The surface level between two boreholes is linearly
interpolated. Alternatively, predefined surfaces can be used to define the top and bottom
surfaces of the soil volume. In order to use predefined surfaces, the Custom option in the
Surfaces tabsheet is selected. An Import (top/bottom) soil surface window appears from
which the desired surface can be selected.
Click the expansion button next to the borehole the stratigraphy of which is to be
created according to the imported data.
• Select the data to be considered from the Field data drop-down menu.
• Select the Interpretation method (currently only the CUR 3 layers option is available.
• Enter the minimal layer thickness (Minimal thickness option) to be considered in
layer generation to avoid multiple very thin layers. The expanded Modify soil layers
window is shown in Figure 4.9.
• Click Apply layers to subdivide the soil cluster into the layers shown in the borehole.
Hint: When soil stratigraphy is generated from CPT logs the created layers will be
applied through the whole model. Boreholes can be used to modify the layer
thicknesses at their location.
» Only one CPT log can be used in a project. Generation of layers resulting
from a new CPT log will overwrite the existing ones.
The geometry of the soil can be imported from predefined files. Instead of using the
Borehole tool, the Import soil tool can be used. Note that the soil material will not be
imported. A detailed description on how geometry is imported is given in Section 5.9.
Geometric entities, structural elements and boundary conditions are defined in the
Structures mode. Click the Structures tab to switch to this mode. The soil layers created
in the Soil mode will appear in the draw area of the Structures mode according to the
default visibilities defined in the Visualization settings option on the Options mode, i.e.
somewhat transparent. Boreholes cannot be modified in the Structures mode.
The assisting tools enable modifying the geometry by changing the location or the
orientation of the object in the model or by creating higher objects with a higher
dimension (from lines to surfaces or from surfaces to volumes).
• Define the coordinates of the rotation point. The default rotation point is the object
reference point. The rotation point can be relocated by dragging it in the draw area.
For cases where the rotation point is not visible, it is possible to drag the rotation
point by clicking the Drag point button in the Rotate object window and dragging it in
the draw area.
• Define the rotation angles. The rotation angles are positive for counterclockwise
rotation. Alternatively, the rotation may be defined by dragging one of the points on
one of the circles.
In the draw area, the circle defined by rotation over the x -axis is marked in red colour.
The circle defined by rotation over the y -axis is marked in green colour. The circle
defined by rotation over the z -axis is marked in blue colour.
Figure 5.2 Rotation point and rotation axis for a selected line
Hint: Note that the options to rotate or pan camera do not change the orientation of
the objects in the model but only the view.
Extrusion
The extrusion of points, lines and surfaces can be defined in the Extrusion tabsheet in
the Extrude object window. As an example, the extrusion of a line is described.
To extrude a line:
• Select a line in the draw area.
Click the Extrude object button. The Extrude line window pops up (Figure 5.3).
• Define the components of the extrusion vector. The length of the extrusion vector is
automatically calculated.
• If the length of the extrusion vector is redefined, the vector components are
automatically recalculated according to their previous values and the new length.
This option enables defining the extrusion by entering the skew components in the first
and second local axis and the height of the extrusion in the direction of the third local axis
(Figure 5.5).
skew (axis 2)
height
skew (axis 1)
Figure 5.5 Skew components
• Select the necessary option in the Shape drop-down menu. Enter the total number
of rows (or columns) after copying and define accordingly the distances between
adjoining copies.
Hint: Note that when Array is used to create copies, besides the geometry entity
itself, the assigned features are copied as well.
The geometric entities are the basic components of the physical model. Features such as
structures, loads etc., can be assigned to geometric entities. The created entities in the
model are automatically named by the program. They can be renamed in the Object
explorers.
A geometric entity can be relocated by selecting it first using the Selection tool in the
toolbar and dragging it to the desired position. Alternatively, the location coordinates can
be changed in the Object explorers.
There are two options of deleting entities in the model. To delete a geometric entity and
all the features assigned to it, right click it and select the Delete option from the appearing
menu. To delete only features assigned, right click it, click on the feature you want to
delete in the appearing menu and click Delete.
Hint: Note that using the <Delete> key, deletes the geometric entity and all the
features assigned to it.
5.2.1 POINTS
One of the basic input items for the creation of a geometry is a point. To create a
point click the Create point button in the side toolbar and then click at the location of
the point. A Point load, Point prescribed displacement and Fixed-end anchor feature can
be assigned to a point.
5.2.2 LINES
Another basic input item for the creation of a geometry is a line. To create a line
click the Create line button in the side toolbar and then click at the location of the
first and the second point of the line in the draw area. Optionally, additional points may be
created resulting in a chain of lines. After creating the last (end) point of the line,
right-click or use the <Esc> key to stop drawing. When a line is created, its first and
second points are also created.
properties of structural elements having orthotropic material behaviour and to show the
forces in the Output program.
The local coordinate system can be modified in the Model explorer. The first local axis
always corresponds to the axial line direction. The Automatic axis function defines
automatically the second axis of the line according to the global coordinate system.
When the Manual option in the AxisFunction drop-down menu in the Model explorer
(Figure 5.7) is selected, the user has to define the direction of the second axis. In both
cases the third axis is determined automatically from the first and the second axis
according to the right hand rotation rule.
A Beam, Embedded pile, Line load, Line prescribed displacement, Node-to-node anchor,
Well and Line drain can be assigned to a line.
5.2.3 POLYCURVES
A polycurve is a collection of lines and arcs that form a single continuous curve. Each
line and arc composing the polycurve is defined as a segment. The starting point of each
new segment is the end point of the previous segment.
To create a polycurve:
Click the Create polycurve button in the side toolbar.
• Click the location where the polycurve is to be inserted. The Shape designer
window is displayed (Figure 5.8).
The Shape designer window consists of a toolbar and two tabsheets where the
properties of the polycurve are defined and a display area showing the created polycurve.
Polycurve layout in the model: The general properties of the polycurve in the model
are defined using the options available in the Shape tabsheet.
• The coordinates of the insertion point can be defined in the Object explorers. By
default the coordinates of the location clicked in the model before the Shape
designer window is displayed are shown. As the coordinates are changed the
polycurve will be relocated.
• The Orientation axis 1 and the Orientation axis 2 define a plane in which the
polycurve is drawn.
• By default the insertion point is the starting point of the first segment of the
polycurve. The Offset to begin point group can be used to specify offset distances
from the orientation axis when insertions point (origin of the orientation axis) and the
starting point of the first segment of the polycurve do not coincide.
the newly added segment. The newly created segment is displayed in the segment list
(Figure 5.10).
Figure 5.10 Segment list displaying the segments composing the polycurve
Radius
Segment angle
Discretization
angle
When specifying a negative radius for an Arc segment, the arc will be mirrored in the line
described by the tangent vector. Note that when the Circular shape option is selected,
the created polycurve consists only of one Arc segment with a Segment angle value of
360◦ . An arc is discretized into line segments (cords). The angle subtended by each cord
along the discretized arc is called the Discretization angle.
Relative start angle To define the value of the angle between the start tangent of the
next segment and the end tangent of the previous segment
(Figure 5.12).
segment 2
segment 1
To extend the polycurve to the symmetry axis click the corresponding button in the
toolbar of the Shape designer window. A new segment is added connecting the last
point of the polycurve to the axis of symmetry.
To close the polycurve by mirroring the existing segments, click the Symmetrically
close polycurve button in the toolbar of the Shape designer window. This is only
possible when one half is completed up to the symmetry axis.
To close the polycurve click the corresponding button in the toolbar of the
Shape designer window. A new segment is created connecting the last point of the
polycurve to the first one.
Modification of polycurve
Modification of a polycurve is possible in the Model explorer after closing the Shape
designer window (Figure 5.13).
The menu appearing when a polycurve is right-clicked either in the model or in the Object
explorers displays the following specific options for polycurves:
• Add segment
• Extend to symmetry axis
• Close symmetrically
• Close Polycurve
• Create surface
5.2.4 SURFACES
The third basic input item for the creation of a geometry is a surface. To create
a surface click Create surface and then on the locations of the points that define the
surface.
The first point defined in the surface is the reference point of the surface. After the first
point is created, the user may draw a line by entering another point. After the first line is
created, the user may draw a surface by entering another point. As the third surface point
is defined, the Surface points window pops up (Figure 5.14).
Hint: The plane of the created polygon is displayed in the model as the third point
is defined. The program allows creation of new points located only in this
plane by clicking.
» To add new points out of the plane, right-click the points in the Surface points
window and select either the Insert point or the Add point option.
The points defining the surface and their coordinates are listed in the Surface points
window. Crossing lines are indicated in the Plane state option. Existence of crossing
lines indicates that no valid surface can be created. When there are no crossing lines, the
condition of the Plane state is OK. The options to define the plane of the surface are:
First points: The first three points in the Surface points window are
considered to define the surface plane. Values in the error
column indicate whether the surface is in or out of the defined
plane.
Best fit: Best fit of surface points is considered to define a plane surface
that fits the best to the defined points. The error column is not
applicable.
The points defining the surface can be projected to the plane of the surface by clicking
Project points to plane in the Surface points window. The creation of a surface is finished
by clicking on the selection tool in the side tool bar.
Hint: Note that the points and the lines defining a surface are not shown in the
Model explorer.
» The points and lines defining the surface can be added to the model by
selecting the Decompose into outlines option from the right mouse button
menu.
Modification of surfaces
The points defining surfaces can be modified using the options available in the menu that
appears when the point ID is right-clicked in the Surface points window. Using these
options a point can be added, inserted or deleted.
The points defining the surface can be relocated by changing their coordinates in the
corresponding cells in the Surface points window.
Note that the options to modify the points defining a surface are available in the Create
surface submenu (Figure 5.16) and they are selectable after the surface to be modified is
double clicked first.
menu in the Model explorer assigns the first axis according to the global coordinate
system. The direction of the first axis can be defined by the user if the Manual option is
selected. For example, in case of a straight tunnel it will be convenient to use the
direction of the tunnel as the direction of the first axis. The direction of the second axis is
determined automatically from the first and third axis according to the right hand rule.
A Plate, Geogrid, Interfaces, Surface load, Surface prescribed displacement, Contraction,
Surface groundwater flow boundary condition and Surface drain can be assigned to a
surface (either created or imported).
Imported surfaces
Predefined surfaces can be imported in the Structures mode (Section 5.9). The geometry
of an imported surface can be modified after being decomposed into outlines.
5.2.5 VOLUMES
In PLAXIS 3D volumes can be created by extruding surfaces (Section 5.1.3). Import
volume feature enables usage of predefined volumes (Section 5.9). Note that the soil
volumes generated according to the soil layers defined in the Modify soil layers window
and the volumes imported in the Soil mode are listed under Soil volumes in the Geometry
subtree of the Model explorer. The volumes created by extruding surfaces or that are
imported in the Structures mode are listed under Volumes in the Geometry subtree.
Hint: The surfaces and the outlines defining a volume created or imported in the
Structures mode can be added to the model by selecting the Decompose into
surfaces or Decompose into outlines option respectively from the right
mouse button menu.
Soil and interfaces material data sets, Volumetric strain and Water conditions can be
assigned to volumes.
Hint: Multiple entities can be selected by keeping the <Ctrl> key pressed while
clicking on the entities.
Hint: Note that Volume strain can be assigned to soil clusters only in the Phase
definition modes.
» If an embedded pile or a node-to node anchor is assigned to a line in the
model, no other features can be assigned to it.
Invert normal
The Invert normal option is available for surfaces in the right mouse button menu. Note
that the location of an interface relative to the surface to which it is assigned depends on
the type of the interface (positive or negative) and the normal of the surface. If the normal
of a surface is inverted, the assigned interface will automatically be applied to the other
side of the surface.
mergeequivalents geometry
The Merge equivalent geometric objects command may also be used to merge objects
that are not exactly located at the same position, but very close to each other. By default,
the distance tolerance is 0.001 length units, equal to the standard snap distance (see
below), but the user can include the distance tolerance in the command, for example:
Snap
Drawing imperfections such as disconnected or overlapping geometric entities can lead
to intersection or meshing difficulties. The Snap feature enables the correction of such
imperfections in a range of 0.001 length units. The default value of snapping distance
can be temporarily modified when snapping is performed using the corresponding
Decompose
The Decompose option enables decomposing the selected geometric entity into lower
entities. A surface is decomposed into outlines and a volume is decomposed into
surfaces. An example of decomposition of a surface into outlines is given in the Tutorial 3
of the Tutorial Manual.
Combine
A multiple selection of geometric entities can be combined (merged/fused together) by
selecting the corresponding option in the right mouse menu. The Figure 5.18 displays
two volumes before and after combining. After combining different entities into one new
entity, the original entities do not exist anymore.
Group
The Group option in the right mouse menu enables creating groups of entities in the
model. The created group will be listed under the Groups subtree in the Model explorer.
Grouping enables modification of the common properties of the group at the same time.
After grouping different entities into a group, the original entities still exist.
5.3 LOADS
Loads are features that can be assigned to geometric entities. Loads can be assigned to
existing geometric entities by right clicking the entity and selecting the corresponding
option in the appearing menu. Alternatively, the options available in the Create load
menu can be used. To assign loads to existing geometry entities:
Select the button corresponding to the load from the menu appearing when the
Create load is selected (Figure 5.19).
• Drag the tool (keep the left mouse button down) and drop it on the geometric entity.
A geometric entity can be created and a load can be assigned to it simultaneously using
the options available in the menu displayed as the Create load button is selected in the
side toolbar, (Figure 5.19). This provides a faster definition of loads. Instead of creating
the geometric entity and then assigning a load to it, the process is completed in one step.
The process is similar to the creation of a geometric entity with the difference that the
load is assigned to it as well.
Although the input values of loads are specified in the Structures mode, the activation,
deactivation or change of loads may be considered in the framework of Staged
construction (Section 7.9.3).
On a part of the geometry where both fixities and loads are applied and active, the fixities
have priority over the loads during calculations. Hence, it is not useful to apply loads on a
fixed geometric entity. However, it is possible to apply a load in the free directions.
The Load subtree in the Object explorers consists of two parts where the static and the
dynamic components of the load can be defined respectively.
Static load
The distribution and the components of the static load are assigned in the first part of the
Load subtree. Note that the Distribution option is not available for Point loads.
Dynamic load
The distribution and the components of the dynamic load are assigned in the second part
of the Load subtree when the subtree of the dynamic load is extended. Note that the
Distribution option is not available for Point loads. Besides the Distribution and the
components of the load, Multipliers can be defined for each of the load components
separately. One of the available multipliers defined in Attributes library (Section 5.5.1)
can be selected from the drop-down menu. The view of the dynamic load in the Object
explorers is given in Figure 5.20.
Prescribed displacements are special conditions that can be imposed to the model to
control the displacements at certain locations. Prescribed displacements can be
assigned to existing geometric entities by right clicking the entity and selecting the
corresponding option in the appearing menu. Alternatively, the options available in the
Create prescribed displacement menu can be used. To assign prescribed displacements
to existing geometry entities:
Select the button corresponding to the prescribed displacement from the menu
appearing when the Create prescribed displacement button is selected (Figure
5.30).
• Drag the tool (keep the left mouse button down) and drop it on the geometric entity.
The creation of a prescribed displacement is similar to the creation of a geometric entity.
The options available in the menu displayed as the Create prescribed displacement
button is selected in the side toolbar, (Figure 5.24), provide a faster definition of
prescribed displacements. Instead of creating the geometric entity and then assigning a
prescribed displacement to it, the process is completed in one step.
Although the input values of prescribed displacements are specified in the Structures
mode, the activation, deactivation or change of values may be considered in the
framework of Staged construction (Section 7.9.1).
On a geometric entity where both prescribed displacements and loads are applied and
active, the prescribed displacement has priority over the load during calculations. If both
prescribed and fixed displacement are assigned to a line, the fixed displacement will be
taken into consideration. However, when not all the displacement directions are fixed, it is
possible to apply a load in the free directions.
The Prescribed displacement subtree in the Object explorers consists of two parts where
the static and the dynamic components of the prescribed displacement can be defined
respectively.
5.4.4 CONTRACTION
Contraction can be used to 'shrink' a surface or to simulate a volume loss around a tunnel
lining. Contraction is specified as a strain without any unit associated.
In order to define a contraction on surfaces, the Create surface contraction
option should be selected from the menu appearing as the Create prescribed
displacement icon is clicked.
The options to define a surface contraction are:
Constant: A constant contraction is defined. The contraction is applied in
the local Axis2-direction.
Axial increment: A linear contraction varying in the local Axis1-direction is
defined. The contraction is applied in the local Axis2-direction.
Hint: When contraction is used to simulate the volume loss around a TBM tunnel,
note that the value of surface contraction to be applied should be half the
volume loss as a percentage of the tunnel volume.
the Edit option is selected in the appearing menu (Figure 5.26). The window consists of
two tabsheets, i.e Displacement multipliers and Load multipliers where the multipliers
applicable to the prescribed displacements and to the loads can be defined respectively.
The buttons available under the multipliers tabs can be used to add new multipliers or to
delete the selected multiplier in the list. When a new multiplier is added, the options to
define it are displayed.
Name The name of the multiplier can be defined.
Signal The type of the signal can be assigned. The options available
are Harmonic and Table. More information on definition of a
harmonic multiplier and definition of the multipliers in a table is
given in Section 5.5.2 and Section 5.8.2 respectively.
Type The options Displacements, Velocities and Accelerations are
available. Note that the Type drop-down menu is available only
for displacement multipliers. For Load multipliers the type does
not need to be specified.
Drift correction It is used to correct the displacement drift. Due to the integration
of the accelerations and velocities, a drift might occur in the
displacements. When selected, the displacement drift is
corrected by applying a low frequency motion from the beginning
of the calculation and by correcting the acceleration accordingly.
F = M̂ F̂ sin (ω t + φ0 )
in which:
M̂ Amplitude multiplier
F̂ Input value of the load
ω 2 π f with f = Frequency in Hz
φ0 Initial phase angle in degrees
(M̂ F̂ is the amplitude of the dynamic load.)
Hint: A dynamic load can also suddenly be applied in a single time step or sub
step (block load). In case of a Harmonic load multiplier, a block load can be
modelled by setting the Amplitude multiplier equal to the magnitude of the
block load, the Frequency to 0Hz and Initial phase angle to 90◦ giving the
relation F = M̂ F̂ . In case of a Load multiplier from data file, a block load
can directly be defined.
Figure 5.28 shows the Multipliers window where a Harmonic signal is defined and plotted.
Figure 5.28 Definition and display of a Harmonic signal in the Dynamic multipliers window
Click the Insert button to insert a new row before the selected row in the table.
Click the Delete button to delete the selected row in the table.
The values can be defined by clicking the cell in the table and by typing the value.
Hint: PLAXIS assumes the data file is located in the current project directory when
no directory is specified in the Dynamic loading window.
A signal, either defined in a table or received from a file and modified, can be
saved using the Save button in the toolbar enabling the usage of the signal in other
projects or validating the effect of the modifications in the current project.
A signal, either defined in a table or received
from a file and modified, can be copied using the Copy button in the toolbar.
Copied data from other applications (using <Ctrl+c>) can be imported by using the
Paste button. The Import data window appears (Figure 5.29). The starting row of
the data to be imported can be defined in the From row cell. The data can be parsed as
Plain text files or Strong motion CD-Rom files (SMC). The data and the plot is displayed
in the Dynamic multipliers window after pressing OK.
ASCII file: An ASCII file can be created by the user with any text editor. In every line a
pair of values (Dynamic time and corresponding multiplier) is defined, leaving at least one
space between them. The time should increase in each new line. It is not necessary to
use constant time intervals.
If the time steps in the dynamic analysis are such that they do not correspond with the
time series given in the file, the multipliers at a given (Dynamic) time will be linearly
interpolated from the data in the file. If the Dynamic time in the calculation is beyond the
last time value in the file a constant value, equal to the last multiplier in the file, will be
used in the calculations.
SMC file: In addition, it is possible to use earthquake records in SMC-format as input for
earthquake loading. The SMC (Strong Motion CD-ROM) format is currently used by the
U.S. Geological Survey National Strong-motion Program to record data of earthquakes
and other strong vibrations. This format uses ASCII character codes and provides text
headers, integer headers, real headers, and comments followed by either digitised
time-series coordinates or response values. The header information is designed to
provide the user with information about the earthquake and the recording instrument.
Most of the SMC files contain accelerations, but they may also contain velocity or
displacement series and response spectra. It is strongly recommended to use corrected
earthquake data records, i.e. time series, that are corrected for final drift and non-zero
final velocities. SMC files should be used in combination with prescribed boundary
displacements at the bottom of a geometry model.
The strong motion data are collected by the U.S. Geological Survey and are available
from the National Geophysical Data Center (NGDC) of the National Oceanic and
Atmospheric Administration. Information on NGDC products is available on the
World-wide Web at http://www.ngdc.noaa.gov/hazard or by writing to:
National Geophysical Data Center NOAA/EGC/1
325 Broadway
Boulder, Colorado 80303
USA
Hint: The time value in dynamic multipliers always refers to the global dynamic
time in the whole list of calculation phases, rather than to the time interval of
an individual phase. This means that in series of continuing dynamic
calculation phases each phase will only use its consecutive part of the
dynamic multiplier.
Note that the assignment of structures is also possible by right-clicking the geometric
entity either in the draw area or in the Object explorers and selecting the corresponding
option in the appearing menu. Assignment of materials to structures is described in
Section 6.9.
Structures can be activated, de-activated or the assigned material data sets can be
modified in calculation phases.
5.6.2 BEAMS
Beams are structural objects to model slender (one-dimensional) structures with a
significant flexural rigidity (bending stiffness) and an axial stiffness. The creation of
a beam is similar to the creation of a geometry line (Section 5.2.2).
Beam elements
After meshing beams are composed of 3-node line elements (beam elements) with six
degrees of freedom per node: three translational degrees of freedom (ux , uy and uz ) and
three rotational degrees of freedom (φx , φy and φz ). The beam elements are based on
Mindlin's beam theory (Bathe (1982)). This theory allows for beam deflections due to
shearing as well as bending. In addition, the element can change length when an axial
force is applied. Note that the element cannot sustain torsion. Element stiffness matrices
are based on the properties as defined in the material data sets and numerically
integrated from the four Gaussian integration points (stress points) along the length of the
element and numerically integrated from the 4 (2x2) Gaussian integration points√over the
cross section of the element. These 4 stress points are √ located at a distance 1/6 3deq3 in
the local 3-direction of the beam and at a distance of /6 3deq2 in the local 2-direction of
1
or the right hand mouse button menu in the draw area (Section 6.9).
Hint: Since installation effects cannot be considered, the embedded piles option
should be primarily used for pile types that cause a limited disturbance of the
surrounding soil during installation. This may include some types of bored
piles, but obviously not driven piles or soil displacement piles.
The creation of a pile is similar to the creation of a geometry line (Section 5.2.2). The
point of connection of the pile to the surrounding geometry and the type of the connection
can be specified in the Object explorers.
In the case of the pile connection point coinciding with soil elements, the connection
behaviour is as follows:
Free The Connection point is free to move and rotate relative to the
surrounding soil. The pile and the surrounding soil interact only
through the special interface elements.
Hinged The displacement at the connection point of the pile is directly
coupled with the displacement of the soil element in which the
pile connection point is located, which means that they undergo
exactly the same displacement. However, their rotation may be
different.
Rigid The displacement and rotation at the pile connection point are
both coupled with the displacement and rotation of the soil
element in which the pile top is located, provided that this
element has rotational degrees of freedom. This option only
applies if the pile connection point coincides with structural
elements like plates or beams.
In the case of the pile connection point coinciding with structural elements (without
interface elements) the connection behaviour is as follows:
Free The Connection point is free to move and rotate relative to the
surrounding soil and the structural element. The pile interacts
with the surrounding soil and the structural elements only
through the special interface elements. Note that the structural
elements are not free to move and rotate relative to the
surrounding soil.
Hinged The displacement at the connection point of the pile is directly
coupled with the displacement of the structural element and/or
the surrounding soil in which the pile connection point is located,
which means that they undergo exactly the same displacement.
Rigid The displacement and rotation at the pile connection point are
both coupled with the displacement and rotation of the structural
and/or soil element in which the pile top is located.
In the case of the pile connection point coinciding with structural elements which
interaction with the surrounding soil is moddelled using interfaces, the connection
behaviour is as follows:
Free The Connection point is free to move and rotate relative to the
surrounding soil and the structural element. The pile interacts
with the surrounding soil and the structural elements only
through the special interface elements. In the same way, the
structural element interacts with the surrounding soil through the
defined interface elements. This is the default option in the
program.
Hinged The displacement at the connection point of the pile is directly
coupled with the displacement of the structural element in which
the pile connection point is located, which means that they
undergo exactly the same displacement. Note that the structural
element interacts with the surrounding soil through the defined
interface elements.
Rigid The displacement and rotation at the pile connection point are
both coupled with the displacement and rotation of the structural
element in which the pile top is located. Note that the structural
element interacts with the surrounding soil through the defined
interface elements.
Note that when an anchor is present in the model at the position of an embedded pile
connection point, the anchor will be automatically connected to the embedded pile
connection point and not to the corresponding point in the soil element. The options
specified for the embedded pile in the Object explorers are valid for the connection of the
anchor-pile connection to the surrounding soil.
Hint: When embedded piles penetrate a volume cluster with linear elastic material
behaviour, the specified value of the shaft resistance is ignored. The reason
for this is that the linear elastic material is not supposed to be soil, but part of
a structure. The connection between the pile and the structure is supposed
to be rigid to avoid, for example, punching of piles through a concrete deck.
After meshing, the beam elements are 3-node line elements with six degrees of freedom
per node: Three translational degrees of freedom (ux , uy and uz ) and three rotational
degrees of freedom (ϕx , ϕy and ϕz ). Element stiffness matrices are numerically
integrated from the four Gaussian integration points (stress points). The element allows
for beam deflections due to shearing as well as bending. In addition, the element can
change length when an axial force is applied.
The special interface elements are different from the regular interface elements as used
along walls or volume piles. Therefore, at the position of the beam element nodes, virtual
nodes are created in the soil volume element from the element shape functions. The
special interface forms a connection between the beam element nodes and these virtual
nodes, and thus with all nodes of the soil volume element.
Pile forces (structural forces) are evaluated at the beam element integration points and
extrapolated to the beam element nodes. These forces can be viewed graphically and
tabulated in the Output program. Details about the embedded pile element formulations
are given in the Scientific Manual.
5.6.5 PLATES
Plates are structural objects used to model thin two-dimensional structures in the
ground with a significant flexural rigidity (bending stiffness). The creation of a plate
is similar to the creation of a geometry surface (Section 5.2.4).
Plate elements
After meshing plates are composed of 6-node triangular plate elements with six degrees
of freedom per node: three translational degrees of freedom (ux , uy and uz ) and three
rotational degrees of freedom (φx , φy and φz ). The plate elements are based on Mindlin's
plate theory (Bathe, 1982). This theory allows for plate deflections due to shearing as
well as bending. In addition, the element can change length when an axial force is
applied. Element stiffness matrices are based on the properties as defined in the material
data sets and numerically integrated from the three pairs of Gaussian integration
√ points
(stress points). Within each pair, stress points are located at a distance 1/6 3d above
and below the plate centre-line. Details about the element formulation are given in the
Scientific Manual.
Plate properties
The material properties of plates are contained in Plates material data sets (Section 6.4)
and can be conveniently assigned using drag-and-drop, the Object explorers or the right
hand mouse button menu in the draw area (Section 6.9).
Structural forces are evaluated at the plate element integration points (Scientific Manual)
and extrapolated to the element nodes. These forces can be viewed graphically and
tabulated in the Output program. Details about the material properties are given in
Section 6.4 and Material Models Manual.
5.6.6 GEOGRIDS
Geogrids are slender structures with an axial stiffness but with no bending stiffness.
Geogrids can only sustain tensile forces and no compression. These objects are
generally used to model soil reinforcements.
Geogrid elements
After meshing, geogrids are composed of 6-node triangular surface elements with three
translational degrees of freedom per node (ux , uy and uz ). Element stiffness matrices are
based on the properties as defined in the material data sets and numerically integrated
from the three Gaussian integration points (stress points). The element can change
length when a tension force is applied. Details about the element formulation are given in
the Scientific Manual.
When a geogrid element is connected to another structural element, they share the
translational degrees of freedom in the connecting node(s).
Geogrid properties
The basic material parameter is the axial stiffness EA. In addition the tension force can
be limited to allow for the simulation of tension failure. Axial forces can be viewed
graphically and tabulated in the Output program. Details about the material properties are
given in Section 6.5 and Material Models Manual.
5.6.7 INTERFACES
Interfaces are joint elements to be added to plates or geogrids to allow for a proper
modelling of soil-structure interaction. Interfaces may be used to simulate, for example,
the thin zone of intensely shearing material at the contact between a plate and the
surrounding soil. Interfaces can be created next to plate or geogrid elements or between
two soil volumes.
An interface can be created using the Create interface option available in
the menu displayed when the Create structure button in the side toolbar is clicked.
The creation of an interface is similar to the creation of a geometry surface (Section
5.2.4). In this case a surface with an interface assigned to it will be created. When
the geometric entity (surface) is already available in the model it is advised to assign an
interface to it without recreating the geometry in order to prevent the model from being
unnecessarily large and unwieldy. Assignment of features to existing geometric entities is
described in Section 5.2.6.
Distinction is made between a positive interface (the side of the surface at the positive
local z -direction) and a negative interface (the side of the surface at the negative local
z -direction).
Interface properties
The interfaces created in the model are listed in the Interfaces subtree in the Object
explorers. Properties such as Material mode, permeability condition and Virtual thickness
factor, available for each interface in the Object explorers can be assigned to an interface
in the model (Figure 5.31).
6.1.4.
Custom A material dataset can be assigned directly to an interface when
the Custom option is selected. Note that the strength reduction
factor of the selected material is by default set to one.
Permeability: By default structural elements such as Plates and Geogrids are
permeable. Impermeability can be introduced by assigning an interface to the geometric
entity, besided the structural element. Note that for an interface a checkbox is available in
the Objects explorers, enabling specifying the permeability condition (either fully
permeable or fully impermeable). In Groundwater flow only calculations no structural
element needs to be assigned to the geometric entity as deformations will not be taken
into account. Hence, it is sufficient to use only interfaces to block the flow.
Virtual interface thickness: Each interface has assigned to it a 'virtual thickness' which
is an imaginary dimension used to define the material properties of the interface. The
higher the virtual thickness is, the more elastic deformations are generated. In general,
interface elements are supposed to generate very little elastic deformations and therefore
the virtual thickness should be small. On the other hand, if the virtual thickness is too
small, numerical ill-conditioning may occur. The virtual thickness is calculated as the
Virtual thickness factor times the global element size. The global element size is
determined by the global coarseness setting for the mesh generation (Section 7.1.1). The
default value of the Virtual thickness factor is 0.1. This value can be changed in the
Object explorers. In general, care should be taken when changing the default factor.
However, if interface elements are subjected to very large normal stresses, it may be
required to reduce the Virtual thickness factor. Further details of the significance of the
virtual thickness are given in Section 6.1.4.
Interface elements
After meshing, interfaces are composed of 12-node interface elements. Interface
elements consist of pairs of nodes, compatible with the 6-noded triangular side of a soil
element or plate element. In some output plots, interface elements are shown to have a
finite thickness, but in the finite element formulation the coordinates of each node pair are
identical, which means that the element has a zero thickness.
The stiffness matrix for triangular interface elements is based on the properties as
defined in the material data sets and obtained by means of Gaussian integration using 6
integration points. The position of these integration points (or stress points) is chosen
such that the numerical integration is exact for linear stress distributions. For more details
about the element formulation reference is made to the Scientific Manual.
At interface ends interface element node pairs are 'degenerated ' to single nodes. Also
when structural elements are connected perpendicular to each other (for example, a plate
is connected to a beam), interface element node pairs are locally 'degenerated' to single
nodes to avoid a disconnection between both structural elements.
non-physical stress oscillations. This problem can be solved by making use of interface
elements as shown in Figure 5.33.
Figure 5.32 Inflexible corner point, causing poor quality stress results
Figure 5.33 shows that the problem of stress oscillation may be prevented by specifying
additional interface elements inside the soil body. These elements will enhance the
flexibility of the finite element mesh and will thus prevent non-physical stress results.
However, these elements should not introduce an unrealistic weakness in the soil.
Therefore special attention should be made to the properties of these interface elements
(Section 6.1.4).
At interface ends, interface element node pairs are 'degenerated' to single nodes. When
an interface exists on only one side of a structure, the interface end node is collapsed to
the corresponding structure node. As a result, there exists only one node at that location.
When interfaces are assigned on both sides of a structure the end nodes of the interfaces
are collapsed, however the collapsed node is not the corresponding structure node. As a
result, there exist two node at that specific location, the structure node and the node
shared by both interfaces. The interactions between the structure and the surrounding
soil at that node is controlled by a combined effect of both interfaces. Also when
structural elements are connected perpendicular to each other (for example, a plate is
connected to a beam), interface element node pairs are locally 'degenerated' to single
nodes to avoid a disconnection between both structural elements, connecting them
rigidly. A description on how to define a non-rigid connection is given in Section 7.9.4.
Additional theoretical details on this special use of interface elements is provided by
Goodman, Taylor & Brekke (1968) and van Langen & Vermeer (1991).
• Drag the tool (keep the left mouse button down) and drop it on the geometric entity.
Note that the assignment of hydraulic conditions is also possible by right-clicking the
geometric entity either in the draw area or in the Object explorers and selecting the
corresponding option in the appearing menu.
Although the input values of hydraulic conditions are specified in the geometry model, the
activation, deactivation or change of values may be considered in the framework of
Staged construction (Section 7.9.1).
In addition to the specific boundary conditions as considered here, global model
conditions for the model as a whole in terms of 'open' and 'closed' boundaries as well as
precipitation conditions can be defined for each calculation phase in the Model conditions
subtree in the Model explorer (Section 7.9.9). Note that specific hydraulic conditions have
priority over global model conditions.
For transient groundwater flow and fully coupled flow-deformation analysis, hydraulic
conditions can be defined as a function of time using so-called Flow functions (Section
5.8).
5.7.1 WELL
In order to define a well, the Create well option should be selected from the menu
appearing as the Create hydraulic conditions button is clicked. The creation of a
well in the geometry model is similar to the creation of a line (Section 5.2.2). Wells are
used to prescribe lines inside the geometry model where a specific flux (discharge) is
extracted from or infiltrated into the soil. The required properties of a well are:
Behaviour The behaviour of the well needs to be specified. The available
options are Extraction (taking water from the ground) and
Infiltration (bringing water into the ground).
|Qwell | The discharge of the well.
hmin The minimum possible head in the well. When the groundwater
head reduces below the hmin level no further extraction will
occur. As a habit, hmin may be set equal to the bottom level of
the well in the ground.
This option is only relevant for groundwater flow calculations and fully coupled
flow-deformation analysis. Wells can be activated or de-activated in calculation phases.
Hint: When a well intersects multiple soil layers the prescribed flux for each soil
layer is a function of the saturated permeability and the intersected depth.
» Note that the saturated permeability is specified in the Flow parameters
tabsheet of the material dataset.
clicked. The creation of a drain in the geometry model is similar to the creation of a
surface (Section 5.2.4). Surface drains are used to prescribe surfaces inside the
geometry model where (excess) pore pressures are reduced.
Together with the creation of a surface drain, the input of a groundwater head (h) is
required. This option is only relevant for consolidation analyses or groundwater flow
calculations. In such calculations, the pore pressure in all nodes of the drain is reduced
such that it is equivalent to the given head. Pore pressures lower than the equivalent to
the given head are not affected by the drain. Drains can be activated or de-activated in
calculation phases.
Figure 5.35 Options for the behaviour of Surface groundwater flow boundary conditions
Hint: The boundary conditions at the extremities of the model can be conveniently
specified in the Model conditions subtree in the Model explorer. By default
the bottom of the model is set to Closed, preventing flow across it, whereas
the three remaining boundaries are set to Open (Seepage).
» In a calculation phase the hydraulic conditions defined using the Surface flow
boundary conditions feature always prevail over the Model conditions. As an
example, if if a surface flow boundary conditions is placed on the top surface
and precipitation is specified in Model conditions only the behaviour specified
for the surface flow boundary condition will be considered, ignoring the
precipitation.
Seepage: A Seepage boundary is a boundary where water can flow in or out freely. It is
generally used at the ground surface above the phreatic level or above the external water
level.
If a surface flow boundary is set to Seepage and is completely above the (external) water
level, then the seepage condition applies to this boundary. This means that water inside
the geometry may flow freely out of this boundary.
If a boundary is Seepage and completely below the (external) water level, the free
boundary condition automatically turns into a groundwater head condition. In that case
the magnitude of the groundwater head in each boundary node is determined by the
vertical distance between the boundary node and the water level.
The pore pressure is zero where a (external) water level crosses a geometry boundary
surface. The part of the geometry surface above the transition line is treated as a
boundary above the water level, whereas the part of the geometry surface below the
transition line is treated as a boundary below the water level. Hence, different conditions
can apply to such a geometry boundary surface. This is possible because, in general, a
geometry surface consists of many nodes and the actual information on boundary
conditions as used by the calculation program is contained in the boundary nodes rather
than in geometry surfaces.
Hint: Note that a boundary that is explicitly set to Seepage will NOT automatically
turn into an Infiltration boundary when precipitation is specified.
Closed: When the Closed option is assigned to a boundary, no flow can occur across it.
Flow, in this case, refers to groundwater flow (in groundwater flow calculations and fully
coupled flow-deformation analysis) as well as dissipation of excess pore pressure (in
consolidation calculations).
Head: In addition to the automatic setting of hydraulic conditions based on the general
phreatic level in Model conditions (Section 7.8.4), a prescribed groundwater head may be
entered manually to the Surface groundwater flow boundary conditions defined by the
user in the model. If a groundwater head is prescribed at an outer geometry boundary,
external water pressures will be generated for that boundary. The deformation analysis
program will treat external water pressures as traction loads and they are taken into
account together with the soil weight and the pore pressures. Figure 5.36 shows the
options available for the definition of head in Object explorers.
The options available for variation of head along the boundary are:
Constant A constant value of the head (href ) is specified along the
boundary.
Linear A linear variation of the head along the boundary is specified.
The parameters required to specify a linear variation are:
href The reference head value at the
specified reference level.
hinc,z Increment value specifying the change
of the head with depth.
The options available for Time dependency of Head assigned to a surface flow boundary
condition are:
Constant The head in the boundary does not change with time
Time dependent The head in the boundary varies with time. This option is only
available for transient flow and fully coupled flow-deformation
analysis. The flow function representing the time dependency
can be selected from the drop-down menu. Information on Flow
functions is available in Section 5.8.
Inflow: Inflow to the model through a boundary can be specified by selecting the
corresponding Behaviour option in the Object explorers. Figure 5.37 shows the options
available for the definition of inflow.
The options available for the distribution of the inflow discharge along the boundary are:
Constant A constant value of the discharge is specified along the
boundary. The value is assigned to the |qref | parameter.
Linear A linear variation of the discharge along the boundary is
specified. The parameters required to specify a linear variation
are:
|qref | The reference head value at the
specified reference level.
|qinc,z | Increment value specifying the rate of
variation of discharge along the
boundary
The options available for variation of discharge with time are (Time dependency):
Constant The discharge does not change with time
Time dependent The discharge in the boundary varies with time. This option is
only available for transient flow and fully coupled
flow-deformation analysis. The flow function representing the
time dependency can be selected from the drop-down menu.
Information on Flow functions is available in Section 5.8.
Infiltration: Apart from the automatic generation of infiltration boundary conditions from
precipitation (Section 7.8.4), infiltration conditions may also be specified manually for
geometry boundaries above the water level. Infiltration boundary condition is a
conditional inflow. Figure 5.39 shows the properties of infiltration in Object explorers. The
parameters required for an infiltration boundary are:
q Recharge (infiltration), specified in the unit of length per unit of
time. Negative values can be used to model evapotranspiration
(evaporation + transpiration).
ψmax Maximum pore pressure head, relative to the elevation of the
boundary, specified in the unit of length (default 0.1 length units).
ψmin Minimum pore pressure head, relative to the elevation of the
boundary, specified in the unit of length (default -1.0 length
units).
The options available for variation of discharge with time are (Time dependency):
Constant The discharge does not change with time
Time dependent The discharge in the boundary varies with time. This option is
only available for transient flow and fully coupled
flow-deformation analysis. The flow function representing the
time dependency can be selected from the drop-down menu.
Information on Flow functions is available in Section 5.8.
Flow functions describing the variation with time of quantities such as Head and
Discharge can be defined in the corresponding subtree under the Attributes library in the
Model explorer (Figure 5.40).
Flow functions can be assigned to Water levels or Surface groundwater flow boundary
conditions. Note that time-dependent conditions can be defined for all behaviour options
except for Seepage and Closed.
Functions can be defined in the corresponding window that pops up when the Flow
functions subtree is right-clicked and the Edit option is selected in the appearing menu.
The window consists of two tabsheets, i.e Head functions and Discharge functions, where
the functions applicable to head and discharge can be defined respectively (Figure 5.41).
The buttons available under the function tabs can be used to add new functions or to
delete the function selected in the list. When a new function is added, the options to
define it are displayed.
Name The name of the flow function can be defined.
Signal The type of the signal can be assigned, The options available
are Linear, Harmonic and Table.
The options available in the Signal drop-down menu are described in the following
sections.
5.8.1 HARMONIC
This option can be used when a condition varies harmonically in time. The harmonic
variation of the water level is generally described as:
in which A is the amplitude (in unit of length), T is the wave period (in unit of time) and φ0
is the initial phase angle. y0 is the centre of the harmonic variation, not necessarily equal
to the water level at the end of the previous phase (y '0 ). In fact the later is equal to:
y '0 = y0 + A sin(ω0 )
Figure 5.42 shows the Flow functions window where a Harmonic signal is defined and
plotted.
Figure 5.42 Definition and display of a Harmonic signal in the Flow functions window
Hint: PLAXIS assumes the data file is located in the current project directory when
no directory is specified.
A signal, either defined in a table or received from a file and modified, can be
saved using the Save button in the toolbar enabling the usage of the signal in other
projects or validating the effect of the modifications in the current project.
A signal, either defined in a table or received
from a file and modified, can be copied using the Copy button in the toolbar.
Copied data from other applications (using <Ctrl+c>) can be imported by using the
Paste button. The Import data window appears (Figure 5.43). The starting row of
the data to be imported can be defined in the From row cell. The data and the plot is
displayed in the Flow functions window after pressing OK.
Clicking the Open .txt file button on the right hand side of the window will open the
Open window where the file can be selected. The file must be an ASCII file that can
be created with any text editor. For every line a pair of values (actual time and
corresponding water level value) must be defined, leaving at least one space between
them. Note that PLAXIS only supports the English notation of decimal numbers using a
dot. The resulting graph of the input data is shown in the Graph tabsheet of the Time
dependent condition window.
5.8.3 LINEAR
This option can be used to describe the increase or decrease of a condition linearly in
time. For a linear variation of groundwater head, the input of the following parameters are
required:
∆Head This parameter, specified in unit of length, represents the
increase or decrease of the water level in the time interval for the
current calculation phase. Hence, together with the time interval
this parameter determines the rate of the water level increase or
decrease.
t This parameter represents the time interval for the calculation
phase, expressed in unit of time. Note that the time interval
Hint: The time value in flow functions always refers to the global time in the whole
list of calculation phases rather than to the time interval of an individual
phase. This means that in a series of consecutive flow calculations each
phase will only use its corresponding part of the flow function.
Figure 5.44 Options for boundary condition behaviour in Phase definition modes
PLAXIS 3D allows to import geometry. This can be done by using the Import soil volume
button in the Soil mode or the Import volume or the top and bottom soil surfaces can be
imported in the Surfaces tabsheet of the Modify soil layers window.
In the Structures mode, geometry can be imported by clicking the Import structures
button available in the side toolbar. Note that only surfaces and volumes can be
imported. When importing a geometry file, a list of possible file extensions are shown.
Most common import file types are:
AutoCAD native (*.DWG) and interchange (*.DXF)
For DWG and DXF files only the 3DFACE entities are supported
(i.e. surfaces/volume boundaries made up of triangles or
quadrilaterals), other entity types are ignored.
3D Studio files (*.3DS) Valid surfaces or volumes need to be defined.
In both cases a file requester pops up from which the desired surface or volume can be
selected.
After selecting the file, the Import structures window pops up (Figure 5.45), where the
geometry can be modified. It is possible to import only surfaces or volumes by keeping
Scaling
The imported geometry can be scaled by defining factors for each of the global
directions. The Keep original aspect ratio option enables applying the same scale to all
the directions. In this case, only the scale value for the x-direction can be changed.
Axis transformation
The orientation of the imported geometry might need to be modified in order to be
properly mounted in the project.
Source vertical axis It enables reorienting the imported geometry such that the
vertical axes of the imported geometry and the program
coincide. The source vertical axis enables axis transformation
according to the right hand rule of rotation. Note that the vertical
direction in PLAXIS 3D is the z-axis.
Mirror y-axis The geometry will be transformed such that the values of y
coordinates will be changed to −y coordinates (effectively in the
xz-plane).
Swap x-axis and y-axis The geometry will be transformed such that the values of x
coordinates will be changed to y coordinates and vice versa.
Insertion point
When a geometric object is imported, its reference point is located at (0 0 0) in the model.
Redefining the location of the insertion point will relocate the reference point.
Bounding box
The coordinates of the box bounding the imported geometry are shown.
Hint: Material properties of volumes or surfaces are not imported and need to be
created in the material data base.
» A description of the feature and an illustration movie is given in the
Knowledge Base section of the PLAXIS website (www.plaxis.nl).
In PLAXIS, soil properties and material properties of structures are stored in material
data sets. There are six different types of material sets grouped as data sets for soil and
interfaces, beams, embedded piles, plates, geogrids anchors. All data sets are stored in
the material database. From the database, the data sets can be assigned to the soil
clusters or to the corresponding structural objects in the geometry model.
The material database can be activated by selecting the option Show materials
in the Soil menu of the Soil mode or in the Structures menu of the Structures mode.
Another way to activate the material database is by clicking on the Show materials button
available in the Soil, Structures and Staged construction modes.
As a result, the Material sets window appears showing the contents of the project
material database. The window can be extended to show the global database by clicking
the Show global button in the upper part of the window. The Material sets window
displaying the material defined in the current project and the ones available in a selected
global database is shown in Figure 6.1.
Figure 6.1 Material sets window showing the project and the global database
The database of a new project is empty. The global database can be used to store
material data sets in a global folder and to exchange data sets between different projects.
At both sides of the window (Project materials and Global materials) there are two
drop-down menus and a tree view. The Set type can be selected from the drop-down
menu on the left hand side. The Set type parameter determines which type of material
data set is displayed in the tree view (Soil and interfaces, Plates, Geogrids, Beams,
Embedded piles, Anchors).
The data sets in the tree view are identified by a user-defined name. The data sets for
Soil and interfaces can be ordered in groups according to the material model, the
material type or the name of the data set by selecting this order in the Group order
drop-down menu. The None option can be used to discard the group ordering.
The small buttons between the two tree views can be used to copy individual data sets
from the project database to the selected global database or vice versa.
To copy the selected project material set to the global database.
To copy all the project material sets of the specified type to the global database.
To copy the selected global material set to the project database.
The location of the selected global database is shown below its tree view. The buttons
below the tree view of the global database enable actions in the global database.
Select To select an existing global database.
Delete To delete a selected material data set from the selected global
database.
By default, the global database for soil and interface data contains the data sets of all the
tutorials and it is contained in the file 'SoilMat.matdb'. This file is compatible with other
PLAXIS database files for soil and interfaces and is stored in the installation folder of
PLAXIS 3D. Material data sets for structural elements will be contained in separate files.
Similarly, the global data bases for plates, geogrids, beams, embedded piles and anchors
are contained in the files 'PlateMat3D.matdb', 'GeogridMat.matdb', 'BeamMat.matdb',
'EmbeddedPileMat.matdb' and 'AnchorMat3D.matdb' respectively.
Note that besides the global material files (*.matdb), it is possible to select project
material files (*.plxmat) and legacy project material files (*.mat) as global database.
In addition, databases with data sets of standard sheet-pile wall profiles are available
from the Plaxis Knowledge Base
(http://kb.plaxis.nl/downloads/material-parameter-datasets-sheetpiles-and-beams).
Hint: A new global database can be created by clicking the Select button, defining
the name of the new global database and clicking Open.
The project data base can be managed using the buttons below the tree view of the
project database.
New To create a new data set in the project. As a result, a new
window appears in which the material properties or model
parameters can be entered. The first item to be entered is
always the Identification, which is the user-defined name of the
data set. After completing a data set, it will appear in the tree
view, indicated by its name as defined by the Identification.
Edit To modify the selected data set in the project material database.
SoilTest To perform standard soil lab tests. A separate window will open
where several basic soil tests can be simulated and the
behaviour of the selected soil material model with the given
material parameters can be checked (Section 6.3).
Copy To create a copy of a selected data set in the project material
database.
Delete To delete a selected material data set from the project material
database.
The material properties and model parameters for soil clusters are entered in material
data sets (Figure 6.2). The properties in the data sets are divided into five tabsheets:
General, Parameters, Flow parameters, Interfaces and Initial.
Material model
Soil and rock tend to behave in a highly non-linear way under load. This non-linear
stress-strain behaviour can be modelled at several levels of sophistication. Clearly, the
number of model parameters increases with the level of sophistication. PLAXIS supports
different models to simulate the behaviour of soil and other continua. The models and
their parameters are described in detail in the Material Models Manual. A short
discussion of the available models is given below:
Linear elastic model: This model represents Hooke's law of isotropic linear elasticity.
The linear elastic model is too limited for the simulation of soil behaviour. It is primarily
used for stiff structures in the soil.
Mohr-Coulomb model (MC): This well-known linear elastic perfectly-plastic model is
used as a first approximation of soil behaviour in general. It is recommended to use this
model for a first analysis of the problem considered. A constant average stiffness is
estimated for the soil layer. Due to this constant stiffness, computations tend to be
relatively fast and a first estimate of deformations can be obtained.
Hardening Soil model (HS): This is an advanced model for the simulation of soil
behaviour. The Hardening Soil model is an elastoplastic type of hyperbolic model,
formulated in the framework of shear hardening plasticity. Moreover, the model involves
compression hardening to simulate irreversible compaction of soil under primary
compression. This second-order model can be used to simulate the behaviour of sands
and gravel as well as softer types of soil such as clays and silts.
Hardening Soil model with small-strain stiffness (HSsmall): This is an elastoplastic
type of hyperbolic model, similar to the Hardening Soil model. Moreover, this model
incorporates strain dependent stiffness moduli, simulating the different reaction of soils
from small strains (for example vibrations with strain levels below 10-5 ) to large strains
(engineering strain levels above 10-3 ).
Soft Soil model (SS): This is a Cam-Clay type model that can be used to simulate the
behaviour of soft soils like normally consolidated clays and peat. The model performs
best in situations of primary compression.
Soft Soil Creep model (SSC): This is a second order model formulated in the
framework of viscoplasticity. The model can be used to simulate the time-dependent
behaviour of soft soils like normally consolidated clays and peat. The model includes
logarithmic primary and secondary compression.
Jointed Rock model (JR): This is an anisotropic elastic-perfectly plastic model where
plastic shearing can only occur in a limited number of shearing directions. This model
can be used to simulate the anisotropic behaviour of stratified or jointed rock.
Modified Cam-Clay model (MCC): This well-known critical state model can be used to
simulate the behaviour of normally consolidated soft soils. The model assumes a
logarithmic relationship between the volumetric strain and the mean effective stress.
NGI-ADP model (NGI-ADP): The NGI-ADP model may be used for capacity,
deformation and soil-structure interaction analysis involving undrained loading of clay.
Distinct anisotropic stress strengths may be defined for different stress paths.
Hoek-Brown model (HB): This well-known elastic perfectly-plastic model is used to
simulate the isotropic behaviour of rock. A constant stiffness is used for the rock mass.
Shear failure and tension failure are described by a non-linear stress curve.
User-defined soil models (UDSM): With this option it is possible to use other
constitutive models than the standard PLAXIS models. For a detailed description of this
facility, reference is made to the Material Models Manual. Links to existing User-defined
soil models are available on the Plaxis Knowledge Base (http://kb.plaxis.nl/models).
Drainage type
In principle, all model parameters in PLAXIS are meant to represent the effective soil
response, i.e. the relationship between the stresses and the strains associated with the
soil skeleton. An important feature of soil is the presence of pore water. Pore pressures
significantly influence the (time-dependent) soil response. PLAXIS offers several options
to enable incorporation of the water-skeleton interaction in the soil response. The most
advanced option is a Fully coupled flow-deformation analysis. However, in many cases it
is sufficient to analyse either the long-term (drained) response or the short-term
(undrained) response without considering the time-dependent development of pore
pressures. In the latter case (undrained), excess pore pressures are generated as a
result of stress changes (loading or unloading). The dissipation of these excess pore
pressures with time can be analysed in a Consolidation calculation.
The simplified water-skeleton interaction, as considered in a Plastic calculation, a Safety
analysis or a Dynamic analysis, is defined by the Drainage type parameter. PLAXIS
offers a choice of different types of drainage:
Drained behaviour: Using this setting no excess pore pressures are generated. This is
clearly the case for dry soils and also for full drainage due to a high permeability (sands)
and/or a low rate of loading. This option may also be used to simulate long-term soil
behaviour without the need to model the precise history of undrained loading and
consolidation.
Undrained behaviour: This setting is used for saturated soils in cases where pore
water cannot freely flow through the soil skeleton. Flow of pore water can sometimes be
neglected due to a low permeability (clays) and/or a high rate of loading. All clusters that
are specified as undrained will indeed behave undrained, even if the cluster or a part of
the cluster is located above the phreatic level.
Distinction is made between three different methods of modelling undrained soil
behaviour. Method A is an undrained effective stress analysis with effective stiffness as
well as effective strength parameters. This method will give a prediction of the pore
pressures and the analysis can be followed by a consolidation analysis. The undrained
shear strength (su ) is a consequence of the model rather than an input parameter. It is
recommended to check this shear strength with known data. To consider this type of
analysis, the Undrained (A) option should be selected in the Drainage type drop-down
menu.
Method B is an undrained effective stress analysis with effective stiffness parameters and
undrained strength parameters. The undrained shear strength su is an input parameter.
This method will give a prediction of pore pressures. However, when followed by a
consolidation analysis, the undrained shear strength (su ) is not updated, since this is an
input parameter. To consider this type of analysis, the Undrained (B) option should be
selected in the Drainage type drop-down menu.
Method C is an undrained total stress analysis with all parameters undrained. This
method will not give a prediction of pore pressures. Therefore it is not useful to perform a
Hint: The Drainage type setting is only considered in a Plastic calculation, a Safety
analysis or a Dynamic analysis. When a Consolidation analysis or a Fully
coupled flow-deformation analysis is performed, the Drainage type is ignored
and the soil response is determined by the Permeability of the material.
C = αM + βK
The parameters α and β are the Rayleigh coefficients and can be specified in the
corresponding cells in the Parameters tabsheet of the Soil window (Figure 6.3).
α is the parameter that determines the influence of mass in the damping of the system.
The higher α is, the more the lower frequencies are damped. β is the parameter that
determines the influence of stiffness in the damping of the system. The higher β is, the
more the higher frequencies are damped. In PLAXIS 3D, these parameters can be
specified for each material data set for soil and interfaces as well as for material data sets
for plates. In this way, the (viscous) damping characteristics can be specified for each
individual material in the finite element model.
Despite the considerable amount of research work in the field of dynamics, little has been
achieved yet for the development of a commonly accepted procedure for damping
parameter identification. Instead, for engineering purposes, some measures are made to
account for material damping. A commonly used engineering parameter is the damping
ratio ξ . The damping ratio is defined as ξ = 1 for critical damping, i.e. exactly the amount
of damping needed to let a single degree-of-freedom system that is released from an
initial excitation u0 , smoothly stop without rebouncing.
Considering Rayleigh damping, a relationship can be established between the damping
ratio ξ and the Rayleigh damping parameters α and β :
α + β ω 2 = 2 ω ξ and ω = 2 π f
Overdamped (ξ > 1)
Critically damped (ξ = 1)
Underdamped (ξ < 1)
Figure 6.4 Role of damping ratio ξ in free vibration of a single degree-of-freedom system
Solving this equation for two different target frequencies and corresponding target
damping ratios gives the required Rayleigh damping coefficients:
ω1 ξ2 − ω2 ξ1 ω1 ξ1 − ω2 ξ2
α = 2ω1 ω2 and β = 2
ω12 − ω22 ω12 − ω22
For example, when it is desired to have a target damping of 8% at the target frequencies
f = 1.5 Hz and 8.0 Hz, the corresponding Rayleigh damping ratios are α = 1.2698 and
β = 0.002681. From Figure 6.5 it can be seen that within the range of frequencies as
defined by the target frequencies the damping is less than the target damping, whereas
outside this range the damping is more than the target damping.
0.5
0.4
1.5 8.0
Damping ratio (-)
0.3
0.2
Damping curve
0.1
8%
Influence of β Influence of α
0
0.1 1 10 100
Frequency (Hz)
The damping parameters (α and β ) can be automatically calculated by the program when
the target damping ratio (ξ ) and the target frequencies (f) are specified in the pane
displayed in the General tabsheet when one of the cells corresponding to the damping
parameters is clicked (Figure 6.6). A graph shows the damping ratio as a function of the
frequency.
Hint: Optional drainage types when the Linear Elastic model is selected are:
Drained, Undrained (A), Undrained (C), and Non-porous.
» In the case of Undrained (A) or Non-porous drainage types, the same
parameters are used as for drained behaviour.
» In the case of Undrained (C) drainage type, an undrained Young's modulus
(Eu ) and undrained Poisson's ratio (νu ) are used.
The model involves two elastic stiffness parameters, namely the effective Young's
modulus E ' and the effective Poisson's ratio ν '.
E' : Effective Young's modulus [kN/m2 ]
ν' : Effective Poisson's ratio [-]
During the input for the Linear Elastic model the values of the shear modulus G and the
oedometer modulus Eoed are presented as auxiliary parameters (alternatives).
G : Shear modulus, where G = E' [kN/m2 ]
2(1 + ν ')
E '(1 − ν ')
Eoed : Oedometer modulus, where Eoed = [kN/m2 ]
(1 + ν ')(1 − 2ν ')
Figure 6.7 Parameters tabsheet for the Linear Elastic model (drained behaviour)
Note that the alternatives are influenced by the input values of E ' and ν '. Entering a
particular value for one of the alternatives G or Eoed results in a change of the Young's
modulus E '.
It is possible for the Linear Elastic model to specify a stiffness that varies linearly with
depth. Therefore, the increment of stiffness per unit of depth, E 'inc , can be defined.
Together with the input of E 'inc the input of zref becomes relevant. For any z -coordinate
above zref the stiffness is equal to E 'ref . For any z -coordinate below zref the stiffness is
given by:
The Linear Elastic model is usually inappropriate to model the highly non-linear behaviour
of soil, but it is of interest to simulate structural behaviour, such as thick concrete walls or
plates, for which strength properties are usually very high compared with those of soil.
For these applications, the Linear Elastic model will often be selected together with
Non-porous type of material behaviour in order to exclude pore pressures from these
structural elements.
Hint: When embedded piles penetrate a volume cluster with linear elastic material
behaviour, the specified value of the shaft resistance is ignored. The reason
for this is that the linear elastic material is not supposed to be soil, but part of
the structure. The connection between the pile and the structure is supposed
to be rigid to avoid, for example, punching of piles through a concrete deck.
Beside the parameters related to strength and stiffness of the soil, the velocities of wave
propagation in soil can be defined in the Parameters tabsheet of the Soil window when
the Dynamics module of the program is available. These velocities are:
p
Vs : Shear wave velocity, where Vs = G/ρ [m/s]
p
Vp : Compression wave velocity, where Vp = Eoed /ρ [m/s]
where ρ = γ/g .
Hint: Note that the wave velocities are influenced by the input values of E ' and ν '.
Entering a particular value for one of the wave velocities results in a change
of the Young's modulus.
» Velocities of wave propagation in soil can be defined only for models with
stress independent stiffness.
Figure 6.8 Parameters tabsheet for the Mohr-Coulomb model (drained behaviour)
for E 'inc which is the increment of stiffness per unit of depth. Together with the input of
E 'inc the input of zref becomes relevant. For any z -coordinate above zref the stiffness is
equal to E 'ref . For any z -coordinate below zref the stiffness is given by:
In some practical problems an area with tensile stresses may develop. This is allowed
when the shear stress is sufficiently small. However, the soil surface near a trench in clay
sometimes shows tensile cracks. This indicates that soil may also fail in tension instead
of in shear. Such behaviour can be included in a PLAXIS analysis by selecting the
Tension cut-off option. When selecting the Tension cut-off option the allowable tensile
strength may be entered. The default value of the tensile strength is zero.
Hardening Soil model (HS): The Parameters tabsheet for the Hardening Soil model is
shown in Figure 6.9.
Hint: Optional drainage types when Hardening Soil model is selected are:
Drained, Undrained (A), and Undrained (B).
» In the case of Undrained (A) drainage type, the same parameters are used
as for drained behaviour.
» In the case of Undrained (B) drainage type, ϕ = ϕu = 0 , ψ = 0 and the
undrained shear strength su is used instead of the effective cohesion (c').
Figure 6.9 Parameters tabsheet for the Hardening Soil model (drained behaviour)
Instead of entering the basic parameters for soil stiffness, alternative parameters can be
entered. These parameters are listed below:
Cc : Compression index [-]
Cs : Swelling index or reloading index [-]
einit : Initial void ratio [-]
In addition, advanced parameters can be defined for stiffness (it is advised to use the
default setting):
νur : Poisson's ratio for unloading-reloading (default ν = 0.2) [-]
pref : Reference stress for stiffnesses (default pref = 100 [kN/m2 ]
kN/m2 )
K0nc : K0 -value for normal consolidation (default K0nc = [-]
1 − sin ϕ)
The strength parameters of the present hardening model coincide with those of the
non-hardening Mohr-Coulomb model:
c 'ref : Effective cohesion [kN/m2 ]
ϕ' : Effective angle of internal friction [◦ ]
ψ : Angle of dilatancy [◦ ]
In addition, advanced parameters can be defined for strength:
c 'inc : As in Mohr-Coulomb model (default cinc = 0) [kN/m3 ]
zref : Reference level [m]
Rf : Failure ratio qf / qa (default Rf = 0.9) [-]
Hint: Optional drainage types when Hardening Soil model with small-strain
stiffness is selected are: Drained, Undrained (A), and Undrained (B).
» In the case of Undrained (A) drainage type, the same parameters are used
as for drained behaviour.
» In the case of Undrained (B) drainage type, ϕ = ϕu = 0, ψ = 0 and the
undrained shear strength su is used instead of the effective cohesion (c').
Figure 6.10 Parameters tabsheet for the HS small model (drained behaviour)
All other parameters, including the alternative stiffness parameters, remain the same as
in the standard Hardening Soil model. In summary, the input stiffness parameters of the
HS small model are listed below:
Parameters for stiffness:
ref
E50 : Secant stiffness in standard drained triaxial test [kN/m2 ]
ref
Eoed : Tangent stiffness for primary oedometer loading [kN/m2 ]
ref
Eur : unloading / reload stiffness at engineering strains [kN/m2 ]
( ≈ 10−3 to 10−2 )
m : Power for stress-level dependency of stiffness [-]
Alternative parameters for stiffness:
Cc : Compression index [-]
Cs : Swelling index or reloading index [-]
einit : Initial void ratio [-]
Advanced parameters for stiffness:
νur : Poisson's ratio for unloading-reloading (default ν = 0.2) [-]
p ref
: Reference stress for stiffnesses (default p ref
= 100 [kN/m2 ]
kN/m2 )
K0nc : K0 -value for normal consolidation (default K0nc = [-]
1 − sin ϕ)
Parameters for strength:
c 'ref : Effective cohesion [kN/m2 ]
ϕ' : Effective angle of internal friction [◦ ]
ψ : Angle of dilatancy [◦ ]
Advanced parameters for strength:
c 'inc : As in Mohr-Coulomb model (default c 'inc = 0) [kN/m3 ]
zref : Reference level [m]
Rf : Failure ratio qf / qa (default Rf = 0.9) [-]
Tension : To be selected when tension cut-off is considered [-]
cut-off
Tensile : The allowable tensile strength [kN/m2 ]
strength
Parameters for small strain stiffness:
γ0.7 : shear strain at which Gs = 0.722G0 [-]
G0ref : reference shear modulus at very small strains [kN/m2 ]
( < 10−6 )
In some practical problems an area with tensile stresses may develop. This is allowed
when the shear stress is sufficiently small. However, the soil surface near a trench in clay
sometimes shows tensile cracks. This indicates that soil may also fail in tension instead
Hysteretic damping
The elastic modulus ratio is plotted as a function of the shear strain (γ ) in a side pane
when specifying the small-strain stiffness parameters (Modulus reduction curve). The HS
small model shows typical hysteretic behaviour when subjected to cyclic shear loading. In
dynamic calculations this leads to hysteretic damping. The damping ratio is plotted as a
function of the cyclic shear strain γc . Details are given in Brinkgreve, Kappert & Bonnier
(2007).
Hint: Note that the Modulus reduction curve and the Damping curve are based on
fully elastic behaviour. Plastic strains as a result of hardening or local failure
may lead to significant lower stiffness and higher damping.
Soft Soil model (SS): The Parameters tabsheet for the Soft Soil model is shown in
Figure 6.12.
Hint: Optional drainage types when Soft Soil model is selected are: Drained and
Undrained (A).
» In the case of Undrained (A) drainage type, the same parameters are used
as for drained behaviour.
Figure 6.12 Parameters tabsheet for the Soft Soil model (drained behaviour)
when the shear stress is sufficiently small. However, the soil surface near a trench in clay
sometimes shows tensile cracks. This indicates that soil may also fail in tension instead
of in shear. Such behaviour can be included in a PLAXIS analysis by selecting the
Tension cut-off option. When selecting the Tension cut-off option the allowable tensile
strength may be entered. The default value of the tensile strength is zero.
Soft Soil Creep model (SSC): The Parameters tabsheet for the Soft Soil Creep model
is shown in Figure 6.13.
Hint: Optional drainage types when Soft Soil Creep model is selected are: Drained
and Undrained (A).
» In the case of Undrained (A) drainage type, the same parameters are used
as for drained behaviour.
Figure 6.13 Parameters tabsheet for the Soft Soil Creep model (drained behaviour)
Hint: Optional drainage types when Jointed Rock model is selected are: Drained
and Non-porous.
» In the case of Non-porous drainage type, the same parameters are used as
for drained behaviour.
Figure 6.14 Parameters tabsheet for the Jointed Rock model (drained behaviour)
ϕi : Friction angle [◦ ]
ψi : Dilatancy angle [◦ ]
σt,i : Tensile strength [kN/m2 ]
Definition of joint directions (Plane i=1, 2, 3):
n : Number of joint directions (1 ≤ n ≤ 3) [-]
α1,i : Dip angle [◦ ]
α2,i : Dip direction [◦ ]
In some practical problems an area with tensile stresses may develop. This is allowed
when the shear stress is sufficiently small. However, the soil surface near a trench in clay
sometimes shows tensile cracks. This indicates that soil may also fail in tension instead
of in shear. Such behaviour can be included in a PLAXIS analysis by selecting the
Tension cut-off option. When selecting the Tension cut-off option the allowable tensile
strength may be entered. The default value of the tensile strength is zero.
Modified Cam-Clay model (MCC): This is a critical state model that can be used to
simulate the behaviour of normally consolidated soft soils. The model assumes a
logarithmic relationship between the volumetric strain and the mean effective stress. The
Parameters tabsheet for the Modified Cam-Clay model is shown in Figure 6.15.
Parameters for stiffness:
λ : Cam-Clay compression index [-]
κ : Cam-Clay swelling index [-]
ν : Poisson’s ratio [-]
Hint: Optional drainage types when Modified Cam-Clay model is selected are:
Drained and Undrained (A).
» In the case of Undrained (A) drainage type, the same parameters are used
as for drained behaviour.
Figure 6.15 Parameters tabsheet for the Modified Cam-Clay model (drained behaviour)
NGI-ADP model (NGI-ADP): The NGI-ADP model may be used for capacity,
deformation and soil-structure interaction analysis involving undrained loading of clay.
The Parameters tabsheet for the NGI-ADP model is shown in Figure 6.16.
Hint: Optional drainage types when NGI-ADP model is selected are: Drained,
Undrained (B) and Undrained (C).
» In the case of Undrained (B) drainage type, the same parameters are used
as for drained behaviour.
A
su,inc : Increase of shear strength with depth [kN/m2 /m]
suP /suA : Ratio of (plane strain) passive shear strength over [-]
(plane strain) active shear strength
τ0 /suA : Initial mobilization (default = 0.7) [-]
suDSS /suA : Ratio of direct simple shear strength over (plane strain) [-]
active shear strength
Advanced parameters:
ν' : Effective Poisson's ratio [-]
νu : Undrained Poisson's ratio [-]
Hoek-Brown model (HB): The Parameters tabsheet for the Hoek-Brown model is
shown in Figure 6.17.
The stiffness parameters of the Hoek-Brown model are:
E : Young's modulus [kN/m2 ]
ν : Poisson's ratio [-]
The Hoek-Brown parameters are:
Hint: Optional drainage types when Hoek-Brown model is selected are: Drained
and Non-porous.
» In the case of Non-porous drainage type, the same parameters are used as
for drained behaviour.
Figure 6.17 Parameters tabsheet for the Hoek-Brown model (drained behaviour)
Hint: Available drainage types when User-defined soil models is selected are:
Drained, Undrained (A) and Non-porous.
Hint: Although the predefined data sets have been created for the convenience of
the user, the user remains at all times responsible for the model parameters
that he/she uses. Note that these predefined data sets have limited accuracy.
Hypres: The Hypres series is an international soil classification system. The hydraulic
models available for Hypres data set are the Van Genuchten model and the Approximate
Van Genuchten (see Sections 16.1 and 16.2 of the Material Models Manual).
A distinction can be made between Topsoil and Subsoil. In general, soils are considered
to be subsoils. The Type drop-down menu for the Hypres data set includes Coarse,
Medium, Medium fine, Fine, Very fine and Organic soils.
Hint: Only soil layers that are located not more than 1 m below the ground surface
are considered to be Upper soils.
The selected soil type and grading (particle fractions) is indicated in the soil texture
triangle. As an alternative, the user can also select the type of soil by clicking one of the
sections in the triangle or by manually specifying the particle fraction values (Figure 6.19).
The predefined parameters for both the Van Genuchten model as well as the
Approximate Van Genuchten model are shown in Table 6.1 and 6.2.
Table 6.1 Hypres series with Van Genuchten parameters
θr (-) θs (-) Ksat (m/day) ga (1/m) gl (-) gn (-)
Topsoil:
coarse 0.025 0.403 0.600 3.83 1.2500 1.3774
medium 0.010 0.439 0.121 3.14 -2.3421 1.1804
medium fine 0.010 0.430 0.0227 0.83 -0.5884 1.2539
fine 0.010 0.520 0.248 3.67 -1.9772 1.1012
very fine 0.010 0.614 0.150 2.65 2.5000 1.1033
Subsoil:
coarse 0.025 0.366 0.700 4.30 1.2500 1.5206
medium 0.010 0.392 0.108 2.49 -0.7437 1.1689
medium fine 0.010 0.412 0.0400 0.82 0.5000 1.2179
fine 0.010 0.481 0.0850 1.98 -3.7124 1.0861
very fine 0.010 0.538 0.0823 1.68 0.0001 1.0730
organic 0.010 0.766 0.0800 1.30 0.4000 1.2039
USDA: The USDA series is another international soil classification system. The
hydraulic models available for USDA data set are the Van Genuchten model and the
Approximate Van Genuchten (see Sections 16.1 and 16.2 of the Material Models
Manual).
The Type drop-down menu for the USDA date set includes Sand, Loamy sand, Sandy
loam, Loam, Silt, Silt loam, Sandy clay loam, Clay loam, Silty clay loam, Sandy clay, Silty
clay and Clay. The selected soil type and grading (particle fractions) are different from
the Hypres data sets and can be visualised in the soil texture triangle. As an alternative,
the user can also select the type of soil by clicking one of the sections in the triangle or by
manually specifying the particle fraction values (Figure 6.20).
The parameters for the Van Genuchten and the Approximate Van Genuchten models are
shown in Table 6.3 and 6.4.
Staring: The Staring series is a soil classification system which is mainly used in The
Netherlands. The hydraulic models available for Staring data set are the Van Genuchten
model and the Approximate Van Genuchten (see Sections 16.1 and 16.2 of the Material
Table 6.3 USDA series with Van Genuchten parameters (gl = 0.5 for all sets)
θr (-) θs (-) Ksat (m/day) ga (1/m) gn (-)
sand 0.045 0.430 7.13 14.5 2.68
loamy sand 0.057 0.410 3.50 12.4 2.28
sandy loam 0.065 0.410 1.06 7.5 1.89
loam 0.078 0.430 0.250 3.6 1.56
silt 0.034 0.460 0.600 1.6 1.37
silty loam 0.067 0.450 0.108 2.0 1.41
sandy clay loam 0.100 0.390 0.314 5.9 1.48
clayey loam 0.095 0.410 0.624 1.9 1.31
silty clayey loam 0.089 0.430 0.168 1.0 1.23
sandy clay 0.100 0.380 0.288 2.7 1.23
silty clay 0.070 0.360 0.00475 0.5 1.09
clay 0.068 0.380 0.0475 0.8 1.09
Models Manual).
A distinction can be made between Topsoil and Subsoil. In general, soils are considered
to be subsoils. The Type drop-down menu for the Staring series (Figure 6.21) contains
the following subsoils: Non-loamy sand (O1), Loamy sand (O2), Very loamy sand (O3),
Extremely loamy sand (O4), Coarse sand (O5), Boulder clay (O6), River loam (O7),
Sandy loam (O8), Silt loam (O9), Clayey loam (O10), Light clay (O11), Heavy clay (O12),
Very heavy clay (O13), Loam (O14), Heavy loam (O15), Oligotrophic peat (O16),
Eutrophic peat (O17) and Peaty layer (O18), and the following topsoils: Non-loamy sand
(B1), Loamy sand (B2), Very loamy sand (B3), Extremely loamy sand (B4), Coarse sand
(B5), Boulder clay (B6), Sandy loam (B7), Silt loam (B8), Clayey loam (B9), Light clay
(B10), Heavy clay (B11), Very heavy clay (B12), Loam (B13), Heavy loam (B14), Peaty
sand (B15), Sandy peat (B16), Peaty clay (B17) and Clayey peat (B18). The selected soil
type and grading (particle fractions) are different from the Hypres and the USDA data
sets. The parameters of the hydraulic model for the selected soil type are displayed in the
Soil tab at the right side of the Flow parameters tabsheet.
Hint: Only soil layers that are located not more than 1 m below the ground surface
are considered to be Upper soils.
User defined: The User defined option enables the user to define both saturated and
unsaturated properties manually. Please note that this option requires adequate
experience with unsaturated groundwater flow modelling. The hydraulic models available
are:
Van Genuchten This well-known and widely accepted model requires direct input
of the residual saturation Sres , the saturation at p = 0, Ssat and
the three fitting parameters gn , ga and gl (see Section 16.1 in the
Material Models Manual).
Spline The Spline function requires direct input of the capillary height ψ
(in unit of length), the relative permeability Kr (-), and the degree
of saturation Sr (-). Data for the Spline function can be entered
by clicking the Table tab. During the calculations, the flow
calculation kernel will use 'smooth' relationships based on a
spline function between the relative permeability and the
capillary height and also between the relative saturation and the
capillary height.
Saturated When the Saturated option is selected, no extra data input is
required. During the calculations, PLAXIS will continuously use
the saturated permeabilities for soil layers where a Saturated
data set was assigned.
permeability and degree of saturation. The negative sign indicates suction. Above the
level of ψunsat , the value of Kr and S remains constant. In this way a minimum degree of
saturation (Smin ) is guaranteed (Figure 6.23). It is used to limit the relative permeability
Kr and degree of saturation for high unsaturated zones.
By default a very large value is assigned to ψunsat (= 104 ). This value is only an indication
that the unsaturated zone is by default unlimited.
Change of permeability (ck ): This advanced feature is to account for the change of
permeability during a consolidation analysis. This can be applied by entering a proper
value for the ck parameter and the void ratio's. On entering a real value, the permeability
will change according to the formula:
k ∆e
log =
k0 ck
where ∆e is the change in void ratio, k is the permeability in the calculation and k0 is the
input value of the permeability in the data set (= kx and kz ). Note that a proper input of
the initial void ratio einit , in the General tabsheet is required. It is recommended to use a
changing permeability only in combination with the Hardening Soil model, Hardening Soil
model with small-strain stiffness, Soft Soil model or the Soft Soil Creep model. In that
case the ck -value is generally in the order of the compression index Cc . For all other
models the ck -value should be left to its default value of 1015 .
effective friction angle ϕ' and the dilatancy angle ψ '. In case of the User-defined soil
models, the tangent stiffness for primary oedometer loading Eoed ref
, the effective cohesion
c 'ref , the effective friction angle ϕ', the dilatancy angle ψ ' and the parameters UD-Power
and UD-Pref are required as interface parameters. For more information on the interface
parameters required for the User-defined soil models, see Section 14.3 in Material
Models Manual.
Interface strength
In case of the Linear Elastic model, the Mohr-Coulomb model, the Hardening Soil model,
the HS small model, the Soft Soil model, the Soft Soil Creep model, the Jointed Rock
model or the Hoek-Brown model, the interface strength is defined by the parameter Rinter .
The interface strength can be set using the following options:
Rigid: This option is used when the interface should not have a reduced strength with
respect to the strength in the surrounding soil. For example, extended interfaces around
corners of structural objects (Figure 5.32) are not intended for soil-structure interaction
and should not have reduced strength properties. The strength of these interfaces should
be assigned as Rigid (which corresponds to Rinter = 1.0). As a result, the interface
properties, including the dilatancy angle ψi , are the same as the soil properties in the
data set, except for Poisson's ratio νi (see further).
Manual: The value of Rinter can be entered manually if the interface strength is set to
Manual. In general, for real soil-structure interaction the interface is weaker and more
flexible than the surrounding soil, which means that the value of Rinter should be less than
1. Suitable values for Rinter for the case of the interaction between various types of soil
and structures in the soil can be found in the literature. In the absence of detailed
information it may be assumed that Rinter is of the order of 2/3. A value of Rinter greater
than 1 would not normally be used.
When the interface is elastic then both slipping (relative movement parallel to the
interface) and gapping or overlapping (i.e. relative displacements perpendicular to the
interface) could be expected to occur.
1 − νi
Eoed,i = 2 Gi
1 − 2 νi
2
Gi = Rinter Gsoil ≤ Gsoil
νi = 0.45
Hint: Note that a reduced value of Rinter not only reduces the interface strength,
but also the interface stiffness.
It is clear from these equations that, if the elastic parameters are set to low values, the
elastic displacements may be excessively large. If the values of the elastic parameters
are too large, however, this can result in numerical ill-conditioning of the stiffness matrix.
The key factor in the stiffness is the virtual thickness. This value is automatically chosen
such that an adequate stiffness is obtained.
Manual with residual strength: When the limit value of the interface strength as
defined by Rinter is reached, the interface strength may soften down to a reduced value
as defined by Rinter ,residual . Definition of the Rinter ,residual is possible when the Manual with
residual strength option is selected for the interface strength.
Interface strength (Rinter ): An elastic-plastic model is used to describe the behaviour of
interfaces for the modelling of soil-structure interaction. The Coulomb criterion is used to
distinguish between elastic behaviour, where small displacements can occur within the
interface, and plastic interface behaviour when permanent slip may occur. For the
interface to remain elastic the shear stress τ is given by:
|τ |= −σn tan ϕi + ci
where ϕi and ci are the friction angle and cohesion of the interface. The strength
properties of interfaces are linked to the strength properties of a soil layer. Each data set
has an associated strength reduction factor for interfaces Rinter . The interface properties
are calculated from the soil properties in the associated data set and the strength
reduction factor by applying the following rules:
ci = Rinter csoil
In addition to Coulomb's shear stress criterion, the tension cut-off criterion, as described
before (see Section 6.1.2), also applies to interfaces (if not deactivated):
Hint: Note that the same values of partial factors in Design approaches are applied
to both interface strength Rinter and residual interface strength Rinter ,residual .
» A reduced residual strength is not recommended to be used in Safety
calculations.
Consider gap closure: When the interface tensile strength is reached a gap may occur
between the structure and the soil. When the load is reversed, the contact between the
structure and the soil needs to be restored before a compressive stress can developed.
This is achieved by selecting the Consider gap closure option in the Interfaces tabsheet
of the Soil window. If the option is NOT selected, contact stresses will immediately
develop upon load reversal, which may not be realistic.
Interfaces using the Hoek-Brown model: When using the Hoek-Brown model as a
continuum model to describe the behaviour of a rock section in which interface elements
are used, equivalent interface strength properties ϕi , ci and σt,i are derived from this
model. The general shear strength criterion for interfaces as well as the tensile strength
criterion are still used in this case:
|τ |≤ −σn tan ϕi + ci
σn ≤ σt,i
Starting point for the calculation of the interface strength properties is the minor principal
effective stress σ '3 in the adjacent continuum element. At this value of confining stress
the tangent to the Hoek-Brown contour is calculated and expressed in terms of ϕ and c :
f'
sin ϕ =
2 + f'
2σ '3 sin ϕ
1 − sin ϕ
c= f+
2 cos ϕ 1 − sinϕ
where
−σ '3
f = σci mb +c a
σci
−σ '3
f ' = amb mb + s a−1
σci
and a, mb , s and ci are the Hoek-Brown model parameters in the corresponding material
data set. The interface friction angle ϕ'i and adhesion c 'i as well as the interface tensile
strength σt,i are now calculated using the interface strength reduction factor Rinter :
tan ϕi = Rinter c
ci = Rinter c
sσci
σt,i = Rinter σt = Rinter
mb
For more information about the Hoek-Brown model and an explanation of its parameters,
reference is made to Chapter 4 of the Material Models Manual.
Interfaces using the Modified Cam-Clay model: If the Modified Cam-Clay model is
selected in the Parameters tabsheet to describe the behaviour of the surrounding soil, the
following parameters are required to model the interface behaviour:
cref : Cohesion of the interface [kN/m2 ]
ϕi : Internal friction angle of the interface [◦ ]
ψi : Dilatancy angle of the interface [◦ ]
When the interface is elastic then both slipping (relative movement parallel to the
interface) and gapping or overlapping (i.e. relative displacements perpendicular to the
interface) could be expected to occur.
The magnitudes of these displacements are:
σn σn ti
Elastic gap displacement = =
Kn Eoed,i
τ τ ti
Elastic slip displacement = =
Ks Gi
where Gi is the shear modulus of the interface, Eoed,i is the one-dimensional
compression modulus of the interface and ti is the virtual thickness of the interface,
generated during the creation of interfaces in the geometry model (Section 5.6.7). KN is
the elastic interface normal stiffness and KS is the elastic interface shear stiffness. The
shear and compression moduli are related by the expressions:
3 (1 − νi ) σn
Eoed,i =
λ (1 + νi ) (1 + e0 )
3(1 − 2νi ) σn
Gi =
2(1 + νi ) λ(1 + e0 )
νi = 0.45
The K0 -values can be defined automatically by selecting the option Automatic in the K0
determination drop-down menu or manually by selecting the option Manual.
K0 -values
In general, two K0 -values can be specified, one for the x -direction and one for the
y -direction:
K0 = 1 − sin ϕ
For advanced models (Hardening Soil model, Hardening Soil model with small-strain
stiffness, Soft Soil model, Soft Soil Creep model, Modified Cam-Clay model) the default
value is based on the K0nc model parameter and is also influenced by the OCR -value and
POP -value in the following way:
νur
K0nc POP − POP
νur
K0,x = K0,y = K0nc OCR − (OCR − 1) + 0 νur
1 −
1 − νur σzz
The POP -value will result in a stress-dependent K0 -value within the layers resulting in
invisible K0 -values.
Be careful with very low or very high K0 -values, since these values might bring the initial
stress in a state of failure. For a cohesionless material it can easily be shown that to
avoid failure, the value of K0 is bounded by:
1 − sin ϕ 1 + sin ϕ
< K0 <
1 + sin ϕ 1 − sin ϕ
The checkbox can be used to set the K0,y value equal to the K0,x value.
σp
OCR=
σ '0zz POP
σ '0zz σp σ '0zz σp
a. Using OCR b. Using POP
Figure 6.26 Illustration of vertical pre-consolidation stress in relation to the in-situ vertical effective
stress
These two ways of specifying the vertical pre-consolidation stress are illustrated in Figure
6.26.
eq
The pre-consolidation stress σp is used to compute pp which determines the initial
eq
position of a cap-type yield surface in the advanced soil models. The calculation of pp is
based on the stress state:
where K0nc is the K0 -value associated with normally consolidated states of stress, which
is based on Jaky's formula, K0nc ≈ 1 − sin ϕ, or it is a direct input parameter for the
advanced soil models.
Hint: The Drainage type setting is only considered in a Plastic calculation, a Safety
analysis or a Dynamic analysis. When a Consolidation analysis or a Fully
coupled flow-deformation analysis is performed, the Drainage type is ignored
and the soil response is determined by the saturated permeability of the
material that is specified in the Flow parameters tabsheet of the material
datase.
Hint: The modelling of undrained soil behaviour is even more complicated than the
modelling of drained behaviour. Therefore, the user is advised to take the
utmost care with the modelling of undrained soil behaviour.
parameters will provide the correct undrained shear strength in the analysis. This is
because the effective stress path that is followed in an undrained analysis may not be the
same as in reality, due to the limitations of the applied soil model.
In order to enable a direct control on the shear strength, PLAXIS allows for an undrained
effective stress analysis with direct input of the undrained shear strength (Undrained (B)).
Hint: For Undrained (B) and Undrained (C) an increased shear strength with depth
can be modelled using the advanced parameter su,inc .
The SoilTest option is a quick and convenient procedure to simulate basic soil tests on
the basis of a single point algorithm, i.e. without the need to create a complete finite
element model. This option can be used to compare the behaviour as defined by the soil
model and the parameters of a soil data set with the results of laboratory test data
obtained from a site investigation. It also offers the possibility to optimise model
parameters such that a best fit is obtained between the model results and the lab test
data. The SoilTest facility works for any soil model, both standard soil models as well as
user-defined models.
The SoilTest option is available from the Material sets window if a soil data set is selected
(see Figure 6.27). Alternatively, the SoilTest option can be reached from the Soil dialog.
Once the SoilTest option has been selected, a separate window will open (Figure 6.28).
This window contains a menu, a toolbar and several smaller sections. The various items
are described in more detail below.
Main menu
The menus available in the menu bar are:
Figure 6.27 Material sets window showing the project and the global database
File To open, save and close a soil test data file (*.vlt).
Test To select the test that will be simulated. The options available
are Triaxial, Oedometer, CRS, DSS (Simple Shear) and General.
Results To select the configuration of diagrams to display.
Toolbar
The toolbar allows for loading, saving and running of soil test results and opening the
PLAXIS SoilTest - Settings window to set the configuration of the results. It also contains
the parameter optimisation feature (Section 6.3.7).
Material properties
The Material properties box displays the name, material model and parameters of the
currently selected data set. Transferring of material parameters to and from the material
database is possible. To copy the modified parameters to the material database:
Click the Copy material button in the Material properties box.
• In the program open the Material sets window and either select the corresponding
material set or click New.
In the Soil window click Paste material button. The parameters will be copied in the
material database. In the same way it is also possible to copy material from material
database to soil test.
Test area
The type of test and the testing conditions are defined in the test area. The test options
available are Triaxial, Oedometer, CRS, DSS and General. As one of these options is
selected by clicking the corresponding tab, the testing conditions can be defined in the
tabsheet. A more detailed description of the tests is given in the following sections.
Run
The Run button starts the currently selected test∗ . Once the calculation has finished, the
results will be shown in the Results window.
Test configurations
The Test configurations button can be used to add and manage different soil test
configurations. A test configuration contains information about the test type and the
values of test input parameters. To save a test configuration select the Save option in the
menu displayed as the Test configuration button is clicked. The Manage option can be
used to manage the test configurations available. When the Manage option is selected,
the Manage configuration window pops up. Note that the name of the window indicates
the test to which the configuration belongs (Figure 6.29).
The name and the location of the configuration file is indicated in the Filename and Path
respectively in the Manage configurations window.
∗ Although the soil test calculation kernel is a reduced version of the finite element calculation kernel, the
implementation of the soil models is identical.
Set as default
The Set as default button saves the current input parameters as the default parameters.
These will be initialised as such the next time the SoilTest window is opened.
Loaded tests
When previously saved tests of the current type have been opened from the File menu,
the Loaded tests window lists all these tests within each tabsheet. The results of all
loaded tests are shown together with the results of the current test. The Delete button
can be used to remove the selected test from the list of loaded tests. It does not remove
the soil test file (*.vlt) from disk.
Results
The results of the test are displayed in the predefined diagrams in the results area.
6.3.2 OEDOMETER
The Oedometer tabsheet contains facilities to define a one-dimensional compression
(oedometer) test. The following settings can be defined:
|vertical precons. stress|
The vertical pre-consolidation pressure to which the soil has
been subjected. If the soil is normally consolidated this value
should be set equal to the initial stress state, i.e. zero. From the
vertical pre-consolidation stress the program calculates the
isotropic pre-consolidation stress based on the K0nc loading path
(see Section 2.8 of the Material Models Manual). This option is
only available for the advanced soil models.
Mobilized relative shear strength
This option is only available for the Hardening Soil model and HS
small model to set the initial shear hardening contour. This value
must be between 0 (isotropic stress state) and 1 (failure state).
Phases Lists the different phases of the oedometer test. Each phase is
defined by a Duration (in units of time), a vertical Stress
increment (in units of stress) and a Number of steps. The initial
state is always assumed to be stress free. The given stress
increment will be reached at the end of the given duration in the
given number of steps. The input values can be changed by
clicking in the table. A negative stress increment implies
6.3.3 CRS
The CRS tabsheet contains facilities to define a constant rate-of-strain compression test.
The following settings can be defined:
|vertical precons. stress|
The vertical pre-consolidation pressure to which the soil has
been subjected. If the soil is normally consolidated this value
should be set equal to the initial stress state, i.e. zero. From the
vertical pre-consolidation stress the program calculates the
isotropic pre-consolidation stress based on the K0nc loading path
(see Section 2.8 of the Material Models Manual). This option is
only available for the advanced soil models.
Mobilized relative shear strength
This option is only available for the Hardening Soil model and HS
small model to set the initial shear hardening contour. This value
must be between 0 (isotropic stress state) and 1 (failure state).
Phases Lists the different phases of the CRS test. Each phase is defined
by a Duration (in units of time), a vertical Strain increment (in %)
and a Number of steps. The initial state is always assumed to be
stress free. The given strain increment will be reached at the end
of the given duration in the given number of steps. The input
values can be changed by clicking in the table. A negative strain
increment implies additional compression, whereas a positive
strain increment implies unloading or tension. If a period of zero
strain is desired, enter the desired duration with a zero strain
increment.
Add Adds a new phase to the end of the Phases list.
Insert Inserts a new phase before the currently selected phase.
Remove Removes the currently selected phase from the Phases list.
6.3.4 DSS
The DSS tabsheet contains facilities to define a direct simple-shear test. Before
specifying the test conditions, a selection can be made between different test options.
DSS - Options
Drained / undrained DSS test
In the latter case, undrained soil conditions and zero drainage
are assumed (similar as when the Drainage type has been set to
Undrained (A) or Undrained (B), see Section 6.2), irrespective of
the drainage type setting in the material data set.
Isotropically consolidated / K0 -consolidated test
In the latter case the K0 -value (ratio of lateral stress over axial
stress) can be specified to set the initial stress state.
DSS - Conditions
The following settings can be defined:
|vertical precons. stress|
The vertical pre-consolidation pressure to which the soil has
been subjected. If the soil is normally consolidated this value
should be set equal to the initial stress state or kept zero. From
the vertical pre-consolidation stress the program calculates the
isotropic pre-consolidation stress based on the K0nc loading path
(see Section 2.8 of the Material Models Manual). This option is
only available for the advanced soil models.
Mobilized relative shear strength
This option is only available for the Hardening Soil model and HS
small model to set the initial shear hardening contour. This value
must be between 0 (isotropic stress state) and 1 (failure state).
Initial stress |σyy | The absolute value of the initial vertical stress at which the
sample is consolidated, entered in units of stress. In the case of
an isotropically consolidated test, the initial lateral stress is equal
to the initial vertical stress. In the case of a K0 -consolidated test,
the initial lateral stress is equal to K0 σyy .
Time ∆t Time increment (only relevant for time-dependent models;
consolidation is not considered).
Number of steps The number of steps that will be used in the calculation.
Maximum shear strain |γxy |
The maximum value of shear strain (entered in %) that will be
reached in the last calculation step.
6.3.5 GENERAL
The General tabsheet contains facilities to define arbitrary stress and strain conditions.
The following settings can be defined:
Type of test The type of the test, whether Drained or Undrained can be
specified.
|vertical precons. stress|
The vertical pre-consolidation pressure to which the soil has
been subjected. If the soil is normally consolidated this value
should be set equal to the initial stress state or kept zero. From
the vertical pre-consolidation stress the program calculates the
isotropic pre-consolidation stress based on the K0nc loading path
6.3.6 RESULTS
The Results window shows several predefined typical diagrams to display the results of
the current test. Double-clicking one of the graphs opens the selected diagram in a larger
window (Figure 6.30). This window shows the selected diagram, the table of the data
points that are used to plot this diagram as well as the tangent and the secant values of
the plot. Note that the point to be taken into consideration for the calculation of the
tangent and the secant values can be determined by clicking on the plot. The diagram or
the data can be copied to the clipboard by selecting the corresponding option in the
drop-down menu displayed when the Copy button is clicked.
The diagram can be zoomed in or out using the mouse by first clicking and holding the
left mouse button in the diagram area and then moving the mouse to a second location
and releasing the mouse button. Moving the mouse from the left upper corner to the right
lower corner zooms the diagram to the selected area, whereas moving the mouse from
the right lower corner to the left upper corner resets the view. The zoom action can also
be undone using the Zoom out option on the toolbar.
The wheel button of the mouse can be used for panning: click and hold the mouse wheel
down and move the diagram to the desired position. When clicking the left mouse button
on a curve in the diagram, the corresponding secant and tangent line through the
selected point are indicated by dashed lines. This enables, for example, the
back-calculation of stiffness parameters from stress-strain diagrams. The corresponding
secant and tangent values are indicated below the table.
Select parameters
The Select parameters tab shows the parameters of the selected material data set that
could participate in the optimisation process. Click on the square in front of the
parameter(s) that need(s) to be optimised (Figure 6.32). The more parameters are
selected, the more time the optimisation process will take. For the selected parameters,
minimum and maximum values need to be specified. The optimisation algorithm will
search for optimum values within this range. If the optimised value turns out to be equal
to the minimum or maximum value, it might be that the best value lies outside the
specified range.
Note that parameters may influence only specific parts of a test. For example, when
considering a triaxial test, the initial part of the test curve is dominated by stiffness
parameters (such as E50 ), whereas the last part of the curve is dominated by strength
parameters (such as ϕ'). In order to obtain a best fit the optimisation should be
performed in separate runs; one for the stiffness parameter using the initial part of the
curve and one for the strength parameter using the last part of the curve, while fixing the
stiffness as the previously optimised value.
Select curves
The Select curves tab enables selection and uploading of real soil lab test data and
corresponding test conditions. Alternatively, synthetic test data may be used in the form
of other PLAXIS soil test results. In this way it is possible to optimise, for example,
parameters of the Mohr-Coulomb model against simulated tests using the Hardening Soil
model.
Initially, the window shows a tree with the five standard test types (Triaxial, Oedometer,
CRS, DSS and General). For each test type, different test conditions can be defined,
which can be taken into account in the optimisation process. By default, the Current
model test is available as test conditions for each test type. The Current model test
contains the test conditions as previously defined for that test (Figure 6.33).
New test conditions can be defined by selecting the New test configuration option
from the tool bar. This will introduce Custom # under the selected test, for which the
test conditions can be defined in the right-hand panel (Figure 6.34).
In both cases (Current model test and Custom test) corresponding test data need
to be selected and uploaded using the Import curve option. Another possibility is to
upload test conditions together with the test data in case it is stored in the format of a
PLAXIS soil test project (<test>.vlt).
Hence, there are different ways to define test conditions and to select the external test
data. The possibilities are summarized below:
• If the test data corresponds to one of the Current model test conditions, the
corresponding line should be selected in the tree and the Import curve option should
be used to upload the test data (Figure 6.35). The test data are assumed to be
Figure 6.33 Selection of the test curves in the Select curves tabsheet
stored in a text file (<data>.txt) and should contain two columns, separated by a
Space, Tab, Comma, Colon (:), Semicolon (;) or arbitrary character. The separator is
to be indicated at the top of the Import test data window. The meaning of the values
in each column has to be selected from the drop down list below the column. Here,
a selection can be made amongst various stress and strain quantities. Moreover, the
basic units of the test data quantities need to be selected from the drop down lists in
the Units group. By pressing OK the data is read and visualised in a diagram, and
the curve is listed in the tree under the Current model (test) conditions.
• If the test data corresponds to other than one of the current model test conditions,
first new Custom test conditions need to be defined. Select the appropriate test type
and click the New test configuration button. The test conditions of the data to be
uploaded can be defined in the right-hand panel. Subsequently, the Import curve
option should be used to upload the test data. The test data are assumed to be
stored in a text file (<data>.txt) and should contain two columns (see explanation
before). The meaning of the values in each column has to be selected from the drop
down list below the column. Moreover, the basic units of the test quantities need to
be selected from the drop down lists in the Units group. By pressing OK the data is
read and a visualised in a diagram, and the curve is listed in the tree under the
Custom (test) conditions (Figure 6.36).
• If the test data together with the test conditions are stored in the format of a PLAXIS
soil test project (<test>.vlt), the Open file option should be used. After selection of a
valid PLAXIS soil test project, the test conditions are listed under the corresponding
test type in the tree, and the available test data curves are listed under the test
conditions (Figure 6.37). This option should typically be used to fit current model
parameters to synthetic data previously produced in the PLAXIS soil test facility and
stored in <test>.vlt format.
Hint: When a line representing test conditions is selected in the tree, the
corresponding test conditions are shown on the panel at the right-hand side.
» When a line representing test data is selected in the tree, the corresponding
curve is visualised in the diagram, and a table of corresponding data points is
shown at the right-hand side of the diagram.
» A sub-set of test data to be used in the optimisation process can be selected
in the table at the right-hand side by clicking on the corresponding cells,
using the standard multi-select convention (using <Shift> for ranges and
<Ctrl> for individual values). The selected values are indicated as 'thick' lines
in the curve whereas non-selected values are indicated as 'thin' lines.
» A line in the tree (either test conditions or test data) can be removed by
selecting that line and clicking the red cross in the toolbar.
All test data to be used in the optimisation process need to be selected in the tree by
clicking the square in front of the corresponding line (if not already selected). The
corresponding test conditions are automatically selected.
Multiple phases
In the case of an Oedometer, CRS or General test, the SoilTest facility allows for multiple
phases. However, the parameter optimisation facility can only deal with one phase at a
time. Therefore, after importing the test data, the desired calculation phase needs to be
selected from the drop down list above the test data curve, together with the
corresponding part of the test data in the column at the right-hand side. In this way it is
possible, for example, to optimise a primary loading stiffness against the first (loading)
phase in an oedometer test and the unloading stiffness against the second (unloading)
phase. Note that the test data needs to be for each phase.
Settings
The Settings tab enables the accuracy selection of the optimisation process (Figure
6.38). Three levels of search intensity are available: Coarse and quick, Moderate,
Thorough. In addition, the relative tolerance of the search algorithm can be selected. The
default value is 1E-3. Note that a more rigorous optimisation may give more accurate
results, but also requires more calculation time. The calculation time also depends on the
number of parameters to be optimised, as selected in the first tab.
Resulting parameters
The Resulting parameters tab shows the optimum values of the parameters used to
obtain the best fit to the selected test data in addition to the minimum and maximum
values and the reference values in the material data set (Figure 6.39). If the optimum
value is equal to the minimum or maximum value, it might be that the best value lies
outside the specified range. Finally, the table shows the sensitivity of the selected
parameters. A sensitivity of 100% means that the parameter has a high influence on the
simulated test results, whereas a low sensitivity means that the parameter has a low
influence on the simulated test results. Note that a low sensitivity also means that the test
may not be suitable to optimise that parameter and, as a result, the suggested optimum
value may not be accurate. Therefore it is better to do separate optimisations for different
parameters based on relevant sections of test data curves rather than one optimisation
with multiple parameters based on the full data curves.
Resulting charts
The Resulting charts tab shows the results of the selected tests (Figure 6.40).
For each test, three curves are visible:
Optimisation target This curve represents the uploaded test data.
Optimisation results This curve represents the simulated test with optimised
parameters.
Reference simulation This curve represents the simulated test with original
parameters. It has no meaning in the optimisation process, but
just shows how good or bad the existing material data set would
fit the uploaded test data for the selected test conditions without
optimisation.
Limitations
The Parameter optimisation facility should be used with care. Note that parameters
optimised for soil lab tests may not be the best parameters for the practical application as
considered in the finite element model. This is because the application may involve
stress levels, stress paths and strain levels which might be significantly different from the
ones that occur in the soil lab tests.
Furthermore, the parameter optimisation facility has the following limitations:
• It is not possible to automatically optimise test data curves that consist of multiple
phases (for example loading and unloading phases). Such curves need to be
uploaded for each phase and the corresponding part of the curves (phase) needs to
In addition to material data sets for soil and interfaces, the material properties and model
parameters for plates are also entered in separate material data sets. A data set for
plates generally represents a certain type of plate material or plate profile, and can be
assigned to the corresponding (group of) plate elements in the geometry model.
2 3
N1 E1 ε1 N2 E2 ε2
Figure 6.42 Definition of positive normal forces (N), shear forces (Q) and bending moments (M) for
a plate based on local system of axes
6.4.2 PROPERTIES
The properties required for plates can be grouped into general properties and stiffness
properties.
General properties
A plate has two general properties:
d: The (equivalent) thickness (in the unit of length) is the material cross section
area of the plate across its major axial direction per 1 m width. For massive
plates without a particular profile this is just the plate thickness, but for plates
that have a certain profile (such as sheet-pile walls or sandwich plates), the
thickness is relatively small and should be properly calculated from the above
definition.
γ: The unit weight is the unit weight of the material from which the plate is
composed. The product γ · d determines the distributed weight of the plate.
Hint: When specifying a unit weight, please consider the fact that the element itself
does not occupy any volume and overlaps with the soil elements. Hence, it
might be considered to subtract the unit soil weight from the real unit weight
of the plate material in order to compensate for the overlap. For partially
overlapping plates the reduction of the unit weight should be proportional.
Stiffness properties
Plate stiffnesses will be linear. PLAXIS 3D allows for orthotropic as well as anisotropic
material behaviour in plates, which is defined by the following parameters:
E1 : Young's modulus in first axial direction.
E2 : Young's modulus in second axial direction.
G12 : In-plane shear modulus.
G13 : Out-of-plane shear modulus related to shear deformation over first direction.
G23 : Out-of-plane shear modulus related to shear deformation over second
direction.
ν12 : Poisson's ratio.
These parameters appear in the following (approximate) relationships for structural
forces:
N1 E1 d ν12 E2 d ε1
=
N2 ν12 E2 d E2 d ε2
Q12 G12 d 0 0 γ12
∗
Q13 = 0 kG13 d 0 γ
13
∗
Q23 0 0 kG23 d γ23
E d3 ν E d3
1 12 2 0
M11 12 12 κ11
ν12 E2 d 3 E2 d 3
M22 = 0 κ
12 12 22
G12 d 3
M12 0 0 κ12
12
Where k represents the modification factor to the shear modulus and is taken equal to 5/6 .
Figure 6.42 visualises the plate's local system of axes and the major quantities. The local
system of axes in a plate element is such that the first and the second local axis lie in the
plane of the plate whereas the third axis is perpendicular to the plane of the plate.
If the Isotropic option is checked the input is limited to E1 and ν12 , where as E2 = E1 and
G12 = G13 = G23 = E/2(1 + ν12 ).
More information about the behaviour and structural forces in plates can be found in
Chapter 15.6 of the Material Models Manual.
In addition to material data sets for soil and interfaces, the material properties and model
parameters for geogrids are also entered in separate material data sets. Geogrids are
flexible elastic elements that represent a grid or sheet of fabric. Geogrids cannot sustain
compressive forces. A data set for geogrids generally represents a certain type of
geogrid material, and can be assigned to the corresponding (group of) geogrid elements
in the geometry model.
Material type:
There are two available options, describing the material type of a plate. These
options are Elastic and Elastoplastic. The availability of the parameters
defined in the Properties box depends on the selected material type.
6.5.2 PROPERTIES
The properties required for geogrids can be grouped into stiffness properties and
strength properties in case of elastoplastic behaviour.
Isotropic
Different stiffnesses in-plane and out-of-plane may be considered. The latter is most
relevant for axisymmetric models when modelling geogrids with an anisotropic pattern. If
this is not the case, the Isotropic option may be selected to ensure that both stiffness are
equal.
Stiffness properties
For elastic behaviour, the axial stiffness EA should be specified. PLAXIS 3D allows for
orthotropic as well as anisotropic material behaviour in geogrids, which is defined by the
following parameters:
EA1 : The normal elastic stiffness in 1-direction (in plane).
EA2 : The normal elastic stiffness in 2-direction (out of plane, anisotropic behaviour).
GA: In-plane shear stiffness (anisotropic behaviour).
The axial stiffness EA is usually provided by the geogrid manufacturer and can be
determined from diagrams in which the elongation of the geogrid is plotted against the
applied force in a longitudinal direction. The axial stiffness is the ratio of the axial force
per unit width and the axial strain (∆l/l where ∆l is the elongation and l is the length):
F
EA =
∆l/l
If the Isotropic option is checked the input is limited to EA1 where as EA1 =EA2 and GA
= EA1 /2.
In addition to material data sets for soil and interfaces, the material properties and model
parameters for beams are also entered in separate material data sets. A data set for
beams generally represents a certain type of beam material or beam profile, and can be
assigned to the corresponding (group of) beam elements in the geometry model.
6.6.2 PROPERTIES
The properties required for beams can be grouped into general properties and stiffness
properties.
General properties
A beam has two general properties:
A: The cross section area is the actual area (in the unit of length squared)
perpendicular to the axial beam direction where beam material is present. For
beams that have a certain profile (such as steel beams), the cross section
area can be found in tables that are provided by steel factories.
γ: The unit weight (in the unit of force per unit of volume) is the unit weight of the
material from which the beam is composed. The product γ · A determines the
distributed beam weight.
Hint: When specifying a unit weight, please consider the fact that the element itself
does not occupy any volume and overlaps with the soil elements. Hence, it
might be considered to subtract the unit soil weight from the real unit weight
of the beam material in order to compensate for the overlap. For partially
overlapping beams the reduction of the unit weight should be proportional.
» Please note that when reducing the unit weight, the axial forces in the beam
may not be realistic.
Stiffness properties
Stiffnesses can only be linear. Beam stiffnesses involve the following properties:
E: Young's modulus.
I3 : Moment of inertia against bending around the third axis.
1
1 1 1
3 3 3 3
2 2 2 2
I3 M3 κ3 I2 M2 κ2 EN
Figure 6.45 Definition of moment of inertia (I ), positive bending moment (M ), positive curvature (κ)
and stiffness (E ) for a horizontal beam based on local system of axes
More information about the behaviour and structural forces in beams can be found in
Chapter 15.2 of the Material Models Manual.
Properties and model parameters for embedded piles are entered in separate material
data sets. A data set for embedded piles generally represents a certain type of pile,
including the pile material and geometric properties, as well as the interaction properties
with the surrounding soil (pile bearing capacity).
Note that the embedded pile material data set does not contain so-called 'p-y curves', nor
equivalent spring constants. In fact, the stiffness response of an embedded pile
subjected to loading is the result of the specified pile length, equivalent radius, spacing,
stiffness and bearing capacity as well as the stiffness of the surrounding soil.
Hint: In contrast to what is common in the Finite Element Method, the bearing
capacity of an embedded pile is considered to be an input parameter rather
than the result of the finite element calculation. The user should realise the
importance of this input parameter. Preferably, the input value of this
parameter should be based on representative pile load test data. Moreover, it
is advised to perform a calibration in which the behaviour of the embedded
pile is compared with the behaviour as measured from the pile load test.
Since embedded piles are used in a row, the group action must be taken into
account when defining the pile bearing capacity.
Geometric properties
An embedded pile requires several geometric parameters used to calculate additional
properties:
Pile type: Either a Predefined or a User defined type can be selected.
Dynamic properties
For dynamic behaviour, two additional parameters can be specified as material
properties:
Rayleigh α:
Rayleigh damping parameter determining the influence of mass in the
damping of the system.
Rayleigh β :
Rayleigh damping parameter determining the influence of the stiffness in the
damping of the system.
For more information on Rayleigh damping, see Page 115.
From the pile geometric properties an equivalent radius for the elastic zone, Req , is
determined:
np o
(2Iavg /A) where Iavg = (I2 + I3 )/2
p
Req = max (A/π),
The definition of various quantities according to the pile's local system of axes are
visualised in Figure 6.47. More information about the behaviour and forces in embedded
piles can be found in Section 15.8 of the Material Models Manual.
2 2 2 2
3 3 3 3
1 1 1 1
I3 M3 κ3 I2 M2 κ2 ENε
Figure 6.47 Definition of moment of inertia (I ), positive bending moment (M ), positive curvature (κ)
and stiffness (E ) for a vertical user-defined pile type (beam profile) based on the local
system of axes
|ts |= Tmax
The input for the shaft resistance is defined by means of the skin resistance at the pile
top, Ttop,max (in force per unit pile length) and the skin resistance at the pile bottom,
Tbot,max (in force per unit pile length). This way of defining the pile skin resistance is
mostly applicable to piles in a homogeneous soil layer. Using this approach the total pile
bearing capacity, Npile , is given by:
1
Npile = Fmax + Lpile Ttop,max + Tbot,max
2
where Lpile is the pile length.
Hint: Note that the length of the embedded pile and the magnitude of the skin
resistance increments are inversely proportional.
In addition to the shaft resistance, the embedded pile has extra bearing capacity at the
base. The base resistance Fmax can be entered directly (in the unit of force) in the
embedded pile material data set window.
Hint: The base resistance is only mobilized when the pile body moves in the
direction of the base (example: with a load on top).
The pile bearing capacities are automatically divided by the pile spacing in order to obtain
the equivalent bearing capacity per unit of width in the out-of-plane direction.
Linear
The first way (Linear) is the easiest way to enter the skin resistance profile. The input is
defined by means of the skin resistance at the pile top, Ttop,max (in force per unit pile
length) and the skin resistance at the pile bottom, Tbot,max (in force per unit pile length).
This way of defining the pile skin resistance is mostly applicable to piles in a
homogeneous soil layer. Using this approach the total pile bearing capacity, Npile , is given
by:
1
Npile = Fmax + Lpile Ttop,max + Tbot,max
2
where Lpile is the pile length.
Hint: Note that the length of the embedded pile and the magnitude of the skin
resistance increments are inversely proportional.
Multi-linear
The second way (Multi-linear) can be used to take into account inhomogeneous or
multiple soil layers with different properties and, as a result, different resistances. The
skin resistance, Tmax , is defined in a table at different positions along the pile, L, where L
is measured from the pile top (L=0) to the bottom of the pile (L = Lpile ). Using this
approach the total pile bearing capacity, Npile , is given by:
n−1
X 1
Npile = Fmax + (Li+1 − Li )(Ti + Ti+1 )
i=1
2
Layer dependent
The third way (Layer dependent) can be used to relate the local skin resistance to the
strength properties (cohesion c and friction angle ϕ) of the soil layer in which the pile is
located, and the interface strength reduction factor, Rinter , as defined in the material data
set of the corresponding soil layers (Section 6.1.4). In this respect the special interface in
the embedded pile behaves similar as an interface along a wall, except that it is a line
interface rather than a sheet. Using this approach the pile bearing capacity is based on
the stress state in the soil, and thus unknown at the start of a calculation. To avoid that
the skin resistance could increase to undesired high values, an overall maximum
resistance (constant value along the pile in force per unit pile length) can be specified in
the embedded pile material data set, which acts as an overall cut-off value.
Hint: The pile-soil interaction parameters in the embedded pile material data set
involve only the pile bearing capacity (skin resistance and base resistance).
Note that the material data set does NOT include the stiffness response of
the pile in the soil (or p − y curve). The stiffness response is the result of the
pile length, equivalent radius, stiffness and bearing capacity as well as the
stiffness of the soil layers in which the pile is located.
In order to ensure that a realistic pile bearing capacity as specified can actually be
reached, a zone in the soil volume elements surrounding the beam is identified where
any kind of soil plasticity is excluded (elastic zone). The size of this zone is determined by
the embedded pile's diameter or equivalent radius Req (Section 6.7). The elastic zone
makes the embedded pile almost behave like a volume pile. However, installation effects
of piles are not taken into account and the pile-soil interaction is modelled at the centre
rather than at the circumference.
In addition to displacement differences and shear forces in axial direction along the pile,
the pile can undergo transverse forces, t⊥ , due to lateral displacements. These
transverse forces are not limited in the special interface element that connects the pile
with the soil, but, in general, they are limited due to failure conditions in the surrounding
soil itself outside the elastic zone. However, embedded piles are not meant to be used as
laterally loaded piles and will therefore not show accurate failure loads when subjected to
transverse forces.
More details about the way the shear and transverse forces are calculated on the basis of
displacement differences between the embedded beam element and the surrounding soil
element are described in the Material Models Manual.
In addition to material data sets for soil and interfaces, the material properties and model
parameters for anchors are also entered in separate material data sets. A material data
set for anchors may contain the properties of node-to-node anchors as well as fixed-end
anchors. In both cases the anchor is just a spring element. A data set for anchors
generally represents a certain type of anchor material, and can be assigned to the
corresponding (group of) anchor elements in the geometry model.
Identification:
A user may specify any identification title for a data set. It is advisable to use a
meaningful name since the data set will appear in the database tree view by its
identification.
Comments:
A user may write down comments related to the material data set.
Colour: Colour can be used as a distinction tool in the model.
Material type:
There are three available options, describing the material type of an anchor.
These options are Elastic and Elastoplastic. The availability of the parameters
defined in the Properties box depends on the selected material type.
6.8.2 PROPERTIES
The properties required for anchors can be grouped into stiffness properties and strength
properties in case of elastoplastic behaviour.
Stiffness properties
An anchor requires only one stiffness parameter:
EA: Axial stiffness, entered per anchor in the unit of force
Figure 6.49 The force-displacement diagram displaying the elastoplastic behaviour of anchors
In the same way as the stiffness, the maximum anchor forces are divided by the
out-of-plane spacing in order to obtain the proper maximum force in a plane strain
analysis.
The Elastoplastic with residual strength option can be used to model anchor failure or
softening behaviour (e.g. buckling of struts). When this option is selected two residual
anchor forces can be specified:
Fresidual,tens :
Residual tension force
Fresidual,comp :
Residual compression force
The Force-displacement diagram displaying the elastoplastic behaviour with residual
strength of the anchors is given in Figure 6.50.
If, during a calculation, the maximum anchor force is reached, the maximum force will
immediately reduce to the residual force. From that point on the anchor force will not
exceed the residual force anymore. Even if the anchor force would intermediately reduce
to lower values, the defined residual force will be its maximum limit.
Note that if the anchor has failed (in tension, compression or both) the residual force will
be valid in the following calculation phases where the anchor is active. If the anchor is
deactivated in a phase and reactivated in the next phase, the maximum anchor force will
be restored, assuming that the anchor is a completely a new one.
Figure 6.50 The force-displacement diagram displaying the elastoplastic behaviour with residual
strength of the anchors
window. The prestress force is not considered to be a material property and is therefore
not included in an anchor data set.
After creating material data sets, the data sets must be assigned to the corresponding
geometry components (soil layers and structures). This can be done in different ways,
which are explained below. The methods described below are primarily meant to assign
properties to the initial geometry. For details on the change of properties during
calculations in the framework of staged construction (Section 7.9.6).
Soil layers
Regarding soil data, material data sets can be assigned to individual soil layers in the
boreholes. Therefore a borehole should be double clicked to open the corresponding
Modify soil layers window. In the Modify soil layers window the Materials button at the
lower right hand side of the window should be clicked to open the material database.
To assign a data set to a particular soil layer, select the desired data set from the material
database tree view (click on the data set and hold the left hand mouse button down),
drag it to the soil column in the borehole window (hold the mouse button down while
moving) and drop it on the desired layer (release the mouse button). The layer should
now show the corresponding material data set colour. The drag and drop procedure
should be repeated until all layers have their appropriate data set. Note that material sets
cannot be dragged directly from the global database tree view and must be copied to the
project database first.
When multiple boreholes are used it should be noted that assigning a data set to a layer
in one particular borehole will also influence the other boreholes, since all layers appear
in all boreholes, except for layers with a zero thickness.
Structures
Regarding structures (fixed-end anchors, beams, node-to-node anchors, embedded
piles, plates, geogrids and interfaces), there are three different methods of assigning
material data sets.
• The first method is based on an open Material sets window, showing the created
material sets in the project database tree view. The desired material set can be
dragged (select it and keep the left mouse button down) to the draw area and
dropped on the desired component. It can be seen from the shape of the cursor
whether or not it is valid to drop the material set.
• The second method is to select the desired structure in the draw area or model
explorer and select the desired material set from the Set material of the right hand
mouse button menu.
• The third method is to select the desired structure in the draw area or model explorer
and use the Material combo box in the explorer to define the material data set.
When the geometry modelling process is complete, one can proceed with the
calculations. This consists of generation of mesh and definition of the construction
stages.
In the engineering practice, a project is divided into project phases. Similarly, a
calculation process in PLAXIS is also divided into calculation phases. Examples of
calculation phases are the activation of a particular loading at a certain time, the
simulation of a construction stage, the introduction of a consolidation period, the
calculation of a safety factor, etc. Each calculation phase is generally divided into a
number of calculation steps. This is necessary because the non-linear behaviour of the
soil requires loadings to be applied in small proportions (called load steps). In most
cases, however, it is sufficient to specify the situation that has to be reached at the end of
a calculation phase. Robust and automatic procedures in PLAXIS will take care of the
sub-division into appropriate load steps.
The construction stages can be defined in the Water levels and Staged construction
modes.The first calculation phase (Initial phase) is always a calculation of the initial
stress field for the initial geometry configuration by means of Gravity loading or K0
procedure. Alternatively, it may be indicated that the calculations only involve
groundwaterflow. After this initial phase, subsequent calculation phases may be defined
by the user. In each phase, the type of calculation must be selected.
For deformation calculations distinction is made between Plastic, Consolidation, Fully
coupled flow-deformation or Safety calculation. The different types of calculations are
explained in Section 7.3.
To perform finite element calculations, the geometry has to be divided into elements. A
composition of finite elements is called a finite element mesh. The mesh is created in the
Mesh mode. The mesh should be sufficiently fine to obtain accurate numerical results.
On the other hand, very fine meshes should be avoided since this will lead to excessive
calculation times. The PLAXIS 3D program allows for a fully automatic generation of
finite element meshes. The mesh generation process takes into account the soil
stratigraphy as well as all structural objects, loads and boundary conditions.
Elements
The basic soil elements of the 3D finite element mesh are the 10-node tetrahedral
elements (Figure 7.1).
In addition to the soil elements, special types of elements are used to model structural
behaviour. For beams, 3-node line elements are used, which are compatible with the
3-node edges of a soil element. In addition, 6-node plate and geogrid elements are used
to simulate the behaviour of plates and geogrids respectively. Moreover, 12-node
interface elements are used to simulate soil-structure interaction behaviour. The element
formulations are given in the Scientific Manual.
η
4
4
10
9 8
1 3
2 51 7 3 ξ
6
ζ 2
Figure 7.1 3D soil elements (10-node tetrahedrons)
To generate the mesh, click the Generate mesh button in the side toolbar
of the Mesh mode or select the corresponding option in the Mesh menu. The Mesh
options window pops up where the general mesh properties can be defined (Figure 7.2).
The mesh is generated by clicking on the OK button in the Mesh options window.
mesh locally.
To locally refine the mesh either click the Finer mesh button and select the
geometry entity (volume, surface, line or point) where a finer mesh is required or
right click the geometry entity and select the Finer mesh option from the appearing menu.
To locally coarsen the mesh click the Coarser mesh button and select the geometry
entity (volume, surface, line or point) where a coarser mesh is required or right click
the geometry entity and select the Finer mesh option from the appearing menu.
To reset the local refinement in a geometry entity right click it either in the draw area
or in Selection explorer and select the corresponding option from the appearing
menu.
Hint: In the Mesh mode, the colours indicating the material data sets are not
displayed. The whole model is shown in gray. However refined and
coarsened objects are displayed in green and yellow respectively.
» The more the object is refined, the lighter is the shadow of green it is
displayed.
» The more the object is coarsened, the lighter is the shadow of yellow it is
displayed.
Hint: Although interface elements have a zero thickness, the interfaces in the
mesh are drawn with a certain thickness to show the connections between
soil elements, structural elements and interface elements. The scale factor
(Section 8.3.7) can be used to reduce the graphical thickness of the
interfaces.
Finite element calculations can be divided into several sequential calculation phases.
Each calculation phase corresponds to a particular loading or construction stage.
The construction stages can be defined in the Staged construction mode. The calculation
phases are listed in the Phases explorer.
Hint: Note that for a new project the initial phase is automatically added. It cannot
be deleted. More information on the initial phase is given in Section 7.3.1.
Toolbar
The buttons in the toolbar enable introducing new phases, removing phases and
accessing the Phases window where the settings of phases are defined.
To introduce a new calculation phase (a child phase), select a reference phase
(parent phase) in the list and click the Add phase button. A new phase is added
Figure 7.4 Extended view of the right mouse menu of Phases explorer
Phase identification
The ID of the phases is displayed in the Phases explorer. The ID of the phase consist of
the caption and the name (part inside the square brackets). The name of the phases is
determined consecutively by the program and it can not be modified by the user. The
user may redefine the caption part of the ID of the phase in the Phases window (Section
7.7.1).
Special cases
In some special cases, the order of calculation phases is not straightforward. Examples
of some cases are:
• The Initial phase may be selected as reference if different loadings or loading
sequences are to be considered separately for the same project.
• For a certain situation, a load is increased until failure to determine the safety
margin. When continuing the construction process, the next phase should start from
the previous construction stage rather than from the failure situation.
• A third example where the phase ordering is not straightforward is in calculations
where safety analysis for intermediate construction stages is considered. The
calculation type in this case is Safety. In general, such a phase results in a state of
failure. When continuing the construction process, the next stage should start from
the previous phase rather than from the results of the safety analysis. Alternatively,
safety analyses for the various construction stages can be performed at the end of
the calculation process. In that case, the reference phase selected in the Start from
phase drop-down menu should refer to the corresponding construction stage.
The Spreadsheet view displays all the calculation information in a spreedsheet form
(Figure 7.6) to facilitate a comparison of parameters for different phases.
In the Phases window, users need to select at least the Calculation type and the Loading
type for each new phase. PLAXIS provides convenient default values for most calculation
control parameters, but the user can change these values. A description of the
calculation types and control parameters is given in the next section.
The first step in a PLAXIS analysis is defining a calculation type of a phase in the
Calculation type drop-down menu in the Phases window. The options available are K0
procedure and Gravity loading for the initial phase, and Plastic, Consolidation, Safety and
Dynamic for other phases.
Hint: As a rule, one should use the K0 procedure only in cases with a horizontal
surface and with all soil layers and phreatic levels parallel to the surface. For
all other cases, Gravity loading should be used.
K0 procedure
K0 procedure is a special calculation method available in PLAXIS to define the
initial stresses for the model, taking into account the loading history of the soil. The
parameters required in the initial stresses development procedures are defined in the
Initial tabsheet of material data sets for soil and interfaces (Section 6.1.5).
Two K0 values can be specified, one for the x -direction and one for the y -direction.
In practice, the value of K0 for a normally consolidated soil is often assumed to be related
to the friction angle by Jaky's empirical expression:
K0 = 1 − sinϕ
Using very low or very high K0 -values in the K0 procedure may lead to stresses that
violate the Mohr-Coulomb failure condition. In this case PLAXIS automatically reduces
the lateral stresses such that the failure condition is obeyed. Hence, these stress points
are in a plastic state and are thus indicated as plastic points. Although the corrected
stress state obeys the failure condition, it may result in a stress field which is not in
equilibrium. It is generally preferable to generate an initial stress field that does not
contain Mohr-Coulomb plastic points.
Hint: The plot of plastic points may be viewed after the presentation of the initial
effective stresses in the Output program by selecting the Plastic points option
from the Stresses menu (see Section 9.3.7).
For a cohesionless material it can easily be shown that to avoid Mohr-Coulomb plasticity,
the value of K0 is bounded by:
1 − sin ϕ 1 + sin ϕ
< K0 <
1 + sin ϕ 1 − sin ϕ
When the K0 procedure is adopted, PLAXIS will generate vertical stresses that are in
equilibrium with the self-weight of the soil. Horizontal stresses, however, are calculated
from the specified value of K0 . Even if the value of K0 is chosen such that plasticity does
not occur, the K0 procedure does not ensure that the complete stress field is in
equilibrium. Full equilibrium is only obtained for a horizontal soil surface with any soil
layers parallel to this surface and a horizontal phreatic level. If the stress field requires
only small equilibrium corrections, then these may be carried out using the calculation
procedures described below. If the stresses are substantially out of equilibrium, then the
K0 procedure should be abandoned in favor of the Gravity loading procedure.
At the end of the K0 procedure, the full soil is weight activated. The soil weight can not be
changed in any other calculation phase.
Gravity loading
Gravity loading is a type of Plastic calculation (Section 7.3.2), in which initial
stresses are generated based on the volumetric weight of the soil. If Gravity loading
is adopted, then the initial stresses are set up by applying the soil self-weight in the first
calculation phase. In this case, when using an elastic perfectly-plastic soil model such as
the Mohr-Coulomb model, the ratio of horizontal effective stress over vertical effective
stress, K0 , depends strongly on the assumed values of Poisson's ratio. It is important to
choose values of Poisson's ratio that give realistic values of K0 . If necessary, separate
material data sets may be used with Poisson's ratio adjusted to provide the proper
K0 -value during gravity loading. These sets may be changed by other material sets in
subsequent calculations (Section 7.9.6). For one-dimensional compression an elastic
computation will give:
ν
K0 =
(1 − ν)
If a value of K0 of 0.5 is required, for example, then it is necessary to specify a value of
Poisson's ratio of 0.333. As Poisson's ratio must be lower than 0.5, it is not
straightforward to generate K0 values larger than 1 using Gravity loading. If K0 values
larger than 1 are desired, it is necessary to simulate the loading history and use different
Poisson's ratio for loading and unloading or use the K0 procedure.
When advanced soil models are used, the resulting K0 -value after gravity loading
corresponds to the K0nc in the material data set.
Hint: To make sure that Gravity loading results in initial effective stresses in
situations where undrained materials are used, the parameter Ignore
undrained behaviour should be selected.
» Once the initial stresses have been set up using Gravity loading, the
displacements should be reset to zero at the start of the next calculation
phase. This removes the effect of the initial stress generation procedure on
the displacements developed during subsequent calculations, whereas the
stresses remain.
In some cases plastic points will be generated during the Gravity loading procedure. For
cohesionless soils in one-dimensional compression, for example, plastic Mohr-Coulomb
points will be generated unless the following inequality is satisfied:
1 − sin ϕ ν
< <1
1 + sin ϕ 1−ν
If there are a small number of plastic points, it is advisable to perform a plastic nil-step.
When using the Hardening Soil model and defining a normally consolidated initial stress
state (OCR = 1.0 and POP = 0.0), the plot of plastic points shows many hardening points.
Users need not be concerned about these plastic points as they just indicate a normally
consolidated stress state.
Plastic nil-step
If the K0 procedure generates an initial stress field that is not in equilibrium or where
Mohr-Coulomb plastic points occur, then a plastic nil-step should be adopted. A plastic
nil-step is a plastic calculation step in which no additional load is applied (Section 7.3.7).
After this step has been completed, the stress field will be in equilibrium and all stresses
will obey the failure condition.
If the original K0 procedure generates a stress field that is far from equilibrium, then the
plastic nil-step may fail to converge. This happens, for example, when the K0 procedure
is applied to problems with very steep slopes. For these problems, the Gravity loading
procedure should be adopted.
It is important to ensure that displacements calculated during a plastic nil-step (if it is
applied immediately after generating the initial stresses) do not affect later calculations.
This is achieved by selecting the Reset displacements to zero parameter in the
subsequent calculation phase (Section 7.9).
give a reasonably accurate prediction of the final situation, although the precise loading
history is not followed and the process of consolidation is not dealt with explicitly.
An elastic-plastic deformation analysis where undrained behaviour (Undrained (A) or
Undrained (B)) is temporarily ignored can be defined by checking the Ignore undr.
behaviour (A, B) parameter. In this case the stiffness of water is not taken into account.
Note that Ignore undrained behaviour does not affect materials of which the drainage
type is set to Undrained (C).
When changing the geometry configuration (Section 7.9) it is also possible (for each
calculation phase) to redefine the water boundary conditions and recalculate the pore
pressures (Section 7.8). For more details on theoretical formulations of a plastic
calculation reference should be made to the Scientific Manual.
In a Plastic calculation loading can be defined in the sense of changing the load
combination, stress state, weight, strength or stiffness of elements, activated by
changing the load and geometry configuration or pore pressure distribution by means of
Staged construction. In this case, the total load level that is to be reached at the end of
the calculation phase is defined by specifying a new geometry and load configuration,
and/or pore pressure distribution, in the Staged construction mode (Section 7.9).
The options for Pore pressure calculation type for a Plastic phase are:
• Phreatic
• Use pressures form previous phase
• Steady state groundwater flow
More information on Pore pressure calculation type is given in Section 7.5.
Hint: In PLAXIS, pore pressures are divided into steady-state pore pressures and
excess pore pressures. Steady state pore pressures are generated
according to the water conditions assigned to the soil layers for each phase,
whereas excess pore pressures are calculated as a result of undrained soil
behaviour (Undrained (A) or Undrained (B)) or consolidation. A Consolidation
calculation in PLAXIS only affects the excess pore pressures.
» A Consolidation calculation does not affect Undrained (C) materials.
available strength
SF = = value of ΣMsf at failure
strength at failure
The ΣMsf -value of a particular calculation step can be found in the Calculation
information window displayed as the corresponding option is selected in the Project
menu of the Output program. It is also recommended to view the development of ΣMsf
for the whole calculation using the Curves option (Chapter 10.2). In this way it can be
checked whether a constant value is obtained while the deformation is continuing; in
other words: whether a failure mechanism has fully developed. If a failure mechanism
has not fully developed, then the calculation must be repeated with a larger number of
additional steps.
To capture the failure of the structure accurately, the use of Arc-length control parameter
is required. The use of a Tolerated error of no more than 1% is also required. Both
requirements are complied with when using the default iteration parameters (Section
7.7.3).
Hint: When performing Safety calculation without Arc-length control, the reduction
factor ΣMsf cannot go down and an overestimation of safety factor can occur.
When using Safety calculation in combination with advanced soil models, these models
will actually behave as a standard Mohr-Coulomb model, since stress-dependent
stiffness behaviour and hardening effects are excluded from the analysis. In that case,
the stiffness is calculated at the beginning of the calculation phase based on the starting
stresses and kept constant until the calculation phase is completed. Note that when using
the Modified Cam-Clay model the strength is not reduced at all since this model does not
have a cohesion or friction angle as model parameter.
Hint: In case of the Jointed Rock model the strength on all the planes will be
reduced by ΣMsf.
» Strength in the Modified Cam-Clay model is not reduced in Safety analysis.
» When using Safety analysis in combination with user-defined soil models,
none of the parameters of these models will be reduced.
with
−σ '3
f σci a
f red = = mb +s
η η σci
and
v
u
u 1
− 1 f '2
u X
2
u
1 X
u Msf
η= Msf 2 − f ' t1 + + f '
u
2 2 − f ' 2
where
−σ '3
∂f
f' = = −amb mb + s a−1
∂σ '3 σci
More details and a derivation of the above equations can be found in Benz, Schwab,
Vermeer & Kauther (2007).
Updated mesh
The geometry of the model considered in a Safety calculation depends on whether the
Updated mesh option is selected or not in the parent phase. If the mesh is updated, the
resulting geometry at the end of the parent phase will be considered in the safety
calculations.
During a safety calculation the mesh is not updated at the beginning of each load step
even if the Updated mesh option is selected for the Safety phase.
The applied dynamic load is the product of the input value of the defined dynamic load
and the corresponding dynamic load multiplier.
Besides the activation of the dynamic load or dynamic prescribed displacement,
absorbent (viscous) boundary conditions can be defined for a Dynamic calculation. For a
more detailed description of the boundary conditions see Section 7.9.9.
Hint: Note that steady-state pore pressures in a Dynamic calculation are always
taken from the steady-state pore pressures generated in the parent phase
(see Section 7.5.1).
» It is possible to calculate excess pore pressures in undrained soil layers in a
dynamic analysis. However, the accuracy at which pore pressures are
generated depends on the capabilities of the soil models being used.
analysis of reinforced soil structures, the analysis of large offshore footing collapse
problems and the study of problems where soils are soft and large deformations occur.
When large deformation theory is included in a finite element program some special
features need to be considered. Firstly it is necessary to include additional terms in the
structure stiffness matrix to model the effects of large structural distortions on the finite
element equations.
Secondly, it is necessary to include a procedure to model correctly the stress changes
that occur when finite material rotations occur. This particular feature of large
displacement theory is usually dealt with by adopting a definition of stress rate that
includes rotation rate terms. Several stress rate definitions have been proposed by
researchers working in this field although none of these are wholly satisfactory. In
PLAXIS the co-rotational rate of Kirchhoff stress (otherwise known as the Hill stress rate)
is adopted. This stress rate would be expected to give accurate results provided that the
shear strains do not become excessive.
Thirdly, it is necessary to update the finite element mesh as the calculation proceeds.
This is done automatically within PLAXIS when the Updated mesh option is selected.
It should be clear from the descriptions given above that the updated mesh procedures
used in PLAXIS involve considerably more than simply updating nodal coordinates as the
calculation proceeds. These calculation procedures are in fact based on an approach
known as an Updated Lagrangian formulation (Bathe, 1982). Implementation of this
formulation within PLAXIS is based on the use of various advanced techniques that are
beyond the scope of this manual (van Langen, 1991).
The three basic types of calculations (Plastic, Consolidation and Safety) can optionally
be performed as an Updated mesh analysis, taking into account the effects of large
deformations. Therefore, the Updated mesh parameter should be selected. The
Updated mesh option cannot be used in a fully coupled flow-deformation analysis or a
dynamic calculation.
Please note that an updated mesh calculation cannot be followed by a 'normal'
calculation. Reversely, a normal calculation can be followed by an updated mesh
calculation, provided that the option Reset displacements to zero is used (Section 7.9).
It should be noted that an updated mesh analysis takes much more time and is less
robust than a normal calculation. Hence, this option should only be used in special cases.
Hint: Note that water conditions at the mesh boundary and pore pressures in
stress points are not updated as the mesh is updated.
Distributed loads
Distributed loads on deformed boundaries are taken into account as if those boundaries
were not deformed. This is to avoid that the total force involved does not change when
the boundary stretches or shrinks. This also applies to axisymmetric applications where
the radius changes as a result of deformation.
Calculation procedures
In order to carry out an updated mesh analysis the Updated mesh checkbox in the
Deformation control parameters subtree in the Phases window should be selected
(Section 7.9).
Updated mesh calculations are carried out using iteration procedures similar to the
conventional calculation options (Plastic or Consolidation) as described in preceding
sections. Therefore an updated mesh analysis uses the same parameters. However,
because of the large deformation effect, the stiffness matrix is always updated at the
beginning of a load step. Due to this procedure and to the additional terms and more
complex formulations, the iterative procedure in an updated mesh analysis is
considerably slower than that for conventional calculations.
Safety calculations
The geometry of the model considered in a Safety calculation depends on whether the
Updated mesh option is selected or not in the parent phase. If the mesh is updated, the
resulting geometry at the end of the parent phase will be considered in the safety
calculations.
During a safety calculation the mesh is not updated at the beginning of each load step
even if the Updated mesh option is selected for the Safety phase.
Practical considerations
Updated mesh analysis tends to require more computer time than an equivalent,
conventional calculation. It is recommended, therefore, that when a new project is under
study a conventional calculation is carried out before an updated mesh analysis is
attempted.
It is not possible to give simple guidelines that may be used to indicate when an updated
mesh analysis is necessary and where a conventional analysis is sufficient. One simple
approach would be to inspect the deformed mesh at the end of a conventional calculation
using the Deformed mesh option in the Output program. If the geometry changes are
large (on a real scale!) then significant importance of geometric effects might be
suspected. In this case the calculation should be repeated using the updated mesh
option. It cannot definitely be decided from the general magnitudes of the deformations
obtained from a conventional plasticity calculation whether geometric effects are
important or not. If the user is in any doubt about whether updated mesh analysis is
necessary then the issue can only be resolved by carrying out the updated mesh analysis
and comparing the results with the equivalent conventional analysis.
In general, it is not appropriate to use an updated mesh calculation for gravity loading to
set up the initial stress field. Displacements resulting from gravity loading are physically
meaningless and should therefore be reset to zero. Resetting displacements to zero is
not possible after an updated mesh analysis. Hence, gravity loading should be applied in
a normal plastic calculation.
Changing from a 'normal' plastic calculation or consolidation analysis to an updated mesh
analysis is only valid when displacements are reset to zero, because a series of updated
mesh analyses must start from an undeformed geometry. Changing from an updated
mesh calculation to a 'normal' plastic calculation or consolidation analysis is not valid,
The loading type is specified in the Loading input drop-down menu in the Phases window.
Only one of the described loading types can be activated in any single calculation phase.
The available loading types depend on the selected type of calculation.
Hint: Note that the Loading type drop-down menu is NOT available for
Groundwater flow only calculation type.
Hint: The Staged construction loading type is available for Plastic, Consolidation
and Dynamic calculations.
Hint: The Minimum pore pressure loading type is available for Consolidation
calculations.
where Pmax is the maximum excess pore pressure reached in the previous phase which
can be found in the the Reached values subtree in the Phases window (Section 7.19).
The calculation stops when the maximum absolute excess pore pressure is below this
calculated value of Minimum pore pressure. Note that the maximum number of additional
steps as defined by the parameter Max steps will not be reached if the Minimum pore
pressure criterion is met before.
available strength X
SF = = value of Msf at failure
strength at failure
The ΣMsf -value of a particular calculation step can be found in the Calculation
Information window of the Output program. It is also recommended to view the
development of ΣMsf for the whole calculation using the Curves option (Chapter 10). In
this way it can be checked whether a constant value is obtained while the deformation is
continuing; in other words: whether a failure mechanism has fully developed. If a failure
mechanism has not fully developed, then the calculation must be repeated with a larger
number of additional steps.
Hint: The Incremental multipliers loading type is available for Safety calculations.
Water pressure can be 'external' water pressure (i.e. 'water load' on model boundaries) or
'internal' water pressure, which is known as pore water pressure. Pore water pressure is
included in what is denoted in PLAXIS as active pore pressure. Active pore pressure is
composed of steady-state pore pressure and excess pore pressure, and may also include
suction (positive pore water stress). In general, the steady-state component of water
pressure (both the 'external' and 'internal' part) is considered to be input data and
supposed to be known at the beginning of a deformation analysis, whereas excess pore
pressure is the result of undrained loading or consolidation. This section deals with the
calculation of steady-state water pressure as input for a deformation analysis, which can
be generated according to the options available in the Phases window.
the global water level, but the 'internal (phreatic) ' part will most likely change.
Hint: If no surface groundwater flow boundary conditions are defined, the hydraulic
boundary conditions for a steady-state groundwater flow calculation are
taken from the Global water level, which means a prescribed hydraulic head
below the water level and seepage above.
» The Steady state groundwater flow option is available for K0 procedure,
Gravity loading, Groundwater flow only, Plastic and Consolidation calculation
types.
Steady-state pore pressure may include suction in the unsaturated zone above the
phreatic surface. If it is desired to exclude suction from the steady-state pore pressure, as
input for a deformation analysis, the option Ignore suction may be used.
When soil plasticity is involved in a finite element calculation the equations become
non-linear, which means that the problem needs to be solved in a series of calculation
steps. An important part of the non-linear solution procedure is the choice of step size
and the solution algorithm to be used.
During each calculation step, the equilibrium errors in the solution are successively
reduced using a series of iterations. The iteration procedure is based on an accelerated
initial stress method. If the calculation step is of a suitable size then the number of
iterations required for equilibrium will be relatively small, usually around ten.
If the step size is too small, then many steps are required to reach the desired load level
and computer time will be excessive. On the other hand, if the step size is too large then
the number of iterations required for equilibrium may become excessive or the solution
procedure may even diverge.
PLAXIS has an automatic load stepping procedure for the solution of non-linear plasticity
problems. Users do not need to worry about the proper selection of load steps and
numerical procedures, since the program will automatically use the most appropriate
procedure to guarantee optimum performance. The parameters for the load stepping
procedures are available in the Numerical control parameters subtree in the Phases
window (Figure 7.7).
• PLAXIS performs a trial calculation step and determines a suitable step size on the
basis of this trial.
• PLAXIS sets the initial load step size to be equal to the final load step size of any
previous calculation.
The first method is generally adopted. The second method would only be used if the
loading applied during the current load step is similar to that applied during the previous
load step, for example if the number of load steps applied in the previous calculation
proved to be insufficient.
In subsequent steps, the automatic load stepping procedures are adopted (Section
7.6.1). If at the end of the calculation, the defined state or load level has been reached,
the calculation is considered to be successful. A successful calculation is indicated by a
check mark in a green circle in the Phases explorer and the Phases window.
If the defined state or load level has NOT been reached, the calculation is considered to
have failed. A failed calculation is indicated by a cross mark in a red circle in the Phases
explorer or the Phases window. A message describing the error is given in the Log info
for the last calculation box in the Phases window:
Prescribed ultimate state not reached; Soil body collapses: A collapse load has
been reached. In this case, the total specified load has not been applied. Collapse is
assumed when the applied load reduces in magnitude in X successive calculation steps
(where X is the maximum number of unloading steps; see Section 7.7.3) and the current
stiffness parameter CSP is less than 0.015 (see Section 7.10.8 for the definition of CSP).
It is also possible that the problem is failing but due to switched-off arc-length control, the
program is not allowed to take negative step sizes. The user should check the output of
the last step and judge whether the project is failing or not. In case of failure,
recalculating the project with a higher Additional steps parameter is useless.
Prescribed ultimate state not reached; load advancement procedure fails. Try
manual control: The load advancement procedure is unable to further increase the
applied load, but the current stiffness parameter CSP is larger than 0.015. In this case
the total load specified has not been applied. The user can now attempt to rerun the
calculation with slight changes to the iterative parameters in Numerical control
parameters subtree in the the Phases window, in particular turning off the Arc-length
control type parameter.
Prescribed ultimate state not reached; Not enough load steps: The maximum
specified number of additional load steps have been applied. In this case, it is likely that
the calculation stops before the total specified load has been applied. It is advised to
recalculate the phase with an increased value of Max steps.
Cancelled by user: This occurs when the calculation process is terminated by clicking
Stop in the Active tasks window.
Prescribed ultimate state not reached; Numerical error: A numerical error has
occurred. In this case, the total specified load has not been applied. There may be
different causes for a numerical error. Most likely, it is related to an input error. Careful
inspection of the input data, the finite element mesh and the defined calculation phase is
suggested.
Severe divergence: This is detected when the global error is increasing and has
reached huge values. This error, for example, can be caused by very small time steps in
a consolidation phase. The program scales down the step size when the tolerated error
cannot be reached, resulting in small time steps. One of the reasons can be that a failure
situation is reached. As for consolidation the arc-length procedure is not used, the
program cannot really detect failure.
File xxxx not found: Such a message appears when a file that ought to exist does not
exist.
Messages may indicate errors related to the iterative solution algorithm or the matrix
condition. In the case of 'floating' elements (insufficient boundary conditions), one could
get a message indicating that the matrix is nearly singular. Checking and improving the
defined calculation phase usually solves the problem.
Another problem related to the solution process can occur due to insufficient direct
computer memory (RAM). In such cases, the iterative solver cannot store the minimum
amount of data necessary to have sufficient accuracy. As a result, the iteration process
converges very slowly, or does not reach the accuracy condition. To solve this problem,
either the problem size should be reduced, or internal memory capacity should be
increased.
The control parameters of a particular calculation phase and the corresponding solution
procedure can be defined in the Phases window.
Phase identification
The ID of the phase consist of the caption and the name (in square brackets). The name
of the phases is determined consecutively by the program and it can not be modified by
the user. The user may redefine the caption of the phase in the Phases window.
Calculation type
The calculation type of the selected phase can be defined in the corresponding
drop-down menu in the Phases window. The options available are:
K0 procedure
Gravity loading
Plastic
Consolidation
Fully coupled analysis
Safety
Dynamic
A description of the different analysis types available in PLAXIS is given in Section 7.3.
Loading type
The options available for the Loading type parameter are:
Staged construction
Incremental multipliers
Target SumMsf
Minimum pore pressure
Degree of consolidation
A description of the different loading types available in PLAXIS is given in Section 7.4.
Σ Mweight
ΣMweight is the total multiplier for the material weight. For ΣMweight = 1, the material
weight as specified by the unit weights in the material data sets is applied. In general,
ΣMweight remains at its default value of 1, but there are situations in which ΣMweight may
be changed:
• In the case of modelling simplified soil tests by means of finite element models, self
weight of the material may be disregarded since the stresses are dominated by
external loads rather than material weight. This can either be achieved by using zero
unit weights in the material data set or by using ΣMweight = 0 during the calculations.
• In the case of over-consolidated material in situations where the K0 procedure
cannot be used (for example slopes and non-horizontal ground surfaces), gravity
loading can be used in the initial phase with an increased value of ΣMweight equal to
the over-consolidation ratio. In this way, the pre-consolidation stress is properly
initialised, at least when using advanced soil models. In the subsequent phase,
ΣMweight should be set back to 1.0 in order to have realistic initial soil stresses whilst
the pre-consolidation stress (in advanced models) will memorize the
over-consolidated stress level.
• In the case of modelling a centrifuge test, ΣMweight may be used to simulate the
increased amount of gravity. Hence, in order to simulate a centrifuge test at 100 g,
ΣMweight should be set to 100.
Time interval
A non-zero value for the time-related parameters is only relevant when a transient
groundwater flow calculation, a consolidation analysis, a fully coupled flow-deformation
analysis or a dynamic analysis is performed or when using time-dependent material
models (such as the Soft Soil Creep model).
Time interval To define the total time period considered in the current
calculation phase, expressed in the unit of time as defined in the
Project properties window.
Dynamic time interval To define the total time period considered in the current dynamic
calculation phase, expressed in seconds [s].
A dynamic analysis uses a different time parameter than other types of calculations. The
time parameter in a dynamic analysis is the Dynamic time, which is always expressed in
seconds [s], regardless of the unit of time as specified in the Project properties window.
In a series of calculation phases in which some of them are dynamic, the Dynamic time is
only increased in the dynamic phases (even non-successive), while the Dynamic time is
kept constant in other types of calculations (whether before, in-between or after the
dynamic phases).
The Dynamic time is not affected by the regular time parameter. Reversely, the regular
time parameter includes the Dynamic time.
Hint: The Ignore undrained behaviour option is not available for a Consolidation
analysis or a Fully coupled flow-deformation analysis, since these calculation
types do not consider the Drainage type as specified in the material data
sets, but use the material permeability instead.
Updated mesh
This option must be selected when the calculation should be performed as a large
deformation analysis according the Updated Lagrange formulation. For a more detailed
Table 7.1 The material models and the corresponding reset parameters
Ignore suction
When steady-state pore pressures are generated by the phreatic level option or by a
groundwater flow calculation, tensile pore water stresses (suction) will be generated
above the phreatic level. Although suction is a realistic phenomenon, the use of suction in
a deformation analysis may lead to an increased shear strength when effective strength
parameters are used for the soil. In order to avoid this, suction can be cut off by selecting
the Ignore suction option. Note that the Ignore suction option does not affect the excess
pore pressures generated in undrained soil layers, except in the case of a fully coupled
flow-deformation analysis. The default setting is such that Ignore suction is selected.
When Ignore suction is selected, soil below the phreatic level is considered to behave
fully saturated, whereas soil above the phreatic level is considered to behave as dry soil,
although excess pore pressures may still occur here. Positive values of steady-state pore
stresses will be set to zero. However, any excess pore pressure above the phreatic level,
both positive and negative, will be taken into account. This requires the effective degree
of saturation, Seff , to be set to 1 (except in the case of a fully coupled flow-deformation
analysis. Hence, positive excess pore stresses will always be taken into account to their
full extent until the cavitation cut-off is reached, except in the case of a fully coupled
flow-deformation analysis.
When Ignore suction is NOT selected, suction is allowed and included in the pore water
stresses (both steady state and excess pore stresses). In this case the effective degree
of saturation, Seff , determines the proportion of the suction that is included in the active
pore stresses. This depends on the soil-water retention curve as defined in the material
data set of the soil layer.
Cavitation cut-off
In case of unloading of undrained materials (undrained A or B) tensile excess pore
stresses may be generated. These excess pore stresses might give rise to tensile pore
water stresses. In case the cavitation cut-off option is activated, excess pore pressures
are limited so that the tensile pore water stress is never larger than the cavitation stress.
By default, the cavitation cut-off option is not activated. If it is activated, the default
cavitation stress is 100 kN/m2 . The Cavitation cut-off option is not available for a fully
coupled flow-deformation analysis.
Max steps
This parameter specifies the maximum number of calculation steps (load steps) that are
performed in a particular calculation phase.
If Plastic, Consolidation or Fully coupled flow-deformation is selected as the calculation
type, then the number of additional steps should be set to an integer number
representing the required number of steps for this calculation phase. In this case the
defined value is an upper bound to the actual number of steps that will be executed. In
general, it is desired that such a calculation is completed within the defined number of
steps and stops when either the prescribed ultimate state is reached or the soil body
collapses. If such a calculation reaches the maximum number of steps, it usually means
that the ultimate level has not been reached. By default, the Max steps parameter is set
to 250, but this number can be changed within the range 1 to 10000.
If Safety or Dynamic is selected as the calculation type, then the number of additional
steps is always exactly executed. In general, it is desired that such a calculation is
completed within the defined number of steps and stops when either the prescribed
ultimate state is reached or the soil body collapses. If such a calculation reaches the
maximum number of steps, it usually means that the ultimate level has not been reached.
By default, the Max steps parameter is set to 100 for Safety and 250 for Dynamic
calculations which is generally sufficient to complete the calculation phase. However, this
number may be changed within the range 1 to 10000.
Figure 7.11 Dialog box displaying the retrieved value for the First time step parameter
used for the Number of sub steps, this parameter is automatically calculated in the kernel.
To be able to see and change these two parameters, the user has the possibility to
retrieve them by pressing the Retrieve button. By pressing the Apply button, these values
will be applied and will be used by kernel without any change. It should be noted that
PLAXIS always tries to find the closest number of steps to the Max steps that the user
specifies (in case of using the default option, the number of steps will be the closest
number to 250 as shown in Figure 7.12).
Figure 7.12 Dialog box displaying the retrieved value for the Number of sub steps parameter
Tolerated error
In any non-linear analysis where a finite number of calculation steps are used there will
be some drift from the exact solution, as shown in Figure 7.13. The purpose of a solution
algorithm is to ensure that the equilibrium errors, both locally and globally, remain within
acceptable bounds (Section 7.10.8). The error limits adopted in PLAXIS are linked
closely to the specified value of the Tolerated error.
Within each step, the calculation program continues to carry out iterations until the
calculated errors are smaller than the specified value. If the tolerated error is set to a high
value then the calculation will be relatively quick but may be inaccurate. If a low tolerated
error is adopted then computer time may become excessive. In general, the standard
setting of 0.01 is suitable for most calculations.
numerical solution
load
exact solution
displacement
Figure 7.13 Computed solution versus exact solution
If a calculation gives failure loads that tend to reduce unexpectedly with increasing
Hint: Be careful when using a tolerated error larger than the default value of 0.01,
as this may give inaccurate results which are not in equilibrium.
displacement, then this is a possible indication of excessive drift of the finite element
results from the exact solution. In these cases, the calculation should be repeated using
a lower value of the tolerated error. For further details of the error checking procedures
used in PLAXIS see Section 7.10.8.
Most Validation & Verification examples have been performed with a tolerated error lower
than the default settings.
Over-relaxation factor
To reduce the number of iterations needed for convergence, PLAXIS makes use of an
over-relaxation procedure as indicated in Figure 7.14. Over-relaxation is a deliberate
over-estimation of the equilibrium error under the assumption that by just solving the
exact unbalance, the next iteration is still far from equilibrium. The purpose of this is to
reduce the number of the iterations in order to speed up the calculation. The parameter
that controls the degree of over-relaxation is the over-relaxation factor. The theoretical
upper bound value is 2.0, but this value should never be used. For low soil friction angles,
for example ϕ < 20◦ , an over-relaxation factor of about 1.5 tends to optimise the iterative
procedure. If the problem contains soil with higher friction angles, however, then a lower
value may be required. The standard setting of 1.2 is acceptable in most calculations.
load
load
displacement displacement
ation
Figure 7.14 Influence of over-relaxation
that take place. This parameter is required only to ensure that computer time does not
become excessive due to errors in the specification of the calculation. The standard
value of Maximum number of iterations is 60, but this number may be changed within the
range 1 to 100.
If the maximum allowable number of iterations is reached in the final step of a calculation
phase, then the final result may be inaccurate. If this is the case then the message
'Maximum iterations reached in final step' is displayed in the Log info for last calculation
box of the Phases window. Such a situation occasionally occurs when the solution
process does not converge. This may have various causes, but it mostly indicates an
input error.
If the soil friction angles are relatively high, or if high-order soil models are used, then it
may be appropriate to increase the desired minimum and maximum from their standard
values to obtain a solution without the use of excessive computer time. In these cases
the following values are suggested:
Arc-length control
The Arc-length control procedure is a method that is by default selected in PLAXIS to
obtain reliable collapse loads for load-controlled calculations. The iterative procedure
adopted when arc-length control is not used is shown in Figure 7.7.3 for the case where a
collapse load is being approached. In the case shown, the algorithm will not converge. If
arc-length control is adopted, however, the program will automatically evaluate the
portion of the external load that must be applied for collapse as shown in Figure 7.7.3.
load load
step 3 step 3 arc
step 2 step 2
step 1 step 1
displacement displacement
a. load control b. arc-length control
Figure 7.15 Influence of Arc-length control
Arc-length control is activated by selecting the corresponding check box in the Numerical
control parameters subtree. The arc-length control procedure should be used for
load-controlled calculations, but it may be deactivated, if desired, for
displacement-controlled calculations.
Setting the Arc-length control type parameter to Auto enables the procedure only when
the CSP (current stiffness parameter) is less than 0.5. In this way arc-length is only used
when the material exhibits significant plastic behaviour.
number of alternative (but still correct) solutions with the same tight tolerance, but with,
for instance, different step size. Hence, is it not really necessary to enforce a tight
tolerance during the whole calculation, as a feasible solution is located in a certain band,
defined by all alternative paths. When the gradual error reduction procedure is enabled,
the solution does not follow one of the paths, but, possibly, jumps from one to another.
However, in the end of the calculation the error is forced to reduce to the levels, as
defined for the calculation phase.
Nonetheless, application of this procedure should be done with caution. If the initial
tolerated error is too high, hypothetically, the solution can switch to a path that is not even
connected to the current initial conditions (hence infeasible). One way to prevent this is to
use much smaller allowable Max load fraction in one step than the default value used in
PLAXIS. The maximum step size is 50% of the whole stage load. The user might
consider reducing the step size to something like 2% or less, hence requiring at least 50
steps (100% / 2%) to obtain the solution of the calculation phase.
Extrapolation
It is a numerical procedure, which is automatically used in PLAXIS if applicable, when a
certain loading that was applied in the previous calculation step is continued in the next
step. In this case, the displacement solution to the previous load increment can be used
as a first estimate of the solution to the new load increment. Although this first estimate is
generally not exact (because of the non-linear soil behaviour), the solution is usually
better than the solution according to the initial stress method (based on the use of the
elastic stiffness matrix) (Figure 7.16).
load
load
displacement displacement
a. without extrapolation b. with extrapolation
Figure 7.16 Difference between elastic prediction (a) and extrapolation from previous step (b)
After the first iteration, subsequent iterations are based on the elastic stiffness matrix, as
in the initial stress method (Zienkiewicz, 1977). Nevertheless, using Extrapolation the
total number of iterations needed to reach equilibrium is less than without extrapolation.
The extrapolation procedure is particularly useful when the soil is highly plastic.
Tolerated error and Over-relaxation factor whereas for a Transient groundwater flow
analysis a higher number of parameters is required to tdefine the iterative procedure.
Max steps
This parameter specifies the maximum number of calculation steps that are performed in
a phase which Pore pressure calculation type is either Steady state groundwater flow or
Transient groundwater flow. By default, the Max steps parameter is set to 1000.
Tolerated error
Similar to the Tolerated error option in deformation analyses (Section 7.7.3), the
Tolerated error option in groundwater flow is used to check if the result is accurate
enough or not. In case of unsaturated groundwater flow, checking the local error (error
which occurs in few elements) is also needed. In this case, due to the fact that the
permeability can change significantly, it is possible to have isolated water in some
elements in unsaturated zones. The specified Tolerated error is used for checking both
global flow error and local errors. PLAXIS uses higher Tolerated error to check the local
flow errors (usually 10 times larger value). If both global and local errors are less than the
Tolerated error, the groundwater flow is terminated.
Over-relaxation factor
Since the unsaturated groundwater flow analyses may be highly non-linear, PLAXIS uses
an over-relaxation factor to speed up the steady state calculations and to decrease the
number of steps needed. This factor is by default 1.5. For confined flow analyses which
are linear a factor of 1 is enough. In case of highly non-linear analyses for which PLAXIS
needs more steps, higher values than 1.5 might result in faster convergence.
For an average acceleration scheme you can use the standard settings (α = 0.25 and
β = 0.5). Using a higher β -value and corresponding α-value results in a damped
Newmark scheme (e.g. α = 0.3025 and β = 0.6).
Reached total time It is the actual accumulated time at the end of a finished
calculation phase.
CSP - Relative stiffness It is the relative stiffness parameter calculated at the end of the
phase, which is a measure for the amount of plasticity that
occurs during the calculation (Section 7.10.8). A value of 1.0
means that the full model is elastic; a value close to zero
indicates that failure is approached.
Reached total force They are the component of the resulting reactions on the nodes
where non-zero prescribed displacement is applied.
Pmax - Reached max pp
It is the value of the maximum pore pressure that has been
actually reached at the end of a calculation phase. This value
corresponds to the maximum excess pore pressure for a
Consolidation analysis, the steady state pore pressure for a
groundwater flow analysis and the maximum active pore
pressure for a Fully coupled flow-deformation analysis.
ΣMstage - Reached phase proportion
It is an indication of the proportion of the unbalance solved in a
phase where the Staged construction loading type is applied
(Section 7.4.1).
ΣMweight - Reached material weight proportion
It is an indication of the total proportion of the material weight
applied in a calculation (Section 7.7.1). A reached value of 1.00
means that the specified weight of soil and structural materials is
fully applied.
ΣMsf - Reached safety factor
It is the resulting value of the ΣMsf parameter in a Safety phase
(Section 7.3.5).
PLAXIS is generally used for effective stress analysis in which total stresses are divided
into effective stresses, σ ', and active pore pressures, pactive .
σ = σ ' + pactive
Active pore pressure (pactive ) is defined as the effective saturation, Seff , times the pore
water pressure, pwater .
Pore water pressure differs from active pore pressure when the degree of saturation is
less than unity. PLAXIS can deal with saturated soil below the phreatic level, as well as
with partially saturated soil above the phreatic level. More details about partially saturated
soil behaviour are provided in Section 6.1.3.
In the pore water pressure a further division is made between steady state pore pressure,
psteady , and excess pore pressure, pexcess .
Excess pore pressures are pore pressures that occur as a result of stress changes in
undrained materials. In this respect, changes in stress may be a result of loading,
unloading, a change in hydraulic conditions or consolidation. Hence, excess pore
pressures are a result of a deformation analysis. In a Plastic calculation, a Safety
analysis or a Dynamic analysis, excess pore pressures can occur in clusters for which the
Drainage type is set to Undrained A or Undrained B. In a Consolidation analysis and a
Fully coupled flow-deformation analysis, excess pore pressures can occur in any material
(except Non-porous materials), depending on the permeability as defined in the
corresponding material data set. In the case of a Fully coupled flow-deformation analysis,
excess pore pressure is calculated by subtracting the steady-state pore pressure, psteady ,
from the pore water pressure, pwater .
Steady-state pore pressures are pore pressures that represent a stable state. Since
steady-state pore pressures are not supposed to change during a deformation analysis,
they are considered input data. There are various ways to define and generate
steady-state pore pressures. The way to indicate how pore pressures are generated for a
particular calculation phase is indicated by means of the Pore pressure calculation type
parameter in the Phases window. In contrast to other deformation calculation types, a
Fully coupled flow-deformation analysis does not allow for the selection of the Pore
pressure calculation type, because in this case total pore water pressures, pwater , are
calculated together with displacements. To enable a division in steady-state pore
pressures and excess pore pressures, the steady-state pore pressures, psteady , are
automatically calculated on the basis of a preliminary steady-state groundwater flow
calculation using the hydraulic boundary conditions at the end of the calculation phase.
This enables the calculation and output of excess pore pressures, pexcess , in all steps:
The definition of water conditions is needed to generate pore pressures in the soil and to
generate external water pressures, if applicable. The definition of water conditions can be
divided into the creation of water levels and assigning water conditions to volume
clusters. These features will be described in more detail below.
The continuity of the flow functions assigned to water levels in consecutive phases can
be maintained by using the options available in the pop up menu displayed when the
water level is right-clicked. These options are:
Create the reached water level
A non-horizontal water level can be created by defining more points in the geometry,
while the Create water level option is still active. The creation of water levels is basically
the same as the creation of surfaces (Section 5.2.4). If necessary, the point coordinates
or the pressure increment with depth (Pinc) may be changed in the table. Unlike
surfaces, user water levels can only be moved in vertical direction by dragging them with
the mouse. A new user water level may also be created by duplicating an existing water
level and modifying it afterwards. To duplicate an existing water level, right click it and
select Duplicate from the appearing menu. User water levels can be deleted by
right-clicking them and selecting the Delete option from the appearing menu.
Time dependancy can be specified for horizontal User water levels. The time dependent
behaviour can be specified by selecting the corresponding flow function.
Global level
When the Global level option is selected, the pore pressure distribution in the soil layer is
generated according to the global water level of the phase considered. The resulting pore
pressure distribution in the soil layer is then hydrostatic. This does not mean that the pore
pressure distribution in the soil as a whole is hydrostatic, since other soil layers may have
other hydraulic conditions.
Custom level
When the Custom level option is selected in the Water conditions drop-down menu, a
further selection should be made from the list of available water levels (Borehole water
levels and User water levels). The resulting pore pressure distribution in the soil layer is
then hydrostatic.
Head
The Head option is a quick alternative for the generation of hydrostatic pore pressures
based on a horizontal water level. After selecting the Head option in the Water conditions
drop-down menu, a reference level (zref ) should be defined, indicating the level where the
pore pressure will be zero.
Hint: The multi-select option may be used to efficiently set water conditions for
multiple soil layers at once.
» When a non-horizontal water level is used to generate a hydrostatic pore
pressure distribution, this distribution may not be fully realistic, since in reality
a non-horizontal water level would be the result of a groundwater flow
process in which case the pore pressure distribution may be non-hydrostatic.
User-defined
If the pore pressure distribution is linear but not hydrostatic, the User-defined option
should be selected in the Water conditions drop-down menu. Subsequently, values for
the reference level (zref ), the reference pressure (pref ) and the pressure increment (pinc )
should be entered. Note that pressure is entered as a negative value.
Interpolate
It is possible to generate pore pressures in a soil volume based on the pore pressures of
the adjacent soil volumes (above and below) by selecting the Interpolate option in the
Condition drop-down menu. This option is, for example, used if a relatively impermeable
layer is located between two permeable layers with a different groundwater head. The
pore pressure distribution in the relatively impermeable layer will not be hydrostatic, so it
cannot be defined by means of a phreatic level.
On selecting the option Interpolate the pore pressure in that soil volume is interpolated
linearly in vertical direction, starting from the value at the bottom of the soil volume above
and ending at the value at the top of the soil volume below. The Interpolate option can be
used repetitively in two or more successive soil volumes (on top of each other). In the
case that a starting value for the vertical interpolation of the pore pressure cannot be
found, then the starting point will be based on the Global water level.
Dry
In addition to the global water pressure distribution it is possible to remove water
pressures from individual groups of volumes in order to make them 'dry'. Deactivation of
water can be done independent from the soil itself. Hence, if the soil is deactivated and
the water level is above the excavation level, then there is still water in the excavated
area. If it is the user's intention to simulate a dry excavation, then the water must be
explicitly deactivated. Please note that the water pressure in adjacent soil volumes is not
affected and may need to be changed manually.
Hint: The Dry clusters behave as non-porous materials. As a result, neither initial
nor excess pore pressures is taken into account and flow is not possible
through the cluster.
Boundary conditions-Phases
Besides the options described in Section 5.7.4 for surface groundwater flow boundary
conditions, two other options are available in Calculation modes. These options are:
Constant value from previous phase
The water level reached at the end of the previous phase will be
kept constant during the current phase.
Maintain function from previous phase
The water level in the current phase will continue to vary
according to the flow function assigned to the boundary in the
previous phase. The accumulated time from previous phases is
used as an offset to enter the flow function for the current phase.
7.8.5 PRECIPITATION
The Precipitation option can be used to specify a general vertical recharge or infiltration
(q ) due to weather conditions. This condition is applied at all boundaries that represent
the ground surface. This option can be selected in the Model conditions subtree in the
Model explorer (Figure 7.22).
Hint: Note that when both Precipitation and a surface groundwater flow boundary
conditions is created at the top boundary of the model, the condition specified
for the boundary condition is taken into account ignoring Precipitation.
In PLAXIS 3D, it is possible to change the geometry and load configuration by activating
or deactivating loads, soil volume clusters or structural objects created in the geometry
input for each calculation phase. The material data sets can be reassigned and the water
pressure distribution in the geometry can be changed as well.
Changes to the geometry configuration generally cause substantial out-of-balance forces.
These out-of-balance forces are stepwise applied to the finite element mesh using the
automatic load stepping procedures (Section 7.6.1).
Hint: The option to activate or deactivate objects should not be confused with the
option to show or hide objects. A hidden object (invisible) may still be active.
» The visibility of deactivated objects is defined in the Visibility tabsheet of the
Visualisation settings window (Section 3.5.3).
compression). Starting with a single layer and adding one layer on top of the first will
give settlements of the top surface. If a third layer is subsequently added to the
second layer, it will be given an initial deformation corresponding to the settlements
of the surface.
The following rules apply for interfaces:
Deactivated interfaces The interfaces can be deactivated when modelling soil-structure
interaction is not desired. Nodes generated for interfaces in
mesh generation process are still there. They have stiff elastic
behaviour; fully permeable (Consolidation analysis).
Activated interfaces Elasto-plastic behaviour; fully impermeable (Consolidation
analysis).
A connection is active when both parts are active in a calculation phase. Modifications
done in connections are valid for all the phases where the connection is active. It is not
possible to assign different properties to a connection in different phases.
If a connection has been defined, by default, the rotation degree of freedom of the 'child's'
node(s) involved in the connection is disconnected (free) from the corresponding degree
of freedom in the 'Parent's' node(s), whereas the translational degrees of freedom are
tied (fixed). This means that the default connection is a hinge. The user can modify a
connection by defining for each 'child' degree of freedom whether it is free or fixed.
Distinction is made between rotation around the first local axis and translation in the three
individual local directions. In order to define a dilation joint, at least one translational
degree of freedom should be set free.
Hint: The order of selection effects the part of the connection assigned to each
plate. The first selected plate will be assigned as the Custom part of the
connection whereas the second selected part will be assigned as the
Reference part of the connection.
» The role of the parts of the connection can be interchanged by right clicking
the connection in the Object explorer and selecting the Invert connected part
option from the appearing menu.
» The order of selection is not important in a plate-beam connection. In that
case the plate will always be the reference part of the connection.
or in the Selection explorer. Values can be defined for the three strain components in x -,
y - and z -direction and for volumetric strain. A positive value of the strain component
represents an expansion, whereas a negative value represents a shrinkage in that
direction.
Deformation
PLAXIS automatically imposes a set of general fixities to the boundaries of the geometry
model. These conditions are generated according to the following rules.
Soil volume:
• Vertical model boundaries with their normal in x -direction (i.e. parallel to the
yz -plane) are fixed in x -direction (ux = 0) and free in y - and z -direction.
• Vertical model boundaries with their normal in y -direction (i.e. parallel to the
xz -plane) are fixed in y -direction (uy = 0) and free in x - and z -direction.
• Vertical model boundaries with their normal neither in x - nor in y -direction are fixed
in x - and y -direction (ux = uy = 0) and free in z -direction.
• The model bottom boundary is fixed in all directions (ux = uy = uz = 0).
• The 'ground surface' is free in all directions.
Structures: Beams and plates that extend to the model boundary (xmin , xmax , ymin ,
ymax ) where at least one displacement direction is fixed obtain at least two fixed rotations
in the points at the boundary.
• At vertical model boundaries with a normal in x -direction: φy = φz = 0 (φx = free).
Hint: Note that the standard displacement fixities (Section 7.9.9) apply to soil
volumes and structures.
» User-defined prescribed displacements always prevail over the automatically
defined fixities. This means that if a prescribed surface displacement is
placed on the top surface and this surface displacement shares an edge with
the side surface of the model (e.g. in the modelling of a symmetrical smooth
footing), there will be no horizontal fixity on the shared edge. In this case,
add a line prescribed displacement along this edge and assign the proper
fixity conditions to it.
Dynamics
A viscous boundary is aimed to absorb the increments of stresses on the boundaries
caused by dynamic loading, that otherwise would be reflected inside the soil body. When
a Dynamic calculation is defined, the boundary conditions can be defined for the phase
by selecting the proper options displayed when the Dynamics subtree is expanded under
Model conditions in Model explorer (Figure 7.27).
The Relaxation coefficients C1 and C2 are used to improve the wave absorption on the
viscous boundaries. C1 corrects the dissipation in the direction normal to the boundary
and C2 does this in the tangential direction. If the boundaries are subjected only to waves
that come in perpendicular to the boundary, relaxation is not necessary (C1 = C2 = 1).
When there are waves in arbitrary direction (which is normally the case), C2 has to be
adjusted to improve the absorption. The standard values are C1 = 1 and C2 = 1.
Such boundary conditions are required for a groundwater flow calculation, a consolidation
analysis or a fully coupled flow-deformation analysis. More information regarding the flow
boundary conditions is available in Section 7.8.4.
Precipitation
The Precipitation subtree in the Model explorer (Figure 7.29) can be used to specify a
general vertical recharge or infiltration (q) due to weather conditions. This condition is
applied at all boundaries that represent the ground surface. More information on
Precipitation is available in Section 7.8.5.
Water
The global water level in a calculation phase can be defined in the Water subtree (Figure
7.30). More information on Global water level is available in Section 7.8.2.
When a calculation phase has been defined, its calculation can be started.
direct visual check before the calculation is started. The active part of the 3D model is
presented in the Output program.
See Chapter 8 for more details on the inspection of 3D models in the Output program.
After the preview, press the Close button to return to the Calculation mode.
Kernel information
Start time The time indicating the start of the calculation is displayed.
Memory used The memory occupied by the calculation process is displayed.
Calculation progress
A small load-displacement curve for one of the pre-selected nodes for curves is shown in
the Calculation progress group box. By default, the curve is shown for the first selected
node. Curve for other pre-selected nodes is shown when the node is selected in the
drop-down menu. The presented graph may be used to roughly evaluate the progress of
the calculation.
Plastic analysis For a plastic analysis the development of the
P
Mstage
parameter is plotted against the displacement.
Consolidation analysis In case of a Consolidation analysis, the maximum excess pore
pressure, Pexcess,max is plotted against the logarithm of time.
Fully coupled flow-deformation analysis
In case of a Fully coupled flow-deformation analysis analysis, the
maximum active pore pressure, Pactive,max , is plotted against the
logarithm of time.
Safety analysis In case of Safety analysis, the development of Msf is plotted
P
against the displacement.
Dynamic analysis In case of Dynamic analysis, the displacement is plotted against
the dynamic time.
Groundwater flow (steady-state)
In case of Groundwater flow (steady-state), the maximum steady
pore pressure, Psteady,max , at steady state is plotted.
Groundwater flow (transient)
In case of Groundwater flow (transient), the maximum steady
pore pressure, Psteady,max , is plotted against the logarithm of
time.
Calculation status
The calculation status indicates what part of the calculation process is currently being
executed. The following processes are indicated:
Reading data. . . Reading input data from the disk.
Renumbering. . . Optimising the node numbering and determining matrix
properties.
Profile. . . Determining the profile of the global stiffness matrix or
pre-conditioner.
Forming matrix. . . Forming the global stiffness matrix.
Forming pre-conditioner. . .
Forming the pre-conditioner for the iterative solution procedure.
Solving equations. . . Solving the global system of equations to obtain the
displacement increments.
Calculating stresses. . . Calculating the strain increments and constitutive stresses.
Reaction forces. . . Calculating the reaction forces and the out-of-balance force.
Writing data. . . Writing output data to disk.
Note that when the calculation of the phase is completed, a warning will appear indicating
that the intermediate results are no longer available. A more detailed description on how
to display the results of a calculated phase is given in Section 7.10.6.
Modification of geometry
When the geometry is slightly modified (small relocation of objects, slight modification of
their geometry or deleting objects) in the Geometry modes, the model is automatically
regenerated in the Calculation modes. The projects should be re-meshed. Note that the
phases in the Phases explorer are marked for calculation automatically. The deleted
geometry entities will not be any more available in the calculation phases. However, the
program preserves the definition of the phases.
If significant changes in the geometry are made then all settings need to be redefined,
since PLAXIS is not able to properly regenerate the settings automatically.
In the determination step by definition the properties of the target feature in the target
phase are included. The properties of the target feature are examined in all child phases
of the target. If the properties of target feature in a given child phase are identical to those
in target phase, those properties will be included in the regeneration, and the process will
repeat for all child phases of the child phase. Otherwise, the properties of target feature
in the child phase will not be included in the regeneration process.
When establishing the new values, the properties of target feature in the parent phase of
the target phase are examined. If the parent phase is the initial phase (which has no
previous phase), the properties as defined in the geometry will be taken as source. These
properties are then applied to everything that was selected in the determination process.
Note that the regeneration process is implicitly applied when switching from the
Geometry modes to the Calculation modes, in order to maintain staged model integrity
even if the geometry was modified. This procedure is mainly useful for quick propagation
of staged construction changes to child phases.
equilibrium stress
stress
B constitutive stress
strain
Figure 7.32 Equilibrium and constitutive stresses
At the end of each iteration, two important values of stress are calculated by PLAXIS. The
first of these, the 'equilibrium stress', is the stress calculated directly from the stiffness
matrix (e.g. point A in Figure 7.32). The second important stress, the 'constitutive stress',
is the value of stress on the material stress-strain curve at the same strain as the
equilibrium stress, i.e. point B in Figure 7.32.
The dashed line in Figure 7.32 indicates the path of the equilibrium stress. In general this
equilibrium stress path depends on the nature of the stress field and the applied loading.
For the case of a soil element obeying the Mohr-Coulomb criterion, the local error for the
particular stress point at the end of the iteration is defined:
kσ e - σ c k
Local error =
Tmax
In this equation the numerator is a norm of the difference between the equilibrium stress
tensor, σ e , and the constitutive stress tensor, σ c . This norm is defined by:
q
kσ e − σ c k =
e − σ c ) 2 + σ e − σ c 2 + (σ e − σ c ) 2 +
(σxx xx yy yy zz zz
e c 2
σxy − σxy + σyz e − σ c 2 + (σ e − σ c ) 2
yz zx zx
The denominator of the equation for the local error is the maximum value of the shear
stress as defined by the Coulomb failure criterion. In case of the Mohr-Coulomb model,
Tmax is defined as:
When the stress point is located in an interface element the following expression is used:
q
(σne − σnc ) 2 + (τ e − τ c ) 2
Local error =
ci − σnc tan ϕi
where σn and τ represent the normal and shear stresses respectively in the interface. To
quantify the local accuracy, the concept of inaccurate plastic points is used. A plastic
point is defined to be inaccurate if the local error exceeds the value of the user specified
tolerated error (Section 7.7.3).
Termination of iterations
For PLAXIS to terminate the iterations in the current load step, all of the following three
error checks must be satisfied. For further details of these error-checking procedures,
see Vermeer & van Langen (1989).
This icon represents the Output program. The main output quantities of a finite
element calculation are the displacements and the stresses. In addition, when a
finite element model involves structural elements, the structural forces in these elements
are calculated. An extensive range of facilities exists within the PLAXIS 3D Output
program to display the results of a finite element analysis. This chapter gives a
description of the features available in the program.
If the Output program is activated by running its executable file or by clicking the Output
program button in the Input program, the user has to select the model and the
appropriate calculation phase or step number for which the results are to be viewed
(Figure 8.1). More options on how to activate the Output are given in Section 8.3.1.
When a particular project is selected, the file requester displays the corresponding list of
calculation phases from which a further selection should be made. If it is desired to select
an intermediate calculation step, then a single mouse click should be given on the plus
icon (+) at the left of the desired phase. As a result, the calculation list expands a list with
all available step numbers for this phase, from which the desired step number can be
selected.
Hint: Please note that the number of the individual steps available depends on the
value assigned to Max steps saved in the Parameters tabsheet of the
Phases window.
Once an output step of a particular project has been opened, the combo box in the
toolbar will contain a list of available output steps, indicated by the step number and
corresponding phase number.
Title bar
The title bar gives information about the project name, the step number and the type of
information/results displayed.
Menu bar
The menu bar contains all output items and operations facilities of the Output program
(Section 8.2).
Toolbars
Buttons for different features in the Output program are located above and at the left side
of the plot area. A hint about the function of each tool is given as the cursor is located on
it.
Plot area
The calculation results are displayed in the Plot area. The results can be displayed in
graphical or tabular form. More information on how to handle the plot is given in Section
8.4.
Status bar
The status bar displays the locations of the cursor and the viewpoint and a hint about the
object in the model and their element numbers.
The menu bar contains drop-down menus covering the options available in the Output
program. The main results from a finite element calculation are deformations and
stresses. Hence, these two aspects form the major part of the Output menu. When
displaying a basic 3D geometry model, the menu consists of the File, View, Project,
Geometry, Mesh, Deformations, Stresses, Tools, Window and Help menus. Note that the
availability of the menus in the bar depends on the type of data that is presented on the
output form.
Use result smoothing To reduce the numerical noise resulting from the extrapolation of
the results obtained in stress points (e.g. stress, force) to nodes.
This option is available for plots and tables. Note that the option
is by default selected in plot presentation of the results.
Title To toggle the display of the title of the active plot in the caption.
Legend To toggle the display of the legend of contours or shadings.
Axes To toggle the display of the global x -, y - and z -axes in the active
plot (displayed in the lower right corner).
Local axes To toggle the display of the local 1-, 2- and 3-axes of the
structures. This option is only available when viewing structures.
Shrink To shrink the soil volumes in the model. This option is available
in the Model view.
Expand To expand the structures in the model. This option is available in
the Model view.
Explode To disassemble the components of the project geometry. This
option is available in the Model view.
Implode To assemble the components of the project geometry. This
option is available in the Model view.
Settings To set various graphical attributes, such as object and
background colours, symbol size, font size and diffuse shading.
Move cross section forward
To move the created cross section through the model enabling a
visual display of the results in the model. This option is available
in the Cross section view.
Move cross section backward
To move the created cross section through the model enabling a
visual display of the results in the model. This option is available
in the Cross section view.
Arrows To display the results as arrows.
Contour lines To display the results as contour lines.
Shadings To display the results as shadings.
Iso surfaces To display the results as iso-areas.
Node labels To display the results at nodes.
Stress point labels To display the results at stress points.
Deformation plane To display the deformed shape of cross sections, geogrids or
plates.
Distribution plane To project the results perpendicularly to the plane creating a
distribution plane for cross sections and plates.
Deformation To display the deformed shape for beams, embedded piles and
anchors.
Vertical cross section To select a user-defined vertical cross section with a distribution
of the presented quantity. The cross section is defined by
specifying two points in the corresponding window. To define the
cut line in the plot, press the left mouse button at one end of the
cross section; hold the mouse button down while moving the
mouse to the other end of the line; release the mouse button.
The cross section is presented on a new form.
Horizontal cross section To select a user-defined horizontal cross section with a
distribution of the presented quantity. The cross section must be
selected by defining the z -coordinate. The cross section is
presented on a new form.
Free cross section To select a user-defined cross section with a distribution of the
presented quantity. The cross section is defined by specifying
three points in the corresponding window. The cross section is
presented on a new form.
Line cross section To define a line on which the distribution of the presented
quantity is displayed. The line is defined by specifying two points
in the corresponding window. The results on the defined line are
presented on a new form.
Cross section curves To display a plot of the results along the line cross sections. The
values in the x-axis in the plot are the distances of the points
from the first point in the cross section.
Hint box To display a hint box with information in individual nodes or
stress points (if nodes or stress points are displayed).
Cross section points To display the points defining the cross section. These points are
displayed as greyed out in the Cross section points window.
Their location can not be modified. This option is valid only when
the Cross section view is active.
Distance measurement To measure the distance between two nodes in the model both
for the original mesh as well as the deformed mesh. This option
is valid in the Model view only when nodes and/or stress points
are displayed in the plot (Section 8.3.14).
Besides displaying the calculation results, the Output program provides tools to handle
the view and enable a better examination of the results. The buttons are grouped in the
toolbar below the menu bar and in the side tool bar. The tools and their functionality is
described in the following sections.
Copy to clipboard
Data as displayed in output forms may be exported to other programs using
the Windows clipboard function. When clicking on the Copy to clipboard button, the
Copy window appears in which selections can be made of the various plot components
that are to be included in the copy (Figure 8.3).
Print
Hardcopies of graphs and tables can be produced by sending the output to
an external printer. When the Print button is clicked or the corresponding option is
selected in the File menu, the Print window appears, in which various plot components
that are to be included in the hardcopy can be selected (Figure 8.4).
When pressing the Setup button, the standard printer setup window is presented in which
specific printer settings can be changed. When the Print button is clicked, the plot is send
to the printer. This process is fully carried out by the Windows® operating system.
Hint: When the Copy to clipboard option or the Print option is used on a plot that
shows a zoomed part of the model, only the part that is currently visible will
be exported to the clipboard or the printer.
Export
Data in output forms may be exported to files. When the Export to file button is
clicked, the Export window appears. Note that a text scaling factor can be defined.
Instead of the PLAXIS logo in the frame, it is also possible to insert a company logo. This
logo has to be provided as a bitmap and can be selected in the Print window after
clicking on the logo.
selected, a window pops up (Figure 8.8) where the factor can be defined. Note that this
option is also available in the right mouse click pop-up menu.
Hint: The default value of the Scale factor depends on the size of the model.
» The Scale factor may be used to increase or reduce the displayed (virtual)
thickness of interfaces in the Connectivity plot.
8.3.8 TABLES
The tabular form of the results given in the plot can be obtained by clicking on the
Table button or by selecting the corresponding option in the menu. Note that this
option is also available in the right mouse click pop-up menu.
Hint: The table of displacements may be used to view the global node numbers
and corresponding coordinates of individual elements.
Displaying of tables
By default, a table is presented in ascending order according to the global element
number and local node or stress point. However, a different ordering may be obtained by
clicking on the small triangle in the column header of the desired quantity on which the
ordering should be based. Another click on the same column header changes the
ordering from ascending to descending.
The options available in the right click pop-up menu are:
Select for curves To select the right clicked point in the table to be considered in
curves.
Align To align the text in the selected part of the table.
Decimal To display data in decimal representation.
Scientific To display data in scientific representation.
Hint: The values in the tables contain the most accurate information, whereas
information in plots can be influenced or be less accurate due to smoothing
or extrapolation of information from stress points to nodes.
The average results are displayed at the center of each soil element. The length of
each line represents the magnitude of the principal quantity (stress or strain) and the
direction indicates the principal direction. Positive direction is indicated by arrows.
Scaling of the results is possible.
The results are displayed in different colours at each stress point of the soil
elements. The length of each line represents the magnitude of the principal quantity
(stress or strain) and the direction indicates the principal direction. Positive direction
is indicated by arrows. Scaling of the results is possible.
The average results are displayed in different colours at the center of each soil
element. The length of each line represents the magnitude of the principal quantity
(stress or strain) and the direction indicates the principal direction. Positive direction
is indicated by arrows. Scaling of the results is possible.
The deformed shape of cross sections, geogrids or plates is displayed. Scaling of
the results is possible.
The distributions of the results in cross sections, plates, geogrids or interfaces is
displayed. Scaling of the results is possible.
The deformed shape of line structures (e.g. beam) is displayed. The relative
deformation is indicated by arrows. Scaling of the results is possible.
The distributions of the results in line structures (e.g. beam) is displayed. Scaling of
the results is possible.
The wireframe distributions of the results in line structures (e.g. beam) is displayed.
Scaling of the results is possible.
The Plastic points option shows the stress points that are in a plastic state,
displayed in a plot of the undeformed geometry (Section 9.3.7). Scaling of the
results is possible. When scaling is used, it is possible to pull the interfaces out of
the plates, however the stress points will remain at their physical locations.
The availability of the display type buttons in the toolbar can be toggled on/off by
selecting the corresponding options in the View menu.
Hint: To quickly view structures inside the 3D model, the Materials option may be
deselected from the Geometry.
Output for structures and interfaces can be obtained by clicking the Select
structures button and then double clicking the desired object in the 3D model. As a
result, a new form is opened on which the selected object appears. At the same time the
menu changes to provide the particular type of output for the selected object.
Visible model components are indicated by an open eye, whereas invisible ones are
indicated by a closed eye. By clicking on the button, the view of the components
(individual and/or groups) can be toggled from being visible to being invisible and vice
versa. A group is expanded by clicking on the + sign in front of the group. Clusters that
have been set inactive in the framework of staged construction are always invisible and
cannot be made visible.
Hint: The cluster numbers are activated by selecting the Cluster numbers option in
the Mesh menu.
The information in the Model explorer can be narrowed according to the filtering criteria
specified at the corresponding cell.
The Model explorer can be fully expanded by selecting the Expand all in menu displayed
when the Model explorer is right-clicked. The displaying menu provides the option
Collapse all that reverts the effect of the Expand all option.
The Show all option will make all the object active in the selected phase visible. The Hide
all will revert the action.The Invert selection option will toggle all visible elements invisible
and all invisible elements visible. The Deselect all button will set all elements to invisible.
On pressing the Close button the Partial geometry window is closed without further
changes.
Apart from the Model explorer, individual volume elements or entire clusters of volume
elements can be made invisible by holding down the <Ctrl> key, the <Shift> key or both
keys at the same time, respectively, while clicking on an element in the 3D model. These
elements can be visible again by clicking the corresponding check boxes in the Model
explorer.
Clicking the Hide soil button in the side bar menu enables hiding parts of the soil.
To hide soil elements, click the Hide soil button first and hold the <Ctrl> key pressed
while clicking on the soil elements. To hide soil clusters, click the Hide soil button first
and hold the <Shift> key pressed while clicking on the soil clusters.
Clicking on the Hide soil in the rectangle button enables hiding the soil
in the rectangle drawn in the model. The drawing order of the rectangle effects the
resulting hidden soil elements.
To hide only the soil elements that fall completely in the defined rectangle, first click the
Hide soil in the rectangle button. In the model, click at the point defining the upper left
corner of the rectangle, drag the mouse to the point defining the lower right corner of the
rectangle and click again.
To hide all the soil elements that are intersected by the defined rectangle, first click the
Hide soil in the rectangle button. In the model, click at the point defining the lower right
corner of the rectangle, drag the mouse to the point defining the upper left corner of the
rectangle and click again.
selection of this option, the Cross section points window pops up in which the three cross
section coordinates can be defined.
A line cross section can be defined by clicking the Line cross section button in
the side button bar or by selecting the Line cross section option in the View menu.
Upon selection of this option, the Cross section points window pops up in which the two
cross section coordinates can be defined.
After the cross section has been selected, a new form is opened in which the distribution
of a quantity is presented on the indicated cross section. At the same time, the menu
changes to allow for the selection of all other quantities that may be viewed on the
indicated cross section.
Multiple cross sections may be drawn in the same geometry. Each cross section will
appear in a different output form. To identify different cross sections, the end points of a
cross section are indicated with characters in alphabetical order. The points defining the
cross section can be viewed by selecting the Cross section points option in the Tools
menu.
In addition to the output quantities that are available for the 3D model, a cross section
allows for the display of cross section stresses, i.e. effective normal stresses σ 'N , total
normal stresses σN , vertical shear stresses τs and horizontal shear stresses τt .
Hint: It is possible to move a cross section in the direction of its normal while the
presentation of results is updated for the new location of the cross section.
• Using the <Ctrl–> and <Ctrl-+> keys will move the cross section 1/100
times the diagonal of the geometry model.
• Using the <Ctrl-Shift–> and <Ctrl-Shift-+> will move the cross section
1/1000 times the diagonal of the geometry model.
Coordinates The Original and Deformed coordinates for the first node/stress
point and the second node/stress point.
∆x The Original and Deformed x-component of the distance
between the points.
∆y The Original and Deformed y-component of the distance
between the points.
∆z The Original and Deformed z-component of the distance
between the points.
Distance The Original (v) and Deformed (v') distance between the points.
Elongation Increase of the distance between the selected points before and
after deformation without considering the rotation of the line
between the two points.
|∆u| The change in the distance between the selected points before
and after deformation.
|∆u|perpendicular
∆u
=v
P'2 '- v
- P'1
P'2 Rotation
v'=
P'1 orientation x
Elongation
- P1
v = P2 P2
P1 x
orientation
Draw scanline
When the Contour lines option is selected for the plot in cross sections or plates, a
distribution of the values can be displayed by clicking on the Draw scanline button in
the side toolbar and drawing a line on the regions of interest. Note that this option is also
available in the right mouse click pop-up menu.
Interactive ruler
The result value at a specific location in a structure or cross section can be
displayed by clicking the Interactive ruler button in the side toolbar and by moving
the cursor to the point of interest. The current value (corresponding to the point on the
cross section line), a minimum value (based on the minimum value in the distribution),
and a maximum value (based on the maximum value in the distribution) are shown along
the ruler. The Interactive ruler is available in the Structure and Cross section views.
The distribution of the results in the model is shown in the display area.
The presence of the legend, title bar, and axes in the draw area is arranged using the
options in the View menu (Section 8.2.2).
Hint: The icon in the title bar indicates the view in which the results are displayed.
A more detailed description on Views is given in Section 8.5.
8.4.1 LEGEND
The Legend is available for the display options where a variation in colour describes the
variation in the displayed result values. It is activated by selecting the corresponding
option in the View menu. When the Legend is double clicked, a window pops up, where
the scaling and the colouring can be defined (Figure 8.13). Note that this option is also
available in the right mouse click pop-up menu.
The distribution of values in the legend can be locked by clicking the Lock the
legend button. When the legend is locked, the value distribution will not change as
the <Ctrl>+<+> or <Ctrl>+<-> keys are used to move the cross section through the model.
Symbol size To modify the size of the symbols in the display for nodes,
forces, etc.
Diffuse shading To make the appearance of the 3D model even more realistic,
the Diffuse shading option may be used. Using this option,
object surfaces that have the same colour by definition (such as
soil elements with the same material data set) appear 'brighter'
or 'darker', depending on their orientation with respect to the
viewer. Object surfaces appear most bright when the normal to
the surface points in the direction of the viewer. The surfaces
become darker the more the normal deviates from this direction.
The contrast can be set to the desired magnitude using the slide
bar.
Anti aliasing To select a convenient anti aliasing method from the options
available in the drop-down menu.
Rendering method To select a convenient rendering method from the options
available in the drop-down menu.
Display Toggle the display of the Cluster borders.
The displaying colours can be arranged in the Colours tabsheet of the Settings window
(Figure 8.15).
The function of the left and the middle mouse buttons can be defined in the Manipulation
tabsheet of the Settings window (Figure 8.16).
The display of particular results can be toggled on/off in the Results tabsheet of the
Settings window (Figure 8.17).
The default visibility of Nodes, Stress points and Disabled structures can be toggled
on/off in the Default visibilities tabsheet of the Settings window (Figure 8.18).
In the Output program the results are displayed in different views. The view type is
indicated by the corresponding icon in the title bar (Figure 8.12). The available views are:
Step 2: Select the phases for which results will be included in the report (Figure 8.20).
Step 3: Select general information sets to be included in the report. Note that the
selection can be saved as a new set besides All and None sets (Figure 8.21).
Step 4: Select model view sets to be displayed in the report. Note that the selection can
be saved as a new set besides All and None sets (Figure 8.22).
Step 5: Select structure view sets to be displayed in the report are selected. Note that the
selection can be saved as a new set besides All and None sets (Figure 8.23).
Step 6: Select saved views to be included in the report (Figure 8.24). For more details on
saved views, see Section 8.3.4.
Step 8: A summary of the number of rows and the number of figures in the report is given
(Figure 8.26). The report is created as the Export button is clicked. A progress bar
appears displaying the number of the remaining rows and images.
The Create animation option is available in the View menu. If the option is
selected, the Create animation window appears (see Figure 8.28). The phases and
calculation steps to be included in the animations can be selected. Note that the
animation will include only available output steps. This depends on the Max. steps stored
parameter as defined for each phase in the Phases window. After selecting the phase(s),
click OK to start the process. The progress of this process is indicated in a separate
window.
If a large number of steps is to be included in the animation, the process may take some
minutes after which the animation is presented. The result is stored in an animation file
(*.AVI) in the project data directory.
A Connectivity plot is a plot of the mesh in which the element connections are clearly
visualised. It is the result of the meshing process. It is available only in the representation
of spatial variation of the results. This plot is particularly of interest when interface
elements are included in the mesh. Interface elements are composed of pairs of nodes in
which the nodes in a pair have the same coordinates. In the Connectivity plot however,
the nodes in a pair are drawn with a certain distance in between so that it is made clear
how nodes are connected to adjacent elements. This option is available from the Mesh
menu.
In the Connectivity plot it can, for example, be seen that when an interface is present
between two soil elements, that the soil elements do not have common nodes and that
the connection is formed by the interface. In a situation where interfaces are placed
along both sides of a plate (Positive interface and Negative interface), the plate and the
adjacent soil elements do not have nodes in common. The connection between the plate
and the soil is formed by the interface. An example of Connectivity plot is given in Figure
9.1.
9.2 DEFORMATIONS
The Deformations menu contains various options to visualise the displacements and
strains in the finite element model. By default, the displayed quantities are scaled
automatically by a factor (1, 2 or 5) ·10n to give a diagram that may be read conveniently.
The scale factor may be changed by clicking the Scale factor button in the toolbar
or by selecting the Scale option from the View menu. The scale factor for strains
refers to a reference value of strain that is drawn as a certain percentage of the geometry
dimensions. To be able to compare plots of different calculation phases or different
projects, the scale factors in the different plots must be made equal.
Elevation
The height of the soil volume is displayed when the Elevation option is selected in the
Deformations menu.
way, the settlement of the last construction layer will be limited and the largest settlement
will most likely occur in the middle of the embankment, as expected. When plotting the
settlements in a vertical cross section through the embankment, the results are
somewhat discontinuous, but the overall settlement profile is more realistic than without
choosing this option (Figure 9.2b). The more construction layers are used, the smoother
the settlement profile is (Figure 9.2c).
a. b. c.
Figure 9.2 Settlement profile of an embankment on a stiff foundation layer: a. Phase displacement
results; b. Sum phase displacements results (5 construction layers); c. Sum phase
displacements results (10 construction layers)
Note that, if displacements are set to zero in intermediate phases, the Sum phase
displacements results may (at some point) be larger than the total displacements.
9.2.6 VELOCITIES
The option Velocities contains the different components of the velocities at the end of the
current calculation step, displayed on a plot of the geometry. This option may be selected
from the Deformations menu. A further selection can be made among the velocity
vectors, |v |, and the individual velocity components, vx , vy and vz . The velocities may be
presented as Arrows, Contour lines, Iso surfaces or Shadings by clicking the appropriate
button in the toolbar (Section 8.3.10).
9.2.7 ACCELERATIONS
The option Accelerations contains the different components of the accelerations at the
end of the current calculation step, displayed on a plot of the geometry. This option may
be selected from the Deformations menu. A further selection can be made among the
acceleration vectors, |a|, and the individual acceleration components, ax , ay and az . The
accelerations may be presented as Arrows, Contour lines, Iso surfaces or Shadings by
clicking the appropriate button in the toolbar (Section 8.3.10).
can be made among the principal strain directions, the individual principal strain
components 1 , 2 , 3 , (1 + 3 )/2, (1 - 3 )/2, the angle, the volumetric strain v , the
deviatoric strain γs and the void ratio e.
• Note that the principal strain components are arranged in algebraic order:
1 > 2 > 3
v = xx + yy + zz + xx yy + xx zz + yy zz + xx yy zz
• The deviatoric strain is calculated as:
s
2 v 2 v 2 v 2 1 2 2 + γ2
γs = xx − + yy − + zz − + γxy + γyz zx
3 3 3 3 2
• The void ratio is calculated as:
e = e0 + (1 + e0 )v
9.3 STRESSES
Various options are available to visualize the stress state in the finite element model.
Hint: By default, the stresses developed in non-porous materials are not displayed
in the plot. To display them select the Show stress for nonporous material
option in the Results tabsheet of the Settings window (Section 8.4.2).
σzz
z
σzy
σzx σyz
σxz
σyy
y σxy σyx
x σxx
Hence, σ '1 is the largest compressive (or smallest tension) principal stress and σ '3 is the
smallest compressive (or largest tension) principal stress.
The Mobilised shear strength τmob is the maximum value of shear stress (i.e. the radius
of the Mohr stress circle or half the maximum principal stress difference).
The Relative shear stress τrel gives an indication of the proximity of the stress point to the
failure envelope, and is defined as:
τmob
τrel =
τmax
where τmax is the maximum value of shear stress for the case where the Mohr's circle is
expanded to touch the Coulomb failure envelope while keeping the center of Mohr's circle
constant.
σ1 ' + σ3 '
τmax = − sin ϕ + c cos ϕ
2
Hint: Particularly when the soil strength has been defined by means of effective
strength parameters (Undrained (A) it is useful to plot the mobilised shear
strength τmob in a vertical cross section and to check this against a known
shear strength profile.
When using the Hoek-Brown model to describe the behaviour of a rock section, the
definition of the maximum shear stress τmax is slightly modified. Starting from the
Hoek-Brown failure criterion:
1 '3
f (σ '3 ) = σci mb −σ
τmax = f (σ '3 ) where σci + s
a
(9.2)
2
The relative shear stress is correspondingly defined by:
σ1 ≤ σ2 ≤ σ3
Hence, σ1 is the largest compressive (or the smallest tension) principal stress and σ3 is
the smallest compressive (or the largest tension) principal stress.
σ = σ ' + pactive
Active pore pressures are defined as the effective saturation, Seff , times the pore water
pressures, pw .
pactive = Seff · pw
Pore water pressure differs from active pore pressure when the degree of saturation is
less than unity, which is usually the case above the phreatic level. Below the phreatic
level pactive and pw are generally equal.
As an alternative to the pore water pressure, pw , the groundwater head, h, can be viewed:
pw
h=z−
γw
where z is the vertical coordinate and γw is the unit weight of water. In the pore water
pressure a further distinction is made between steady state pore pressure, psteady , and
excess pore pressure, pexcess .
pw = psteady + pexcess
where steady-state pore pressure is the steady-state or long term part of pore pressure,
which is supposed to be input data for a deformation analysis.
Excess pore pressure is the result from undrained behaviour (Undrained (A) or (B), or
low-permeable materials), and is affected by stress changes due to loading or unloading,
a (sudden) change in hydraulic conditions and consolidation.
Below an overview is given of the various pore pressure related quantities that are
available in the Stresses menu of the PLAXIS Output program. Most quantities can be
presented as Contour lines, Shadings or Iso surfaces. Although pore pressures do not
have principal directions, the Principal stresses presentation can be useful to view pore
pressures inside the model. In that case the colour of the lines represents the magnitude
of the pore pressure and the directions coincide with the x-, y- and z-axis.
Groundwater head: The groundwater head is an alternative quantity for the pore water
pressure, equivalent to the top level of a free water column.
Active pore pressure: The part of total stress in the pores of the soil. In saturated soil,
active pore pressure is equal to the pore water pressure. In unsaturated soil, active pore
pressure is the effective saturation times the pore water pressure.
Pore water pressure: The water pressure in the pores of the soil consisting of
steady-state pore pressure and excess pore pressure. In most cases, pore water
pressure is negative (pressure). However, due to capillary action or undrained unloading,
pore water pressure may become positive (suction).
Steady-state pore pressures: Long-term part of pore water pressure which is
supposed to be input data for a deformation analysis. Steady-state pore pressure is the
result of a direct pore pressure generation based on phreatic levels and cluster-related
pore pressure definition, or the result of a steady-state groundwater flow calculation. In
the case of a fully coupled flow-deformation analysis, steady-state pore pressure is
calculated from a steady-state groundwater flow calculation using the hydraulic boundary
conditions at the end of the calculation phase. When the option Ignore suction has been
used, the steady-state pore pressure only includes negative values.
Excess pore pressure: Instationary part of pore water pressure that is a result of
undrained behaviour (Undrained (A) or (B), or low-permeable materials), affected by
stress changes due to loading or unloading, a (sudden) change in hydraulic conditions
and consolidation. In the case of a fully coupled flow-deformation analysis, excess pore
pressure is the difference between the calculated pore water pressure and the
steady-state pore pressure. The latter is calculated from a steady-state groundwater flow
calculation using the hydraulic boundary conditions at the end of the calculation phase.
Excess pore pressure may show positive values (suction) as a result of undrained
unloading, even when the option Ignore suction has been selected.
Extreme excess pore pressure: Maximum and minimum values of excess pore
pressure over the whole calculation phase up to the selected calculation step.
Change in pore pressure per phase: The change in active pore pressure from the
beginning of a phase up to the selected calculation step.
Suction: Any positive value (tension) in the pore water pressure. Suction can be a result
of capillary action (included in the steady-state pore pressure if the option Ignore suction
has NOT been used) or undrained unloading (included in the excess pore pressure).
Effective suction: Any positive value (tension) in the active pore pressure. When
multiplied by the tangent of the friction angle, the effective suction gives a kind of’'artificial
cohesion' in the soil.
The failure points are particularly useful to check whether the size of the mesh is
sufficient. If the zone of plasticity reaches a mesh boundary (excluding the symmetry
plane in a symmetric model) then this suggests that the size of the mesh may be too
small. In this case the calculation should be repeated with a larger model.
When Plastic points is selected in the Stresses menu the Plastic points dialog is shown
(Figure 9.4). Here the user can select which types of plastic points are displayed. When
the Stress points option is selected, all other stress points are indicated by a purple
diamond shape (). For details of the use of advanced soil models, the user is referred to
the Material Models Manual.
By default both accurate and inaccurate plastic points are displayed in the model. Only
the inaccurate plastic points are displayed as the corresponding check box is selected in
the Plastic points window. Inaccurate plastic points are points where the local error is
larger then the tolerated error (Section 7.10.8).
Hint: The Plastic point history option in the Stresses menu enables displaying in
the model all the points that have been plastic in any step in the calculation
history (depending on the specified criteria, Failure, Tension cut-off, etc.)
generated up to the current calculation phase.
9.3.10 WELLS
When Wells is selected in the Stresses menu a table appears displaying the wells
available in the model, the nodes representing the well and their location, the discharge
of the well and the defined minimum groundwater head.
9.3.11 DRAINS
When Drains is selected in the Stresses menu a table appears displaying the drains
available in the model, the nodes representing the drain and their location, the total
discharge and the defined groundwater head of the drain.
By default, structures (i.e. anchors, beams, embedded piles, geogrids and plates) and
interfaces are displayed in the geometry. Otherwise, these objects may be displayed by
selecting the Structures or Interfaces option from the Geometry menu. To quickly view
structures inside the 3D model, the Materials option may be deselected from the
Geometry menu.
Output for structures and interfaces can be obtained by clicking the Select
structures button and then double clicking the desired object in the 3D model. As a
result, a new form is opened on which the selected object appears. At the same time the
menu changes to provide the particular type of output for the selected object.
All objects of the same type with the same local coordinate system are automatically
selected. When multiple objects or multiple groups of objects of the same type need to
be selected, the <Shift> key should be used while selecting the objects. The last object to
be included in the plot should then be double clicked. When all objects of the same type
are to be selected, select one of the objects while pressing <Ctrl-A> simultaneously. If it
is desired to select one or more individual elements from a group, the <Ctrl> key should
be used while selecting the desired element.
Another option of selecting structural elements in the output is by clicking the Drag
a window to select structures button and drawing a rectangle in the model. As a
results, the structures in the rectangle will be selected.
Figure 9.5 Sign convention for axial forces in beams and plates
The Structure axes option from the View menu may be used to display the beam's local
system of axes (1,2,3). The first direction is always the axial direction. The second and
third axis are always perpendicular to the beam axis.
The Axial force N is the force in the first (axial) direction (Figure 9.6b). An axial force is
positive when it generates tensile stresses, as indicated in Figure 9.5.
3
N ε Q12 Q13
13 2 13 13
a. Local Axes b. Axial force N c. Shear force Q12 d. Shear force Q13
Figure 9.6 Positive axial force and shear forces in beams
1
1
3
3
2
2
I3 M3 κ3 I2 M2 κ2
κ3
a. Bending moment M3 b. Bending moment M2
Figure 9.7 Positive bending moments in beams
The Shear force Q12 is the shear force over the second beam axis (Figure 9.6c), whereas
the Shear force Q13 is the shear force over the third beam axis (Figure 9.6d).
The Bending moment M3 is the bending moment due to bending around the third axis
(Figure 9.7a), whereas the Bending moment M2 is the bending moment due to bending
around the second axis (Figure 9.7b).
2 3
ε2 ε2 N1 E1 ε1 N2 E2 ε2
a. Local plate directions b. Axial force N1 c. Axial force N2
Figure 9.8 Positive axial forces in plates and geogrids
Hint: The Axes option from the View menu may be used to display the pile's local
system of axes (1,2,3). The first direction is always the axial direction. The
second direction is perpendicular to the pile axis with a component in the
global x-direction and the third direction is also perpendicular to the pile with
a component in the global z-direction.
The pile-soil interaction forces are obtained from the special interface that is
automatically applied between the embedded beam elements and the surrounding soil
volume elements. The Skin force Tskin , expressed in the unit of force per unit of pile
length, is the force related to the relative displacement in the pile's first direction (axial
direction). This force is limited by the skin resistance as defined in the embedded pile
material data set (Section 6.7).
The interaction force T2 relates to the relative displacement perpendicular to the pile in
the pile's second direction whereas the interaction force T3 relates to the relative
displacement perpendicular to the pile in the pile's third direction. These quantities are
expressed in the unit of force per unit of pile length. Note that T2 and T3 are not limited.
In fact, when these forces become very large, plasticity will occur outside the elastic zone
in the surrounding soil volume elements.
The maximum shear stress Tmax is the limit defined for the material dataset. The relative
shear stress Trel gives an indication of the proximity of the stress point to the failure
envelope.
The pile foot force Ffoot , expressed in the unit of force, is obtained from the relative
displacement in the axial pile direction between the foot or tip of the pile and the
surrounding soil. The foot force is shown in the plot of the Axial force N. The foot force is
limited by the base resistance as defined in the embedded pile material data set (Section
6.7).
10 CURVES
The location in the model where the variation of results through calculation steps is to be
analyzed is specified by selecting nodes or stress points in the model. The selection of
points should be done preferably before but may also be done after calculating the
project.
In order to specify points to be considered in curves, the Select points for curves
option should be selected. In the Input program, this button is available in the side
toolbar of the Mesh and Staged construction modes. Selecting this option will open the
Output program displaying the Connectivity plot and the Select points window.
Nodes and stress points can be selected in the Output program either by clicking the
Select points for curves button in the side toolbar or by selecting the corresponding option
in the Tools menu. More information on selecting procedure is given in Section 10.1.1.
It is important to consider the differences in selecting the points before or after starting
the calculation process. A more detailed description is given in Section 10.1.2 and
Section 10.1.3.
Hint: When the Select points for curves option is selected but the Select points
window is closed, it can be displayed by selecting the Mesh point selection
option in the Tools menu.
If the finite element mesh is regenerated (after being refined or modified), the position of
nodes and stress points will change. As a result, previously selected nodes and stress
points may appear in completely different positions. Therefore nodes and stress points
should be reselected after regeneration of the mesh.
Hint: Pre-calculation points provide detailed information related to stress and strain
at those points. However information about structural forces and state
parameters is not provided.
The information available for selected points (nodes or stress point) depends on the view
in which they have been selected in the Output program.
The points selected in the Model view, can be used to generate curves related to
displacements, stresses, strains and state parameters in soil elements. The Model
view is the default view in the Output program.
The points selected in the Structure view, can be used to generate curves related
to resulting structural forces. The points should be selected after selecting the
structure first (Section 8.3.11). The Structure view is displayed when structures are
selected and double clicked.
Hint: The type of the active view is indicated by the corresponding icon under the
plot.
To generate curves, the Curves manager option should be selected from the Tools
menu or the corresponding button in the toolbar should be clicked. As a result, the
Curves manager window appears with three tabsheets named Charts, Curve points and
Select points.
The Charts tabsheet contains the saved charts that were previously generated for the
current project. The Curve points tabsheet gives an overview of the nodes and stress
points that were selected for the generation of curves, with an indication of their
coordinates. The list includes the points selected before the calculation (pre-calc) as well
as the points selected after the calculation (post-calc) (Figure 10.2). For points that are
part of a structure further information is given in the list about the type of structure and
the corresponding structure element number.The Select points window is described in
Section 10.1.1
As a next step to generate curves, the New button should be pressed while the Charts
tabsheet is active. As a result, the Curve generation window appears, as presented in
Figure 10.3.
Two similar groups with various items are shown, one for the x -axis and one for the
y -axis of the curve. The x -axis corresponds to the horizontal axis and the y -axis
corresponds to the vertical axis. For each axis, a combination of selections should be
made to define which quantity is plotted on that axis. First, for each axis a selection
should be made whether the data to be shown is related to the general project (Project)
or a particular selected node or stress point. The tree in the Curve generation window will
then show all quantities which are available depending for this type of data. The tree can
be expanded by clicking the + sign in front of a group. The Invert sign option may be
selected to multiply all values of the x -quantity or the y -quantity by -1. When both
quantities have been defined and the OK button is pressed, the curve is generated and
presented in a chart window.
The combination of the step-dependent values of the x -quantity and the y -quantity forms
the points of the curve to be plotted. The number of curve points corresponds to the
available calculation step numbers plus one. The first curve point (corresponding to step
0) is numbered as 1.
Hint: When curves are generated from points selected after the calculation, only
information of saved steps can be considered. The number of the saved
steps for each calculation phase is defined by the Maximum number of steps
stored option in the Parameters tabsheet of the Phases window (Section
7.7.3).
» All the calculation results are available for the pre-selected points.
load applied in current phase - Total load applied in previous phase). Also other types of
curves can be generated.
The selection of Displacement must be completed with the selection of a node in the
drop-down menu and the selection of a displacement component in the Deformations
subtree. The type of displacement can be either the length of the displacement vector
(|u |) or one of the individual displacement components (ux , uy or uz ). The displacements
are expressed in the unit of length, as specified in the Project properties window of the
Input program.
To define a multiplier on the y -axis, first the Project option should be selected as the
activation of a load system is not related to a particular point in the geometry. The
selection must be completed with the selection of the desired load system, represented
by the corresponding multiplier in the Multiplier subtree. Note that the 'load' is not
expressed in units of stress or force but in a multiplier value without unit. To obtain the
actual load, the presented value should be multiplied by the input load as specified by
means of staged construction.
Another quantity that can be presented in a curve is the Pore pressure. This quantity is
available for selected nodes as well as stress points. In the Pore pressures subtree of the
Stresses tree pactive , psteady or pexcess can be selected. Pore pressures are expressed in
the unit of stress.
When non-zero prescribed displacements are activated in a calculation, the reaction
forces against the prescribed displacements in the x - and y -direction are calculated and
stored as output parameters. These force components can also be used in the
load-displacement curves by selecting the option Project and then selecting one of the
forces in the Forces subtree.
The Amplification tabsheet enables obtaining the plot which shows the ratio of the
acceleration response of any point (Top) to the acceleration response of another point
(Bottom) which is preferably the point where input load is applied (Figure 10.6). This
gives the magnification of the response at one point with respect to given excitation.
For the curves created in the PSA tabsheet of the Curve generation window, the
Displacement response factor can be selected (Figure 10.8), to display the variation of
displacement with frequency.
Hint: The Settings window is displayed by right clicking the chart and selecting the
corresponding option in the appearing menu or by selecting the option in the
Format menu.
Figure 10.8 Displacement response factor option available for PSA curves
Once a curve has been generated, a new chart window is opened in which the generated
curve is presented. The quantities used to generate the curve are plotted along the x -
and y -axis. By default, a legend is presented at the right hand side of the chart. For all
curves in a chart, the legend contains the Curve title, which is automatically generated
with the curve. An example of the curves in Output program is given in Figure 10.9.
File menu
The File menu is basically the same with the one available in the Output program. For a
more detailed description see Section 8.2.1.
Edit menu
Note that Edit menu is only available when the curves are displayed. The options
available can be used to include curves in the current chart. These options are:
Copy To export the chart to other programs using the Windows
clipboard function. This feature is described in detail in Section
8.3.2.
Add curve from current project
To add a new curve to the active chart from the current project.
Add curve from another project
To add a new curve to the active chart from another project.
Add curve from clipboard
To add a new curve to the active chart from clipboard.
Hint: The added curves are redefined using the data from either the current
project, another project or clipboard. It is not possible to mount a generated
curve to the current chart.
» It is possible to add a curve to the active chart using the Add curve option in
the corresponding option in the right mouse click pop-up menu.
View menu
The display of the results in the window is arranged using the options available in the
View menu. These options are:
Reset view To reset a zoomed view.
Hint: For a more detailed view of particular regions in curves, press the left mouse
button at a corner of the zoom area; hold the mouse button down and move
the mouse to the opposite corner of the zoom area; then release the button.
The program will zoom into the selected area. The zoom option may be used
repetitively.
The zoomed view can be reset by clicking the corresponding button in the
toolbar as well.
Format menu
The Format menu contains the Settings option, selecting which displays the
corresponding window where the layout of the chart and curves can be modified.
The layout and presentation of charts can be modified by clicking the Settings
button available in the toolbar or by selecting the corresponding option in the
Format menu. Alternatively, the Settings option can be selected from the Format menu of
the right mouse button menu. As a result, the Settings window will appear. Distinction is
made between the chart settings displayed on the first tabsheet and the curve settings
displayed on a separate tabsheet for each curve. The options available in the Chart
tabsheet can be used to customize the frame and axes of the chart (Section 10.4.1). The
options available in the tabsheets of the curves can be used to customize the plot
(Section 10.4.2).
If the correct settings are defined, the OK button may be pressed to activate the settings
and to close the window. Alternatively, the Apply button may be pressed to activate the
settings, keeping the Settings window active. The changes to the settings can be ignored
by pressing the Cancel button.
Titles By default, a title is given to the x -axis and the y -axis, based on
the quantity that is selected for the curve generation. However,
this title may be changed in the Title edit boxes of the
corresponding axis group. In addition, a title may be given to the
full chart, which can be entered in the Chart name edit box. This
title should not be confused with the Curve title as described in
above.
Scaling of x- and y-axis By default, the range of values indicated on the x - and y -axis is
scaled automatically, but the user can select the Manual option
and enter the desired range in the Minimum and Maximum edit
boxes. As a result, data outside this range will not appear in the
The Settings window contains for each of the curves in the current chart a tabsheet with
the same options (Figure 10.11).
Title A default title is given to any curve during its generation. This
title may be changed in the Curve title edit box. When a legend
is presented for the active chart in the main window, the Curve
title appears in the legend.
Show curve When multiple curves are present within one chart, it may be
useful to hide temporarily one or more curves to focus attention
on the others. The Show curve option may be deselected for this
purpose.
Phases The Phases button may be used to select for which calculation
phases the curve has to be generated. This option is useful
when not all calculation phases should be included in the curve.
Fitting To draw a smooth curve, the user can select the Fitting item.
When doing so, the type of fitting can be selected from the Type
combo box. The Spline fitting generally gives the most
satisfactory results, but, as an alternative, a curve can be fitted to
a polynomial using the least squares method.
Line and marker presentation
Various options are available to customise the appearance of the
curve lines and markers.
Arrow buttons The arrow buttons can be used to change the order of the curves
in the legend.
Regenerate The Regenerate button may be used to regenerate a previously
generated curve to comply with new data (Section 10.5).
Add curve The Add curve button may be used to add new curves to the
current chart (Section 10.6).
Delete When multiple curves are present within one chart, the Delete
button may be used to erase a curve.
If, for any reason, a calculation process is repeated or extended with new calculation
phases, it is generally desirable to update existing curves to comply with the new data.
This can be done by means of the Regenerate facility. This facility is available in the
Settings tabsheet (Figure 10.10), which can be opened by selecting the Settings option
from the Format menu. When clicking on the Regenerate button, the Curve generation
window appears, showing the existing setting for x - and y -axis. Pressing the OK button
is sufficient to regenerate the curve to include the new data. Another OK closes the
Settings window and displays the newly generated curve.
When multiple curves are used in one chart, the Regenerate facility should be used for
each curve individually. The Regenerate facility may also be used to change the quantity
that is plotted on the x - or y -axis.
It is often useful to compare similar curves for different points in a geometry, or even in
different geometries or projects. Therefore PLAXIS allows for the generation of more
than one curve in the same chart. Once a single curve has been generated, the Add
curve options in the Edit menu can be used to generate a new curve in the current chart.
As an alternative, the Add curve option from the Settings window or from the right mouse
button menu can be used. Distinction is made between a new curve from the current
project, a new curve from another project or curves available on the clipboard.
The Add curve procedure is similar to the generation of a new curve (Section 10.2).
However, when it comes to the actual generation of the curve, the program imposes
some restrictions on the selection of data to be presented on the x - and the y -axis. This
is to ensure that the new data are consistent with the data of the existing curve.
When the Add curve option is used, the current chart is modified. In order to preserve the
current chart, a copy of it can be created by selecting it first in the list and then by clicking
the Copy button in the Curves manager window.
11 REFERENCES
[1] Bathe, K.J. (1982). Finite element analysis in engineering analysis. Prentice-Hall,
New Jersey.
[2] Benz, T., Schwab, R., Vermeer, P.A., Kauther, R.A. (2007). A Hoek-Brown criterion
with intrinsic material strength factorization. Int. J. of Rock Mechanics and Mining
Sci., 45(2), 210–222.
[3] Brinkgreve, R.B.J., Bakker, H.L. (1991). Non-linear finite element analysis of safety
factors. In Proc. 7th Int. Conf. on Comp. Methods and Advances in Geomechanics.
Cairns, Australia, 1117–1122.
[4] Brinkgreve, R.B.J., Kappert, M.H., Bonnier, P.G. (2007). Hysteretic damping in
small-strain stiffness model.
[5] CUR (2004). Geotechnical exchange format for cpt-data. Technical report, CUR.
[6] Goodman, R.E., Taylor, R.L., Brekke, T.L. (1968). A model for mechanics of jointed
rock. Journal of the Soil Mechanics and Foundations Division, 94, 19–43.
[7] Schank, O., Gärtner, K. (2006). On fast factorization pivoting methods for symmetric
indefinite systems. Electronic Transactions on Numerical Analysis, 23, 158–179.
[8] Schank, O., Wächter, A., Hagemann, M. (2007). Matching−based preprocessing
algorithms to the solution of saddle−point problems in large−scale nonconvex
interior−point optimization. Computational Optimization and Applications, 36 (2-3),
321–341.
[9] van Langen, H. (1991). Numerical analysis of soil structure interaction. Ph.D.
thesis, Delft University of Technology, The Netherlands.
[10] van Langen, H., Vermeer, P.A. (1991). Interface elements for singular plasticity
points. Int. J. Num. Analyt. Meth. in Geomech., 15, 301–315.
[11] Vermeer, P.A., van Langen, H. (1989). Soil collapse computations with finite
elements. In Ingenieur-Archive 59. 221–236.
[12] Zienkiewicz, O.C. (1977). The Finite Element Method. McGraw-Hill, London.
INDEX
A output · 259
Cluster · 268
Accelerations · 284
Command
g · 284
line · 23
Anchor
prestressing · 239 Connections · 235
properties · 178 Connectivity plot · 258, 261
Arc-length control · 195 Copy object
Array array · 60
copy objects · 60 Copy to clipboard
Assigning data sets · 179 Input · 24
Automatic Output · 262
error checks · 250 Coulomb point · 191
step size · 207 Create animation · 255
stepsize · 220 Cross section
Avoid predeforming · 282 free · 260
horizontal · 260
B line · 260
Beams output · 270
output · 294 vertical · 260
Boundary conditions Curve
adjustments during calculation · 249 generation · 301
groundwater head · 229 regeneration · 312
submerged boundaries · 198 settings · 310
viscous · 240 Curves manager · 263
C D
Calculation Damping
abort · 243 Rayleigh · 115
Advanced mode · 194 Data sets
automatic step size · 220 assigning · 179
phase · 249 Default iterative procedure · 216
plastic · 192 deformation
staged construction · 201 control parameters · 212
type Deformations · 281
Consolidation · 193 Displacement
Dynamic · 197 incremental · 283
Gravity · 191 phase · 282
K0 procedure · 189 reset · 213
Plastic · 192 total · 282
Plastic nil-step · 192 Distributed load
Safety · 195 change · 235
Calculation steps Distributed prescribed displacement
max number of steps stored · 216 change · 238
max steps saved · 300 Drain
CamClay · 112 line · 94
CavitationCutoff · 214 surface · 95
Clipboard Drained behaviour · 113
APPENDIX A - THE FORMAT REQUIRED FOR THE FILES TO IMPORT CPT DATA
The file must contain a small header specifying the contained data. If it is not available,
the field data cannot be recognized and the importing of the file will be aborted. The
header must be the first content in the file in order to prevent that arbitrary data of the file
is misinterpreted. The header must have the following format:
X[m] # optional comment
Y # optional comment
Z # optional comment
d[m] q f x x u # optional comment
# optional comment
X , Y and Z denote a number specifying the location and elevation of the starting point of
the CPT. Missing coordinates are set automatically to zero. Note that PLAXIS does not
use this information for modelling purposes.
Following the coordinate specification, specifiers for the column content must be
specified. The following specifiers are supported:
d penetration depth
q tip resistance
f sleeve friction
u pore pressure
x columns of data to be skipped
Except for the 'x' character, column specifiers can occur only once. If a column specifier
occurs more than once, the double specifier and all successive specifiers are skipped.
Units can be introduces inside square brackets next to the column specifier. Comments
can be added behind a '#'.
The first number containing a decimal separator determines whether European or
English settings are considered in a file.
• 1. or .23 or 1.23 or 1.23E-5 renders English format settings
• 1, or ,23 or 1,23 or 1,23E-5 renders European format settings
• 1E-5 or 125 assumes English format settings, the next number containing a decimal
separator decides on the format settings.
• Thousand separators are not recognized.
The user can assign the extension *.cpt to a file. In this case the program will recognize
the content as CPT data. No header needs to be specified because the file extension
indicates its content. If no header is specified or its data is incomplete, the following
values are assumed:
• Missing coordinates are set to 0.0
The program accepts only [m] as length unit. The conversion of the data units to the
project units shown in Table A.1 is possible for pressure.
Table A.1 Convertable units
Data unit Project unit
2
1 Pa 1 N/m
1 kPa 1000 N/m2 1 kN/m2
1 MPa 1 · 106 N/m2 1000kN/m2
2 6 2
1 N/mm 1 · 10 N/m 1000 kN/m2
1 kN/m2 1000 N/m2 1 kN/m2
† For hydrostatic cases, the situations different from the Head are considered. If the phreatic level in hydrostatic
conditions is equal to Head, no extra water level will be generated.
‡ The number of generated water levels for cases where more than one soil layer is available can be derived
accordingly.
As of 2013, PLAXIS 3D can handle unsaturated soils. The unsaturated soil behaviour
and suction are taken into account using Bishop' s stress, whereas conventional
calculations (without suction) are performed using Terzaghi' s effective stress. The way to
distinguish between conventional calculations (without suction) and calculations in which
unsaturated soil behaviour is taken into account (with suction) is by means of the
calculation option Ignore suction. This option is available for all calculation types, i.e.
Plastic calculations, Safety analysis, Dynamic analysis, Consolidation analysis and Fully
coupled flow-deformation analysis.
Before defining the possibilities and limitations of both options, a definition is given of the
various stresses and related quantities involved. Note that compression is negative and
tension (or suction) are positive.
Total stresses are divided into effective stresses and pore pressures. Pore pressures are
composed of steady-state pore pressures and excess pore pressures, including suction.
Since decades, PLAXIS uses the term 'active pore pressure' to denote the contribution of
the pore stresses in the total stresses. In the context of unsaturated soil behaviour, active
pore pressure is defined as:
Active pore pressure: pactive = Seff (psteady + pexcess ) = Seff pw
Pore water pressure: pwater = psteady + pexcess
Note that for fully saturated soil the 'active pore pressure' is equal to the 'pore water
pressure', but there is a clear difference in the unsaturated zone.
Steady-state pore pressure: psteady
• Input for a deformation analysis
• Direct generation based on phreatic levels and cluster-related pore pressure
distribution
• Or, calculated from steady-state groundwater flow calculation
Excess pore pressure: pexcess
• Result from undrained behaviour (Undrained (A) or (B)) and Kw /n from Table C.1
• Affected by loading, a (sudden) change in hydraulic conditions and consolidation.
• In the case of a fully coupled flow-deformation analysis, psteady is calculated from a
preliminary steady-state groundwater flow calculation using the hydraulic boundary
conditions at the end of the calculation phase. This enables the calculation and
output of pexcess for a fully coupled analysis in all steps. pexcess = pwater − psteady .
Effective degree of saturation: Seff = (S − Sres )/(Ssat − Sres )
Degree of saturation: S
Residual degree of saturation: Sres
Saturated degree of saturation: Ssat (usually 1.0)
When the calculation option Ignore suction is selected, soil below the phreatic level is
considered to be fully saturated (S = 1), whereas soil above the phreatic level is
considered to be ideally unsaturated (S = 0). The phreatic level itself is defined directly
by the user or calculated as a result of a groundwater flow calculation or a fully coupled
flow-deformation analysis. Positive steady-state pore stresses will be set to zero.
However, any excess pore pressure above and below the phreatic level, both positive and
negative, will be taken into account. This requires Seff to be set to 1, regardless of S , Sres
and Ssat .
When selecting the calculation option Ignore suction, the following conditions apply:
• Any positive value of psteady is cut at 0.
• On or below the phreatic level (psteady ≤ 0)1 : S = 1; Sres = 0; Ssat = 1 > Seff = 1
• In drained and undrained (A,B) materials2 : Seff = 1 (over-rules previous value)
1. For a fully coupled flow-deformation analysis the position of the phreatic level is
based on pwater instead of psteady .
2. In contrast to other calculation types, Seff (representing χ) is NOT overruled in a fully
coupled flow-deformation analysis. Hence, Seff and pactive are 0 above the phreatic
level.
With this, the stresses and related quantities as defined previously fully apply. In
principle, all calculation types (Plastic, Consolidation, Fully coupled flow-deformation,
Dynamic, Safety) are available, but some specific features or combinations are not
available. Limitations are described below (see Section C.6).
Note that it is still possible to have suction (pore water tension) as a result of unloading in
Undrained (A,B) materials. Hence, the option Ignore suction will NOT ignore suction as a
result of excess pore tension.
When the calculation option Ignore suction is NOT selected, suction is allowed and
included in the active pore pressure and pore water pressure. The saturation of the soil
depends on the soil-water retention curve as defined in the corresponding material data
sets.
With this, the stresses and related quantities as defined previously fully apply. In
principle, all calculation types (Plastic, Consolidation, Fully coupled flow-deformation,
Dynamic, Safety) are available, but some specific features or combinations are not
available. Limitations are described below (see Section C.6).
Positive pore stresses may also occur as a result of unloading in undrained materials.
However, there is a limit to the amount of positive pore stresses (pore water tension) that
can occur, which is the cavitation stress. This limit can be taken into account in PLAXIS
by selecting the Cavitation cut-off option and setting the cavitation stress. As a result, the
pore water pressure, pwater , cannot be higher than the cavitation stress, pcavitation . As
soon as the cavitation stress is reached, the excess pore pressure is manipulated:
At cavitation pexcess = pcavitation − psteady
Although a proper definition of the unit weight of soil would involve the degree of
saturation, PLAXIS applies a unit weight based on the current position of the phreatic
level in all types of calculations.
There are limitations in the use of specific features in the various calculations types in
PLAXIS depending on whether suction is ignored or allowed. An overview of the
limitations is given below.
C.6.3 K0-PROCEDURE
• Only possible as first (initial) calculation phase.
• Does not consider external loads; only material weight.
• Recommended for situations involving a horizontal ground surface and horizontal
soil layering without excavations or structures.
Drainage type No pore pressures generated in soil clusters of the types
Undrained (C) and Non-porous; no generation of excess
pore pressures
Steady-state groundwater flow Possible, but not recommended since non-hydrostatic pore
pressures may lead to non-equilibrium.
Transient groundwater flow Not possible
Suction Included in steady-state pore pressures when suction is
allowed
Cavitation cut-off Not possible, but also not relevant; undrained behaviour is
ignored
Updated mesh Not possible (deformations are not considered)
Updated water pressures Not possible (deformations are not considered)
3 1 − 2νu 1 − 2ν0
Cluster just Kw = 0 Kw = 0 Kw = Kwsat · 10−8 Not relevant Not relevant
been activated
Non-porous or Kw = 0 Kw = 0 Kw = 0 Kw = 0 Kw = 0
dry cluster
APPENDIX C - CONVENTIONAL AND UNSATURATED SOIL BEHAVIOUR IN PLAXIS
The Python wrapper for the PLAXIS HTTP REST API hides the technicalities of dealing
with the HTTP communications behind a comfortable object-oriented wrapper.
Prerequisites:
• Install the latest Python 3.3.x release compatible with your system from
http://python.org/download/. The installer comes with a code editor called IDLE.
There are many alternative editors available which may be more to your liking, such
as PyScripter (https://code.google.com/p/pyscripter/). The rest of the document will
however assume IDLE as the editor of choice.
• At least a very rudimentary knowledge of the Python language is recommended, as
this document makes no attempt to teach it. Good resources for this purpose are:
• Dive Into Python 3 (http://www.diveintopython3.net/)
• After Hours Programming
(http://www.afterhoursprogramming.com/tutorial/Python/Overview/)
• Think Python: How to Think Like a Computer Scientist
(http://greenteapress.com/thinkpython/html/index.html)
• Your firewall must not block the PLAXIS application from accessing the internet, nor
must it block other applications (in particular the python.exe executable) from talking
to the remote scripting server embedded inside the PLAXIS application.
• Ensure you find an available port and then start the server. Information about ports
currently in use can be retrieved in a Windows Command Prompt, with the following
When using the remote scripting in an automated process, you can start the server
without manual interaction by launching your PLAXIS application with the AppServerPort
command line parameter. For example in order to start the server in PLAXIS 3D on port
21403:
c:\Program Files (x86)\Plaxis\PLAXIS 3D\Plaxis3DInput.exe --AppServerPort=21403
When the server is running, the main window of your PLAXIS application will reflect this
information:
• Start IDLE from the Windows Start menu. This will present you with a Python Shell.
• Define a variable indicating where your PLAXIS application is installed. For example:
>>> plaxis_path = r'c:\Program Files (x86)\Plaxis\PLAXIS 3D'
• Make sure to not include the trailing backslash. Note that the leading prompt ">>>"
is already present in the Python prompt, you do not have to type it yourself.
• Import the scripting library:
>>> import imp
>>> found_module = imp.find_module('plxscripting', [plaxis_path])
>>> plxscripting = imp.load_module('plxscripting', *found_module)
>>> from plxscripting.easy import *
• Connect to the PLAXIS application (21403 is the port number, substitute yours as
applicable):
>>> s, g = new_server('localhost', 21403)
• You should get a reply: 'OK'. If instead an error is returned, check the following for
the PLAXIS application you are trying to control:
• Is the application running?
• Does its title bar indicate that the server is active?
• Is the server active on the port you have specified when calling new_server?
You can also try a different port.
• Is your firewall blocking communications to the application?
Connecting to a PLAXIS application running on a computer other than your own requires
you to know the name or IP address of that machine. The IP address can be found by
typing in a Windows command prompt the following:
ipconfig | find "IPv4"
The result will give a set of four numbers separated by dots. E.g.:
IPv4 Address. . . . . . . . . . . : 10.4.2.178
In order to connect to this machine from within the same network, you can simply replace
'localhost' in the boilerplate connection code with the IP address you found, such as
'10.4.2.178'. For example:
>>> s, g = new_server('10.4.2.178', 21403)
If you want to connect to the same machine over the internet, you will need to set up
appropriate port forwarding rules in your router. Note that opening up your computer to
access from the network or even the internet carries significant security risks - do so only
if you have considered these risks.
All commands that are available in the PLAXIS command line are also available via the
Python wrapper. For more information, see the commands reference provided with your
PLAXIS application.
• Start by connecting from Python to PLAXIS 3D as described above.
• Execute the code displayed below. For the sake of completeness, the replies
returned by the different statements are included, but some of them are
session-specific and will be different for you. It is helpful while executing the Python
code to keep an eye on the command line in the PLAXIS application, in order to gain
an understanding of how Python code maps to commands. Keep also in mind that
In the examples above, we have looked at possibilities for using the API interactively from
a Python prompt. While this is very useful for testing and finding out how things work,
soon you will need the ability to save code and run it again without having to retype it.
• In IDLE make a new text file (File > New Window).
• Start the file with the boilerplate code to get access to the remote scripting libraries
(change the path to fit your installation)
plaxis_path = r'c:\Program Files (x86)\Plaxis\PLAXIS 3D' # no trailing backslash!
import imp
found_module = imp.find_module('plxscripting', [plaxis_path])
plxscripting = imp.load_module('plxscripting', *found_module)
After that, you can type the code, save and run it using Run > Run Module. A Python
script equivalent to the interactive session presented above looks as follows:
plaxis_path = r'c:\Program Files (x86)\Plaxis\PLAXIS 3D' # no trailing backslash!
import imp
found_module = imp.find_module('plxscripting', [plaxis_path])
plxscripting = imp.load_module('plxscripting', *found_module)
from plxscripting.easy import *
s, g = new_server('localhost', 21403)
s.new()
g.gotostructures()
g.lineload((3, 5, 0), (7, 5, 0))
g.gotomesh()
g.mesh(0.2)
g.gotostages()
phase1 = g.phase(g.Phases[0])
g.LineLoads[0].Active[phase1] = True
g.calculate()
g.save(r'c:\data\scripting_sample')
If you wish to reduce the boilerplate code required to get access to the plxscripting
module, see http://docs.python.org/3.3/using/cmdline.html#envvar-PYTHONPATH.
However, if you have multiple versions/applications that support remote scripting, not all
of them may be compatible with one particular version of the remote scripting API.
For some projects the overhead of updating the user interface may slow down the
execution of the script quite significantly. If you want maximum execution speed, launch
your PLAXIS application with the command parameter --NO_CONTROLLERS. For example:
Plaxis3DInput.exe --AppServerPort=21403 --NO_CONTROLLERS
The downside of this approach is that you will be unable to monitor progress or diagnose
problems visually. Instead, you will have to:
1. Stop the running script
2. Save the current project
3. Restart the application and open the saved project
While using Python scripts to manipulate the state of your PLAXIS application, it is
technically possible to interact with the same application in other ways as well (e.g. by
changing a project manually, or firing commands separately in an interpreter while the
If you encounter such an error, simply stop the execution of your script and restart it.
Opening a project, removing all lines containing beams and saving again
s.open(r'c:\data\lines_with_beams.p3d')
lines = [b.Parent for b in g.Beams]
g.delete(*lines)
# alternatively (look at the difference in the command line):
# g.delete(lines)
g.save()
Reporting all points that are not used by lines and that have no features
points_in_lines = set()
for line in g.Lines[:]:
for point in line:
points_in_lines.add(point)
points = set(g.Points)
points_no_lines = points - points_in_lines
points_to_report = []
for p in points_no_lines:
if len(p.UserFeatures.value) == 0:
points_to_report.append(p)
for p in points_to_report:
print(p.Name)
g.pointload(diagonal_points)
def get_bottom_point(line):
first, second = line.First, line.Second
delta = second.z - first.z
if abs(delta) < 1E-4: # horizontal line
return None
if delta > 0: # Second point is higher -> move First down
return first
else: # First point is higher -> move Second down
return second
def lengthen_embedded_piles(extra_depth):
# start by selecting the lines which have embedded piles
pile_lines = [pile.Parent for pile in g.EmbeddedPiles[:]]
lengthen_embedded_piles(2.5)
def phase_to_string(phase):
return "{} [{}]".format(phase.Identification, phase.Name)
def phase_to_string_with_error_code(phase):
return "{} [{}]. LogInfo: {}".format(
phase.Identification, phase.Name, phase.LogInfo)
def report_phases(phases_list):
successful_phases = []
failed_phases = []
uncalculated_phases = []
for phase in phases_list:
if phase.ShouldCalculate:
uncalculated_phases.append(phase)
elif phase.LogInfo == '0':
successful_phases.append(phase)
else:
failed_phases.append(phase)
report_lines = []
report_phases_by_kind(report_lines, successful_phases,
'Successful', phase_to_string)
report_phases_by_kind(report_lines, failed_phases,
'Failed', phase_to_string_with_error_code)
report_phases_by_kind(report_lines, uncalculated_phases,
'Skipped', phase_to_string)
return '\n'.join(report_lines)
def calculate_with_report():
phases_to_calculate = []
res = g.calculate()
report = [
'Title: {}'.format(g.Project.Title),
'Location: {}'.format(g.Project.Filename),
'',
report_phases(phases_to_calculate),
'',
'Command feedback:',
res]
return '\n'.join(report)
message = calculate_with_report()
# Run one command at a time and have the wrapper handle the response
# for you - not much extra control.
for command in commands:
print(s.call_and_handle_command(command))
# Run one command at a time and handle the response yourself. The
# response now consists of Python primitives representing the
# JSON communication rather than comfortable objects.
for command in commands:
print(s.call_commands(command))
# input parameters
nr_piles = 4 # must be > 1
total_angle_deg = 180 # over which to distribute the piles, must be > 0 and <= 360
pile_length = 10
pile_inclination_deg = 15
array_center = (5, 5, 0)
array_radius = 10
# calculations
def rel_coords_to_absolute(rel_coords, radius, offset):
coords = [radius * coord for coord in rel_coords]
return [p[0] + p[1] for p in zip(coords, offset)]
pile_inclination_rad = radians(pile_inclination_deg)
array_bottom_radius = array_radius + \
pile_length * sin(pile_inclination_rad)
array_bottom_center = (array_center[0], array_center[1],
array_center[2] - pile_length * cos(pile_inclination_rad))
pile_coords = []
for i in range(nr_piles):
angle = radians(total_angle_deg * i / (nr_piles - 1))
rel_coords = [cos(angle), sin(angle), 0]
top_coords = rel_coords_to_absolute(
rel_coords, array_radius, array_center)
bottom_coords = rel_coords_to_absolute(
rel_coords, array_bottom_radius, array_bottom_center)
pile_coords.append((top_coords, bottom_coords))
G.1 OVERVIEW
PLAXIS provides a HTTP based API. This API allows a developer or user to do normal
PLAXIS tasks such as creating projects, firing regular commands as found on the existing
command line, performing meshing and calculations. In addition, unique representations
of Plaxis objects can be retrieved, kept, and then used in later operations, such as the
fetching of their properties and methods, as well as calling methods. This provides
additional flexibility and power to the existing command line.
There are several categories of actions that can be performed via the HTTP API. Each of
these have their own URL path, as follows:
/environment Start, close, recover or open projects
/commands Fire regular PLAXIS command line commands and fetch any
created objects
/members Fetch the members of objects
/propertyvalues Fetch the property values of objects
/namedobjects Fetch existing objects by their command line names
/list Fetch items from listable objects
/enumeration Fetch all possible values of enumeration objects
Resources which create or retrieve objects return representations of those objects. These
representations comprise a Globally Unique Identifier (GUID), a type name, and a flag
indicating whether the object is listable or not. The JSON structure for this representation
is as follows:
"guid": <string>
"islistable": <boolean>
"ownerguid": <string>
"type": <string>
The server can recognise and retrieve information about objects contained within the
PLAXIS project based on that object's GUID. This means that, for example, a particular
PLAXIS object can be created using the /commands resource, and then later the values
of that object's properties can be retrieved by sending that object's GUID in a request to
the /propertyvalues resource. The same GUID can be used in the creation of additional
objects by sending it in a request to the /commands resource. If the object is listable,
sending the guid as part of a request to the /list resource will enable access to the
object's properties by index, and so on. The ownerguid is only returned for objects that
are intrinsic properties. It is the guid of the owning object.
It is possible to perform a subset of the actions that are described below by using a
browser address bar. This can be useful for quickly viewing the data representations that
G.4.1 /environment
This resource is for managing projects. Creating, restoring, closing or opening a project
are permitted. If a project is opened, then a filename string is supplied for locating an
existing project locally, or otherwise becomes the name of a new project.
Note that it is not possible to save projects via this resource. Instead, fire a save
command using the /commands resource.
Method: POST
Data parameters
{
"action":
{
"name": *1 <string> from {"new", "open", "close", "restore"}*,
"filename": *1 <string>*
}
}
Success responses
Code Reason
200 OK
Error responses
Code Reason
404 STATUS NOT FOUND "Incorrect or unspecified action: [action-name]"
404 STATUS NOT FOUND "Incorrect or unspecified file name, or not a valid project: [filename]"
400 BAD REQUEST "Syntax of request could not be understood by the server."
415 UNSUPPORTED MEDIA TYPE "Unsupported content-encoding of [encoding type]"
G.4.2 /commands
This resource allows all standard commands (see command reference pages) to be run
on an active PLAXIS project. Any objects created as the result of commands are
represented in the response.
The server can recognise and "translate" GUIDs that refer to existing objects, as part of
the command string. This means that it is possible to use GUIDs in the place of regular
command line names to refer to other objects. This is beneficial, because the names of
objects may be changed over the lifetime of a project (e.g. using the rename command),
but the GUID remains fixed.
For example, it is possible to create a line between two existing points by referring to the
GUIDs of those points in the command string of a request:
line <guid1> <guid2>
Note that a command that fails on the command line will not result in an error response
from the server. These are treated as successful requests to the API. The server
response indicates whether a command failed or not in the 'success' field. (See example
below.) Whether the command succeeded or not, the 'additionalinfo' field indicates the
response from that action as seen on the regular command line.
Method: POST
Data parameters
{
"action":
{
"commands":
[
*{1, ...}[<string>]*
]
}
}
Success responses
Code Reason
200 OK
Data
{
"commands":
[
{
"feedback":
{
"extrainfo": "Added Point_1",
"returnedobjects":
[
{
"islistable": true,
"type": "Point",
"guid": "{860C1FC5-BDE1-4887-98A2-8FE20B61FCD5}"
}
],
"debuginfo": "",
"success": true,
"errorpos": -1
},
"command": "point 1 2 3"
}
]
}
{
"commands": [
{
"feedback": {
"extrainfo": "Added Point_3\r\nAdded Line_1\r\nAdded EmbeddedPile_1",
"returnedobjects": [
{
"islistable": true,
"type": "Point",
"guid": "{92851C68-561A-4C21-9ACD-46F50905FD4E}"
},{
"islistable": true,
"type": "Line",
"guid": "{7DB7DF8E-1FE4-4624-A3D2-F682DB9C8F18}"
},{
"islistable": false,
"type": "EmbeddedPile",
"guid": "{A678F396-C8F3-41EE-9A8E-5798BB0B6054}"
}],
"debuginfo": "",
"success": true,
"errorpos": -1
},
"command": "embeddedpile 4 4 4 Point_1"
}]
}
Code Reason
200 OK
Data
{
"commands":
[
{
"feedback":
{
"extrainfo": "Command \"bad\" is not recognized as a global command.
Check the spelling or try specifying a target object.",
"debuginfo": "",
"success": false,
"errorpos": 1
},
"command": "bad"
}
]
}
Error responses
Code Reason
400 BAD REQUEST "Syntax of request could not be understood by the server."
404 STATUS NOT FOUND "Commands not specified"
405 METHOD NOT ALLOWED "No active project"
415 UNSUPPORTED MEDIA TYPE "Unsupported content-encoding of [encoding type]"
G.4.3 /members
This resource fetches the members (properties and methods) of an existing object within
PLAXIS, identified by its GUID.
Property members are represented as objects.The type names of these objects are
simplified in comparison to their type names within PLAXIS. This means that all 'Text'
subtypes for instance will be simply described as having a type of 'Text'. The same
principle applies to numerical, boolean, enumeration and object types.
Method members are represented as names only.
Method: POST
Data Parameters
{
"action":
{
"members":
[
*{1,...}[GUID string]*
]
}
}
Success responses
Code Reason
200 OK
Data
{
"queries": {
"{68A2C1A7-924B-46C3-89DF-F597DAADC359}": {
"extrainfo": "",
"success": true,
"properties": {
"z": {
"islistable": false,
"ownerguid": "{68A2C1A7-924B-46C3-89DF-F597DAADC359}",
"type": "Number",
"guid": "{5FCCF1CC-7576-4565-B931-DE54241EA21B}"
},
"y": {
"islistable": false,
"ownerguid": "{68A2C1A7-924B-46C3-89DF-F597DAADC359}",
"type": "Number",
"guid": "{D52CBED3-8C59-4967-A5BA-89ECE6E9BD94}"
},
"x": {
"islistable": false,
"ownerguid": "{68A2C1A7-924B-46C3-89DF-F597DAADC359}",
"type": "Number",
"guid": "{25DA6A75-4F25-42F3-9654-3672A93FCE1F}"
},
"Name": {
"islistable": false,
"ownerguid": "{68A2C1A7-924B-46C3-89DF-F597DAADC359}",
"type": "Text",
"guid": "{A4209590-9487-4B98-9966-34A6A58DBE2D}"
},
"Comments": {
"islistable": false,
"ownerguid": "{68A2C1A7-924B-46C3-89DF-F597DAADC359}",
"type": "Text",
"guid": "{9430592C-7D5D-4735-B572-2B29C87FBF1F}"
}
},
"commands": [
"echo","__dump","commands","rename","set","info","__observers","setproperties","move","rotate","rotateline"],
"commandlinename": "Point_1"
}
}
}
{
"queries": {
"{0D04D0CB-86BF-4430-882C-FA56E45DF7AB}": {
"extrainfo": "GUID does not refer to object in registry: {0D04D0CB-86BF-4430-882C-FA56E45DF7AB}",
"success": false
}
}
}
{
"queries": {
"bad": {
"extrainfo": "Supplied string is not a valid GUID: bad",
"success": false
}
}
}
Error responses
Code Reason
400 BAD REQUEST "Syntax of request could not be understood by the server."
404 STATUS NOT FOUND "No GUIDs specified"
405 METHOD NOT ALLOWED "No active project"
415 UNSUPPORTED MEDIA TYPE "Unsupported content-encoding of [encoding type]"
G.4.4 /propertyvalues
This resource provides the current values of particular properties of PLAXIS objects. If a
property is itself an object (such as the point properties of a line), then it is represented as
a full object. If a property is a primitive value type, then the relevant primitive value is
returned.
There are two ways of calling this resource. One is to supply just a guid of the owner
object, in which case all properties of this object are represented in the response. The
other way is to supply not just the owner guid, but also a property name and an optional
phase parameter. Such a structure is treated as a request for one single property. If that
property is staged, and the supplied phase guid is valid, then the value for that property
and phase is returned.
Method: POST
Data Parameters
{
"action":
{
"propertyvalues":
[*{1} GUID string*]
}
}
{
"action":
{
"propertyvalues":
{
"owner": *{1} GUID string*
("phaseguid: *{1} GUID string*)
}
}
}
Success responses
{
"queries": {
"{D4981EA8-EC07-471D-A8C3-7510B4F74F25}": {
"extrainfo": "",
"success": true,
"properties": {
"Second": {
"islistable": true,
"type": "Point",
"guid": "{A67959C9-C12D-41D6-8136-52C1E5711165}"
},
"Length": 6.92820323027551,
"First": {
"islistable": true,
"type": "Point",
"guid": "{4B889C3F-7F5C-4611-8D31-F4A39C71B11A}"
},
"AxisFunction": 0,
"Name": "Line_1",
"AxisVectorZ": 0,
"AxisVectorY": 0,
"AxisVectorX": 0,
"Comments": ""
}
}
}
}
{
"queries":
{
"{A10835D7-AE38-449C-BF67-6C9CE22705BD}":
{
"extrainfo": "",
"properties": {"Comments": ""},
"success": true
}
}
}
Error responses
Code Reason
400 BAD REQUEST "Syntax of request could not be understood by the server."
404 STATUS NOT FOUND "No GUIDs specified"
405 METHOD NOT ALLOWED "No active project"
415 UNSUPPORTED MEDIA TYPE "Unsupported content-encoding of [encoding type]"
G.4.5 /namedobjects
Supplies a representation of objects within PLAXIS which are identified by their
command line name. This is useful when accessing internal lists of objects within Plaxis
with known names, such as lists of particular user features or geometric objects. This
resource can also be used to retrieve unique representations of objects for which the
command line name is known, such as 'Point_1'.
Method: POST
Data Parameters
Example Request:
{
"action":
{
"namedobjects":
[
"namedobjects": ["Points"]
]
}
}
{
"action":
{
"namedobjects":
[
"namedobjects": ["Beam_1"]
]
}
}
Success Responses
Code Reason
200 OK
Data
{
"namedobjects": {
"Points": {
"extrainfo": "",
"success": true,
"returnedobject": {
"islistable": true,
"type": "ModelGroup",
"guid": "{32871FF1-809C-470A-94B0-4B56B28BC67B}"
}
}
}
}
Code Reason
200 OK
Data
{
"namedobjects": {
"Beam_1": {
"extrainfo": "",
"success": true,
"returnedobject": {
"islistable": false,
"type": "Beam",
"guid": "{C0C4EDA6-DFF7-4593-A280-A1EC20D5EADD}"
}
}
}
}
Error responses
Code Reason
400 BAD REQUEST "Syntax of request could not be understood by the server."
404 STATUS NOT FOUND "No object names specified"
405 METHOD NOT ALLOWED "No active project"
415 UNSUPPORTED MEDIA TYPE "Unsupported content-encoding of [encoding type]"
G.4.6 /list
Provides the ability to perform read operations on lists within PLAXIS, such as getting
values at particular indices, finding the number of objects in the list, and filtering.
Method: POST
Data Parameters
Example Request:
{"action": {"listqueries": [{"guid": "{CF1DECEB-A28D-4609-B38C-8D3FF5E573A1}",
"method": "sublist",
"startindex": 0,
"stopindex": 1]}}
Success Responses
{"listqueries": [{"extrainfo": "",
"guid": "{CF1DECEB-A28D-4609-B38C-8D3FF5E573A1}",
"methodname": "sublist",
"outputdata": [{"guid": "{81A754C8-1393-48BA-AD6E-5F6953F9FBE7}",
"islistable": true,
"type": "Point"}],
"startindex": 0,
"stopindex": 1,
"success": true}]}
Error responses
Code Reason
400 BAD REQUEST "Syntax of request could not be understood by the server."
405 METHOD NOT ALLOWED "No active project"
415 UNSUPPORTED MEDIA TYPE "Unsupported content-encoding of [encoding type]"
G.4.7 /enumeration
Provides a list of all possible enumeration values of one or more enumeration objects,
when supplied with their GUIDs.
Success Responses
Code Reason
200 OK
Data
{
"queries": {
"{31AFFD41-0EB4-4A67-94F4-298642E17507}": {
"extrainfo": "",
"success": true,
"enumvalues": {
"hinged": 1,
"free": 2,
"rigid": 0
}
}
}
}
Code Reason
200 OK
Data
{
"queries": {
"{0B4774CE-3307-4B86-B644-BFDE82CF3AAA}": {
"extrainfo": "GUID does not refer to object in registry:
{0B4774CE-3307-4B86-B644-BFDE82CF3AAA}",
"success": false
}
}
}
Code Reason
200 OK
Data
{
"queries": {
"{C896C54E-1CBA-4B0F-9987-E1FF57A9131E}": {
Error responses
Code Reason
400 BAD REQUEST "Syntax of request could not be understood by the server."
405 METHOD NOT ALLOWED "No active project"
415 UNSUPPORTED MEDIA TYPE "Unsupported content-encoding of [encoding type]"