有限元基础理论教程 lecture05幻灯片.pptVIP

2.4 单元刚度矩阵 2.5单元刚度矩阵的性质与物理意义 假设单元的结点位移如下： 单元刚度矩阵的性质： 1）对称性 2）奇异性 3）对角线上的主元恒正 2.6整体分析 整体分析包括以下4个步骤： 1）建立整体刚度矩阵， 2）根据约束条件修改整体刚度矩阵， 3）解方程组，求出结点的位移， 4）根据结点位移，求出单元的应变和应力。 建立整体刚度矩阵的基本方法是刚度集成法，即整体刚度矩阵是单元刚度矩阵的集成。 编程实现整体刚度矩阵集成 NEC(i,1)，表示单元(i)中I结点的整体编号； NEC(i,2)，表示单元(i)中J结点的整体编号； NEC(i,3)，表示单元(i)中M结点的整体编号。 分块矩阵[Krs]中的系数按以下地址加入整体刚度矩阵： 在编制有限元程序时，可以用一个二维数组来保存单元结点的局部编号与该结点在整体结构中的总体编号之间的关系。假定最大单元数目为NELEM，用Fortran语言可以定义以下数组， DIMENSION NEC(NELEM,3) 数组的行与单元对应，数组的第1列为单元中I结点的整体编号，第2列为单元中J结点的整体编号，第一列为单元中M结点的整体编号。例如， NEC(2,1)=5，表示单元(2)中I结点的整体编号为5； NEC(2,2)=2，表示单元(2)中J结点的整体编号为2； NEC(2,3)=4，表示单元(2