-多自由度系统的数值计算方法.pptVIP

基地建设目标和总体思路 * * * * 返回首页 Theory of Vibration with Applications 瑞利(Rayleigh)能量法 李兹(Ritz)法 子空间迭代法 多自由度系统 多自由度系统的数值计算方法 返回首页 Theory of Vibration with Applications 在求解多自由度系统的固有频率和主振型的问题时，随着系统自由度数目的增加，这种求解计算工作量也随之加大。因此，通常要借助计算机进行数值计算。 常用的数值计算方法有： 瑞利法 李兹法 子空间迭代法 下面介绍这几种常用的数值计算方法及计算机的应用。 多自由度系统 多自由度系统的数值计算方法 返回首页 Theory of Vibration with Applications 多自由度系统 多自由度系统的数值计算方法—瑞利法 单自由度系统 运动微分方程 位移函数 速度函数 瑞利商： 返回首页 Theory of Vibration with Applications 设A为振型矢量，对于简谐振动，其最大动能和最大势能为 对于保守系统，由能量守恒，则有 若A是系统的第i阶主振型A(i)，则得相应的主频率的平方 若A是任意的n维矢量，则可得 称为瑞利商为了区别用位移方程求得的值，又称之为瑞利第一商。 多自由度系统 多自由度系统的数值计算方法—瑞利法 返回首页 Theory of Vibration with Applications 瑞利第一商值是否为系统某一主频率的平方，则决定于所取矢量A。如果A与某一主振型矢量接近，则所得瑞利商是相应的固有频率的近似值。实际上，对高阶振型很难做出合理的假设，而对于第一阶主振型则比较容易估计，所以此方法常用于求基频，现推证如下。 按照振型叠加的原理，系统的任何可能位移，包括假设振型，都可以描述为各阶主振型的线性组合。现取假设振型A是正则振型矢量的线性组合，即 多自由度系统 多自由度系统的数值计算方法—瑞利法 返回首页 Theory of Vibration with Applications 现取假设振型A是正则振型矢量的线性组合，即 是组合系数的列矩阵，且为非全为零的常数 Ci可用振型的正交条件求出。即 前乘 多自由度系统 多自由度系统的数值计算方法—瑞利法 返回首页 Theory of Vibration with Applications 代入 多自由度系统 多自由度系统的数值计算方法—瑞利法 返回首页 Theory of Vibration with Applications 由此可见，瑞利商的平方根是基频ω1的近似值。假设振型越接近于真实的第一阶振型，则结果越准确。通常，以系统的静变形作为假设振型，可以得到较满意的结果。 ＜＜1 由于假设振型A接近于第一阶主振型，所以有， 多自由度系统 多自由度系统的数值计算方法—瑞利法 返回首页 Theory of Vibration with Applications 可以看出，用瑞利法求出的基频近似值大于实际的基频ω1 。这是由于假设振型偏离了第一阶振型，相当于给系统增加了约束，因而增加了刚度，使求得的结果高于真实的值。 由于 ＞1 多自由度系统 多自由度系统的数值计算方法—瑞利法 返回首页 Theory of Vibration with Applications 如果采用位移方程描述系统的运动微分方程，即 前乘以 同理，若A是任意的n矢量，则有 称为瑞利第二商 若假设振型接近于第一阶主振型时，则 是基频 的近似值 给出同样假设振型的同一振动系统，用瑞利第二商计算的结果，要比用瑞利第一商计算的结果更精确一些。 多自由度系统 多自由度系统的数值计算方法—瑞利法 返回首页 Theory of Vibration with Applications 例1 用瑞利法求图示三自由度扭转系统的第一阶固有频率的估值。已知k1=k2=k3=k；I1=I2=I3=I。 解：系统的质量矩阵和刚度矩阵为 逆矩阵 计算得 求第一阶固有频率的估值，取假设振型 多自由度系统 多自由度系统的数值计算方法—瑞利法 返回首页 Theory of Vibration with Applications 在上面的计算中，假设振型比较“粗糙”，与该系统的第一阶固有频率 ，精确到第四位值的比较误差较大。 多自由度系统 多自由度系统的数值计算方法—瑞利法 返回首页 Theory of Vibration with Ap