[工学]平面有限元.ppt

* 组集整体刚度矩阵 由于[Krs]=[Ksr]T，又单元1和单元2的节点号按123对应341，则可得： 按刚度集成法可得整体刚度矩阵为： * * 所以组集的整体刚度矩阵为： * 先求出各单元的应力矩阵[S]1、[S]2，然后再求得各单元的应力分量： 6. 计算各单元应力矩阵，求出各单元应力 单元应力可看作是单元形心处的应力值。 * 7. 引入约束条件，修改刚度方程并求解 根据约束条件：u1 =v1=0；v2=0；u4=0和等效节点力列阵： ，并代入刚度方程： ，划去[K]中与0位移相对应的1，2，4，7的行 和列，则刚度方程变为： 求解上面方程组可得出节点位移为： 所以 * 例4-2 图4-16所示为一平面应力问题离散化以后的结构图，其中图（a）为离散化后的总体结构，图（b）为单元1，2，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 图 4-16 计算实例2的结构图 * 　　首先求确定各单元刚度所需的系数 　 及面积A，对于单元1，2，4有： 解： 对于单元3有: * 其次，求出各单元的单元刚度矩阵。对于1，2，4单元，其单元刚度矩阵为： i j m i j m * 各单元的节点编号与总体结构的总编号之间的对应关系见表4-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 表4-2 * 将各单元刚度矩阵按节点总数及相应的节点号关系扩充成12*12矩阵,分别如下: * * * * 将扩充后的各单元刚度矩阵相加,得总体刚度矩阵K，即： * 所以结构总方程为： 其中 考虑到边界条件： * 用对角元乘大数法消除奇异性后的结构总体方程为： * 由以上方程解得的各节点的位移为： * 然后将相应的节点位移代入公式，可分别求得各单元的应 变和应力。 对于单元1： * 对于单元2： * 对于单元3： * 对于单元4： * ⒊ 体积力的等效载荷列阵{P} 　　与表面力的情况类似，体积力的等效载荷列阵也是由单元体积力的等效节点力在各节点处合成以后，按节点号码顺序排列而成，即 ( j ) 式中单元e的体积力的等效节点力为 ( k ) * y 0 x i j m 设在单元ijk上受有分布体力。取微元 体 ，则此微元体可看成一个集 中力： 体力等效节点力： 图 4-8 由（4-47）式，并对三 角形单元面积积分，即 可得到体力位置到节点 上的等效节点力列阵： * 当单元体是均质、等厚、比重为 时，则x=0，y=- 故有： ? 单元面积 * 矩形单元也是一种常用的单元，它采用了比常应变三角形单元次数更高的位移模式，因而可以更好地反映弹性体中的位移状态和应力状态。 矩形单元1234如图4-9所示，其边长分别为2a和2b，两边分别平行于x、y轴。若取该矩形的四个角点为节点，因每个节点位移有两个分量，所以矩形单元共有8个自由度。采用§4-2节中的方法，同样可以完成对这种单元的力学特性分析。然而，如果我们引入一个局部坐标系?、?，那么就可以推出比较简洁的结果。 第六节 矩形单元 图4-9 矩形单元12