[2018年最新整理]4-有限差分法-20161119.ppt

（5）按照公式（2－1）计算应力的分量。 说明： 1．以上是针对单连体导出的结果。对于多连体，情况就不象这样简单。 2. 如果一部分边界是曲线的，或是不与坐标轴正交，则边界附近将出现不规则的内结点。对于这样的结点，差分方程（2－2）必须加以修正。 4.6 例 深梁的应力函数差分解 现以如图所示的混凝土深梁为例，应用应力函数的差分解求出应力分量。已知混凝土深梁上边受有均布向下的铅直荷载q，并由下角点处的反力维持平衡。 (1）计算边界上各结点的φ 、 和 值。 取A为基点，且 由上面公式所得的计算结果见下表。 （2）计算边界外以行各虚结点处的值。 由式（2-6）及前表可得 (3）边界内各结点的差分方程,由式（2-2）可知 联立求解上式，可得（以qh2为单位 ) （4）计算结点外一行各结点处的值。由(a)、(b)、(c)可得 （5）计算应力。对于结点M，由式（1）可知 同理可得 沿着梁的中线MA，的变化如下图和右图所示。 差分方法的优缺点 优点： 1.差分法是解微分方程边值问题和弹性力学问题的有效方法； 2.差分法简便易行； 3.对于某些结构，为了更精确地分析局部的应力状态，可以用差分法进行分析； 缺点： 1.对于曲线边界和斜边界等产生的不等间距网格的处理，比较麻烦和易于出错； 2.比较比较适用于求解二维问题或平面问题； 3.比较适用于等间距网格，对于应力变化较为剧烈时，需要采用二次网格进行计算。 课堂作业 用差分法计算下图中A和B点的应力分量。 F F B A a a 2 3 4 1 6 5 X Y ·7 * 考察 Richardson 格式的稳定性。 用 表示计算 所产生的误差，如果右端 无误差存在， 则 满足： 取 （14） 假设k - 1层之前无误差存在。即 ，而在第k层产生了 误差。 ，这一层其它点也无误差，而且在计算过程 中不再产生新的误差，利用（14）式算出误差? 的传播如下表： r = ? 时 Richardson 格式的误差传播 j j0 – 4 j0 – 3 j0 – 2 j0 –1 j0 j0 +1 j0 +2 j0 +3 j0 +4 k ?? ? -2? ? ? -4? 7? 4? ? ? -6? 17? -24? 17? -6? ? ? -8? 31? -68? 89? -68? 31? -8? ? -10? 49? -144? 277? -388? 277? -144? 49? -10? 71? -260? 641? -109? 1311? -109? 641? -260? 71? r≤ 1/2 时古典显格式的误差传播 j j0 – 4 j0 – 3 j0 – 2 j0 –1 j0 j0 +1 j0 +2 j0 +3 j0 +4 k ?? 0.5? 0 0.5? 0.25? 0 0.5? 0 0.25 ? 0.125? 0 0. 375? 0 0.375 ? 0 0.125 ? 0.0625? 0 0.25? 0 0.375? 0 0.25? 0 0.0625? 如果选用 r = ? 时的古典显格式，误差方程为： 差分格式关于初值稳定的实际含义是：如果其解 在某一层存在误差，则由它引起的以后各层上的误差不 超过原始误差的M倍（M为与? 无关的常数）。 因此，在稳定的条件下，只要初始误差足够小， 以后各层的误差也能足够小。 以上构造的几种差分格式中， 古典显格式：r ≤1/2时稳定 古典隐格式：绝对稳定 Richardson格式：绝对不稳定 六点对称格式：绝对稳定。 稳定性概念： 4.4 双曲型方程的差分解法 一阶线性双曲型方程最简单的形式为 （4.4.1) 当给定初始条件 (4.4.2) 以后，容易验证，双曲型方程(8.4.1)的解为： (4.4.3) 也就是说，在平面 x t上，沿着 （k 是常数） 这样的直线，u 的值保持不变。这种直线叫做特征