采用黏弹性人工边界时显式算法稳定性条件

黏弹性人工边界是处理无限域波动问题常用的数值模拟方法。采用显式时域逐步积分算法进行计算时，受黏弹性人工边界的阻尼、刚度等影响，人工边界区的稳定性比内部计算域的更严格，尚无明确、实用的稳定性判别准则，这限制了黏弹性人工边界在显式动力分析中的应用。针对二维黏弹性人工边界，利用基于局部子系统的稳定性分析方法和基于传递矩阵谱半径的稳定性判别准则，给出了可代表整体模型局部特征的不同边界子系统的稳定性条件解析解。通过对比分析不同计算区域的稳定性条件及其影响因素，证明了整体模型的稳定性由角点子系统控制。在此基础上，获得了含黏弹性人工边界的整体模型在显示动力计算中的统一稳定性判别准则和简化实用计算方法。在实际应用中，令积分时间步长满足稳定性条件，即可顺利完成整体模型的动力计算。以上研究可为将黏弹性人工边界应用于显式动力计算时积分时间步长的合理选取提供参考。

Stability conditions of explicit algorithms when using viscoelastic artificial boundaries

Viscoelastic artificial boundary is a commonly used numerical simulation method to deal with the wave propagation problems in an infinite domain. When the explicit time-domain stepwise integration algorithm is adopted for such numerical analysis, the stability conditions of the artificial boundary area are more stringent than those of the internal domain due to the influence of the damping and stiffness of the viscoelastic artificial boundary. However, there is currently no clear and practical stability criterion for this problem, which affects the reasonable selection of the integral time step when using the viscoelastic artificial boundaries, and further restricts the application of viscoelastic artificial boundary in the explicit dynamic analysis. Aiming at the two-dimensional (2D) viscoelastic artificial boundary, two typical types of boundary subsystem that can represent the localized characteristics of the overall numerical model, namely the edge boundary subsystem and the corner boundary subsystem, were established and their motion equations as well as the transfer matrixes were obtained according to the stability analysis method based on the local subsystem. Then through the stability criteria based on the spectral radius of the transfer matrix, the analytical solutions of the stability conditions of different local subsystems were derived. Through the comparative analysis of the stability conditions of different calculation areas and their influencing factors, it is found that the stability of the overall model is controlled by the corner boundary subsystem. On that basis, a uniform stability criterion and a simplified practical calculation method of the stability condition for the overall model with 2D viscoelastic artificial boundary in explicit dynamic calculations were proposed. In practical applications, the dynamic calculation of the overall model can be successfully completed once the integral time step meets the proposed stability condition of the numerical system. This study provides theoretical guidance for the reasonable selection of the integral time step when applying 2D viscoelastic artificial boundaries in explicit dynamic calculations.