结构动力学理论的和其在地震工程的中地应用的

《结构动力学理论的和其在地震工程的中地应用的》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、5章动力反应的数值计算如果激励[作用力或地面加速度]是随时间任意变化的，或者体系是非线性的，那么对单自由度体系的运动方程进行解析求解通常是不可能的。这类问题可以通过数值时间步进法对微分方程进行积分来处理。在应用力学广阔的学科领域中，有关各种类型微分方程数值求解方法的文献（包括几部著作中的主要章节）浩如烟海，这些文献包括这些方法的数学进展以及它们的精度、收敛性、稳定性和计算机实现等问题。然而，本章仅对在单自由度体系动力反应分析中特别有用的很少几种方法进行简要介绍，这些介绍仅提供这些方法的基本概念和计算算法。尽管这些对许多实际问题和应用研究已经足够了，但是读者应该明白，有关这个主题存在大量

2、的知识。5.1时间步进法对于一个非弹性体系，欲采用数值求解的运动方程为或者(5.1.1)初始条件假定体系具有线性粘滞阻尼，不过，也可以考虑其他形式的阻尼（包括非线性阻尼），后面会明显看到这一点。然而由于缺乏阻尼信息．因此很少这样做，特别是在大振幅运动时。作用力由一系列离散值给出:,到N。时间间隔（5.1.2）图5.1.1时间步进法的记号通常取为常数，尽管这不是必需的。在离散时刻（表示为时刻）确定反应，单自由度体系的位移、速度和加速度分别为、和。假定这些值是已知的，它们在时刻满足方程（5.1.3）式中，是时刻的抗力，对于线弹性体系，，但是如果体系是非弹性的，那么它会依赖于时刻以前的位移时

3、程和速度。将要介绍的数值方法将使我们能够确定+1时刻满足方程(5.1.1)的反应、和，即在+1时刻（5.1.4）对于=0，1，2，3，…，连续使用时间步进法，即可给出=0，l，2，3，…所有瞬时所需的反应。已知的初始条件和提供了起动该方法的必要信息。从时刻到+1时刻的步进一般不是精确的方法，许多在数值上可以实现的近似方法是可能的。对于数值方法，有三个重要的要求：(1)收敛性一随着时间步长的减少，数值解应逼近精确解；(2)稳定性一在存在数值舍入误差的情况下，数值解应是稳定的；(3)精度一数值方法应提供与精确解足够接近的结果。这些重要的问题在本书中均作简要的讨论，全面的论述可在着重微分方程

4、数值解法的书中找到。本章介绍三种类型的时间步进法：(1)基于激励函数插值的方法；(2)基于速度和加速度有限差分表达的方法；(3)基于假设加速度变化的方法。前两类中各只介绍一种方法，第三类中介绍两种方法。5.2基于激励插值的方法对于线性体系，通过在每个时间间隔里对激励进行插值，并利用第4章的方法进行精确求解，能推导出一种非常有效的数值方法。如果时间间隔较短，则线性插值是令人满意的。图5.2.1所示的时间间隔，激励函数为（5.2.1a）其中（5.2.1b）时间变量从0到变化。为数学上简单起见，我们首先考虑无阻尼体系，后面再将该方法扩展到有阻尼体系。待求解的方程为（5.2.2）在时间间隔内，

5、反应为三部分之和：(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递推公式

6、中的系数（）因为递推公式是从运动方程的精确解推导出的，因此对时问步长大小的唯一限制条件是，允许它对于激励函数有一个接近的逼近，并以较密的时间间隔提供反应结果，以使反应峰值不会被漏掉。这种数值方法对于激励由紧密的时间间隔定义的情况（例如对于地震地面加速度的情况）特别有用，从而使得线性捕值即可得到较完美的结果。如果时间步长是常数，则系数A，B，…，仅需计算一次。这种数值方法所要求的运动方程精确解答仪对线性体系是可行的。如上所述，这种方法用于单白由度体系较便利，但是对于多自由度体系则是不切实际的，除非它们的反应由振型反应的叠加（第12章和第13章）来获得。5.3中心差分法这种方法是基于对位移

7、时间导数（即速度和加速度）的有限差分近似进行的。步长，则时刻的速度和加速度的中心差分表达式为(5.3.1)将速度和加速度的这些近似表达式代人方程(5.1.3)中，对线弹性体系，得(5.3.2)在这个方程中，和假定是已知的（来自于前面时间步内方法的执行）。将这些已知量移到右侧，导得(5.3.3)或写成(5.3.4)其中（5.3.5）（5.3.6）则未知的由下式给出（5.3.7）+l时刻的解答是根据时刻的平衡条件即方程(5.1.3)确定的，而不是以