《有限差分法求解偏微分方程》由会员分享，可在线阅读，更多相关《有限差分法求解偏微分方程（18页珍藏版）》请在金锄头文库上搜索。
1、word有限差分法求解偏微分方程摘要：本文主要使用有限差分法求解计算力学中的系统数学模型，推导了有限差分法的理论根底，并在此根底上给出了局部有限差分法求解偏微分方程的算例验证了推导的正确性与操作可行性。关键词：计算力学，偏微分方程，有限差分法Abstract：This dissertation mainly focuses on solving the mathematic model of putation mechanics with finite-difference method. The theoretical basis of finite-difference is derive
2、d in the second part of the dissertation, and then I use MATLAB to program the algorithms to solve some partial differential equations to confirm the correctness of the derivation and the feasibility of the method.Key words：putation Mechanics, Partial Differential Equations, Finite-Difference Method
3、1 引言机械系统设计常常需要从力学观点进展结构设计以与结构分析，而这些分析的前提就是建立工程问题的数学模型。通过对机械系统应用自然的根本定律和原理得到带有相关边界条件和初始条件的微分积分方程，这些微分积分方程构成了系统的数学模型。求解这些数学模型的方法大致分为解析法和数值法两种，而解析法的局限性众所周知，当系统的边界条件和受载情况复杂一点，往往求不出问题的解析解或近似解。另一方面，计算机技术的开展使得计算更准确、更迅速。因此，对于绝大多数工程问题，研究其数值解法更具有实用价值。对于微分方程而言，主要分为差分法和积分法两种，本论文主要讨论差分法。2 有限差分法理论根底2.1 有限差分法的根本思想
4、当系统的数学模型建立后，我们面对的主要问题就是微分积分方程的求解。根本思想是用离散的只含有限个未知量的差分方程组去近似地代替连续变量的微分方程和定解条件，并把差分方程组的解作为微分方程定解问题的近似解。将原方程与边界条件中的微分用差分来近似，对于方程中的积分用求和或与机械求积公式来近似代替，从而把原微分积分方程和边界条件转化成差分方程组。有限差分法求解偏微分方程的步骤主要有以下几步：n 区域离散，即把所给偏微分方程的求解区域细分成由有限个格点组成的网格，这些离散点称作网格的节点；n 近似替代，即采用有限差分公式替代每一个格点的导数；n 逼近求解，换而言之，这一过程可以看作是用一个插值多项式与其
5、微分来代替偏微分方程的解的过程。从原如此上说，这种方法仍然可以达到任意满意的计算精度。因为方程的连续数值解可以通过减小独立变量离散取值的间格，或者通过离散点上的函数值进展插值计算来近似得到。理论上，当网格步长趋近于零时，差分方程组的解应该收敛于准确解，但由于机器字节的限制，网格步长不可能也没有必要取得无限小，那么差分法的收敛性或者说算法的稳定性就显得至关重要。因此，在运用有限差分法时，除了要保证精度外，还必须要保证其收敛性。2.2 系统微分方程的一般形式1由于大多数工程问题都是二维问题，所以得到的微分方程一般都是偏微分方程，对于一维问题得到的是常微分方程，解法与偏微分方程类似，故为了不是一般性
6、，这里只讨论偏微分方程。由于工程中高阶偏微分较少出现，所以本文仅仅给出二阶偏微分方程的一般形式，对于高阶的偏微分，可进展类似地推广。二阶偏微分方程的一般形式如下：其中，为弹性体上的某一特征物理量连续函数。当A、B、C都是常数时，1式称为准线性，有三种准线性方程形式：n 如果，如此称为椭圆型方程；n 如果，如此称为抛物型方程；n 如果，如此称为双曲型方程。椭圆型方程主要用来处理稳态或静态问题，如热传导等问题；抛物线方程主要用来处理瞬态问题，如渗透、扩散等问题；双曲型方程主要用来处理振动问题，如玄震动、薄膜震动等问题。除了上述微分方程外，必须给出定解条件，通常有如下三类：n 第一类边界条件Diri
7、chlet条件：；n 第二类边界条件Neumann条件：；n 第三类边界条件Robin条件：；其中，为求解域的边界，为的单位外法矢， 。第二类和第三类边界条件统称为导数边界条件。2.3 有限差分方程的数学根底2.3.1 一元函数导数的差分公式一个函数在x点上的导数，可以近似地用它所临近的两点上的函数值的差分来表示。函数在处的泰勒展式如下：23对一个单变量函数，以步长将区间等距划分，我们得到一系列节点：，然后求出在这些节点上的近似值。与节点相邻的节点有和，因此在点处可以构造如下形式的展开式：4由式3和式4可得到：5n 一阶向前差分：6n 一阶向后差分：7n 一阶中心差分：不妨，记，如此式5、6、
8、7分别简写为：8n 一阶向前差分：9n 一阶向后差分：10n 一阶中心差分：根据式8、式9和式10，可得二阶差分：11n 二阶向前差分：12n 二阶向后差分：13n 二阶中心差分：差分公式13是以相隔2h的两结点处的函数值来表示中间结点处的一阶导数值，可称为中点导数公式。式11和式12是以相邻三结点处的函数值来表示一个端点处的一阶导数值，可称为端点导数公式。应当指出：中点导数公式与端点导数公式相比，精度较高。因为前者反映了结点两边的函数变化，而后者却只反映了结点一边的函数变化。因此，我们总是尽可能应用前者，而只有在无法应用前者时才不得不应用后者。14但是，由于式11中的各阶导数均使用的是向前差
9、分，导致用到的节点不相邻，同时为了均衡误差，将节点处用到的一阶差分换成向后差分，如此式11修正为：同理，根据上述推导过程，可得到任意阶的差分公式：15n n阶向前差分：16n n阶向后差分：17n n阶中心差分：说明，上述公式中各节点处前一阶导数的代入可能存在不一致，可能是向前差分、向后差分或者中心差分，从而使最终的公式在系数上存在差异。当然，也可以对各相邻节点进展需要阶数的泰勒展开，从而建立方程组直接求各阶导数。2.3.2 微分方程转化为线性方程18由于三种类型的微分方程解法类似，故这里仅以椭圆型微分方程为例，将微分方程转化为代数方程，对于双曲型和抛物型方程依次类推即可。不妨记：称为拉普拉斯
10、算子，和是求解域上的连续函数。假设求解区域为：，将求解区域划分成个网格，其中：，如图1所示。记，如此根据式14可得到：图1 五点差分公式式18也称为五点差分公式，同理根据式12和式13可分别得到向前差分公式19和向后差分公式20，如图2所示。19n 向前差分20n 向后差分图2 向前差分左和向后差分右图3 中心差分、向前差分和向后差分的拉普拉斯算子表示21利用中心差分公式18，由于式18在点处具有二阶精度，所以式18可近似改写成下式：根据椭圆方程的具体形式可以将其分为以下三种形式：n 拉普拉斯Laplace方程：n 泊松Poison方程：n 赫耳墨次Helmholtz方程：根据式21，可建立三
11、种不同形式椭圆方程的代数方程如下：n 拉普拉斯方程：化简后得到拉普拉斯方程的计算公式：2223n 泊松方程：n 赫耳墨次方程：242.3.3 建立有限差分方程组根据式2224建立方程组，但是需要知道对应的边界条件才能使方程组存在定解，根据2.2中可知，边界条件一般分为狄利克雷边界条件和导数边界条件两种，下面分别给出这两种边界条件的有限差分方程组的建立过程：n 狄利克雷边界条件：对于狄氏条件而言，给出了边界上各节点出的函数计算公式，直接代入节点值计算即可，如下所示为矩形区域的边界点计算：左边界25右边界下边界上边界n 导数边界条件：26以右边界点为例，对于右边界点，根据Neumann条件可得下式
12、：对于拉普拉斯方程，根据计算公式22，对于边界上的点可得：2728显然，上式中的在求解域外，是未知量。根据中心差分公式10可得到：根据式28可得到逼近表示：，并且具有2阶逼近精度，代入式27可得下式：29同理，对于其它边界可获得如下边界方程：30下边界上边界左边界图4 Neumann条件算子对于泊松方程和赫耳墨次方程同样根据上述方法，获得边界条件的线性方程，然后将这些方程添加到式2224所建立的方程组中，从而建立起个元的线性方程组，解该方程组即可获得各节点的函数值。对于上述过程建立的线性方程组的求解，可采用多种方法，比如Jacob迭代法、Gauss-Seidel迭代法、超松弛迭代法SOR法、高
13、斯消元法等方法求解。2.4 有限差分法的收敛性和稳定性由于迭代法必须保证收敛性，所以在解有限差分方程组时还应保证其收敛性，也就是通常所说的算法稳定性。有限差分法的算法稳定性可以通过特征值方法、傅里叶变换冯诺依曼条件以与能量估计等方法来判断，下面给出常用的冯诺依曼条件：n 向前差分：，绝对收敛；n 向后差分：，绝对收敛；n 中心差分：对任何的对不收敛；假设求解域内方向网格划分的步长为，方向网格划分的步长为，将偏微分方程化为标准形式，具体来说标准形式如下：31n 双曲方程：对于式31所示的双曲方程，冯诺依曼条件为：.32n 抛物方程：对于式32所示的抛物方程，冯诺依曼条件为：.33n 椭圆方程：对
14、于式33所示的椭圆方程，冯诺依曼条件为：.34为了使算法在任何情况下都能保持稳定性，去掉对网格划分的冯诺依曼条件，通常采用隐式方案，对五点差分公式中的节点所在的行做差分，然后把这些差分的加权作为中心点的差分值，如此拉普拉斯算子可修正为：利用式34进展计算时，稳定性没有任何限制。取不同的值得到不同的差分公式，通常取.35为了提高计算精度，很明显的一个措施就是网格细分，但是由于随着网格步长的减小，未知量的数目将会呈指数增长，网格划分太细会导致计算量过于庞大而无法计算。通常，我们可以通过提高逼近的精度，采用更高精度的差分公式，例如对公式21进展修改，可得到九点差分公式：图5 九点差分公式3 有限差分法求解实例根据上述推导有限差分法理论，对于不同类型的偏微分方程建立有限差分方程组，采用mat lab编程给出一些计算实例如下：1. 椭圆型方程n 拉普拉斯方程：；求解域：下面分别给出拉普拉斯方程在不同的边界条件下的解。a) 狄利克雷边界条件：下边界：上边界：左边界：右边界：图6 狄利克雷边界条件下拉普拉斯方程
