07213第四章__多自由度系統振动(4-4).docVIP

§4-4 多自由度系统的数值方法 从以上各章节的分析中可以看到，求解振动系统的固有频率和主振型，是研究各类振动问题的重要内容。当系统自由度较少时，我们可以从系统的特征方程求出其特征值（即固有频率的平方值），从而直接解出固有频率的精确值，然后再求出特征向量（即主振型）。 由此可见，寻找多自由度系统特征方程的特征值和特征向量的问题，一般都要用计算机来处理。但对特征值问题的各种近似解解法，在实际工作还常常是有用的。这里选择了几种常用的近似解法来加以介绍。 一、瑞利法 在第二章中曾介绍过，对单自由度保守系统可以用最大动能与最大势能相等的原理来求出它的固有频率。瑞利法指出，只要对位移有合理分布的假设，这种方法也可以用于求多自由度系统的最低阶固有频率（基频）。 设系统的质量矩阵、刚度矩阵和位移列阵、速度列阵分别为。则多自由度系统的动能T与势能U的表达式为： （4.100） （4.101） 当系统以第r阶固有频率和相应的主振型作第r阶主振动时，则： （4.102） 故 （4.103） 将（4.102）及（4.103）式分别代入（4.100）及（4.101）式，可得到系统在作主振动时的最大动能和势能： （4.104） （4.105） 根据机械能守恒原理，使上两式相等，即可解得系统第r阶固有频率的平方值： （4.106） 由于未求出系统的各阶固有频率之前，我们无法知道系统的各阶主振型。因此，我们不可能应用（4.106）式来计算系统的各阶固有频率。但是，我们已经知道，系统的任意一个位移振幅向量，可以用系统的各阶正则振型的线性组合来表示，即； （4.107） 所以我们可以先假设一个系统的位移振幅向量，并对它进行以下运算： （4.108） （4.109） 在进行以上运算时，我们应用了正则振型的性质，即（4.86）和（4.95）式，因此，式中的就是系统的第r阶固有频率。 将（4.108）、（4.109）式，按（4.106）式形式进行计算，可得以下关系式： （4.110） 显然，按上式所计算出的值，并不是系统的任何一阶固有频率，而是系统所有固有频率的加权平均值。我们把这值称之为瑞利商。 在瑞利商中，权因子为，而则决定了各阶主振型的影响大小。若第s阶主振型的权因子比其它各阶权因子都大得多，则显然，瑞利商值将接近于第s阶固有频率的平方值。而且，当很接近于时，则，且 瑞利法一般均用于计算系统的最低阶固有频率的近似值。因为系统的第一阶主振型 比较容易近似地确定，因而我们所选择的比较容易接近于第一阶主振型。此时， 若对（4.110）式的分子和分母项均除以，则瑞利商的计算式可写成以下形式： （4.111） 由此可见，当选择的接近于时，则按（4.110）式计算出的瑞利商可作为系统最低阶固有频率的平方的近似值。而且，这一近似值总比的精确值大。当我们所假论的越接近时，则值越接近。只有当选得恰好是时，算出的瑞利商才正好是的精确值。 二、邓柯莱法 邓柯莱法与瑞利法一样，均只用于求系统的最低阶固有频率。但瑞利法给出了基本频率的上限，而邓柯莱法给出了基本频率的下限。 在本章第一节，曾给出以柔度矩阵表示的多自由度无阻尼系统的运动位移方程为： 若系统作自由振动，则上式变为： （4.112） 设上式有简谐振动解，即 将上式代入系统的运动位移方程，并将全式除以，则得： 假设系统的质量矩阵[m]为对角阵，则根据上式可得出系统的的特征方程为： （4.113） 将上式写成展开式为： （4.114） 再将（4.114）式展开，得： （4.115） 根据代数方程的理论，上列多项式的n个根之和应等于该多项式中第二项的系数的负值，即： 因为，所以上式等式左边各项中，第一项最大。故将这一项保留，而略去其余项，则有： （4.116） 由于柔度影响系数和刚度影响系数对单自由度系统来说，就是系统的柔度和刚度，它们互为倒数。因而， （4.117） 式中，可以理解为当多自由度系统只有质量存在而变成单自由度系统时的固有频率。 将（4.117）式代入（4.116）式，得： （4.118） 上式表明，系统最低阶固有频率平方值的倒数，近似等于各质量（i=1，2，…，n）单独存在时所得固有频率平方值的倒数之和。由于（4.116）和（4.118）式的左边略去了，所以求出的就大于它的真实值，从而计算出的系统最低固有频率将小于它的精确值。 （4.118）式称为邓柯莱（Dunkerley）公式，这是邓柯莱用实验确定多圆盘的轴横向振动的固有频率时建立的。 若以表示单独安装各圆盘时在该点轴的静挠度，则： 将以上各式代入（4.118）式，可得：