第三章势能极小.ppt

第三章 势能极小原理的有限元解法 有限元法的出现是从实际结构的直接离散化开始的，但很快被人们认识到，它是变分问题里兹解法的一种。由此而巩固了有限元法的数学基础。 了解有限元法的数学基础 有助于理解方法的数学实质； 有助于扩大可用的函数形式； 有助于扩大方法的应用范围； 有助于估计近似计算得误差。 概念 当弹性体受到外力作用时，它就发生变形，因而外力在变形位移方向上对弹性体做功。 如果不计及弹性体在加载和卸载时能量的损失，即当结构系统是保守系统，则对这样的物体在变形时所做的功，可以看成是储存在物体中的能量，称为应变能。因此，应变能可以看成是弹性体变形时，它所吸收的能量。 总位能 在讨论弹性系统的能量时，不仅要考虑外力所做的功（对应于系统位能的，称为外力位能），还要考虑和变形相对应的应力（内力）所做的功，即应变能（对于系统位能的，有时称为内力位能）。 这样，当系统能量以位能表述时，则总位能包括外力位能和应变能。 应变能 应变能只和最终的应力状态有关，和应力作用的先后次序无关。 作用在物体上的外力可分为两类 表面力—分布在物体表面上的力，如加在梁上的载荷、流体静压力等都是表面力。 体积力—分布在物体体积中的力，如重力、离心力、惯性力等都是体积力。 如果结构系统仅有体积力作用，则外力位能为（平面问题） 如果结构系统仅有表面力作用，则外力位能为（平面问题） 当体积力和表面力同时作用时，则外力位能为 总位能 按照力学的一般说法，任何一个实际状态的弹性结构系统的总位能是这个系统从实际状态运动到某一参考状态（通常取结构的卸载状态作为参考状态）时它的所有作用力所做的功。 结构系统的作用力包括外力（体积力、表面力等）和内力（与结构在外力作下产生的变形相应的应力）。内力位能就是已阐明的应变能U。因为在卸载时，应力总是和应变同方向的，所以应变能是正值。外力位能是结构从它的最终位置（实际状态或受载状态）恢复到它的初始状态（卸载状态或参考状态）时，结构上的每一个外力所做的功W，显然它是负值。 这样，弹性结构在外力作用下的总位能可以写成 上式说明总位能п是一个函数的函数，即泛函。因为п是u，v的函数，而u，v又是x，y的函数。 上式表述了弹性体变形时的能量关系。对它应用变分法可以导出能量原理中的一些重要原理。例如，有限元法中很有用的最小位能原理等。 泛函 设有函数族y(x)，它的每一个函数值都有一个п值与之对应，则变量п叫做依赖于自变函数y(x)的泛函，记为п = п [y(x)] 。由于y(x)是自变量x的函数，所以泛函是函数的函数。 变分法 变分问题是研究泛函的极大值或极小值（简称极值）的问题。 它的解法非常类似于求解函数的极大值和极小值的方法。 即变分在泛函研究中所起的作用和微分在函数研究中所起的作用一样。 例如 若可微函数f(x)在内点x=x0达到极值，则在这点处有df=0； 而对变分问题来说，若具有变分的泛函п [y(x)] 在y= y0(x)上达到极值，则在y= y0(x)上有δп =0。 步骤 首先建立所讨论问题的泛函式п [y(x)] 继而进行δп的变分运算 自极值条件δп =0获得所论问题的微分方程 然后对该微分方程边值问题求解，得到函数y(x) ，它就是所论问题的正确解。 泛函变分的近似计算法 变分问题的微分方程在许多情况下不便于进行积分求解。 因而就产生“泛函变分的近似计算法”或称“变分问题中的直接法”。里兹法就是其中的一种。 里兹法 里兹法是就所讨论问题的整个区域来设定可能解的线性组合式的，因而有时比较困难、复杂。 和里兹法一样，有限元法也是求解泛函极值问题一种近似方法。它的方法步骤也和里兹法相似。 不同的是，有限元法是把所讨论问题的整个区域剖分成许多子区域，即通常所说的单元，然后在各个单元内采用和里兹法相似的办法求解。 最小位能原理 变形弹性体的总位能п，其中的位移u，v是满足连续条件的任意单值连续函数，它们将是泛函п的容许函数。 就是说， u，v可以具有很多组不同的值。 当然，不同组的u，v可以得出不同的泛函п值。 但是其中总会有一组u，v的值可使泛函п取最小值。 根据变分原理，使泛函取极值的那一组u，v的值就是问题的正确解。它反映在外力作用下，弹性体的真实变形。 物理意义 若有梁受外力{F}作用而变形，可以设想在满足变形连续条件的前提下，它有多种变形曲线的可能，如图中的虚线和实线所示。但是总有一条变形曲线是真正的变形曲线。例如，图中实线所示的变形曲线。其它的满足变形连续条件的可能变形曲线则是它的接近曲线 如果从变分原理来考虑，这条真实的变形曲线是从泛函п取极值，即从δп =0的条件计算出来的。 若从能量原量来考虑，则变形弹性体在受外力作用而处于平衡状态时，在很多可能的变形曲线（容许曲线）中，使总位能最小的那一条曲线是真正的变