结构力学第8章307999658
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
（3）当i、j为相关结点时，副子块Kij就等于连接ij的杆 单元中相应的子块；若i、j不相关，则Kij为零子块。
（4）总刚度矩阵为对称矩阵。
（5）总刚度矩阵为稀疏带状矩阵。愈是大型结构，带 状分布规律就愈明显。
（6）总刚度矩阵主对角元素都大于零。通常是主对角 元素占优势的矩阵，因源自文库，线形方程组的解有较好的稳定 性。
3. 弹性支座的处理 通常用主对角元素叠加法处理弹性支座。如果结构的 第j个自由度是弹性约束，那么，把弹性支座的刚度系数 叠加到原始刚度矩阵主对角线的第j个元素上即可得到经 约束处理后的总刚度方程。
4. 总刚度方程和总刚度矩阵的性质与特点 总刚度方程为整体结构的结点荷载与结点位移之间的 关系式，是结构应满足的平衡条件。无论何种结构，其总 刚度方程都具有统一的形式:
K=P
式中K为总刚度矩阵，为结构的结点位移列向量，P
为结点力列向量。
总刚度矩阵K反应了整个结构的刚度，是描述结点 力与结点位移之间关系的系数矩阵。其矩阵的性质与 特点：
( 1 )元素kij的物理意义为：当△j=1而其他位移分量为零 时产生在△i方向的杆端力。
（2）主子块Kii是由结点i的相关单元中与结点i相应的主 子块叠加而得。
一、基本概念
结构矩阵分析是采用矩阵方法分析结构力学问题的一种 方法。与传统的力法、位移法相对应，在结构矩阵分析中 也有矩阵力法和矩阵位移法，或柔度法与刚度法。矩阵位 移法易于实现计算过程程序化而被广泛应用。
矩阵位移法是结构力学中的位移法加上矩阵方法。矩阵 位移法的基本未知量也是结点位移——独立的线位移和转 角。但由于有时考虑杆件的轴向变形，且把杆件铰结端的 转角也作为基本未知量，因此，基本未知量数目比传统位 移法的基本未知量多一些。
设结点位移向量中第r个位移等于零, 即r=0 ，则在结
构的原始刚度矩阵k中的第r行第r列中主对角元素krr改为1 其余元素改为零。同时将结点结点荷载列向量P中的第r个
分量也改为零。 即
krr 1 krs ksr 0 Pr 0
(s r)
对于支座位移等于给定值时，采用“乘大数法”。设
结点位移向量，中第r个位移等于d0，在矩阵K与向量P中， 主对角元素krr 改为Gkrr，将Pr改为d0Gkrr，其中G为一大 数通常取108～1010 。
单元定位向量：按单元连接结点编号顺序由结点未知 位移编号组成的向量。
（2）边界条件处理。对于刚性支座，其位移总码均编 为零。对于支座位移等于给定值时，通常也将其位移总码 均编为零，将支座结点位移的影响转换成单元非结点荷载, 即，将支座结点位移转换成与该支座结点位移连接的各单 元在单元坐标系中的杆端位移，求出由此给定的杆端位移 产生的单元固端力，然后转换成等效结点荷载。
36
9 12 15 18 21 24 27
deMA M X I1 N
如果结构内部存在组合结点，并采用先处理法，则不 能用上述公式计算总刚度矩阵的最大半带宽，而应按照单 元编，利用单元定位向量求出总刚度矩阵的最大半带宽。 设用MAX表示单元(e)定位向量中的最大分量，MIN表示 单元(e)定位向量中的最小分量，则
矩阵位移法的基本思路是:
(1) 先把结构离散成单元，进行单元分析,建立单元杆 端力与杆端位移之间的关系；
(2)在单元分析的基础上，考虑结构的几何条件和平衡条 件，将这些离散单元组合成原来的结构，进行整体分析， 建立结构的结点力与结点位移之间的关系，即结构的总刚 度方程，进而求解结构的结点位移和单元杆端力。
例如图示刚架，按图a 编码，d=3×(9+1)=30 ，而按 b 图编码，d=3×(3+1)=12 。
(a) 1
2
3
4
56 7 89
10 11 12 13 14 15 16 17 18
19 20 21 22 23 24 25 26 27
(b) 1
4
7
10 13 16 19 22 25
25
8 11 14 17 20 23 26
对于每个结点位移分量数相同的结构，原始刚度矩阵 的阶数为结构的总结点数乘以结点位移分量的数目，例 如，每个结点位移分量数为3的平面刚架，结构原始刚 度矩阵的阶数为3n×3n 。
（2）边界条件处理
对于刚性支座，用划行划列法处理刚性支座，即直接 划去原始刚度方程中与零位移对应的行和列。这样做有 时要改变原方程的排列顺序,会给编程带来麻烦。为了不 改变原方程的排列顺序，同时又要引入边界条件,采用 “主一副零”法。
5. 总刚度矩阵的最大半带宽
总刚度矩阵的上三角部分，从某行的主对角元素到该行 最末一个非零元素所具有的元素的个数称为该行的半带宽。 各行半带宽的最大值称为总刚度矩阵的最大半带宽。
对应于后处理法，结构内部不存在组合结点时最大半带 宽的计算公式为:d=(b+1)c ，其中b为单元两端结点编码的 最大差; c为结构中一个结点的位移分量数，显然，最大半 带宽与结构的结点编码的顺序有关。通常应使相邻结点编 码的最大差值为最小，即d 值为最小。
在从单元分析到整体分析的计算过程中，全部采用矩阵 运算。
二、总刚度矩阵的集成及约束处理
集成总刚度矩阵最常用的方法是直接刚度法，即由单元 刚度矩阵直接集成结构刚度矩阵，又可分为后处理法和先 处理法。
1. 后处理法
(1) 集成。对所有单元不做边界条件处理，均采用自由 式的单元刚度矩阵，按单元的结点编号将单元刚度矩阵分 为四个子块（阶数相同），逐块地将结点所对应的子块在 结构的原始刚度矩阵中对号入座，形成结构的原始刚度矩 阵。由于结点位移分量中包括了非自由结点的已知位移， 原始刚度矩阵为奇异的，需进行边界条件处理，才能求解 自由结点位移。由于原始刚度矩阵的阶数较高，所以后处 理法的主要缺点是占用较多的计算机内存。
2. 先处理法
(1) 集成。将单元刚度矩阵先按边界条件进行处理,然后 按照单元连接结点的总位移编号将单元刚度矩阵的元素在 结构的刚度矩阵中对号入座，形成总刚后即可进行求解。 上述过程可通过引入定位向量来实现。在单元定位向量中 考虑边界条件，凡给定的结点位移分量，其位移总码均编 为零，与总码编为零相应的行、列元素在集成总刚时被屏 弃在外。
（4）总刚度矩阵为对称矩阵。
（5）总刚度矩阵为稀疏带状矩阵。愈是大型结构，带 状分布规律就愈明显。
（6）总刚度矩阵主对角元素都大于零。通常是主对角 元素占优势的矩阵，因源自文库，线形方程组的解有较好的稳定 性。
3. 弹性支座的处理 通常用主对角元素叠加法处理弹性支座。如果结构的 第j个自由度是弹性约束，那么，把弹性支座的刚度系数 叠加到原始刚度矩阵主对角线的第j个元素上即可得到经 约束处理后的总刚度方程。
4. 总刚度方程和总刚度矩阵的性质与特点 总刚度方程为整体结构的结点荷载与结点位移之间的 关系式，是结构应满足的平衡条件。无论何种结构，其总 刚度方程都具有统一的形式:
K=P
式中K为总刚度矩阵，为结构的结点位移列向量，P
为结点力列向量。
总刚度矩阵K反应了整个结构的刚度，是描述结点 力与结点位移之间关系的系数矩阵。其矩阵的性质与 特点：
( 1 )元素kij的物理意义为：当△j=1而其他位移分量为零 时产生在△i方向的杆端力。
（2）主子块Kii是由结点i的相关单元中与结点i相应的主 子块叠加而得。
一、基本概念
结构矩阵分析是采用矩阵方法分析结构力学问题的一种 方法。与传统的力法、位移法相对应，在结构矩阵分析中 也有矩阵力法和矩阵位移法，或柔度法与刚度法。矩阵位 移法易于实现计算过程程序化而被广泛应用。
矩阵位移法是结构力学中的位移法加上矩阵方法。矩阵 位移法的基本未知量也是结点位移——独立的线位移和转 角。但由于有时考虑杆件的轴向变形，且把杆件铰结端的 转角也作为基本未知量，因此，基本未知量数目比传统位 移法的基本未知量多一些。
设结点位移向量中第r个位移等于零, 即r=0 ，则在结
构的原始刚度矩阵k中的第r行第r列中主对角元素krr改为1 其余元素改为零。同时将结点结点荷载列向量P中的第r个
分量也改为零。 即
krr 1 krs ksr 0 Pr 0
(s r)
对于支座位移等于给定值时，采用“乘大数法”。设
结点位移向量，中第r个位移等于d0，在矩阵K与向量P中， 主对角元素krr 改为Gkrr，将Pr改为d0Gkrr，其中G为一大 数通常取108～1010 。
单元定位向量：按单元连接结点编号顺序由结点未知 位移编号组成的向量。
（2）边界条件处理。对于刚性支座，其位移总码均编 为零。对于支座位移等于给定值时，通常也将其位移总码 均编为零，将支座结点位移的影响转换成单元非结点荷载, 即，将支座结点位移转换成与该支座结点位移连接的各单 元在单元坐标系中的杆端位移，求出由此给定的杆端位移 产生的单元固端力，然后转换成等效结点荷载。
36
9 12 15 18 21 24 27
deMA M X I1 N
如果结构内部存在组合结点，并采用先处理法，则不 能用上述公式计算总刚度矩阵的最大半带宽，而应按照单 元编，利用单元定位向量求出总刚度矩阵的最大半带宽。 设用MAX表示单元(e)定位向量中的最大分量，MIN表示 单元(e)定位向量中的最小分量，则
矩阵位移法的基本思路是:
(1) 先把结构离散成单元，进行单元分析,建立单元杆 端力与杆端位移之间的关系；
(2)在单元分析的基础上，考虑结构的几何条件和平衡条 件，将这些离散单元组合成原来的结构，进行整体分析， 建立结构的结点力与结点位移之间的关系，即结构的总刚 度方程，进而求解结构的结点位移和单元杆端力。
例如图示刚架，按图a 编码，d=3×(9+1)=30 ，而按 b 图编码，d=3×(3+1)=12 。
(a) 1
2
3
4
56 7 89
10 11 12 13 14 15 16 17 18
19 20 21 22 23 24 25 26 27
(b) 1
4
7
10 13 16 19 22 25
25
8 11 14 17 20 23 26
对于每个结点位移分量数相同的结构，原始刚度矩阵 的阶数为结构的总结点数乘以结点位移分量的数目，例 如，每个结点位移分量数为3的平面刚架，结构原始刚 度矩阵的阶数为3n×3n 。
（2）边界条件处理
对于刚性支座，用划行划列法处理刚性支座，即直接 划去原始刚度方程中与零位移对应的行和列。这样做有 时要改变原方程的排列顺序,会给编程带来麻烦。为了不 改变原方程的排列顺序，同时又要引入边界条件,采用 “主一副零”法。
5. 总刚度矩阵的最大半带宽
总刚度矩阵的上三角部分，从某行的主对角元素到该行 最末一个非零元素所具有的元素的个数称为该行的半带宽。 各行半带宽的最大值称为总刚度矩阵的最大半带宽。
对应于后处理法，结构内部不存在组合结点时最大半带 宽的计算公式为:d=(b+1)c ，其中b为单元两端结点编码的 最大差; c为结构中一个结点的位移分量数，显然，最大半 带宽与结构的结点编码的顺序有关。通常应使相邻结点编 码的最大差值为最小，即d 值为最小。
在从单元分析到整体分析的计算过程中，全部采用矩阵 运算。
二、总刚度矩阵的集成及约束处理
集成总刚度矩阵最常用的方法是直接刚度法，即由单元 刚度矩阵直接集成结构刚度矩阵，又可分为后处理法和先 处理法。
1. 后处理法
(1) 集成。对所有单元不做边界条件处理，均采用自由 式的单元刚度矩阵，按单元的结点编号将单元刚度矩阵分 为四个子块（阶数相同），逐块地将结点所对应的子块在 结构的原始刚度矩阵中对号入座，形成结构的原始刚度矩 阵。由于结点位移分量中包括了非自由结点的已知位移， 原始刚度矩阵为奇异的，需进行边界条件处理，才能求解 自由结点位移。由于原始刚度矩阵的阶数较高，所以后处 理法的主要缺点是占用较多的计算机内存。
2. 先处理法
(1) 集成。将单元刚度矩阵先按边界条件进行处理,然后 按照单元连接结点的总位移编号将单元刚度矩阵的元素在 结构的刚度矩阵中对号入座，形成总刚后即可进行求解。 上述过程可通过引入定位向量来实现。在单元定位向量中 考虑边界条件，凡给定的结点位移分量，其位移总码均编 为零，与总码编为零相应的行、列元素在集成总刚时被屏 弃在外。