python 仿真积分模块_形变仿真中的积分方法-CSDN博客

形变仿真中，积分方法主要有显式积分方法(Explicit time integration scheme)、隐式时间积分方法(Implicit time integration scheme)等。

下面介绍中心插值法、Houbolt 方法、Wilson-$\theta$ 方法、以及 Newmark 方法。(这些方法应该也是包含在显式、隐式方法中的)

0、变量定义

假定 $t$ 时刻，质点的位置、速度、加速度分别为 $\boldsymbol{u}_t$ 、$\dot{\boldsymbol{u}}_t$ 、 $\ddot{\boldsymbol{u}}_t$ 。仿真过程中，时间间隔为 $$

1、中心插值法

在中心插值法中，是按中心差分将速度和加速度矢量离散化为：

\[\begin{aligned}

\dot{\boldsymbol{u}}_t = \frac{1}{2\Delta t}(\boldsymbol{u}_{t + \Delta t} - \boldsymbol{u}_{t - \Delta t}) \\

\ddot{\boldsymbol{u}}_t = \frac{1}{2\Delta t^2}(\boldsymbol{u}_{t + \Delta t} - 2\boldsymbol{u}_{t} + \boldsymbol{u}_{t - \Delta t})

\end{aligned}

在中心插值法中，将 $t$ 时刻的速度和加速度用相邻时刻的位移来表示。

2、线性加速度法和 Wilson-$\theta$ 法

线性加速度法和 Wilson-$\theta$ 法，都是属于逐步积分法。

(略)

3、Newmark 方法

Newmark 在 1959 年提出的逐步积分格式，称为 Newmark 方法。它的基本假定是

\[\begin{aligned}

\dot{\boldsymbol{u}}_{t + \Delta t} = \dot{\boldsymbol{u}}_{t} + [(1-\delta)\ddot{\boldsymbol{u}}_{t} + \delta \ddot{\boldsymbol{u}}_{t + \Delta t}] \Delta t \\

\boldsymbol{u}_{t + \Delta t} = \boldsymbol{u}_{t} + \dot{\boldsymbol{u}}_{t} \Delta t + [(-\frac{1}{2} - \alpha)\ddot{\boldsymbol{u}}_{t} + \alpha \ddot{\boldsymbol{u}}_{t + \Delta t}] \Delta t^2

\end{aligned}

其中，$\delta$ 和 $\alpha$ 是按积分的精度和稳定性要求可以调整的参数。当 $\delta = \frac{1}{2}$ ， $\alpha = \frac{1}{6}$ 时，它就是线性加速算法。所以，Newmark 方法也可以理解为线性加速法俄一个小延伸。

Newmark 法最初提出作为无条件稳定的一种积分格式是常平均加速度法，即假定从 $t$ 到 $t + \Delta t$ 时刻，加速度不变，取为常数 $\frac{1}{2}(\ddot{\boldsymbol{u}}_{t} + \ddot{\boldsymbol{u}}_{t + \Delta t})$ 。此时，取 $\delta = \frac{1}{2}$ ， $\alpha = \frac{1}{4}$ 。常平均加速度法是应用的最为广泛的逐步积分方法之一。研究表明，当 $\delta \ge 0.5$ ， $\alpha \ge 0.25(0.5 + \delta)^2$ 时，Newmark 方法是无条件稳定的。

4、Houbolt 方法

(略)

小结