第五讲结构动力学中的常用数值方法.doc

第五章 结构动力学中的常用数值方法 5．1．结构动力响应的数值算法 当c为比例阻尼、线性问题模态叠加最常用。但当C无法解耦，有非线性存在，有冲击作用（激起高阶模态，此时模态叠加法中的高阶模态不可以忽略）。此时就要借助数值积分方法，在结构动力学问题中，有一类方法称为直接积分方法最为常用。所识直接是为模态叠加法相对照来说，模态叠加法在求解之前，需要对原方程进行解耦处理，而本节的方法不用作解耦的处理，直接求解。（由以力学，工程中的力学问题为主要研究对象的学者发展出来的） 中心差分法的解题步骤 初始值计算 形成刚度矩阵K，质量矩阵M和阻尼矩阵C。 定初始值，，。 选择时间步长，使它满足，并计算 ，， 计算 形成等效质量阵 对阵进行三角分解 2．对每一时间步长 （1） 计算时刻t的等效载荷 （2） 求解时刻的位移 （3） 如需要计算时刻t的速度和加速度值，则 若系统的质量矩阵和阻尼矩阵为对角阵时，则计算可进一步简化。 纽马克法的解题步骤 1.初始值计算 （1）形成系统刚度矩阵K，质量矩阵M和阻尼矩阵C （2）定初始值，，。 （3）选择时间步长，参数、。并计算积分常数 ， ，， ，， ， （4）形成等效刚度矩阵 （5）矩阵进行三角分解 对第一时间步长 （1）计算时刻的等效载荷 （2）求解时刻的位移 （3）计算时刻的加速度和速度 威尔逊-法的解题步骤 初始值计算 （1）形成系统刚度矩阵K，质量矩阵M和阻尼矩阵C （2）定初始值，，。 （3）选择时间步长，并计算积分常数 ，， ，， ，， （4）形成等效刚度 （5）将等效刚度进行三角分解 2.对每一个时间步长 （1）计算时刻的等效载荷 （2）求解时刻的位移 （3）计算在时刻的加速度、速度和位移 5．2 结构动力响应数值算法性能分析 对公式（5.1）描述的线性系统结构动力学问题，已经有证明对整个多自由度的积分，等价于将模态分解后对单自由度的积分的结果进行模态叠加，因此可以通过对单自由度问题的分析，来说明算法的特性，其中阻尼均假设为比例阻尼，这样，模态分解后的单自由度结构动力学方程为： （5-29） 以下算法的性能分析，均将算法用于这个方程。 5．2．1 算法用于结构动力学方程的有限差分表示 将数值计算方法应用于（5-29）, 即分别在相邻的不同时刻应用算法可得如下一般形式 （5-30） A为放大矩阵或称逼近算子，为载荷逼近算子。 例如将Newmak方法应用于方程（5-29）有： （5-33） ， 矩阵A的特征多项式为 （5-34） 其中A1，A为该矩阵的两个特征向量，分别为矩阵的迹的一半和矩阵的行列式 （5-35） （5-36） （5-37） 对Newmak方法有： , （5-38） 其中 h为时间步长，。Newmak方法放大矩阵的规模是二维的，因此特征值也只有两个，可以根据它们进行分析。 有的算法放大矩阵是三维的，例如Wilson-方法,在无阻尼情况下放大矩阵为： （5-39） 放大矩阵A的特征多项式为： （5－40） 其中A1，A2，A3为该矩阵的三个特征向量，分别为矩阵的迹的一半、各阶主子式的和以及矩阵的行列式，对Wilson-方法有 （5-41） 此外，在几个不同时刻应用数值算法，然后将方程中的速度和加速度项消去，可得数值算法关于位移的差分方程，例如Newmak方法，有 （5-42） 很显然，其特征方程与其放大矩阵A的特征方程是相同的，使用关于位移的线性多步方式和放大矩阵来说明算法性能是一样的，只不过各有方便之处。 5．2．2 算法的稳定性分析 设为放大矩阵A的特征值，则定义为A的谱半径，若特征值互异，则的算法是稳定的，但若有重特征根，则要求。如果算法的稳定性要求对步长的选取有限制，称算法是有条件稳定的，反之为无条件稳定的。放大矩阵的谱半径小于等于1成立的充分