[工学]第五章多自由自由度系统的数值解法.pptVIP

第一节： Rayleigh 法 第二节：Dunkerley 法 第三节：矩阵迭代法 第四节 传递矩阵法 第五节 梁 将（1）（2）（3）（4）写成矩阵形式： （5） 弹性元件的场 矩阵＝〔F〕 3.质量元件的点矩阵 4.分段传递矩阵 （5）代入（6）得： （7） 分段传递矩阵 （6） 点传递矩阵＝〔P〕 5.系统传递矩阵 利用方程（7）可以从系统的左边边界－－位置0算起，计算到右边边界－－位置n，可得： （8） 系统传递矩阵 关于梁，通常的边界条件有： 简支 自由 固定 利用边界条件可 解方程式（8） 多自由度系统的数值方法 第5章 一、刚度法 个自由度无阻尼系统特征值问题方程： (1) 即有： (2) 方程(2)两边前乘 得： (3) 上式：分子与第r阶势能有关，分母与第r阶模态动能有关如果是 任一向量 ，则 (4) 称为Rayleigh商。若 是 ，则 即为 系统特征向量 ，形成n维空间中一组线性 独立的完备系。 同一空间中的任一向量{u}，可用特征向量的线性组合来表示 ， 即： (5) 式中： 是常数 将(5)代入(4)，同时考虑到： 得： (6) 由(6)知： 是特征值 加权平均。 如果 与第r阶特征向量 很接近，则意味着 与 相比较是很小的， 有： (8) (7) 式中： (9) 的误差为二阶。 即表明：Rayleigh商在特征值向量的领域中有稳定的值。 二、用Rayleigh法计算系统第一阶固频 当r=1时，由方程（8）可得： （10） （11） Rayleigh商大于真值 注意：计算 模态向量真值 就可以求得比较精确的勤 值。 n个自由度特征值问题方程为： (1) (2) 以两个自由度系统为例来说明二个自由度系统，其特征值问题方程为： (3) (4) (5) 对比(4)(5)两式可知： (6) (7) (8) (9) (10) 例：质量为 ，长为 的均质悬臂梁，其端部有集中质量 ，试确定 。假设梁的弯曲刚度为EI。 EI 解：对于端部带有质量m，而略去梁本身质量的悬臂梁，其固频为： 均质悬臂梁的第一阶固频为： 一、特征值问题 n个自由度无阻尼系统特征问题方程为： 可表示为： (1) (2) 叫动力矩阵。 为了求得系统特征值和特征向量，可用(2)式进行迭代。 二、迭代原理： (3) (4) 式中： 且随着r的增加而减小。 试探向量， 相比，高阶模态在 中的作用在减小。 是我们为了得到系统的特 征向量 而构造的 可以认为是 的一个改进。 为了改进 ，令： (5) 显然 比 和 更好地试探向向量，重复上述过程，则有： (6) 如果迭代数S有足够的大： (7) 由方程(7)知：当S足够大时，级数(6)的第一项是决定性的，级数将收敛于 ，即 和 都可认为是 满足方程(2)。 和 相互成比例，比例常数为 ，由此得： 和 三、迭代法求高阶固有模态的特征值和特征向量 由变换 由上式得主坐标{p}与广义坐标的关系： (8) 如果能从系统特征值问题方程(2)中消去第一阶固有模态的影响，则通过迭代将得到第二阶固有模态的特征值和特征向量。消去第一阶固有模态的影响，就是使第一阶主坐标 由方程(8)得： (9) 以三个自由度系统为例来进行讨论：(9)式为： 将上式写成矩阵形式： 上式可变换为： (10) (11) 清型矩阵 令： (12) 由于矩阵[G]已被矩阵 修改，第一阶固有模态已被消除。 代入(2)进行迭代，可以得到第二阶固有模态的特征值和特征向量。 四、 的求法 由(11)知，要求得 ，关键是求得 而 的解不可能由 求得 。 外振型矩阵[u]中的其它列向量仍是未知 的。 利用特征向量的正交关系 ： 上式右乘以 (13) (14) ∵ 为对角阵， ∴方程(13)的第一行为： 为得第三阶和更高阶固有模态的特征值和特征向量，可用类似方法，依次建立起相应的矩阵 ，得到修改后的 ，依次代入方程(2)进行迭代。 例：有一个三自由度系统，运动方程为： 用迭代法求其固频和特征