8.有限差分法基础.pdfVIP

有限差分法基础 成都大学 张永强 1 有 限 元 思 想 2 工程中常用的数值解法：有限单元法，差分法。 有限单元法: 是以有限个单元的集合体来代替 连续体，属于物理上的近似。 差分法: 是把弹性力学的基本方程和边界条件 （一般均为微分方程）近似地改用差分方程 （代 数方程）来表示，把求解微分方程的问题改换成 为求解代数方程的问题，属于数学上的近似。 3 1.有限差分法简介 有限差分法（FDM）是计算机数值模拟最早采用的方法， 至今仍被广泛运用。 该方法将解域划分为差分网格，用有限个网络节点代替 连续的求解域。 有限差分法通过泰勒级数展开等方法，把控制方程中的 导数用网格节点上的函数值得差商代替进行离散，从而 建立以网格节点上的值为未知数的代数方程组。 该方法是一种直接将微分问题变为代数问题的近似值解 法，数学概念直观，表达简单，是发展较早且比较成熟 的数值方法。 4 2.有限差分法的数学基础 有限差分法的数学基础是用差分代替微分 用差商代替微商的意义：是用函数在某区域内的平均 变化率来代替函数的真是变化率。 而根据泰勒级数展开可以看出，用差商代替微商必然 会带来阶段误差，相应的用差分方程代替微分方程也会 带来误差，因此，在应用有限差分法进行计算的时候， 必须注意差分方程的形式，建立方法及由此产生的误差 如何减小误差？ 5 3.有限差分解题基本步骤  1）建立微分方程,根据问题的性质选择计算区域，建立 微分方程式，写出初始条件和边界条件。 (物理问题) 2）构建差分格式,首先对求解域进行离散化，确定计算 节点，选择网格布局，差分形式和步长；然后以有限差 分代替无限微分，以差商代替微商，以差分方程代替微 分方程及边界条件。 3）求解差分方程,差分方程通常是一组数量较多的线性 代数方程，其求解方法主要包括两种:精确法和近似法。 精确法又称直接发，包括矩阵法，高斯消元法及主元素 消元法等；近似法又称间接法，以迭代法为主，主要包 括直接迭代法，间接迭代法以及超松弛迭代法。 4）精度分析和检验,对所得到的数值进行精度与收敛性 分析和检验 6 4.商用有限差分软件简介 商用有限差分软件主要包括FLAC、UDEC/3DEC和PFC程序， （1） 是一个基于显式有限差分法的连续介质程序，主要用来进行土 （2）UDEC/3DEC是针对岩体不连续问题开发，用于模拟非连续介质 （3）PFC是利用显式差分算法和离散元理论开发的微、细观力学程 以及颗粒的流动(大位移）问题。 7 本节主要内容 1.差分原理 2.逼近误差 3.差分方程 4.截断误差 5.相容性 6.收敛性 7.稳定性 8.Lax等价定理