第六章杆系单元选读.ppt

6.4.3 整体刚度矩阵的性质 整体刚度矩阵 中位于主对角线上的子块 ，称为主子块，其余 为副子块。 a. 中主子块 由结点i的各相关单元的主子块扩展之后叠加求得，即 b. 当结点i、 j为单元e的相关结点时， 中副子块 为该单元e相应的副子块，即 。 c. 当结点i、 j为非相关结点时， 中副子块 为零子块，即 。 d. 仅与各单元的几何特性、材料特性，即A、I、l、E等因素有关。 e. 为对称方阵， f. 为奇异矩阵，其逆矩阵不存在，因为建立整体刚度矩阵时没有考虑结构的边界约束条件。 6.4 整体刚度矩阵 g. 为稀疏矩阵，整体刚度矩阵中的非零元素分布区域的宽度与结点编号有关，非零元素分布在以对角线为中心的带状区域内，称为带状分布规律，见图6-10(a)。在包括对角线元素在内的区域中，每行所具有的元素个数叫做把半带宽，以d表示。 最大半带宽等于相邻结点号的最大差值加 1 与结点自由度数的乘积，结点号差越大半带宽也就越大。计算机以半带宽方式存储，见图3-10(b)。半带宽越窄，计算机的存储量就越少，而且可以大幅度减少求解方程所需的运算次数。其效果对大型结构显得尤为突出。 图6-10 整体刚度矩阵存储方法 h. 整体刚度矩阵稀疏阵。 故整体刚度矩阵不能求逆，必须作约束处理方能正确地将结点位移求出，进而求出结构的应力场。 (a) 带状分布规律 (b) 带状存储 6.4 整体刚度矩阵 6.5 约束处理及求解 6.5.1 约束处理的必要性 建立结构原始平衡方程式 时，并未考虑支承条件（约束），也就是说，将原始结构处理成一个自由悬空的、存在刚体位移的几何可变结构。整体刚度矩阵是奇异矩阵，因此，无法求解。可以参照第 2 章的原则，结合实际工程结构引入支承条件，即对结构原始平衡方程式 做约束处理。 约束处理后的方程称为基本平衡方程。 统一记为 3.5.2 约束处理方法 约束处理常用方法有填0置1法和乘大数法。采用这两种方法不会破坏整体刚度矩阵的对称性、稀疏性及带状分布等特性。 下面以图6-11所示刚架结构为例，解释如何进行约束处理。对于下图所示刚架结构 设结点位移列向量为 设结点载荷列向量为 固定支座 (b) 支座强迫位移已知 图6-11 结构约束 6.5 约束处理及求解 * 其原始平衡方程式为 按照每个结点的位移分量将上式展开为 6.5 约束处理及求解 对于如图6-11(a)所示，结构约束（支座）位移全部为零，此时做约束处理时，采用填0置1法比较适宜。 对于如图6-11(b)所示，某约束（支座）位移为给定的强迫值，此时做约束处理时，采用乘大数法比较适宜。 (1) 填0置1法 如右图所示结点1、3处为固定支座，可知 将整体刚度矩阵中与之相对应的主对角元素全部置换成1， 相应行和列上的其它元素均改为0。 同时，所在同一行上的载荷分量替换成0，则有 6.5 约束处理及求解 则 也可简便地采用划行划列的办法。在整体刚度矩阵中将与约束位移为 0 的行和列划掉，包括相关的所在行的位移和载荷向量。 6.5 约束处理及求解 处理后得基本平衡方程 (2) 乘大数法 右图所示刚架，结点1为固定支座，结点3处在方向的约束为已知强迫位移。即 将整体刚度矩阵中与之相对应的主对角元素全部乘以一个大数N，一般取 。同时，将相应同一行上的载荷分量替换成 N 乘以其主对角刚度系数和给定的强迫位移（包括零位移）。 6.5 约束处理及求解 NJUT zhuang7802@163.com * * * * NJUT 南京工业大学 第六章 杆系有限单元法 庄海洋 南京工业大学交通学院 6.1 结构离散与向量表示 工