能量原理与变分法精选PPT.ppt

《能量原理与变分法精选PPT.ppt》由会员分享，可在线阅读，更多相关《能量原理与变分法精选PPT.ppt（123页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、关于能量原理与变分法关于能量原理与变分法第1页，讲稿共123张，创作于星期三11-1 11-1 弹性体的形变势能弹性体的形变势能主主主主 要要要要 内内内内 容容容容 11-2 11-2 位移变分方程位移变分方程11-3 11-3 位移变分法位移变分法11-4 11-4 位移变分法应用于平面问题位移变分法应用于平面问题11-5 11-5 应力变分方程应力变分方程11-6 11-6 应力变分法应力变分法11-7 11-7 应力变分法应于平面问题应力变分法应于平面问题11-8 11-8 应力变分法应于扭转问题应力变分法应于扭转问题11-9 11-9 解答的唯一性解答的唯一性11-10 11-10

2、功的互等定理功的互等定理第2页，讲稿共123张，创作于星期三11-0 11-0 引引 言言1.弹性力学问题的弹性力学问题的微分提法微分提法及其及其解法解法：（1）平衡微分方程）平衡微分方程（2）几何方程）几何方程（3）物理方程）物理方程（4）边界条件）边界条件应力边界条件；应力边界条件；位移边界条件；位移边界条件；定定解解问问题题求解方法求解方法：（1）按位移求解）按位移求解基本方程：基本方程：（a）以位移为基本未知量）以位移为基本未知量的的平衡微分方程平衡微分方程;（2）按应力求解）按应力求解基本方程：基本方程：（a）平衡微分方程；）平衡微分方程；（b）边界条件。）边界条件。（b）相容方程；

3、相容方程；（c）边界条件。边界条件。（a）归结为归结为求解联立的微分方程求解联立的微分方程组组；求解特点：求解特点：（b）难以求得难以求得解析解解析解。从研究从研究微小单元微小单元体入手，考察其体入手，考察其平衡平衡、变形变形、材料性质材料性质，建立基本方程：，建立基本方程：第3页，讲稿共123张，创作于星期三2.弹性力学问题的弹性力学问题的变分提法变分提法及其及其解法解法：基本思想基本思想：在在所有可能的解所有可能的解中，求出最接近于精确解的解；中，求出最接近于精确解的解；将定解问题转变为将定解问题转变为求解线性方程组求解线性方程组。弹性力学中的变分原理弹性力学中的变分原理 能量原理能量原理

4、 直接处理直接处理整个弹性系统整个弹性系统，考虑系统的，考虑系统的能量关系能量关系，建立一些泛函的，建立一些泛函的变分方程变分方程，将弹性力学问题归结为将弹性力学问题归结为在给定约束条件下求泛函极（驻）值的变分问题在给定约束条件下求泛函极（驻）值的变分问题。（变分解法也称（变分解法也称能量法能量法）（a）以）以位移位移为基本未知量，为基本未知量，得到得到最小势（位）能原理最小势（位）能原理等。等。（b）以）以应力应力为基本未知量，为基本未知量，得到得到最小余能原理最小余能原理等。等。（c）同时以）同时以位移、应力、应变位移、应力、应变为未知量，为未知量，得到得到广义（约束）变分原理。广义（约束

5、）变分原理。位移法位移法 力法力法 混合法混合法 有限单元法有限单元法、边界元法边界元法、离散元法离散元法 等等数值解法数值解法的理论基础。的理论基础。求解方法求解方法：里兹（里兹（Ritz）法，）法，伽辽金（伽辽金（Galerkin）法，）法，加权残值（加权残值（余量）法等。余量）法等。第4页，讲稿共123张，创作于星期三3.弹性力学问题的弹性力学问题的数值解法数值解法：（a）直接求解联立的微分方程组（弹性力学的基本方程）直接求解联立的微分方程组（弹性力学的基本方程）有限差分法；有限差分法；基本思想：基本思想：将将导数导数运算近似地用运算近似地用差分差分运算代替；运算代替；将定解问题转变为将

6、定解问题转变为求解线性方程组求解线性方程组。典型软件：典型软件：FLAC实质：实质：将将变量离散变量离散。（b）对）对变分方程变分方程进行数值求解进行数值求解 有限单元法有限单元法、边界元法边界元法、离散元法离散元法 等等典型软件：典型软件：ANSYS，MARC，ADINA，SAP，NASTRAN，ABAQUS 等；等；基于有限元法的分析软件；基于有限元法的分析软件；UDEC 基于离散元法的分析软件；基于离散元法的分析软件；基本思想：基本思想：将求解将求解区域离散区域离散，离散成有限个小区域（离散成有限个小区域（单元单元），），在小区域（单元）上假设在小区域（单元）上假设可能解可能解，最后由能

