可恶！又被COMSOL炫到了！

流固耦合

以青岛某地铁站为地质背景，在COMSOL Multiphysics上建立不同倾角断层隧道开挖流固耦合模型，并对隧道开挖卸荷诱发断层活化引起突水过程中围岩应力场、渗流场进行分析。

建立模型

首先建立固体力学模型，获得隧道开挖过程中围岩应力场演化云图，而后将所得结果导入达西流动模型中进行二次运算，得出隧道开挖过程中渗流场演化规律云图，从而建立起固流耦合模型。

1、坐标及范围

本次建模局部坐标系为solid.sx、solid.sy和solid.sz  ，隧道开挖方向为 ｚ 轴负方向，延伸长度为24ｍ ；垂直隧道开挖方向向右、向上分别为 ｘ 轴正方向、ｙ 轴正方向，延伸长度分别为72ｍ 、45ｍ 。

模型中隧道断面垂高为７ｍ ，隧道地表水压约为0.1ＭＰａ ，断层上下盘落差为1.7ｍ ，模型中断层倾角取55° 、 65° 、75°三个值。

2、网格划分

根据地质背景拱顶覆岩分布特征，采用1:1的比例建立模型，设定所建的模型均为线弹性，且是各向同性的材料。

3、测点布置

隧道开挖卸荷效应首先波及隧道围岩，使围岩应力峰值点及集中区向深部转移。为有效获取隧道开挖卸荷及断层活化突水过程中拱顶覆岩中近断层侧围岩应力场、渗流场的演化规律云图，在隧道拱顶与断层上盘相交处距拱顶以上 0ｍ 、 5ｍ 、 10ｍ 处分别布设3个应力监测点，实时采集捕捉隧道开挖推进过程中对应点位处应力、渗流等物理参量的数值。结合对测点数据的分析结果，辅助数值模拟结果的分析和表征。

4、岩层组分参数

参考工程现场隧道拱顶覆岩不同岩层岩性以及数值模拟所需相关岩石物理力学性质的相似原理，整个模型纵向上划分为４种不同的岩层组分。

模拟分析

1、围岩应力场演化规律

隧道开挖至9.6ｍ 时，断层已开始活化，与隧道围岩应力卸荷及重分布形成联动互馈；随着断层倾角的增大，断层对应力的吸收阻隔作用与断层倾角呈正相关；倾角越大，断层对开挖卸荷扰动敏感程度越低。

2、围岩渗流场演化规律

随开挖推进，覆岩渗流场呈现向隧道围岩运动的趋势，水体流速和流量加大，隧道围岩低水压区逐渐向拱底集中；断层倾角愈大，围岩破碎程度越高，断层对水的汇集导通作用越强。

投稿专栏

“COMSOL仿真交流”微信公众号自创办以来得到了广大科研工作者和研究生的广泛关注和支持。为更好地服务仿真学习与研究，本公众号现因业务需要长期招聘供稿作者。您可发送简历到likeapoem@163.com，邮件主题请注明“姓名+供稿作者”。

欢迎“光声力热”类相关专业的科研爱好者加入“COMSOL仿真交流”团队。

来源： COMSOL仿真交流