Multi-Directional Viscous Damping Absorbing Boundary in Numerical Simulation of Elastic Wave Dynamic Response

被引：0

作者：

Zhao, Jianguo ^{[1
]}

Yu, Yang ^{[2
]}

Xu, Hao ^{[1
,3
]}

Zhang, Rongtang ^{[1
]}

Ma, Yuxi ^{[1
]}

Li, Jialiang ^{[1
]}

机构：

[1] Wuhan Polytech Univ, Sch Civil Engn & Architecture, Wuhan 430023, Peoples R China

[2] Univ New South Wales, Ctr Infrastructure Engn & Safety, Sch Civil & Environm Engn, Sydney, NSW 2052, Australia

[3] Engn Res Ctr Key Technol Phosphogypsum Rd Construc, Xiaogan 420900, Hubei, Peoples R China

来源：

APPLIED SCIENCES-BASEL | 2024年 / 14卷 / 05期

关键词：

finite element method; numerical simulation; seismic wave field; absorbing boundary; semi-infinite space; structural non-destructive testing; PERFECTLY MATCHED LAYER; FINITE-ELEMENT-ANALYSIS; TIME-STEP CONSTRAINTS; PML; PROPAGATION; DOMAIN; IMPLEMENTATION; ELASTODYNAMICS;

D O I：

10.3390/app14052183

中图分类号：

O6 [化学];

学科分类号：

0703 ;

摘要：

Numerical seismic wave field simulation is essential for studying the dynamic responses in semi-infinite space, and the absorbing boundary setting is critical for simulation accuracy. This study addresses spherical waves incident from the free boundary by applying dynamic equations and Rayleigh damping. A new multi-directional viscous damping absorbing boundary (MVDB) method is proposed based on regional attenuation. An approximate formula for the damping value is established, which can achieve absorbing the boundary setting by only solving the mass damping coefficients without increasing the absorbing region grid cells or depending on the spatial and temporal walking distance. The validity and stability of the proposed method are proven through numerical calculations with seismic sources incident from different angles. Meanwhile, the key parameters affecting the absorption of the MVDB are analyzed, and the best implementation scheme is provided. In order to meet the requirements of mediums with different elastic parameters for boundary absorption and ensure the high efficiency of numerical calculations, the damping amplitude control coefficients k can be set between 1.02 and 1.12, the thickness of the absorbing region L is set to 2-3 times of the wavelength of the incident transverse wave, and the thickness of the single absorbing layer is set to the size of the discrete mesh of the model Delta l.

引用

页数：22