flac3d 双线隧道开挖和基坑开挖。 临近既有隧道基坑开挖。 首先进行隧道开挖，考虑应力释放...
flac3d 双线隧道开挖和基坑开挖。 临近既有隧道基坑开挖。 首先进行隧道开挖，考虑应力释放，使用反力支撑法，使用shell壳单元支护。 然后进行基坑开挖，使用地连墙和对撑支护。 分三层开挖。
最近手头有个挺有意思的FLAC3D项目，既要搞双线隧道开挖，又要在既有隧道旁边挖基坑。这活儿就像在豆腐块旁边切蛋糕，既要保证新开挖的稳定，又不能影响原来的结构。今天咱们就唠唠这个项目的实现思路，顺带分享些关键代码。
先说说隧道开挖部分。这里用反力支撑法控制应力释放，相当于让岩土体自己慢慢卸力。支护方面选了shell壳单元，毕竟要模拟喷射混凝土这类薄壳结构。看这段核心代码：
;隧道开挖循环 def tunnel_excavate loop n (1,3) command model null range cyl end1 0 0 0 end2 50 0 0 rad 5 shell cmodel elastic thick 0.2 range position-z 0 10 struct node fix velocity range cyl end1 0 0 0 end2 50 0 0 rad 5.5 endcommand step 2000 endloop end tunnel_excavate这段代码有几个关键点：
- 分三步开挖（loop循环），每次挖掉指定圆柱范围的岩体
- 立马给隧道壁挂上0.2米厚的弹性壳单元（模拟喷射混凝土）
- 设置结构节点固定，相当于反力支撑的边界条件
- 每步开挖后跑2000步计算确保收敛
重点来了！开挖后的应力云图显示，最大主应力集中在隧道拱肩位置，这和现场监测数据吻合。shell单元应力输出值要控制在<2MPa，不然后续施工得调整支护参数。
接下来是紧邻隧道的基坑开挖。这里用了地连墙+对撑的组合拳，分三层开挖。看这段关键实现：
;地连墙创建 struct wall create ... prop young 30e9 poisson 0.2 ... ;基坑开挖函数 fish define dig_pit loop layer (1,3) excavate_pit = excavate_pit + 5 ;每层挖深5米 command model null range excavate_pit struct prop wall stress=0 ;激活地连墙 install strut id=layer ;安装对撑 endcommand step 3000 endloop end这里埋了个坑：地连墙的接触面刚度参数设置不对的话，墙后土体位移会异常增大。记得用struct node fix displacement锁定地连墙底端。三层开挖时，每层都要及时安装水平对撑，代码里用install strut按层号自动部署。
监测点数据用fish函数实时抓取：
fish_history tunnel_displacement = struct.node.disp(tunnel_node) fish_history pit_wall_stress = struct.wall.stress(wall_segment)项目收尾时发现个有趣现象：当基坑开挖到第二层时，既有隧道拱顶位移突然增大0.5mm。排查发现是地连墙接缝处刚度突变导致的，后来在模型里加了接触面摩擦系数就稳了。这个案例告诉我们，临近施工时结构连接处的细部处理才是真·魔鬼细节。
最后说个实用技巧：用plot bitmap命令导出云图时，记得把色标范围固定，不然不同施工阶段的云图对比会看得你怀疑人生。现场实测数据导入可以用table命令直接读取csv，比手动输数据靠谱多了。