第9章-能量法.pptVIP

第9章 能量原理及变分法 求解弹性力学问题，可以归结为求解偏微分方程的边值问题，要在给定的边界上，求得微分方程的精确解，在数学上往往比较困难，有时甚至是不可能的。为此，必须寻求能提供足够精确的近似解法，其中基于能量原理的变分法以及以变分法为基础的有限单元法，才使得弹性力学在工程中获得广泛应用成为可能。 * * 能量原理的各个判据跟能量泛函的极值相联系。泛函就是函数的函数。如应变能密度W是σij或εij的函数，而 又是x,y,z的函数，所以W称为泛函。求泛函极值的问题，称为变分问题。因此，用能量原理的一些定则或判据求解弹性力学问题的方法，称为变分法。 §9.1 虚位移原理 则在体积内有 图9-1 (给定位移Su) (给定面力Sσ) (9.1) 在边界上 (9.2) 现设想一个处于平衡状态的物体，由于某种原因发生了约束所允许的任意微小的虚位移 。实际的力系在虚位移上所作的功称为虚功。 即弹性体的外力在约束所允许的虚位移上所作的功为， (9.3a) 或 (9.3b) 在物体产生微小虚变形的过程中，该物体的总虚应变能为 或 (9.4b) 于是虚位移原理表达为 (9.5a) 或 (9.5b) 或 虚位移原理表述为：外力在约束所允许的虚位移上所作的虚功等于相应应力在虚应变上所产生的总虚应变能。式(9-5)又称为拉格朗日变分方程。 现在给出式(9-5)的证明： (9.6) 由奥-高公式： (9.7) 将式(9-7)代入式(9-6)得： (9.8) 由平衡条件： 和几何方程， (9.9) 即 (9.10) 于是得到： 故 式(9-5)得证。 然后再来证明虚位移原理等价于平衡条件和应力边界条件。 注意式(9.10)，则(9-5b)式左端为： (9.11) 根据高斯的分部积分公式： (9.12) 将式(9.12)代入式(9.11)，再将式(9.11)代入(9-5b)式，即为 移项整理得： (9.13) 由于 的任意性，有物体的平衡方程： 和应力边界条件： 可见虚位移原理等价于平衡条件和应力边界条件. 另外，应当指出，式(9.5b)等号左边表示由于产生虚位移 ，而引起的物体内的总虚应变能。这种虚位移实际上应理解为真实位移的变分，而不是其它随便的一种位移函数，而且式中(9-5c)的 ，也不是别的什么虚应变，而是由于 引起的。 §9.2 最小势能原理 当存在应变能密度函数W时，虚功方程可写为 将上式中的变分符号移至积分号以外。 于是有， (9.15) 其中： (9.16) 附加条件为： （在 上） 可以证明，总势能的二阶变分恒为正，即 ，于是最小势能原理表述为：弹性体在所有满足给定几何边界条件的位移场中，真实的位移场使物体的总势能取最小值。 最小势能原理说明：真实的位移场除了满足几何边界条件外，还必须满足最小势能原理的变分方程。而它又完全等价于平衡方程和应力边界条件。二者之间只是形式上的不同，这种形式上的变更扩大了弹性力学的解题的范围。另外，虚功方程和最小势能原理是等价的，它们都是弹性体在实际平衡状态的位移发生虚位移时，能量原理的具体应用，只是表达方式不同而已。 §9.3 位移变分法 其基本思路是：设定满足位移边界条件的含有待定系数的位移分量的表达式，然后由位移变分方程(9-5)或(9-15)得出含有待定系数的线性方程组，解方程得系数，从而得到位移分量的具体表达式。 选取几何允许的位移函数为 (9.17) 其中 为三个待定系数，位移函数应满足下列要求： (1). 满足位移边界条件，即， （在 上） (2). 为坐标x,y,z的线性独立的设定函数，且满足 （在 上） 对位移取变分，由于设定的函数 仅为坐标的函数，不参加变分，所以位移的变分是由对 系数的变分实现的，于是得， (9.18) 应变能的变分为 (9.19) 将式(9.18)、(9.19)代入(9.5)式或(9.15)式得： (9.20) 由于变分 的任意性，可得， (9.21) 由式(4.56)及式(9.5)可见，应变能应是系数Am、Bm、Cm的二次函数，因而方程(9.6)是这些系数的一次方程，方程数目为3n个，与未知数数目相等，由式(9-14)可以确定全部系数，从而求得位移分量。设个方法称为李兹(Rit)法。 方程(9.21)也可以采用另外的形式，总势能Π可以为位移分量的函