有限差分法的介绍及简单应用_初学者.ppt

第十章 有限差分法 弹性力学的经典解法存在一定的局限性，当弹性体的边界条件和受载情况复杂一点，往往无法求得偏微分方程的边值问题的解析解。因此，各种数值解法便具有重要的实际意义。;第一节 差分方程; 工程中常用得数值解法有有限单元法和差分法。; 我们在弹性体上，用相隔等间距h而平行于坐标轴的两组平行线织成正方形网格，Δx=Δy=h，如图。; 我们将只考虑离开结点０充分近的那些结点，即（x-x0）充分小。于是可不计（x-x0）的三次及更高次幂的各项，则上式简写为：;联立(c),(d),解得差分公式： ; 以上(１-１)~(１-４)是基本差分公式,从而可导出其它的差分公式如下：; 差分公式(１-１)及(１-３)是以相隔2h的两结点处的函数值来表示中间结点处的一阶导数值，可称为中点导数公式。;第二节 应力函数的差分解; 一旦求得弹性体全部节点??φ值后，就可按应力分量差分公式（对节点0）算得弹性体各节点的应力。; 可见，用差分法解平面问题，共有两大任务： 一、建立差分方程 将（1-6）代入双调和方程; 二、联立求解这些线性代数方程，就能求得各内结点处的值。;代入上式，即得： ;关于边界上任一点处; 将此式亦从A点到B点沿s进行积分，就得到边界上任一点B处的φ值。为此利用分部积分法，得： ;将式(b),(c)代入，整理得：; 从图易看出，式（2－3）右边的积分式表示A与B之间的，x方向的面力之和；式（2－4）右边的积分式表示A与B之间的，y方向的面力之和；式（2－5）右边的积分式表示A与B之间的面力对于B点的矩。; 至此，我们解决了怎样 计算边界上各结点;所以有 ; 然后由面力的矩及面力之和算出边界上所有各结点处φ 的值，以及所必需的一些 及 值，即垂直于边界方向的导数值。;（5）按照公式（2－1）计算应力的分量。;第三节 例 深梁的应力函数差分解;(1）计算边界上各结点的φ 、 和 值。 取A为基点，且;（2）计算边界外以行各虚结点处的值。 由式（2-6）及前表可得;(3）边界内各结点的差分方程,由式（2-2）可知;联立求解上式，可得（以qh2为单位 );同理可得