第五章_有限差分法(2008_11改).ppt

所以有 当求出全部结点上的φ值以后，我们就可按应力分量的差分公式（2－1）计算应力分量。 用差分法解弹性平面问题时，可按下列步骤进行： （1） 在边界上任意选定一个结点作为基点A， 取 然后由面力的矩及面力之和算出边界上所有各结点处φ 的值，以及所必需的一些 及 值，即垂直于边界方向的导数值。 （2）应用公式（2－6），将边界外一行虚结点处的φ值用边界内的相应结点处的φ值来表示。 （3）对边界内的各结点建立差分方程（2－2），联立求解这些结点处的值。 （4）按照公式（2－6），算出边界外一行的各虚结点处的φ值。 （5）按照公式（2－1）计算应力的分量。 说明： 1．以上是针对单连体导出的结果。对于多连体，情况就不象这样简单。 2. 如果一部分边界是曲线的，或是不与坐标轴正交，则边界附近将出现不规则的内结点。对于这样的结点，差分方程（2－2）必须加以修正。 第三节 例 深梁的应力函数差分解 现以如图所示的混凝土深梁为例，应用应力函数的差分解求出应力分量。已知混凝土深梁上边受有均布向下的铅直荷载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 第五章 有限差分法 从弹性力学的基本方程建立以来，这些方程在各种问题的边界条件下如何求解，一直是很多数学工作者和力学工作者研究的内容。即弹性力学的经典解法存在一定的局限性，当弹性体的边界条件和受载情况复杂一点，往往无法求得偏微分方程的边值问题的解析解，许多工程重要问题，不能够得出函数式的解答。 因此，弹性力学问题的各种数值解法便具有重要的实际意义。 工程中常用的数值解法有有限单元法和差分法。 有限单元法 是以有限个单元的集合体来代替连续体，属于物理上的近似。 差分法 是把弹性力学的基本方程和边界条件（一般均为微分方程）近似地改用差分方程（代数方程）来表示，把求解微分方程的问题改换成为求解代数方程的问题，属于数学上的近似。 第一节 差分方程 第二节 应力函数的差分解 第三节 深梁应力函数的差分解 第一节 差分方程 差分法是沿用已久的一种数值解法。随着计算机的普及和相应的软件发展，此法成为解弹性力学问题的一种有效的方法。 我们在弹性体上，用相隔等间距h而平行于坐标轴的两组平行线织成正方形网格，Δx=Δy=h，如图。 设f=f(x,y)为弹性体内的某一个连续函数。该函数在平行于x轴的一根网线上，如在３-０-１上,它只随x坐标的改变而变化。在邻近结点０处，函数f可展为泰勒级数如下： 我们将只考虑离开结点０充分近的那些结点，即（x-x0）充分小。于是可不计（x-x0）的三次及更高次幂的各项，则上式简写为： 在结点３，x=x0-h, 在结点1, x=x0+h，代入(b) 得： 联立(c),(d),解得差分公式： 同理，在网线４-０-２上可得到差分公式 差分公式(１-１)及(１-３)是以相隔2h的两结点处的函数值来表示中间结点处的一阶导数值，可称为中点导数公式。 以相邻三结点处的函数值来表示一个端点处的一阶导数值，可称为端点导数公式。 应当指出：中点导数公式与端点导数公式相比，精度较高。因为前者反映了结点两边的函数变化，而后者却只反映了结点一边的函数变化。因此，我们总是尽可能应用前者，而只有在无法应用前者时才不得不应用后者。 以上(１-１)~(１-４)