第三章 平面问题有限单元法.docVIP

第三章 平面问题有限单元法 3.1平面问题的单元划分 1. 问题的提出 （1）连续体上无自然的结点和单元分界面，必须人为地用若干个离算点将连续体划分为有限个单元 （2）没有通用的公式能描述节点位移和单元内任一点的位移及应力关系，应该选择适当的位移函数来建立这种关系 2. 要求 （1）几何近似：计算模型应在几何上与原型和结构近似 （2）物理近似：离算的单元的物理力学性质与原型近似 （3）边界条件近似 3. 离算化过程及内容 （1）确定计算模型几何形状和尺寸，包括边界约束条件和荷载情况； （2）选择单元类型，用有限个单元划分求解区域； （3）单元和结点编号 为保证单元面积不为负，单元结点应按逆时针编号 4. 定义单元 （1）定义每个单元的节点编号：i,j,k,…; （2）定义每个单元结点在总体坐标系中的坐标 i(x,y) j(x,y) k(x,y) （3）给出每个单元的材料性质： 用于平面问题离散化单元主要有：三角形三结点单元（简称三角形单元）、三角形6节点单元、四节点四边形单元、八节点四边形单元、任意四边形单元 5. 划分求解域注意事项 （1）曲线边界，可用多段直线模拟； （2）每个三角形（四边形）单元，边长比不应过于悬殊，以免计算误差过大或者奇异； （3）任一单元的顶点必须是相邻单元的顶点，不能是内点； （4）重点区域加密处理； （5）荷载突变处（均布荷载及集中荷载）加密处理，突变处需设置结点； （6）材料性质差异悬殊处应分别划分单元； （7）充分利用对称性进行简化处理； （8）兼顾精度和经济性。 3.2单元的位移函数及插值函数 根据有限元法的基本思路，将弹性体离散成有限个单元体的组合，以结点的位移作为未知量。弹性体内实际的位移分布可以用单元内的位移分布函数来分块近似地表示。在单元内的位移变化可以假定一个函数来表示，这个函数称为单元位移函数、或单元位移模式。 对于弹性力学平面问题，单元位移函数可以用多项式表示， （3.1） 多项式中包含的项数越多，就越接近实际的位移分布，越精确。具体取多少项，由单元形式来确定。即以结点位移来确定位移函数中的待定系数。 图3.1 如图3.1所示的3结点三角形单元，结点I、J、M的坐标分别为、、，结点位移分别为、、、、、。六个节点位移只能确定六个多项式的系数，所以3结点三角形单元的位移函数如下， （3.2） 将3个结点上的坐标和位移分量代入公式（3.2）就可以将六个待定系数用结点坐标和位移分量表示出来。 将水平位移分量和结点坐标代入（3.2）中的第一式， 写成矩阵形式， （3.3） 令 ， 则有 （3.4） ，A为三角形单元的面积。 [T]的伴随矩阵为， （3.5） 令 （3.6） 则 （3.7） 同样，将垂直位移分量与结点坐标代入公式（3.2）中的第二式，可得， （3.8） 将（3.7）、（3.8）代回（3.2）整理后可得， 令 （下标i，j，m轮换） 可得 （3.9） 单元内的位移记为 单元的结点位移记为 单元内的位移函数可以简写成， （3.10） 把[N]称为形态矩阵，Ni称为形态函数。 选择单元位移函数应满足以下条件： 反映单元的刚体位移与常量应变。 2）相邻单元在公共边界上的位移连续，即单元之间不能重叠，也不能脱离。 由（3.2）可以将单元位移表示成以下的形式， 反映了刚体位移和常应变。 单元位移函数是线性插值函数，因此单元边界上各点的位移可以由两个结点的位移完全确定。两个单元的边界共用两个结点，所以边界上的位移连续。 形态函数Ni具有以下性质： 在单元结点上形态函数的值为1或为0。 2）在单元中的任意一点上，三个形态函数之和等于1。 用来计算三角形面积时，要注意单元结点的排列顺序，当三个结点i，j，m取逆时针顺序时，；当三个结点i，j，m取顺时针顺序时，。 例题：如图3.2所示等腰三角形单元，求其形态矩阵[N]。 解： 由 在公式中轮换下标可以计算得 ，， ，， ， 三角形积为 形态函数为 形态矩阵为 三角形面积的计算公式可得， 如果把三个结点按顺时针方向排列，即i（a，0），j（0，0），m（0，a） 3.3单元载荷移置 有限元法的求解对象是单元的组合体，因此作用在弹性体上的外力，需要移置到相应的结点上成为结点载荷。载荷移置要满足静力等效原则。静力等效是指原载荷与结点载荷在任意虚位移上做的虚功相等。 单元的虚位移可以用结点的虚位移表示为， （3.11） 令