Taichi模拟金属cube与凝胶的相互作用出现数值爆炸的问题

我想用Taichi模拟金属cube与长方体凝胶的相互作用。凝胶下表面固定,金属cube紧贴凝胶上表面并给一个朝凝胶运动的初速度。

通过查表定义真实符合物理的材料参数
凝胶材料参数设置为:E0, nu0 = 1e4, 0.45 (杨氏模量单位Pa)
金属材料参数设置为:E1, nu1 = 2.06e7, 0.25

发现数值爆炸。如果调小二者的杨氏模量到10^2可以正常模拟,但是这样材料就不符合物理了。请问数值爆炸是什么原因造成的呢?为什么很多Taichi代码里给的杨氏模量大多在10^3Pa这个量级?Taichi是不擅长模拟杨氏模量非常大的物体吗?谢谢