《弹性力学教学课件》差分法简介.pptVIP

* 弹性力学的经典解法存在一定的局限性，当弹性体的边界条件和受载情况复杂一点，往往无法求得偏微分方程的边值问题的解析解。因此，各种数值解法便具有重要的实际意义。差分法就是数值解法的一种。 所谓差分法，是把基本方程和边界条件（一般均为微分方程）近似地改用差分方程（代数方程）来表示，把求解微分方程的问题改换成为求解代数方程的问题。 * 差分公式的推导 我们在弹性体上，用相隔等间距h而平行于坐标轴的两组平行线织成正方形网格，如图。 设f=f(x,y)为弹性体内的某一个连续函数。该函数在平行于x轴的一根网线上，例如在3-０-１上，它只随x坐标的改变而变化。在邻近结点０处，函数f可展为泰勒级数如下： * 我们将只考虑离开结点０充分近的那些结点，即（x-x0）充分小。于是可不计（x-x0）的三次及更高次幂的各项，则上式简写为： 在结点３，x=x0-h；在结点1， x=x0+h。代入(b) 得： 联立（c）、（d），解得差分公式： * 同理，在网线4-0-2上可得到差分公式： 以上（１）—（４）是基本差分公式，从而可导出其它的差分公式如下： 图7-2 * 差分公式(１)及(３)是以相隔2h的两结点处的函数值来表示中间结点处的一阶导数值，可称为中点导数公式。 以相邻三结点处的函数值来表示一个端点处的一阶导数值，可称