7、量原理最后由能量原理（变分原理）确定其最优解。（变分原理）确定其最优解。将问题转变为将问题转变为求解求解大型大型的线性方程组的线性方程组。第5页，讲稿共123张，创作于星期三11-1 11-1 弹性体的形变势能弹性体的形变势能1.形变势能的一般表达式形变势能的一般表达式Pxl0l单向拉伸：单向拉伸：PlOPl外力所做的功：外力所做的功：由于在静载（缓慢加载）条件下，由于在静载（缓慢加载）条件下，其它能量损失很小，所外力功全部转化其它能量损失很小，所外力功全部转化杆件的杆件的形变势能（形变势能（变形能变形能）U：杆件的体积杆件的体积令：令：单位体积的变形能，单位体积的变形能，称为称为比能比能。三

8、向应力状态：三向应力状态：一点的应力状态：一点的应力状态：xyz第6页，讲稿共123张，创作于星期三三向应力状态：三向应力状态：一点的应力状态：一点的应力状态：xyz 由能量守恒原理，形变势能的值与弹性体受力的由能量守恒原理，形变势能的值与弹性体受力的次序无关次序无关，只，只取决于最终的状态。取决于最终的状态。假定所有应力分量与应变分量全部按同样的比例增加，此时，单元假定所有应力分量与应变分量全部按同样的比例增加，此时，单元体的体的形变比能形变比能：（a）整个弹性体的整个弹性体的形变势能：形变势能：（b）（c）若用张量表示：若用张量表示：形变比能：形变比能：整体形变势能：整体形变势能：第7页，

9、讲稿共123张，创作于星期三2.形变势能的应力分量表示形变势能的应力分量表示在线弹性的情况下，由物理方程（在线弹性的情况下，由物理方程（8-17）：代入式（代入式（a），整理得形变势能的表达式：），整理得形变势能的表达式：（d）（e）代入式（代入式（b），有：），有：（11-1）将式（将式（e）分别对）分别对6 个应力分量求导，并将其结果与物理方程比较，得：个应力分量求导，并将其结果与物理方程比较，得：第8页，讲稿共123张，创作于星期三（11-2）表明：表明：弹性体的比能对于任一应力分量的改变率，就等于相应的形变分量。弹性体的比能对于任一应力分量的改变率，就等于相应的形变分量。3.形变势能的

10、应变分量表示形变势能的应变分量表示用应变表示的物理方程（用应变表示的物理方程（8-19）：）：（f）第9页，讲稿共123张，创作于星期三或：或：代入式（代入式（a）：）：（a）并整理可得：并整理可得：（g）（11-3）0 1/2，U 0 即弹性体的形变势能是非负的量。即弹性体的形变势能是非负的量。将上式对将上式对6个应变分量分别求导，再与应力表示的物理方程（个应变分量分别求导，再与应力表示的物理方程（8-17）比较，）比较，可得：可得：第10页，讲稿共123张，创作于星期三（11-4）将几何方程（将几何方程（8-9）代入上式，得：）代入上式，得：弹性体的比能对于任一应变分量的改变率，就等于相应

11、的应力分量。弹性体的比能对于任一应变分量的改变率，就等于相应的应力分量。格林公式格林公式4.形变势能的位移分量表示形变势能的位移分量表示表明：表明：（11-5）第11页，讲稿共123张，创作于星期三11-2 11-2 位移变分方程位移变分方程1.泛函与变分的概念泛函与变分的概念（1）泛函的概念）泛函的概念函数：函数：x 自变量；自变量；y 因变量，或称自变量因变量，或称自变量 x 的函数。的函数。泛函：泛函：x 自变量；自变量；y 为一变函数；为一变函数；F 为函数为函数 y 的函数，的函数，称为称为泛函泛函。例例1：P1 弯矩方程弯矩方程梁的形变势能：梁的形变势能：ABlx 泛函泛函例例2：

12、第12页，讲稿共123张，创作于星期三例例2：因为因为所以，所以，U 被称为被称为形变势能泛函形变势能泛函。（2）变分与变分法）变分与变分法设：设：当自变量当自变量 x 有一增量：有一增量：函数函数 y 也有一增量：也有一增量：dy 与与 dx，分别称为自变量，分别称为自变量 x 与函数与函数 y 的的 微分。微分。微分问题微分问题P1ABlx设：设：函数函数 y 有一增量：有一增量：泛函泛函 U 也有一增量：也有一增量：第13页，讲稿共123张，创作于星期三P1ABlx设：设：函数函数 y 也有一增量：也有一增量：泛函泛函 U 也有一增量：也有一增量：函数的增量函数的增量y、泛函的增量、泛函

