第3章-平面问题的有限元法

第3章平面问题的有限元法3.1结构的离散化3.2三角形常应变单元的位移模式和形函数3.5整体分析3.6等效节点载荷计算3.8有限元分析的实例3.3单元刚度矩阵3.4单元位移函数的选择原则3.7约束条件的处理将连续体变换为离散化结构将连续体划分为有限多个、有限大小的单元，并使这些单元仅在节点处连结起来，构成所谓“离散化结构”。(c)深梁（离散化结构）3.1结构的离散化离散化要注意:1.单元形状的选择:平面问题的单元，按其几何特性可分为两类：以三节点三角形为基础；以任意四边形为基础。较高精度的三角形等参数单元；运用非常广泛的四边形等参数单元。这两类都可以增加节点也构成一系列单元：首选三角形单元和等参数单元。2.对称性的利用利用结构和载荷的对称性：如结构和载荷都对于某轴对称，可以取一半来分析；若对于x轴和y轴都对称，可以取四分之一来分析。3.单元的划分原则通常集中载荷的作用点、分布载荷强度的突变点、分布载荷与自由边界的分界点，支承点都应取为节点单元的形状和尺寸可以根据要求进行调整。对于重要或应力变化急剧的部位，单元应划分得小些；对于次要和应力变化缓慢的部位，单元可划分得大些；中间地带以大小逐渐变化的单元来过渡。单元的划分原则单元数量要根据计算精度和计算机的容量来决定。在保证精度的前提下，尽可能减少单元数量。不要把不同厚度或不同材料的区域划分在一个单元里。单元的划分原则根据误差分析，应力及位移的误差都和单元的最小内角正弦成反比，所以单元的边长力求接近相等。即单元的三（四）条边长尽量不要悬殊太大。4.节点的编号应尽量使同一单元的节点编号相差小些，以减少整体刚度矩阵的半带宽，节约计算机存储。上图，节点顺短边编号为好。3.2三角形常应变单元的位移模式和形函数首先以平面单元中最基本的三节点三角形单元为例，介绍有限元法。单元分析的步骤可表示如下：节点位移内部各点位移→应变应力→→节点力单元分析分为四步求出相邻各量之间的转换关系,综合起来,得出由节点位移求节点力的转换关系:eeekFek单元刚度矩阵位移模式1.位移模式Tiiiewvu...单元的若干个节点有基本未知量，即位移模式:单元内任一点的位移表达式，假定为坐标的简单函数。反映单元的位移分布形态，是单元内的插值函数。在节点处等于该节点位移。位移模式可表示为：eNfN为形态矩阵（形函数矩阵）平面问题每个节点位移分量有两个，所以整个单元有6个节点位移分量，即6个自由度。单元节点位移列阵:TmmjjiiTmTjTievuvuvu三角形单元有6个自由度，内部任一点的位移是由6个节点位移分量完全确定的，位移模式中应含有6个待定系数，所以位移模式可取为：ayxvyxu。654321,位移函数一般用多项式来构造。位移模式:单元内任一点的位移表达式，假定为坐标的简单函数。反映单元的位移分布形态。在弹性体内，位移变化非常复杂。有限元法将整个弹性体分割成许多小单元，在每个单元内采用简单的函数来近似表达单元的真实位移，将各单元连接起来，便可近似表达整个弹性体的真实位移函数。这种化整为零、化繁为简的方法，正是有限元法的精华。假设节点i，j，m的坐标分别为(xi,yi),(xj,yj),(xm,ym)2.形函数联立求解左边3个方程，得：其中A为三角形单元的面积注意：为了使得出的面积值不为负值，节点i,j,m的次序必须是逆时针。至于将那个节点作为起始点i则没有关系。同理，求解右边的三个方程，得到a4，a5，a6，解得：mmmmjjjjiiiivycxbavycxbavycxbaAv21mjimjimmjjixxcyybyxyxa11,11,式中：i，j，m轮换mmmmjjjjiiiiuycxbauycxbauycxbaAu21整理后得：)m,,(A21轮换令jiycxbaNiiiiiimmjjiiiimjjiivNvNvNvNvuNuNuNuNum其中Ni，Nj，Nm是坐标的线性函数，反应了单元的位移形态，称为形（状）函数。eekjiNINININ写成矩阵形式式中：I二阶单位阵，N形函数矩阵f3.三角形面积坐标定义：在三角形内任一点P，向三个角点（节点）连线，将原三角形分割成三个子三角形，设子三角形的面积分别是：Ai，Aj，Am，则：AALiiAALjjAALmm即面积坐标定义为子三角形与原三角形面积之比；记为：P（Li，Lj，Lm）。面积坐标的性质：AAAAmji1.1mjiLLLLi，Lj，Lm中只有两个是独立的。2.三角形三个角点处)0,0,1(i)1,0,0(m)0,1,0(j3.三条边上i-j:Lm=0j-m:Li=0m-i:Lj=0形心处：31mjiLLL推论：三角形内一条平行于三角形任一边的直线上的各点，具有相同的与该边对应的坐标值。iiiHhL面积坐标与直角坐标的转换：)(2111121ycxbayxyxyxAiiimmjji（i,j,m))(21ycxbaAAALiiiii(i,j,m)因此：iiNLmmNLjjNL即三角形面积坐标就是三角形相应的形函数。mmjjiiyxyxyxA11121所以，位移模式也可以用面积坐标表示为：iimmjjiiiimjjiivLvLvLvLvuLuLuLuLum)(21ycxbaAAALiiiii(i,j,m)将面积坐标的表达式：写成矩阵形式：yxcbacbacbaALLLmmmjjjiiimji121求逆得：mjimjimjiLLLyyyxxxyx1111第1行展开为面积坐标性质1，第2行和第3行展开即为局部的面积坐标和整体直角坐标的关系：iimmjjiiiimjjiiyLyLyLyLvxLxLxLxLxm例题下图为一平面应力的直角三角形单元，直角边长均为a,厚度为t，弹性模量为E,泊松比μ=0.3,求形函数。1.单元应变emmjjiimjimjixyyxbcbcbccccbbbAxvyuyvxu00000021eB),,(0021mjibccbABiiiiimjiBBBB应变矩阵应变矩阵为常量，单元内应变是常数3.3单元刚度矩阵eeSDBDmjimjiSSSBBBDDBS2.单元应力S——称为应力转换矩阵)(),,(2121)1(22fmjibμcμcμbμcbAμEiiiiii。iiDBS应用平面应力问题的弹性矩阵：2100010112ED阵：平面应力问题，弹性矩应变矩阵为常量，单元内应力也是常数，相邻单元的应变与应力将产生突变，但位移是连续的。)1(2210001-10-112-11-1））（（）（阵：平面应变问题，弹性矩ED能量转换与守恒定律，是自然界基本的运动规律之一。实功原理：处于平衡状态的可变形固体，在受外力作用而变形时外力对其相应的位移所做的功（实功），等于积蓄在物体中的应变能（实应变能）。能量法的优点：与坐标系的选择无关，因而应用极为广泛。能量法与数学工具—变分法的结合，导出虚位移（虚功）原理，使得用数学分析的方法解决力学问题的理论得到发展而更趋完善。3.虚位移（功）原理TmmjjiievuvuvuTymxmyjxjyixieFFFFFFF单元节点力列阵：单元节点虚位移列阵：节点力在虚位移所做的功：ymmyjjxjjyiixiiFuFvFuFvFuWeTeFW简写为：4.单元刚度矩阵eB单元虚应变：单元内应力在虚应变上所做的功（虚应变能）：tdxdyUxyxyyyAxx)(其中：t为单元厚度eeSDBD单元应力：eTATeTAtdxdyDBBtdxdyUW虚功原理eeeeTekFtdxdyDBBFtdxdyDBBkTe为单元面积A单元刚度矩阵ke取决于单元的大小、方向和弹性常数，而与单元的位置无关，即不随单元或坐标轴的平行移动而改变。对于三角形常应变单元：tADBBkTe单元刚度矩阵为对称矩阵。mmmjmijmjjjiimijiimjiTmTjTiekkkkkkkkkBBBDBBBksrsrsrsrsrsrsrsrrsbbcccbbcbccbccbbAEtk21212121)1(42对于平面应力问题：mjismjir,,;,,其中例题下图为一平面应力的直角三角形单元，直角边长均为a,厚度为t，弹性模量为E,泊松比μ=0.3,求单元刚度矩阵。理论力学中质点、质点系（刚体）的虚位移原理；材料力学中杆件的虚位移原理。弹性力学中的虚位移（虚功）原理：在外力作用下处于平衡状态的变形体，当给与该物体微小位移时，外力总虚功在数值上等于变形体的总虚应变能。虚：微小的、任意的、可能的，变分的思路实功是力在自己产生位移上所做的功，虚功是力在别的（人为的）因素产生的位移上做的功。所谓”虚“并不是虚无，而是可能、虚设的意思。“虚”的表达：δ虚位移（虚功）原理：3.4单元位移函数的选择原则三角形常应变单元简单，精度较差，要提高精度：1.增加单元数目和节点数目；2.采用更高精度的单元。FEM中的一系列工作，都是以位移模式为基础的。所以当单元趋于很小时，即△x,△y→0时，为了使FEM之解逼近于真解，即为了保证FEM收敛性,位移模式应满足下列条件：1.位移模式必须能反映单元的刚体位移。单元位移包含两部分：本单元的形变引起的位移；其他单元的形变引起的位移，即刚体位移。在位移函数中，常数项即提供刚体位移。2.位移模式必须能反映单元的常量应变。单元应变包含两部分：变量应变和常量应变。位移函数的一次项提供常量应变。当单元→0时，单元中的位移和应变都趋近于基本量—刚体位移和常量位移。3.位移模式应尽可能反映位移的连续性使相邻单元之间的位移保持连续，即受力后，相邻单元在公共边界上，即既不互相脱离，也不互相嵌入。使相邻单元在公共节点处具有相同的位移。使单元内部的位移保持连续。位移函数取坐标的单值连续函数。满足条件1、2的单元，称为完备单元；满足条件3的单元，称为协调单元。常采用“帕斯卡三角形”来选取位移模式代数多项式的形式。按照帕斯卡三角形选择位移模式的原则：1.多项式的阶次及项数，由单元的节点数目和自由度数目来决定。保证多项式中的待定系数同单元的自由度数目相一致，以避免在确定待定系数时增加困难。2.当高次多项式只选取一部分项时，应遵循“对称性”原则，即取其最高次中的位置对称的相应项，以保证在各坐标轴方向上具有相同的精度。3.应满足完备性和协调性要求。3节点三角形单元：yaxaavyaxaau6543216节点三角形单元：2121121098726524321yaxyaxayaxaavyaxyaxayaxaau4节点四边形单元：xyayaxaavxyayaxaau876543213.5整体分析结构的整体分析是将离散后的所有单元通过节点连接成原结构，进行分析。分析过程是将所有单元平衡方程组集成整体平衡方程，引进边界条件后求解整体节点位移向量。整体平衡方程：F=KδK为整体刚度矩阵TTnTTnTTnTTnFFFF21