Please wait a minute...
Frontiers of Structural and Civil Engineering

ISSN 2095-2430

ISSN 2095-2449(Online)

CN 10-1023/X

Postal Subscription Code 80-968

2018 Impact Factor: 1.272

Front. Struct. Civ. Eng.    2020, Vol. 14 Issue (4) : 998-1011    https://doi.org/10.1007/s11709-020-0621-8
RESEARCH ARTICLE
Three-dimensional finite difference analysis of shallow sprayed concrete tunnels crossing a reverse fault or a normal fault: A parametric study
Masoud RANJBARNIA1(), Milad ZAHERI1, Daniel DIAS2,3
1. Department of Geotechnical Engineering, Faculty of Civil Engineering, University of Tabriz, Tabriz, Iran
2. School of Automotive and Transportation Engineering, Hefei University of Technology, Hefei 230009, China
3. Université Grenoble Alpes, Grenoble F-38000, France
 Download: PDF(4334 KB)   HTML
 Export: BibTeX | EndNote | Reference Manager | ProCite | RefWorks
Abstract

Urban tunnels crossing faults are always at the risk of severe damages. In this paper, the effects of a reverse and a normal fault movement on a transversely crossing shallow shotcreted tunnel are investigated by 3D finite difference analysis. After verifying the accuracy of the numerical simulation predictions with the centrifuge physical model results, a parametric study is then conducted. That is, the effects of various parameters such as the sprayed concrete thickness, the geo-mechanical properties of soil, the tunnel depth, and the fault plane dip angle are studied on the displacements of the ground surface and the tunnel structure, and on the plastic strains of the soil mass around tunnel. The results of each case of reverse and normal faulting are independently discussed and then compared with each other. It is obtained that deeper tunnels show greater displacements for both types of faulting.

Keywords urban tunnel      sprayed concrete      reverse fault      normal fault      finite difference analysis     
Corresponding Author(s): Masoud RANJBARNIA   
Just Accepted Date: 07 May 2020   Online First Date: 10 July 2020    Issue Date: 27 August 2020
 Cite this article:   
Masoud RANJBARNIA,Milad ZAHERI,Daniel DIAS. Three-dimensional finite difference analysis of shallow sprayed concrete tunnels crossing a reverse fault or a normal fault: A parametric study[J]. Front. Struct. Civ. Eng., 2020, 14(4): 998-1011.
 URL:  
