(新)第三章 平面问题的有限单元法——精品.ppt

由于单元2若按341对应单元1的123排码时,则这两个单元刚度矩阵内容完全一样，故有: 3 4 1 3 4 1 返回 组集整体刚度矩阵 由于[Krs]=[Ksr]T，又单元1和单元2的结点号按123对应341，则可得： 按刚度集成法可得整体刚度矩阵为： 返回 返回 所以组集的整体刚度矩阵为： 返回 6. 引入约束条件，修改刚度方程并求解 根据约束条件：u1 =v1=0；v2=0；u4=0和等效结点力列阵： ，并代入刚度方程： ，划去[K]中与0位移相对应的1，2，4，7的行 和列，则刚度方程变为： 求解上面方程组可得出结点位移为： 所以 先求出各单元的应力矩阵[S]1、[S]2，然后再求得各单元的应力分量： 7. 计算各单元应力矩阵，求出各单元应力 单元应力可看作是单元形心处的应力值。 返回 例3-2 图3-16所示为一平面应力问题离散化以后的结构图，其中图（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 图 3-16 计算实例2的结构图 返回 　　首先求确定各单元刚度所需的系数 　 及面积A，对于单元1，2，4有： 解： 对于单元3有: 返回 其次，求出各单元的单元刚度矩阵。对于1，2，4单元，其单元刚度矩阵为： i j m i j m 返回 各单元的结点编号与总体结构的总编号之间的对应关系见表3-2。 对于单元3，其单元刚度矩阵为： i j m i j m 返回 各单元结点号与总体结点号对应表 单元号 1 2 3 4 结点号 节 点 总 编 号 I 1 2 2 3 j 2 4 5 5 m 3 5 3 6 表3-2 返回 将各单元刚度矩阵按结点总数及相应的结点号关系扩充成12*12矩阵,分别如下: 返回 返回 返回 返回 将扩充后的各单元刚度矩阵相加,得总体刚度矩阵K，即： 返回 所以结构总方程为： 其中 考虑到边界条件： 返回 用对角元乘大数法消除奇异性后的结构总体方程为： 返回 由以上方程解得的各结点的位移为： 返回 然后将相应的结点位移代入公式，可分别求得各单元的应 变和应力。 对于单元1： 返回 对于单元2： 返回 对于单元3： 返回 对于单元4： 返回 3.8.1 六结点三角形单元 一、位移模式与形函数 取三角形顶点和边中点作结点，位移模式为： 六结点三角形单元 用面积坐标表示的形函数为： 3.8 平面问题的高次单元 二、应变 六结点三角形单元 3.8.2 矩形双线性单元 矩形单元 矩形单元结点位移、结点力列阵 一、位移模式与形函数 正方形规则单元 正方形单元与矩形单元的关系 形函数的性质：本点处值为1，它点处值为0 4、特殊载荷的移置： y 0 x i j m 已知在ij边受有面力q，则移置到i、j 结点上的等效结点力为： y 0 x i j m 当某一边上有三角形分布的面力时， 可由刚体静力等效直接写出 　　如果任意三角形单元ijk的重心c上受有自重 则按刚体静力等效原理可把W直接移置到i，j，m三个结点上而组成 ： 等效结点力、载荷列阵小结 例:图示三角形单元的i j m的jm边上作用有如图所示线性均匀分布面载荷，三角形单元厚度