zhenbo

ISSN 2096-7780 CN 10-1665/P

段雪良, 周正华, 卞祝, 韩轶, 赵玲, 李政, 廖成亮, 贺家聪, 刘伟. 插值算法对透射边界数值模拟结果的影响分析[J]. 地震科学进展 . DOI: 10.19987/j.dzkxjz.2024-046
引用本文: 段雪良, 周正华, 卞祝, 韩轶, 赵玲, 李政, 廖成亮, 贺家聪, 刘伟. 插值算法对透射边界数值模拟结果的影响分析[J]. 地震科学进展 . DOI: 10.19987/j.dzkxjz.2024-046
Duan Xueliang, Zhou Zhenghua, Bian Zhu, Han Yi, Zhao Ling, Li Zheng, Liao Chengliang, He Jiacong, Liu Wei. Analysis of the influence of interpolation algorithm on the numerical simulation results of transmission boundary[J]. Progress in Earthquake Sciences. DOI: 10.19987/j.dzkxjz.2024-046
Citation: Duan Xueliang, Zhou Zhenghua, Bian Zhu, Han Yi, Zhao Ling, Li Zheng, Liao Chengliang, He Jiacong, Liu Wei. Analysis of the influence of interpolation algorithm on the numerical simulation results of transmission boundary[J]. Progress in Earthquake Sciences. DOI: 10.19987/j.dzkxjz.2024-046

插值算法对透射边界数值模拟结果的影响分析

Analysis of the influence of interpolation algorithm on the numerical simulation results of transmission boundary

  • 摘要: 一个时间步内外形波以人工波速走过的距离与网格尺寸往往不一致,因此在应用透射边界公式时,一般通过时间插值或空间插值的方法将透射边界的计算结点用有限元结点表示。本文导出了拉格朗日插值算法、埃尔米特插值算法以及三次样条插值算法的透射边界一次透射插值公式,探讨不同插值算法对透射边界数值模拟结果的影响。首先将弹性半空间中入射P/SV波自由场问题的数值解与解析解对比,讨论三种插值算法对计算精度的影响,然后以半圆河谷地形入射P/SV波散射问题为例,对比不同插值算法计算得到的结果。数值结果表明:不同插值算法对透射边界数值模拟结果有显著影响,三种插值算法中,三次样条插值算法计算精度最高,其次是拉格朗插值算法,相比之下,埃尔米特插值算法计算得到的精度较低。相比于自由场,散射问题的计算结果受三种插值算法的影响较大,拉格朗日插值算法与三次样条插值算法计算结果相近,而埃尔米特插值算法的计算结果差异较大。本研究成果有助于为在应用透射边界模拟场地动力反应过程时,选择合适的插值算法提供科学依据。

     

    Abstract: The distance traveled by contour wave with artificial wave velocity in a step-time is often inconsistent with the mesh size. Therefore, when the transmission boundary formula is applied, the calculation nodes of transmission boundary are usually represented by finite element nodes through the method of time interpolation or spatial interpolation. In this paper, the transmission boundary interpolation formulas of Lagrange interpolation algorithm, Hermite interpolation algorithm and cubic spline interpolation algorithm are derived, and the influence of different interpolation algorithms on the numerical simulation results of transmission boundary is discussed. Firstly, the numerical and analytical solutions of the incident P/SV wave free field problem in elastic half space are compared, and the influence of three interpolation algorithms on the calculation accuracy is discussed. Then, taking the scattering problem of incident P/SV wave in semi-circular valley terrain as an example, the results obtained by different interpolation algorithms are compared. The numerical results show that: Different interpolation algorithms have significant influence on the numerical simulation results of transmission boundary. Among the three interpolation algorithms, the cubic spline interpolation algorithm has the highest accuracy, followed by Lagrange interpolation algorithm, and the Hermite interpolation algorithm has lower accuracy. Compared with the free field, the results of scattering problem are more affected by three interpolation algorithms, and the results of Lagrange interpolation algorithm and cubic spline interpolation algorithm are similar, while the results of Hermite interpolation algorithm are different. The findings of this study can serve as a scientific foundation for the selection of an appropriate interpolation algorithm when utilizing transmission boundaries to simulate the dynamic reaction process of the site.

     

/

返回文章
返回