王圣柱, 周建科, 张奎华, 王艺豪. 2015: 有限元三角网格中波动的频散分析及数值仿真. 地震学报, 37(3): 493-507. DOI: 10.11939/jass.2015.03.012
引用本文: 王圣柱, 周建科, 张奎华, 王艺豪. 2015: 有限元三角网格中波动的频散分析及数值仿真. 地震学报, 37(3): 493-507. DOI: 10.11939/jass.2015.03.012
Wang Shengzhu, Zhou Jianke, Zhang Kuihua, Wang Yihao. 2015: Dispersion analysis and numerical simulation of wave motion in finite element algorithm with triangular meshes. Acta Seismologica Sinica, 37(3): 493-507. DOI: 10.11939/jass.2015.03.012
Citation: Wang Shengzhu, Zhou Jianke, Zhang Kuihua, Wang Yihao. 2015: Dispersion analysis and numerical simulation of wave motion in finite element algorithm with triangular meshes. Acta Seismologica Sinica, 37(3): 493-507. DOI: 10.11939/jass.2015.03.012

有限元三角网格中波动的频散分析及数值仿真

Dispersion analysis and numerical simulation of wave motion in finite element algorithm with triangular meshes

  • 摘要: 波动问题有限元离散后会引起数值误差, 数值频散的本质就是数值误差传播引起的非物理解. 数值频散不仅没有实际意义, 而且还会影响对真实波动现象的认识. 为厘清有限元三角网格中波动数值频散的影响因素, 本文推导了集中质量矩阵和一致质量矩阵的频散函数, 同时给出了组合质量矩阵的频散函数, 并对不同质量矩阵的数值频散进行了对比研究. 理论分析和数值计算结果表明: 有限元三角网格中波动的数值频散受网格布局、 波传播方向、 单元网格纵横比以及质量矩阵的影响; 一致质量矩阵的数值频散比集中质量矩阵更易受到波传播方向的影响; 不合理的三角网格单元会对数值相速度(数值频散)产生不良影响; 正三角网格中波动的数值频散几乎不受波传播方向的影响; 一致质量矩阵与集中质量矩阵的线性组合能够有效地压制数值频散.

     

    Abstract: The finite element discretization of wave motion can cause numerical error. The numerical dispersion is in essence the non-physical solution caused by the propagation of numerical error. Not only has the numerical dispersion no practical meaning, but also it can affect the understanding of real fluctuations. In order to clarify the influence factors of numerical dispersion in finite element algorithm with triangular meshes, the dispersion functions of lumped mass matrix and consistent mass matrix are derived respectively, and the dispersion function of combined mass matrix is also given. And then the numerical dispersion of different mass matrices are compared. The results of theoretical analysis and numerical simulation indicate: ① The numerical dispersion in finite element algorithm with triangular meshes depend on the mesh layout, direction of wave propagation, the ratio of vertical length to horizontal length and the mass matrix; ② The numerical dispersion of consistent mass marix is more easily affected by wave propagation direction than that of lumped mass matrix; ③ The irrational triangular meshes have a bad effect on numerical phase velocity (numerical dispersion); ④ Equilateral triangular meshes can provide superior results regardless of the propagation direction; ⑤ The linear combination of lumped mass matrix and consistent mass matrix can effectively suppress numerical dispersion.

     

/

返回文章
返回