Improvement research on the stability of explicit integration algorithms with 3D viscoelastic artificial boundary elements

Abstract

PurposeWhen explicit integral analysis is performed on a numerical model with viscoelastic artificial boundary elements, an instability phenomenon is likely to occur in the boundary area, reducing the computational efficiency of the numerical calculation and limiting the use of viscoelastic artificial boundary elements in the explicit dynamic analysis of large-scale engineering sites. The main purpose of this study is to improve the stability condition of viscoelastic artificial boundary elements.Design/methodology/approachA stability analysis method based on local subsystems was adopted to analyze and improve the stability conditions of three-dimensional (3D) viscoelastic artificial boundary elements. Typical boundary subsystems that can represent the localized characteristics of the overall model were established, and their analytical stability conditions were derived with an analysis based on the spectral radius of the transfer matrix. Then, after analyzing the influence of each physical parameter on the analytical-stability conditions, a method for improving the stability condition of the explicit algorithm by increasing the mass density of the artificial boundary elements was proposed.FindingsNumerical wave propagation simulations in uniform and layered half-space models show that, on the premise of ensuring the accuracy of the viscoelastic artificial boundary, the proposed method can effectively improve the numerical stability and the efficiency of the explicit dynamic calculations for the overall system.Originality/valueThe stability improvement method proposed in this study are significant for improving the applicability of viscoelastic artificial boundary elements in explicit dynamic calculations and the calculation efficiency of wave analysis at large-scale engineering sites.

Study on the Influence of a Rubber-Modified Soil Isolation Layer on the Isolation Performance of Frame Structures with Different Foundation Forms

“…The shock absorbing layer of 50 mm thickness is better than the shock absorbing layer of 150 mm and 100mm thickness. Li et al [17] proposed a numerical stability analysis method based on the subsystem of a three-dimensional viscoelastic artificial boundary element explicit algorithm. This method identified the subsystem that governs the stability of the overall value system of artificial boundaries.…”