《多自由度系统的数值计算方法》由会员分享，可在线阅读，更多相关《多自由度系统的数值计算方法（37页珍藏版）》请在金锄头文库上搜索。
1、 返回首页Theory of Vibration with Applications瑞利瑞利瑞利瑞利( (RayleighRayleigh) )能量法能量法能量法能量法李兹李兹李兹李兹(Ritz)(Ritz)法法法法子空间迭代法子空间迭代法子空间迭代法子空间迭代法1 返回首页Theory of Vibration with Applications 在求解多自由度系统的固有频率和主振型的问题时，在求解多自由度系统的固有频率和主振型的问题时，随着系统自由度数目的增加，这种求解计算工作量也随之随着系统自由度数目的增加，这种求解计算工作量也随之加大。因此，通常要借助计算机进行数值计算。加大。因此，通
2、常要借助计算机进行数值计算。 常用的数值计算方法有：常用的数值计算方法有： 瑞利法瑞利法 李兹法李兹法 子空间迭代法子空间迭代法 下面介绍这几种常用的数值计算方法及计算机的应用。下面介绍这几种常用的数值计算方法及计算机的应用。2 返回首页Theory of Vibration with Applications单自由度系统单自由度系统运动微分方程运动微分方程位移函数位移函数速度函数速度函数瑞利商：瑞利商：3 返回首页Theory of Vibration with Applications设A为振型矢量，对于简谐振动，其最大动能和最大势能为对于保守系统，由能量守恒，则有若A是系统的第i阶主振型
3、A(i)，则得相应的主频率的平方若A是任意的n维矢量，则可得称为瑞利商为了区别用位移方程求得的值，又称之为瑞利第一商。 4 返回首页Theory of Vibration with Applications瑞利第一商值是否为系统某一主频率的平方，则决定于所取矢量A。如果A与某一主振型矢量接近，则所得瑞利商是相应的固有频率的近似值。实际上，对高阶振型很难做出合理的假设，而对于第一阶主振型则比较容易估计，所以此方法常用于求基频，现推证如下。按照振型叠加的原理，系统的任何可能位移，包括假设振型，都可以描述为各阶主振型的线性组合。现取假设振型A是正则振型矢量的线性组合，即5 返回首页Theory of
4、 Vibration with Applications现取假设振型A是正则振型矢量的线性组合，即是组合系数的列矩阵，且为非全为零的常数Ci可用振型的正交条件求出。即前乘6 返回首页Theory of Vibration with Applications代入7 返回首页Theory of Vibration with Applications由此可见，瑞利商的平方根是基频1的近似值。假设振型越接近于真实的第一阶振型，则结果越准确。通常，以系统的静变形作为假设振型，可以得到较满意的结果。 1由于假设振型A接近于第一阶主振型，所以有，8 返回首页Theory of Vibration with
5、Applications可以看出，用瑞利法求出的基频近似值大于实际的基频1 。这是由于假设振型偏离了第一阶振型，相当于给系统增加了约束，因而增加了刚度，使求得的结果高于真实的值。 由于 19 返回首页Theory of Vibration with Applications如果采用位移方程描述系统的运动微分方程，即前乘以同理，若A是任意的n矢量，则有称为瑞利第二商若假设振型接近于第一阶主振型时，则 是基频 的近似值 给出同样假设振型的同一振动系统，用瑞利第二商计算的结果，要比用瑞利第一商计算的结果更精确一些。 10 返回首页Theory of Vibration with Applicatio
6、ns例1 用瑞利法求图示三自由度扭转系统的第一阶固有频率的估值。已知k1=k2=k3=k；I1=I2=I3=I。 解：系统的质量矩阵和刚度矩阵为逆矩阵计算得求第一阶固有频率的估值，取假设振型11 返回首页Theory of Vibration with Applications在上面的计算中，假设振型比较“粗糙”，与该系统的第一阶固有频率 ，精确到第四位值的比较误差较大。12 返回首页Theory of Vibration with Applications如果进一步改进假设振型，即以静变形曲线为假设振型，如设显然，在工程上，若以静变形曲线作为假设振型，可以得到很好的第一阶固有频率的近似值。1
7、3 返回首页Theory of Vibration with Applications用瑞利法估算的基频的精度取决于假设的振型对第一阶主振型的近似程度，而且其值总是精确值的上限。李兹法对近似振型给出更合理的假设，从而使算出的基频值进一步下降，并且还可得系统较低的前几阶固有频率及相应的主振型在李兹法中，系统的近似主振型假设为是选取的s个线性独立的假设振型矩阵s维待定系数14 返回首页Theory of Vibration with Applications由于 在系统的真实主振型处取驻值，这些驻值即相应的各阶固有频率的平方 ，所以a的各元素由下式确定15 返回首页Theory of Vibrat
8、ion with Applicationsn个自由度缩减至s 自由度。刚度矩阵质量矩阵16 返回首页Theory of Vibration with Applications李兹法是一种缩减系统自由度数的近似方法。频率方程求出s个固有频率，即n自由度系统的前s阶固有频率。解出其相应的特征矢量求出n自由度系统的前s阶主振型正交性17 返回首页Theory of Vibration with Applications对于瑞利第二商利用驻值条件可得s个方程，将其写成矩阵形式特征方程求出s个固有频率，即n自由度系统的前s阶固有频率。解出其相应的特征矢量求出n自由度系统的前s阶主振型18 返回首页The
9、ory of Vibration with Applications例2 用李兹法求图示四自由度振动系统的前二阶固有频率及主振型。 解：由条件可求出系统的质量矩阵、刚度矩阵和柔度矩阵设振型19 返回首页Theory of Vibration with Applications求出求出2个固有频率，即4自由度系统的前2阶固有频率。20 返回首页Theory of Vibration with Applications求出系统的前二阶主振型21 返回首页Theory of Vibration with Applications将将矩矩阵阵迭迭代代法法与与李李兹兹法法结结合合起起来来，可可以以得得到
10、到一一种种新新的的计计算算方法，即子空间迭代法。方法，即子空间迭代法。它它对对求求解解自自由由度度数数较较大大系系统统的的较较低低的的前前若若干干阶阶固固有有频频率率及及主振型非常有效。主振型非常有效。22 返回首页Theory of Vibration with Applications计算系统的前P阶固有频率和主振型，按照李兹法，可假设s个振型且sP。将这些假设振型排列成ns阶矩阵，即其中每个 都包含有前P阶振型的成分，也含有高阶振型的成分。23 返回首页Theory of Vibration with Applications为了提高李兹法求得的振型和频率的精确度，将A0代入动力矩阵中进
11、行迭代，并对各列阵分别归一化后得这样做的目的是使 比A0含有较强的低阶振型成分，缩小高阶成分。但如果继续用 进行迭代，所有各阶振型即 的各列都将趋于A(1)。 为了避免这一点，可以在迭代过程中进行振型的正交化。用李兹法进行振型正交化具有收敛快的特点。因为它是利用瑞利取驻值的条件，寻求s2个aij的系数，使得 的每一列都成为相对应振型A(i)的最佳近似。24 返回首页Theory of Vibration with Applications所以用 作为假设振型，再按李兹法求解，即设可求得广义质量矩阵和广义刚度矩阵ss阶待定系数方阵得到s个值 及对应的特征矢量 再由李兹法特征值问题，即求解方程从而
12、求出25 返回首页Theory of Vibration with Applications然后，以求出的 作为假设振型进行迭代，可求得与李兹法特征值问题，解出 。由李兹法，即不断地重复矩阵迭代和李兹法的过程，就可以得到所需精度的振型和固有频率。26迭代的功能是使这s个矢量的低阶成分不断地相对放大，即向 张成的子空间靠拢。 返回首页Theory of Vibration with Applications 子空间迭代法是对一组假设振型反复地使用迭代法和李兹法的运算。 从几何观点上看，原n阶特征值系统有n个线性无关的特征矢量，它们之间是正交的，张成一个n维空间。而假设的s个线性无关的n维矢量 张
13、成一个s 维子空间，27 如果只迭代不进行正交化，最后这s个矢量将指向同一方向，即A(1)的方向。 由于用李兹法作了正交处理，则这些矢量不断旋转，最后分别指向前s个特征值的方向。 返回首页Theory of Vibration with Applications即由张成的一个s 维子空间，经反复地迭代正交化的旋转而逼近于由所张成的子空间。28 返回首页Theory of Vibration with Applications 在实践中发现，最低的几阶振型一般收敛很快，经过二至三次迭代便已稳定在某一数值。 在以后的迭代中不能使这几个低阶振型值的精度进一步提高，只是随着迭代次数的增加，将有越来越多
14、的低阶振型值稳定下来。 所以，在计算时要多取几个假设振型，如果所需求的是P个振型，则假设振型个数s一般应在2P与2P+8之间取值。29 返回首页Theory of Vibration with Applications 子空间迭代法有很大的优点，它可以有效地克服由于等固有频率或几个频率非常接近时收敛速度慢的困难。 同时，在大型复杂结构的振动分析中，系统的自由度数目可达几百甚至上千，但是，实际需用的固有频率与主振型只是最低的三、四十个，通常对此系统要进行坐标缩聚。 与其它方法相比，子空间迭代法具有精度高和可靠的优点。因此，它已成为大型复杂结构振动分析的最有效的方法之一。 30 返回首页Theor
15、y of Vibration with Applications例3 用子空间迭代法求例2中所示系统的前二阶固有频率及振型。解：系统的质量矩阵、刚度矩阵和柔度矩阵已由例2求出。现取假设振型由动力矩阵迭代得到31 返回首页Theory of Vibration with Applications将各列分别归一化得求得再由李兹法特征值问题为其中 32 返回首页Theory of Vibration with Applications由上述方程有非零解的条件，得频率方程为各列分别归一化后，得33 返回首页Theory of Vibration with Applications重复上述过程进行第二次迭代，由归一化后得则有由34 返回首页Theory of Vibration with Applications解得得频率方程为由于 近似于单位矩阵，所以有35 返回首页Theory of Vibration with Applications由于 近似于单位矩阵，所以有结束迭代，求得系统的前二阶固有频率及相应的主振型为36Theory of Vibration with Applications习题习题4.1137
