第三章 平面问题的有限元.pptVIP

机械有限元分析 Mechanical Finite Element Analysis and Its Applications 主讲：苏建新 第三章 平面问题的有限元 知识总结 离散过程： 建立坐标系，选择单元类型，划分网格，单元结点编号，给出结点坐标，每一单元结点编号，将位移约束及载荷表示在结点上。 离散应注意问题： 1、对同一问题，采用不同单元精度不同，结点越多，精 度越高，但计算量越大 2、在边界曲折、应力集中处单元的尺寸要小些，同一问题最大与最小单元尺寸倍数不宜过大 3、在集中力作用点以及分布力突变点最好选为结点 4、厚度变化、物性变化也应在划分单元时区分开来 5、单元形状影响计算精度，60度较好 对应的单元1写出相应的单元刚度方程，即单元节点平衡方程 —1号单元第i点（i=1，2，3）节点所受力 为便于组装整体刚度矩阵，将上式以整体节点位移δ表示 K (1) —1号单元的扩大刚度矩阵或称为单元贡献矩阵。 同理，对其他各单元也可以写出单元贡献矩阵，对单元（2） K (2) —2号单元的扩大刚度矩阵或称为单元贡献矩阵。 —2号单元第i点（i=1，2，3）节点所受力 对单元（3） K (3) —3号单元的扩大刚度矩阵或称为单元贡献矩阵。 —3号单元第i点（i=1，2，3）节点所受力 对单元（4） K (4) —3号单元的扩大刚度矩阵或称为单元贡献矩阵。 —3号单元第i点（i=1，2，3）节点所受力 对于任意一个节点，可能承受两种力的作用，一种是其他单元给予该节点的反作用力；另一种是作用在节点上的等效节点力。对整体而言，前者属于内力，后者属于外力，每个节点在两种力的作用下处于平衡。将各节点所受力相加，由于对整体而言，单元给予节点的反作用力属于内力，在相加过程中相互抵消，所以各节点所受力相加的结果只有外力，即等效结点力，从而得到整体节点载荷列阵，即 代入，得到总体刚度矩阵，即 K为总体刚度矩阵，即 可得到整体节点载荷与整体节点位移之间的关系式，即整体有限元方程，即 总体刚度矩阵组装的基本步骤： 1）将单元刚度矩阵中的每个子块放到在整体刚度矩阵中的对应位置上，得到单元的扩大刚度矩阵。 注意对于单元刚度矩阵是按局部编码排列的，即对应单元刚度矩阵中的I,j,m；对于整体刚度矩阵是按整体编码排列的，即按节点号码以从小到大的顺序排列。组装中，必须知道单元节点的局部编码与该节点在整体结构中的整体编码之间的关系，才能得到单元刚度矩阵中的每个子块在整体刚得矩阵中的位置。将单元刚度矩阵中的每个子块按总体编码顺序重新排列后，可以得到单元的扩大矩阵。 2）将全部单元的扩大矩阵相加得到整体刚度矩阵 例3.4 如图所示有限元模型，弹性模量为E，厚度为t，为简化计算取u=0,求整体刚度矩阵 解：该模型中共有6个节点，4个单元，各单元的信息如表所示 K63(4) K65(4) K66(4) K25(3) K23(3) K22(3) K52(2) K54(2) K55(2) K31(1) K32(1) K33(1) K53(4) K55(4) K56(4) K35(3) K33(3) K32(3) K42(2) K44(2) K45(2) K21(1) K22(1) K23(1) K33(4) K35(4) K36(4) K55(3) K53(3) K52(3) K22(2) K24(2) K25(2) K11(1) K12(1) K13(1) 以整体编码表示的单元刚度矩阵子块 i、j、m i、j、m i、j、m i、j、m 局部编码 3、5、6 5、3、2 2、4、5 1、2、3 整体编码 （4） （3） （2） （1） 单元编号 单元刚度矩阵的特点 用有限元方法分析复杂工程问题时，节点的数目比较多，整体刚度矩阵的阶数通常也是很高的。计算量大，浪费资源、降低计算效率，如果能根据整体刚度矩阵的特点进行编程，可减低资源消耗，提高计算效率。必须理解整体刚度矩阵的特点: 1)K是对称矩阵 由整体刚度矩阵的对称性和整体刚度矩阵的组装过程，可知整体刚度矩阵必为对称矩阵。利用对称性，在编写程序时，只存储整体刚度矩阵的上三角或下三角部分即可。 2)K中主对角元素总是正的 主对角元素Knn表示节点在某一方向上产生单位位移，而其他位移均为0时，在节点的这一方向上必须施加的力。