梁振动方程的一个稳定的有限差分近似.pdfVIP

维普资讯 第 24巷 第 1期 林 化 工 学 院 学 报 V0I．24 No．1 2007{F2月 JOURNALOFJlLI51INSn．rII．I．EOFCHEMICALTECHNOLOGY Feb． 20 文章编号 ：1007-2853(2007)01-0079-03 梁振动方程的一个稳定的有限差分近似 许士菊’，王长华 (1．吉林大学数学学院，吉林 长春 130012；2．吉化九中 数学教研室，吉林 吉林 132022) 摘要：针对以梁振动问题为背景的一个四阶线性偏微分方程的初边值问题，给出了一个有限差分格式． 其想法是使用适当的非对称格式，实现了问题的稳定的显式计算． 关 键 词：有限差分法；四阶方程；非对称法 中图分类号：O175．2 文献标识码：A 对于大型科学与工程计算问题，算法的稳定可 精度的数值解，差分步长必须很小，此时，提高计 靠与简便是提高计算效率的关键．构造高效率的数 算效率是个重要的问题．通常的显式差分格式虽 值方法一直是人们关心的问题．众所周知，用有限 然具有计算简便的优点，但由于有较严格网比限 差分方法近似求解偏微分方程问题具有重要的理 制(稳定性要求)，计算效率并不理想；而隐式格 论意义和应用价值 ¨。J．但是，对于高阶方程，有限 式又大大增加了计算的复杂性．本文借助于人们 差分方法具有相当的困难，其表现为：显格式的稳 处理抛物方程的非对称差分格式的思想，给出了 定性条件比较苛刻，很难放大时间方向的差分步 一 个求解问题(1)～(3)的一个稳定的可显式计 长；而其隐格式要大规模求解具有较大带宽的线性 算的差分格式，以提高计算的整体效率 代数方程组，其系数矩阵的条件数随着阶数的增加 急剧变坏．因此有必要寻找有效的方法解决这类问 2 问题的差分逼近 题．本文对于以梁振动为背景的一个四阶方程构造 了一个有限差分算法，基本思想是使用非对称格 为了得到稳定的计算方法，首先对问题 (1) 式，实现了近似解的稳定的显式计算． ～ (3)进行转换．引入函数 Ou tl，= k—Ou tt，， == —一— tl， = — 1 问 题 at a 则问题 (1)～(3)转化为如下初边值问题： kOw 0 l 在某些梁振动问题中，需要考虑如下四阶方程： = 一 ， t0 (4) u O +k2u+- — 了