結构动力学理论及其在地震工程中的应用.docxVIP

5章 动力反应的数值计算 如果激励[作用力或地面加速度]是随时间任意变化的，或者体系是非线性的，那么对单自由度体系的运动方程进行解析求解通常是不可能的。这类问题可以通过数值时间步进法对微分方程进行积分来处理。在应用力学广阔的学科领域中，有关各种类型微分方程数值求解方法的文献（包括几部著作中的主要章节）浩如烟海，这些文献包括这些方法的数学进展以及它们的精度、收敛性、稳定性和计算机实现等问题。 然而，本章仅对在单自由度体系动力反应分析中特别有用的很少几种方法进行简要介绍，这些介绍仅提供这些方法的基本概念和计算算法。尽管这些对许多实际问题和应用研究已经足够了，但是读者应该明白，有关这个主题存在大量的知识。 5.1 时间步进法 对于一个非弹性体系，欲采用数值求解的运动方程为 或者 (5.1.1) 初始条件 假定体系具有线性粘滞阻尼，不过，也可以考虑其他形式的阻尼（包括非线性阻尼），后面会明显看到这一点。然而由于缺乏阻尼信息．因此很少这样做，特别是在大振幅运动时。作用力由一系列离散值给出: , 到N。时间间隔 （5.1.2） 图5.1.1 时间步进法的记号 通常取为常数，尽管这不是必需的。在离散时刻（表示为时刻）确定反应，单自由度体系的位移、速度和加速度分别为、和。假定这些值是已知的，它们在时刻满足方程 （5.1.3） 式中，是时刻的抗力，对于线弹性体系，，但是如果体系是非弹性的，那么它会依赖于时刻以前的位移时程和速度。将要介绍的数值方法将使我们能够确定+1时刻满足方程(5.1.1)的反应、和，即在+1时刻 （5.1.4） 对于=0，1，2，3，…，连续使用时间步进法，即可给出=0，l，2，3，…所有瞬时所需的反应。已知的初始条件和提供了起动该方法的必要信息。 从时刻到+1时刻的步进一般不是精确的方法，许多在数值上可以实现的近似方法是可能的。对于数值方法，有三个重要的要求：(1)收敛性一随着时间步长的减少，数值解应逼近精确解；(2)稳定性一在存在数值舍入误差的情况下，数值解应是稳定的；(3)精度一数值方法应提供与精确解足够接近的结果。这些重要的问题在本书中均作简要的讨论，全面的论述可在着重微分方程数值解法的书中找到。 本章介绍三种类型的时间步进法：(1)基于激励函数插值的方法；(2)基于速度和加速度有限差分表达的方法；(3)基于假设加速度变化的方法。前两类中各只介绍一种方法，第三类中介绍两种方法。 5.2 基于激励插值的方法 对于线性体系，通过在每个时间间隔里对激励进行插值，并利用第4章的方法进行精确求解，能推导出一种非常有效的数值方法。如果时间间隔较短，则线性插值是令人满意的。图5. 2.1所示的时间间隔，激励函数为 （5.2.1a） 其中 （5.2.1b） 时间变量从0到变化。为数学上简单起见，我们首先考虑无阻尼体系，后面再将该方法扩展到有阻尼体系。待求解的方程为 （5.2.2） 在时间间隔内，反应为三部分之和：(1) =0时刻的初位移和初速度引起的自由振动；(2)零初始条件下对阶跃力的反应；(3)零初始条件下对斜坡力的反应。对这三种情况分别采用来自§2.1、§4.3和§4.4中已有的解答，得 （5.2.3a） （5.2.3b） 计算时的这些等式，得+1时刻的位移和速度： (5.2.4a) (5.2.4b) 将式（5.2.1b）代入后，可将这些等式重写为如下的递推公式： (5.2.5a) (5.2.5b) 对于欠临界阻尼体系（即），重复上面的推导，表明式(5.2.5)也适用于有阻尼体系，系数A，B，…,的表达式由表5.2.1给出，这些系数取决于体系的参数、k和ζ以及时间间隔。 表5.2.1 递推公式中的系数（）