结构力学第8章.pptVIP

烟台大学 第四章 静定结构的位移计算 烟台大学 第8章 矩阵位移法 返回 开篇 帮助 一、基本概念 结构矩阵分析是采用矩阵方法分析结构力学问题的一种方法。与传统的力法、位移法相对应，在结构矩阵分析中也有矩阵力法和矩阵位移法，或柔度法与刚度法。矩阵位移法易于实现计算过程程序化而被广泛应用。 矩阵位移法是结构力学中的位移法加上矩阵方法。矩阵位移法的基本未知量也是结点位移——独立的线位移和转角。但由于有时考虑杆件的轴向变形，且把杆件铰结端的转角也作为基本未知量，因此，基本未知量数目比传统位移法的基本未知量多一些。 矩阵位移法的基本思路是: (1) 先把结构离散成单元，进行单元分析,建立单元杆端力与杆端位移之间的关系； (2)在单元分析的基础上，考虑结构的几何条件和平衡条件，将这些离散单元组合成原来的结构，进行整体分析，建立结构的结点力与结点位移之间的关系，即结构的总刚度方程，进而求解结构的结点位移和单元杆端力。 在从单元分析到整体分析的计算过程中，全部采用矩阵运算。 集成总刚度矩阵最常用的方法是直接刚度法，即由单元刚度矩阵直接集成结构刚度矩阵，又可分为后处理法和先处理法。 1. 后处理法 (1) 集成。对所有单元不做边界条件处理，均采用自由式的单元刚度矩阵，按单元的结点编号将单元刚度矩阵分为四个子块（阶数相同），逐块地将结点所对应的子块在结构的原始刚度矩阵中对号入座，形成结构的原始刚度矩阵。由于结点位移分量中包括了非自由结点的已知位移，原始刚度矩阵为奇异的，需进行边界条件处理，才能求解自由结点位移。由于原始刚度矩阵的阶数较高，所以后处理法的主要缺点是占用较多的计算机内存。 二、总刚度矩阵的集成及约束处理 对于每个结点位移分量数相同的结构，原始刚度矩阵的阶数为结构的总结点数乘以结点位移分量的数目，例如，每个结点位移分量数为3的平面刚架，结构原始刚度矩阵的阶数为3n×3n 。 对于刚性支座，用划行划列法处理刚性支座，即直接划去原始刚度方程中与零位移对应的行和列。这样做有时要改变原方程的排列顺序,会给编程带来麻烦。为了不改变原方程的排列顺序，同时又要引入边界条件,采用“主一副零”法。 （2）边界条件处理 设结点位移向量中第r个位移等于零, 即?r=0 ，则在结构的原始刚度矩阵k中的第r行第r列中主对角元素krr改为1其余元素改为零。同时将结点结点荷载列向量P中的第r个分量也改为零。 即 对于支座位移等于给定值时，采用“乘大数法”。设结点位移向量中第r个位移等于d0，在矩阵K与向量P中，主对角元素krr 改为G?krr，将Pr改为d0?G?krr，其中G为一大数通常取108～1010 。 ， 单元定位向量：按单元连接结点编号顺序由结点未知位移编号组成的向量。 2. 先处理法 (1) 集成。将单元刚度矩阵先按边界条件进行处理,然后按照单元连接结点的总位移编号将单元刚度矩阵的元素在结构的刚度矩阵中对号入座，形成总刚后即可进行求解。上述过程可通过引入定位向量来实现。在单元定位向量中考虑边界条件，凡给定的结点位移分量，其位移总码均编为零，与总码编为零相应的行、列元素在集成总刚时被屏弃在外。 （2）边界条件处理。对于刚性支座，其位移总码均编为零。对于支座位移等于给定值时，通常也将其位移总码均编为零，将支座结点位移的影响转换成单元非结点荷载,即，将支座结点位移转换成与该支座结点位移连接的各单元在单元坐标系中的杆端位移，求出由此给定的杆端位移产生的单元固端力，然后转换成等效结点荷载。 通常用主对角元素叠加法处理弹性支座。如果结构的第j个自由度是弹性约束，那么，把弹性支座的刚度系数叠加到原始刚度矩阵主对角线的第j个元素上即可得到经约束处理后的总刚度方程。 3. 弹性支座的处理 总刚度方程为整体结构的结点荷载与结点位移之间的关系式，是结构应满足的平衡条件。无论何种结构，其总刚度方程都具有统一的形式: 4. 总刚度方程和总刚度矩阵的性质与特点 K?=P 式中K为总刚度矩阵，?为结构的结点位移列向量，P为结点力列向量。 总刚度矩阵K反应了整个结构的刚度，是描述结点力与结点位移之间关系的系数矩阵。其矩阵的性质与特点： ( 1 )元素kij的物理意义为：当△j=1而其他位移分量为零时产生在△i方向的杆端力。 （2）主子块Kii是由结点i的相关单元中与结点i相应的主子块叠加而得。 （3）当i、j为相关结点时，副子块Kij就等于连接ij的杆单元中相应的子块；若i、j不相关，则Kij为零子块。 （4）总刚度矩阵为对称矩阵。 （5）总刚度矩阵为稀疏带状矩阵