弹性力学变分解法

1、关于弹性力学的变分解法第一张，PPT共七十六页，创作于2022年6月泛函的定义及举例函数：对于变量x的某一变域中的每一个值，y都有唯一一个值与之相对应，那么变量y称作变量x的函数。记为： y=f (x)x称为函数的自变量。泛函：对于某一类函数y()中的每一个函数y(x)，变量J都有一个值与之相对应，那么变量J称作依赖于函数y(x)的泛函。记为： J=J y(x)y(x)称为泛函的宗量。第二张，PPT共七十六页，创作于2022年6月例子 第三张，PPT共七十六页，创作于2022年6月第四张，PPT共七十六页，创作于2022年6月第五张，PPT共七十六页，创作于2022年6月第六张，PPT共七十六

2、页，创作于2022年6月第七张，PPT共七十六页，创作于2022年6月第八张，PPT共七十六页，创作于2022年6月泛函的极值泛函极值定理： 若可微泛函Jy(x)在y0(x)上达到极值，则在y= y0(x)上的变分为零。即第九张，PPT共七十六页，创作于2022年6月14.1 弹性体的虚功原理第十张，PPT共七十六页，创作于2022年6月第十一张，PPT共七十六页，创作于2022年6月第十二张，PPT共七十六页，创作于2022年6月静力可能的应力与几何可能的位移第十三张，PPT共七十六页，创作于2022年6月第十四张，PPT共七十六页，创作于2022年6月高斯公式第十五张，PPT共七十六页，创

3、作于2022年6月第十六张，PPT共七十六页，创作于2022年6月14.2 贝蒂互换定理第十七张，PPT共七十六页，创作于2022年6月第十八张，PPT共七十六页，创作于2022年6月第十九张，PPT共七十六页，创作于2022年6月第二十张，PPT共七十六页，创作于2022年6月14.3 位移变分方程 最小势能原理第二十一张，PPT共七十六页，创作于2022年6月第二十二张，PPT共七十六页，创作于2022年6月第二十三张，PPT共七十六页，创作于2022年6月第二十四张，PPT共七十六页，创作于2022年6月第二十五张，PPT共七十六页，创作于2022年6月第二十六张，PPT共七十六页，创作

4、于2022年6月第二十七张，PPT共七十六页，创作于2022年6月第二十八张，PPT共七十六页，创作于2022年6月第二十九张，PPT共七十六页，创作于2022年6月14.4 最小势能原理推导以位移表示的平衡微分方程及边界条件第三十张，PPT共七十六页，创作于2022年6月第三十一张，PPT共七十六页，创作于2022年6月第三十二张，PPT共七十六页，创作于2022年6月第三十三张，PPT共七十六页，创作于2022年6月第三十四张，PPT共七十六页，创作于2022年6月第三十五张，PPT共七十六页，创作于2022年6月第三十六张，PPT共七十六页，创作于2022年6月 一、里兹（Ritz）法基

5、本思想：设定位移函数的表达形式，使其满足位移边界条件，其中含有若干待定常数，然后利用最小势能原理(位移变分方程)确定这些常数，即得位移解。设选取的位移表达式如下：（a）其中：为互不相关的 3m 个系数；为设定的函数，且在边界上有：为边界上取零值的设定函数 显然，上述函数满足位移边界条件。此时，位移的变分由系数 Am、Bm、 Cm的变分来实现。与变分无关。14.5 基于最小势能原理的近似计算方法第三十七张，PPT共七十六页，创作于2022年6月（b）位移的变分：形变势能的变分：由应变能计算式可知：（c）根据最小原理，有：第三十八张，PPT共七十六页，创作于2022年6月将上式整理、移项、合并，可

6、得：完全任意，且互相独立，要使上式成立，则须有：第三十九张，PPT共七十六页，创作于2022年6月 Ritz 法方程或称 Rayleigh- Ritz 法方程说明：（1）由 U 的表达式（4-20）可知，U 是系数的二次函数，因而，上式为各系数的线性方程 组。互不相关，因而，总可以求出全部的系数。（2）求出了系数就可求得其它量，如位移、应力等（3）在假定位移函数时，须保证其满足全部位移边界条件。第四十张，PPT共七十六页，创作于2022年6月第四十一张，PPT共七十六页，创作于2022年6月第四十二张，PPT共七十六页，创作于2022年6月第四十三张，PPT共七十六页，创作于2022年6月第四

7、十四张，PPT共七十六页，创作于2022年6月第四十五张，PPT共七十六页，创作于2022年6月第四十六张，PPT共七十六页，创作于2022年6月基本思想构造位移试函数满足位移(面力)边界条件RayleighRitz(瑞利里兹)法(伽辽金)法 通过能量变分，偏微分方程边值问题转化为线性代数方程组。位移边界条件位移与面力边界条件第四十七张，PPT共七十六页，创作于2022年6月解：用瑞利里兹法位移试函数 例1： 两端简支的等截面梁，受均匀分布载荷q作用如图所示,不计体力。试求解梁的挠度w(x) 满足梁的位移边界条件：在x=0，l处，w=0 简支梁的形变势能为： 第四十八张，PPT共七十六页，创作

