Python球物理模拟

Python球物理模拟

2D杆具有3个自由度（2个速度/位置+ 1个旋转/角频率）。
我将代表中心位置，该位置已通过力以通常的方式进行了修改，并使用旋转（为简单起见，围绕系统中心）变量来计算粒子的位置。
旋转由力修改

ang_accel = F * r * sin (angle(F,r)) / (2*M * r^2)

哪里

ang_accel

是角加速度

F

是作用在特定球上的力，因此有2个扭矩*相加，因为有两个力相加（矢量方向）以更新中心的位置。

r

是长度的一半，

angle(F,r)

是力矢量和半径矢量之间的角度（从中心到受力的粒子），

因此，这

F * r * sin (angle(F,r))

是围绕中心的转矩，并且

2*M * r^2

是围绕中心的两点的系统的惯性矩。