哈工大研究生课程高等结构动力学第五章概论.pptVIP

第五章 结构动力学中常用的数值解法 响应求解方法的选择取决的因素有：载荷、结构、精度要求、非线性影响程度、方法的稳定性等。 作业：5-1（2） 作业：5-1（1）、（3）- 5 ;5-4 * * §5.1概述 数值分析技术为结构的动态分析提供了有力的保障，为工程结构在各种复杂的动力学环境下的模拟和仿真提供了有效工具。 工程结构的动态分析主要包括两个方面：结构的动态特性分析和结构动态响应分析 标准特征值问题 1 雅可比方法（Jacobi）、 2.Rayleigh-Ritz 3.子空间迭代法 4. 行列式搜索法 行列式搜索法是求解大型特征值问题的另一种方法。它的特点是综合运用多项式加速割线迭代，移轴向量逆迭代，Sturm序列的性质以及Gram－Schmidt正交化过程，直接计算所需要的任意特征对，通常是计算最小的部分特征值及相应的特征向量。 因此，它是一种计算部分特征对的特殊求解方法。此方法具有计算速度快，精度高，灵活等优点。 响应数值分析： 1.中心差分法 2.Wilson-θ法 3.Newmark法 ? 综合各方面的因素，比较、权衡，才能判定所应采取的方法；有时为了互相验证，也可以同时采取两种以上的方法来处理动响应分析 ? 对结构过于复杂的情况，宜采用直接积分法，结构较简单的情况可采用模态迭加法。 ? 对精度要求较低的初步设计阶段，可采用取少数模态的模态迭加法。对精度要求较高的最后设计阶段，宜采用直接积分法 §5.2 求解系统固有频率主振型的近似解法 1.邓克利法：是邓克利首先通过实验方法建立起来的一个计 算公式，后来才得到完整的数学证明。 设质量矩阵 ，柔度矩阵为 则有 1894年邓克利：提出一种近似计算多圆盘轴横向振动基频的实用方法（偏小） 设系统作j阶主振动,则有： 代入得特征方程： 有 假设质量矩阵为对角阵，展开得： 根据多项式的根与系数之间的关系 的n个根 ， 之和为 由于二阶频率往往比基频高得多 忽略 设 如例题 2.雅可比（Jacobi）法求特征方程 即可断定[D]的n个对角元素就是[A]的n个特征值，而[S]的第i列就是[D]中第i个对角元素所对应的特征向量， [S]为坐标变换矩阵。 在[A]中非对角线元素中选取一个绝对值最大的元素，设为 ，利用平面旋转矩阵 对 进行正交变换： 其中 用雅可比法求n阶对称矩阵[A]的特征值和特征向量的步骤 设 为单位矩阵 ②在A中选取非对角线元素中绝对值最大的元素 ④Vb实现 ③ 3.瑞利（Rayleigh）法和里兹法 已知系统的刚度[K]，质量[M]，并设定系统的j阶主振型为 对于作简谐运动的多自由度系统，其动能T与势能V 系统作j阶主振动时 速度及加速度 里兹法：是瑞利法的改进 将瑞利法使用的单个假设模态改进为若干个独立的假设模态 的线性组合 瑞利商在真实模态处取驻值 （参见 刘延柱振动力学 107页和谢官模振动力学244页） 得 问题又归结为矩阵的本征值问题，但与原系统的本征值比，矩阵的阶数r小于原系统的阶数n. 4.矩阵迭代法 标准化 整理得： 对于正定系统： 可写为 迭代 假设一个初值 解： 取初值： 振型值趋于稳定。 迭代终止： 几种数值算法的比较： 根据瑞利商的性质，原则上可用瑞利商计算任意阶固有频率，但由于高阶主振型很难合理假设，所以瑞利商一般用于求基频，瑞利商求的基频是真实值的上限，这是因为假设的一阶主振型与真实振型的偏差，相当于对系统附加了某些约束，从而提高了系统的刚度，使基频有所提高。 矩阵迭代法只能求前几阶频率，求较高阶的频率和模态需要把前几阶的模态剔除才能收敛。 邓柯莱公式建立的是质量矩阵为对角阵的情况，而且第二阶及第二阶以上的固有频率通常远大于基频，邓柯莱公式计算出的基频显然是精确值的下限。 5.子空间迭代法 子空间迭代法实质就是对一组试验向量反复地使用里 兹法和矩阵迭代法，它将矩阵迭代法每次迭代一个假设模态，发展为同时迭代系统的前r阶假设模态，因而提高了计算效率，迭代过程中各阶假设模态的正交性由里兹法保证。 设系统的前r阶模态 构成全部n阶模态所张成的线性空间的一个子空间，任 选r个独立的向量 作为子空间的假设模态,组成 阶矩阵 各个假设模态 总能表示成真实模态的线性组合。 ②式 代入②式两边乘以 阵 如此迭代k次后 第二项包含的高于r阶模态的成分比第一个求和式更快趋近 于零。 子空间迭代的步骤 假设初始模态矩阵 将上式左乘矩阵[?D] 利用 各列的线性组合表示子空间基的一次近似 以 为假设模态，进行里兹法计算，以确定系数矩阵 先作出以下r阶方阵 瑞利商取驻值 解此本征值，得到 个本征值和 个本征向量 代入式 中则子空间的一次近似 完全确定，构成的各个模态满足正交性条件，至此完成第一次迭代， 将 代替 进