1、,天津科技大学,有限元分析,Finite Element Analysis,内容 平面矩形单元与平面等参单元1 平面矩形单元分析2 等参单元的概念3 平面四节点等参单元分析要求 理解：平面矩形单元分析流程平面矩形单元与三角形单元的比较等参单元的目的和实现方法掌握： 等参单元分析的基本思想课后作业 收集、阅读平面8节点等参元分析资料,回顾,连续体有限元分析的基本流程,整体 离散,单元分析,单元 组装,引入边界条件整体解算,连续体结构,回顾,连续体有限元分析中的几次近似,1. 逼近性离散（网格剖分）,2. 单元分片插值（单元分析）,几何形状，特别是边界的逼近程度。,插值多项式的项数的截取。,3.
2、其它,回顾,连续体有限元分析中的几次近似,1. 逼近性离散（网格剖分）,提高精度的办法：,减小单元尺寸,改变单元形状,回顾,连续体有限元分析中的几次近似,2. 单元分片插值（单元分析）,提高精度的办法：,增加插值多项式的项数,线性插值函数,如何增加？,增加节点,一、平面矩形单元,矩形单元也是常用的单元之一，由于采用了比常应变三角形单元更高次数的位移模式，故可以更好地反映弹性体的位移状态和应力状态。,如图所示四节点矩形单元，记单元的节点位移向量 和节点力向量 为：,为了能推导出简洁的结果，在这里引入无量纲坐标：,单元位移场,由图可以看出，节点条件共有8个，因此，x和y方向的位移场可以各有4个待定
3、系数，可以取以下多项式作为单元的位移场模式：,它们是具有完全一次项的非完全二次项，其中以上两式中右端的第四项 是考虑到x和y方向的对称性而取的。,由节点条件，在 处，有,回代，可以求解出待定系数 ，然后整理可得,其中，N为单元的形函数矩阵，,如以无量纲坐标系来表达，则上式可以写成,其中：,单元应变场,根据单元的位移场函数式，由几何方程可以得到单元的应变场表达式，,记为：,这里，B矩阵称为几何矩阵。B矩阵可以表示为分块矩阵的形式,其中,注意：矩形单元的应变场为一次线性函数。,单元应力场,由物理方程及应变矩阵，可以得到单元的应力场表达式，,其中 为应力矩阵，D称为弹性矩阵，对于平面应力问题，,注意
4、：矩形单元的应力场为一次线性函数。,将应力矩阵表示为分块矩阵的形式,其中：,对于平面应变问题，只需将E换为 ， 换为 。,单元刚度矩阵,和三角形单元一样，可以根据虚功理导出节点位移向量和节点力向量之间关系，即单元的刚度矩阵 ，可以将其写成分块的形式。,其中,对于平面应力问题，如果单元厚度t为常数，则可得刚度矩阵的显式形式：,积分得：,平面矩形单元小结 优点：矩形单元的应力、应变为一次线性函数，精度要比三角形三节点高； 不足：实际问题很难用4节点矩形单元划分，特别是边界适合性不强； 问题：能否构造一种任意四边形单元，则在提高精度得前提下，边界适应性还强？,等参单元（iso-parametric
5、element）的概念,二、平面等参单元,问题：能否利用规则的平面矩形单元的结果来研究不规则的任意四边形单元的计算公式？,思路：任意直四边形可看成是正四边形(常称为母元)的变形，由于正四边形（母元）的位移函数、单刚矩阵均已得到，则可利用正四边形单元的结果研究任意四边形。,重点：1）构造任意四边形与母元间的坐标（形状）变换关系,2）利用坐标变换关系和母元的计算公式，推导任意四边形的单刚矩阵,（包括母元位移函数、应变矩阵、刚度矩阵转换过程中的导数、积分计算）,等参单元分析范例平面4节点等参单元,1、等参变换（坐标映射）,目的：建立矩形母单元与任意四边形单元的坐标映射关系,求：,解法：插值,代入4个
6、角点坐标，确定系数。,求出待定系数，得,其中：,同矩形单元位移形函数,将四角点的局部坐标代入,2、等参单元位移函数,从坐标变换可知,等参单元位移与母元间位移仅相差坐标变换式,而母元单元内任意点P的位移函数,Ni同矩形单元位移形函数，即与坐标变换形函数相同，故得名等参单元。,3、等参单元应变矩阵,由几何方程，得,新问题：形函数是局部坐标的函数，而局部坐标又是整体坐标的函数，故：,3、等参单元应变矩阵,称为雅克比矩阵，且,3、等参单元应变矩阵,4、等参单元应力矩阵,5、等参单元刚度矩阵,求微小平行四边形面积,注：等参单元的刚度积分一般很难有解析式，必须进行数值积分，目前普遍采用高斯数值积分法（略）
7、。,d，d在（x，y）系中的分量为,等参单元小结,1、等参单元存在的充要条件是|J|0 为了保证能进行等参变换(即总体坐标与局部坐标一一对应)，通常要求总体坐标系下的单元为凸，即不能有内角大于或等于或接近180度情况。,等参单元小结,2、等参单元的优点是当单元边界呈二次以上的曲线时，容易用很少的单元去逼近曲线边界。,3、前述推导要求：保持坐标变换中几何模式阶次与描述单元位移函数中形函数的阶次相同，故被称为等参元。如取坐标变换的几何模式阶次较单元的位移函数的阶次高，则称此单元为超参元，反之，为亚参元。,上机实验安排,时间：第8周 星期三 第12节（08:0009:40） 地点：C302 上机实验题目：结构、连续体Ansys有限元分析（2） 实验目的：Ansys结构有限元分析 实验内容：开始完成课程大作业,作 业,阅读教材平面8节点等参元的分析流程,再 见,
