《平面分析-整体刚度矩阵》.pptVIP

5 平面问题有限元分析整体刚度矩阵; 前文对单元体进行了分折，得到了单元刚度方程 ，但要解决问题，还必须进一步建立整个计算模型的整体刚度方程。完成这一步的关键，在于怎样将单元的刚度矩阵和节点荷载列阵，分别“组装”成整体刚度矩阵和整体节点荷载列阵。 这里通过研究任意节点的平衡来建立整体刚度矩阵，该方法不但比较直观、易懂，而且对怎样编写计算机程序是很有帮助的。; 整体节点载荷列阵：由各节点所受等效节点力按节点号码以从小到大的顺序??列组成的列阵。等效节点力是由集中力、表面力和体积力共同移置构成的，其中集中力包括直接作用在弹性体上的外力和边界约束力，如支座反力。;式中： ，; 对每个单元都可以写出相应的单元刚度方程 ，即单元节点平衡方程。例如，对①号单元，有 ; 同理，对于②单元，有 ;对于3单元，有;对于④单元，有; 对于任意一个节点，可能承受两种力的作用，一种是其它单元给予该节点的反作用力；另一种是作用在节点上的等效节点力。对整体而言，前者属于内力，后者属于外力，每个节点在两种力的作用下处于平衡。 将各单元刚度方程左边相加，即将各节点所受力相加，由于对于整体而言，单元给予节点的反作用力属于内力，在相加过程中相互抵消，所以各节点所受力相加的结果只有外力，即等效节点力，从而得到整体节点荷载列阵，如下; 将各单元刚度方程右边相加，从而得到整体刚度矩阵，如下; 通过以上分析得，整体节点载荷与整体节点位移之间的关系式，即结构整体有限元方程，如下; 整体刚度矩阵组装的基本步骤： 1）将单元刚度矩阵中的每个子块放到在整体刚度矩阵中的对应位置上，得到单元的扩大刚度矩阵。注意对于单元刚度矩阵是按照局部编码排列的，即对应单元刚度矩阵中的i、j、m；对于整体刚度矩阵是按照整体编码排列的，即按节点号码以从小到大的顺序排列。在组装过程中，必须知道单元节点的局部编码与该节点在整体结构中的整体编码之间的关系，才能得到单元刚度矩阵中的每个子块在整体刚度矩阵中的位置。将单元刚度矩阵中的每个子块按总体编码顺序重新排列后，可以得到单元的扩大矩阵。例如在图中，单元②的局部编码为i、j、m，对应整体编码为1、3、4，然后将单元②刚度矩阵中的每个子块按总体编码顺序重新排列后，可以得到单元的扩大矩阵。注意有些书籍中将局部编码表示为1、2、3或1，2，3等； 2）将全部单元的扩大矩阵相加得到整体刚度矩阵。;通过以上组装过程可以得到组装整体刚度矩阵的一般规则： 1）结构中的等效节点力是相关单元结点力的叠加，整体刚度矩阵的子矩阵是相关单元的单元刚度矩阵子矩阵的集成；; 3）当 中 时，若rs边是组合体的内边，则 就是共用该边的两相邻单元刚度矩阵中的子矩阵 的相加。如13边为单元①和②的共用边，则 ; 上述组装基本步骤和规则具有普遍意义 对于不同类型、不同形式的单元，只是相应节点的子矩阵的阶数（节点自由度×节点自由度）可能不同，至于组装整体刚度矩阵的规律仍是相同的。正是因为有了这种组装规律，使得有限元法能够很方便地应用电子计算机进行计算。;;单元编号;同上例类似的分析，得;;，由图形可知，25边为单元②和③的共用边，则;;;cumtfem@126.com;;4. 是奇异矩阵，在排除刚体位移后，它是正定阵 ;1. 是对称矩阵;3. 是稀疏矩阵，非零元素呈带状分布;5.2 整体刚度矩阵的特点;; 显然，带状刚度矩阵的带宽取决于单元网格中相邻节点号码的最大差值D。把半个斜带形区域中各行所具有的非零元素的最大个数叫做刚度矩阵的半带宽（包括主对角元），用B表示，如下 B=2(D+1) 通常的有限元程序，一般都利用刚度矩阵的对称和稀疏带状的特点，在计算求解中，只存储上半带的元素，即所谓的半带存储。因此，在划分完有限元网格进行节点编号时，要采用合理的编码方式，使同一单元中相邻两节点的号码差尽可能小，以便节省存储空间、提高计算效率。; 对于同样的有限元单元网格，按照图（a）的结点编码，最大的半带宽为 按照图（b）的结点编码，最大的半带宽为B 按照图（c）的结点编码，最大的半带宽为B ;