有限元 2-弹性力学平面问题(2.1三角形单元,2.2问题讨论).pptVIP

有限单元法 第2章 弹性力学平面问题 有限单元法 材料力学、结构力学与弹性力学的比较 弹性力学中的基本假定。 弹性力学平面问题 平面应力问题 平面应变问题 位移与应变之间的几何方程 平面应力状态下的应力应变关系 平面应变状态下的应力应变关系 常用的平面单元 第2章 弹性力学平面问题有限单元法 2.1 三角形单元 2.2 三角形单元中几个问题的讨论 2.3 平面问题有限元程序设计? 2.4 矩形单元 2.5 六结点三角形单元 2.6 四结点四边形单元 2.7 八结点曲线四边形等参元? 2.8 几个问题的补充 (1)介质的离散化 (2)单元特性计算 (3)单元组合体的结构分析 第2章 弹性力学平面问题有限单元法 2.1 三角形单元 2.2 三角形单元中几个问题的讨论 2.3 平面问题有限元程序设计? 2.4 矩形单元 2.5 六结点三角形单元 2.6 四结点四边形单元 2.7 八结点曲线四边形等参元? 2.8 几个问题的补充 2.1 三角形单元 三角形单元是有限元分析中的常见单元形式之一，它的优点是: ①对边界形状的适应性较好 ②单刚形式及其推导比较简单 这种单元的计算精度较低，使用的时候必须进行精细的网格划分，但它仍然是一种常用的单元。 一、结点位移和结点力列阵 从某一结构中取出一典型三角形单元，三个节点按逆时针编号为i,j,m，在平面应力问题中，单元的每个结点上有沿x、y两个方向的力和位移。 单元的结点位移列阵 结点力列阵 二、单元位移函数和形状函数 有限单元法是一种近似方法，在单元分析中，首先要求假定（构造）一组在单元内有定义的位移函数作为近似计算的基础。即以结点位移为已知量，假定一个能表示单元内部（包括边界）任意点位移变化规律的函数。 构造位移函数的方法是：以结点(i,j,m)为定点。以位移(ui ,vi , uj,vj,um,vm )为定点上的函数值，利用普通的函数插值法构造出一个单元位移函数。 在平面应力问题中，有u,v两个方向的位移,若假定单元位移函数是线性的，则可表示成: 式中的6个待定常数α1 ,…, α6 可由已知的6个结点位移分量(3个结点的坐标)确定。将3个结点坐标(xi,yi ),(xj,yj ),(xm,ym )代入上式得如下两组线性方程: 联立求解左边的三个方程，可以求出待定常数α1 、α2 、α3 只要A不为0，则可由上式解出 代回（2-1-2)，得 将三角形单元的位移函数用矩阵表示： 三、单元的应变和应力 １、由节点位移求应变──几何矩阵[B] 由弹性力学知，弹性力学平面问题的几何方程： 将 {f}=[N]{d}e代入: ２、应力矩阵[S] 由弹性力学知，由应力求应变的物理方程是: 将 {ε}=[B]{d}e代入上式得: 四、单元刚度矩阵 将几何矩阵[B]和弹性矩阵[D] 代入单元刚度矩阵的一般表达式: 子块: 对此作些改变，便会获得不同的单元刚度矩阵形式，如将位移列阵的排列改为: {d}e=[ui uj um vi vj vm ]T。 （可自行推导，此种单元刚度矩阵的显式可参见《Finite Element Analysis Fundamentals》） 作业3： 1．求形函数Ni(x,y)在三角形形心（xc,yc）上的函数 值。 2．设图中i点有水平位移ui=1，试由单元刚度方程 写出各链杆的反力；并证明各水平、竖向反力之和为0。 3．求图示单元1、2的单元刚度矩阵和应力矩阵(已知E，m，t)。 第2章 弹性力学平面问题有限单元法 2.1三角形单元 2.2 三角形单元中几个问题的讨论 2.3平面问题有限元程序设计? 2.4 矩形单元 2.5 六结点三角形单元 2.6 四结点四边形单元 2.7八结点曲线四边形等参元? 2.8 几个问题的补充 2.2 三角形单元中几个问题的讨论 一、形函数的物理意义 故形函数Ni(x,y)表示当结点i发生单位位移时，在单元内部产生的位移分布状态，形函数Nj ，Nm 亦具有类似的性质。因此，Ni ,Nj ,Nm 称为位移的形状(态)函数，简称形函数，[N]即称形函数矩阵。 二、形函数的几何意义 其中形函数: 若从单元内任意点p(x,y)向各顶点引连线,将其分成三个小三角形，则上式中的行列式恰好是小三角形Pjm面积的2倍 面积坐标 三角形单元形函数的性质 三、有限元的收敛性 解的收敛性也可理解为一个问题的解的精度，较粗的分，影响一个实际问题的解的精度可分成三个方面： 实际物理问题→①理想化力学模