第5章差分法及变分法解平面问题概论.ppt

§5-1 差分公式的推导 §5-3 应力函数差分解的实例 * * 第一节 差分公式的推导 第二节 应力函数的差分解 第三节 应力函数差分解的实例 第四节 弹性体的形变势能 第五节 位移变分方程 第六节 位移变分法 第七节 位移变分法应用于平面问题 图5-1 工程中许多重要的实际问题，并不能得出函数式解答，必须进行数值计算法求解，差分法是数值解法之一。 差分法是把基本方程和边界条件（一般都是微分方程）近似地改用差分方程（代数方程）来表示，把求解微分方程的问题转化为求解代数方程的问题。 首先把弹性体用间距为h的两组平行线织成网格，如图5-1。设函数 是弹性体内的某一连续函数，将函数在平行于轴的网线上，如在3—0—1上，在临近结点0处，展开成台勒级数： （a） 在结点3及结点1，x分别等于 及 ，即将 及 代入（b）式有： （b） （5-2） （5-1） 同理，可以得到： （5-3） （5-4） 公式（5-1）至（5-4）称为基本差分公式，其它差分公式可以由基本差分公式得到： （5-5） 同理，可以得到： （5-6） §5-2 应力函数的差分解 1. 应力函数的差分解 当不计体力时，平面问题中的应力分量为： ， ， （a） 把任一结点0处的应力分量表示成： （5-11） 可见，若已知了结点处的F值，就可以求出各结点处的应力值。 将差分公式（5-6）代入相容方程，即 即可得出： （5-12） 对于弹性体内的每一点，都可建立一个这样的差分方程，但对边界内一行的结点（距边界为h），建立的差分方程中还将包含边界上各结点处的F值，并包含边界外一行的虚结点处的F值。 由应力边界条件（2-18），有： ， 利用式（a），将其转化为： (b) 图5-2 由图（5-2）可知， ， 将（b）式改写为： ， （c） 将（c）式从A点到B点对s积分，得： ， （d） 因为 ，故由分步积分法可得： 将（c）式代入得： 再将式（d）代入，得： （e） 由式（e）及式（d）可见，若已知 、 、 ，即可由面力分量 、 求得于 、 、 。取 = = ，于是式（e）及式（d）简化为： （5-13） （5-14） （5-15） 对多连体情况，只能直接应用式（e）及式（d），不能应用简化式（5-13）至（5-15）。这就使得应力函数的差分解在多连体问题中应用不方便。 式（5-13）右边的积分表示A与B之间的x方向面力之和；式（5-14）右边的积分表示A与B之间的y方向面力之和；式（5-15）右边的积分表示A与B之间的面力对B点的矩（在x轴向右y轴向下的坐标系下，力矩以顺时针转向为正）。 边界外一行的（距边界为h的）虚结点处的F值，可用边界内一行结点处的F值及边界上的导数值来表示。如虚结点13及14 ： ， ， （5-16） 2. 差分法的解题步骤 （1）在边界上任意选定一个结点作为基点A ，取 ，再由面力的矩及面力之和算出边界上所有各结点处的F值，以及应用式（5-16）时所必需的一些 值及 值。 （2）应用式（5-16），将边界外一行各虚结点处的F值用边界内的相应结点处的F值来表示。 （3）对边界内的各结点建立差分方程（5-12），联立求解这些结点处的F值。 （4）按式（5-16），算出边界外一行各虚结点处的F值。 （5）按式（5-11）计算应力分量。 设有正方形的砼深梁，图5-3，上面受均布压力q，由下角点处的力平衡，用应力函数的差分法求出应力分量。 取坐标如图，网格划分h=1/6边长，由于对称，只计算一半梁。 图5-3 （1）选基点A，取 ，计算边界各结点处的F值及 所需的 值及 值。 0 0 0 0 0 0 0 0 / / / 0 0 / / / / / / 0 M L K J E、F、G、H、I D B、C A 结点 （2）将边界外一行各虚结点处的F值用边界内的相应结点处的F值表示出来。 上、下边， =0，有： （a） 左边界，