5 平面分析-整体刚度矩阵

《5 平面分析-整体刚度矩阵》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、2021/6/13平面问题有限元分析-总刚15平面问题有限元分析整体刚度矩阵曹国华5.1整体刚度矩阵5.2整体刚度矩阵的特点5.3边界条件5.4计算结果整理5.5收敛准则前文对单元体进行了分折，得到了单元刚度方程，但要解决问题，还必须进一步建立整个计算模型的整体刚度方程。完成这一步的关键，在于怎样将单元的刚度矩阵和节点荷载列阵，分别“组装”成整体刚度矩阵和整体节点荷载列阵。这里通过研究任意节点的平衡来建立整体刚度矩阵，该方法不但比较直观、易懂，而且对怎样编写计算机程序是很有帮助的。2021/6/132平面问题有限元分析-总刚5.1整体刚度矩

2、阵整体节点载荷列阵：由各节点所受等效节点力按节点号码以从小到大的顺序排列组成的列阵。等效节点力是由集中力、表面力和体积力共同移置构成的，其中集中力包括直接作用在弹性体上的外力和边界约束力，如支座反力。为了研究整体刚度矩阵的组装过程，先引入两个概念。整体节点位移列阵：由各节点位移按节点号码以从小到大的顺序排列组成的列阵。2021/6/133平面问题有限元分析-总刚5.1整体刚度矩阵式中：，不失一般性，仅考虑计算模型中有4个单元，如图所示。四个单元的整体节点位移列阵为2021/6/134平面问题有限元分析-总刚5.1整体刚度矩阵对每个单元都可以

3、写出相应的单元刚度方程，即单元节点平衡方程。例如，对①号单元，有式中：——①号单元中第i（i=1,2,3）节点所受力。为了便于组装整体刚度矩阵，将上式以整体节点位移表示，即——①号单元的扩大刚度矩阵或称为单元贡献矩阵。2021/6/135平面问题有限元分析-总刚5.1整体刚度矩阵同理，对于②单元，有式中：——②号单元中第i（i=1,3,4）节点所受力；——②号单元的扩大刚度矩阵。2021/6/136平面问题有限元分析-总刚5.1整体刚度矩阵对于3单元，有式中：——3号单元中第i（i=3,4,5）节点所受力；——3号单元的扩大刚度矩阵。202

4、1/6/137平面问题有限元分析-总刚5.1整体刚度矩阵对于④单元，有式中：——④号单元中第i（i=3,4,5）节点所受力；——④号单元的扩大刚度矩阵。2021/6/138平面问题有限元分析-总刚5.1整体刚度矩阵对于任意一个节点，可能承受两种力的作用，一种是其它单元给予该节点的反作用力；另一种是作用在节点上的等效节点力。对整体而言，前者属于内力，后者属于外力，每个节点在两种力的作用下处于平衡。将各单元刚度方程左边相加，即将各节点所受力相加，由于对于整体而言，单元给予节点的反作用力属于内力，在相加过程中相互抵消，所以各节点所受力相加的结果只

5、有外力，即等效节点力，从而得到整体节点荷载列阵，如下2021/6/139平面问题有限元分析-总刚5.1整体刚度矩阵将各单元刚度方程右边相加，从而得到整体刚度矩阵，如下2021/6/1310平面问题有限元分析-总刚5.1整体刚度矩阵通过以上分析得，整体节点载荷与整体节点位移之间的关系式，即结构整体有限元方程，如下式中：——整体刚度矩阵。2021/6/1311平面问题有限元分析-总刚5.1整体刚度矩阵整体刚度矩阵组装的基本步骤：1）将单元刚度矩阵中的每个子块放到在整体刚度矩阵中的对应位置上，得到单元的扩大刚度矩阵。注意对于单元刚度矩阵是按照局部

6、编码排列的，即对应单元刚度矩阵中的i、j、m；对于整体刚度矩阵是按照整体编码排列的，即按节点号码以从小到大的顺序排列。在组装过程中，必须知道单元节点的局部编码与该节点在整体结构中的整体编码之间的关系，才能得到单元刚度矩阵中的每个子块在整体刚度矩阵中的位置。将单元刚度矩阵中的每个子块按总体编码顺序重新排列后，可以得到单元的扩大矩阵。例如在图中，单元②的局部编码为i、j、m，对应整体编码为1、3、4，然后将单元②刚度矩阵中的每个子块按总体编码顺序重新排列后，可以得到单元的扩大矩阵。注意有些书籍中将局部编码表示为1、2、3或1，2，3等；2）将全

7、部单元的扩大矩阵相加得到整体刚度矩阵。2021/6/1312平面问题有限元分析-总刚5.1整体刚度矩阵通过以上组装过程可以得到组装整体刚度矩阵的一般规则：1）结构中的等效节点力是相关单元结点力的叠加，整体刚度矩阵的子矩阵是相关单元的单元刚度矩阵子矩阵的集成；2）当整体刚度矩阵中的子矩阵中r=s时，该节点（节点r或s）被哪几个单元所共有，则就是这几个单元的刚度矩阵中的子矩阵的相加。如应该是单元①－④中对应子矩阵的集成，即2021/6/1313平面问题有限元分析-总刚5.1整体刚度矩阵3）当中时，若rs边是组合体的内边，则就是共用该边的两相邻单

8、元刚度矩阵中的子矩阵的相加。如13边为单元①和②的共用边，则4）当中r和s不同属于任何单元时，则=0。如节点r＝1和s＝5不同属于任何单元，此时=0。2021/6/1314平面问