Estimation of Surface Fault Displacement by High Performance Computing

摘要

Reliable estimation of surface fault displacements is crucial to the safety of nuclear power plant facilities. It is necessary to develop a numerical method for the estimation. In the study, we develop a finite element method in which the following two functions are implemented: (1) a symplectic time integration of an explicit scheme to properly conserve the energy of the system; and (2) rigorously formulated joint elements of high order. The infinite element method is enhanced with parallel computing capability. We apply the developed method to solve simple three-dimensional models of faults embedded in a rock mass. It includes a comparison of results from quasi-static and dynamic simulations and investigation of the sensitivity of results to the shear stiffness on faults. In the study, we propose capacity computing with a quasi-static simulation for uncertainty quantification.