有限元分析第4章平面问题有限单元法.ppt
,李建宇天津科技大学,有限元分析,Finite Element Analysis,内容 Chp.4 平面问题的有限元法 2 4.3 形函数的性质 4.4 整体刚度矩阵的组装及其性质 4.5 等效节点载荷向量要求 理解：连续体有限元分片插值的含义；形函数的功能及其性质；基于虚功方程的整体刚度矩阵的组装；基于虚功方程的等效节点载荷向量的生成 掌握：常应变三角形单元组装技术；常应变三角形单元等效节点载荷向量的生成 课后作业 推导常应变三角形单元总刚和等效载荷向量,有限元分析43,4,5,回顾,连续体有限元分析的基本流程,整体离散,单元分析,单元组装,整体解算,连续体结构,离散化方法,回顾,几何实体的逼近性离散,三角形单元分析,回顾,目标：对三角形单元，建立节点位移与等效节点力之间的转换关系。,？,回顾,三角形单元分析,单元分析的流程,解决办法：插值(分片插值的提法),形函数矩阵,形函数（shape function）及其性质,回顾,下标i,j,m 轮换,性质：,其中：,（a）Ni(x,y)在i点的节点位移为1，其它节点为0。,（b）单元中任一点各形函数的和为1。,力学意义：固定 j，m节点，使i节点发生位移1，则单元内各点的位移场为Ni(x,y)。,力学意义：令单元发生刚体位移u0，则单元内各点的位移均为u0，即,单元分析流程,回顾,解决办法：弹性力学几何方程,B矩阵称为应变矩阵,得,代入,该单元为常应变单元,单元分析流程,回顾,解决办法：弹性力学物理方程,得,代入,S矩阵称为应力矩阵。,例：对于平面应力问题,代入,得,其中,单元分析流程,回顾,解决办法：单元平衡分析,平面问题虚功原理,代入,得,内力虚功,外力虚功,外力虚功内力虚功,单元刚度方程,单元刚度方程建立了单元的节点力与节点位移之间的关系，称为单元刚度矩阵。它是66矩阵，其元素表示该单元的各节点沿坐标方向发生单位位移时引起的节点力，它决定于该单元的形状、大小、方位和弹性常数，而与单元的位置无关，即不随单元或坐标轴的平行移动而改变。,一、单元刚度矩阵的物理意义及其性质,已经求出了下列关系,一、单元刚度矩阵的物理意义及其性质,结点力和节点位移的关系：(以简单平面桁架为例)平面问题中，离散化的单元组合体极为相似，单元组合体在节点载荷的作用下，节点对单元、单元对节点都有作用力与反作用力存在，大小相等方向相反，统称为节点力。节点力和节点位移的关系为单元刚度方程：,一、单元刚度矩阵的物理意义及其性质,单元刚度矩阵的物理意义：将 写成分块矩阵 写成普通方程 其中 表示节点s(s=i,j,m)产生单位位移时，在节点r(r=i,j,m)上所需要施加的节点力的大小。,一、单元刚度矩阵的物理意义及其性质,单元刚度矩阵的物理意义：将节点力列向量 与节点位移列向量 均扩展成(61)阶列矩阵，单元刚度矩阵相应地展开成(66)阶方阵：元素K的脚码，标有“-”的表示水平方向，没有标“-”的表示垂直方向。,一、单元刚度矩阵的物理意义及其性质,单元刚度矩阵的物理意义：单元刚度矩阵的各元素的物理意义：表示节点s(s=i,j,m)在水平方向、垂直方向产生单位位移时，在节点r(r=i,j,m)上分别所要施加的水平节点力和垂直节点力的大小。例如 表示节点j在垂直方向产生单位位移时，在节点i所需要施加的水平节点力的大小。,一、单元刚度矩阵的物理意义及其性质,单元刚度矩阵的性质：1)对称性：是对称矩阵 2)奇异性：是奇异矩阵 单元刚度矩阵所有奇数行的对应元素之和为零，所有偶数行的对应元素之和也为零。由此可见，单元刚度矩阵各列元素的总和为零。由对称性可知，各行元素的总和也为零。,一、单元刚度矩阵的物理意义及其性质,单元刚度矩阵的性质：例题：求下图所示单元的刚度矩阵，设,1、求B2、求D3、求S4、求,回顾,连续体有限元分析的基本流程,整体离散,单元分析,单元组装,整体解算,连续体结构,二、单元组装技术,位移协调条件：各单元共享节点的位移相等,节点平衡条件：各单元内力与节点外力构成平衡力系,二、单元组装技术,以i点为例，利用虚功原理建立平衡方程，设虚位移,各单元i节点等效内力的虚功为:,各单元i节点等效外力的虚功为:,二、单元组装技术,由单元刚度方程：,对每一个单元，将i节点虚位移扩展到单元全部节点，如，,则,Te表示节点位移提取矩阵,对i节点，有,同理，可建立其它节点的平衡方程。,二、单元组装技术,对i节点，有,将每一个单元的节点位移（包括虚位移）扩展到全部节点位移向量，如,则，,Te表示单元节点位移提取矩阵,二、单元组装技术,对所有节点，有,由虚位移的任意性，整理得,总体刚度矩阵,总体刚度方程,二、单元组装技术,总体刚度方程,等效节点外载荷向量。,三、单元等效节点外载荷向量,连续弹性体离散为单元组合体时，为简化受力情况，需把弹性体承受的任意分布的载荷都向节点移置(分解)，而成为节点载荷。将载荷移置到节点上，必须遵循静力等效的原则。静力等效是指原载荷与节点载荷在任意虚位移上做的虚功相等。在一定的位移模式下，移置结果是唯一的，且总能符合静力等效原则。,三、单元等效节点外载荷向量,在线性位移模式下，对于常见的一些载荷，可以通过简单的虚功计算，得出所需的载荷列矩阵。,均质等厚度的三角形单元所受的重力，把1/3的重力移到每个节点,三、单元等效节点外载荷向量,例：,总载荷的2/3移置到节点i，1/3移置到节点j，与原载荷同向,三、单元等效节点外载荷向量,载荷向节点的移置，可以用普遍公式来表示。体力的移置分布面力的移置在线性位移模式下，用直接计算法简单；非线性模式下，要用普遍公式计算。,课后作业,如图（a）所示一高深悬臂梁，在右端部受集中力F作用，材料弹性模量E、泊松比v1/3，悬臂梁的厚度（板厚）为t，如图（b）所示有限元模型，试按平面应力问题，求结构总体刚度矩阵和总体载荷向量。,再 见,
