4-有限差分法20161116.docxVIP

第4章 有限差分法第4章有限差分法刚刚介绍的变分法和加权余量法，都是将实际问题的控制方程转化为积分形式进行求解，而有限差分方法不同，它是用网格覆盖求解区域和时间域，用差分近似替代控制方程中的微分，把求解微分方程的问题改换成求解代数方程问题，进行近似求解。有限差分法是一种经典数值方法，当网格细化或节点较多时，近似解的精度可以得到改进，对于几何形状规则且边界条件简单的问题，差分法有很高的求解精度和收敛性，但用于求解几何形状和边界条件复杂的问题时，它的精度将会降低，甚至产生困难。有限差分法在流体力学和爆炸力学中得到广泛应用，也曾用以求解弹性力学边值问题。目前，弹性静力学一般都不用差分法进行求解了，但是对于弹性动力学等初值问题，仍然要借助于有限差分法予以解决。有限差分法的关键是构造差分格式，保证解的收敛性和解的稳定性。这里简要介绍椭圆型方程、抛物型方程和双曲型方程的差分格式构造方法。4.1有限差分法的原理差分法是微分方程的一种近似数值解法。它不是去寻求函数式的解答，而是去求出函数在一些网格结点上的数值。具体地讲，差分法就是先将求解域离散化，在这些离散点上用差商近似地代替导数，从而把基本方程和边界条件近似地改用差分方程来表示，把求解微分方程的问题改换为求解代数方程的问题。考虑一维热传导方程的第一类初边值问题（4-1）其中，，a为正常数。假定和在相应区域光滑，并且在满足相容条件，使上述问题有