An end-to-end 3D seismic simulation of underground structures due to point dislocation source by using an FK-FEM hybrid approach

Abstract

Based on the domain reduction idea and artificial boundary substructure method, this paper proposes an FK-FEM hybrid approach by integrating the advantages of FK and FEM (i.e., FK can efficiently generate high-frequency three translational motion, while FEM has rich elements types and constitutive models). An advantage of this approach is that it realizes the entire process simulation from point dislocation source to underground structure. Compared with the plane wave field input method, the FK-FEM hybrid approach can reflect the spatial variability of seismic motion and the influence of source and propagation path. This approach can provide an effective solution for seismic analysis of underground structures under scenario of earthquake in regions where strong earthquakes may occur but are not recorded, especially when active faults, crustal, and soil parameters are available. Taking Daikai subway station as an example, the seismic response of the underground structure is simulated after verifying the correctness of the approach and the effects of crustal velocity structure and source parameters on the seismic response of Daikai station are discussed. In this example, the influence of velocity structure on the maximum interlayer displacement angle of underground structure is 96.5% and the change of source parameters can lead to the change of structural failure direction.

Graphical abstract

Keywords

source-to-structure simulation / FK-FEM hybrid approach / underground structures / point dislocation source

References

[1] Dong Y H, Peng F L, Guo T F. Quantitative assessment method on urban vitality of metro-led underground space based on multi-source data: A case study of Shanghai Inner Ring area. Tunnelling and Underground Space Technology, 2021, 116: 104108

[2] Xie H, Zhang Y, Chen Y, Peng Q, Liao Z, Zhu J. A case study of development and utilization of urban underground space in Shenzhen and the Guangdong-Hong Kong-Macao Greater Bay Area. Tunnelling and Underground Space Technology, 2021, 107: 103651

[3] Guo D, Zhu X, Xie J, Zhang C, Zhao Z, Zheng L. Undergraduate program for urban underground space engineering in China: Exploration and practice. Tunnelling and Underground Space Technology, 2021, 116: 104084

[4] Zaini F, Hussin K, Raid M M. Legal considerations for urban underground space development in Malaysia. Underground Space, 2017, 2(4): 234–245

[5] Yang J, Zhuang H, Wang W, Zhou Z, Chen G. Seismic performance and effective isolation of a large multilayered underground subway station. Soil Dynamics and Earthquake Engineering, 2021, 142: 106560

[6] Iida H, Hiroto T, Yoshida N, Iwafuji M. Damage to Daikai subway station. Soils and Foundations, 1996, 36(Special): 283–300

[7] Samata S, Ohuchi H, Matsuda T. A study of the damage of subway structures during the 1995 Hanshin-Awaji earthquake. Cement and Concrete Composites, 1997, 19(3): 223–239

[8] Yamaguchi A, Mori T, Kazama M, Yoshida N. Liquefaction in Tohoku district during the 2011 off the Pacific Coast of Tohoku Earthquake. Soil and Foundation, 2012, 52(5): 811–829

[9] Xu Z, Du X, Xu C, Hao H, Bi K, Jiang J. Numerical research on seismic response characteristics of shallow buried rectangular underground structure. Soil Dynamics and Earthquake Engineering, 2019, 116: 242–252

[10] Li W, Chen Q. Effect of vertical ground motions and overburden depth on the seismic responses of large underground structures. Engineering Structures, 2020, 205: 110073

[11] Hashash Y M A, Hook J J, Schmidt B, I-Chiang Yao J. Seismic design and analysis of underground structures. Tunnelling and Underground Space Technology, 2001, 16(4): 247–293

[12] Kuesel T R. Earthquake design criteria for subways. Journal of the Structural Division, 1969, 95(6): 1213–1231

[13] Liu J, Wang W, Dasgupta G. Pushover analysis of underground structures: method and application. Science China. Technological Sciences, 2014, 57(2): 423–437

[14] Penzien J. Seismically induced racking of tunnel linings. Earthquake Engineering & Structural Dynamics, 2000, 29(5): 683–691

[15] TateishiA. A study on loading method of seismic deformation method. Proceedings of the Civil Society, 1992, 1992(441): 157–166 (in Japanese)

[16] Tateishi A. A study on seismic analysis methods in the cross section of underground structures using static finite element method. Structural Engineering/Earthquake Engineering, 2005, 22(1): 41s–54s

[17] WangJ N. Seismic Design of Tunnels: A Simple State-of-the-Art Design Approach. New York: Brinckerhoff Parsons Inc., 1993

[18] Zhao M, Gao Z, Du X, Wang J, Zhong Z. Response spectrum method for seismic soil-structure interaction analysis of underground structure. Bulletin of Earthquake Engineering, 2019, 17(9): 5339–5363

[19] Xu C, Zhang Z, Li Y, Du X. Validation of a numerical model based on dynamic centrifuge tests and studies on the earthquake damage mechanism of underground frame structures. Tunnelling and Underground Space Technology, 2020, 104: 103538

[20] Zhu T, Wang R, Zhang J M. Evaluation of various seismic response analysis methods for underground structures in saturated sand. Tunnelling and Underground Space Technology, 2021, 110: 103803

[21] IwatateTKobayashi YKusuHRinK. Investigation and shaking table tests of subway structures of the Hyogoken-Nanbu earthquake. In: 12th World Conference on Earthquake Engineering. Auckland: New Zealand Society for Earthquake Engineering Upper Hutt, 2020

