Seismic response of three-dimensional rockfill dams using the Indirect Boundary Element Method

The Indirect Boundary Element Method (IBEM) is used to compute the seismic response of a three-dimensional rockfill dam model. The IBEM is based on a single layer integral representation of elastic fields in terms of the full-space Green function, or fundamental solution of the equations of dynamic elasticity, and the associated force densities along the boundaries. The method has been applied to simulate the ground motion in several configurations of surface geology. Moreover, the IBEM has been used as benchmark to test other procedures. We compute the seismic response of a three-dimensional rockfill dam model placed within a canyon that constitutes an irregularity on the surface of an elastic half-space. The rockfill is also assumed elastic with hysteretic damping to account for energy dissipation. Various types of incident waves are considered to analyze the physical characteristics of the response: symmetries, amplifications, impulse response and the like. Computations are performed in the frequency domain and lead to time response using Fourier analysis. In the present implementation a symmetrical model is used to test symmetries. The boundaries of each region are discretized into boundary elements whose size depends on the shortest wavelength, typically, six boundary segments per wavelength. Usually, the seismic response of rockfill dams is simulated using either finite elements (FEM) or finite differences (FDM). In most applications, commercial tools that combine features of these methods are used to assess the seismic response of the system for a given motion at the base of model. However, in order to consider realistic excitation of seismic waves with different incidence angles and azimuth we explore the IBEM.