有限元分析及数学求解原理.ppt

对于该问题的最小势能原理，其数学变分问题为： 设有满足位移边界条件的许可位移场函数，其中真实的一组使得 以下泛函取极小值，即 计算剪支梁势能的变分（复合函数的求导）有 考虑到剪力、弯矩和位移之间的关系，以及边界条件等，势能的变分 可以变为 可以看出，要使得势能的变分等于0，必须要满足以下几个条件： 1、两端的弯矩等于0 2、力平衡方程 对于能量泛函求二次变分有 可以知道，能量泛函的极值为极小值。 其中又有 由变分方法，对泛函取极值，故有 故有 以上的证明说明，总能够找到满足位移边界条件的试函数，即许可 位移场，在满足几何条件和物理方程的条件下，当势能取最小值时， 其结果精度可以满足剩下的平衡方程以及力边界条件。 但是在实际上，由于选择许可位移场的局限性和盲目性，一般很难 将真正精确的位移场包含在许可位移场中，这样，就不可能由最小势 能原理来求出精确解，只能在所选择的试函数（即许可位移场）范围 内，通过最小势能原理求出最好的一组解，它在加权残值最小的意义 下，对平衡方程和力边界条件进行最佳逼近。 由以上的情况可以看出，弹性问题的最小势能原理的特点为： 1、试函数为许可位移场，即只需要满足位移边界条件，而不用满足力边界条件 2、积分中试函数的最高阶导数教低（对梁的弯曲问题导数为2阶，对于一般问题导数为1阶） 3、整个计算方法是计算以个全场的积分 4、由取积分问题的最小值，可以将原方程的求解化为线性方程组的求解 5、整个方法的处理流程比较规范 几种方法之间的关系 1、变分方法是从纯数学的角度来描述，只是用方程、泛函、极值等术语来表达和推导，而没有考虑背后的物理意义 2、最小势能原理和虚功原理的提法为物理提法，他们有相应的物理表达和物理意义。 3、加权残值发是一种数学提法，它基于试函数来求取对原始方程的最佳逼近。 在某种意义上，它们可以互相转化： 变分法的泛函极值——最小势能原理 变分过程——加权残值意义下对平衡方程和力边界条件的逼近。实际上都是变分方法的几种变换形式。 4 各种求解方法的特点及比较 从以上图表可以看出，微分形式的求解方法和积分形式的求解方 法存在着本质上的不同，正式引入了试函数，使得求解的难度大大降 低。 由于工程问题非常复杂，要求所采用的方法具有教好的规范性， 较低的难度，较低的函数连续性要求，较明确的物理概念，较好的通 用性。 而基于最小势能原理的求解方法具有比较明显的综合优势，因 此，可以在该原理的基础上发展出能广泛适用于工程中任意复杂问题 的求解方法，但是必须处理的关键问题包括：复杂物体的几何描述、试函数的确定与选取、全场试函数的表达。 基本方程——偏微分方程的边值问题 变分原理 ——偏微分方程的边值问题转换为代数方程 有限元原理——数值分析方法 工程应用广泛，理论基础——变分原理。 为什么变分原理在工程上的应用有限，而有限元原理却应用广泛。 有限元原理与经典变分原理的差别 位移试函数 位移边界Su 位移试函数构造困难 有限元不是整体选取试函数 而是在弹性体内分区（单元）完成的 试函数形式简单统一 近年来，随着现代科学技术的发展，特别是计算机技术的迅速发展和广泛应用，使得以有限元方法为代表的计算力学的迅速发展，改变了弹性力学理论在工程应用领域的处境。 有限元方法将计算数学与工程分析相结合，极大地扩展和延伸了力学理论与方法，取得了当代力学理论应用的高度成就。 大型通用有限元程序的广泛应用，使得有限元成为结构分析工具。 以此为基础，CAD, CAE等技术的应用使得计算机不仅成为数值分析的工具，而且成为设计分析的工具。 x方向平衡 y方向平衡 。 求解方程 得 其中c1,c2,c3是待定系数。最后可得 讨论 应变能 外力功 势能 由于单元有四个位移分量，可设梁单元的位移模式v(x)为包含4个待定常数的三次多项式： 根据边界条件可以确定待定系数，将其进一步回代，可以得到用节点位移表示的梁单元位移。 式中 根据梁的平面假定可知梁单元的轴向应变为： 这里利用平面假设（变形后横截面仍保持平面，与纵线正交）如图： 从而可以由单向虎克定律得出单元的轴向应力： 由虚功原理可以推得 2 弹性力学问题近似求解的加权残值法 直接针对原始三大类方程边界条件下求解三大类变量往往是非常 困难的，尤其是当几何形状和边界条件比较复杂时，一般求不出相应 的解析解。 如果事先假定满足一定边界条件的试函数，再在此基础上进行近 似求解，则可以大大降低求解难度。这种试函数方法可以使