8、于2022年6月积分后可得:外力势能为:当m为奇数时。第四十九张，PPT共七十六页，创作于2022年6月有：所以回代 第五十张，PPT共七十六页，创作于2022年6月挠曲线表达式是无穷级数精确解这个级数收敛很快，只要取少数几项就可以得到足够的精度。如果取一项 这一结果与精确值十分接近 。第五十一张，PPT共七十六页，创作于2022年6月例2:如图所示简支梁，中点处承受有集中P，试求的梁的挠曲线方程。 PABlxy解:（1）假设位移试探函数（必须满足位移边界条件）设位移试探函数为（取一项）：式中：a 为待定常数。（2）计算形变势能 U：（ a）（ b）显然，式（a）满足端点的位移边界条件:（3）

9、代入Ritz 法方程，求解（ c）（ d）第五十二张，PPT共七十六页，创作于2022年6月PABlxy讨论：（1）中点的挠度：（ e）而材料力学的结果：两者比较：式（a）的结果偏小1.46%。如果取如下位移函数：式中项数 m 取得越多，则求得精度就越高。（2）所取的位移函数必须满足位移边界条件。（3）位移函数选取不是唯一的，如：第五十三张，PPT共七十六页，创作于2022年6月PABlxy例3:如图所示简支梁，中点处承受有集中P，试求的梁的挠曲线方程。 解:（1）假设位移试探函数式中：A1、A2 为待定常数。显然，式（a）满足端点的位移边界条件:（2）计算：梁的形变势能:（3）代入 Ritz

10、 法方程:第五十四张，PPT共七十六页，创作于2022年6月PABlxy例3:如图所示简支梁，中点处承受有集中P，试求的梁的挠曲线方程。 解:位移函数（a）（3）代入 Ritz 法方程:所求挠曲线方程 :第五十五张，PPT共七十六页，创作于2022年6月PABlxy所求挠曲线方程:中点挠度:而材料力学的结果：第五十六张，PPT共七十六页，创作于2022年6月解：位移试函数 例4：矩形薄板，四边固定，受平行于板面的体力作用。设坐标轴如图所示，试用RayleighRitz法求解。 m和n为正整数在边界x=0，a，和y=0，b上，u=v=0，所以试函数满足位移边界条件。 第五十七张，PPT共七十六页

11、，创作于2022年6月平面应力问题 因此 第五十八张，PPT共七十六页，创作于2022年6月将位移试函数代入求导数后再积分 因此 如果体力已知，积分可求待定系数Amn和Bmn第五十九张，PPT共七十六页，创作于2022年6月例5：图示薄板，宽为 a，高度为 b，左边和下边受连杆支承，右边和上边分别受有均布压力 q1和 q2 作用，不计体力。试求薄板的位移。解：（1）假设位移函数（a）满足边界条件：在式（a）中u,v 各取一项 ，即（b）（2）计算形变势能 U将式（b）代入平面应力情形下形变势能公式，有积分得：（c）第六十张，PPT共七十六页，创作于2022年6月（c）（3）代入Ritz 法方程

12、求解体力有在右边界：在上边界：于是有：将式（c）代入，得（11-15）第六十一张，PPT共七十六页，创作于2022年6月联立求解，得：（f）代入位移表达式（b），得：（g）讨论：（1）如果在位移式（a）中再多取一些系数如：A2、B2等，但是经计算，这些系数全为零。（2）位移解（g）满足几何方程、平衡方程和边界条件。表明：位移解（g）为问题的精确解。第六十二张，PPT共七十六页，创作于2022年6月Ritz 法解题步骤：（1）假设位移函数，使其满足边界条件；（2） 计算形变势能 U ；（3）代入Ritz 法方程求解待定系数；（4）回代求解位移、应力等。第六十三张，PPT共七十六页，创作于2022

13、年6月例6:图示矩形薄板，宽为2 a，高度为2 b，左右两边和下边均被固定，而上边的给定位移为：（a）不计体力。试求薄板的位移和应力。解:（1）假设位移函数只取一项，即 m =1, 将位移分量设为：（b）显然，可满足位移边界条件：第六十四张，PPT共七十六页，创作于2022年6月（2）代入Galerkin 法方程求解该问题中无应力边界条件，式（b）满足全部条件。可用伽辽金（Galerkin）法求解。X = Y = 0，m = 1，伽辽金法方程变为：（c）第六十五张，PPT共七十六页，创作于2022年6月 将其代入伽辽金方程（c）, 可求得：第六十六张，PPT共七十六页，创作于2022年6月代回位移表达式（b）, 得位移解答：当 b = a，取 = 0.2时，上述解答成为：第六十七张，PPT共七十六页，创作于2022年6月（3）求应力分量应用几何方程及物理方程，可求得应力为：第六十八张，PPT共七十六页，创作于2022年6月第六十九张，PPT共七十六页，创作于2022年6月由虚功原理：应力变分方程（虚应力方程）第七十张，PPT共七十六页，创作于2022年6月由广义虚功原理：第七十一张，PPT共七十六页，创作于2022年6月表明