13、的增量 U 等称为变等称为变分分。研究研究自变函数的增量自变函数的增量与与泛函的增量泛函的增量 间关系称为间关系称为变变分问题分问题。例如：例如：Pcr（1）压杆稳定问题）压杆稳定问题 寻求压杆形变势能寻求压杆形变势能 U 达到达到最大值最大值时的压力时的压力 P 值。值。（2）球下落问题）球下落问题12 球从球从位置位置1下落下落至至位置位置2，所需时间，所需时间为为T，当当 最速下降问题最速下降问题 泛函的变分问题泛函的变分问题第14页，讲稿共123张，创作于星期三（3）变分及其性质）变分及其性质定义：定义：泛函泛函增量：增量：函数函数连续性：连续性：称函数称函数 z 在在 x0 点连续。

14、点连续。当当有有称泛函称泛函 U 在在 y0(x)处零阶接近。处零阶接近。当当有有称泛函称泛函 U 在在 y0(x)处一阶接近。处一阶接近。当当有有称泛函称泛函 U 在在 y0(x)处二阶接近。处二阶接近。第15页，讲稿共123张，创作于星期三泛函泛函函数函数微分：微分：当当x0时，时，0，则，则 z 可用其可用其线性主部表示其微分。即线性主部表示其微分。即 U 增量的线性主部增量的线性主部变分：变分：当当 max|y|0时，时，max 0，则，则 U 可用其线性主部表示可用其线性主部表示,即即极值：极值：若若在在 x0 处有极值，处有极值，则有：则有：若若 Uy(x)在在 y0(x)处有极值

15、，处有极值，条件：条件：一阶变分为零一阶变分为零。当当取得极值取得极值 称为称为强极值强极值当当取得极值取得极值 称为称为弱极值弱极值极值：极值：第16页，讲稿共123张，创作于星期三（4）变分的运算）变分的运算变分与微分运算：变分与微分运算：变分运算与微分运算互相交换变分运算与微分运算互相交换。变分与积分运算：变分与积分运算：变分运算与积分运算互相交换变分运算与积分运算互相交换。复合函数的变分：复合函数的变分：其中：其中：一阶变分：一阶变分：第17页，讲稿共123张，创作于星期三复合函数的变分：复合函数的变分：其中：其中：一阶变分：一阶变分：自变量自变量 x 的变分的变分 x 0二阶变分：二

16、阶变分：二阶变分用于判别驻值点是取得二阶变分用于判别驻值点是取得极大值极大值还是还是极小值极小值。第18页，讲稿共123张，创作于星期三2.位移变分方程位移变分方程建立：弹性体的建立：弹性体的形变势能形变势能与与位移位移间间变分变分关系关系 位移变分方程位移变分方程qP应力边界应力边界 S位移边界位移边界 Su设弹性体在外力作用下，处于平衡状态。设弹性体在外力作用下，处于平衡状态。边界：边界：位移场：位移场：应力场：应力场：满足：平衡方程、几何满足：平衡方程、几何方程、物理方方程、物理方程、边界条件。程、边界条件。称为称为真实解真实解（1）任给弹性体一微小的位移变化：）任给弹性体一微小的位移变

17、化：满足两个条件：满足两个条件：（1）不破坏平衡状态；）不破坏平衡状态；（2）不破坏约束条件，即为约束所允许。）不破坏约束条件，即为约束所允许。第19页，讲稿共123张，创作于星期三任给弹性体一微小的位移变化：任给弹性体一微小的位移变化：满足两个条件：满足两个条件：（1）不破坏平衡状态；）不破坏平衡状态；（2）不破坏约束条件，即为）不破坏约束条件，即为约束所允许。约束所允许。qP应力边界应力边界 S位移边界位移边界 Su变化后的位移状态：变化后的位移状态：称为称为位移的变分位移的变分，或，或虚位移虚位移。（2）考察弹性体的能量变化）考察弹性体的能量变化：由能量守恒原理：由能量守恒原理：弹性体变

18、形势能的增加，等于外力势能的减少。弹性体变形势能的增加，等于外力势能的减少。（在没有温度改变、动能改变的情况下）（在没有温度改变、动能改变的情况下）设：设：表示弹性变形势能的增量；表示弹性变形势能的增量；表示外力在虚位移上所做的功，它在数值上等于外力势能的表示外力在虚位移上所做的功，它在数值上等于外力势能的减少。减少。则有：则有：第20页，讲稿共123张，创作于星期三外力的虚功：外力的虚功：体力：体力：面力：面力：外力外力代入前式：代入前式：（11-6）表明：表明：物体形变势能的变分，等于外力在虚位移上所做的虚功。物体形变势能的变分，等于外力在虚位移上所做的虚功。式（式（11-6）称为）称为位

19、移变分方程位移变分方程，也称，也称 Lagrange 变分方程变分方程。第21页，讲稿共123张，创作于星期三3.虚功方程虚功方程由式（由式（b）:两边求变分两边求变分:将将 U1 视为应变分量视为应变分量的函数的函数由格林公式由格林公式:（11-4）第22页，讲稿共123张，创作于星期三表示表示：实际应力在虚应变上所做的虚功实际应力在虚应变上所做的虚功 内力的虚功内力的虚功将上式代入位移变分方程（将上式代入位移变分方程（11-6），有），有（11-7）虚功方程虚功方程表明表明表明表明：如果如果在虚位移发生前，弹性体处于在虚位移发生前，弹性体处于平衡平衡状态，状态，则则在虚位移发生过程中，在虚

