第二章 有限元法的基本概念.pptVIP

平面问题的应变协调方程 第二章 有限元法的基本概念 2.1 有限元法的产生及应用 有限元法(Finite Element Method，简称FEM)，又称 有限单元法、有限元素法，它是随着大型高速数字电子计算机的出现而发展起来的一种有效的数值计算方法。 传统的材料力学给出了某些计算位移和应力的解析公 式。但它只能用于几种简单形体，而且在推导那些公式时 用到一些基本假设。例如，等截面杆、简支梁、悬臂梁、 轴等。 有些用材料力学不能计算的问题如压力集中、深梁、 板弯曲、厚壁圆筒等，可以用弹性力学理论来解释。如对 于平板问题来说，矩形板、圆形板都能推导出计算公式， 不规则形状的板就不能。 2.1 有限元法的产生及应用 有限元法用离散化(discretization)的方法彻底地解决了 这一问题。不管多么复杂的结构，总可以用一个离散化的模 型来代替。整个结构被划分成有限个小的区域，称为单元， 所划分的每个单元的几何形状是有一定规则的。先用弹性力 学理论计算出每个单元的特性，然后再把各单元集合在一起 形成完整的结构。用适当的方法对这个离散化的模型求解， 就可以代替对原有结构的求解。 2.2 有限元法采用的变分原理 在固体力学中，有限元法的应用都是根据变分原理来推 导单元特性和有限元方程的。最常用的三种变分原理是：最 小势能原理、余能原理和雷纳斯原理。采用不同的变分原理 会得到不同的未知场变量。 当采用势能原理时，必须假设单元内位移场函数的形 式，这种有限元分析方法称为位移法或协调法； 当采用余能原理时，必须假设单元内应力场的形式，这 种方法称为力法或平衡法； 当采用雷纳斯原理时，就必须同时假设某些位移和某些 应力，这种方法称为混合法。 2.2 有限元法采用的变分原理 关于变分概念 微分是变量的增量，变分是函数的增量，通 常用δ表示。 2.3 单元与离散化模型 单元是有限元法的基本环节，任何结构都被离散化成单 元，所谓离散化，就是将弹性连续体用线或面分成许多互不 重迭的有限多个有限大小的构件，这些有限大小的构件就称 为有限单元，也简称单元。单元与单元之间在一些点上相连 接，称为“结点”或“节点”(node，joint)。 单元可分为一维单元(one-dimensional element)、二维单 元(two-dimensional element)和三维单元(three-dimensional element)。我们分别简单介绍一下各种单元。 2.3.1 一维单元 一维单元主要有杆单元和梁单元，杆单元与梁单元都是断面尺寸远小于长度尺寸的构件，也就是说它的基本形状是细长形的。有限元法中认为，相邻的单元在端部结点处互相连接。对于杆单元是互相铰接的，即假设杆只受轴向力，两杆的夹角可以改变，杆与杆之间只传递力而不能传递力矩。对于梁单元是刚接的(或称固接)，即假设梁与梁之间既能传递力又能传递力矩。杆单元的变形只是伸长与缩短，不可能发生弯曲而梁单元除了长度方向上变形之外，还可能产生剪切、扭转和弯曲变形。 在有限元模型网格图中，只表现它的长度而无另外两维几何尺寸。它的断面参数将和材料、温度等参数一起作为单元的性质另外提供。 例子： 2.3.3 三维单元 不能简化为二维问题的三维连续体也需要人为地确定单元划分的分 界面，如图所示的轴承座结构，可离散化为由三维单元所组成的模型。 三维单元的形状是个三维实体，所以不需要任何其它数据来描述几何形 状。单元与单元在结点处作空间铰接(球面铰接)，单元间只传递任意方向 的力而不能传递力矩。三维单元通常划分为六面体。 2.4 刚度矩阵和柔度矩阵 弹簧的刚度和柔度 一个弹簧受有力P作用，使之伸长为u，在弹簧极限范围内，拉力 愈大则延伸愈长，其关系如下： c表示施加在弹簧上的单位力所引起的弹簧变形，称为柔度。 k是使弹簧产生单位位移(变形)需要加在弹簧上的内力，称为刚度。 设有一弹性体，在其上作用有广义力P1、P2、…、Pi、…、Pn，作 用点分别为1、2、…、i、…、n。设在支座的约束下，弹性体不能发生 刚体运动，仅产生弹性变形。其相应的广义位移为δ1、δ2、…、 δi、…、δn，如以节点i的位移δ为例，广义位移是弹性体受这一组广义力P1、P2、…、Pi、…、Pn共同影响的结果。 如果用位移来表示载荷，则方程组有如下的形式 矩阵[k]称为刚度矩阵。而每个系数kij表示节点j有单位位移(δj=1)而在节点