fs

fD
f D / x c tg
[C ], [ K ] 为时变矩阵 非线性问题：
x (t )

x (t )
fD
fs
x (t )
x (t )
2.增量平衡方程
f I (t ) f D (t ) fs (t ) P(t )
t t 时刻：
（1）
f I (t t ) f D (t t ) f s (t t ) P(t t ) f I (t ) f D (t ) f s (t ) P(t )
fI (t ) f I (t t ) f I (t ) [M ]x(t t ) x(t ) [M ]x(t )
f D (t ) f D (t t ) f D (t ) [C (t )]x(t )
----增量方程(3)
•
•
显式方法(explicit)在方程求解过程中只涉及到历史的第i和i－1步的信 息，而当前的第i＋1步的信息（比如空间上的其他点）不会涉及到,而 隐式方法(implicit)在求解当前点(第i＋1步)时，会涉及到其他已知点 的第i+1步信息,所以需要迭代。
显式求解与隐式在数学上说主要是在求解的递推公式一个是用显式方 程表示，一个是用影视方程来表示。比如a(n)=a(n-1)+b(n-1),后一次 迭代可以由前一次直接求解，这就是显示方程，如果a(n)=a(n1)+f[a(n)]，f[a(n))为a(n)的函数，此时a(n)不能用方程显示表示， 及数学上的隐函数，一般很难直接求解，多用迭代试算法间接求解。
f D (t ) f D (t t ) f D (t ) [C (t )]x(t )
df c t D dx t
f D (t t )
fD
斜率c(t )
f D (t )
x
x (t )
f D
df D x dx
f s (t ) f s (t t ) f s (t ) [ K (t )]x(t )
4.2 时域分析法——直接积分法 （时程分析）
4.2.1 概述
• 数值积分法是根据已知的初始时刻的位移、速度、加速
度和荷载条件，计算下一时刻振动响应的方法。 • 数值积分法属 于初值问题， 也称步步积分 法或时程分析 法。
u
ui+1 ui
O
ti
t
ti 1
时程分析法基本概念
因此，时程分析法是对结构物的运动微分方程直接进行 逐步积分求解的一种动力分析方法。由时程分析可得到各 质点随时间变化的位移、速度和加速度动力反应，并进而 可计算出构件内力的时程变化关系。由于此法是对运动方 程直接求解，又称直接动力分析法。 直接动力分析包括确定性动力分析与非确定性动力分 析两大类，即确定性动力分析中的时程分析法与非确定性 分析的随机振动分析法，这里主要介绍时程分析法。 《抗震规范》规定，重要的工程结构，例如：大跨桥 梁，特别不规则建筑、甲类建筑，高度超出规定范围的高 层建筑应采用时程分析法进行补充计算。
显式积分、隐式积分
• 显式积分是在第i步计算中状态ti满足运动方程式的
计算方法。
ui1 ui Δ t f ti , ui
•当前时刻的位移只与前一时刻的加速度和位移有关，这就意味着当前 时刻的位移求解无需迭代过程。 另外，只要将运动方程中的质量矩阵 和阻尼矩阵对角化（线性无关），前一时刻的加速度求解无需解联立 方程组，从而使问题大大简化，这就是所谓的显式求解法。
df k t s dx t
x (t ) x(t t )
df s x dx
f s 斜率k (t )
f s (t t ) f s (t )
x
x (t )
f s
结构在t时刻的刚度矩阵 由t时刻结构各构件的切线刚度确定
x(t t )
x (t )
[M ]x(t ) [C(t )]x(t ) [ K (t )]x(t ) P(t )
1.运动方程
[M ]x(t ) [C(t )]x(t ) [ K (t )]x(t ) [ M ]1 xg
f I (t ) f D (t ) fs (t ) P(t )
[C ], [ K ] 为常数矩阵 线性问题：
（1 ）
f s / x k tg
（2 ）
将(1),(2)两式相减：
令
x t x t t x t x t x t t x t
{x t } {x(t t )} {x(t )}
P(t ) P(t t ) P(t )
df c t D dx
f D (t t )
fD
斜率c(t )
f D (t )
x
x (t )
f D
df D x dx
x (t ) x(t t )
f I (t ) f D (t ) f s (t ) P(t ) fI (t ) f I (t t ) f I (t ) [M ]x(t t ) x(t ) [M ]x(t )
• 用显式积分法计算结构响应时，为了提高计算精
度，时间间隔Δt必须十分小，否则，随着计算步 数的增加，误差不断积累、出现发散的现象。
显式积分、隐式积分
• 隐式积分是满足ห้องสมุดไป่ตู้i+Δt时刻运动方程式的计算方法。
这种算法由于函数f包含未知的结构响应（求解 ti+Δt时刻的反应，但函数f却含有待求的ti+Δt时刻 的结构响应），因此，需要通过反复迭代的方法 计算ti+Δt时刻的响应，它是一个非线性计算的过 程。隐式积分归结为求解非线性方程组。