20、位移发生过程中，外力外力在虚位移上所做的虚功，等于应力在虚应变上所做的虚功。在虚位移上所做的虚功，等于应力在虚应变上所做的虚功。虚功方程虚功方程 是是有限单元法有限单元法的理论基础，也是许多的理论基础，也是许多变分原理变分原理的基础。的基础。第23页，讲稿共123张，创作于星期三4.最小势能原理最小势能原理 也是位移变分方程的一个应用也是位移变分方程的一个应用由位移变分方程：由位移变分方程：由于虚位移为由于虚位移为微小的微小的、为约束所允许为约束所允许的，所以，可认为在虚位移发生过程中，外力的，所以，可认为在虚位移发生过程中，外力的大小和方向都不变，只是作用点位置有微小变化。的大小和方向都不变

21、，只是作用点位置有微小变化。于是，有：于是，有：以以为零势能状态，为零势能状态，并用并用 V 表示任意状态的外力势能，则表示任意状态的外力势能，则外力在外力在可能位移可能位移上所做的功上所做的功W，即，即代入前式，有代入前式，有第24页，讲稿共123张，创作于星期三其中：其中：形变势能与外力势能的总和，形变势能与外力势能的总和，称为称为系统的总势能系统的总势能表明：表明：表明：表明：在给定的外力作用下，实际存在的位移应使系统的总势能的变分为零。在给定的外力作用下，实际存在的位移应使系统的总势能的变分为零。等价于总势能等价于总势能 U+V 取驻值。取驻值。极值势能原理极值势能原理平衡状态：平衡状

22、态：（1）稳定平衡状态；）稳定平衡状态；（2）不稳定平衡状态；）不稳定平衡状态；（3）随宜平衡状态；）随宜平衡状态；稳定平衡稳定平衡不稳定平衡不稳定平衡随宜平衡随宜平衡 势能取势能取极小值极小值 势能取势能取极大值极大值 不定不定最小势能原理最小势能原理最小势能原理最小势能原理：在给定的外力作用下，满足位移边界条件的在给定的外力作用下，满足位移边界条件的各组位移中，实际存在的位移，应使系统的总势各组位移中，实际存在的位移，应使系统的总势能成为驻值。当系统处于能成为驻值。当系统处于稳定平衡稳定平衡时，总势能取时，总势能取极小值，通常也为最小值。极小值，通常也为最小值。第25页，讲稿共123张，创

23、作于星期三实际存在的位实际存在的位移应满足：移应满足：（1）位移边界条件；）位移边界条件；（2）平衡方程（位移形式）；）平衡方程（位移形式）；（3）应力边界条件。）应力边界条件。（1）位移边界条件；）位移边界条件；（2）位移变分方程。）位移变分方程。因而，有：因而，有：位移变分方程位移变分方程（1）平衡方程；）平衡方程；（2）应力边界条件。）应力边界条件。（可互相导出）（可互相导出）（最小势能原理）（最小势能原理）5.伽辽金变分方程伽辽金变分方程 由虚功方程建立当位移分量满足：由虚功方程建立当位移分量满足：位移边界条件位移边界条件、应力边界条件应力边界条件时，弹性时，弹性体的体的位移变分应满足

24、的条件。位移变分应满足的条件。将将虚应变虚应变用虚位移表示：用虚位移表示：（c）将其代入虚功方程：将其代入虚功方程：第26页，讲稿共123张，创作于星期三（11-7）第27页，讲稿共123张，创作于星期三 同理，可得到其余各项的结果：同理，可得到其余各项的结果：将其代入虚功方程左边，有：将其代入虚功方程左边，有：（11-7）第28页，讲稿共123张，创作于星期三 将其代入虚功方程，并整理有：将其代入虚功方程，并整理有：当应力边界条件满足时，当应力边界条件满足时，000 上式可简化为：上式可简化为：（11-7）第29页，讲稿共123张，创作于星期三（10-8）伽辽金（伽辽金（Galerkin）变

25、分方程）变分方程表明：表明：表明：表明：当所取位移分量当所取位移分量同时满足同时满足：位移边界条件、应力边界条件时，：位移边界条件、应力边界条件时，其位移变分需满足的方程。其位移变分需满足的方程。第30页，讲稿共123张，创作于星期三（11-6）（1）位移变分方程）位移变分方程（2）虚功方程）虚功方程（11-7）位移变分方程小结：位移变分方程小结：也称也称 Lagrange 变分方程变分方程：（3）最小势能原理）最小势能原理说明：说明：说明：说明：（1）只要求：可能（虚）位移满足位移边界条件；）只要求：可能（虚）位移满足位移边界条件；（2）对虚功方程，也适用）对虚功方程，也适用各种材料的物理方

