有限差分法弹性力学

有限差分法弹性力学第一页，共三十页，编辑于2023年，星期三

工程中常用的数值解法有有限单元法和差分法。

有限单元法

是以有限个单元的集合体来代替连续体，属于物理上的近似。

差分法

是把弹性力学的基本方程和边界条件（一般均为微分方程）近似地改用差分方程（代数方程）来表示，把求解微分方程的问题改换成为求解代数方程的问题，属于数学上的近似。第二页，共三十页，编辑于2023年，星期三

第一节差分方程

第二节应力函数的差分解

第三节深梁应力函数的差分解第三页，共三十页，编辑于2023年，星期三第一节差分方程

差分法是沿用已久的一种数值解法。随着计算机的普及和相应的软件发展，此法成为解弹性力学问题的一种有效的方法。第四页，共三十页，编辑于2023年，星期三

我们在弹性体上，用相隔等间距h而平行于坐标轴的两组平行线织成正方形网格，Δx=Δy=h，如图。

设f=f(x,y)为弹性体内的某一个连续函数。该函数在平行于x轴的一根网线上，如在３-０-１上,它只随x坐标的改变而变化。在邻近结点０处，函数f可展为泰勒级数如下：第五页，共三十页，编辑于2023年，星期三

我们将只考虑离开结点０充分近的那些结点，即（x-x0）充分小。于是可不计（x-x0）的三次及更高次幂的各项，则上式简写为：在结点３，x=x0-h,在结点1,x=x0+h，代入(b)得：第六页，共三十页，编辑于2023年，星期三联立(c),(d),解得差分公式：

同理，在网线４-０-２上可得到差分公式第七页，共三十页，编辑于2023年，星期三

差分公式(１-１)及(１-３)是以相隔2h的两结点处的函数值来表示中间结点处的一阶导数值，可称为中点导数公式。

以相邻三结点处的函数值来表示一个端点处的一阶导数值，可称为端点导数公式。

应当指出：中点导数公式与端点导数公式相比，精度较高。因为前者反映了结点两边的函数变化，而后者却只反映了结点一边的函数变化。因此，我们总是尽可能应用前者，而只有在无法应用前者时才不得不应用后者。第八页，共三十页，编辑于2023年，星期三

以上(１-１)~(１-４)是基本差分公式,从而可导出其它的差分公式如下：第九页，共三十页，编辑于2023年，星期三第二节应力函数的差分解

当不计体力时，我们已把弹性力学平面问题归结为在给定边界条件下求解双调和方程的问题。用差分法解平面问题，就应先将双调和方程变换为差分方程，而后求解之。第十页，共三十页，编辑于2023年，星期三

一旦求得弹性体全部节点的φ值后，就可按应力分量差分公式（对节点0）算得弹性体各节点的应力。第十一页，共三十页，编辑于2023年，星期三

可见，用差分法解平面问题，共有两大任务：一、建立差分方程将（1-6~8）代入双调和方程

对于弹性体边界以内的每一结点，都可以建立这样一个差分方程。整理即得第十二页，共三十页，编辑于2023年，星期三

二、联立求解这些线性代数方程，就能求得各内结点处的值。

为了求得边界上各结点处的φ值，须要应用应力边界条件，即：

一般建立和求解差分方程，在数学上不会遇到很大困难。但是，当对于边界内一行的（距边界为h的）结点，建立的差分方程还将涉及边界上各结点处的φ值，并包含边界外一行的虚结点处的φ值。第十三页，共三十页，编辑于2023年，星期三代入上式，即得：

l1=cos(N,x)=cosα=dy/ds,l2=cos(N,y)=sinα=-dx/ds,于是，式（a）可改写为：由右图可见，第十四页，共三十页，编辑于2023年，星期三关于边界上任一点处由此得：

的值，可将（b）式从A点到B点对s积分得到：第十五页，共三十页，编辑于2023年，星期三

将此式亦从A点到B点沿s进行积分，就得到边界上任一点B处的φ值。为此利用分部积分法，得：

由高等数学可知，第十六页，共三十页，编辑于2023年，星期三将式(b),(c)代入，整理得：

由前知，把应力函数加上一个线性函数，并不影响应力。因此，可设想把应力函数加上a+bx+cy，然后调整a，b，c三个数值，使得由式(d)及式(c)可见，设即可根据面力分量及导数求得为已知，第十七页，共三十页，编辑于2023年，星期三

从图易看出，式（2－3）右边的积分式表示A与B之间的x方向的面力之和；式（2－4）右边的积分式表示A与B之间的y方向的面力之和；式（2－5）右边的积分式表示A与B之间的面力对于B点的矩。于是式(d)，式(c)即简化为：第十八页，共三十页，编辑于2023年，星期三

至此，我们解决了怎样计算边界上各结点的值的问题。

至于边界外一行虚结点处的值，则可用边界上结点处的或值和边界内一行相应结点处的值来表示。例如，对于图1中的虚结点14，因为有

第十九页，共三十页，编辑于2023年，星期三所以有

当求出全部结点上的φ值以后，我们就可按应力分量的差分公式（2－1）计算应力分量。

用差分法解弹性平面问题时，可按下列步骤进行：（1）在边界上任意选定一个结点作为基点A，取第二十页，共三十页，编辑于2023年，星期三

然后由面力的矩及面力之和算出边界上所有各结点处φ

的值，以及所必需的一些及值，即垂直于边界方向的导数值。

（2）应用公式（2－6），将边界外一行虚结点处的φ值用边界内的相应结点处的φ值来表示。

（3）对边界内的各结点建立差分方程（2－2），联立求解这些结点处的值。

（4）按照公式（2－6），算出边界外一行的各虚结点处的φ值。第二十一页，共三十页，编辑于2023年，星期三（5）按照公式（2－1）计算应力的分量。说明：

1．以上是针对单连体导出的结果。对于多连体，情况就不象这样简单。

2.

如果一部分边界是曲线的，或是不与坐标轴正交，则边界附近将出现不规则的内结点。对于这样的结点，差分方程（2－2）必须加以修正。第二十二页，共三十页，编辑于2023年，星期三第三节例深梁的应力函数差分解

现以如图所示的混凝土深梁为例，应用应力函数的差分解求出应力分量。已知混凝土深梁上边受有均布向下的铅直荷载q，并由下角点处的反力维持平衡。第二十三页，共三十页，编辑于2023年，星期三(1）计算边界上各结点的φ

、和值。取A为基点，且由上面公式所得的计算结果见下表。第二十四页，共三十页，编辑于2023年，星期三（2）计算边界外以行各虚结点处的值。由式（2-6）及前表可得第二十五页，共三十页，编辑于2023年，星期三(3）边界内各结点的差分方程,由式（2-2）可知第二十六页，共三十页，编辑于2023年，星期三联立求解上式，可得（以qh2为单位)

（4）计算结点外一行各结点处的值。由(a)、(b)、(c)可得（5）计算应力。对于结点M，由式（1）可知

第二十七页，共三十页，编辑于2023年，星期三同理可得

沿着梁的中线MA，的变化如下图和右图所示。

第二十八页，共三十页，编辑于2023年，星期三差分方法的优缺点优点：

1.差分法是解微分方程边值问题和弹性力学问题的有效方法；

2.差分法简便易行；

3.对于某些结构，为了更精确地分析局部的应力状态，可以用差分法进行分析；缺点：

1.对于曲线边界和斜边界等产生的不等间距网格的