MatDEM应用：水城滑坡

通过改进MatDEM的3DSlope示例代码，基于滑坡前后地面数字高程数据，建立边坡的三维离散元模型。

滑坡前后未移动的部分设为固定单元，以减小计算量。经优化切割，最终模型包含64万单元，其中14万活动单元，单元平均半径为1.75米。数值模拟耗时22小时，模拟真实世界时间70秒。

图1 水城滑坡模拟图

7月26日无人机拍摄的贵州省六盘水市水城县鸡场镇坪地村岔沟组山体滑坡现场。

图2 水城滑坡现场      新华社 图

从下图中的Energy Conversion图中可以看到，30~40秒时动能达到最大值，相当于滑坡速度达到最大值。70秒时滑坡的动能趋于零，而重力势能还在减少，说明滑坡表面颗粒还在局部运动。

图3 滑坡过程中的能量转化图

感谢中国地震局地质研究所提供滑坡高程数据，3DSlope示例代码包含在MatDEM软件包中，具体使用可参见教学视频和专著。