一维变带宽有限元刚度矩阵C++实现(系列二)
发布时间
阅读量:
阅读量
在上一篇中,我们获得了辅助数组pDiag中存储的半带宽值。在节点自由度重新编号后,在节点自由度重新编号后,在节点自由度重新编号后,在节点自由度重新编号后,在节点自由度重新编号后,在节点自由度重新编号后,在节点自由度重新编号后,在节点自由度重新编号后,在节点自由度重新编号后,在节点自由度重新编号后,在节点自由度重新编号后,在节点自由度重新编号后的总刚度矩阵[K]呈现为:

此时,所有的非零元素都集中在带内。
一维变带宽存储的思路为:
用于将总刚度矩阵[K]其下三角区域中的每一行数据从各自带内起始位置依次排列形成一个序列,并将其存储在一维数组pGK[L]中。其中L表示一维数组中元素的数量,即总刚度矩阵下三角区域的所有带内元素数量。
但是为了明确确定pGK中的元素在[K]中的行列编号,则需将[K]中各对角线上的元素被存储到一维数组pGK中的特定位置(称为对角元地址),并将这些值存入另一辅助数组pDiag[n]中。(其中n表示总刚矩阵[K]的阶数)举例说明这种存储方法:设有某个系数矩阵:
全部评论 (0)
还没有任何评论哟~
