时域分析法直接积分法.pptVIP

代入结构运力方程式 得到计算 时刻位移的方程组 得到位移后进一步得到t+τ加速度，然后内插得到t+Δt加速度，进一步计算位移和速度。 当θ≧1.37时Wilsonθ法为无条件稳定的计算方法 wilsonθ法的计算过程 1. 计算常数 2. 等效刚度计算 3. 对每一时间步计算等效荷载增量 4. 解方程计算位移 5.计算时刻t+Δt的响应 用增量形式表示的Wilsonθ法（推导省略） 计算公式 增量方程 得到 时刻的响应， 再转变成t+Δt的响应： 直接积分法的补充说明 一种算法很难同时兼顾稳定性和精度，稳定与精度往往具有相反的倾向，稳定性好的计算方法精度相对比较差一些。一般而言，显式积分的稳定性差一些，而且时间间隔的取值对计算稳定性的影响很大。 稳定性好的计算方法并不意味可以用任意大的时间步长进行积分计算，无条件稳定的计算方法（比如Wilsonθ法）虽然不会发生数学上的发散现象，但是容易出现早期结果偏大、后期出现异常振动的计算结果，而且也有可能导致高频振动的计算结果失真的现象，因此，这种算法一般不太合适于带有冲击响应的结构计算。对于高次振动成分比较重要的计算，用非常小的时间步长按显式积分较多。 W=15kN x(t) 例： 求位移时程曲线， 恢复力时程曲线， 最大位移，最大恢复力， 开始时静止。 t(s) 0.1 0.8 2.5 4 3.5 2.5 1.5 1 0.5 fs 3kN 0.05 fD 0.05 解： 1.确定步长 计算步骤： 1.确定积分步长?t 2.确定当前积分步长内结构的 质量，刚度和阻尼矩阵以及 阻尼力和恢复力 3.计算初始加速度 4.确定等效刚度K*和等效荷载 矩阵P* 5.计算 6.根据公式(11)计算结构的位移， 速度和加速度响应 P(t)(kN) t(s) 0.1 0.8 2.5 4 3.5 2.5 1.5 1 0.5 fs x(m) 3kN 0.05 1. t=0 2. t=0.1s 弹性阶段 计算步骤： 1.确定积分步长?t 2.确定当前积分步长内结构的 质量，刚度和阻尼矩阵以及 阻尼力和恢复力 3.计算初始加速度 4.确定等效刚度K*和等效荷载 矩阵P* 5.计算 6.根据公式(11)计算结构的位移， 速度和加速度响应 P(t)(kN) t(s) 0.1 0.8 2.5 4 3.5 2.5 1.5 1 0.5 fs x(m) 3kN 0.05 2. t=0.1s 弹性阶段 计算步骤： 1.确定积分步长?t 2.确定当前积分步长内结构的 质量，刚度和阻尼矩阵以及 阻尼力和恢复力 3.计算初始加速度 4.确定等效刚度K*和等效荷载 矩阵P* 5.计算 6.根据公式(11)计算结构的位移， 速度和加速度响应 P(t)(kN) t(s) 0.1 0.8 2.5 4 3.5 2.5 1.5 1 0.5 fs x(m) 3kN 0.05 3. t=0.2s 4. t=0.3s 5. t=0.4s 屈服 计算步骤： 1.确定积分步长?t 2.确定当前积分步长内结构的 质量，刚度和阻尼矩阵以及 阻尼力和恢复力 3.计算初始加速度 4.确定等效刚度K*和等效荷载 矩阵P* 5.计算 6.根据公式(11)计算结构的位移， 速度和加速度响应 6. t=0.5s 7. t=0.6s 8. t=0.7s 塑性发展 塑性发展 位移减少，恢复弹性 P(t)(kN) t(s) 0.1 0.8 2.5 4 3.5 2.5 1.5 1 0.5 fs x(m) 3kN 0.05 计算步骤： 1.确定积分步长?t 2.确定当前积分步长内结构的 质量，刚度和阻尼矩阵以及 阻尼力和恢复力 3.计算初始加速度 4.确定等效刚度K*和等效荷载 矩阵P* 5.计算 6.根据公式(11)计算结构的位移， 速度和加速度响应 8. t=0.7s 位移减少，恢复弹性 fs x(m) 3kN 0.05 x(0.7) x(0.7) 9. t=0.8s … … 计算步骤： 1.确定积分步长?t 2.确定当前积分步长内结构的 质量，刚度和阻尼矩阵以及 阻尼力和恢复力 3.计算初始加速度 4.确定等效刚度K*和等效荷载 矩阵P* 5.计算 6.根据公式(11)计算结构的位移， 速度和加速度响应 t(s) 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 fs(kN) 0 0.1488 1.069 2.86 3 3 3 2.934 … … … … … x(cm) 0 0.248 1.78 4.77 7.91 10.19 11.43 (max) 11.31 … … … … … x(t)(cm) t(s) 0.3 0.6 0.9 11.43 3 fs(t)(kN) t