03 3D Tunnel Connection
03 3D Tunnel Connection
03 3D Tunnel Connection
Chapter 3
3D Tunnel Connection
Workflow
: Geometry > Point & Curve > Rectangle 7
Chapter 3
3D Tunnel Connection
3D Tunnel Connection
1.1 Learning Purpose
Section 1
Overview The purpose of this analysis is to see the effect of tunnel construction on surrounding facilities and to
review the stability of support structures in advance. It is possible to choose the solution types depend on
whether it's 2D analysis or 3D analysis. If you are processing 2D analysis, you have to search for the
method to consider 3D ground behavior due to tunnel excavation. Sometimes it is very difficult to process
the analysis in 2D, for example junctions which connect tunnel branches, ventilating holes and vertical
shafts together need to be modeled in 3D. In this case you should process 3D analysis to analyze the
ground behavior.
Modeling tunnel
connection
Tunnel modeling
Boolean geometry
Generate and manage mesh sets
Setting construction stages
Results Plotting graphs and tables for construction stage analysis results
This model has various layers of strata with vertical shaft, main tunnel and passageway. Since this model is
symmetry, generate half of the entire model. Then divide it properly and process construction stage
analysis. Define construction stages to Vertical shaft excavation -> Passageway installation -> Main tunnel
excavation.
Modeling diagram
Analysis Setting
[Unit : kN, m]
Table. Ground material Name Soft rock Weathered rock Weathering soil
Material Isotropic Isotropic Isotropic
Model type Mohr-Coulomb Mohr-Coulomb Mohr-Coulomb
General
Elastic Modulus(E) 2.0E+06 5.0E+05 5.0E+04
Poissons Ratio(v) 0.25 0.30 0.30
Unit Weight(r) 25 23 18
Ko 1.0 0.7 0.5
Porous
Unit Weight(Saturated) 25 23 18
Initial Void Ratio(eo) 0.5 0.5 0.5
Drained Parameters Drained Drained Drained
Non-Linear
Cohesion(C) 200 20 20
Friction Angle 35 33 30
[Unit : kN, m]
Table. Structure Name Conc Panel Shotcrete rockbolt
material Material Isotropic Isotropic Isotropic
Model type Elastic Elastic Elastic
Elastic Modulus(E) 2.0E+07 1.5E+07 2.10E+08
Poissons Ratio(v) 0.2 0.2 0.3
Unit Weight(r) 25 24 78.5
Define Structure
material
Properties represent physical attributes of the meshes and will be assigned to mesh sets during mesh
generation. While defining ground and structure properties, firstly choose the material to be used. And for
structure properties, structure types and cross-section shapes (cross-section stiffness) should be further
defined.
Use 'Plate' element for continuous wall in 3D model. And for rockbolt which is 1D element, use 'Embedded
Truss' element. Such as 'Truss' element, 'Embedded Truss' element only bears axial force, but there is a
convenience to use it because you don't need to create shared nodes with the ground. Since it has to be
located inside of the ground, it is mainly used to model rockbolts, nails, and anchors etc.
[Unit : KN, m]
Table. Ground property Name Soft rock Weathered rock Weathering soil
Type 3D 3D 3D
Material Soft rock Weathered rock Weathering soil
The properties for each structure property are defined as below. After setting the section shape, the
section rigidly is automatically calculated.
[Unit : kN, m]
Structure property Name Conc Panel Main_Passageway S/C Vertical shaft S/C Rockbolt
Type 2D 2D 2D 1D
Model type Shell Shell Shell Embedded Truss
Material C/P S/C S/C Rockbolt
Spacing - -
Section shape - Solid Round
Section Thickness TH=0.3 TH=0.15 TH=0.2 D=0.025
Define structure
property-Conc Panel
Define structure
property-
Main/Passageway S/C
Define structure
property-Vertical shaft
S/C
Define structure
property (rockbolt)
Define structure
property(Section shape)
This tutorial focus on tunnel geometry modeling, generating and managing mesh set through Boolean
Section 4 operation, setting construction stages, and analyzing results. You can start the tutorial by opening the start
Modeling file in which basic materials and properties have already been predefined.
In GTS NX, 2 types of coordinate system are used: Global coordinate system (GCS) and Work plane
coordinate system (WCS).
GCS is a fixed coordinate system indicated in red(X axis), green(Y axis), and blue (Z axis) colors at the right
bottom corner on the screen.
WCS is a coordinate system which locates at center of the working window and moves with the work plane.
So if the work plane changes, WCS will also change. 3D absolute coordinates are necessary to locate
geometry in the space, but practically relative coordinates are commonly used, for example to indicate
length. You can process modeling by inputting 2D coordinates (XY plane in WCS) after moving the work
plane to proper location.
Keep in mind that in case of inputting direction for extruding geometry or defining load/boundary
conditions, it always follows Global coordinate system.
Global coordinate
system(GCS) and
Workplane coordinate
system(WCS)
Select tunnel type to '3 Center Circle', and section type to 'Full'.
Input 5.5, 90, 9, 25 at R1, A1, R2, A2.
Check [Include Rock Bolts] option. And enter '13','4' in number, length of rock bolts.
Input '1.5' at Arrangement > Tangential pitch.
Check [Make wire] option. Name geometry set name to [Main tunnel].
Click [OK].
Tunnel section(Main)
Generated tunnel
section shape
To prevent the numerical error when generating tunnel shape at the origin, create tunnel shape at the
random location and then transfer to the origin.
Double click Work Tree > Work Plane > 'YZ (1, 0, 0)'.
Click [Normal] in the view toolbar.
Select tunnel type to '3 Center Circle', and section type to 'Full'.
Input '4, 90, 9, 20 at R1, A1, R2, and A2.
Check [Include Rock Bolts] option. And enter '10','4' in number, length of rock bolts.
Input '1.5' at Arrangement > Tangential pitch.
Check [Make wire] option. Name geometry set name to [Passageway].
Click [OK].
Modeling tunnel
section(Passageway)
Transfer passageway
Transferred tunnel
section
You can select the object shape at the generated geometry set in the Work Tree.
Embedded solid(Main
tunnel/Vertical shaft)
To primarily filter certain type of object which will be selected, choose the object type from [Selection
Filter].
[Selection Filter] consists of object types registered in the software depending on current working status
and command type. The user can select specific object types from this composition.
Double click Work Tree > Work Plane > YZ (1, 0, 0).
Click [Normal] in the view toolbar.
Check [Make Face] option.
Select the geometry set to Passageway.
Make a rectangle which is little bit bigger than passageway tunnel section.
Click [OK].
Divide solid(Passageway)
Double click Work Tree > Work Plane > YX (0, 0, 1).
Click [Normal] in the view toolbar.
Check [Make Face] option.
Select the geometry set to Vertical shaft.
Make a rectangle which is little bit bigger than vertical shaft tunnel section.
Click [OK].
Divide solid(Vertical
shaft)
Divide solid(Ground)
Double click Work Tree > Work Plane > XY (0, 0,1).
Input the start point '25, 0'. And select 'Vertical shaft' geometry set.
Select [Apply]. And put end location to '4, 0.
Click [OK].
Generated
rockbolt(Vertical shaft)
As the mesh size gets smaller, the result comes out closer to the true answer. But small mesh size will
also lead to longer analysis time. So only generate small and dense meshes around the part (around
tunnel) where we want to see the results.
Generate mesh
(Area A~C)
C
Hide all the meshes and geometries except for solids of Main tunnel, Passageway, Vertical shaft
in the Work Tree.
Select [Front] in the view Toolbar.
In the Mesh > Element > Extract, select type to 'Face'.
Check [Skip Duplicated Faces] and [Register Based-on Owner Shape] option.
Select the property to '5:Main_connection S/C'.
Click [OK].
Uncheck to hide Main tunnel, Passageway, Vertical solids in the Work Tree > Geometry.
Select [Right] in the View Toolbar.
Refer to image below, drag & select elements of area A, B, C, D.
Click [OK].
Delete mesh
Generated mesh(Plate)
Generated mesh(Main
tunnel rockbolt)
Delete
mesh(Passageway
rockbolt)
Generate
mesh(Passageway
rockbolt)
Select rockbolt edges (6) of the tunnel registered in Work Tree > geometry > Vertical shaft > Line.
Set [Size Method] to Size and input 2.
Check if [Attribute ID] set as 7.
Select [Preview] button to see the locations of generated nodes.
Click [OK].
Select the entire 'Main tunnel' in the Work Tree > Mesh.
Set Global Rectangular for coordinate in [Sorting Order]. And select 1st to 'Y'.
Set the [Order] to 'Ascending' and enter Name to 'Main tunnel#'. Input '1' in the [Starting suffix
Number].
Click [Apply].
Select the entire 'Passageway' in the Work Tree > Mesh.
Set Global Rectangular for coordinate in [Sorting Order]. And select 1st to 'X'.
Set the [Order] to 'Descending' and enter Name to 'Passageway#'. Input '1' in the [Starting suffix
Number].
Click [Apply].
Select the entire 'Vertical shaft' in the Work Tree > Mesh.
Set Global Rectangular for coordinate in [Sorting Order]. And select 1st to 'Z'.
Set the [Order] to 'Descending' and enter Name to 'Vertical shaft#'. Input '1' in the [Starting suffix
Number].
Click [Apply].
In the same way rename Shotcrete and Rockbolt. Rename the ground mesh with [F2] key to 'Soft
rock', 'Weathered rock', 'Weathering soil'.
Through View Toolbar > Display mode(Mesh)( ) > Property color or Material color, you can easily
check out assigned property.
Set the [Name] to Gravity-1 and enter 'Self weight' as the [Load set] name.
Put -1 in Gz this defines the gravity accelerations direction.
Click [OK].
In the [Auto] tab, set name and boundary set name to Constraint-1, Constraint.
In 3D model analysis, displacement constraints are applied in x direction for left/right, y direction for
front/back, xyz direction for the bottom part. GTS NX automatically recognizes and sets the boundary
conditions for the model.
In GTS NX, it is possible to create several construction stage sets to perform different analysis.
Table. Stage definition Set Name A/R End Postfix Start Stage Stage Inc.
wizard Vertical shaft#- R 1 1 1
Vertical shaft S/C#- A 1 2 1
Vertical shaft R/B#- A 1 2 1
Passageway#- R 1 8 1
Passageway S/C#- A 1 9 1
Passageway R/B#- A 1 9 1
Main tunnel#- R 1 12 1
Main tunnel S/C#- A 1 13 1
Main tunnel R/B#- A 1 13 1
Apply Assignment
Rules
Select Main tunnel#-, Passageway#-, Vertical shaft#-, Soft rock, weathered rock, weathering soil.
Drag selected item to [I.S.] of [Element, Boundary, and Load Activation Status].
Drag boundary condition, load to [I.S.] as following image.
Define construction
stage(Initial step)
Analysis control
Selected the analysis cases that you want to perform and click [OK].
During the analysis, you can check the calculation process in real-time. Messages such as
whether the results converge or not, warnings and errors can be checked through [Output
Window].
The model is automatically saved before being analyzed. The information about analysis is saved
in *.OUT file in the same place where the model file is saved.
When analyzing a model, the results can be controlled by Analysis > Analysis Case > General >
Output Control option. And if you set the output option as [Binary and Text], the results of nodes
and elements after the analysis are also plotted in *.OUT file.
You can verify displacement, stress, member force etc. in the Result Tree after the analysis. All the results
Section 6 can be plotted as contour, table, and graph. In this tutorial, the main result items which need to be checked
Results are listed below.
Verify by 'Displacement' of the work tree. TX, TY, TZ signify displacements in directions of X, Y, Z. Since
gravity direction of this model is Z d, tunnel crown displacement is shown in TZ TRANSLATION, and tunnel
internal displacement can be verified in TX TRANSLATION, TY TRANSLATION.
'(V)' refers to the result item which can represent both in contour and vector formats at the same time. In
GTS NX, it is possible to show contour/vector simultaneously for displacement and stress results.
Select the last stage (S27) to check the result in the Work Tree > Results > 3D Tunnel Connection
Analysis, select Displacement > TZ TRANSLATION (V).
Select [Front] in View Toolbar. And verify the results in the front view.
Crown displacement
(Isometric1)
Crown displacement
(Front view)
Click Result > Advance > Probe to see the value of selected node. Also by this function you can
find the locations of maximum, minimum results.
Probe
Simulation bar at the bottom of the work window helps you to view the evolution of results in each
construction stage.
Select the stage (S27) in the Work Tree > Results > 3D Tunnel Connection Analysis, and select
Displacement > TX TRANSLATION (V), TX TRANSLATION to see the displacement results in X, Y
directions.
X direction
displacement
Y direction
displacement
The Segment stress can be verified by 'Solid Stresses' in the Results Tree. S-XX, S-YY, S-ZZ represent the
stresses in each direction. You can see the maximum principal stresses in 'S-PRINCIPAL A (V)', and
minimum principal stress in 'S-PRINCIPAL C (V)'.
Select the stage(S27) in the Work Tree > Results > 3D Tunnel Connection Analysis, and select the
Solid Stresses > S-PRINCIPAL A(V), S-PRINCIPAL C(V) to verify the maximum principal stresses
and minimum principal stresses.
Maximum principal
stress
Minimum principal
stress
Deselect Result > General > Contour, select vector. In the property window at the left bottom of
the work window, check [Only Free Face] option of vector. Select [Apply].
Property window(Vector)
Minimum principal
stress(Vector)
Verify stresses of Shotcrete and axis force of Rockbolt. Since shotcrete is shell element, results about X, Y,
Z axis of element axis are plotted. Shell is a plate element which has a thickness, so results plotted in TOP,
MID, and BOTTOM part of the plate.
Select a stage (S27) to check the result in the Work Tree > Results > 3D Tunnel Connection
Analysis, select Shell Element Stresses > S-XX TOP.
X direction stress of
Shotcrete
Exclude no result
entity
Select a stage (S27) to check the result in the Work Tree > Results > 3D Tunnel Connection
Analysis, select Truss Element Forces > AXIAL FORCE.
In case of 3D model, you can slice the model with specific plane and check the results at the cutting plane.
In GTS NX, by slicing the model with cutting planes you can easily check the results of the specific points
located at inner part of the model.
Select a stage (S27) to check the result in the Work Tree > Results > 3D Tunnel Connection
Analysis, select Displacement > TOTAL TRANSLATION(V).
Select [Clipping Plane] in Advanced View Toolbar. In [Define Plane] window, select plane
direction to 'X' and direction to '0m'. Click [Add] button to create 'plane1'.
Sliced by plane1
Sliced by plane1&2
Select [3D->2D] in Result > Advanced > Others. [3D->2D] wizard is for verifying results by tagging
the result when 3D model is sliced by specific plane. Check [Show Points] to see the points in the
slice plane.
Probe
Probe at slice plane