Efficiency of Higher Order Finite Elements for the Analysis of Seismic Wave Propagation

The analysis of wave propagation problems in linear damped media must take into account both propagation features and attenuation process. To perform accurate numerical investigations by the finite differences or finite element method, one must consider a specific problem known as the numerical dispersion of waves. Numerical dispersion may increase the numerical error during the propagation process as the wave velocity (phase and group) depends on the features of the numerical model. In this paper, the numerical modelling of wave propagation by the finite element method is thus analyzed and dis-cussed for linear constitutive laws. Numerical dispersion is analyzed herein through 1D computations investigating the accuracy of higher order 15-node finite elements towards numerical dispersion. Concerning the numerical analy-sis of wave attenuation, a rheological interpretation of the classical Rayleigh assumption has for instance been previously proposed in this journal

Seismic site effects in a deep alluvial basin: numerical analysis by the boundary element method

The main purpose of the paper is the numerical analysis of seismic site effects in Caracas (Venezuela). The analysis is performed considering the boundary element method in the frequency domain. A numerical model including a part of the local topography is considered, it involves a deep alluvial deposit on an elastic bedrock. The amplification of seismic motion (SH-waves, weak motion) is analyzed in terms of level, occurring frequency and location. In this specific site of Caracas, the amplification factor is found to reach a maximum value of 25. Site effects occur in the thickest part of the basin for low frequencies (below 1.0 Hz) and in two intermediate thinner areas for frequencies above 1.0 Hz. The influence of both incidence and shear wave velocities is also investigated. A comparison with microtremor recordings is presented afterwards. The results of both numerical and experimental approaches are in good agreement in terms of fundamental frequencies in the deepest part of the basin. The boundary element method appears to be a reliable and efficient approach for the analysis of seismic site effects

Numerical analysis of seismic wave amplification in Nice (France) and comparisons with experiments

The analysis of site effects is very important since the amplification of seismic motion in some specific areas can be very strong. In this paper, the site considered is located in the centre of Nice on the French Riviera. Site effects are investigated considering a numerical approach (Boundary Element Method) and are compared to experimental results (weak motion and microtremors). The investigation of seismic site effects through numerical approaches is interesting because it shows the dependency of the amplification level on such parameters as wave velocity in surface soil layers, velocity contrast with deep layers, seismic wave type, incidence and damping. In this specific area of Nice, a one-dimensional (1D) analytical analysis of amplification does not give a satisfactory estimation of the maximum reached levels. A boundary element model is then proposed considering different wave types (SH, P, SV) as the seismic loading. The alluvial basin is successively assumed as an isotropic linear elastic medium and an isotropic linear viscoelastic solid (standard solid). The thickness of the surface layer, its mechanical properties, its general shape as well as the seismic wave type involved have a great influence on the maximum amplification and the frequency for which it occurs. For real earthquakes, the numerical results are in very good agreement with experimental measurements for each motion component. Two-dimensional basin effects are found to be very strong and are well reproduced numerically

A Simple Multi-Directional Absorbing Layer Method to Simulate Elastic Wave Propagation in Unbounded Domains

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)

Modelling strong seismic ground motion: three-dimensional loading path versus wavefield polarization

Seismic 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

A Simple and Efficient Regularization Method for 3D BEM: Application to Frequency-Domain Elastodynamics

An efficient and easy-to-implement method is proposed to regularize integral equations in the 3D boundary element method (BEM). The method takes advantage of an assumed three-noded triangle discretization of the boundary surfaces. The method is based on the derivation of analytical expressions of singular integrals. To demonstrate the accuracy of the method, three elastodynamic problems are numerically worked out in the frequency domain: a cavity under harmonic pressure, diffraction of a plane wave by a spherical cavity, and amplification of seismic waves in a semispherical alluvial basin (the second one is also investigated in the time domain). The numerical results are compared to semi-analytical solutions; a close agreement is found for all problems, showing the accuracy of the proposed method

La motivation des élèves de formation professionnelle à poursuivre leurs études dans un programme technique harmonisé

Cette étude vise à mieux comprendre la motivation des élèves de formation professionnelle à poursuivre ou non leurs études dans un programme technique harmonisé en prenant en compte deux modalités d’harmonisation différentes. Des entrevues réalisées auprès des élèves de neuf centres de formation professionnelle ainsi qu’auprès des directions et des enseignants des centres et des cégeps concernés nous apprennent que ceux qui envisagent de poursuivre manifestent souvent de l’intérêt pour les études collégiales avant leur inscription en formation professionnelle, que leurs buts de formation au cégep sont clairs et qu’ils ont reçu une information personnalisée au sujet du programme technique harmonisé.The aim of this study is to better understand the motivation of students in professional training to continue or not their studies in a technical training programme while taking into account two different harmonization modalities. Interviews were conducted with students in nine professional training centres, and the directors and teachers of both centres and colleges (Cegep) involved. The findings show that those who intend to continue often show interest in collegial studies even before their registration in professional training, that they have clear objectives regarding college training, and that they received personalized information regarding a harmonized technical programme.El presente estudio tiene por objetivo entender mejor la motivación de los alumnos de formación profesional para proseguir, o no, con sus estudios en un programa técnico armonizado tomando en cuenta dos modalidades de armonización diferentes. Entrevistas realizadas con alumnos de nueve centros de formación profesional así como direcciones y docentes de estos centros y escuelas preparatorias nos indican que los que planean proseguir manifiestan con frecuencia interés para los estudios colegiales antes de su inscripción en formación profesional, que sus objetivos de formación en la preparatoria son claros y que se les ha proporcionado una información personalizada acerca del programa técnico armonizado