多自由度体系近似计算方法.ppt

误差 瑞利商 设 误差 瑞利商 3-4 里兹(Ritz)法 对于复杂工程问题，动力分析需要计算系统的前几阶固有频率及相应的主振型 Ritz 法对 Rayleigh 法进行了修正，以实现计算低阶固有频率与振型的目的 Ritz 法是一种减缩系统自由度的近似计算方法 Ritz 法对系统的近似振型 X 给出更合理的假设 为选取的 k 个线性无关的假设振型 待定常数向量 代入瑞利商 瑞利商成为a的函数 利用瑞利商在真实主振型处取驻值的性质，由极值条件 特征值问题的阶数 k n Ritz 法实际是一种减缩系统自由度数求解固有振动的近似计算方法 Ritz 法的基本思想 利用 k 个线性无关的假设振型 为基底 在 n 维振型空间中构成一个 k 维子空间 确定瑞利商在该子空间的 k 个极值 将所得 k 个极值作为原系统前 k 阶固有频率平方的近似值 n自由度系统的固有频率 系统的前 k 阶主振型 证明 所得近似主振型关于 M 和 K 具有正交性 Ritz 法的一些性质 ? 若假设振型恰好是主振型 Ritz 法求出的 就是系统的前 k 阶固有频率的精确值 * * * * 多自由度体系近似计算方法 3-1 邓柯莱(Dunkerly)法 邓柯莱（Dunkerly法） 迹法 确定系统基频的估算公式 方法特点：简单实用 定义 系统的动力矩阵为 n个自由度系统的特征值问题为 标准特征值问题 若将特征值按降序排列 系统的基频为 标准特征值问题的特征行列式为 动力矩阵的对角线元素 由代数方程理论，多项式根与系数关系的韦达定理 动力矩阵A的迹 若质量矩阵M为对角阵，动力矩阵的迹为 对角线元素 M 对角线元素 1 设弹性系统只保留第 i 个质量 mi 及相应的弹簧δii ,则系统视为单自由度系统的固有频率为 ? 邓柯莱法计算系统的基频为精确解的下限 ? 只有当 时，迹法可给出比较准确的基频估算值 ? 算例表明，梁结构通常具有以上的特点 举例 三自由度梁弯曲的固有频率与主振型 m 2m m 系统的质量矩阵与柔度矩阵 举例 均质等直梁，试估算梁中央附加集中质量M时的基频 M m EJ 均质简支梁的基频 记简支梁的基频为 不计简支梁质量时系统的固有频率为 均质梁中央附加集中质量M时的基频 M=m Dunkerly法 Rayleigh法 精确解 3-2 矩阵迭代法 工程中的振动问题的响应分析中，系统的低阶固有频率及主振型占有重要地位 矩阵迭代法是求解系统低阶固有频率和主振型的一种简单实用的方法 ? 第一阶固有频率及主振型 向量 向量 给定一个初始迭代向量 x1，由展开定理 x1 与 Φ(1) 不正交 所占比重增加 所占比重减少 动力矩阵迭代一次后，扩大了第一阶主振型在迭代向量中的优势 第一阶主振型在迭代向量中的优势继续扩大 随着迭代次数的增加，第一阶主振型的优势越来越大。当迭代次数充分大时，可近似地得到 迭代后的新向量与原向量个对应元素间仅相差一常数倍 λ1 ? 迭代过程中应对迭代向量作归一化处理 ? 迭代过程收敛速度取决于比值 趋于零的速度 ? 迭代次数取决于系统本身的物理参数和试算向量的选取 举例 矩阵迭代法计算系统的基频及主振型 m m 2m k k 2k x1 x2 x3 系统质量矩阵和刚度矩阵 系统动力矩阵 选取初始迭代向量 r xr 1 2 8 3 7.25 4 7.189655 5 7.184652 6 7.184245 7 7.184210 系统的第一阶固有频率和主振型 r xr 1 2 7 3 7.228572 4 7.189723 5 7.184718 6 7.184253 7 7.184214 系统的第一阶固有频率和主振型 试算向量取系统静载作用时的静变形 ? 较高阶固有频率及主振型 采用动力矩阵迭代的过程，总是不断扩大第一阶主振型的比重。 能否求出第二阶以上的系统固有频率和主振型？ 对于试算初始向量 左乘 动力矩阵迭代 取 不包含有 的成分 由于计算过程中的舍入误差，x2内仍有可能存在 的残余成分 b1尽管很小，但若直接用动力矩阵A继续迭代仍然会不断扩大 的比重。 必须继续剔出它！ 设 于是有 只要在迭代计算中用矩阵A1取代A，迭代的结果会收敛到 和 且矩阵A1的特征值为λ2 ，特征向量为 而相应于 的特征值变为零 证明 当 i=1 时 即 从以上的分析，若已知系统的特征值 相应的特征 向量 。欲求出第 l+1 阶特征值