第3章实用计算方法.pdf

第3章实用计算方法 在本科学习基础上补充和加深 第3章实用计算方法 目录 §3.1 多自由度体系自由度的确定 §3.2 瑞利-里兹法 §3.3 子空间迭代法 §3.1 多自由度体系自由度的确定 在本科内容复习中已经指出体系的自由度取决于确定 质量空间（或平面）位置的独立坐标数。对于图示20层 框架结构，总计有640个结点，静力计算时每结点6个自 由度，因此共有3840个未知位移。 当然，做动力计算时可取与静力分 析相同的模型(例如用单元一致质量矩 阵建立质量矩阵) ，但这是完全没有必 要的。 由于楼板在自身平面内刚度“无穷 大” ，除非整体属于细长的建筑，一般 不考虑轴向变形影响，此外忽略转动 惯量的惯性作用 ，这样每层只有3个自 由度，总计60个动力自由度即可。 §3.1 多自由度体系自由度的确定 在采用集中质量法的情况下，一般都不考虑质量的转 动惯性作用(也即取转动惯量为零) ，这样静力分析的全 部未知位移就可以分成两部分：有和无惯性作用，如果 将位移编号先编有惯性作用后编无惯性作用，则无阻尼 运动方程可写作   t  M 0 u K K u P         s s ss s s s           0 0 u K K u P t    s            式中下标s表示有惯性作用，表示无惯性作用。由此可 进行如下“静力凝聚” ，由矩阵方程第二式可得 u K 1 P t K u       s s  再将此式代回矩阵方程第一式，得只有未知量s的运动方 程  1 1 M u  K K K K u P t K P t s s  ss s  s  s s       这就是缩减自由度后的运动方程。 §3.2 瑞利-里兹法 虽然用静力凝聚可以大量縮减自由度，但缩减自由度 数目和寻求低阶固有频率及振型近似解的最通用技术是 瑞利-里兹法，下面分六部分来介绍它。 3.2.1 瑞利-里兹法 3.2.2 “最佳”近似 3.2.3 近似频率和振型 3.2.4 近似振型的正交性 3.2.5 里兹向量的选择 3.2.6 力相关的里兹向量 3.2.1 瑞利-里兹法 像弹性力学中介绍的类似，设系统位移为 J u(t)= ψ z t Ψz t