3D Tunnel Simulation Using Core Replacement 1.0 Introduction
3D Tunnel Simulation Using Core Replacement 1.0 Introduction
The goal of this tutorial is to demonstrate how to model the tunnel deformation prior to
support installation using the core replacement (material softening) approach.
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 1/27
25/8/2019 RS2: 3D Tunnel Simulation
3. Build a model that relaxes the boundary to the calculated amount in step 2. Add
the support and determine whether a) the tunnel is stable, b) the tunnel wall
deformation meets the specified requirements, and c) the tunnel lining meets
certain factor of safety requirements. If any of these conditions are not met,
choose a different support system and run the analysis again (Core Replacement
Part 2 Model)
Under the General tab, define the units as “Metric, stress as MPa.”
Under the Stages tab, change the number of stages to 9 (see following figure). Fill
in the stage names as shown below. Close the dialog by clicking OK.
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 2/27
25/8/2019 RS2: 3D Tunnel Simulation
2.2 Geometry
1. Right-click the mouse and select the Circle option from the popup menu.
2. Select the Center and radius option, enter Radius = 4 and enter Number of
Segments = 96 and select OK.
3. Enter the circle center: Enter (0,0) in the prompt line and the circular excavation
will be created.
Select Zoom All (or press the F2 function key) to zoom the excavation to the center of
the view.
Now let’s create the external boundary. In RS2, the external boundary may be
automatically generated, or user-defined. Let’s use one of the ‘automatic’ options:
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 3/27
25/8/2019 RS2: 3D Tunnel Simulation
RS2 Select: Boundaries > Add External Click here to see this page in full context
In the Create External Boundary dialog, use set Boundary Type = Box and Expansion
Factor = 5. Select OK, and the external boundary will be automatically created.
2.3 Materials
For Material 1, under the Strength tab, change the Failure Criterion to
Generalized Hoek-Brown and the Material Type to Plastic.
Now define the strength parameters and the Young’s Modulus using the GSI
calculator. Press the GSI calculator button:
In the GSI calculator dialog, set the uniaxal compressive strength of the
intact rock equal to 50 MPa, the GSI equal to 50 and mi equal to 10. To
compute the rock mass deformation modulus, set the modulus ratio (MR) to
400.
Press the OK button. The material properties dialog should now be updated
with the new strength and modulus values.
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 4/27
25/8/2019 RS2: 3D Tunnel Simulation
Click on the Material 2 tab and change the name to E=3000. Change the Initial
Element Loading to None. In the Stiffness tab, change the Young’s Modulus to
3000 MPa.
Now follow the same procedure and set the Young’s modulus of Materials 3
through 8 to 1000, 250, 100, 50, 20, and 10 MPa respectively. Change the names
to reflect the value of the modulus. Make sure that the Initial Element Loading for
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 5/27
25/8/2019 RS2: 3D Tunnel Simulation
The first material, with modulus 6143 MPa and Generalized Hoek-Brown failure
criterion, is the in-situ rock mass. Materials 2 through 8 will be used inside the
excavation (excavation core). The core material is progressively replaced over several
stages. This replacement, along with the modulus reduction, allows the boundary to
progressively deform. In each of the eight stages, the material inside the excavation is
replaced by a material with zero internal stress (i.e. Initial Element Loading = None)
and with a lower modulus than the proceeding stage. In the final stage, the material
inside the excavation is removed. This process models the advancement of the tunnel
face. Each stage (and corresponding core modulus) represents some distance from the
tunnel face, either in front of or behind the face. The final excavated stage represents
the deformed state far away from the tunnel face, at a distance where the face has no
influence on stresses or displacements. What’s left is determining the correspondence
between core modulus and distance from the tunnel face, specifically, the modulus
sequence that yields the deformation at the support installation distance. The support
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 6/27
25/8/2019 RS2: 3D Tunnel Simulation
installation distance being the distance between the tunnel face and where the support
RS2 Click here to see this page in full context
is installed.
To determine the correspondence between core modulus and distance from the tunnel
face, the relationship between tunnel wall deformation and distance from the tunnel
face must be known. Knowing the relationship between tunnel wall displacement and
distance from the tunnel face and knowing the relationship between core modulus and
tunnel wall displacement, the relationship between core modulus and distance from the
tunnel face can then be determined. Knowing this relationship allows the calculation of
the modulus reduction sequence that gives the tunnel wall displacement prior to
support installation.
1. Make sure the Stage 2 tab, E=3000, is selected (at the bottom left of the view).
4. Change to Stage 3, E=1000, by clicking the stage tab at the bottom of the screen.
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 7/27
25/8/2019 RS2: 3D Tunnel Simulation
23. Select the “Excavate” button at the bottom of the Assign dialog.
24. Click the left mouse button inside the tunnel. The material inside the excavation
should now be removed.
25. Close the Assign dialog by clicking on the X in the upper right corner of the dialog.
Now select Stage 1 – the in-situ condition stage. Turn on the minimum data tips mode
using the following command.
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 8/27
25/8/2019 RS2: 3D Tunnel Simulation
Hover the mouse inside the excavation. After a second, a data tip should appear:
RS2 Click here to see this page in full context
Notice that the data tip shows all the materials inside the excavation as a function of
stage.
2.5 Mesh
In the mesh setup dialog, ensure the Element Type is 6 Noded Triangles.
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 9/27
25/8/2019 RS2: 3D Tunnel Simulation
Click the Discretize button and then the Mesh button. Click OK to close the dialog.
Field Stress determines the initial in-situ stress conditions, prior to excavation.
Enter Sigma 1 = 30, Sigma 3 = 15, Sigma Z = 25, Angle = 0, and select OK.
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 10/27
25/8/2019 RS2: 3D Tunnel Simulation
3.0 Results
RS2 and Discussion Click here to see this page in full context
The maximum stress, Sigma 1 for Stage 1 will be displayed. Notice that there is no
variation of stress and that the stress (30 MPa) is equal to the major in-situ field stress.
This is expected since in the first stage the material inside and outside the tunnel
boundary is the in-situ E=6143 material.
Change the contours to plot Total Displacement using the pull-down menu in the
toolbar. The model for Stage 1 will look like this:
There is essentially no displacement in the first stage. Now click through the stages. An
increase in deformation is visible around the tunnel as the core material is replaced and
softened (modulus reduced).
3.1 Step 1: Computing tunnel deformation before support installation using the
Vlachopoulos and Diederichs method
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 11/27
25/8/2019 RS2: 3D Tunnel Simulation
To compute the tunnel deformation at the point of support installation, use the empirical
RS2 Click here to see this page in full context
relationship developed by Vlachopoulos and Diederichs. To use this method,
information from the finite element analysis is required: a) the maximum tunnel wall
displacement far from the tunnel face, and b) the radius of the plastic zone far from the
tunnel face.
Both values can be computed from a plane strain analysis with zero internal pressure
inside the excavation. In the model constructed in this tutorial, the results from stage 9
are used since the material inside the excavation is completely removed in this stage.
Switch to the last stage (Stage 9). Look at the bottom of the legend. The maximum
displacement for this stage is approximately 0.05 m. This is the value of maximum wall
displacement far from the tunnel face. The location of this displacement is in the roof
and floor of the excavation. The location of this displacement is important since any
comparisons of displacement for various core moduli must be made at the same
location.
Several crosses are visible and represent elements in the finite element analysis that
have failed. Zoom Out so that the entire extent of failed points is visible (see below).
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 12/27
25/8/2019 RS2: 3D Tunnel Simulation
The extent of this failed zone represents the extent of the plastic zone around the
RS2 Click here to see this page in full context
tunnel. To determine the radius of the plastic zone, use either the measuring tool or the
dimensioning tool to measure the distance from the center of the tunnel to the
perimeter of the yielded/plastic zone. This tutorial uses the measuring tool:
Select (0, 0) as the location to measure from. Use the mouse to extend the measuring
line vertically until the edge of the yield zone is reached. Click here to select.
The following plot was created using the Vlachopoulos and Diederichs equations
(Vlachopoulos and Diederichs, 2009). The equations can also be found in the Kersten
Lecture, appendix 1 (Hoek et. al., 2008). Using this plot, it is possible to estimate the
amount of closure prior to support installation if the plastic radius and displacement far
from the tunnel face are known.
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 13/27
25/8/2019 RS2: 3D Tunnel Simulation
For our problem, Rp=8.8m, Rt=4m, X=2m, and umax=0.069m. The Distance from
tunnel face/tunnel radius = 2/4 = 0.5. The Plastic zone radius/tunnel radius = 8.8/4 =
2.2. From the above plot this gives Closure/max closure approximately equal to 0.44.
Therefore the closure equals (0.44)*(0.069) = 0.030 m.
As computed above, the tunnel roof displaces 0.030m before the support is installed.
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 14/27
25/8/2019 RS2: 3D Tunnel Simulation
The next step is to determine the core modulus that yields a displacement of 0.030m in
RS2 Click here to see this page in full context
the roof of the tunnel. It is important to maintain the same location as is used to
determine umax, since the location of maximum displacement can change depending
on the magnitude of the internal pressure. This can be seen in this model as larger core
moduli produce larger displacement in the sidewall while smaller core moduli produce
larger displacements in the roof and floor.
To determine the internal pressure that yields a 0.030m roof displacement, plot the
displacement versus stage for a point on the roof of the excavation.
1. When asked to enter a vertex, type in the value 0, 4 for the location and press
Enter. This is a point on the roof of the excavation.
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 15/27
25/8/2019 RS2: 3D Tunnel Simulation
3. Press the Plot button. The following figure shows the plot generated by the
program. This is a plot of displacement versus stage for a point in the roof of the
tunnel.
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 16/27
25/8/2019 RS2: 3D Tunnel Simulation
Right-click in the plot and choose Sampler > Show Sample. Move the sampler by
moving the mouse with the left mouse button. Move the sampler until the displacement
value on the right side of the plot is equal to 0.030m.
This plot shows that in stage 4, the wall displacement in the roof of the tunnel is
approximately 0.030m. This represents a 3-stage material replacement and reduction
of core modulus from E=6143(insitu), to 3000, 1000 and finally 250 MPa.
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 17/27
25/8/2019 RS2: 3D Tunnel Simulation
Excel will launch with a plot of stage number versus displacement. The plot can be
easily changed to show the stage number data to the core modulus. A sample of the
Excel file for this example is included in the Tutorials folder with the RS2 data files.
The following image shows the convergence-confinement plot in Excel for this example.
The plot shows that modulus reduction to 250MPa yields the tunnel wall displacement
computed above for the point of support installation (0.030m).
Steps 1 and 2 as defined in the Problem section at the beginning of this tutorial have
been completed; lets proceed by defining the support system.
Note: see the note at the end of this tutorial about how to carry out the analysis if the required
modulus value lies between two values in the initial modulus reduction sequence.
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 18/27
25/8/2019 RS2: 3D Tunnel Simulation
We will now use the 9-stage model created above and modify it to create the support
RS2 Click here to see this page in full context
design.
Select: Analysis > Project Settings and select the Stages Tab.
It is important that we keep all the core softening stages up to the stage that represents
support installation. This is because the replacement and softening of the core material
in stages 2 and 3 affect the final displacement result. These stages directly influence
the stress path and displacement of the material around the excavation.
Make sure the Stage 5, Support Installed stage tab is selected. Click the Zoom
Excavation button on the toolbar.
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 19/27
25/8/2019 RS2: 3D Tunnel Simulation
Now define the liner properties. The properties correspond to a 200 mm thick layer of
concrete reinforced with W150X18 I-beams spaced at 2-meter intervals along the
tunnel axis.
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 20/27
25/8/2019 RS2: 3D Tunnel Simulation
4. In the Define Liner Properties dialog, for the Reinforcement, enter a spacing of
2m.
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 21/27
25/8/2019 RS2: 3D Tunnel Simulation
Let’s line the tunnel with the liner defined above. First make sure that Stage 5, the
RS2 Click here to see this page in full context
Support Installed stage, is selected.
1. In the Add Liner dialog, select Tunnel Liner as the Liner Property and set the liner
to install at stage 5. Select OK.
2. Click and hold the left mouse button and drag a selection window which encloses
the entire excavation. Release the left mouse button. Notice that all excavation
line segments are selected.
3. Right-click the mouse and select Done Selection, or just press the Enter key. The
entire tunnel will now be lined, as indicated by the thick blue line segments
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 22/27
25/8/2019 RS2: 3D Tunnel Simulation
Click through the stages. Notice how the color of the liner changes from light blue in
stages 1 through 4 to dark blue in stage 5. This indicates that the liner is being installed
in stage 5.
5.0 Compute
Before analyzing the model, save it as a new file called 3D Tunnel Simulation using
Core Replacement (Part 2).fez (make sure to select Save As and not Save, or it will
overwrite the internal pressure reduction file).
Save the file as 3D Tunnel Simulation using Core Replacement (Part 2).fez.
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 23/27
25/8/2019 RS2: 3D Tunnel Simulation
If any other files are loaded in the Interpret program (i.e. the CoreSoftening.fez file),
close them. Click on the tab at the bottom of the program window associated with the
file and use the File > Close menu option to close the file.
Make sure the Stage 5 tab is selected and zoom in on the excavation.
Support capacity diagrams give the engineer a method for determining the factor of
safety of a reinforced concrete liner. For a given factor of safety, capacity envelopes are
plotted in axial force versus moment space and axial force versus shear force space.
Values of axial force, moment, and shear force for the liner are then compared to the
capacity envelopes. If the computed liner values fall inside an envelope, they have a
factor of safety greater than the envelope value. So, if all the computed liner values fall
inside the design factor of safety capacity envelope, the factor of safety of the liner
exceeds the design factor of safety.
The Support Capacity Plot dialog allows the user to choose the support element (i.e.
liner type), the number of envelopes, and the stages from which the liner data is taken.
Use the spin control to increase the number of envelopes to 3. The dialog should look
like:
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 24/27
25/8/2019 RS2: 3D Tunnel Simulation
Press OK.
The following plot is generated. The dark red lines represent the capacity envelopes for
the 3 factors of safety (1, 1.2, 1.4).
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 25/27
25/8/2019 RS2: 3D Tunnel Simulation
Notice that all the data points fall within the factor of safety=1.4 envelope, on all four
plots. This means that the support system chosen has a factor of safety greater than
1.4.
In this example, the required core modulus, which gives the displacement required at
the point of support installation, happens to be exactly equal to one of the original
modulus values chosen for the initial reduction sequence (i.e. 250 MPa). In general,
this will not be the case. That is, the required core modulus will probably lie between
two of the values chosen for the initial modulus reduction sequence. If this occurs:
2. Then either insert a new stage of core replacement, with the required modulus
value, or simply use the nearest stage with a HIGHER modulus value and lower
the material modulus at this stage to the required value (e.g. if the required
modulus is 350 MPa, but the initial sequence goes from 500 to 250, then change
the 500 to 350).
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 26/27
25/8/2019 RS2: 3D Tunnel Simulation
3. Re-run the analysis and check if the new modulus value does in fact give the
RS2 Click here to see this page in full context
desired displacement at the point of support installation. It should be close. If not,
then repeat steps 1 to 3 until the required modulus value is determined.
7.0 References
Hoek, E., Carranza-Torres, C., Diederichs, M.S. and Corkum, B. (2008). Integration of
geotechnical and structural design in tunnelling – 2008 Kersten Lecture. Proceedings
University of Minnesota 56th Annual Geotechnical Engineering Conference.
Minneapolis, 29 February 2008, 1-53.
https://www.rocscience.com/help/rs2/tutorials/RS2__3D_Tunnel_Simulation.htm 27/27