课件：弹性力学与有限元完整版.ppt

所以组集的整体刚度矩阵为： 5. 引入约束条件，修改刚度方程并求解 根据约束条件：u1 =v1=0；v2=0；u4=0和等效节点力列阵： ，并代入刚度方程： ，划去[K]中与0位移相对应的1，2，4，7的行 和列，则刚度方程变为： 求解上面方程组可得出节点位移为： 所以 先求出各单元的应力矩阵[S]1、[S]2，然后再求得各单元的应力分量： 6. 计算各单元应力矩阵，求出各单元应力 单元应力可看作是单元形心处的应力值。 例2 图中所示为一平面应力问题离散化以后的结构图，其中图（a）为离散化后的总体结构，图（b）为单元1，2，3，4的结构，图（c）为单元3的结构。用有限元法计算节点位移、单元应变及单元应力（为简便起见，取泊松比 ，单元厚度t=1）。 x y 1 2 3 4 6 5 1 2 3 4 a 3 i j m a a a 1,2,4 i j m 图 计算实例2的结构图 各单元节点号与总体节点号对应表 单元号 1 2 3 4 节点号 节 点 总 编 号 I 1 2 2 3 j 2 4 5 5 m 3 5 3 6 解： 　　1、求确定各单元刚度所需的系数 　 及面积A，对于单元1，2，4有： 对于单元3有: 2、求出各单元的单元刚度矩阵。对于1，2，4单元，其单元刚度矩阵为： i j m i j m 各单元的节点编号与总体结构的总编号之间的对应关系见表2。 对于单元3，其单元刚度矩阵为： i j m i j m 3、总刚 将各单元刚度矩阵按节点总数及相应的节点号关系扩充成12*12矩阵,分别如下: 将扩充后的各单元刚度矩阵相加,得总体刚度矩阵K，即： 4、结构总方程： 其中 考虑到边界条件： 5、位移约束的处理 用对角元乘大数法消除奇异性后的结构总体方程为： 6、节点的位移 由以上方程解得的各节点的位移为： 7、单元应力计算 然后将相应的节点位移代入公式，可分别求得各单元的应 变和应力。 对于单元1： 对于单元2： 对于单元3： 对于单元4： 6.3 整体刚度矩阵形成方法 1 2 3 4 2 1 q 图 5 组装总刚[k]的一般规则： 1. 当[krs]中r=s时，该点被哪几个单元所共有，则总刚子矩阵[krs]就是这几个单元的刚度矩阵子矩阵[krs]e的相加。 2. 当[krs]中r s时，若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，3 3，4，1 或： 单元 1 ： 1，2，3 1，2，3 单元 2 ： 1，2，3 1，3，4 单元e的刚度矩阵分块形式为： 整体刚度矩阵分块形式为： 其中每个子块是按照节点总码排列的。 通常，采用刚度集成法或直接刚度法来