桥梁结构震反应分析.ppt

4.2桥梁结构地震反应分析方法 结构抗震设计理论发展过程主要经历三个阶段——静力法、动力反应谱法和动态时程分析法 对上式求导可得振子相对速度与绝对加速度反应积分公式： 基于上面公式，用数值积分的方法，可得出各反应的时程曲线。以不同单自由度体系的周期为横坐标，以不同阻尼比为参数，绘出最大相对位移、最大相对速度和最大绝对加速度的谱曲线，分别称为相对位移反应谱、拟相对速度反应谱和拟加速度反应谱 （分别简称为位移反应谱、速度反应谱和加速度反应谱），并用符号记为SD、PSV和PSA 忽略小阻尼比的影响，有： 动态时程分析法 显式积分是在第i步计算中状态ti满足运动方程式的计算方法。 用显式积分法计算结构响应时，为了提高计算精度，时间间隔Δt必须十分小，否则，随着计算步数的增加，误差不断积累、出现发散的现象。 隐式积分是满足ti+Δt时刻运动方程式的计算方法。这种算法由于函数f包含未知的结构响应，因此，需要通过反复迭代的方法计算ti+Δt时刻的响应，它是一个非线性计算的过程。 常用的隐式积分法 隐式积分归结为求解非线性方程组。但是从计算效率考虑，在结构动力响应分析中一般对它的计算过程加以修正提高计算收敛性和效率，各种改进算法中Newmark β法和Wilson θ法是结构振动响应分析中最常用的两种方法。 5.2 Newmark β法 Newmark β法是一种加速度法，它是根据时间增量内假定的加速度变化规律计算结构动力响应的方法。由于时间增量内加速度变化规律的假定形式是任意的，因此Newmark β法有多种形式的计算公式。为了方便理解起见，以下通过几种特殊情况的加速度算法来介绍Newmark β法。 线性加速度法 假定时刻ti到ti+1 (=ti+Δt)之间加速度线性变化 在时刻ti+1结构振动响应应满足运动方程： 代入方程，得到ti+1时刻的加速度： 进一步计算ti+1时刻的速度、位移。 Newmark β法有很多的表示形式，也可以表示成直接计算位移的形式。与直接计算加速度响应的计算方法相比，直接计算位移响应的计算方法更加常用。关于直接计算位移的方法后面再介绍。 2. 平均加速度法 假定加速度在ti- ti+1区间内为平均值： 速度、位移为： 在时刻ti+1结构振动响应应满足运动方程： 得到ti+1时刻的加速度为： 进一步计算ti+1时刻的速度、位移。 3. Newmark β法统一的表达式 用增量形式表示 根据上式，得到速度和加速度的增量 Newmark β法的收敛性 可以证明， Newmark β法当β≥1/4时，计算是无条件收敛的。 Newmark β法是工程计算中最常用的方法， Newmark β法计算过程 5.3 wilsonθ法 同Newmark β法一样，wilsonθ法也是结构振动分析的常用计算方法之一。 wilsonθ法的基本假定是在时间间隔θ Δt（θ≧1.0）内加速度响应线形变化 。 因此，在间隔内任意时刻τ的加速度根据线性内插可以表示为： 代入结构运力方程式 得到位移后进一步得到t+τ加速度，然后内插得到t+Δt加速度，进一步计算位移和速度。 wilsonθ法的计算过程 1. 计算常数 2. 等效刚度计算 3. 对每一时间步计算等效荷载增量 4. 解方程计算位移 5.计算时刻t+Δt的响应 用增量形式表示的Wilsonθ法（推导省略） 计算公式 得到 时刻的响应， 再转变成t+Δt的响应： 直接积分法的补充说明 一种算法很难同时兼顾稳定性和精度，稳定与精度往往具有相反的倾向，稳定性好的计算方法精度相对比较差一些。一般而言，显式积分的稳定性差一些，而且时间间隔的取值对计算稳定性的影响很大。 稳定性好的计算方法并不意味可以用任意大的时间步长进行积分计算，无条件稳定的计算方法（比如Wilsonθ法）虽然不会发生数学上的发散现象，但是容易出现早期结果偏大、后期出现异常振动的计算结果，而且也有可能导致高频振动的计算结果失真的现象，因此，这种算法一般不太合适于带有冲击响应的结构计算。对于高次振动成分比较重要的计算，用非常小的时间步长按显式积分较多。 二、桥梁结构动力计算模型 在利用计算机进行桥梁结构计算分析时，均要根据结构形式、构造、受力特点、计算机容量、要求的精度等各种因素，选择既能较真实地描述结构内力变形性质，又能使用简便的力学计算模型，即结构的计算简图。 在按动态时程法分析桥梁结构地震反应时，由于计算工作量大，在尽量较真实地描述结构反应特点的前提下，宜尽可能简化。 在时程法分析中，最常用的计算模型有堆聚质量模型、脊骨