[22] Zhang Z, Li Y, Xu C, Du X, Dou P, Yan G Y. Study on seismic failure mechanism of shallow buried underground frame structures based on dynamic centrifuge tests. Soil Dynamics and Earthquake Engineering, 2021, 150: 106938

[23] Trifonov O V, Cherniy V P. A semi-analytical approach to a nonlinear stress–strain analysis of buried steel pipelines crossing active faults. Soil Dynamics and Earthquake Engineering, 2010, 30(11): 1298–1308

[24] Wang X, Chen J, Xiao M. Seismic responses of an underground powerhouse structure subjected to oblique incidence SV and P waves. Soil Dynamics and Earthquake Engineering, 2019, 119: 130–143

[25] Wang X, Xiong Q, Zhou H, Chen J, Xiao M. Three-dimensional (3D) dynamic finite element modeling of the effects of a geological fault on the seismic response of underground caverns. Tunnelling and Underground Space Technology, 2020, 96: 103210

[26] Jiang J, Nggar H M, Xu C, Zhong Z, Du X. Effect of ground motion characteristics on seismic fragility of subway station. Soil Dynamics and Earthquake Engineering, 2021, 143: 106618

[27] Zhu Z, Tang Y, Ba Z, Wang K, Gong W. Seismic analysis of high-speed railway irregular bridge–track system considering V-shaped canyon effect. Railway Engineering Science, 2022, 30(1): 57–70

[28] Zhang L, Wang J, Xu Y, He C, Zhang C. A procedure for 3D seismic simulation from rupture to structures by coupling SEM and FEM. Bulletin of the Seismological Society of America, 2020, 110(3): 1134–1148

[29] He C H, Wang J T, Zhang C H, Jin F. Simulation of broadband seismic ground motions at dam canyons by using a deterministic numerical approach. Soil Dynamics and Earthquake Engineering, 2015, 76: 136–144

[30] Mccallen D, Petersson N, Rodgers A, Pitarka A, Miah M, Petrone F, Sjogreen B, Abrahamson N, Tang H. EQSIM—A multidisciplinary framework for fault-to-structure earthquake simulations on exascale computers part I: Computational models and workflow. Earthquake Spectra, 2021, 37(2): 707–735

[31] Yu H T, Yang Y S, Yuan Y, Duan K P, Gu Q. A comparison between vibration and wave methods in seismic analysis of underground structures. China Earthquake Engineering Journal, 2019, 41(04): 845–852

[32] Bielak J, Loukakis K, Hisada Y, Yoshimura C. Domain reduction method for three-dimensional earthquake modeling in localized regions, Part I: Theory. Bulletin of the Seismological Society of America, 2003, 93(2): 817–824

[33] BielakJ. Reply to “Comment on ‘Domain reduction method for three-dimensional earthquake modeling in localized regions, Part I: Theory,’ by J. Bielak, K. Loukakis, Y. Hisada, and C. Yoshimura, and ‘Part II: Verification and applications,’ by C. Yoshimura, J. Bielak, Y. Hisada, and A. Fernandez,” by E. Faccioli, M. Vanini, R. Paolucci, and M. Stupazzini. Bulletin of the Seismological Society of America, 2005, 95(2): 770−773

[34] Yoshimura C. Domain reduction method for three-dimensional earthquake modeling in localized regions, Part II: Verification and applications. Bulletin of the Seismological Society of America, 2003, 93(2): 825–841

[35] Li W, Assimaki D. Simulating soil stiffness degradation in transient site response predictions. Soil Dynamics and Earthquake Engineering, 2010, 30(5): 299–309

[36] Crempien J G F, Archuleta R J. UCSB method for simulation of broadband ground motion from kinematic earthquake sources. Seismological Research Letters, 2015, 86(1): 61–67

[37] Brun M, De Martin F, Richart N. Hybrid asynchronous SEM/FEM co-simulation for seismic nonlinear analysis of concrete gravity dams. Computers & Structures, 2021, 245: 106459

[38] Casadei F, Gabellini E, Fotia G, Maggio F, Quarteroni A. A mortar spectral/finite element method for complex 2D and 3D elastodynamic problems. Computer Methods in Applied Mechanics and Engineering, 2002, 191(45): 5119–5148

[39] Liu J, Tan H, Bao X, Wang D, Li S. Seismic wave input method for three-dimensional soil-structure dynamic interaction analysis based on the substructure of artificial boundaries. Earthquake Engineering and Engineering Vibration, 2019, 18(4): 747–758

[40] Liu J, Bao X, Wang D, Wang P. Seismic response analysis of the reef-seawater system under incident SV wave. Ocean Engineering, 2019, 180: 199–210

[41] AkiKRichardsP G. Quantitative Seismology. 2nd ed. Sausalito: University Science Books, 2002

[42] Ba Z, Sang Q, Wu M, Liang J. The revised direct stiffness matrix method for seismogram synthesis due to dislocations: From crustal to geotechnical scale. Geophysical Journal International, 2021, 227(1): 717–734

[43] Liu J, Du Y, Du X, Wang Z, Wu J. 3D viscous-spring artificial boundary in time domain. Earthquake Engineering and Engineering Vibration, 2006, 5(1): 93–102

[44] Fukuyama E, Ishida M, Dreger D, Kawai H. Automated seismic moment tensor determination by using on-line broadband seismic waveforms. Journal of the Seismological Society of Japan, 1998, 51(1): 149–156

Acknowledgements

This work was supported by Open Foundation of National Engineering Laboratory for High Speed Railway Construction (No. HSR202006) and National Natural Science Foundation of China (Grant Nos. 52178495, 52078498).

RIGHTS & PERMISSIONS

2022 Higher Education Press