单元形函数的讨论解析.ppt

4.6 根据形函数基本原理进行三维实体分析 按以上分析，可以求得4-12节点的位移列阵 即 结果为 1.0e-006 * -0.1160 -0.0204 -0.1927 0.1780 0.0084 0.8733 0.1182 -0.0199 0.1287 -0.1627 0.0124 0.0470 -0.0816 -0.0206 0.0966 0.1007 -0.0179 0.1426 0.0849 0.0371 0.1780 -0.0984 -0.0470 0.0112 4.6 根据形函数基本原理进行三维实体分析 （7）求各节点外力、应力和应变 求出了各节点位移，紧接着可以求出各节点力。利用几何方程求出应变，利用物理方程求出应立，从而完成静力学分析。 4.2 形函数的性质 类似地有 (4.34) 第二，在单元的任一节点上，三个形函数之和等于1，即 (4.34) 简记为 (4.36) 这说明，三个形函数中只有二个是独立的 4.2 形函数的性质 第三，三角形单元任意一条边上的形函数，仅与该边的两端节点坐标有关、而与其它节点坐标 无关。 例如，在2-3 边上有 (4.36) 根据形函数的这一性质可以证明，相邻单元的位移分别进行线性插值之后，在其公共边上将是连续的。例如，单元1-2-3和1-2-4具有公共边1-2。由上式可知，在1-2边上两个单元的第三个形函数都等于0，即 (4.37) 4.2 形函数的性质 不论按哪个单元来计算，公共边1-2上的位移均由下式表示 (4.38) 可见，在公共边上的位移u、v 将完全由公共边上的两个节点1、2的位移所确定，因而相邻单元的位移是保持连续 的 为了能够更好地理解形函数的概念，这里引入面积坐标。在如图4-13所示的三角形单元ijm中，任意一点P(x , y)的位置可以用以下三个比值来确定。 4.3 用面积坐标表示的形函数 图4-13平面三角形单元的面积坐标 (4.40) 式中，?——三角形单元ijm的面积，?i 、?j 、?m 为三角形Pjm、Pmi、Pij的面积。Li，Lj，Lm叫做P点的面积坐标。显然，这三个面积坐标不是完全独立的，这是由于 4.3 用面积坐标表示的形函数 ?i +?j +?m =? (4.41) 所以有 Li +Lj +Lm =1 (4.42) 对于三角形Pjm，其面积为 (4.43) 故有 (4.44) 类似地有 (4.44) (4.46) 4.3 用面积坐标表示的形函数 可见，前面讲述的平面三角形单元的形函数Ni 、Nj 、Nm 等于面积坐标Li 、Lj 、Lm 。 容易看出，单元三个节点的面积坐标分别为 节点 i： Li =1 Lj =0 Lm =0 节点 j： Li =0 Lj =1 Lm =0 节点m： Li =0 Lj =0 Lm =1 根据面积坐标的定义，平行于jm边的某一直线上的所有各点都有相同的坐标Li，并且等于该直线至jm边的距离与节点i至jm边的距离之比，图4-13中给出了Li的一些等值线。平行于其它边的直线也有类似的情况。 4.3 用面积坐标表示的形函数 (4.47) 当面积坐标的函数对直角坐标求导时，有下列公式 (4.48) 不难验证，面积坐标与直角坐标之间还存在以下变换关系： 4.3 用面积坐标表示的形函数 求面积坐标的幂函数在三角形单元上的积分时，有 (4.49) 式中, ?、?、? 为整常数。求面积坐标的幂函数在三角形某一边上的积分值时，有 (4.40) 式中, l为该边的长度。 4.4 有限元的收敛准则 对于一个数值计算方法，一般总希望随着网格的逐步细分所得到的解答能够收敛于问题的精确解。根据前面的分析，在有限元中，一旦确定了单元的形状，位移模式的选择将是非常关键的。由于载荷的移置、应力矩阵和刚度矩阵的建立都依赖于单元的位移模式，所以，如果所选择的位移模式与真实的位移分布有很大的差别，会将很难获得良好的数值解。 可以证明，对于一个给定的位移模式，其刚度系数的数值比精确值要大。所以，在给定的载荷之下，有限元计算模型的变形将比实际结构的变形小。因此细分单元网格，位移近似解将由下方收敛于精确解，即得到真实解的下界。 4.4 有限元的收敛准则 ⑴ 位移模式必须包含单元的刚体位移。也就是，当节点位移由某个刚体位移引起时，弹性体内将不会产生应变。所以，位移模式不但要具有描述单元本身形变的能