26、程各种材料的物理方程。如：塑性材料、非线性弹性材料等。如：塑性材料、非线性弹性材料等。第31页，讲稿共123张，创作于星期三（4）伽辽金（）伽辽金（Galerkin）变分方程）变分方程要求要求：可能（虚）位移满足：可能（虚）位移满足：（1）位移边界条件；）位移边界条件；（2）应力边界条件。）应力边界条件。（10-8）第32页，讲稿共123张，创作于星期三11-3 11-3 位移变分法位移变分法1.里兹里兹（Ritz）法法基本思想：基本思想：设定位移函数设定位移函数的表达形式，使其的表达形式，使其满足位移边界条件满足位移边界条件，其中含有若，其中含有若干待定常数，然后干待定常数，然后利用位移变分

27、方程确定这些常数利用位移变分方程确定这些常数，即得位移解。，即得位移解。设取位移的表达式如下：设取位移的表达式如下：（11-9）其中：其中：为互不相关的为互不相关的 3m 个系数；个系数；为设定的函数，且在边界上有：为设定的函数，且在边界上有：为为边界上为零边界上为零的设定函数的设定函数 显然，上述显然，上述函数满足边界函数满足边界条件条件。此时，位移的变分此时，位移的变分只能由系数只能由系数 Am、Bm、Cm的变分来实现。的变分来实现。与变分无关。与变分无关。第33页，讲稿共123张，创作于星期三（a）位移的变分：位移的变分：形变势能的变分：形变势能的变分：由式（由式（11-5），可知：），

28、可知：（b）将式（将式（a）、（）、（b）代入位移变分方程，有：）代入位移变分方程，有：（11-5）第34页，讲稿共123张，创作于星期三将上式整理、移项、合并，可得：将上式整理、移项、合并，可得：完全任意，且互相独立，完全任意，且互相独立，要使上式成立，则须有：要使上式成立，则须有：第35页，讲稿共123张，创作于星期三（11-10）Ritz 法方程法方程或称或称 Rayleigh-Ritz 法方程法方程说明：说明：（1）由由 U 的位移表达式（的位移表达式（11-5）可知，）可知，U 是系数是系数的二次函数，的二次函数，因而，方程（因而，方程（11-10）为各系数的）为各系数的线性方程线性

29、方程 组组。互不相关，因而，总可以求出全部的系数。互不相关，因而，总可以求出全部的系数。（2）求出了系数求出了系数就可求得其它量，如位移、应力等就可求得其它量，如位移、应力等（3）在假定位移函数时，须保证其在假定位移函数时，须保证其满足全部位移边界条件满足全部位移边界条件。第36页，讲稿共123张，创作于星期三2.伽辽金伽辽金（Galerkin）法法设取位移的表达式如下：设取位移的表达式如下：（11-9）同时满足：同时满足：（1）位移边界条件；）位移边界条件；（2）应力边界条件；）应力边界条件；位移的变分：位移的变分：将其代入伽辽金变分方程（将其代入伽辽金变分方程（10-8）：）：得到：得到：

30、第37页，讲稿共123张，创作于星期三完全任意，且互相独立，完全任意，且互相独立，要使上式成立，则须有：要使上式成立，则须有：第38页，讲稿共123张，创作于星期三将物理方程和几何方程代入，有将物理方程和几何方程代入，有（11-11）伽辽金（伽辽金（Galerkin）法方程）法方程说明：说明：（1）与与 Ritz 法类似，得法类似，得 3m 阶的线方程组，可求出阶的线方程组，可求出3m个系数。个系数。（2）伽辽金（伽辽金（Galerkin）法与）法与 Ritz 法的区别：在于设位移函数时，前者要求法的区别：在于设位移函数时，前者要求同时满足应力、位移边界条件同时满足应力、位移边界条件，而后者，

31、而后者只要求满足位移边界条件只要求满足位移边界条件。第39页，讲稿共123张，创作于星期三位移变分法的应用：位移变分法的应用：（1）求解弹性体的近似解；）求解弹性体的近似解；（2）推导弹性体的平衡微分方程与自然（力）边界条件；）推导弹性体的平衡微分方程与自然（力）边界条件；第40页，讲稿共123张，创作于星期三11-4 11-4 位移变分法应用于平面问题位移变分法应用于平面问题1.形变势能表达式形变势能表达式对于平面应变问题：对于平面应变问题：且且由式（由式（11-5）（11-12）对于平面应力问题：对于平面应力问题：（11-13）（11-5）第41页，讲稿共123张，创作于星期三2.位移函数

