弹性力学的有限元求法.ppt

§6.1 有限元网格(Finite Element Network) 1. 常用单元 (1) 线性单元(linear elements)：三角形、矩形或其他四边形。形函数是线性函数，即单元内任一点的坐标可用单元节点坐标的线性函数来表示。 (2) 抛物形单元(parabola elements)：除了角上有节点（主外节点），边缘上也有节点（副外节点）的单元。有时在内部也有节点（内节点），直边或曲边均可。单元内任一点的坐标，可用一个抛物线内插法来求得。 2. 单元划分要注意的几个问题 (1) 相邻两个单元的节点要与节点重合（外节点与外节点、内节点 与内节点），不能与无节点边重合。 (2) 单元不必是相同尺寸，应力有突变的地方，单元划分应较小。 (3) 任何一个单元必须只能在一种材料区，即它不能跨越两种材料 区。 (4) 同一单元的各个边长，一般不要相差太大。 位移函数构造和收敛性要求 单元中的位移模式采用待定系数的有限多项式做为近似函数，有限多项式选取的原则： －待定系数是由结点位移条件确定的，三角形6结点位移，矩形单元待定系数8结点位移。 －在选择多项式时，必须要选择常数项和完备的一次项。因为这两项可以反映单元刚体位移和常应变的特性。 －选择多项式应由低阶到高阶，并具有坐标对称性。 6.2单元位移函数 弹性体内的实际位移分布可以用单元内的位移分布函数来分片近似地表示。在单元内的位移变化可以假定一个函数来表示，这个函数称为单元位移函数、或单元位移模式。单元位移函数可以用多项式表示， 形态函数Ni具有以下性质： 1）在单元结点上形态函数的值为1或为0。 2）在单元中的任意一点上，三个形态函数之和等于1。 形态函数的几何意义 2.3单元载荷移置 有限元法的求解对象是单元的组合体，待求解的未知变量都定义在单元的结点上，因此作用在弹性体上的外力，需要移置到相应的单元结点上成为结点载荷。 载荷移置要满足静力等效原则。静力等效是指原载荷与结点载荷在任意虚位移上做的虚功(virtual work)相等。 4.4 单元刚度矩阵的求解过程 4.4.1 单元刚度矩阵 * * 图5.3 线形单元及节点位置 图5.4 抛物线形单元及节点位置 1 常数项1 x y 线性项3 x2 xy y2 二次项6 x3 x2y xy2 y3 立方项10 x4 x3y x2y2 xy3 y4 四次项 15 3结点三角形单元内位移由结点的6个位移分量来确定。六个位移分量只能确定六个多项式的系数，所以3结点三角形单元的位移函数如下， (2-6) 将3个结点上的坐标和位移分量代入公式（2-6）就可以将六个待定系数用结点坐标和位移分量表示出来。首先计算位移分量u的系数， (2-7) A为三角形单元的面积 (2-8) [T]的伴随矩阵为， (2-9) (2-10) i,j,m坐标轮换 (2-11) (2-12) 同样，将垂直位移分量与结点坐标代入公式（2-6）中的第二式，可得， 将（2-11）、（2-12）代回（2-6）整理后可得， （下标i，j，m轮换） 令 (2-13) 单元内的位移记为 单元的结点位移记为 单元内的位移函数可以简写成， (2-14) 把[N]称为形态矩阵，Ni称为形态函数。 例2.1、3结点三角形单元如图所示，求其形态矩阵[N]。 结点坐标为， 三角形积为 形态函数为 形态矩阵为 用来计算三角形面积时，要注意单元结点的排列顺序。当三个结点i，j，m取逆时针顺序时，用行列式计算出的三角形面积为正值。 如果把三个结点按顺时针方向排列，即i(a,0)、j(0,0)、m(0,a)，行列式的计算结果为负值。 例2.2、3结点三角形单元如图所示。结点位移如下， ui=0.1a, vi=0.05a, uj=0.15a, vj=0.1a, um=0.05a, vm=0。单元内一点p的坐标为(0.25a,0.5a)，求p点的位移。 对于任意一点p(x,y)， 假定结点位移： 按照上面给出的结点位移，单元实际上产生了刚体位移，对于单元中的任意一个点u(x,y)=1。所以， 任意一点P的形态函数Ni是点P与结点I的对边所构成的三角形面积与整个单元面积之比。 结点的虚位移 单元的虚位移 令结点载荷为 处于平衡状态的弹性体，真实位移