The numerical analysis of elastic wave propagation in unbounded media may be difficult due to spurious waves reflected at the model artificial boundaries. This point is critical for the analysis of wave propagation in heterogeneous or layered solids. Various techniques such as Absorbing Boundary Conditions, infinite elements or Absorbing Boundary Layers (e.g. Perfectly Matched Layers) lead to an important reduction of such spurious reflections. In this paper, a simple absorbing layer method is proposed: it is based on a Rayleigh/Caughey damping formulation which is often already available in existing Finite Element softwares. The principle of the Caughey Absorbing Layer Method is first presented (including a rheological interpretation). The efficiency of the method is then shown through 1D Finite Element simulations considering homogeneous and heterogeneous damping in the absorbing layer. 2D models are considered afterwards to assess the efficiency of the absorbing layer method for various wave types and incidences. A comparison with the PML method is first performed for pure P‐waves and the method is shown to be reliable in a more complex 2D case involving various wave types and incidences. It may thus be used for various types of problems involving elastic waves (e.g. machine vibrations, seismic waves, etc.). Copyright © 2010 John Wiley & Sons, Ltd.

No abstract

[1] An integrated numerical and geophysical approach was used to investigate the role of local seismic amplification in the reactivation of the Salcito landslide (Italy) after the Molise earthquake of 31 October 2002. Numerical stress-strain analysis with FLAC 5.0 FDM software, performed in dynamic configuration, showed that the 1 Hz frequency was consistent with landslide instability conditions. The Fourier spectrum of the triggering earthquake showed two main peaks occurring at 1 and 2 Hz, respectively, which could be related to local effects. The analysis of H/V spectral ratios of ambient noise records obtained in the landslide area, evidenced amplification effects in the 1-3 Hz range. On the basis of an engineering geology model, numerical modeling of both 1-D and 2-D seismic wave propagation was conducted using linear and nonlinear solutions. The simulation outputs showed (1) a 1 Hz amplification ascribable to 2-D effects due to a landslide mass lying within a basin-like geological structure; (2) a double-peak amplification at about 2 and 3 Hz, respectively, ascribable to 1-D resonance of the landslide mass; and (3) 1-D plus lateral wave effects within the landslide mass in the 2.5-3.5 Hz frequency range. These findings suggest that local amplification of ground motion by the Salcito landslide mass may generate a self-excitation process responsible for its reactivation.

International audienceSeismic waves due to strong earthquakes propagating in surficial soil layers may both reduce soil stiffness and increase the energy dissipation into the soil. To investigate seismic wave amplification in such cases, past studies have been devoted to one-directional shear wave propagation in a soil column (1D-propagation) considering one motion component only (1C-polarization). Three independent purely 1C computations may be performed ('1D-1C' approach) and directly superimposed in the case of weak motions (linear behaviour). This research aims at studying local site effects by considering seismic wave propagation in a 1-D soil profile accounting for the influence of the 3-D loading path and non-linear hysteretic behaviour of the soil. In the proposed '1D-3C' approach, the three components (3C-polarization) of the incident wave are simultaneously propagated into a horizontal multilayered soil. A 3-D non-linear constitutive relation for the soil is implemented in the framework of the Finite Element Method in the time domain. The complex rheology of soils is modelled by mean of a multisurface cyclic plasticity model of the Masing-Prandtl-Ishlinskii-Iwan type. The great advantage of this choice is that the only data needed to describe the model is the modulus reduction curve. A parametric study is carried out to characterize the changes in the seismic motion of the surficial layers due to both incident wavefield properties and soil non-linearities. The numerical simulations show a seismic response depending on several parameters such as polarization of seismic waves, material elastic and dynamic properties, as well as on the impedance contrast between layers and frequency content and oscillatory character of the input motion. The 3-D loading path due to the 3C-polarization leads to multi-axial stress interaction that reduces soil strength and increases non-linear effects. The non-linear behaviour of the soil may have beneficial or detrimental effects on the seismic response at the free surface, depending on the energy dissipation rate. Free surface time histories, stress-strain hysteresis loops and in-depth profiles of octahedral stress and strain are estimated for each soil column. The combination of three separate 1D-1C non-linear analyses is compared to the proposed 1D-3C approach, evidencing the influence of the 3C-polarization and the 3-D loading path on strong seismic motions

No abstract