32、设定位移函数设定 由于，两种平面问题中，都不必考虑由于，两种平面问题中，都不必考虑 z 方向的位移方向的位移w，所以位移分量可设为：，所以位移分量可设为：（11-14）式中：各系数的含义和以前相同。式中：各系数的含义和以前相同。3.变分法方程变分法方程Ritz 法方程：法方程：（在（在 z 方向取单位长度）方向取单位长度）（11-15）Galerkin 法方程：法方程：第42页，讲稿共123张，创作于星期三Galerkin 法方程：法方程：（11-16）适用于适用于平面应变问题平面应变问题式中：式中：对于对于平面应力问题平面应力问题：（11-17）第43页，讲稿共123张，创作于星期三例：例：

33、图示薄板，宽为图示薄板，宽为 a，高度为，高度为 b，左边和下边受连杆支，左边和下边受连杆支承，右边和上边分别受有均布压力承，右边和上边分别受有均布压力 q1和和 q2 作用，不计体作用，不计体力。试求薄板的位移。力。试求薄板的位移。解：解：（1）假设位移函数）假设位移函数（a）满足边界条件：满足边界条件：试在式（试在式（a）中只取两个系数：）中只取两个系数：A1、B1，即，即（b）（2）计算形变势能）计算形变势能 U将式（将式（b）代入（）代入（11-13），有），有（平面应力情形下形变势能公式）（平面应力情形下形变势能公式）积分得：积分得：（c）第44页，讲稿共123张，创作于星期三（c）

34、（3）代入）代入Ritz 法方程求解法方程求解体力体力有有在右边界：在右边界：在上边界：在上边界：于是有：于是有：将式（将式（c）代入，得）代入，得（11-15）第45页，讲稿共123张，创作于星期三联立求解，得：联立求解，得：（f）代入位移表达式（代入位移表达式（b），得：），得：（g）讨论：讨论：（1）如果在位移式（如果在位移式（a）中再多取一此系数）中再多取一此系数如：如：A2、B2等，但是经计算，这些系等，但是经计算，这些系数全为零。数全为零。（2）位移解（位移解（g）满足几何方程、平衡方程）满足几何方程、平衡方程和边界条件。和边界条件。表明：表明：位移解（位移解（g）为问题的）为问题

35、的精确解精确解。第46页，讲稿共123张，创作于星期三例：例：图示矩形薄板，宽为图示矩形薄板，宽为2 a，高度为，高度为2 b，左右两边和，左右两边和下边均被固定，而上边的给定位移为：下边均被固定，而上边的给定位移为：（h）不计体力。试求薄板的位移和应力。不计体力。试求薄板的位移和应力。解：解：（1）假设位移函数）假设位移函数取取 m=1,将位移分量设为：将位移分量设为：（i）显然，可满足位移边界条件：显然，可满足位移边界条件：第47页，讲稿共123张，创作于星期三（2）代入）代入Galerkin 法方程求解法方程求解该问题中，无应力边界条件，式（该问题中，无应力边界条件，式（i）满足全部条件

36、。）满足全部条件。可用伽辽金（可用伽辽金（Galerkin）法求解。）法求解。X=Y=0，m=1，伽辽金法方程变为：伽辽金法方程变为：（11-17）（j）第48页，讲稿共123张，创作于星期三 将其代入伽辽金方程（将其代入伽辽金方程（j）,可求得：可求得：代回位移式（代回位移式（h）,有：有：第49页，讲稿共123张，创作于星期三代回位移表达式（代回位移表达式（h）,得位移解答：得位移解答：当当 b=a，取，取 =0.2时，上述解答成为时，上述解答成为：第50页，讲稿共123张，创作于星期三（3）求应力分量）求应力分量应用几何方程及物理方程，可求得应力为：应用几何方程及物理方程，可求得应力为：

37、第51页，讲稿共123张，创作于星期三在在处，处，相应的面力：相应的面力：第52页，讲稿共123张，创作于星期三例：例：如图所示简支梁，中点处承受有集中如图所示简支梁，中点处承受有集中P，试，试求的梁的挠曲线方程。求的梁的挠曲线方程。PABlxy解：解：（1）假设位移试探函数）假设位移试探函数（必须满足位移边界条件）（必须满足位移边界条件）设位移试探函数为（取一项）：设位移试探函数为（取一项）：式中：式中：a 为待定常数。为待定常数。（2）计算：）计算：（a）（b）显然，式（显然，式（a）满足端点的位移边界条件）满足端点的位移边界条件:（3）代入）代入Ritz 法方程，求解法方程，求解（c）（

38、d）第53页，讲稿共123张，创作于星期三PABlxy讨论：讨论：（1）中点的挠度：中点的挠度：（e）而材料力学的结果：而材料力学的结果：两者比较：两者比较：式（式（a）的结果偏小）的结果偏小2%。如果取如下位移函数：如果取如下位移函数：式中项数式中项数 m 取得越多，则求得精度就越高。取得越多，则求得精度就越高。（2）所取的位移函数所取的位移函数必须满足位移边界条件必须满足位移边界条件。（3）位移函数选取不是唯一的，如：位移函数选取不是唯一的，如：第54页，讲稿共123张，创作于星期三PABlxy例例：如图所示简支梁，中点处承受有集中如图所示简支梁，中点处承受有集中P，试求的梁的挠曲线方程。

