如何在体素地形中实现应力/应变力学？

我正在开发一个类似于Minecraft的世界，该世界将地形划分为体素。但是，我也希望景观的不稳定配置可以预期地崩溃。例如，悬垂太重的物体会在“高应力”体素处断裂并折断，而具有不对称腐蚀基础的柱体也会破裂。我正在考虑向地形中的每个体素添加一个“应力矢量”字段，并执行以下操作（伪代码）：

foreach voxel in terrain
    foreach neighbor in voxel.neighbor_voxels()
        if magnitude(voxel.stress - neighbor.stress) > stressThreshold
            detach_voxels(voxel, neighbor)

但是问题是我不知道如何计算这些单独的压力。是否存在某种专门用于体素离散化的基于FEA的算法，可以用来计算每个体素的压力？还是我将这一切都弄错了，还有其他方法可以做到这一点？

您可以考虑对所有已连接体素进行DFS搜索，并确定最依赖或“分支”单个体素的体素。然后运行您的漂亮的伪代码。:)