J4 ›› 2012, Vol. 25 ›› Issue (3): 80-87.doi: 10.3976/j.issn.1002-4026.2012.03.017
马燕,王兆清,唐炳涛
MA Yan, WANG Zhao-Qing, TANG Bing-Chao
摘要:
提出一种数值求解波动问题的高精度重心有理插值配点法。对于给定的时间和空间上的计算节点,采用重心有理插值近似未知函数,建立未知函数关于时间和空间变量导数的微分矩阵。将未知函数的重心有理插值近似函数代入波动问题的控制方程,得到波动问题方程和定解条件的离散代数方程组。利用微分矩阵的记号,将离散后的代数方程组写成简洁的矩阵形式。通过置换法施加边界条件和初始条件,求解代数方程组,得到波动问题在计算节点处的位移值。数值算例表明,重心有理插值配点法具有计算公式简单、计算节点适应性好、程序实施方便和计算精度高的优点。
中图分类号:
TB115.1
开放获取 本文遵循知识共享-署名-非商业性4.0国际许可协议(CC BY-NC 4.0),允许第三方对本刊发表的论文自由共享(即在任何媒介以任何形式复制、发行原文)、演绎(即修改、转换或以原文为基础进行创作),必须给出适当的署名,提供指向本文许可协议的链接,同时表明是否对原文作了修改,不得将本文用于商业目的。CC BY-NC 4.0许可协议详情请访问 https://creativecommons.org/licenses/by-nc/4.0