多自由系统数值计算方法.pptxVIP

第5章多自由度系统的数值计算方法 天津大学 5.1 瑞利(Rayleigh)能量法 5.2 李兹(Ritz)法 5.3 邓克莱(Dunkerley)法 5.4 矩阵迭代法 5.5 子空间迭代法 5.6 传递矩阵法 Mechanical and Structural Vibration 第5章多自由度系统的数值计算方法 5.1 瑞利(Rayleigh)能量法Mechanical and Structural Vibration 5.1 瑞利(Rayleigh)能量法5.1.1瑞利第一商5.1.2瑞利第二商 Mechanical and Structural Vibration 若A是系统的第i阶主振型A(i)，则得相应的主频率的平方 5.1 瑞利(Rayleigh)能量法5.1.1瑞利第一商设A为振型矢量，对于简谐振动，其最大动能和最大势能为对于保守系统，由能量守恒，则有若A是任意的n维矢量，则可得称为瑞利商为了区别用位移方程求得的值，又称之为瑞利第一商。 Mechanical and Structural Vibration 由于 ＞1 5.1 瑞利(Rayleigh)能量法5.1.2瑞利第二商 瑞利商的平方根是基频p1的近似值。假设振型越接近于真实的第一阶振型，则结果越准确。通常，以系统的静变形作为假设振型，可以得到较满意的结果。 用瑞利法求出的基频近似值大于实际的基频p1 。这是由于假设振型偏离了第一阶振型，相当于给系统增加了约束，因而增加了刚度，使求得的结果高于真实的值。 Mechanical and Structural Vibration 前乘以若假设振型接近于第一阶主振型时，则 是基频 的近似值 5.1 瑞利(Rayleigh)能量法5.1.2瑞利第二商 如果采用位移方程描述系统的运动微分方程，即同理，若A是任意的n矢量，则有称为瑞利第二商给出同样假设振型的同一振动系统，用瑞利第二商计算的结果，要比用瑞利第一商计算的结果更精确一些。 Mechanical and Structural Vibration 逆矩阵 5.1 瑞利(Rayleigh)能量法5.1.2瑞利第二商 例5-1 用瑞利法求图示三自由度扭转系统的第一阶固有频率的估值。已知k1=k2=k3=k；I1=I2=I3=I。 解：系统的质量矩阵和刚度矩阵为求第一阶固有频率的估值，取假设振型计算得Mechanical and Structural Vibration 在上面的计算中，假设振型比较“粗糙”，与该系统的第一阶固有频率 ，精确到第四位值的比较误差较大。 5.1 瑞利(Rayleigh)能量法5.1.2瑞利第二商 Mechanical and Structural Vibration 5.1 瑞利(Rayleigh)能量法5.1.2瑞利第二商 如果进一步改进假设振型，即以静变形曲线为假设振型， 设显然，在工程上，若以静变形曲线作为假设振型，可以得到很好的第一阶固有频率的近似值。Mechanical and Structural Vibration 第5章多自由度系统的数值计算方法 5.2 李兹(Ritz)法Mechanical and Structural Vibration s维待定系数矩阵 5.2 李兹(Ritz)法用瑞利法估算的基频的精度取决于假设的振型对第一阶主振型的近似程度，而且其值总是精确值的上限。李兹法对近似振型给出更合理的假设，从而使算出的基频值进一步下降，并且还可得系统较低的前几阶固有频率及相应的主振型在李兹法中，系统的近似主振型假设为是选取的s个线性独立的假设振型Mechanical and Structural Vibration 由于 在系统的真实主振型处取驻值，这些驻值即相应的各阶固有频率的平方 ，所以a的各元素由下式确定 5.2 李兹(Ritz)法Mechanical and Structural Vibration 5.2 李兹(Ritz)法刚度矩阵质量矩阵n个自由度缩减至s 自由度。Mechanical and Structural Vibration 5.2 李兹(Ritz)法李兹法是一种缩减系统自由度数的近似方法。频率方程求出s个固有频率，即n自由度系统的前s阶固有频率。解出其相应的特征矢量求出n自由度系统的前s阶主振型正交性Mechanical and Structural Vibration 5.2 李兹(Ritz)法对于瑞利第二商利用驻值条件可得s个方程，将其写成矩阵形式特征方程求出s个固有频率，即n自由度系统的前s阶固有频率。解出其相应的特征矢量求出n自由度系统