15第13章_线性动力_李永双试题.doc

第13章 线性动力分析 上一章我们已经谈到使用反应谱分析方法计算结构地震作用具有很多优点，但它本质上仍然是一种拟动力方法，并不是结构真实的动力响应分析，只是对于结构动力响应最大值进行估算的近似方法。但是如果仅限于线弹性范围内，反应谱分析法仍然是被认为高效而且合理的方法。 时程分析是结构更加真实的动力分析，而且随着计算机速度的提高，最初时程分析在计算时间方面的困难已经被解决。另外，最主要的是结构设计领域的不断发展，结构分析早已经超出了线弹性的范围，面对更多非线性问题，反应谱法动力分析是没有任何意义的，这时就更需要动力时程分析的发展和应用。 SAP2000为我们提供了线性时程分析和非线性时程分析功能，而且已经可以使用线性时程分析结果直接对结构构件进行设计。本章主要向大家介绍线性时程分析的内容。 13.1 线性时程分析的基本理论 线性时程分析本质上仍然是通过对结构基本动力微分方程（12.2）的求解，来得到结构在动力荷载作用下结构的基本响应大小。与反应谱分析方法不同的是，时程分析方法是将动力作用以时间函数的形式引入微分方程，并通过相应的积分方式得到结构每一个时刻的响应以及响应的变化情况。因此时程分析的积分方式是此类型分析的关键问题之一。 在SAP2000程序中，时程分析可以使用的积分方式为模态积分和直接积分。模态积分方式是以结构的模态分析结果为基础，通过结构不同模态的积分求解来得到结构总的响应值，在程序中如果选择使用模态积分求解方式，需要选择积分求解所基于的模态分析工况，因此在定义该工况前需要先定义相关模态分析工况。一般情况下，如果结构需要考虑高阶振型的影响，就应该在该模态分析工况中考虑更多的振型数目。SAP2000中所使用的模态积分方法又被称为FNA（快速非线性分析方法），我们将在后续章节非线性分析中介绍这一内容，本章主要介绍直接积分方法。 13.1.1时间积分方式 直接积分法又被称为逐步求解法，它是最常规的动力分析求解方法，其本质是在一系列时间间隔范围内求解平衡方程。在SAP2000中时间积分方法的选择是通过其它参数区域时间积分参数的定义来完成的，这里提供了Newmark法、Wilson法、排列法、Hiber-Huges-Taytor法和Chung and Hulbert法等方法可供选择，这些方法都是隐式方法，并且大部分都是无条件稳定的方法。正如Wilson教授所提出的：“只有单步的、隐式的、无条件稳定的方法才应该用于实际结构的逐步地震分析。” 由于本书的篇幅所限，本节只介绍在SAP2000中推荐使用的几种直接积分方法。 1. Newmark法 Nemark法是由Nemark在1959年提出的，并且在随后的发展过程中，许多其他研究者对其也进行了修改和改进。为了说明Newmark积分方法，考虑动力微分方程（11.2）式写成如下的形式： (13.1) 直接使用泰勒级数获得下列两个额外方程： (13.2a) (13.2b) Newmark截短了这些方程并以下列形式来表达这些方程： (13.3c) (13.3d) 如果假定加速度在时间步长内是线性的，可得下列方程： (13.4) 将方程(13.4)代入方程(13.3c和13.3d)，可产生标准形式的Newmark方程： (13.5a) (13.5b) Newmark为结构系统中每个位移DOF的每个时间步长，通过迭代来求解方程 (13.5a、13.5b和13.1)。通过用DOF相关的质量除以方程从方程(13.1)中获得项。 在SAP2000现在的版本中所使用的Newmark方法并不是上面所述，而是1962年Wilson改进的Newmark方法。Wilson将Newmark法写成矩阵形式，并添加了刚度与质量比例阻尼，同时通过在每个时间步长对方程直接求解来消除了对迭代的需要。这要求方程 (13.5a和 13.5b) 以下列形式重写： (13.6a) (13.6b) 其中常量到为：；；；；；。 将方程(13.6a和13.6b)代入方程(13.1)，可在时间“t”处按未知节点位移ut的形式写出系统动力平衡方程。即： (13.7) 2. Wilson法 Newmark方法本身并不是无条件稳定的，1973年Wilson在常规的Newmark方法基础上通过引入一个系数θ达到了无条件的稳定，这就是现在SAP2000程序中所使用的Wilson积分方法。 系数θ的目的是通过修改时间步长，来消除真实解附近不稳定解的振荡。如果θ等于1.0，就是Newmark方法。Wilson方法先将系数θ引入到时间步长公式： (13.8a) 以及一个由以下定义的荷载：