振动力学3new.pptVIP

第三章 线性振动近似计算方法 第一节  邓克利法 第二节  瑞利法和里茨法 作业： P123-5.3 第三节  矩阵迭代法 作业： P123-5.6 第四节 子空间迭代法 选作： P123-5.8 第五节  传递矩阵法 选作： P 123-5.9 , 5.10 第六节广义本征值问题的其它解法 小结： 邓克利法 瑞利法与里茨法 矩阵迭代法 子空间迭代法 传递矩阵法 广义本征值问题的解法 解：两端无约束的边界条件为 令 代入各单元状态变量的第一个元素，得到模态 2.梁的横向振动 用传递矩阵法可以分析梁的横向振动。 梁的支座、梁上的集中质量以及联结支座或集中质量的梁段，将第 和第 集中质量以及联结二者的第 梁段作为第 单元，梁段质量不计。 设支座或集中质量处梁的位移 和截面转角 ,以及弯矩 和剪力 组成状态变量X 第 单元的集中质量为 、梁的长度和抗弯刚度分别为 和 ，下角标表示集中质量和梁段的编号，上角标L和R表示左侧和右侧截面. 根据受力状况，第 质量两侧的状态变量满足 以下关系 设梁的横向振动固有频率 将 代入，导出第 质量左右两侧的传递关系 点传递矩阵为 第 梁段上力的平衡得 梁段两端的位移和转角的关系 导出第 梁段左右两端状态变量之间的传递关系 场传递矩阵为 得出从第 质量右侧至第 质量右侧的状态变量传递关系 其中 将式 代入 单元传递矩阵为 对于带n个集中质量的梁，总能利用各单元传递矩阵的连乘积导出梁的最左端和最右端的传递关系 利用两端边界条件可以确定固有频率和模态。 例题2 用传递矩阵法计算带集中质量简支梁的固有频率。 解：引入量纲一的量 量纲一的状态变量为 点传递矩阵和场传递矩阵为 两端支座之间的传递矩阵 根据两端支座的边界条件 利用此齐次线性代数方程的非零解条件导出本征方程 因此实际应用传递矩阵法时并不进行符号推演，而是先假设一系列试算频率代入各传递矩阵进行数值运算作出 曲线，能满足 的试算频率即固有频率。在图中固有频率由 曲线与横坐标的交点确定。 第二次迭代后得到 归一化后为 精确解 终止迭代， 为第一阶模态 。利用 和 的最后一个元素计算基频,得 归一化后 2.高阶模态及固有频率 可利用同样的方法计算第二阶模态 和频率 。 只要在每次迭代时在试算的假设模态 中将含 的部分剔除，则迭代将收敛到第二阶模态 。利用模态正交性，将式 左乘 ,将 写作 ，推出 在第一次迭代的公式 中减去第一阶模态成分，将式 代入，得 其中矩阵 进行k次迭代后得 迭代公式 可推知，继续迭代的结果必收敛到第二阶模态和频率。 注意： 用同样方法原则上还可以计算更高阶的模态和频率，但由于计算误差的积累，只有前几阶模态和频率有足够的精度。 与公式 比较， 例题2 用矩阵迭代法计算系统的第二阶模态和固有频率。 解：从例题1已解出 第一次迭代后 选取 终止迭代， 为第一阶模态 。利用 和 的最后一个元素计算第二阶频率,得 归一化后 对于高阶频率的迭代与基频情形相比，收敛速度要慢得多。 对于频率方程有重根情形，设 ，则经过k次迭代后式 迭代得到的一阶模态 的线性组合，记作 将 代替式 中的 继续迭代，最终得到 与正交的模态 。 重本征值情形的模态不唯一， 与 是与同一本征值 相对应的两组不同的模态。 3.半正定系统的矩阵迭代法 半