阻尼摆运动的系统动力学VENSIM模拟仿真

阻尼摆运动

以一个在粘稠液体内摆动的摆球为例,示意图如下图。

系统动力学流图如下：

具体方程如下：

(01) AV角速度=

(M质量*G重力/B阻尼)*SIN正弦

Units: **undefined**

(02) A摆角= INTEG (

AV角速度,

1)

Units: **undefined**

(03) B阻尼=

20

Units: **undefined**

(04) FINAL TIME = 400

Units: Month

The final time for the simulation.

(05) G重力=

32

Units: **undefined**

(06) INITIAL TIME = 0

Units: Month

The initial time for the simulation.

(07) M质量=

1

Units: **undefined**

(08) SAVEPER =

TIME STEP

Units: Month [0,?]

The frequency with which output is stored.

(09) SIN正弦=

sin(A摆角*6.28/360)

Units: **undefined**

(10) TIME STEP = 1

Units: Month [0,?]

The time step for the simulation.

摆角如下：

摆角趋于180度，不会发生过摆现象。

角速度会逐渐增大，到逐渐减小。

详细请参考：

王其藩，《系统动力学》，上海财经大学，pp:110-112