MatDEM流固耦合：岩浆侵入数值模拟

基于user_HydraulicBlock建立二维岩浆侵入模型。模型包含3.1万个固体单元，4.7万流体单元。中下部建立一个软弱夹层，其同时具有较高的渗透率。下部挖取一个圆形区域，作为“岩浆房”。

初始时，单元相互胶结，设定胶结断裂后的渗透率增长为100倍。岩层上部施加恒定5MPa压力，在“岩浆房”中施加恒定100MPa的压力。模拟裂隙形成和流体运移。采用笔记本电脑进行计算，模拟耗时2小时，模拟真实世界时间26.4秒。

孔隙压力场动图（红色代表较高压力）

速度场动图（红色代表较高单元速度，尖端破裂时有较高速度）

位移场动图（红色代表较高位移）

模拟多次喷出过程