An unsplit-field viscoelastic complex-frequency-shifted perfectly matched layer for analysis of transient waves in heterogeneous media based on an efficient voxel element method

被引:1
作者
Naeini, Ehsan Fakharzadeh [1 ]
Khaji, Naser [2 ]
Hori, Muneo [3 ]
机构
[1] Tarbiat Modares Univ, Fac Civil & Environm Engn, POB 14115-397, Tehran, Iran
[2] Hiroshima Univ, Grad Sch Adv Sci & Engn, Civil & Environm Engn Program, 1-4-1 Kagamiyama, Higashihiroshima, Hiroshima 7398527, Japan
[3] Japan Agcy Marine Earth Sci & Technol, Res Inst Value Added Informat Generat, Yokohama, Japan
基金
日本学术振兴会;
关键词
Rayleigh damping; Voxel element method; Perfectly matched layer; Transient wave propagation; Viscoelastic; Unsplit-field viscoelastic CFS-PML; DIAGONAL COEFFICIENT MATRICES; ABSORBING BOUNDARY-CONDITION; SOIL DYNAMIC IMPEDANCE; ELASTODYNAMIC PROBLEMS; SEMIANALYTICAL METHOD; GRAZING-INCIDENCE; CONVOLUTION PML; DOMAIN ANALYSIS; PILE GROUPS; PROPAGATION;
D O I
10.1016/j.soildyn.2023.108404
中图分类号
P5 [地质学];
学科分类号
0709 ; 081803 ;
摘要
Wave propagation analysis of heterogeneous unbounded media employing domain discretization methods necessitates truncation of the finite computational domain at its truncation boundaries. The Perfectly Matched Layer (PML) is a robust absorbing boundary condition, which is extensively employed to favorably absorb outgoing scattered waves from the interior computational domain, regardless of frequencies and directions of outgoing waves. Nonetheless, most available PMLs are developed for elastic domains, while their developments in domain discretization methods, which consider damping using the Rayleigh damping formulation, are insufficient. In this study, a viscoelastic PML formulation is developed that accommodates the Rayleigh-type damping. The PMLs are basically developed from a coordinate transformation concept using a complex stretching function. Suitable stretching functions can efficiently control the absorptive properties and the spectral character of the PML. Besides, the stretching functions play a crucial role in the long-time stability of transient wave simulations. The standard stretching function usually results in spurious reflections, especially when propagating waves incident to the truncation boundary at near-grazing incidence. On the other hand, the complex-frequency-shifted (CFS) stretching function has been shown to mitigate long-time instability. Therefore, the proposed viscoelastic PML in the time-domain and frequency-domain are derived using the CFS stretching function for the first time. Comprehensive formulations of the proposed viscoelastic CFS-PML are given, where the displacement field is represented as the sole unknown variable of the problem. Hence, the proposed CFS-PML can be incorporated into any domain discretization method. However, the proposed CFS-PML is further developed in the context of the voxel element method, where element-by-element iterative strategies are implemented to solve the governing equations without storing global matrices in the computer memory. Comprehensive numerical experiments are conducted to demonstrate the validity and efficiency of the proposed viscoelastic CFS-PML formulation in analyzing soil profiles with various damping ratios and properties. Moreover, long-time stability for problems involving waveguides and near-grazing wave incidence is studied, where no instabilities are detected.
引用
收藏
页数:26
相关论文
共 86 条
[1]  
Aslmand M., 2018, Soil Dyn Earthq Eng, V115, P531, DOI DOI 10.1016/J.SOILDYN.2018.07.044
[2]   Decoupled scaled boundary finite element method for analysing dam-reservoir dynamic interaction [J].
Babaee, Reza ;
Khaji, Naser .
INTERNATIONAL JOURNAL OF COMPUTER MATHEMATICS, 2020, 97 (08) :1725-1743
[3]   Perfectly matched layers for transient elastodynamics of unbounded domains [J].
Basu, U ;
Chopra, AK .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 2004, 59 (08) :1039-1074
[4]   Perfectly matched layers for time-harmonic elastodynamics of unbounded domains: theory and finite-element implementation [J].
Basu, U ;
Chopra, AK .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2003, 192 (11-12) :1337-1375
[5]   Free and forced vibrations of a segmented bar by a meshless local Petrov-Galerkin (MLPG) formulation [J].
Batra, R. C. ;
Porfiri, M. ;
Spinello, D. .
COMPUTATIONAL MECHANICS, 2008, 41 (04) :473-491
[6]  
Bécache E, 2004, COMPUTAT GEOSCI, V8, P255
[7]   On the long-time behavior of unsplit perfectly matched layers [J].
Bécache, E ;
Petropoulos, PG ;
Gedney, SD .
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, 2004, 52 (05) :1335-1342
[8]   Dynamic elastic analysis of 3-D flexible pavements under moving vehicles: A unified FEM treatment [J].
Beskou, Niki D. ;
Tsinopoulos, Stephanos V. ;
Theodorakopoulos, Dimitrios D. .
SOIL DYNAMICS AND EARTHQUAKE ENGINEERING, 2016, 82 :63-72
[9]   A high-order doubly asymptotic continued-fraction solution for frequency-domain analysis of vector wave propagation in unbounded domains [J].
Chen, Denghong ;
Wang, Min ;
Dai, Shangqiu ;
Peng, Gang .
SOIL DYNAMICS AND EARTHQUAKE ENGINEERING, 2018, 113 :230-240
[10]   Numerical modelling of wave propagation in anisotropic soil using a displacement unit-impulse-response-based formulation of the scaled boundary finite element method [J].
Chen, Xiaojun ;
Birk, Carolin ;
Song, Chongmin .
SOIL DYNAMICS AND EARTHQUAKE ENGINEERING, 2014, 65 :243-255