39、，试求的梁的挠曲线方程。解解：（1）假设位移试探函数）假设位移试探函数式中：式中：A1、A2 为待定常数。为待定常数。显然，式（显然，式（a）满足端点的位移边界条件）满足端点的位移边界条件:（2）计算：）计算：梁的形变势能梁的形变势能:（3）代入）代入 Ritz 法方程法方程:第55页，讲稿共123张，创作于星期三PABlxy例例：如图所示简支梁，中点处承受有集中如图所示简支梁，中点处承受有集中P，试求的梁的挠曲线方程。，试求的梁的挠曲线方程。解解：位移函数位移函数（a）（3）代入）代入 Ritz 法方程法方程:所求挠曲线方程所求挠曲线方程:第56页，讲稿共123张，创作于星期三PABlxy所

40、求挠曲线方程所求挠曲线方程:中点挠度中点挠度:而材料力学的结果：而材料力学的结果：说明：说明：（1）设定的待定系数个数与所得的线性方程数相同；）设定的待定系数个数与所得的线性方程数相同；（2）亦可用亦可用最小势能原理最小势能原理求解上述问题。求解上述问题。第57页，讲稿共123张，创作于星期三例：例：如图所示简支梁，中点处承受有集中如图所示简支梁，中点处承受有集中P，试，试求的梁的挠曲线方程。求的梁的挠曲线方程。PABlxy解：解：（1）假设位移试探函数）假设位移试探函数（必须满足位移边界条件）（必须满足位移边界条件）设位移试探函数为：设位移试探函数为：式中：式中：a 为待定常数。为待定常数。

41、（2）求系统的总势能）求系统的总势能（a）（b）（c）将式（将式（a）代入，计算得）代入，计算得显然，式（显然，式（a）满足端点的位移边界条件）满足端点的位移边界条件:（3）由最小势能原理确定常数）由最小势能原理确定常数（d）第58页，讲稿共123张，创作于星期三PABlxy说明：说明：（1）（e）与与 Ritz 法结果相同；法结果相同；（2）所取的位移函数所取的位移函数必须满足位移边界条件必须满足位移边界条件。（3）位移函数选取不是唯一的，如：位移函数选取不是唯一的，如：第59页，讲稿共123张，创作于星期三例：例：如图所示，一端固定，另一端有弹性支承的梁，如图所示，一端固定，另一端有弹性支

42、承的梁，跨度为跨度为 l，抗弯刚度为，抗弯刚度为EI，弹簧的刚度为，弹簧的刚度为 k，梁上，梁上作用有分布载荷作用有分布载荷q(x)，试用最小势能原理导出梁，试用最小势能原理导出梁的弯曲微分方程和边界条件。的弯曲微分方程和边界条件。解：解：（1）求系统的总势能）求系统的总势能系统的总势能系统的总势能=梁的弯曲变形能梁的弯曲变形能+弹簧的变形能弹簧的变形能+外力势能外力势能(a)式中：式中：w 为梁的挠度。为梁的挠度。由最小势能原理：由最小势能原理：(b)分部积分：分部积分：（2）对总势能求变分）对总势能求变分第60页，讲稿共123张，创作于星期三将其代入式（将其代入式（b），有），有第61页，

43、讲稿共123张，创作于星期三梁的左端固定，梁的左端固定，有有代入上式，有：代入上式，有：的任意性与相互独立性，的任意性与相互独立性，有有（3）利用位移边界条件和变分的任意性确定所需的结）利用位移边界条件和变分的任意性确定所需的结果。果。第62页，讲稿共123张，创作于星期三 弯曲微分方程弯曲微分方程 力的边界条件力的边界条件表明：表明：最小势能原理最小势能原理等价于等价于平衡微分方程和力的边界条件；平衡微分方程和力的边界条件；第63页，讲稿共123张，创作于星期三Ritz 法解题步骤：法解题步骤：（1）假设位移函数，使其）假设位移函数，使其满足边界条件满足边界条件；（2）计算形变势能计算形变势

44、能 U；（3）代入）代入Ritz 法方程求解待定系数；法方程求解待定系数；（4）回代求解位移、应力等。）回代求解位移、应力等。用最小势能原理解题步骤：用最小势能原理解题步骤：（1）假设位移函数，使其）假设位移函数，使其满足边界条件满足边界条件；（2）计算系统的总势能计算系统的总势能 ；（3）由最小势能原理：由最小势能原理：=0，确定待定系数；，确定待定系数；（4）回代求解位移、应力等。）回代求解位移、应力等。第64页，讲稿共123张，创作于星期三图示简支梁，两端受轴向压力图示简支梁，两端受轴向压力P 作用，在距左端距离作用，在距左端距离 c处受集中力偶处受集中力偶 M 作用，梁的作用，梁的跨度