https://academic.hep.com.cn/fsce/EN/10.1007/s11709-020-0621-8
https://academic.hep.com.cn/fsce/EN/Y2020/V14/I4/998
Fig.1  Hanging wall and foot wall: (a) in the reverse fault; (b) in the normal fault.
Fig.2  Definition of W parameter.
parameter unit value
elastic modulus of soil MPa 20
friction angle of soil - 37
soil density kg/m3 1630
Poisson’s ratio of soil - 0.3
cohesion of soil kN/m2 0
segment density kg/m3 2600
Poisson’s ratio of segment - 0.28
elastic modulus of segment GPa 20
Tab.1  Properties of the soil and of the tunnel in the centrifuge test level [16]
fault type fault angle height of overburden of soil (m) the length of the tunnel affected by normal fault movement (m) maximum displacement of tunnel (m)
reverse 60° 4.42 - 0.550
reverse 75° 4.42 - 0.256
normal 60° 4.42 9.44 2.500
normal 75° 5.90 11.5 2.500
Tab.2  The displacement of tunnel obtained by the centrifuge physical test model level [16]
Fig.3  Effect of mesh sizes on the tunnel deflection values (the first and second numbers in the legend of figure represent the number of meshes in the Y and Z directions, respectively).
Fig.4  Perspective view of the model and its meshes: (a) general model; (b) half of model along the tunnel axis.
Fig.5  (a) Friction and (b) dilation angle parameters versus plastic shear strain used for the dense soil.
parameter unit dense soil loose soil
elastic modulus MPa 34 15
peak friction angle - 34 30
peak dilation angle - 4 a 0
residual friction angle - 30 30
residual dilation angle - 0 a 0
density kg/m3 1800 b 1500 b
Poisson’s ratio - 0.3 0.3
cohesion kN/m2 0 0
Tab.3  Properties of the soil [51,52]
parameter unit value
density kg/m3 2600
Poisson’s ratio - 0.28
thickness m 0.3
elastic modulus GPa 34
Tab.4  Properties of the shotcrete
Fig.6  Effect of the tunnel shotcrete thickness on the deformed mesh due to the reverse faulting (the properties of the dense soil are mentioned in Table 3). (a) Shotcrete thickness= 0.30 m; (b) shotcrete thickness= 0.35 m.
Fig.7  Effect of the tunnel shotcrete thickness on the deformed mesh due to the normal faulting (the properties of the dense soil are mentioned in Table 3). (a) Shotcrete thickness= 0.30 m; (b) shotcrete thickness= 0.35 m.
Fig.8  Effect of shotcrete thickness on the settlement profile along the tunnel crossing a reverse fault or a normal fault.
Fig.9  Effect of the shotcrete thickness on the maximum vertical displacement of a tunnel crossing a reverse fault or a normal fault.
peak friction angle
(°)
peak dilation angle
(°)
soil density (kg/m3) elastic modulus of soil (MPa) maximum vertical tunnel displacement for reverse faulting (m) maximum vertical tunnel displacement for normal faulting (m)
34 4 1800 34 1.14 2.16
36 6 1800 34 1.17 2.16
38 8 1800 34 1.19 2.16
30 - 1500 15 0.62 2.18
28 - 1500 15 0.59 2.18
24 - 1500 15 0.54 2.18
Tab.5  Effect of the soil geo-mechanical properties on the maximum tunnel vertical displacement in the reverse and normal faulting
Fig.10  Effect of the soil mechanical properties on the deformed mesh due to the reverse faulting. (a) φ= 38°, ψ= 8°, soil density= 1800 kg/m3, elastic modulus of soil=24 MPa; (b) φ= 30°, ψ= 0°, soil density= 1500 kg/m3, elastic modulus of soil=15 MPa.
Fig.11  Effect of the soil mechanical properties on the deformed mesh due to the normal faulting. a) φ= 38°, ψ= 8°, soil density= 1800 kg/m3, elastic modulus of soil=24 MPa; (b) φ= 30°, ψ= 0°, soil density= 1500 kg/m3, elastic modulus of soil=15 MPa.
Fig.12  Effect of the soil geo-mechanical properties on the settlement profile along the tunnel crossing a reverse fault or a normal fault.
Fig.13  Effect of the tunnel depth on the deformed mesh due to the reverse faulting. (a) Tunnel depth= 7.4 m; (b) tunnel depth= 14 m.
Fig.14  Effect of the tunnel depth on the deformed mesh due to the normal faulting. (a) Tunnel depth= 7.4 m; (b) tunnel depth= 14 m.
Fig.15  Effect of the tunnel depth on the settlement profile along the tunnel crossing a reverse fault or a normal fault.
Fig.16  Effect of the tunnel depth on the maximum displacement of a tunnel crossing a reverse fault or a normal fault.
Fig.17  Effect of the reverse fault dip angle on the deformed mesh. (a) Fault dip angle= 75°; (b) fault dip angle= 90°.
Fig.18  Effect of the normal fault dip angle on the deformed mesh. (a) Fault dip angle= 75°; (b) fault dip angle= 90°.
Fig.19  Effect of the fault dip angle on the settlement profile along the tunnel crossing a reverse fault or a normal fault.
Fig.20  Effect of the fault dip angle on the maximum displacement of a tunnel crossing a reverse fault or a normal fault.
the cases point of the maximum gradient in the reverse fault (m) Y maximum gradient in the reverse fault point of the maximum gradient in the normal fault (m) Y maximum gradient in the normal fault
the base example case ?8 −0.09 ?−9 0.10
shotcrete thickness= 0.25 m ?8 −0.1? ?−9 0.10
shotcrete thickness= 0.35 m ?8 −0.09 ?−9 0.10
dense sand, φ = 36°, Es = 24 MPa ?8 −0.11 −10 0.11
dense sand, φ = 38°, Es = 24 MPa 10 −0.14 −12 0.12
loose sand, φ = 30°, Es = 15 MPa 10 −0.07 ??11 0.12
loose sand, φ = 28°, Es = 15 MPa ?7 ?−0.065 ??11 0.12
loose sand, φ = 24°, Es = 15 MPa 10 −0.06 ??11 0.12
tunnel depth= 7.4 m −15?? −0.10 −28 0.16
tunnel depth= 14 m ?5 −0.14 ???8 0.18
fault dip angle= 75° ?2 −0.09 ???5 0.09
fault dip angle= 90° −2? −0.09 ???0 0.06
Tab.6  The maximum gradient of settlement and corresponding location for the different cases (reverse and normal faults)
the cases tunnel deformation in the horizontal direction (%) tunnel deformation in the vertical direction (%)
the base example case 1.64 1.66
shotcrete thickness= 0.25 m 2.91 2.96
shotcrete thickness= 0.35 m 0.86 0.91
dense sand, φ = 36°, Es = 24 MPa 1.56 1.59
dense sand, φ = 38°, Es = 24 MPa 1.50 1.54
loose sand, φ = 30°, Es = 15 MPa 1.40 1.44
loose sand, φ = 28°, Es = 15 MPa 1.33 1.35
loose sand, φ = 24°, Es = 15 MPa 1.18 1.18
tunnel depth= 7.4 m 1.74 1.74
tunnel depth= 14 m 1.79 1.88
fault dip angle= 75° 2.20 2.20
fault dip angle= 90° 2.40 2.44
Tab.7  Tunnel deformation in the horizontal and vertical directions in the reverse fault movement
1 M H Baziar, A Nabizadeh, C J Lee, W Y Hung. Centrifuge modeling of interaction between reverse faulting and tunnel. Soil Dynamics and Earthquake Engineering, 2014, 65: 151–164
https://doi.org/10.1016/j.soildyn.2014.04.008
2 M H Baziar, A Nabizadeh, R Mehrabi, C J Lee, W Y Hung. Evaluation of underground tunnel response to reverse fault rupture using numerical approach. Soil Dynamics and Earthquake Engineering, 2016, 83: 1–17
https://doi.org/10.1016/j.soildyn.2015.11.005
3 M Kiani, T Akhlaghi, A Ghalandarzadeh. Experimental modeling of segmental shallow tunnels in alluvial affected by normal faults. Tunnelling and Underground Space Technology, 2016, 51: 108–119
https://doi.org/10.1016/j.tust.2015.10.005
4 M L Lin, C F Chung, F S Jeng, T C Yao. The deformation of overburden soil induced by thrust faulting and its impact on underground tunnels. Engineering Geology, 2007, 92(3–4): 110–132
https://doi.org/10.1016/j.enggeo.2007.03.008
5 Z Wang, Z Zhang, B Gao. Seismic behavior of the tunnel across active fault. In: The 15th World Conference on Earthquake Engineering. Lisbon, 2012
6 M L Lin, C F Chung, F S Jeng. Deformation of overburden soil induced by thrust fault slip. Engineering Geology, 2006, 88(1–2): 70–89
https://doi.org/10.1016/j.enggeo.2006.08.004
7 D Loukidis, G D Bouckovalas, A G Papadimitriou. Analysis of fault rupture propagation through uniform soil cover. Soil Dynamics and Earthquake Engineering, 2009, 29(11–12): 1389–1404
https://doi.org/10.1016/j.soildyn.2009.04.003
8 M Mortazavi Zanjani, A Soroush. Numerical modelling of fault rupture propagation through layered sands. European Journal of Environmental and Civil Engineering, 2017, 23(9): 1139–1155
9 M Hazeghian, A Soroush. DEM simulations to study the effects of the ground surface geometry on dip-slip faulting through granular soils. European Journal of Environmental and Civil Engineering, 2020, 24(7): 861–879
https://doi.org/10.1080/19648189.2018.1428225
10 I Anastasopoulos, A Callerio, M F Bransby, M C R Davies, A E Nahas, E Faccioli, G Gazetas, A Masella, R Paolucci, A Pecker, E Rossignol. Numerical analyses of fault-foundation interaction. Bulletin of Earthquake Engineering, 2008, 6(4): 645–675
https://doi.org/10.1007/s10518-008-9078-1
11 I Anastasopoulos, G Gazetas. Foundation-structure systems over a rupturing normal fault: Part II. Analysis of the Kocaeli case histories. Bulletin of Earthquake Engineering, 2007, 5(3): 277–301
https://doi.org/10.1007/s10518-007-9030-9
12 M Bransby, M Davies, A El Nahas, S Nagaoka. Centrifuge modelling of reverse fault-foundation interaction. Bulletin of Earthquake Engineering, 2008, 6(4): 607–628
https://doi.org/10.1007/s10518-008-9080-7
13 M Bransby, M Davies, A E Nahas. Centrifuge modelling of normal fault-foundation interaction. Bulletin of Earthquake Engineering, 2008, 6(4): 585–605
https://doi.org/10.1007/s10518-008-9079-0
14 C W W Ng, M A Soomro, Y Hong. Three-dimensional centrifuge modelling of pile group responses to side-by-side twin tunnelling. Tunnelling and Underground Space Technology, 2014, 43: 350–361
https://doi.org/10.1016/j.tust.2014.05.002
15 M A Soomro, C W W Ng, K Liu, N A Memon. Pile responses to side-by-side twin tunnelling in stiff clay: Effects of different tunnel depths relative to pile. Computers and Geotechnics, 2017, 84: 101–116
https://doi.org/10.1016/j.compgeo.2016.11.011
16 M Kiani. Effects of Surface Fault Rupture on Shallow Segmental Soil Tunnels-Centrifuge Modeling. Tabriz: University of Tabriz, 2016 (in Persian)
17 N A Do, D Dias. A comparison of 2D and 3D numerical simulations of tunnelling in soft soils. Environmental Earth Sciences, 2017, 76(3): 102
https://doi.org/10.1007/s12665-017-6425-z
18 N A Do, D Dias, P Oreste. 2D seismic numerical analysis of segmental tunnel lining behaviour. Bulletin of the New Zealand Society for Earthquake Engineering, 2014, 47(3): 206–216
https://doi.org/10.5459/bnzsee.47.3.206-216
19 N A Do, D Dias, P Oreste. Numerical investigation of segmental tunnel linings-comparison between the hyperstatic reaction method and a 3D numerical model. Geomechanics and Engineering, 2018, 14(3): 293–299
20 N A Do, D Dias, P Oreste, I Djeran-Maigre. 2D numerical investigations of twin tunnel interaction. Geomechanics & Engineering, 2014, 6(3): 263–275
21 N A Do, D Dias, P Oreste, I Djeran-Maigre. The behaviour of the segmental tunnel lining studied by the hyperstatic reaction method. European Journal of Environmental and Civil Engineering. 2014, 18(4): 489–510
https://doi.org/https://10.1080/19648189.2013.872583
22 N A Do, D Dias, P Oreste. 3D numerical investigation on the interaction between mechanized twin tunnels in soft ground. Environmental Earth Sciences, 2015, 73(5): 2101–2113
https://doi.org/10.1007/s12665-014-3561-6
23 N A Do, D Dias, P Oreste, I Djeran-Maigre. 2D numerical investigation of segmental tunnel lining behavior. Tunnelling and Underground Space Technology, 2013, 37: 115–127
https://doi.org/10.1016/j.tust.2013.03.008
24 N A Do, D Dias, P Oreste, I Djeran-Maigre. 2D tunnel numerical investigation: The influence of the simplified excavation method on tunnel behaviour. Geotechnical and Geological Engineering, 2014, 32(1): 43–58
https://doi.org/10.1007/s10706-013-9690-y
25 N A Do, D Dias, P Oreste, I Djeran-Maigre. Three-dimensional numerical simulation for mechanized tunnelling in soft ground: The influence of the joint pattern. Acta Geotechnica, 2014, 9(4): 673–694
https://doi.org/10.1007/s11440-013-0279-7
26 N A Do, D Dias, P Oreste, I Djeran-Maigre. Three-dimensional numerical simulation of a mechanized twin tunnels in soft ground. Tunnelling and Underground Space Technology, 2014, 42: 40–51
https://doi.org/10.1016/j.tust.2014.02.001
27 N A Do, D Dias, P Oreste, I Djeran-Maigre. Behaviour of segmental tunnel linings under seismic loads studied with the hyperstatic reaction method. Soil Dynamics and Earthquake Engineering, 2015, 79: 108–117
https://doi.org/10.1016/j.soildyn.2015.09.007
28 T Rabczuk, H Ren, X. Zhuang A nonlocal operator method for partial differential equations with application to electromagnetic waveguide problem. Computers, Materials & Continua, 2019, 59(1): 31–35
https://doi.org/https://doi.org/10.32604/cmc.2019.04567
29 H Guo, X Zhuang, T Rabczuk. A deep collocation method for the bending analysis of Kirchhoff plate. Computers, Materials & Continua, 2019, 59(2): 433–456
30 C Anitescu, E Atroshchenko, N Alajlan, T Rabczuk. Artificial neural network methods for the solution of second order boundary value problems. Computers, Materials & Continua, 2019, 59(1): 345–359
31 Itasca. Fast Lagrangian Analysis of Continua in 3-Dimension (FLAC3D V3.1). Itasca Consulting Group, 2007
32 P Areias, T Rabczuk. Steiner-point free edge cutting of tetrahedral meshes with applications in fracture. Finite Elements in Analysis and Design, 2017, 132: 27–41
https://doi.org/10.1016/j.finel.2017.05.001
33 G Liu, X Zhuang, Z Cui. Three-dimensional slope stability analysis using independent cover based numerical manifold and vector method. Engineering Geology, 2017, 225: 83–95
https://doi.org/10.1016/j.enggeo.2017.02.022
34 Y Zhang, R Lackner, M Zeiml, H A Mang. Strong discontinuity embedded approach with standard SOS formulation: Element formulation, energy-based crack-tracking strategy, and validations. Computer Methods in Applied Mechanics and Engineering, 2015, 287: 335–366
https://doi.org/10.1016/j.cma.2015.02.001
35 Y Zhang, X Zhuang, R Lackner. Stability analysis of shotcrete supported crown of NATM tunnels with discontinuity layout optimization. International Journal for Numerical and Analytical Methods in Geomechanics, 2018, 42(11): 1199–1216
https://doi.org/10.1002/nag.2775
36 H Ren, X Zhuang, T Rabczuk. Dual-horizon peridynamics: A stable solution to varying horizons. Computer Methods in Applied Mechanics and Engineering, 2017, 318: 762–782
https://doi.org/10.1016/j.cma.2016.12.031
37 Y Zhang, X Zhuang. Cracking elements: A self-propagating strong discontinuity embedded approach for quasi-brittle fracture. Finite Elements in Analysis and Design, 2018, 144: 84–100
https://doi.org/10.1016/j.finel.2017.10.007
38 T Rabczuk, G Zi, S Bordas, H Nguyen-Xuan. A simple and robust three-dimensional cracking-particle method without enrichment. Computer Methods in Applied Mechanics and Engineering, 2010, 199(37–40): 2437–2455
https://doi.org/10.1016/j.cma.2010.03.031
39 Y Zhang, X Zhuang. Cracking elements method for dynamic brittle fracture. Theoretical and Applied Fracture Mechanics, 2019, 102: 1–9
https://doi.org/10.1016/j.tafmec.2018.09.015
40 Y Zhang. Multi-slicing strategy for the three-dimensional discontinuity layout optimization (3D DLO). International Journal for Numerical and Analytical Methods in Geomechanics, 2017, 41(4): 488–507
https://doi.org/10.1002/nag.2566
41 S W Sloan. Upper bound limit analysis using finite elements and linear programming. International Journal for Numerical and Analytical Methods in Geomechanics, 1989, 13(3): 263–282
https://doi.org/10.1002/nag.1610130304
42 P Areias, J Reinoso, P P Camanho, J César de Sá, T Rabczuk. Effective 2D and 3D crack propagation with local mesh refinement and the screened Poisson equation. Engineering Fracture Mechanics, 2018, 189: 339–360
https://doi.org/10.1016/j.engfracmech.2017.11.017
43 T Rabczuk, T Belytschko. Cracking particles: A simplified meshfree method for arbitrary evolving cracks. International Journal for Numerical Methods in Engineering, 2004, 61(13): 2316–2343
https://doi.org/10.1002/nme.1151
44 H Ren, X Zhuang, Y Cai, T Rabczuk. Dual-horizon peridynamics. International Journal for Numerical Methods in Engineering, 2016, 108(12): 1451–1476
https://doi.org/10.1002/nme.5257
45 S A Martínez-Galván, M P Romo. Assessment of an alternative to deep foundations in compressible clays: The structural cell foundation. Frontiers of Structural and Civil Engineering, 2018, 12(1): 67–80
https://doi.org/10.1007/s11709-017-0399-5
46 X Yang, J Han, R L Parsons, D Leshchinsky. Three-dimensional numerical modeling of single geocell-reinforced sand. Frontiers of Architecture and Civil Engineering in China, 2010, 4(2): 233–240
https://doi.org/10.1007/s11709-010-0020-7
47 J Atkinson. The Mechanics of Soils and Foundations. London: Taylor & Francis, 2007
48 M A Hicks, K Samy. Influence of heterogeneity on undrained clay slope stability. Quarterly Journal of Engineering Geology and Hydrogeology, 2002, 35(1): 41–49
https://doi.org/10.1144/qjegh.35.1.41
49 K M Hamdia, M Silani, X Zhuang, P He, T Rabczuk. Stochastic analysis of the fracture toughness of polymeric nanoparticle composites using polynomial chaos expansions. International Journal of Fracture, 2017, 206(2): 215–227
https://doi.org/10.1007/s10704-017-0210-6
50 N Vu-Bac, T Lahmer, X Zhuang, T Nguyen-Thoi, T Rabczuk. A software framework for probabilistic sensitivity analysis for computationally expensive models. Advances in Engineering Software, 2016, 100: 19–31
https://doi.org/10.1016/j.advengsoft.2016.06.005
51 B M Das, K Sobhan. Principles of Geotechnical Engineering. California: Thomson Learning, 2013
52 J M R Ortiz, J S Gesta, C O Mazo. Infrastructure Application Course. Madrid: Publishing Office of the Official Institute of Architects, 1986 (in Spanish)
[1] Mehdi SABAGH, Abbas GHALANDARZADEH. Centrifuge experiments for shallow tunnels at active reverse fault intersection[J]. Front. Struct. Civ. Eng., 2020, 14(3): 731-745.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed