第四章实用计算方法要点分析.ppt

第四章实用计算方法1第四章实用计算方法§4.2矩阵特征值问题及解法§4.3结构动力响应的数值解法2§4.1能量法求自振频率弹簧变形能：质量块动能：§4.1能量法求自振频率一、瑞利能量法Rayleigh法的理论基础为能量守恒定律。即认为如果没有阻尼力消耗能量的话，在自由振动体系中，能量应该保持常量。最大动能等于最大位能：这个表达式和以前所述的一样，但现在它是从最大变形能应等于最大动能的Rayleigh法概念而得。例子：简支梁，认为是无限自由度体系变形能：最大值：体系动能：由Rayleigh法：最大值：→k*→m*此即为瑞利商例子：简支梁，认为是无限自由度振动形状的选取假定振型为抛物线：能量守恒：假定振型为正弦曲线：能量守恒：原则上，只要满足梁的几何边界条件，形状函数可任意选取，亦即形状函数仅需和具体的支承条件一致。但是，对不是真实振型的任意形状函数，为了保持平衡就必须有附加的外部约束作用，这些附加约束将会使体系变得刚硬，从而使计算频率增大。Rayleigh法计算的频率中，最低的一个，总是最好的近似值！Question:如何确定合理的挠曲形状？Solution:自由振动的位移是由惯性力作用引起的；惯性力正比于质量×加速度（质量分布及位移幅值）因此：正确的振动形式为正比于m(x)的荷载所引起的挠曲线。能量守恒：频率计算将根据静止重量荷载所引起的挠曲线vd(t)进行。此时，体系的变形能必然等于重量荷载所做的功。例假定变形曲线最大位能最大动能Finish?计算频率:例：试用瑞利法求图示楔形悬臂梁的基本频率。宽度b=1。解：设形状函数为满足位移边界条件。精确解为例：试求图示对称刚架的基本频率。解：柱的最大动能梁的最大动能刚架的最大动能例：试求图示等截面悬臂梁的基本频率。解：1.设形状函数为余弦曲线精确解为2.设形状函数为抛物线3.设形状函数为重力引起的位移曲线形状函数为重力引起的位移曲线时20二、李兹能量法用瑞利法求解结构的自振频率的精度取决于假设振型的精度，由于难以估计高阶振型的形状，所以一般情况下，瑞利法只能求得振动基频的上限。李兹发展了瑞利能量法，使求得的最低频率更接近于精确解，且可以求较高阶频率。李兹能量法给出振型的级数形式：（1）将（1）代人瑞利商：得：21引入符号：（2）（3）则（2）式为：（4）2223即：（5）（6）（7）例：试用李兹法求图示楔形悬臂梁的基本频率。宽度b=1。解：设位移函数为例：试求图示楔形悬臂梁的基本频率。宽度b=1。解：精确解为瑞利法解误差：0.07%误差：3%§4.2矩阵特征值问题及解法对多自由度体系的自由振动来说，采用数值分析的方法，求解其固有频率和振型最终可以变成数学上的固有特征值问题（eigenvalueproblem），即26求固有频率和主振型固有特征值方程的解:1)给出了结构的自振频率和主振型，2)使结构在动力荷载作用下的运动方程解耦，即振型叠加法。振型叠加法:只要频率和振型确定以后就可以得到线性结构的动力响应，对于一般的工程问题，只需取前几个振型就能获得具有相当精度的响应解。§4.2矩阵特征值问题及解法1、化广义特征值问题为标准特征值问题27广义特征值问题转化为标准特征值问题进行求解。§4.2矩阵特征值问题及解法将质量矩阵[M]进行乔累斯基（Cholesky）分解，即将其分解为[L]与[L]T的乘积:28广义特征值问题[L]为对角元素不为零的下三角矩阵。标准特征值问题§4.2矩阵特征值问题及解法[K]是对称的，矩阵[A]也具有对称性。所有对称矩阵特征值问题的算法均可以得到利用。如果矩阵[K]是正定的，也可将其进行乔累斯基分解，得到类似于方程的标准特征值问题。29§4.2矩阵特征值问题及解法30Cholesky分解:§4.2矩阵特征值问题及解法则其下三角i=1,2,…,n行的元素为:31§4.2矩阵特征值问题及解法2、标准特征值问题解法特征值问题:一是求解它的全部特征值问题，即所有的特征值和对应的特征向量；另一是求解它的部分特征值问题，即部分（通常是最小或最大的一部分）特征值和对应的特征向量。32在结构动力学中，往往矩阵的阶数都很高，有时不可能，而且也没必要求解全部特征值和特征向量。求解方法分为两大类:一类是变换法,另一类是向量迭代法。雅可比法幂法子空间迭代法§4.2矩阵特征值问题及解法1）雅可比法雅可比法是求解实对称矩阵全部特征值和特征向量的简单有效方法。33[A]为n×n阶的实对称矩阵。根据线性代数理论，任何一个n×n实对称矩阵[A]，可通过一个的正交矩阵[S]，经相似变换化为一对角矩阵，即矩阵[D]的n个对角元素就是[A]的n个特征值，而[S]的第i列，就是[D]中第i个元素所对应的特征向量。§4.2矩阵特征值问题及解法正交矩阵[S]实际上是一个坐标变换矩阵34功能:实现坐标系的旋