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