「Python与地震工程」单自由度体系求解之Newmark-β法

「Python与地震工程」单自由度体系求解之Newmark-β法

原理

Newmark-β法是地震工程领域最经典的逐步积分算法。

推导过程请查阅结构动力学或地震工程学教材，此处仅简单列出逐步递推公式。

已知第 $i$ 步响应，则第 $i+1$ 步位移响应可按下式计算：

$$u_{i+1}=\frac{{\hat{p}}_{i+1}}{\hat{k}}$$

其中

$$\hat{k}=k+\frac{\gamma }{\beta \Delta t}c+\frac{1}{\beta {\left(\Delta t\right)}^2}m$$

p ^ i + 1 = p i + 1 + [ 1 β ( Δ t ) 2 m + γ β Δ t c ] u i + [ 1 β Δ t m + ( γ β − 1 ) c ] u