Tutorial 07 Finite Element Groundwater Seepage
Tutorial 07 Finite Element Groundwater Seepage
Tutorial 07 Finite Element Groundwater Seepage
7-1
Introduction
Within the Slide program, Slide has the capability to carry out a finite element groundwater seepage analysis for steady state or transient conditions. This tutorial will demonstrate steady state seepage. Transient groundwater seepage is discussed in later tutorials. Finite element groundwater analysis in Slide allows you to define and analyze a groundwater problem using the same model as for the slope stability problem. The boundaries of the problem only need to be defined once, and will be used for both the groundwater analysis and the slope stability analysis. After a groundwater seepage analysis is performed, the results (pore pressures), can be automatically utilized by the slope stability analysis engine in Slide. NOTE: The groundwater analysis capability in Slide can be considered a completely self contained groundwater analysis program, and can be used independently of the slope stability functionality of Slide. You may perform a groundwater analysis in Slide, without necessarily performing a slope stability analysis.
Slide v.6.0
Tutorial Manual
7-2
Although the Slide groundwater analysis is geared towards the calculation of pore pressures for slope stability problems, it is not restricted to slope geometry configurations. The groundwater modeling and analysis capabilities in Slide can be used to analyze an arbitrary, 2-dimensional groundwater problem, for saturated / unsaturated flow conditions.
Groundwater Modeling
The groundwater modeling options in Slide are all contained within the Slide Model program. In order to enable steady state groundwater modeling, it is first necessary to set the Groundwater Method in Project Settings to Steady State FEA. When you do this: an Analysis Mode option will be available in the toolbar, which allows you to select either Slope Stability analysis mode, or Steady State Groundwater mode. When you are in Groundwater analysis mode, the menu and toolbar will present all of the necessary groundwater modeling options.
The following general procedure is required, in order to use the Slide program for a groundwater analysis.
Project Settings
In order to perform a steady state groundwater analysis, the first thing you must do is set the Groundwater Method in Project Settings: 1. Select Project Settings from the toolbar or the Analysis menu. 2. Select the Groundwater page, and set the Groundwater Method = Steady State FEA (note: FEA stands for Finite Element Analysis). 3. You may configure the Groundwater Analysis parameters as necessary (i.e. Tolerance or Maximum Number of Iterations) by selecting the FEA Options button. 4. Select OK. The Slide toolbar will now display an Analysis Mode option, as described below.
Slide v.6.0
Tutorial Manual
7-3
The Slope Stability analysis mode allows you to define the elements of your model which are relevant to the slope stability analysis. The Groundwater analysis mode allows you to define the elements of your model which are relevant to the groundwater analysis.
Alternatively, you can change the analysis mode by selecting the tabs which appear at the lower left of the application window, as shown below.
When you switch between Slope Stability analysis mode and Groundwater mode, you will notice that the menus and toolbars are automatically updated.
Slide v.6.0
Tutorial Manual
7-4
Boundaries
The SAME boundaries are used in both the groundwater and slope stability analysis. However, the model boundaries can only be defined when the Analysis Mode = Slope Stability. You MUST create the model boundaries, when Analysis Mode = Slope Stability. You CANNOT create or edit the model boundaries, when Analysis Mode = Groundwater.
For details about defining boundaries in Slide, see the previous tutorials, or the Slide Help system.
Meshing
The groundwater analysis in Slide is a finite element analysis, and therefore a finite element mesh is required in order to solve the problem. To create the finite element mesh: You can do this with a single mouse click, by selecting the Discretize and Mesh option. This will automatically create a graded finite element mesh within your model boundaries. If the mesh requires customization, several options are available in the Mesh menu, which allow you to customize the mesh as necessary.
Slide v.6.0
Tutorial Manual
7-5
Boundary Conditions
Once the finite element mesh is satisfactory, the user must then set up the boundary conditions which define the groundwater problem. This is done with the Set Boundary Conditions option, which allows you to define all of the necessary pressure and flow boundary conditions along the model boundaries.
Hydraulic Properties
The permeability (hydraulic conductivity) characteristics of each material are defined with the Define Hydraulic Properties option. You may define a saturated permeability for each material. In addition, various models are available for defining the unsaturated permeability, or you may create a user-defined permeability function.
Slide v.6.0
Tutorial Manual
7-6
Groundwater Compute
When all of your groundwater boundary conditions and material properties have been defined, then you are ready to run the groundwater analysis. The groundwater analysis engine in Slide is a separate program from the slope stability analysis engine. It allows you to run the groundwater analysis independently of the slope stability analysis. To run the analysis: select the Compute (Groundwater) option, from the Analysis menu or the toolbar. This will run the Slide Groundwater Analysis engine.
You will see a Compute dialog while the analysis is running. When the analysis is finished, you will be returned to the Slide Model program.
After you have computed the groundwater analysis, you should view the results of the groundwater analysis, by selecting the Interpret (groundwater) option. If the groundwater analysis results are satisfactory, then you can return to the Slide Model program, switch the Analysis Mode = Slope Stability, and proceed with your slope stability analysis. When you select the slope stability Compute option, the slope stability analysis will automatically use the pore pressures calculated from the groundwater analysis.
Groundwater Interpret
The results of a groundwater analysis are viewed with the Slide Interpret program, using the options in the Groundwater menu. NOTE: Groundwater analysis results can be viewed simultaneously with the slope stability analysis results. The display of groundwater results (i.e. pore pressure contours), and the display of slope stability results (i.e. slip surfaces, safety factors etc) are fully integrated within the Slide Interpret program.
Slide v.6.0
Tutorial Manual
7-7
If you only wish to view groundwater results without the slope stability results, or vice versa, then you can easily turn display options on or off, as necessary.
To run the Slide Interpret program, select the Interpret (groundwater) option, from the Analysis menu or the toolbar, after you have computed the groundwater analysis with the Compute (Groundwater) option.
After you perform a groundwater analysis with Slide, it is always a good idea to first use the Interpret program, to check that the analysis results are reasonable. If not, then you should go back to the Slide Model program, and check that you have defined your model correctly. A simple introduction to groundwater modeling and data interpretation using Slide is found in the next tutorial.
Slide v.6.0
Tutorial Manual
7-8
This tutorial will demonstrate a simple groundwater seepage analysis using Slide. We will begin with the same model used for Tutorial 05 (the Water Pressure Grid tutorial). However, rather than read in a Water Pressure Grid from a file, we will carry out a seepage analysis in order to determine the pore pressures within the slope. We will then re-run the slope stability analysis, and compare the results with Tutorial 05. MODEL FEATURES: homogeneous, single material slope steady state finite element seepage total head boundary conditions
The finished product of this tutorial can be found in the Tutorial 07 Finite Element Groundwater Seepage.slim data file. All tutorial files installed with Slide 6.0 can be accessed by selecting File > Recent Folders > Tutorials Folder from the Slide main menu.
Slide v.6.0
Tutorial Manual
7-9
Model
First, lets read in the Tutorial 05 file. From the Slide main menu, select File > Recent Folders > Tutorials Folder. You will see the Open File dialog. Open the Tutorial 05 Water Pressure Grid.slim file.
Project Settings
In order to perform a finite element groundwater analysis with Slide, the first thing you must always do, is set the Groundwater Method in Project Settings.
You will see a warning message that the water pressure grid and water table will be deleted. Select Yes. Note: The Water Pressure Grid (the grid of blue triangles) has disappeared from the model. Since we will be obtaining pore pressures from the finite element analysis, the Water Pressure Grid will not be used, and has been deleted. The Water Table which was used for defining the ponded water in Tutorial 5, has also been deleted. When we define the groundwater boundary conditions, ponded water will automatically be created based on the total head boundary conditions.
Slide v.6.0
Tutorial Manual
7-10
The Analysis Mode option now appears in the Slide toolbar. The Analysis Mode allows you to switch between Slope Stability Analysis Mode, and Groundwater Analysis Mode.
To begin with, we will use Slope Stability Analysis mode, since we have to do a bit of editing to the boundaries. (Boundary defining and editing can only be done in Slope Stability Analysis Mode).
Boundary Editing
Before we proceed with the groundwater modeling, we have to make one small modification to the external boundary. To correctly define the boundary conditions for the groundwater problem, we need to add two vertices to the external boundary. Select the Add Vertices option from the toolbar.
The two vertices you have added define the level of the water table (phreatic surface) at the slope face, and at the right edge of the model. These vertices will be necessary in order to correctly assign the total head boundary conditions.
Slide v.6.0
Tutorial Manual
7-11
When the Analysis Mode = Groundwater, you will notice that the menus and toolbars are updated: The modeling options relevant to the Groundwater Analysis are now available. The modeling options which are relevant to the slope stability analysis, are now hidden.
Also notice that the Search Grid (for the circular surface Grid Search) and the Slope Limits symbols, are not displayed, when you are in Groundwater Analysis Mode. In general: Modeling entities which are only applicable to the SLOPE STABILITY ANALYSIS are NOT displayed when ANALYSIS MODE = GROUNDWATER. Modeling entities which are only applicable to the GROUNDWATER ANALYSIS are NOT displayed when ANALYSIS MODE = SLOPE STABILITY.
Slide v.6.0
Tutorial Manual
7-12
Meshing
The groundwater analysis in Slide is a finite element analysis, and therefore a finite element mesh is required in order to solve the problem. You can create the finite element mesh, with a single mouse click, by selecting the Discretize and Mesh option. In one step, this will automatically: 1. Discretize the model boundaries (the discretization of the boundaries forms the framework of the finite element mesh) 2. Generate a graded finite element mesh within the model boundaries.
Slide v.6.0
Tutorial Manual
7-13
However, note that Slide allows total user control over the generation and customization of the mesh. If the mesh requires customization, many different options are available in the Mesh menu, which allow the user to customize the mesh as necessary. Experimenting with the meshing options is beyond the scope of this tutorial. The user is encouraged to experiment with the mesh options, after completing this tutorial. For details about the use of the mesh options, refer to the Slide Help system.
Boundary Conditions
After the finite element mesh has been generated, you must define the boundary conditions which define the groundwater problem you wish to solve. This is done with the Set Boundary Conditions option, which allows you to define all of the necessary pressure and flow boundary conditions along the model boundaries. Notice that by default, when the mesh is generated: the slope surface is given an Unknown boundary condition (P=0 or Q= 0) the bottom and side edges of the external boundary are given the default Zero Nodal Flow boundary condition (no symbol is displayed).
Slide v.6.0
Tutorial Manual
7-14
1. First select Zoom All (you can use the F2 function key), to make sure the model is fully zoomed in the view. 2. Now, make sure the Total Head boundary condition option is selected in the Set Boundary Conditions dialog, as shown above. 3. In the dialog, enter a Total Head Value = 26 meters. Also make sure the Selection Mode = Boundary Segments. 4. Now you must select the desired boundary segments, by clicking on them with the mouse. 5. Click on the THREE segments of the external boundary indicated in Figure 7-3. (i.e. the left edge of the external boundary, and the two segments at the toe of the slope).
Slide v.6.0
Tutorial Manual
7-15
Slide v.6.0
Tutorial Manual
7-16
Hydraulic Properties
The last thing we have to do to complete the groundwater model, is define the hydraulic properties (permeability) of the slope material.
In the Define Hydraulic Properties dialog, enter a saturated permeability Ks = 5e-8. Select OK. NOTE: the model used for this tutorial and the Water Pressure Grid tutorial, is based on a model from the 1989 ACADS soil slope stability programs review. (Giam, P.S.K. & I.B. Donald 1989, Example problems for testing soil slope stability programs, Civil Engineering Research Report No. 8 / 1989, Monash University, Australia). The value of Ks comes from this reference. ALSO NOTE: since we are dealing with a single material model, and since you entered properties with the first (default) material selected, you do not have to Assign these properties to the model. Slide automatically assigns the default properties (i.e. the properties of the first material in the Define Hydraulic Properties dialog) for you. We are now finished with the groundwater modeling, and can proceed to run the groundwater analysis.
Slide v.6.0
Tutorial Manual
7-17
Compute (groundwater)
Before you analyze your model, save it as a file called gw1.slim. (Slide model files have a .slim filename extension).
Slide v.6.0
Tutorial Manual
7-18
Interpret (groundwater)
We can now view the results of the groundwater analysis in the Slide Interpret program.
Slide v.6.0
Tutorial Manual
7-19
Water Table
You will notice on the plot, a pink line which is displayed on the model. This line highlights the location of the Pressure Head = 0 contour boundary. By definition, a Water Table is defined by the location of the Pressure Head = 0 contour boundary. Therefore, for a slope model such as this, this line represents the position of the Water Table (phreatic surface) determined from the finite element analysis. The display of the Water Table can be turned on or off using the toolbar shortcut, the Display Options dialog, or the right-click shortcut (rightclick ON the Water Table and select Hide Water Table). As an optional exercise, you can compare: the Water Table determined from the groundwater analysis (in this tutorial) with the Water Table which we entered in Tutorial 5, for the slope stability analysis.
You will see that the location of the Water Table is very similar in both models. Also notice that the contours of Pressure Head, above the Water Table, have negative values. The negative pressure head calculated above the water table is commonly referred to as the matric suction in the unsaturated zone. This is discussed later in the tutorial. To view contours of other data (Total Head, Pressure, or Discharge Velocity), simply use the mouse to select from the drop-list in the toolbar. Select Total Head contours from the drop-list.
Slide v.6.0
Tutorial Manual
7-20
Flow Vectors
Right click the mouse and select Display Options. Select the Groundwater tab. Toggle ON the Flow Vectors option. Toggle OFF all of the Boundary Condition options. Select Done. (Flow Vectors and other Display Options can also be toggled on or off with shortcut buttons in the toolbar.)
Slide v.6.0
Tutorial Manual
7-21
Flow Lines
Select Total Head contours again. We can also add Flow Lines to the plot. Flow lines can be added individually, with the Add Flow Line option. Or multiple flow lines can be automatically generated with the Add Multiple Flow Lines option.
The generation of the flow lines may take a few seconds. Your screen should look as follows.
Slide v.6.0
Tutorial Manual
7-22
Slide v.6.0
Tutorial Manual
7-23
Iso-Lines
An iso-line is a line of constant contour value, displayed on a contour plot. As we discussed earlier, the pink line which is displayed on the model, represents the Water Table determined by the groundwater analysis. By definition, the Water Table represents an iso-line of zero pressure head. Lets verify that the displayed Water Table does in fact represent a line of zero pressure (P = 0 iso-line), by adding an iso-line to the plot. First, make sure you select Pressure Head contours.
The dialog will display the exact value (Pressure Head) of the location at which you clicked. It may not be exactly zero, so enter zero in the dialog, and select the Add button. An Iso-Line of zero pressure head will be added to the model. It overlaps the displayed Water Table exactly, verifying that the Water Table is the P = 0 line. (If you hide the display of the Water Table, you will see the iso-line in the same location).
Slide v.6.0
Tutorial Manual
7-24
Queries
Another useful feature of the Interpret program, is the ability to create a Query, to obtain detailed analysis results. For groundwater results, a Query allows the user to add a line or polyline, anywhere on the model contours. The Query can then be used to graph values of the contoured data along the Query line or polyline. Lets demonstrate this now. We will create a query which consists of a single vertical line segment, from the vertex at the crest of the slope, to the bottom of the external boundary.
4. Enter a value of 20 in the edit box. Select OK. 5. The query will be created, as you will see by the vertical line segment, and the display of interpolated values at the 20 points along the line segment. 6. Zoom in to the query, so that you can read the values. 7. Now we can graph this data with the Graph Material Queries option (in the Groundwater > Query menu or the toolbar). 8. A shortcut to graph data for a single query, is to right click the mouse ON the Query line. Do this now, and select Graph Data from the popup menu. A graph will be immediately generated, as shown in Figure 7-8.
Slide v.6.0
Tutorial Manual
7-25
3. Now hover the mouse over the contours on the model. As you move the mouse over the contours, the exact interpolated contour value and x-y location, is displayed in a popup data tip. 4. This is a convenient, interactive and graphical way of examining contour values at any point on the model.
Slide v.6.0
Tutorial Manual
7-26
That wraps up our demonstration of viewing groundwater analysis results with the Slide Interpret program. Now, lets get back to the slope stability analysis, so that we can run the analysis using the pore pressures generated by the groundwater analysis.
Slide v.6.0
Tutorial Manual
7-27
Model
In the Slide Model program, change the Analysis Mode from Groundwater to Slope Stability, by selecting from the toolbar, or the tab at the lower left of the screen.
Notice that the modeling entities which are applicable to the slope stability analysis are once again displayed (i.e. Search Grid, Slope Limits, etc). The modeling entities which are applicable to the groundwater analysis (e.g. mesh and boundary conditions), are now hidden.
Compute
Now lets run the slope stability analysis.
Slide v.6.0
Tutorial Manual
7-28
Interpret
To view the results of the analysis:
Figure 7-10: Slope stability analysis using pore pressures from groundwater analysis.
As you can see, after you have performed BOTH: A finite element groundwater analysis, and A slope stability analysis
with Slide, the results of BOTH analyses can be displayed simultaneously in the Slide Interpret program. You can view the contours of groundwater analysis results on the model and also view the slope stability results (slip surfaces etc). Notice that TWO legends are displayed one legend is for the slope stability (safety factor) results, and the other legend is for the groundwater contours. Since we have already viewed and discussed the groundwater analysis results, lets turn OFF the groundwater display options, so that we will only view the slope stability analysis results.
Slide v.6.0
Tutorial Manual
7-29
Use the Mode option in the View menu to quickly toggle the display of Slope Stability or Groundwater results.
TIP: a convenient shortcut for quickly toggling the display of ALL groundwater (or slope stability) results ON or OFF is to use the Mode options in the View menu. However, we still want to display the Water Table from the groundwater analysis, so we will do the following. 1. Right-click and select Display Options. Select the Groundwater tab, and clear ALL of the checkboxes for the groundwater display options, EXCEPT the FEA Water Table and Ponded Water (Hatch) checkboxes. Select Done. 2. Right-click and select Contour Options (groundwater). Set the contour mode = Materials. Select Done. 3. Right-click on the groundwater Legend. Select Hide Legend from the popup menu. Now we are only viewing the slope stability results in Slide Interpret. Lets now compare the results with the Tutorial 5 results.
Slide v.6.0
Tutorial Manual
7-30
Slide v.6.0
Tutorial Manual
7-31
Now lets compare the pore pressures along the Global Minimum slip surface for each file. This way we can directly compare the pore pressures calculated from the pore pressure grid (Tutorial 5) and the groundwater analysis (Tutorial 7). A quick shortcut for graphing data for a slip surface, is to right click on the slip surface (NOTE: you can click on the surface, or on the radial lines joining the slip center to the endpoints of the slip surface). In each view: 1. Right click on the Global Minimum, and select Add Query and Graph from the popup menu. 2. In the Graph Slice Data dialog, select Pore Pressure from the Primary Data list, and select the Create Plot button. 3. Now minimize the model views of each file, and select the Tile Vertically option. You should see the following graphs.
Slide v.6.0
Tutorial Manual
7-32
(The influence of MATRIC SUCTION on the safety factor, is discussed in the last part of this tutorial). To better compare the positive pore pressures, calculated below the water table, do the following: 1. Right click on the pore pressure graph, for the groundwater analysis file, and select Chart Properties from the popup menu. 2. Enter a Minimum Vertical Axis value = 0. Select OK. Your view should now appear as follows:
Slide v.6.0
Tutorial Manual
7-33
tan b
where: = matric suction (POSITIVE value), and b = unsaturated shear strength angle. (NOTE: by convention, the term matric suction implies the POSITIVE, or absolute value of the negative pore pressures calculated in the unsaturated zone). To demonstrate this, try the following: 1. Return to the Slide Model program. 2. Select the groundwater analysis file (GW1.slim). 3. Select Define Material Properties. Notice that, for a groundwater analysis file in Slide, you may specify an Unsaturated Shear Strength Angle. 4. By default, the Unsaturated Shear Strength Angle = 0. This means that matric suction in the unsaturated zone, WILL NOT have any effect on the shear strength or safety factor. 5. However, if you enter a non-zero value for Unsaturated Shear Strength, then slip surfaces which pass through a material in the unsaturated zone, will have INCREASED shear strength, and safety factor. 6. For example, enter an Unsaturated Shear Strength Angle = 15 degrees. Select OK. 7. Now re-run the slope stability analysis. 8. Examine the Global Minimum slip surface. It should now have a safety factor = 1.556 (Bishop analysis).
Slide v.6.0
Tutorial Manual
7-34
9. Specifying an Unsaturated Shear Strength Angle, has increased the Global Minimum safety factor. The Unsaturated Shear Strength Angle is usually not a well known quantity. To obtain an appreciation of its importance, a parametric study can be carried out, in which the Unsaturated Shear Strength angle is varied between 0 and the Friction Angle of the material. This can easily be done using the Sensitivity Analysis option in Slide. See the Slide Help system or Tutorial 09 for more information about Sensitivity Analysis. Unsaturated shear strength can, in some cases, be a critical factor in a slope stability analysis. It has been observed, in some cases, that slopes which are near critical equilibrium (safety factor just over 1), would not be stable without including the effect of matric suction on the shear strength. That concludes this introduction to finite element groundwater analysis using Slide. Further groundwater examples are discussed in the next section.
Slide v.6.0
Tutorial Manual
7-35
Slide v.6.0
Tutorial Manual