The Effects of Boundary Conditions on Accuracy and Efficiency of Finite Element Analysis for Fault Dislocation within Homogeneous Elastic Half-Space

Document Type : Original Research

Authors
Arak University
Abstract
Surface deformation of the earth due to earthquake fault dislocation is very important for predicting ground motions. There are many studies on kinematic modeling of earthquake faults in both analytical and numerical methods. However, suitable investigations on improving usefulness and efficiency of those numerical methods are still necessary for relevant researchers. In this paper, displacement fields for free surface of the earth due to fault dislocation in homogeneous elastic half-space have been investigated by finite element method. Boundary conditions have significant effects on the results of finite element method, especially when the domain of the problem has infinite boundaries (half-space). Therefore, appropriate cares should be taken to increase the efficiency and accuracy of this method. In order to achieve a comprehensive study on this topic, boundaries have been modeled by two approaches here. The first approach uses common elements for infinite boundaries, while the second one uses infinite elements for those boundaries. To verify the results, each problem has been examined by several meshes and numerical solutions have been compared to Okada’s analytical solutions. In addition to the effects of the boundary modeling, the discretization effects have been investigated in order to find a suitable approach to reduce computational efforts and to increase the accuracy and efficiency of finite element method. In the modeling process, contact elements have been employed to impose fault dislocation. Three numerical examples have been provided for these finite element analyses. Each example includes four analyses without infinite elements and three analyses with infinite elements which are also compared together. The results show that not only infinite elements are necessary for quasi-static fault dislocation problems, but also they improve the performance of finite element method, so that with finer meshes and smaller dimensions of a domain, analytical solutions can be captured by numerical solutions with suitable accuracy.

Keywords

Subjects


[1] Zwieten G.J., Hanssen R.F., Gutiérrez M.A. 2013 Overview of a range of solution methods for elastic dislocation problems in geophysics. Journal of Geophysical Research: Solid Earth, 118(4), pp. 1721-1732.
[2] Zakian P., Khaji N. 2016 Spectral finite element simulation of seismic wave propagation and fault dislocation in elastic media. Asian Journal of Civil Engineering, 17(8), pp. 1189-1213.
[3] Zakian P., Khaji N. 2019 A stochastic spectral finite element method for solution of faulting-induced wave propagation in materially random continua without explicitly modeled discontinuities. Computational Mechanics, 64(4), pp. 1017–1048.
[4] Okada Y. 1992 Internal deformation due to shear and tensile faults in a half-space. Bulletin of the Seismological Society of America, 82(2), pp. 1018-1040.
[5] Okada Y. 1985 Surface deformation due to shear and tensile faults in a half-space. Bulletin of the Seismological Society of America, 75(4), pp. 1135-1154.
[6] Zakian P., Khaji N. 2016 Development of higher-order stochastic spectral finite element method for uncertainty analysis of 2D continua. Modares Mechanical Engineering, 16(7), pp. 51-60. ‘’(Persian)’’
[7] Melosh H., Raefsky A. 1981 A simple and efficient method for introducing faults into finite element computations. Bulletin of the Seismological Society of America, 71(5), pp. 1391-1400.
[8] Mirashrafi R., Khaji N. 2016 Determination of slippage field of strike-slip faults using inverse solution methods. Journal of the Earth and Space Physics, 42(2), pp. 233-245. ‘’(Persian)’’
[9] Megna A., Barba S., Santini S. 2005 Normal-fault stress and displacement through finite-element analysis. Annals of Geophysics, 48(6), pp. 1009-1016.
[10] Megna A., Barba S., Santini S., 2008 Dragoni M. Effects of geological complexities on coseismic displacement: hints from 2D numerical modelling. Terra Nova, 20(3), pp. 173-179.
[11] Lavecchia G., Castaldo R., Nardis R. and et al. Ground deformation and source geometry of the 24 August 2016 Amatrice earthquake (Central Italy) investigated through analytical and numerical modeling of DInSAR measurements and structural‐geological data. Geophysical Research Letters, 43(24) (2016).
[12] Zakian P., Khaji N., Soltani N. 2017 A Monte Carlo adapted finite element method for dislocation simulation of faults with uncertain geometry. Journal of Earth System Science, 126(7), pp. 105.
[13] Candela S., Mazzoli S., Megna A., Santini S. 2015 Finite element modelling of stress field perturbations and interseismic crustal deformation in the Val d'Agri region southern Apennines, Italy. Tectonophysics, 657, pp. 245-259.
[14] Cattin R., Briole P., Lyon-Caen H., Bernard P., Pinettes P. 1999 Effects of superficial layers on coseismic displacements for a dip-slip fault and geophysical implications. Geophysical Journal International, 137(1), pp. 149-158.
[15] Rezaei m.h., and Khaji N. 2018 Simulation of 2009 L’aquila Earthquake using Specific barrier model. Modares Civil Engineering journal, 18(2), p. 101-112. ‘’(Persian)’’
[16] Nabavi S.T., Alavi S.A., Mohammadi S., Ghassemi M.R. 2018 Mechanical evolution of transpression zones affected by fault interactions: Insights from 3D elasto-plastic finite element models. Journal of Structural Geology, 106, pp. 19-40.
[17] Castaldo R., De Novellis V., Solaro G. and et al. 2017 Finite element modelling of the 2015 Gorkha earthquake through the joint exploitation of DInSAR measurements and geologic-structural information. Tectonophysics, 714, pp. 125-132.
[18] Monshizadeh Naiin A., Seyedi Hosseininia E. 2018 Numerical Investigation of Reverse Faulting Effect on the Response of Buried Pipes. Modares Civil Engineering journal, 18(2), pp. 233-244.
[19] Liu G.R., Quek S.S., Michael E., 2003 The Finite Element Method: A Practical Course. London, Butterworth-Heinemann.