弹性力学及有限元(5-1).pptVIP

有限差分法概述 * 第五章 用差分法和变分法解平面问题 自从弹性力学基本方程建立以后，这些方程在各种问题的边界条件下如何求解，曾经是很多数学家和力学家研究的内容。 但是，对于工程上的许多重要问题，由于边界条件较为复杂，并不能给出函数式解答（即解析解）。因此，弹性力学问题的各种数值解法（近似解法）具有重要的实际意义。差分法和变分法是沿用较久的两种数值解法。 平面问题的差分法 有限差分法概述 差分公式的推导 应力函数的差分解 应力函数差分解的实例 离散数学中相当于微分学的一个分支。它是一种利用差分算子（运算）之间各种形式的关系进行内插的方法。其中差分算子利用函数在等距（节）点上的值来定义。差分法是一种常用的数值方法。 差分法将微分用有限差分代替，把导数用有限差商代替，从而把基本方程和边界条件（微分方程）近似地改用差分方程（代数方程）来表示，把求解微分方程的问题改换为求解代数方程的问题。 差分法在工具书中的解释 差分法 差分法全称有限差分法（Finite Difference Method）， 它是微分方程的一种近似数值解法。 差分法的优点 有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格结点代替连续的求解域。 有限差分法以Taylor级 数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。 差分法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等。 其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 有限差分法的格式 从格式的精度来划分，有：一阶格式、二阶格式和高阶格式。 从差分的空间形式来划分，分为：中心格式和逆风格式。 考虑时间因子的影响，差分格式还可分为：显格式、隐格式、显隐交替格式等。 目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。 构造差分的方法 差分法的应用 求解海洋地质问题——马在田院士被誉为“给地球做CT的人” 应用于气候研究的海冰动力学模式 在地震资料处理中的应用 传热问题 流体力学 结构分析 土的固结方程 弹性地基梁 土木工程 地质勘探 解电磁场边值问题 物理学 稳定温度场和不稳定温度场 气象学 地震工程 当⊿x为微小值时，上式中的约等号是成立的。 §5-1 差分公式的推导 设f(x,y)为弹性体内的某一个连续函数，它可能是某个应力分量或位移分量，也可能是应力函数，等等。 在右图所示的域上，在结点0的近处，函数f(x,y)可以展为Taylor级数如下： 设y=f(x)为自变量x的连续函数，即： 在结点3及结点1，x分别等于(x0-h)及(x0+h)。代入式（a）,得 假定网格间距h是充分小的，因而可以不计它的三次幂及更高次幂的各项，则式（b）及式（c）简化为 联合求解 及 ，得差分公式： 同样可以得到 公式（5-1）至（5-4）是基本差分公式，可导出其他差分公式。 利用公式（5-1）及（5-3），可以导出混合二阶导数的差分公式： 同样可以导出 以上在导出基本差分公式（5-1）至（5-4）过程中，在式（a ）中略去了( x- x0)的三次幂及更高次幂的各项。这样，就把函数f简化为x的二次函数，即在连续两段网格间距之内，把f看作按抛物线变化。 基本差分公式（5-1）至（5-4）常称为抛物线差分公式。 §5-2 应力函数的差分解 由§2-9可知，在不计体力的情况下，平面问题中的应力分量可用应力函数的二阶导数表示如下： 如果在弹性体上织成如图5-1所示的网格，应用差分公式(5-4)，（5-2）,（5-5），则任一结点0处的应力分量可表示成为 应力分量的差分解 为了求得弹性体边界以内各结点处的?值，须利用应力函数的重调和方程,但须首先把它变换为差分方程。为此，要把差分公式（5-6）至（5-8）代入 这样就得出 重调和方程对应的差分方程 即 式(5-10)即为应力函数的差分方程。 边界条件的差分解 为了求得弹性体边界上各结点处的?值，需要应用边界条件： 利用式(a)，上式变换为 由右图可见： 因此，式(b)可改写为 或表达为 将上列二式对s积分，从A点至B点，得 或 注意 或将式（c）代入而得 对s积分，从A点到B点，由分部积分得 即 再将式（d）代入，得 或 式（d）及式（ｅ）表明：若已知固定基点Ａ的 即可根据面力分量求得边界ｓ上任一