弹性力学用差分法与变分法解平面问题.pptVIP

《弹性力学》 第五章 用差分法和变分法解平面问题  第五章用差分法和变分法解平面问题 内容提要 、差分公式的推导 弹性体的形变势能和外力势能 三、位移变分方程 四、位移变分法 律怎 芝09五、位移变分法例题 弹性力学简明教程 弹性力学简明教程(第三版  差分公式的推导 弹性力学的基本解法是,根据静力平衡条件,形变与位移 之间的几何条件和形变与应力之间的物理条件,建立微分方程 和边界条件 因此,弹性力学问题属于微分方程的边值问题。通过求解, 得出函数表示的精确解答。 对于工程实际问题,由于荷载和边界较复杂,难以求出函 数式的解答。为此,人们探讨弹性力学的各种近似解法,主要 有差分法、变分法和有限单元法。 弹性力学 用差分法和变分法解平面问题  差分公式的推导 差分法是微分方程的一种数值解法,它不是去求解函数f(x), 而是求函数在一些结点上的值f,F2 将微分用有限差分来代替 dx→)△x=x f(x) df→yf=f2-f 将导数用有限差商来代替 , J df、f2-f1 dx△xxn-x O 将微分方程用差分方程(代数方程)代替,求解微分方程问 题化为求解差分方程问题。 弹性力学 用差分法和变分法解平面问题  差分公式的推导 在弹性体上,用相隔等间距h而平 行于坐标轴的两组平行线织成正方形 网格。 设f=f(x,y)为弹性体内的某一个连续 函数,该函数在平行于x轴的一根网线 上,如在3-0-1上,它只随x坐标的改 变而变化。在邻近结点0处,函数f可 展为泰勒级数如下 f=后+(x-x0)+ 2(-(x-5) 弹性力学 用差分法和变分法解平面问题 5  差分公式的推导 6f6-4(f+f3)+(f+f1 ax 4 14f6-2(f1+f2+3+f4)+(f+f6+f1+f 6f6-4(2+f4)+(1o+f2 差分公式(1)及(3)是以相隔2h的两结点处的函数值来表示 中间结点处的一阶导数值,可称为中点导数公式。 以相邻三结点处的函数值来表示一个端点处的一阶导数值, 可称为端点导数公式 应当指出:中点导数公式与端点导数公式相比,精度较高。 因为前者反映了结点两边的函数变化,而后者却只反映了结点 边的函数变化。因此,我们总是尽可能应用前者,而只有在 无法应用前者时才不得不应用后者 弹性力学 用差分法和变分法解平面问题  弹性体的形变势能和外力势能 变分法,是研究泛函及其极值的求解方法。 泛函一一是以函数为自变量的一类函数。 弹性力学变分法,因其泛函就是弹性体的能量(如形变势能、 外力势能),又称为能量法。 弹性力学变分法,是区别于微分方程边值问题的另一种独 立解法,分为 位移变分法:取位移函数为自变量,并以势能极小值条 件导出变分方程。 应力变分法:取应力函数为自变量,并以余能极小值条 件导出变分方程 弹性力学 用差分法和变分法解平面问题 10