Various software packages are available to conduct one-dimensional (1D) and twodimensional (2D) site response analyses (SRAs). In this article, a finite element program is tested with the purpose of assessing the importance of several aspects on the obtained results and verifying the software. Abaqus (Dassault Systèmes) is compared with the 1D SRA software STRATA (Kottke & Rathje, 2008) for simple 1D models to understand the influence of the boundary conditions; as being a 1D SRA program, STRATA does not require vertical boundary conditions. For Abaqus the subroutine by Nielsen (2006Nielsen ( , 2014 is used to implement free-field boundary conditions.In addition, we test the influence of mesh dimension and Rayleigh damping as well as the importance of buffer-zone width. 2D SRAs with Abaqus and FLAC3D (Itasca Consulting Group, 2012), which is commonly used for geotechnical analyses, are compared as part of the assessment. Similar results are obtained from the two programs but Abaqus is preferred as it is more efficient for linear elastic analyses than FLAC3D, which, on the contrary, performs well for soil presenting strongly non-linear behaviour and effective stress. We demonstrate that reliable results can be achieved, not only for simple uniform sites but also for complex sites with multiple layers and dipping stratigraphy.
KeywordsOne-dimensional, two-dimensional, site response analysis, free-field boundary conditions, damping
1.IntroductionAn important part of geotechnical earthquake engineering is the study of the response of the ground under earthquake excitation by means of site response analysis (SRA). Based on the characteristics of the problem studied, several approaches can be used: from one-dimensional (1D) to threedimensional (3D) analyses and from linear-elastic to fully non-linear soil behaviour. 1D linearelastic SRA is often a good starting point for all such studies, even though it is unrealistic in many cases (e.g. when the site is subject to high-amplitude shaking, for which a non-linear analysis would provide more accurate results, and in the case of sedimentary valleys and basins, where 3D effects are pronounced). A complete SRA should include a consideration of the 3D geometry and nonlinear soil behaviour and the consideration of effective stress, which is important, for example, for the study of soil liquefaction. The seismic waves, under these conditions, can generate an increase in the pore water pressure producing a reduction of stiffness and strength of the soil. However, it is not common to perform such analyses, because the more complex a model is, the more input parameters are needed, which means that a full characterization of the site must be available to the engineer. In addition, 3D and non-linear analyses require long computational times, sophisticated software and much experience in conducting the analyses and in interpreting the results.In this work, we test several approaches, but we will limit ourselves to 1D and 2D SRA and to linear-elastic behaviour with viscous damping. A ...