第7章逐步法——对一般动力荷载的反应资料课件

高等结构动力学第七章第七章高等结构动力学对一般动力荷载的反应逐步法高等结构动力学7.1 7.1 一般概念一般概念7.2 7.2 分段精确方法分段精确方法7.3 7.3 数值近似方法数值近似方法一般注释一般注释7.3 7.3 二阶中心差分列式二阶中心差分列式7.3 7.3 积分法积分法7.3 7.3 非线性分析的增量列式非线性分析的增量列式7.3 7.3 线加速度法步骤概要线加速度法步骤概要第七章第七章 对一般动力荷载的反应对一般动力荷载的反应逐步法逐步法高等结构动力学 分析承受任意动力荷载的线性结构，分析承受任意动力荷载的线性结构，Duhamel积积分或频域分析，提供了最方便的解法。分或频域分析，提供了最方便的解法。这两种方法的推导过程中都使用了叠加原理，只这两种方法的推导过程中都使用了叠加原理，只能适用于线性体系，即反应过程中体系的特性保持不变。能适用于线性体系，即反应过程中体系的特性保持不变。另一方面，有许多种重要的结构动力学问题，体系另一方面，有许多种重要的结构动力学问题，体系不能视作线性的。如：足以引起严重破坏的地震运动下不能视作线性的。如：足以引起严重破坏的地震运动下的建筑物反应等等。因此，还需要发展适用于非线性体的建筑物反应等等。因此，还需要发展适用于非线性体系的其它分析方法。系的其它分析方法。7.17.17.17.1 一般概念一般概念一般概念一般概念 7.17.17.17.1 一般概念一般概念一般概念一般概念高等结构动力学动力反应分析的方法动力反应分析的方法1、叠加法、叠加法线性体系，即反应过程中体系的特性保线性体系，即反应过程中体系的特性保持不变；持不变；2、逐步法、逐步法体系不能视作线性的，要发展适用于非体系不能视作线性的，要发展适用于非线性体系的方法。线性体系的方法。高等结构动力学逐步法的思想将荷载和反应历程分成一系列的时间间隔或步；将荷载和反应历程分成一系列的时间间隔或步；每步期间结构特性保持常数；每步期间结构特性保持常数；每步的反应为此步开始时的初始条件（位移及速每步的反应为此步开始时的初始条件（位移及速度）和该步期间的荷载历程引起；度）和该步期间的荷载历程引起；是一个分段线性化的系统。是一个分段线性化的系统。7.17.17.17.1 一般概念一般概念一般概念一般概念高等结构动力学近似的方法方程的近似方程求解方法的近似 1）数值微分 2）数值积分 7.17.17.17.1 一般概念一般概念一般概念一般概念高等结构动力学7.2 7.2 7.2 7.2 分段精确方法分段精确方法分段精确方法分段精确方法7.2 7.2 7.2 7.2 分段精确方法分段精确方法分段精确方法分段精确方法高等结构动力学(7-1)(7-1)(7-2)(7-2)(7-3)(7-3)(7-4)(7-4)(7-5)(7-5)7.2 7.2 7.2 7.2 分段精确方法分段精确方法分段精确方法分段精确方法高等结构动力学其中其中同样地，可获得时间步长期间的速度为同样地，可获得时间步长期间的速度为(7-6)(7-6)7.2 7.2 7.2 7.2 分段精确方法分段精确方法分段精确方法分段精确方法高等结构动力学7.2 7.2 7.2 7.2 分段精确方法分段精确方法分段精确方法分段精确方法高等结构动力学图图E7-2 分段精确计算的反应分段精确计算的反应7.2 7.2 7.2 7.2 分段精确方法分段精确方法分段精确方法分段精确方法高等结构动力学7.3 7.3 7.3 7.3 数值近似方法数值近似方法数值近似方法数值近似方法一般注释一般注释一般注释一般注释数值方法数值方法数值微分、数值积分数值微分、数值积分近似逐步法的要点：近似逐步法的要点：1 1）列式可以为显式亦可为隐式；）列式可以为显式亦可为隐式；2 2）效率是重要的，关系到达到精度的工作量，任何情）效率是重要的，关系到达到精度的工作量，任何情况下步长必须短到足以提供荷载和反应历程足够的精度；况下步长必须短到足以提供荷载和反应历程足够的精度；3 3）产生误差的技术原因）产生误差的技术原因舍入、不稳定性、截断；舍入、不稳定性、截断；4 4）产生误差的自身原因）产生误差的自身原因相位的漂移或频率的显著相位的漂移或频率的显著改变、人工阻尼。改变、人工阻尼。7.3 7.3 7.3 7.3 数值近似方法数值近似方法数值近似方法数值近似方法一般注释一般注释一般注释一般注释高等结构动力学7.4 7.4 7.4 7.4 二阶中心差分列式二阶中心差分列式二阶中心差分列式二阶中心差分列式7.4 7.4 7.4 7.4 二阶中心差分列式二阶中心差分列式二阶中心差分列式二阶中心差分列式高等结构动力学（7-77-7）（7-87-8）7.4 7.4 7.4 7.4 二阶中心差分列式二阶中心差分列式二阶中心差分列式二阶中心差分列式使用中心差分使用中心差分使用中心差分使用中心差分:高等结构动力学（7-97-9）（7-107-10）7.4 7.4 7.4 7.4 二阶中心差分列式二阶中心差分列式二阶中心差分列式二阶中心差分列式高等结构动力学（7-117-11）（7-127-12）7.4 7.4 7.4 7.4 二阶中心差分列式二阶中心差分列式二阶中心差分列式二阶中心差分列式使用平均速度的概念使用平均速度的概念使用平均速度的概念使用平均速度的概念,可以得到可以得到可以得到可以得到;高等结构动力学7.5 7.5 7.5 7.5 积分法积分法积分法积分法 另一类一般性的逐步进行动力反应分析的数值方法是，对每一时间步，从初始到最终条件应用积分前进一步。这个基本概念可用如下式子表示：（7-13a7-13a）（7-13b7-13b）7.5 7.5 7.5 7.5 积分法积分法积分法积分法高等结构动力学7.5 7.5 7.5 7.5 积分法积分法积分法积分法 最终速度和位移依据这些值的初始值加一个积分表达式。速度的变化依赖于加速度历程的积分，而位移的变化依赖于相应的速度积分。为了进行这类分析，首先需要假设在时间步的持续时间内加速度是如何变化的；加速度的假设也控制了速度的变化，因而可以由这一步向前获得下一时间步。高等结构动力学7.5 7.5 7.5 7.5 积分法积分法积分法积分法Euler-Gauss方法 最简单是假设加速度在时间步持续时间内为常数，结果是在持续时间内速度为线性，位移为二次曲线 著名的Euler-Gauss方法。高等结构动力学 列式特性的样式，假设常量加速度是由初值及步长持续时间内所获得的最终加速度的平均。在此图中也显示了速度和位移的表达式，它们是对此加速度在这步持续时间内任意时刻 由逐次积分所获得的，把 代入这些表达式而获得最终速度和位移。7.5 7.5 7.5 7.5 积分法积分法积分法积分法高等结构动力学图7-3 基于常平均加速度的运动h 加速度（常数）速度（线性）位移（二次的）高等结构动力学为了对任意步开始这种分析，首先需要计算初始加速度时刻式（7-7）所示的动力平衡表达式获得。另外，最终加速度需要应用隐式列式，它的值可以由迭代获得。对，这可以解它的值可以由迭代获得。对开始时用任意假设的值，再用图7-3所列式（a）和（b）得到 和的值。然后，用与式（7-7）相当的表达式从动力平衡方程计算时刻值的一个改进，由此再导得速度和和和位移的改进值，最后，迭代收敛于这时间步最终加速度的一个固定的值，这个过程可以前进一步到下一时间步。7.5 7.5 7.5 7.5 积分法积分法积分法积分法高等结构动力学迭代的列式开始时用任意假设的值 ，再用图7-3所列式（a）7.5 7.5 7.5 7.5 积分法积分法积分法积分法高等结构动力学 常平均加速度法的主要优点：是无条件稳定的。也就是说，从一步到下一步不管时间步长选得如何长，误差不会放大。因此时间步长的选择只需要考虑所定义动力激励和结构的振动反映特性。高等结构动力学Newmark法 一种更一般的逐步列式是由Newmark提出的，前面的方法可以作为它的特殊情况。但是也可以在其他一些形式下应用。在Newmark列式中，对最终速度和位移的基本积分式（7-13）如下所示：7.5 7.5 7.5 7.5 积分法积分法积分法积分法高等结构动力学（7-14a7-14a）（7-14b7-14b）由式（7-14a）显然可见，系数 提供了在初始和最终加速度改变影响之间的线性变化的权重，类似地，系数提供了在这些初始和最终加速度对位移改变贡献的权重。高等结构动力学 从该列式性能的研究发现，系数控制了由这个逐步法导致的人工阻尼量；如果=1/2，方法是无人工阻尼的，因此这个值被推荐用于标准的单自由度分析。在式（7-14a）和式（7-14b）中令系数=1/2和=1/4，此时可以看到，Newmark列式直接退化为图7-3所示最终速度和位移的表达式。因此，Newmark=1/4法也可以归诸于常平均加速度法。7.5 7.5 7.5 7.5 积分法积分法积分法积分法高等结构动力学 另一方面，如果取作1/6（用=1/2），最终速度和位移的表达式成为7.5 7.5 7.5 7.5 积分法积分法积分法积分法（7-15a7-15a）（7-15b7-15b）高等结构动力学 这些结果也可以如图7-4所示，由假设在时间步持续期间加速度在 和 的初始到最后值之间线性变化来得到；因此=1/6的Newmark法也称为线加速度法。像常平均加速度法一样，此法在实际中也是广为应用的。但是与=1/4方法对比，线加速度法仅是条件稳定的。可是，与二阶中心差分法一样，在单自由度体系分析中这个限制并不重要，因为要获得动力荷载和反应的满意表示，必须取比这一限制更短的时间步长。7.5 7.5 7.5 7.5 积分法积分法积分法积分法高等结构动力学变换到显式公式 法的隐式列式是不方便应用的，因为每一时间步内为了确定此步终点加速度需要进行迭代。因此，通常被修改为显式形式，目的是最终加速度用其他反应量表示选择一个基本未知量（位移较好）。7.5 7.5 7.5 7.5 积分法积分法积分法积分法高等结构动力学 再代入图7-3式（a）中，获得最终速度表达式为（7-16b7-16b）7.5 7.5 7.5 7.5 积分法积分法积分法积分法（7-16a7-16a）根据图7-3式（b）对最终加速度求解可得高等结构动力学在t1时刻写出动力平衡方程并将式（7-16a）和（7-16b）代入上式，则可导得仅含时间步终点未知位移v1的表达式。经适当归并同类项，此式可写为（7-177-17）这是一个静力平衡方程的形式，它包含等效刚度（7-17a7-17a）7.5 7.5 7.5 7.5 积分法积分法积分法积分法高等结构动力学和等效荷载（7-17b7-17b）在式（7-17）里下标c用以标记常平均加速度法。7.5 7.5 7.5 7.5 积分法积分法积分法积分法高等结构动力学而不是从式（7-16a）求，因而保留了平衡条件。7.5 7.5 7.5 7.5 积分法积分法积分法积分法使用这个显式公式，时间步终点位移v1可直接由式（7-17）计算，所要用到的仅是时间步开始时的数据。然后，此时刻的速度 可用式（7-16b）计算，最后，此时间步终点的加速度 由求解该时刻的动力平衡方程而得 高等结构动力学 采用同样的方法，使用图7-4中的式（a）和（b），也可以类似地将线性加速度法转换为显式形式，这些列式的位移差别就是等效刚度，等效荷载及最终速度的表达式不同。对线加速度分析来说，等效静力平衡方程为 （7-187-18）7.5 7.5 7.5 7.5 积分法积分法积分法积分法高等结构动力学图7-4 基于线性变化加速度的运动 加速度（线性）速度（二次的）位移（三次的）h7.5 7.5 7.5 7.5 积分法积分法积分法积分法高等结构动力学（7-18b7-18b）当位移v 1由式（7-18）计算时，同时刻的速度可由如下表达式给出相当于式（7-16b）:（7-18c7-18c）7.5 7.5 7.5 7.5 积分法积分法积分法积分法其中下标d表示线加速度法。等效刚度和荷载分别为 （7-18a7-18a）高等结构动力学 线加速度法仅仅是条件稳定，但如前面所述，对于单自由度体系分析，这一点并不重要。另一方面，假设每个步长持续时间内加速度线性变化，要比连续用常加速度法能获得真实特性的更好近似。实际上，数值实验结果也证明了线加速度法结果比用常加速度步所得结果优越。基于此理由，对单自由度体系的分析推荐使用线加速度（=1/6）法。7.5 7.5 7.5 7.5 积分法积分法积分法积分法高等结构动力学7.6 7.6 7.6 7.6 非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式上述分析中，体系的动力特性保持不变，仅可用于线性体上述分析中，体系的动力特性保持不变，仅可用于线性体系；在步长系；在步长t t足够小时，可以认为体系的动力特性是常足够小时，可以认为体系的动力特性是常数；数；采用一系列短时间增量采用一系列短时间增量t t计算反应，为了方便取计算反应，为了方便取t t为等为等步长步长；在在步长步长的起点和终点建立动力平的起点和终点建立动力平衡衡条件，并以一个条件，并以一个假设的反应机理为根据，近似地计算在时间增量范围内体假设的反应机理为根据，近似地计算在时间增量范围内体系的运动系的运动(通常忽略去在时间间隔内可能产生的不平衡通常忽略去在时间间隔内可能产生的不平衡)；体系的非线性特性可用每个时间增量起点所求得的当前变体系的非线性特性可用每个时间增量起点所求得的当前变形状态的特性来说明。利用本计算时间区间终点的速度和形状态的特性来说明。利用本计算时间区间终点的速度和位移作为下一计算时间区间的初始条件从而可得到整个反位移作为下一计算时间区间的初始条件从而可得到整个反应应；7.6 7.6 7.6 7.6 非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式高等结构动力学这个过程可以逐步地从加荷开始时起进行到任何所要求的这个过程可以逐步地从加荷开始时起进行到任何所要求的时间，而非线性特性则可用一系列相继改变的线性体系来时间，而非线性特性则可用一系列相继改变的线性体系来逼近。逼近。对于非线性分析，最有效的方法是逐步积分法对于非线性分析，最有效的方法是逐步积分法；7.6 7.6 7.6 7.6 非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式高等结构动力学图图 7-5 非线性动力体系的定义：非线性动力体系的定义：(a)基本单自由度结构；基本单自由度结构；(b)力的平衡；力的平衡；(c)非线性阻尼；非线性阻尼；(d)非线性刚度；非线性刚度；(e)作用荷载作用荷载(a)(b)(c)(d)(e)7.6 7.6 7.6 7.6 非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式高等结构动力学 考虑的结构为图考虑的结构为图7-5(a)所示的单自由度体系，体系的特性所示的单自由度体系，体系的特性m，k，c和和p(t)可以理解为可以理解为2-5节所讨论的广义量，而并不只局限于节所讨论的广义量，而并不只局限于图面上所示的简单情况。作用于体系质量图面上所示的简单情况。作用于体系质量m上的力如图上的力如图7-5(b)所所示，弹簧力和阻尼力的一般非线性性质分别绘于图示，弹簧力和阻尼力的一般非线性性质分别绘于图7-5(c)和和(d)中，中，任意作用荷载则绘于图任意作用荷载则绘于图7-5(e)中。中。在任一瞬间在任一瞬间t，作用于质量，作用于质量m上的力系的平衡要求：上的力系的平衡要求：（7-19a7-19a）而在暂短时间而在暂短时间t以后，平衡方程将为以后，平衡方程将为（7-19b7-19b）从方程从方程(7-19b)减去方程减去方程(7-19a)则可得时间间隔则可得时间间隔t的运动方程的增量形式：的运动方程的增量形式：（7-207-20）7.6 7.6 7.6 7.6 非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式高等结构动力学 这个方程中的增量力可表示如下：这个方程中的增量力可表示如下：（7-21a7-21a）其中，质量其中，质量m不言而喻被假设为常量，而不言而喻被假设为常量，而c(t)和和k(t)项则表示项则表示与时间间隔内速度和位移相应的阻尼和刚度特性，分别如图与时间间隔内速度和位移相应的阻尼和刚度特性，分别如图7-5(c)和和(d)所示。实际上，因为在时间增量末端的速度和位栘将所示。实际上，因为在时间增量末端的速度和位栘将依赖于这些特性，故所示割线斜率只能用依赖于这些特性，故所示割线斜率只能用迭迭代法进行计算。代法进行计算。因因此此，通常用时间间隔起点所定义的切线斜率来代替。，通常用时间间隔起点所定义的切线斜率来代替。（7-227-22）（7-21b7-21b）（7-21c7-21c）（7-21d7-21d）7.6 7.6 7.6 7.6 非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式高等结构动力学 将方程将方程(7-21)代入方程代入方程(7-20)，可导得时间，可导得时间t的增量平衡方程的增量平衡方程的最终形式为：的最终形式为：（7-237-23）在这种类型的分析中，所讨论的材料特性可以包括任何的非在这种类型的分析中，所讨论的材料特性可以包括任何的非线性形式。因此，没有必要规定弹簧力线性形式。因此，没有必要规定弹簧力fs像非线性弹性材那样仅像非线性弹性材那样仅仅依赖于位移。同样也可以很好地说明非线性滞变材料，在这仅依赖于位移。同样也可以很好地说明非线性滞变材料，在这种材料中，力依赖于变形的过去时程以及位移的当前值，唯一种材料中，力依赖于变形的过去时程以及位移的当前值，唯一的要求是刚度特性必须完全由变形的过去和目前状态所确定。的要求是刚度特性必须完全由变形的过去和目前状态所确定。此外，隐含的质量不变的假定显然是可以改动的：质量也可以此外，隐含的质量不变的假定显然是可以改动的：质量也可以表示为随时间变化的量。表示为随时间变化的量。7.6 7.6 7.6 7.6 非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式高等结构动力学 有很多方法可以用来进行方程有很多方法可以用来进行方程(7-23)的数值积分。的数值积分。引入假定：在每个时间增量内加速度线性变化，而引入假定：在每个时间增量内加速度线性变化，而且体系的特性在这个间隔内保持为常量。且体系的特性在这个间隔内保持为常量。质量在时间间隔内的运动用图形的形式绘于图质量在时间间隔内的运动用图形的形式绘于图7-6中，与所中，与所假定的加速度为线性变化的方程一起，图中还分别绘出了相应的假定的加速度为线性变化的方程一起，图中还分别绘出了相应的速度二次变化、位移三次变化的图形。计算在间隔终点速度二次变化、位移三次变化的图形。计算在间隔终点(t)时时后的值，导得速度和位移的增量方程如下：后的值，导得速度和位移的增量方程如下：（7-6a*7-6a*）（7-6b*7-6b*）7.6 7.6 7.6 7.6 非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式高等结构动力学图图 7-6 基于线性变化加速基于线性变化加速 度的增量运动度的增量运动7.6 7.6 7.6 7.6 非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式高等结构动力学问题：由已知的由已知的 titi时刻的反应时刻的反应 ，求，求 ti+ti+t t时刻的反应时刻的反应？增量的表达方式增量的表达方式三量相关，用哪个做基本的未知量？三量相关，用哪个做基本的未知量？7.6 7.6 7.6 7.6 非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式高等结构动力学线性加速度的假定：线性加速度的假定：在在 t=ti,ti+t=ti,ti+t t的间隔内，的间隔内，=0-=0-t t，t=ti+t=ti+7.6 7.6 7.6 7.6 非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式高等结构动力学7.6 7.6 7.6 7.6 非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式高等结构动力学 位移增量作为分析的基本变量，解式位移增量作为分析的基本变量，解式(7-6a*)可得加速度增量，并将这个可得加速度增量，并将这个表达式代入方程表达式代入方程(7-6b*)，得：，得：（7-7a*7-7a*）（7-7b*7-7b*）获得了用位移增量表达的速度和加速度增量。获得了用位移增量表达的速度和加速度增量。把方程把方程(7-7*)代入方程代入方程(7-23),可导得下列的运动方程：可导得下列的运动方程：（7-8*7-8*）7.6 7.6 7.6 7.6 非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式高等结构动力学最终把所有含已知初始条件的各项移到右边，给出最终把所有含已知初始条件的各项移到右边，给出（7-9a*7-9a*）（7-9b*7-9b*）方程方程(7-8*)相当于静力增量平衡关系，即可将荷载增量除以刚相当于静力增量平衡关系，即可将荷载增量除以刚度求得位移增量。在等效荷载和刚度项中包含惯性和阻尼的项反度求得位移增量。在等效荷载和刚度项中包含惯性和阻尼的项反映了动力特性。解方程（映了动力特性。解方程（7-8*），得出位移增量后，将此值代入），得出位移增量后，将此值代入方程方程(7-7b*)即可获得速度增量。即可获得速度增量。下一时段的初始条件由该时段起点速度和位移值加上这些增下一时段的初始条件由该时段起点速度和位移值加上这些增量值得到。量值得到。（7-8*7-8*）其中：其中：7.6 7.6 7.6 7.6 非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式高等结构动力学 分析包含两个重要假定：分析包含两个重要假定：(1)加速度为线性变化；加速度为线性变化；(2)质量、阻尼和刚度特性在时间步长内保持常量。质量、阻尼和刚度特性在时间步长内保持常量。虽然时间步长很短时误差甚小，但这两个假定都不是完全正虽然时间步长很短时误差甚小，但这两个假定都不是完全正确的。因此，误差在增量平衡关系中出现，并积累。为了避免误确的。因此，误差在增量平衡关系中出现，并积累。为了避免误差的积累，在分析的每一步中利用总的平衡条件消除误差。差的积累，在分析的每一步中利用总的平衡条件消除误差。只要从总外荷载中减去总阻尼力和弹性力以表示时间步长起只要从总外荷载中减去总阻尼力和弹性力以表示时间步长起点的加速度就行了。点的加速度就行了。7.6 7.6 7.6 7.6 非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式非线性分析的增量列式高等结构动力学7.7 7.7 7.7 7.7 线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要 对任一给定的时间增量，按如下程序进行运算：对任一给定的时间增量，按如下程序进行运算：(1)初始速度和位移值是已知的。它们或是前一增量的终点初始速度和位移值是已知的。它们或是前一增量的终点值或是问题的初始条件值；值或是问题的初始条件值；(2)利用这些值及结构特定的非线性特性，可找出时间间隔利用这些值及结构特定的非线性特性，可找出时间间隔内的阻尼内的阻尼c(t)、刚度、刚度k(t)以及阻尼力以及阻尼力fD(t)和弹性力和弹性力fS(t)的当前值，的当前值，例如图例如图7-5(c)和和7-5(d)中所示；中所示；(3)初始加速度由下式给出：初始加速度由下式给出：（7-257-25）时间时间t时平衡方程的重新排列；时平衡方程的重新排列；(4)等效荷载增量等效荷载增量 和等效刚度和等效刚度 按方程按方程(7-9*)计算；计算；7.7 7.7 7.7 7.7 线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要高等结构动力学 (5)由方程（由方程（7-8*）可求得位移增量，从而由方程）可求得位移增量，从而由方程(7-7b*)可求可求得速度增量；得速度增量；(6)最后，获得时段终点速度和位移：最后，获得时段终点速度和位移：（7-26a7-26a）（7-26b7-26b）7.7 7.7 7.7 7.7 线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要高等结构动力学 当第当第6步运算完成时，这个时段的分析结束；下一个时段的步运算完成时，这个时段的分析结束；下一个时段的分析，只需将上述整个程序重复进行即可。显然，连续进行上述分析，只需将上述整个程序重复进行即可。显然，连续进行上述运算，就可进行任意多个时间增量的分析，这样就能算出具有任运算，就可进行任意多个时间增量的分析，这样就能算出具有任何非线性性质的单自由度体系的全部反应时程。线性体系自然也何非线性性质的单自由度体系的全部反应时程。线性体系自然也可以用同样的方法进行处理。此时阻尼和刚度特性保持不变，因可以用同样的方法进行处理。此时阻尼和刚度特性保持不变，因此分析的过程将稍为简单一些。此分析的过程将稍为简单一些。7.7 7.7 7.7 7.7 线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要高等结构动力学 逐步积分法的精度，依赖于时间增量逐步积分法的精度，依赖于时间增量t的长的长度，选取时注意三个因素：度，选取时注意三个因素：(1)作用荷载作用荷载p(t)的变化速率；的变化速率；(2)非线性阻尼和刚度特性的复杂性；非线性阻尼和刚度特性的复杂性；(3)结构振动周期结构振动周期T。7.7 7.7 7.7 7.7 线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要高等结构动力学 为可靠地反映这些因素，时间增量必须足为可靠地反映这些因素，时间增量必须足够的短，其中最后一个因素是和体系的自由振动够的短，其中最后一个因素是和体系的自由振动特性联系在一起的。一般来说，材料特性的变化特性联系在一起的。一般来说，材料特性的变化不是关键性的因素，如果发生一个重大的突然变不是关键性的因素，如果发生一个重大的突然变化，例如一个弹塑性弹簧的屈服，这时可以引入化，例如一个弹塑性弹簧的屈服，这时可以引入一个特殊再细分的时间增量来精确地处理这个影一个特殊再细分的时间增量来精确地处理这个影响。同时，要估算能恰当地模拟动力荷载主要形响。同时，要估算能恰当地模拟动力荷载主要形状的时间增量也并不困难。状的时间增量也并不困难。7.7 7.7 7.7 7.7 线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要高等结构动力学 因此，如果加载过程比较简单，时间间隔的选取主要依赖于因此，如果加载过程比较简单，时间间隔的选取主要依赖于结构的振动周期。这个线性加速度方法只是有条件稳定的，如果结构的振动周期。这个线性加速度方法只是有条件稳定的，如果时间增量大于振动周期的一半左右，则将给出扩散的解。但是，时间增量大于振动周期的一半左右，则将给出扩散的解。但是，为了提供适当的精度，时间增量必须比这个短的多，因此不稳定为了提供适当的精度，时间增量必须比这个短的多，因此不稳定是不会成为问题的。一般来说，按照经验，如果增量是不会成为问题的。一般来说，按照经验，如果增量一一周期比周期比t/T1/10，则可获得可靠的结果。如果对于得出的解有任则可获得可靠的结果。如果对于得出的解有任何怀疑，则第二次分析时取时间增值的一半进行，何怀疑，则第二次分析时取时间增值的一半进行，如果在第如果在第二次分析中反应没有明显的变化，二次分析中反应没有明显的变化，则可以认为则可以认为数值积分所产生的误差是可以忽略不计的。数值积分所产生的误差是可以忽略不计的。7.7 7.7 7.7 7.7 线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要高等结构动力学例题例题 E7-2 为了说明应用上述线性加速度逐步积分法的手算方法为了说明应用上述线性加速度逐步积分法的手算方法,我们我们来计算图来计算图E7-3所给出的弹塑性单自由度刚架在所给的加荷过程下的反应所给出的弹塑性单自由度刚架在所给的加荷过程下的反应.在分析中在分析中,时间步长取为时间步长取为0.1秒秒,它比要得到精度较高的结果所需的步长较它比要得到精度较高的结果所需的步长较长一些长一些,但为了说明这个手算方法但为了说明这个手算方法,这个步长已经足够了这个步长已经足够了.在这个结构里在这个结构里,假定阻尼系数保持常量假定阻尼系数保持常量,因此非线性影响仅由于发因此非线性影响仅由于发生屈服时刚度将改变才发生的生屈服时刚度将改变才发生的.此时等效刚度可以被表达为此时等效刚度可以被表达为参见方程参见方程(8-9c)其中其中k(t)或是或是5千磅千磅/英寸或是零英寸或是零,它们分别相应于刚架的弹性和屈服阶段它们分别相应于刚架的弹性和屈服阶段,等效增量荷载由下式给出等效增量荷载由下式给出参见方程参见方程(7-9b*):由方程由方程(7-7b*)给出的速度增量成为给出的速度增量成为:7.7 7.7 7.7 7.7 线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要高等结构动力学便于手算的表格示于表便于手算的表格示于表E7-1图图 E7-3 弹塑性刚架和动力荷载弹塑性刚架和动力荷载7.7 7.7 7.7 7.7 线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要高等结构动力学7.7 7.7 7.7 7.7 线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要高等结构动力学手算方法演示手算方法演示:首先,在t=0时刻(第1步),初始荷载对应的阻尼系数c=0.2千磅秒/英寸（恒为常量），刚度k(0)=5千磅/英寸，7.7 7.7 7.7 7.7 线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要高等结构动力学t=0.1(20)时刻,荷载7.7 7.7 7.7 7.7 线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要高等结构动力学如此重复进行,直至加载结束.当计算至t=0.3时刻,7.7 7.7 7.7 7.7 线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要高等结构动力学此时刚架发生屈服(v12英寸),故刚度k（0.3）=0。当t=0.7秒时，后面的计算重复前述步骤即可。7.7 7.7 7.7 7.7 线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要高等结构动力学图图 E7-4 弹塑性和弹性反应的比较（图弹塑性和弹性反应的比较（图E7-3的刚架）的刚架）7.7 7.7 7.7 7.7 线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要线性加速度法步骤概要高等结构动力学人有了知识，就会具备各种分析能力，人有了知识，就会具备各种分析能力，明辨是非的能力。明辨是非的能力。所以我们要勤恳读书，广泛阅读，所以我们要勤恳读书，广泛阅读，古人说古人说“书中自有黄金屋。书中自有黄金屋。”通过阅读科技书籍，我们能丰富知识，通过阅读科技书籍，我们能丰富知识，培养逻辑思维能力；培养逻辑思维能力；通过阅读文学作品，我们能提高文学鉴赏水平，通过阅读文学作品，我们能提高文学鉴赏水平，培养文学情趣；培养文学情趣；通过阅读报刊，我们能增长见识，扩大自己的知识面。通过阅读报刊，我们能增长见识，扩大自己的知识面。有许多书籍还能培养我们的道德情操，有许多书籍还能培养我们的道德情操，给我们巨大的精神力量，给我们巨大的精神力量，鼓舞我们前进鼓舞我们前进。