内域波动数值模拟的显式方法

刘恒

(中国地震局工程力学研究所，哈尔滨 150080)

内域波动数值模拟的显式方法

刘恒

(中国地震局工程力学研究所，哈尔滨 150080)

如何更精确高效地模拟大型、复杂系统内域的波动是发展和完善近场波动数值模拟技术的一个重要研究课题。内域波动的数值模拟通常采用计算量较小的显式方法，但现有的时空解耦显式有限元方法的精度只有二阶;低精度不仅影响数值模拟的精度，而且制约着计算效率的提高。鉴于此，本文旨在探索内域波动数值模拟具有更高精度且稳定的显式方法。作者发展了内域波动数值模拟的现有显式有限元解耦技术，提出了一种高精度且稳定的显式数值模拟方法。

(1)依据波速有限的概念，从波动方程的精确解出发提出了一种新的显式数值模拟方法。此方法与现有有限元技术相似之处在于:适于处理非规则网格节点，且节点递推公式是具有显式和时空解耦特征。其不同之处在于:新方法可给出时空离散精度皆为2M阶的稳定格式，M为正整数。本文通过一维模型详细论述了这一方法的可行性:建立了非规则网格节点递推公式，详细分析了均匀网格标量波动数值模拟的精度和稳定性，提出了构建时空精度皆为2M阶(M为正整数)的稳定递推公式的技术途径，并以构建二阶(M=1)和四阶(M=2)公式为例予以说明。

(2)将一维情形研究结果推广到高维情形，分别建立了二维、三维非规则网格节点的递推公式。针对二维正方形均匀网格详细论述了时空离散精度皆为2M阶的稳定递推公式的构建方法，给出了二阶和四阶递推公式的具体结果;并对三维立方体均匀网格模型作了简要的讨论。

(3)本文通过一维、二维模型算例检验了波动方程数值模拟显式方法的精度和稳定性等理论结果，指出了高阶公式对提高计算效率的价值。

(4)从空间解耦有限元常微分方程组(结构动力学方程)出发，通过被积函数的拉格朗日多项式内插和分部积分导出了一组具有高阶精度的时域显式积分格式，并将其推广为适用于一般数学物理方程的显式积分格式。作者以一个简单的线性时不变系统为例，初步考察了此积分格式的稳定性。

内域波动;显式方法;高精度;稳定性;有限元

TU435;

A;

10.3969/j.issn.0235-4975.2010.01.014