弹性力学与有限元完整版.pptVIP

三角形单元的位移模式假定为位移模式:位移模式矩阵表达:位移模式通式——单元内任一点的位移；——单元的结点位移向量；——单元的形函数矩阵。形函数的性质面积坐标形函数与面积坐标的关系三角形的面积位移模式反映了单元内任意一点的位移与结点位移之间的关系。是有限元计算精度的关键。§4单元分析4.1单元上任意一点的应变4.2单元上任意一点的应力4.3单元的能量4.4单元刚度矩阵的性质4.1单元上任意一点的应变几何方程或写成通式[B]矩阵叫做单元几何矩阵,反映了单元内任意一点的应变分量与结点位移之间的关系几何矩阵[B]中的每个元素,均为常数,它们由结点坐标确定。单元内任意一点P的应变分量与坐标（x,y）无关,说明单元中应变是常量。4.2单元上任意一点的应力物理方程[D]——弹性矩阵平面应力问题平面应变问题[S]叫做应力矩阵4.3单元的能量1.单元的应变能一维问题应变能密度为平面问题应变能密度为称为单元刚度矩阵,简称单刚,它反映了单元应变能与单元结点向量之间的关系。2.外力势能(1)、体力势能(2)、面力势能(3)、集中力势能3、单元的总势能4.4单元刚度矩阵的性质某单元的刚度矩阵，仔细看看，会发现该矩阵有哪些特点？4.4单元刚度矩阵的性质1、对称性单元刚度矩阵是对称方阵，其元素都对称于主对角线。2、奇异性单元刚度矩阵中任意一行或列元素之和为零。其物理意义是在没有给单元施加任何约束时，单元可有刚体运动，位移不能唯一的确定。3、主对角线元素恒为正值主对角线元素是正值说明结点位移方向与施加结点荷载的方向是一致的。4.单元刚度矩阵与单元位置无关单元刚度矩阵与单元位置无关，也就是单元在平移时，[K]e不变；单元结点排列顺序不同时，[K]e中元素大小不变，而排列顺序相应改变。§5等效节点载荷弹性体所受外力包括体积力、表面力、集中力。分别作用在弹性体内部、物体表面上、物体的一个点上。载荷列阵{R},是由弹性体的全部单元的等效节点力集合而成,是将全部载荷转移到单元的节点上,它们的作用位置发生了变化——载荷移置。它们的作用效果是等效的,故称等效节点力向量{R}e。各种载荷分别移置到节点上,再逐点加以合成求得单元的等效结点载荷。1.体力等效结点载荷自重情况下:y0xijmy0xijm2.面力等效结点载荷6.1结构的结点位移向量假设弹性体被划分为N个单元和n个节点,整个弹性体的节点位移向量{?}2n×1整个弹性体的载荷列阵{R}2n×1§6整体分析矢量——有方向性,外力、应力,不能直接相加标量——没有方向,只有大小,可以相加。弹性体的能量是标量,可以直接相加。6.2结构的总势能单刚的扩充为了实现上述运算扩展——结构的总刚度矩阵——结构总的体力列阵——结构总的面力列阵——结构总的集中力列阵结构的总势能6.3整体刚度矩阵形成方法123421q图5组装总刚[k]的一般规则:1.当[krs]中r=s时,该点被哪几个单元所共有,则总刚子矩阵[krs]就是这几个单元的刚度矩阵子矩阵[krs]e的相加。2.当[krs]中rs时，若rs边是组合体的内边，则总体刚度矩阵[krs]就是共用该边的两相邻单元单刚子矩阵[krs]e的相加。3.当[krs]中r和s不同属于任何单元时,则总体刚度矩阵[krs]=[0]。下面，我们考查一个组装总刚的实例:1.整体刚度矩阵及载荷列阵的组集根据叠加原理,整体结构的各个刚度矩阵的元素显然是由有关单元的单元刚度矩阵的元素组集而成的,为了便于理解,现结合图5说明组集过程。图中有两种编码：一是节点总码：1、2、3.4；二是节点局部码，是每个单元的三个节点按逆时针方向的顺序各自编码为1，2，3。图中两个单元的局部码与总码的对应关系为:单元1：1，2，3 1，2，3单元2：1，2，33，4，1或：单元1：1，2，3 1，2，3单元2：1，2，31，3，4单元e的刚度矩阵分块形式为:整体刚