Logo image
Back
Time-domain PML formulation for modeling viscoelastic waves with Rayleigh-type damping in an unbounded domain: Theory and application in ABAQUS
Journal article   Open access   Peer reviewed

Time-domain PML formulation for modeling viscoelastic waves with Rayleigh-type damping in an unbounded domain: Theory and application in ABAQUS

Masoud Khazaei Poul and Aspasia Zerva
Finite elements in analysis and design, v 152, pp 1-16
01 Dec 2018
DOI: https://doi.org/10.1016/j.finel.2018.08.004
url
https://doi.org/10.1016/j.finel.2018.08.004View
Accepted (AM)Maybe Open Access (Publisher Bronze),  Open

Abstract

Mathematics Mathematics, Applied Mechanics Physical Sciences Science & Technology Technology
The Perfectly Matched Layer (PML) method is a powerful approach to absorb outgoing waves from a finite computational domain in various media, such as elastic, poroelastic, anisotropic, and viscoelastic, by means of layers of artificial material placed at the finite domain boundaries. However, its use in most finite-element method (FEM) codes, which take into account viscous damping by employing only the Rayleigh and Lindsday damping formulation, is inadequate. This paper introduces a viscoelastic Perfectly Matched Layer (PML) formulation with Rayleigh-type damping for finite-element time-domain analyses. The PML formulations in the frequency- and time-domain are derived using a two-parameter complex coordinate stretching function. The displacement field is the only unknown variable in the formulation, so that the approach can be readily implemented in general finite-element codes. The weak form formulation in the time-domain is implemented in ABAQUS/Standard by merging it with a user-defined element (UEL) subroutine in Fortran90. The UEL is developed for 2D plane-strain (PE4ML) and plane-stress (PS4ML) four-node, linear, isoparametric, quadrilateral elements. The Hilber-Hughes-Taylor implicit time integration scheme is used to evaluate the unknown displacements, velocities, and accelerations. The validity and efficiency of the viscoelastic PML formulation and the UEL subroutine are examined with two numerical examples, namely a single-layer and a multi-layered soil profile with different damping ratios and properties. The results also highlight the poor performance of elastic PML layers, when they are placed at the boundaries of a viscoelastic interior domain, because, as the damping ratio increases, the reflection at the interface of the interior domain and the elastic PML domain increases. Long-time stability of the model is also examined, and no instabilities are observed.

Metrics

16 Record Views
27 citations in Web of Science
27 citations in Scopus

Details

InCites Highlights

Data related to this publication, from InCites Benchmarking & Analytics tool:

Web of Science research areas
Mathematics, Applied
Mechanics
Logo image