Examples • Example Applications
Swelling of a Fully Wetted Slope (FLAC2D)
Problem Statement
Note
This project reproduces Example 19 from FLAC 8.1. The project file for this example is available to be viewed/run in FLAC2D.[1] The project’s main data files are shown at the end of this example.
This example illustrates the estimation of volume change (swelling deformation) of a fill slope subjected to post-construction wetting. Swelling deformations are calculated using the swell constitutive model. The example is based upon an actual field case of a fill slope presented by Noorany et al. (1999).
The problem conditions correspond to a fully wetted fill slope. The slope is composed of three layers of fill material with different elastic, plastic, and wetting properties. The slope height is approximately 22 m, and the slope angle is roughly 24 degrees. The slope crest is at elevation 304.8 m. The top layer (soil 3) extends to elevation 293.0 m, the second layer (soil 4) extends to elevation 284.0 m, and the third layer (soil 5) extends approximately 3.0m to a slightly angled rock foundation. The slope face contains small 1.8-m wide benches at elevations 298.0 m and 284.0 m.
Relative horizontal movements between two points at two different sections on the fill slope crest were monitored when the slope was fully wetted. The monitoring points were approximately 26 m apart. The average swelling movement between points was found to be 103 mm.
Modeling Procedure
Model Geometry
The model zones are created in a sketch using the Slope Wizard tool. After opening the slope wizard, fill in the form as shown in Figure 1. Note that you don’t need to fill in the slope angles. Slope angle is automatically calculated in the dialog based on the other dimensions. Make sure you click in the slope angle boxes after filling in rise and run for each level and ensure that the numbers match those shown below. Also make sure to check the box at the bottom to Mirror layout and uncheck the box to Create internal edges.
The resulting slope is shown in Figure 2.
Next, delete the bottom-right point (at x=89, y = 279). Draw a line to join the bottom-left point and bottom-right point.
Next we need to define the material layers. Draw horizontal lines at y=284 and y=293. The model should now appear as shown in Figure 3.
Next we mesh the model. Set the zone size to 1 and mesh all polygons. The model should appear as in Figure 4.
Finally, we can define the material layers. Click on the top layer and set the group name to Soil3
. Set the group name of the middle layer to Soil4
and the bottom to Soil5
. Save the project and the model state. You may also wish to save the commands used to create the model. Do do this, click on the State Record tab at the bottom, right-click where you see the commands, and select .
Model Properties
The slope fill is modeled as Swell material. The constitutive model and properties can be set in the model pane, or by using commands. The assigned properties are shown in Table 1. Be careful to remain consistent with units. Stiffnesses are shown as kPa in Table 1 but Pa units are actually used in the example.
Soil |
3 |
4 |
5 |
---|---|---|---|
mass density (\(kg/m^3\)) |
1728.3 |
1575.4 |
1978.5 |
bulk modulus (kPa) |
15400 |
32500 |
29300 |
shear modulus (kPa) |
7120 |
15000 |
22000 |
friction angle (deg) |
30 |
30 |
40 |
cohesion (kPa) |
4.8 |
4.8 |
1.44 |
swelling type |
logarithmic |
logarithmic |
linear |
a1 |
1.533 |
0.694 |
0.0004 |
c1 |
-0.0187 |
-0.0211 |
0 |
a3 |
0.436 |
0.468 |
0.00047 |
c3 |
-0.0215 |
-0.023 |
0.001 |
m1 |
0.023 |
0.028 |
0.0 |
m3 |
0.036 |
0.036 |
0.0 |
Initially, the slope is considered to be un-wetted. The swell model requires an initial stress for each zone, so first the model must be solved with a mohr-coulomb material. We initially set the density, moduli and strength properties as shown in Table 1. The model is then solved.
The model is then changed to swell
and all properties are assigned. For this example, principal
swelling directions match the global coordinate system (angle = 0). Atmospheric pressure is set to
101.3 kPa, and swelling stresses are introduced over 20,000 steps.
Displacements and velocities for all gridpoints in the model are reset to zero, and the calculation is continued to introduce swelling.
Simulation Results
Figure 5 shows the vectors of the displacements induced by full wetting of the slope. Figure 6 shows the history plots of x-displacement and y-displacement at the crest edge. The maximum horizontal movement (to the right) obtained at the edge of the crest is ~131 mm. Swell (vertical displacement) at the same location is ~84 mm.
It is possible to compare results of the FLAC2D simulation with field data (see Noorany et al., 1999). The relative horizontal movement obtained from the FLAC2D analysis between two points on top of the slope (x=13.5 and the crest, x=39.6) is found to be 103 mm. This shows good agreement with field data. The calculated movement is compared to the measured value in Figure 7.
Reference
Noorany, I., S. Frydman and C. Detournay. “Prediction of Soil Slope Deformation Due toWetting”, in FLAC and Numerical Modeling in Geomechanics, Proceedings of the First International FLAC Symposium (Minneapolis, September 1999), pp. 101-107, C. Detournay and R. Hart, Eds. Rotterdam: Balkema (1999).
Data Files
initial.dat
model new
program call 'sketch'
; assign Mohr-coulomb properties. Swell properties are assigned later
; units are Pa
zone cmodel assign mohr-coulomb
zone property density 1728.3 bulk 15.4e6 shear 7.12e6 friction 30 ...
cohesion 4800 range group 'Soil3'
zone property density 1575.4 bulk 32.5e6 shear 15.0e6 friction 30 ...
cohesion 4800 range group 'Soil4'
zone property density 1978.5 bulk 29.3e6 shear 22.0e6 friction 40 ...
cohesion 1440 range group 'Soil5'
; identify boundaries
zone face skin
; apply boundary conditions
zone face apply velocity-x 0 range group 'East' or 'West'
zone face apply velocity 0 0 range group 'Bottom'
; turn on gravity and solve
model gravity 9.81
model large-strain off
model solve elastic
model save 'initial'
swell.dat
model restore 'initial'
; reset displacements and velocities
zone gridpoint initialize displacement 0 0
zone gridpoint initialize velocity 0 0
; assign swell properties.
zone cmodel assign swell
; MC properties
zone property density 1728.3 bulk 15.4e6 shear 7.12e6 friction 30 ...
cohesion 4800 range group 'Soil3'
zone property density 1575.4 bulk 32.5e6 shear 15.0e6 friction 30 ...
cohesion 4800 range group 'Soil4'
zone property density 1978.5 bulk 29.3e6 shear 22.0e6 friction 40 ...
cohesion 1440 range group 'Soil5'
; swell properties
; these for all layers
zone property normal (0,1) pressure-reference 101.3e3 number-start 20000
; these are layer-specific
zone property constant-a-1 1.533 constant-c-1 -0.0187 ...
constant-a-3 0.436 constant-c-3 -0.0215 ...
constant-m-1 0.023 constant-m-3 0.036 ...
range group 'Soil3'
zone property constant-a-1 0.694 constant-c-1 -0.0211 ...
constant-a-3 0.468 constant-c-3 -0.023 ...
constant-m-1 0.028 constant-m-3 0.036 ...
range group 'Soil4'
zone property constant-a-1 0.0004 constant-c-1 0 ...
constant-a-3 0.00047 constant-c-3 0.001 ...
flag-swell 2 ...
range group 'Soil5'
; take histories at the crest
zone history name 'disp-x' displacement-x position 40.2,304.8
zone history name 'disp-y' displacement-y position 40.2,304.8
; get relative displacement between two locations
fish define xdisp_rel
zone.field.name = 'displacement-x'
; displacements in mm
local xdisp_right = zone.field.get(39.6,304.8)*1000
local xdisp_left = zone.field.get(13.5,304.8)*1000
xdisp_rel = xdisp_right - xdisp_left
xdisp_meas = 102 ; field data
end
fish history name 'FLAC2D' xdisp_rel
fish history name 'Measured' xdisp_meas
model cycle 25000
model save 'swell'
Endnote
⇐ Pile-Supported Highway Embankment | Undrained Cylindrical Cavity Expansion in a Cam-Clay Medium (FLAC3D) ⇒
Was this helpful? ... | Itasca Software © 2024, Itasca | Updated: Sep 15, 2024 |