有限元-第3讲补充-平面问题-整体刚度矩阵.ppt

《有限元-第3讲补充-平面问题-整体刚度矩阵.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、15平面问题有限元分析整体刚度矩阵曹国华整体刚度矩阵整体刚度矩阵的特点边界条件处理计算结果整理前文对单元体进行了分折，得到了单元刚度方程，但要解决问题，还必须进一步建立整个计算模型的整体刚度方程。完成这一步的关键，在于怎样将单元的刚度矩阵和节点荷载列阵，分别“组装”成整体刚度矩阵和整体节点荷载列阵。通过研究任意节点的平衡来建立整体刚度矩阵，该方法不但比较直观、易懂，而且对怎样编写计算机程序是很有帮助。2整体刚度矩阵整体节点载荷列阵：由各节点所受等效节点力按节点号码以从小到大的顺序排列组成的列阵。等效节点力是由集中力、表面力和体积力共同移置构成的，其中集中力包括直接作用在弹性体上的外力

2、和边界约束力，如支座反力。为了研究整体刚度矩阵的组装过程，先引入两个概念。整体节点位移列阵：由各节点位移按节点号码以从小到大的顺序排列组成的列阵。3整体刚度矩阵式中：，不失一般性，仅考虑计算模型中有4个单元，如图所示。四个单元的整体节点位移列阵为4整体刚度矩阵对每个单元都可以写出相应的单元刚度方程，即单元节点平衡方程。例如，对①号单元，有式中：——①号单元中第i（i=1,2,3）节点所受力。为了便于组装整体刚度矩阵，将上式以整体节点位移表示，即——①号单元的扩大刚度矩阵或称为单元贡献矩阵。5整体刚度矩阵同理，对于②单元，有式中：——②号单元中第i（i=1,3,4）节点所受力；——②号

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

4、力，即等效节点力，从而得到整体节点荷载列阵，如下9整体刚度矩阵将各单元刚度方程右边相加，从而得到整体刚度矩阵，如下10整体刚度矩阵通过以上分析得，整体节点载荷与整体节点位移之间的关系式，即结构整体有限元方程，如下式中：——整体刚度矩阵。11整体刚度矩阵整体刚度矩阵组装的基本步骤：1）将单元刚度矩阵中的每个子块放到在整体刚度矩阵中的对应位置上，得到单元的扩大刚度矩阵。注意对于单元刚度矩阵是按照局部编码排列的，即对应单元刚度矩阵中的i、j、m；对于整体刚度矩阵是按照整体编码排列的，即按节点号码以从小到大的顺序排列。在组装过程中，必须知道单元节点的局部编码与该节点在整体结构中的整体编码之间

5、的关系，才能得到单元刚度矩阵中的每个子块在整体刚度矩阵中的位置。将单元刚度矩阵中的每个子块按总体编码顺序重新排列后，可以得到单元的扩大矩阵。例如在图中，单元②的局部编码为i、j、m，对应整体编码为1、3、4，然后将单元②刚度矩阵中的每个子块按总体编码顺序重新排列后，可以得到单元的扩大矩阵。注意有些书籍中将局部编码表示为1、2、3或1，2，3等；2）将全部单元的扩大矩阵相加得到整体刚度矩阵。12整体刚度矩阵通过以上组装过程可以得到组装整体刚度矩阵的一般规则：1）结构中的等效节点力是相关单元结点力的叠加，整体刚度矩阵的子矩阵是相关单元的单元刚度矩阵子矩阵的集成；2）当整体刚度矩阵中的子矩

6、阵中r=s时，该节点（节点r或s）被哪几个单元所共有，则就是这几个单元的刚度矩阵中的子矩阵的相加。如应该是单元①－④中对应子矩阵的集成，即13整体刚度矩阵3）当中时，若rs边是组合体的内边，则就是共用该边的两相邻单元刚度矩阵中的子矩阵的相加。如13边为单元①和②的共用边，则4）当中r和s不同属于任何单元时，则=0。如节点r＝1和s＝5不同属于任何单元，此时=0。14整体刚度矩阵上述组装基本步骤和规则具有普遍意义对于不同类型、不同形式的单元，只是相应节点的子矩阵的阶数（节点自由度×节点自由度）可能不同，至于组装整体刚度矩阵的规律仍是相同的。正是因为有了这种组装规律，使得有限元法能够很方

7、便地应用电子计算机进行计算。16例：如图所示有限元模型，弹性模量为，厚度为，为简化计算取，求整体刚度矩阵。17单元编号①②③④整体编码1、2、32、4、55、3、23、5、6局部编码i、j、mi、j、mi、j、mi、j、m以整体编码表示的单元刚度矩阵子块解：该模型中共有6个节点，4个单元，各单元的信息如表所示。各单元信息18同上例类似的分析，得根据单元刚度矩阵的性质，可知，若3单元5,3,2，则整体刚度矩阵中的各子块是对所有单元相应的子块求和得到的（实际只