岩土数值分析2-1.docVIP

有限差分法 在岩土工程中常常会遇到大量的偏微分方程的求解问题，许多具体过程都可以用偏微分方程的定解问题来描述，例如渗流问题??、土的固结问题等。除了一些特殊问题外，要直接得到偏微分方程的精确解很困难，不但受求解域几何形状、初始条件和边界条件不规则性的限制，而且由于方程中的参数往往依赖于方程的解，使得求解更加困难。随着科学技术的发展，特别市电子计算机的诞生及其运算速度的飞速提高，偏微分方程的近似解法（数值解）在理论上和方法上都有了很大的进展。过去许多难一解决的实际问题现在也得以解决。在偏微分方程的数值解法中应用最广的是有限元法和有限差分法。有限差分法是最早被应用于工程科学中的数学方法之一，作为一种解决偏微分方程定解近似方法，在岩土工程中有着广泛的应用。 2.1有限差分法的基本概念 差分法的基本思想就是要把要求解问题的微分方程及其边界条件用离散的、只含有限个未知数的差分方程（代数方程组）来表示，把求解微分方程的问题转化为求解代数方程的问题，并用代数方程的解作为微分方程的近似解。具体的做法是用差分网格离散求解域，用差商近似代替导数的表达式，得到差分方程组并求解得到差分解，原来微分方程的解可用此差分解来代替。网格划分越细，差分解就越逼近精确解。 2.1.1 有限差分网格的剖分 有限差分法邱界偏微分方程组是先要把连续问题离散化，即把连续的求解区域作网格划分。下面以二维问题为例来说明网格划分。假设所研究的问题是关于空间变量x和时间变量t的偏微分方程组，而研究的区域是xa,b],t[0,T]，如图2-1所示。在x-t平面上画两族平行于坐标轴的直线，在上述区域划分为矩形网格，这些直线的交点称为网格点或节点。一般来说，等距的网格划分较为常见（当然，不等距的网格亦可划分）。设空间方向的距离为，记为h,称其为空间步长；时间方向的步长为，记为。为了研究方便起见，网格划分中的每一个节点（x,t）简记为（i,j）. 经过网格剖分，把连续的区域离散为以下区域（离散点的集合） D= (2-1) 图2-1 差分法网格划分 2.1.2 差分方式 用差商来近似代替导数可得到差分方式。设f(x,，t)为所要求解的某一连续函数。函数在平行于x轴的一根网线上（如直线(i-1,j)-(i,j)-(i+1,j)上），只随x坐标的变化而变化，可将函数展开为以下泰勒级数形式： (2-2) 节点（i-1,j）处有 (2-3) 可得 （向前差分） 也可利用节点（i+1, j）得到上述偏导数的另一种差商形式 （向后差分） 以上两种差商的计算中略去了步长h的二次幂及其以后各项，在连续的某一网格区间内，把函数f简化为按直线变化。向前差分和向后差分被称为偏心差分，是一种最基本、最简单的构造差商的方法，常用来对非对称性变量（例如时间变量等）进行差分计算。用差商来近似代替导数的前提条件是步长h充分小，即网格划分的越细，差分法计算结果越接近精确解答。但是，受计算机存储量和计算速度的限制，网格不可能划分得太细，也就是说单靠细化网格来提高速度不现实。在同样的网格划分下为了得到精度较高的解答，可以通过提高差分公式精度的方法来实现.即在泰勒级数展开式中多取几项,例如分别在节点(i-1,j)和节点(i+1,j)处多去一项得到 (2-6) (2-7) 联立求解可得 (2-8) (2-9) 这种差分公式叫中心差分,其特点是在连续的两段网格区间内,把函数f简化为x的二次函数,看作按抛物线变化.如果要得到更高精度的差分方程,只要在函数的级数展开式中多取几项即可.由于高精度的差分格式中涉及的节点数目太多,应用不便,因而较少采用高精度的差分公式来解决实际问题. 表2-1 常用差分公式表 点导数 向前差分 向后差分 中心差分 同理可得沿t方向的差分公式及x和t的混合导数的差分公式(以中心差分为例) (2-10) (2-