弹性力学简明_第四版_徐芝纶第五章详解.ppt

§5-1 差分公式的推导 差分法是微分方程的一种数值解法。 它不是去求解函数 ，而是求函数在一些结点上的值 。 差分法的内容是： 导数差分公式的导出： 应用泰勒级数导出差分公式，可得出统一的格式，避免任意性，并可估计其误差量级，式 b 的误差为 。 对结点3， 得： 式 d 称为向后差分公式。 稳定温度场中的温度场函数T x, y 应满足下列方程和边界条件： （在 A 中）, a （在 上）, b （在 上）.　 c 现在我们将式 a 、 b 、 c 转化为差分形式。应用图5－1网格，和抛物线差分公式， （1）将　　　　化为差分公式，得 （2）若x边界516上为第一类边界条件，则 已知。 （3）若y边界627上为第二类边界条件，已 知　　，则 ?由于 所以得 这时，边界点2的　是未知的，对2点 须列出式（d）的方程。此方程涉及到 值，可将式（e）代入。 例2 稳定温度场问题的差 分解。设图中的矩形 域为6m×4m ，取网 格间距为h 2m，布 置网格如图，各边界 点的已知温度值如图 所示，试求内结点a 、b的稳定温度值。 解：对a、b 列出方程如下： 解出 1.比较导数的抛物线差分公式和线性差分公式的区别。 2.应用抛物线差分公式 5-2 ，试导出三阶导数 的差分公式。 §5-2 应力函数的差分解 （3）求出 后，由下式求应力（假设无体力）: 对边界内一行结点列式 e 方程时，需要求出边界点和边界外一行结点（虚结点）的 值。 为了求虚结点的 值，需要求出边界点 的 、 值。 3.应用应力边界条件 b ，求出边界点的 、 、 值。 ⑶∵A为定点， 、 和 、 、 均为常数，而式 h 中，加减x，y的一次式不影响应力，∴可取 故边界结点的 和导数值，由式 g 、 h 简化为 式 i 的物理意义是： 第一式表示从A到B边界上x向面力的主矢量； 第二式表示从A到B边界上y向面力的主矢量 改号; 第三式表示从A到B边界上面力对B点的力距， 图中以顺时针向为正。 因此，可以按物理意义直接求 ⑷ 由式 i 的第三式，可求出边界点的 值; 由式 i 的前两式，可求出边界点 的 、 值，然后再求出边 界外一行虚结点的 值。 思考题 1、将应力函数　看成是覆盖于区域A和边 界s上的一个曲面，则在边界上，各点 的　值与从 A（基点）到B面力的合力 距有关， 的一阶导数值与A到B的面力 的合力（主矢量）有关；而在区域内， 应力分量与　曲面的曲率、扭率有关。 §5－3应力函数差分解的实例 首先考虑对称性，可以减少未知值数目，并大量减少计算工作量。 按照物理意义，求出边界点上的 和其导数值（如书中所示）： 比较： 材料力学解─AM上 为直线分布， 弹性力学解─AM上 为曲线分布， 由此又说明，材料力学解法只适用于杆件。 （1）差分法是解微分方程边值问题和弹性 力学问题的有效方法。 （2）差分法简便易行，且总能求出解答。 （3）差分法可配合材料力学、结构力学解 法，精确地分析结构的局部应力状 态。 （1）对于曲线边界和不等间距网格的计算 较麻烦。 （2）差分法比较适用于平面问题或二维问 题。 （3）凡是近似解，在求导运算时会降低精 度。如 的误差为 ，则应力 的误差为 。 思考题： 1.试用线性向前或向后差分公式，导出 的 差分方程。 §5－4 弹性体的形变势能 外力势能 （3）对于平面应力问题 或平面应变问题 单元体积上应力所做的功都是 （4）假设没有转化为非机械能和动能，则 应力所做的功全部转化为弹性体的 内力势能，又称为形变势能，或应变 能， 存贮于物体内部。 ─单位体积的形变势 能（形变势能密度）。 （5）整个弹性体的形变势能是 3.形变势能 的性质 （1） 是应变或位移的二次泛函， 故不能应用叠加原理。 （2）应变或位移发生时， 总是正的，即 （3） 的大小与受力次序无关。 （4）对应变的导数，等于对应的应力： 4.弹性体的总势能，是外力势能和内力 （形变）势能之和， 1.试证明在线性的应力与应变关系下， 。 2. 试由式 e 导出式 g 。 3. 试列出极坐标系中平面应力问题的形变势能公式，并与式 d 、 e 和 f 相比较。 §5－5　位移变分方程 在位移变分法中，所取泛函为总势能 ，其自变量（宗量）为位移状态函数 ， 。 现在来导出位移变分方程。 2.虚位移状态 ⑴ 虚位移（数学上称为位移变分） ， 　表示在约束条件允许下，平衡状态附近的微小位移增量，如图所示。 虚位移应满足 上的