45、为跨度为 l。试用最小势能原理求的梁的挠曲线方程。试用最小势能原理求的梁的挠曲线方程。例：例：设梁的挠曲线方程可设为：设梁的挠曲线方程可设为：解：解：设定梁的挠曲线函数求系统的总势：设定梁的挠曲线函数求系统的总势：第65页，讲稿共123张，创作于星期三代入总势能计算公式：代入总势能计算公式：由最小势能原理求出待定系数由最小势能原理求出待定系数：由于，由于，Am不能等于零，可求得：不能等于零，可求得：第66页，讲稿共123张，创作于星期三梁的挠曲线方程为：梁的挠曲线方程为：梁的最小失稳载荷为：梁的最小失稳载荷为：第67页，讲稿共123张，创作于星期三11-5 11-5 应力变分方程应力变分方程1

46、.形变余能形变余能Pxl0lO（1）单向应力状态）单向应力状态设：设：一般的应力应变关系一般的应力应变关系形变势能（比能）：形变势能（比能）：dd00 单位体积的形变势能（比能）单位体积的形变势能（比能）形变余能（比能）：形变余能（比能）：单位体积的形变余能单位体积的形变余能（比余能）（比余能）对线弹性体，显然有：对线弹性体，显然有：形变形变势能势能（比能）等于形（比能）等于形变变余能余能（比余能）（比余能）表明：形变表明：形变比余能比余能在数值上等于图中矩形面积减去在数值上等于图中矩形面积减去 U1 后余下的面积。后余下的面积。第68页，讲稿共123张，创作于星期三（2）三向应力状态）三向应

47、力状态对线弹性体，有：对线弹性体，有：物体形变余能：物体形变余能：对线弹性体对线弹性体：物体形变余能常用应力表示：物体形变余能常用应力表示：第69页，讲稿共123张，创作于星期三（3）形变余能的变分）形变余能的变分对照形变余比能的表达式，有：对照形变余比能的表达式，有：由应力表示的卡氏（由应力表示的卡氏（Castigliano）定理）定理代入形变余比能的变分表达式，有：代入形变余比能的变分表达式，有：第70页，讲稿共123张，创作于星期三若将上式中应变分量利用几何方程表示成位移形式，有：若将上式中应变分量利用几何方程表示成位移形式，有：代入形变余比能的变分表达式，有：代入形变余比能的变分表达式

48、，有：第71页，讲稿共123张，创作于星期三2.应力变分方程应力变分方程设有任一弹性体，在外力的作用下处于平衡状态。其应力和位移分别为：设有任一弹性体，在外力的作用下处于平衡状态。其应力和位移分别为：实际的应力和位移实际的应力和位移建立：物体建立：物体形变余能的变化形变余能的变化与与应力变分应力变分之间的之间的关系关系。（1）应力的变分）应力的变分假设：作用于物体的体力不变，而应力分量发生如下变分：假设：作用于物体的体力不变，而应力分量发生如下变分：常称为常称为虚应力虚应力满足：满足：（1）平衡微分方程；）平衡微分方程；（2）应力边界条件）应力边界条件（即：在应力边界上变分应为零）。（即：在应

49、力边界上变分应为零）。变化后应力状态：变化后应力状态：（2）应力变分方程）应力变分方程都满足平衡方程都满足平衡方程并作用于同样的体力，并作用于同样的体力，将其分别代入平衡微分方程，并进行比较，应有：将其分别代入平衡微分方程，并进行比较，应有：第72页，讲稿共123张，创作于星期三（a）张量表示张量表示在在位移给定位移给定的边界上，的边界上，由于应力的变分必然引起该边界上面力的变分：由于应力的变分必然引起该边界上面力的变分：由边界上应力与边界面间关系，在位移给定边界上，应有：由边界上应力与边界面间关系，在位移给定边界上，应有：（b）张量表示张量表示第73页，讲稿共123张，创作于星期三由形变余能

50、的变分：由形变余能的变分：利用奥利用奥-高公式，将上式每一项作变换，如：高公式，将上式每一项作变换，如：将其代入应变余能的变分，并整理有：将其代入应变余能的变分，并整理有：第74页，讲稿共123张，创作于星期三000得到：得到：（11-18）上式表明：上式表明：由于应力的变分，由于应力的变分，形变余能的变分形变余能的变分等于等于面力的变分在实际位移上所做的功面力的变分在实际位移上所做的功（虚功）。（虚功）。应力变分方程，应力变分方程，也称也称Castigliano变分方程。变分方程。第75页，讲稿共123张，创作于星期三说明：说明：（1）要求应力的变分满足：要求应力的变分满足：平衡微分方程；平