有限元基础理论教程.pptVIP

第1页，共24页，星期日，2025年，2月5日先考虑结点n有水平方向位移约束，与n结点水平方向对应的平衡方程为：在[K]矩阵中，第2n-1行的对角线元素改为1，该行中全部非对角线元素改为0；在{P}中，第2n-1个元素改为0。为了保持[K]矩阵的对称性，将第2n-1列的全部非对角元素也改为0。应该换成下面的方程：第2页，共24页，星期日，2025年，2月5日如果结点n在垂直方向有位移约束，则（2-29）中的第2n个方程修改为，在[K]矩阵中，第2n行的对角线元素改为1，该行中全部非对角线元素改为0；在{P}中，第2n个元素改为0。为了保持[K]矩阵的对称性，将第2n列的全部非对角元素也改为0。第3页，共24页，星期日，2025年，2月5日例2.9、结构的位移约束条件如图所示，结构平衡的方程组如下。修改整体刚度矩阵。第4页，共24页，星期日，2025年，2月5日根据结点4的位移约束，修改该整体刚度矩阵的系数。第5页，共24页，星期日，2025年，2月5日根据结点1和结点6的位移约束条件继续修改整体刚度矩阵可以得到以下的形式，第6页，共24页，星期日，2025年，2月5日如果结点n处存在一个已知非零的水平方向位移，这时的约束条件为，在[K]矩阵中，第2n-1行的对角线元素乘上一个大数A，将方程修改为，第7页，共24页，星期日，2025年，2月5日处理约束条件的意义在于，强迫单元内部的分片近似位移场满足整体结构的位移边界条件，得到整个结构的近似解。在定义位移约束时，要消除结构的刚体位移。第8页，共24页，星期日，2025年，2月5日1）用相关单元的单元刚度矩阵计算结点力，再由力的平衡关系得到约束反力。与位移约束对应的约束反力如何计算？2）给矩阵相应的对角元素加上一个大数，将载荷列阵的对应元素置为零。相当于用一个刚度很大的弹簧代替位移约束。n结点在x方向位移为零，用刚度系数为A的弹簧代替，由于刚度很大，un是一个很小的值。方程修改为，第9页，共24页，星期日，2025年，2月5日2.8整体刚度矩阵的特点与存储方法 整体刚度矩阵具有以下几个显著的特点：对称性，稀疏性，非零系数带形分布。1）对称性 由单元刚度矩阵的对称性和整体刚度矩阵的集成规则，可知整体刚度矩阵必为对称矩阵。利用对称性，只保存整体矩阵上三角部分的系数即可。2）稀疏性 单元刚度矩阵的多数元素为零，非零元素的个数只占较小的部分。第10页，共24页，星期日，2025年，2月5日3）非零元素带形分布整体刚度矩阵的非零元素分布在以对角线为中心的带形区域内，这种矩阵称为带形矩阵。在包括对角线元素的半个带形区域内，每行具有的元素个数叫做半带宽。最大半带宽用d表示，例2.9所示结构的最大半带宽为，第11页，共24页，星期日，2025年，2月5日二维等带宽存储 设整体刚度矩阵[K]为一个n行、n列的矩阵，最大半带宽为d。利用带形矩阵的特点和对称性，只需要保存以d为固定带宽的上半带的元素，称为二维等带宽存储。 整体刚度矩阵[K]每行中的上半带元素取出，保存在另一个矩阵[K*]的对应行中，得到一个n行、d列矩阵[K*]。 把元素在[K]矩阵中的行、列编码记为r、s，在矩阵[K*]中的行、列编码记为r*、s*，对应关系如下： r*=r s*=s-r+1第12页，共24页，星期日，2025年，2月5日[K]矩阵中的对角线元素保存在新矩阵中的第1列中，[K]矩阵中的r行元素仍然保存在新矩阵的r行中，[K]矩阵中的s列元素则按照新的列编码保存在新矩阵的不同列中。仍然会保存一些零元素，但是元素寻址很方便。整体刚度矩阵[K]等带宽矩阵[K*]第13页，共24页，星期日，2025年，2月5日二维等带宽存储所存储的元素数量取决于最大半带宽d的值，d的值则由单元结点的编码方式决定。相同的有限单元网格按照图2.13(a)的结点编码，最大的半带宽为14；按照图2.13(b)的结点编码，最大的半带宽为18。图2.13(a)图2.13(b)第14页，共24页，星期日，2025年，2月5日一维变带宽存储按行的一维变带宽存储，按照每一行的半带宽把半带宽内的元素保存到一维数组中。使用辅助数组定义对角元素在一维数组中的位置。需要存储的矩阵元素最少。解方程时的地址计算比较复杂，会带来一些额外的计算量。第15页，共24页，星期日，2025年，2月5日2.9线性方程组解法线性方程组的解法：直接解法 包括高斯消去法、等带宽高斯消去法、三角分解法，以及适用于大型方程组求解的分块算法和波前法等。迭代解法 雅可比迭代法、高斯-赛德尔迭代法、超松弛迭代法和共轭梯度法等。第16页，共24页，星期