PLAXIS - 2D - CEV21 - Tutorial - 08 - Construction of A Road Embankment
PLAXIS - 2D - CEV21 - Tutorial - 08 - Construction of A Road Embankment
PLAXIS - 2D - CEV21 - Tutorial - 08 - Construction of A Road Embankment
• Consolidation analysis
• Modelling drains
• Change of permeability during consolidation
• Safety analysis (strength reduction)
• Updated mesh analysis (large deformations)
Geometry
The embankment is 16 m wide and 4 m high. The slopes have an inclination of 1:3. The problem is symmetric, so
only one half is modelled (in this case the right half is chosen). The embankment itself is composed of loose
sandy soil. The subsoil consists of 6 m of soft soil. The upper 3 m is peat and the lower 3 m is clay. The phreatic
level is located 1 m below the original ground surface. Under the soft soil layers there is a dense sand layer of
which 4 m are considered in the model.
12 m 16 m 12 m
Road embankment 4m
Peat 3m
Clay 3m
Dense sand
A number of material set are needed for this tutorial. The properties of the materials are as follows:
Table 1: Material properties of the sand and clay layer and the interfaces
General
Hardening Hardening
Material model - Soft soil Soft soil -
soil soil
Undrained Undrained
Drainage type - Drained Drained -
(A) (A)
Parameters
Unloading / reloading
Eurref 75· 103 105 · 103 - - kN/m2
stiffness
Dilatancy angle ψ 0 3 0 0 °
Groundwater
Interfaces
Initial
1.
Click the Materials button to open the Material sets window.
2. Create soil material data sets according to Table 20 (on page 4) and assign them to the corresponding layers
in the borehole (see Figure 91 (on page 3)).
3. Close the Modify soil layers window and proceed to the Structures mode to define the embankment and
drains.
Note: The initial void ratio (einit) and the change in permeability (ck) should be defined to enable the modelling
of a change in the permeability in a consolidation analysis due to compression of the soil. This option is
recommended when using advanced models.
The embankment and the drains are defined in the Structures mode.
1. Click the Structures tab to proceed with the input of structural elements in the Structures mode.
1. Click the Create soil polygon button in the side toolbar and select the Create soil polygon option.
2. Define the embankment in the drawing area by clicking on (0 0), (0 4), (8 4) and (20 0).
3. Right-click the created polygon and assign the Embankment data set to the soil polygon.
In this project the effect of the drains on the consolidation time will be investigated by comparing the results
with a case without drains. Drains will only be active for the calculation phases in the case with drains.
1. Click the Create hydraulic conditions button in the side toolbar and select the Create drain option in the
appearing menu.
2. Drains are defined in the soft layers (clay and peat; y = 0 m to y = -6 m). The distance between two
consecutive drains is 2 m. Considering the symmetry, the first drain is located at 1 m distance from the model
boundary. 10 drains will be created in total. The head is defined at 0.0 m.
Note:
The modelling of drains in a plane strain model actually involves the use of an equivalent (lateral) permeability
in the surrounding soil based on the drain pattern. The latter has been omitted in this simplified example. More
information can be in found in literature1.
1 Achtergronden bij numerieke modellering van geotechnische constructies, deel 2. CUR 191. Stichting CUR,
Gouda Indraratna, B.N., Redana, I.W., Salim, W. (2000), Predicted and observed behaviour of soft clay
foundations stabilised with vertical drains. Proc. GeoEng. 2000, Melbourne.
The initial water pressures are fully hydrostatic and based on a general phreatic level located at y = -1 m. Note
that a phreatic level is automatically created at y = -1 m, according to the value specified for Head in the
borehole. In addition to the phreatic level, attention must be paid to the boundary conditions for the
consolidation analysis that will be performed during the calculation process. Without giving any additional
input, all boundaries except for the bottom boundary are draining so that water can freely flow out of these
boundaries and excess pore pressures can dissipate. In the current situation, however, the left vertical boundary
must be closed because this is a line of symmetry, so horizontal flow should not occur. The remaining
boundaries are open because the excess pore pressures can be dissipated through these boundaries. In order to
define the appropriate consolidation boundary conditions, follow these steps:
Consolidation analysis
A consolidation analysis introduces the dimension of time in the calculations. In order to correctly perform a
consolidation analysis a proper time step must be selected. The use of time steps that are smaller than a critical
minimum value can result in stress oscillations.
The consolidation option in PLAXIS 2D allows for a fully automatic time stepping procedure that takes this
critical time step into account. Within this procedure there are three main possibilities:
Consolidate for a predefined period, including the effects of changes to the active geometry (Staged
construction).
Consolidate until all excess pore pressures in the geometry have reduced to a predefined minimum
value (Minimum excess pore pressure).
The first two possibilities will be used in this exercise. To define the calculation phases, follow these steps:
2.
In the Phases window select the Consolidation option from the Calculation type drop-down menu in
the General subtree.
3.
Make sure that the Staged construction option is selected for the Loading type.
4. Enter a Time interval of 1 day. The default values of the remaining parameters will be used.
5. In the Staged construction mode activate the second part of the embankment.
Safety analysis
In the design of an embankment it is important to consider not only the final stability, but also the stability
during the construction. It is clear from the output results that a failure mechanism starts to develop after the
second construction phase.
It is interesting to evaluate a global safety factor at this stage of the problem, and also for other stages of
construction.
To calculate the global safety factor for the road embankment at different stages of construction, follow these
steps:
1. Select Phase 1 in the Phases explorer.
2.
Add a new calculation phase.
3. Double-click on the new phase to open the Phases window.
4. In the Phases window the selected phase is automatically selected in the Start from phase drop-down
menu.
5.
In the General subtree, select Safety as calculation type.
6.
The Incremental multipliers option is already selected in the Loading input box. The first increment of
the multiplier that controls the strength reduction process, Msf, is set to 0.1.
7. In order to exclude existing deformations from the resulting failure mechanism, select the Reset
displacements to zero option in the Deformation control parameters subtree.
8. The first safety calculation has now been defined.
9. Follow the same steps to create new calculation phases that analyse the stability at the end of each
consolidation phase.
Calculate
Before starting the calculation it is suggested that you select nodes or stress points for a later generation of load-
displacement curves or stress and strain diagrams. To do this, follow the steps given below.
1. Click the Select points for curves button in the side toolbar.
2. As the first point, select the toe of the embankmentat (20 0).
3. The second point will be used to plot the development (and decay) of excess pore pressures. To this end, a
point somewhere in the middle of the soft soil layers at the left side of the model is needed, hence underneath
the middle of the embankment. For instance at (0 -3).
4. Click the Calculate button to calculate the project.
During a consolidation analysis the development of time can be viewed in the upper part of the calculation info
window.
In addition to the multipliers, a parameter Pexcess,max occurs, which indicates the current maximum excess pore
pressure. This parameter is of interest in the case of a Minimum excess pore pressure consolidation analysis,
where all pore pressures are specified to reduce below a predefined value.
Results
After the calculation has finished, select the third phase and click the View calculation results button .
The Output window now shows the deformed mesh after the undrained construction of the final part of the
embankment. Considering the results of the third phase, the deformed mesh shows the uplift of the embankment
toe and hinterland due to the undrained behaviour.
1. Press <Ctrl + 7> to display the developed excess pore pressures (see Appendix C of the Reference Manual for
more shortcuts). They can also be displayed by selecting the menu Stresses > Pore pressures > P excess.
2.
Click the Center principal directions button . The principal directions of excess pressures are displayed
at the center of each soil element. The results are displayed in the following figure.
It is clear that the highest excess pore pressure occurs under the embankment centre.
1. Select Phase 4 in the drop down menu.
2. Click the Contour lines button in the toolbar to display the results as contours.
3.
Use the Draw scanline button or the corresponding option in the View menu to define the position of
the contour line labels.
It can be seen that the settlement of the original soil surface and the embankment increases considerably during
the fourth phase. This is due to the dissipation of the excess pore pressures (= consolidation), which causes
further settlement of the soil. The figure below shows the remaining excess pore pressure distribution after
consolidation. Check that the maximum value is below 1.0 kN/m2.
Figure 15: Excess pore pressure contours after consolidation to Pexcess < 1.0 kN/m2
The Curves manager can be used to view the development, with time, of the excess pore pressure under the
embankment. In order to create such a curve, follow these steps:
1. Create a new curve by clicking the Curves manager button .
2. For the x-axis, select the Project option from the drop-down menu and select Time in the tree.
3. For the y-axis select the point in the middle of the soft soil layers (Point B) from the drop-down menu. In the
tree select Stresses > Pore pressure > p excess.
4. Select the Invert sign option for the y-axis.
5. Click OK.
6. Open the Curve settings (F3) and go to the second tabsheet.
7. In the Show box click the Phases button. By default all phases are selected to show in the curve. For the
clarity of the curve, hide the Safety phases (phases 5 - 8).
8. Click OK to close the Curve settings window.
A curve similar to the following one should appear:
The figure above clearly shows the four calculation phases. During the construction phases the excess pore
pressure increases with a small increase in time while during the consolidation periods the excess pore pressure
decreases with time. In fact, consolidation already occurs during construction of the embankment, as this
involves a small time interval. From the curve it can be seen that more than 50 days are needed to reach full
consolidation.
Save the chart before closing the Output program.
Additional displacements are generated during a Safety calculation. The total displacements do not have a
physical meaning, but the incremental displacements in the final step (at failure) give an indication of the likely
failure mechanism.
In order to view the mechanisms in the three different stages of the embankment construction:
1.
Select one of these phases and click the View calculation results button .
2. Select the menu Deformations > Incremental displacements > |Δu|.
3. Change the presentation from Arrows to Shadings . The resulting plots give a good impression of the
failure mechanisms. The magnitude of the displacement increments is not relevant.
Figure 17: Shadings of the total displacement increments indicating the most applicable failure mechanism of the
embankment in the final stage
The safety factor can be obtained from the Calculation info option of the Project menu. The Multipliers
tabsheet of the Calculation information window represents the actual values of the load multipliers. The value
of ΣMsf represents the safety factor, provided that this value is indeed more or less constant during the previous
few steps.
The best way to evaluate the safety factor, however, is to plot a curve in which the parameter ΣMsf is plotted
against the displacements of a certain node. Although the displacements are not relevant, they indicate whether
or not a failure mechanism has developed.
In order to evaluate the safety factors for the three situations in this way, follow these steps:
1. Click the Curves manager button in the toolbar.
2. Click New in the Charts tabsheet.
3. In the Curve generation window, select the embankment toe (Point A) for the x-axis. Select Deformations >
Total displacements > |u|.
4. For the y-axis, select Project > Multipliers > ΣMsf. The Safety phases are considered in the chart.
5. Right-click on the chart and select the Settings option in the appearing menu. The Settings window pops up.
6. In the tabsheet corresponding to the curve click the Phases button.
7. In the Select phases window select Phase 5:
11. In the Settings window click the Chart tab to open the corresponding tabsheet.
12. In the Chart tabsheet specify the chart name.
13. Set the scaling of the x-axis to Manual and set the value of Maximum to 1:
Using drains
In this section the effect of the drains in the project will be investigated. Four new phases will be introduced
having the same properties as the first four consolidation phases. The first of these new phases should start from
the initial phase. The differences in the new phases are:
• The drains should be active in all the new phases. Activate them in the Staged construction mode.
• The Time interval in the first three of the consolidation phases (9 to 11) is 1 day. The last phase is set to
Minimum excess pore pressure and a value of 1.0 kN/m2 is assigned to the minimum pressure (|P-stop|).
Follow these steps:
1. After the calculation is finished, save the project, then select the last phase and click the View calculation
results button . The Output window now shows the deformed mesh after the drained construction of the
final part of the embankment. In order to compare the effect of the drains, the excess pore pressure
dissipation in the second point can be used.
2. Click the Curves manager button to open the Curves manager.
3. In the Chart tabsheet double-click Chart 1 (pexcess of the second point at (0 -3) versus time). The chart is
displayed. Close the Curves manager.
4. Double-click the curve in the legend at the right of the chart. The Settings window pops up.
5. Click the Add curve button and select the From current project ... option in the appearing menu. The Curve
generation window pops up.
6. Select the Invert sign option for y-axis and click OK to accept the selected options.
7. In the chart a new curve is added and a new tabsheet corresponding to it is opened in the Settings window.
Click the Phases button. From the displayed window select the Initial phase and the last four phases
(drains) and click OK.
8. In the Settings window change the titles of the curves in the corresponding tabsheets.
9. In the Chart tabsheet specify the chart name.
10. Click Apply to preview the generated curve and click OK to close the Settings window. The chart gives a
clear view of the effect of drains in the time required for the excess pore pressures to dissipate:
Note:
Instead of adding a new curve, the existing curve can be regenerated using the corresponding button in the
Curves settings window.
Figure 22: Effect of updated mesh and water pressures analysis on resulting settlements
It can be seen that the settlements are less when the Updated mesh and Updated water pressures options are
used (red curve). This is partly because the Updated mesh procedure includes second order deformation effects
by which changes of the geometry are taken into account, and partly because the Updated water pressures
procedure results in smaller effective weights of the embankment. This last effect is caused by the buoyancy of
the soil settling below the (constant) phreatic level. The use of these procedures allows for a realistic analysis of
settlements, taking into account the positive effects of large deformations.