FLAC3D 6.0命令流实现隧道台阶法开挖模拟实战解析

用 FLAC3D 6.0 做隧道开挖模拟，绕不开两个难题：一是新版命令体系跟老版本差异很大，二是台阶法这种“边挖边支、分步推进”的施工工序，怎么用命令流正确表达出来。很多人第一次用命令流做台阶法隧道施工，算出来的位移不是偏大就是偏小，问题往往不在力学参数，而在开挖顺序、支护时机和地应力平衡这三件事没理顺。

这篇分享就是围绕这套流程来做梳理：从模型分组设计、初始地应力平衡，到上台阶开挖、初期支护、下台阶跟进，再到结果提取和常见坑。内容以 FLAC3D 6.0 的命令流为主线，适合正在用命令流做隧道施工模拟、或者打算从全断面开挖转向分步工序模拟的朋友参考。我尽量把每一步前后逻辑说透，而不是只丢一串代码。

1. 台阶法模拟前，先理清FLAC3D 6.0的命令流骨架

1.1 我为什么坚持全命令流做隧道施工模拟

我在FLAC3D里做隧道模拟，基本不靠GUI拉模型。不是GUI不好用，而是隧道施工模拟有个特点：一个完整项目要跑几十个工况，台阶位置、进尺长度、支护参数经常要调。用命令流把模型和工序组织好，改一个参数、调一个进尺长度，重新跑一遍就行。尤其是台阶法，多了上台阶、下台阶、核心土的先后关系，如果不用分组和脚本，纯靠GUI点鼠标，工作量大到没法维护。

我自己的习惯是：先画一张施工流程图，把每一个循环的“挖哪块、支哪块、算多少步”写清楚，再对应写命令流。这张图不画，直接抄网上的命令片段，最后必然乱。FLAC3D 6.0的命令流本质上就是按施工顺序执行的，先建模、再平衡地应力、然后分步开挖和支护，最后求解和提取数据。你把它当一段“施工日志”来写，而不是当代码堆，思路会顺很多。

1.2 6.0命令体系里和台阶法模拟强相关的几组关键词

FLAC3D从5.0版本开始命令体系大改，6.0延续了这套体系。如果你是从4.0甚至更早版本过来的，第一件事就是别拿旧示例硬套。拿台阶法模拟来说，核心命令其实就那么几组：

• model new / model title：新建工程和命名。
• zone create / zone cmodel assign / zone property：建网格、定本构、给参数。
• zone group：给区域分组，这是台阶法模拟的生命线。
• zone face apply / zone gridpoint fix：施加边界条件。
• zone cmodel assign null：把指定区域“挖空”，模拟开挖。
• sel shell create / sel cable create：施作初期支护和锚杆。
• model solve / model step / model save：求解控制。

我的理解是，这些命令可以映射到施工动作上：zone group 相当于“备料”，zone cmodel assign null 相当于“开挖”，sel shell 相当于“喷混凝土”，model step 相当于“等围岩稳定一段时间”。命令流不是把命令堆在一起，而是把工序还原成显式求解步骤。只要想明白每个施工动作对应哪个命令，剩下的就是排序问题。

2. 建模阶段的分组设计：台阶法能不能算准，从这一步就决定了

2.1 模型尺寸怎么取：边界效应与计算量的平衡

模型范围取多少，直接影响计算精度和时间。我这次按一个典型的双车道隧道来举例：隧道跨度12m，高10m，埋深20m。模型沿隧道轴向取60m，横向从隧道中线向两侧各取24m，底部取到洞底以下约18m，顶部以地表为界。整体尺寸大致满足3到5倍洞径的影响范围，边界效应可以接受。

网格方面，隧道周边加密到0.5到1m，远处适当放宽。台阶法后处理要看的是围岩扰动和支护受力，开挖区附近的网格太乱，云图会全是毛刺，结果很难说服自己。FLAC3D 6.0虽然自带网格生成工具，但复杂断面我更推荐外部网格软件建模后导入。本文为了讲命令流逻辑，先用规则矩形断面近似台阶形洞室，实际工程中这种做法也够用，关键是把分组方式理解透。

2.2 把隧道断面拆成“可开挖单元”：Group的命名和划分

台阶法模拟的第一步其实是分组设计。很多人建模时不分上下台阶，等计算时再临时选range，结果常常选错区域或者漏掉核心土。我的做法是：模型生成后，立刻把整个开挖范围按施工顺序切成不同的group。

以双台阶预留核心土为例，我把上台阶、下台阶、核心土分别设组，再按沿隧道轴向的3m进尺分段。隧道轴向为x，洞轴线y=12m，隧道轮廓简化为矩形：上台阶z=14到20m，下台阶z=10到14m，跨度y=6到18m。第一循环进尺取x=33到36m。对应命令写成：

; 作用：按开挖顺序对模型分组 zone group 'upper_1' range position-x 33 36 position-y 6 18 position-z 14 20 zone group 'lower_1' range position-x 33 36 position-y 6 18 position-z 10 14 zone group 'core_1' range position-x 33 36 position-y 8 16 position-z 10 14

组名规则我建议按“部位_循环号”来。比如 upper_1 表示上台阶第1循环，lower_2 表示下台阶第2循环，core_1 表示核心土第1循环。后面写循环开挖命令时，只需要不断替换循环号，脚本结构完全一致。

2.3 边界条件、重力与初始地应力三步法

边界条件这块基本是流水作业：四周法向约束、底部固定、顶部自由。FLAC3D 6.0里可以用zone face apply或者zone gridpoint fix。我习惯把x和y方向四周全部锁法向速度，底部锁z向，顶部不做约束。

初始地应力我推荐用三步法，很多新手容易在这一步埋雷。第一步，全局用弹性模型，只加重力sovle；第二步，平衡后把位移和速度清零；第三步，把围岩改成摩尔库伦模型，重新solve。为什么这么做？因为如果直接给摩尔库伦模型，初始地应力场不容易稳定，尤其是节理裂隙发育的岩体，一上来就大概率算不到平衡。

; 作用：边界约束+重力平衡 zone gridpoint fix velocity-z range position-z 0 zone face apply velocity-normal 0 range position-x 0 zone face apply velocity-normal 0 range position-x 64 zone face apply velocity-normal 0 range position-y 0 zone face apply velocity-normal 0 range position-y 24 model gravity 0 0 -9.81 zone cmodel assign elastic model solve ; 位移清零，准备进入弹塑性计算 zone gridpoint initialize displacement zone gridpoint initialize velocity zone cmodel assign mohr-coulomb range group '围岩' zone property density 2200 bulk 3e8 shear 1.5e8 friction 28 cohesion 0.15e6 tension 0 model solve

如果埋深更深或者有构造应力，可以用更精细的侧压力系数初始化方法，比如按K0系数给定水平应力。但对常规隧道模拟，自重平衡三步法简单可靠，现场也够用。

2.4 网格质量检查：face skin 和模型体检

网格建好之后别急着开挖。我会先执行zone face skin，把外表面提取出来，后续做初期支护壳单元时需要用到开挖临空面的surface信息。如果用的是外部网格导入，建议在建模软件里就把surface面处理好，否则到了工序阶段再补，非常被动。

另外跑计算前必须做一遍“模型体检”。用list zone group查看分组数量是否和设计一致，用plot zone查看网格有没有飞点、重叠。这一步看着不起眼，但能省下后面排查问题的几个小时。我遇到过几次，因为网格局部没连上，shell单元生成后和zone脱开，结果位移全乱。

3. 台阶法施工工序的命令流编排：从开挖到初支再到仰拱

3.1 开挖一循环：null本构与应力释放的先后手

台阶法施工和全断面最大的区别，就是“开挖—支护—再开挖”的节奏。比如上台阶第1循环开挖，我先把对应的 group 设为 null 本构，相当于把这一块岩体从计算中移除。然后不要急着加支护，先跑几百步，让围岩应力释放。喷混凝土从开挖到闭合有滞后，这个滞后反映在数值里，就是无支护状态的少量计算步。

; 上台阶第1循环开挖：先挖空，再让应力释放 zone cmodel assign null range group 'upper_1' model step 300 ; 释放步之后，施作上台阶初期支护 sel shell create surface 'tunnel_face_1' thickness 0.25 sel property young 23e9 poisson 0.2 density 2400

这里再强调一下顺序：先null、再step、最后shell。顺序一旦颠倒，等于把支护提前到变形之前，结果会偏刚，拱顶沉降明显偏小，造成假安全。释放步数取多少合适？我一般按工程经验控制30%到60%的应力释放率，先试算一个循环，对比监测点位移，再定参数。

3.2 初期支护的命令处理：shell模拟喷射混凝土、cable模拟锚杆

喷射混凝土我用shell单元，厚度按实际初支厚度给，比如25cm。FLAC3D 6.0中shell创建前要确保开挖面surface已经存在，毕竟壳单元是依附在临空面上的。锚杆用cable单元，长度3.5m，按环向间距和纵向间距折算布置；循环进尺内我一般布置4到6根简化锚杆，重点观察轴力变化趋势，而不是追求逐根精确建模。

结构单元生成后一定要可视化检查一遍。FLAC3D里结构单元节点和zone网格点位置接近时会自动连接，但偶尔也会出现没接上的情况，尤其是复杂曲面导入模型。plot sel 看一眼，发现问题及时补连接或重新生成，否则支护相当于悬空，结果完全失真。

3.3 下台阶跟进与“封闭成环”的次序控制

台阶法的核心就是下台阶及时跟上。上台阶支护完成，不能一直把断面拉成大跨度空腔。我的做法是上下台阶错开1到2个循环，比如上台阶开挖到第3循环时，开始挖下台阶第1循环。命令和上台阶一致，只是换一下group名称：

; 下台阶第1循环开挖：上台阶已经推进2个循环后再跟进 zone cmodel assign null range group 'lower_1' model step 300 sel shell create surface 'lower_face_1' thickness 0.25 sel property young 23e9 poisson 0.2 density 2400

初支和下台阶做完后，仰拱尽量尽早闭合。所谓“封闭成环”，在数值上就是让初期支护形成一个连续环路。如果上台阶不断往前挖、下台阶长期不跟进，shell单元在下部没有约束，拱脚位置位移会异常偏大，这在实际施工中也是大忌。因此，命令流里的开挖顺序必须和现场施工组织顺序对齐。

3.4 用FISH封装循环步进：批量开挖命令流的思路

循环多了，一条条手写命令不仅烦，还容易漏。我习惯把“开挖—释放—支护”封装成FISH函数，循环调用。

; FISH示意：把单循环台阶开挖封装成函数 fish define do_bench(cycle) zone cmodel assign null range group string.build('upper_%1', cycle) model step 300 sel shell create surface string.build('tunnel_face_%1', cycle) thickness 0.25 ... end do_bench(1)

这样改循环进尺、改支护参数，只需要改函数内部的少数行。前提是建模阶段分组足够规范，group名称带循环号且格式统一。如果分组是随手命名的，FISH函数写的再漂亮也无从下手。

4. 结果检查配命令：位移、塑性区、支护受力这么看

4.1 先布监测点再开挖：gridpoint history的正确用法

开挖之前，一定要先把monitor布好。FLAC3D 6.0里可以用zone gridpoint history把拱顶正上方、台阶拱脚、地表位置下沉量记录下来。每个监测点起一个通俗的名字，后面通过plot history直接看曲线。我习惯把历史曲线按施工顺序叠加：某一次开挖突变多大、支护后收敛趋势如何，一目了然。如果不提前布置，开挖完再补，前面的数据就永远补不回来了。

4.2 塑性区扩展怎么看工序合不合理

塑性区分布是最直观的工序反馈。命令上可以按开挖阶段保存多个model save文件，然后用plot zone state 查看塑性区。注意区分shear-n和tension-n，前者代表剪切破坏，后者代表拉破坏，隧道围岩里拉破坏往往更值得警惕。

经验上，如果上台阶开挖后塑性区集中在台阶拱脚位置，说明台阶高度或者错距设计偏大；如果下台阶开挖后塑性区沿洞周扩展过深，要么初支滞后，要么一次开挖段太长。这些从云图上一眼能看出来，也是我判断模型结果是否合理的第一个抓手。

4.3 从支护内力判断初支是否“吃力”

初期支护不能只看位移，还要看结构单元受力。shell单元看轴力、弯矩；cable单元看轴力。命令上用sel list shell、sel list cable一类命令可以输出这些数据。某一段shell轴力远超混凝土允许值，说明这一段初支太吃力，需要考虑缩短循环进尺、增加锚杆，或者调整台阶错距。用数值模拟的真正价值就在这里，不只是看位移够不够，而是辅助判断支护方案是否合理。

5. 实测踩坑记录：五个值得你注意的细节

5.1 zone delete 和 zone cmodel assign null：选不好就翻车

FLAC3D模拟开挖有两种常见做法：zone delete直接删除网格，或者把区域设为null本构。很多人看到网上的老命令喜欢用delete，但在台阶法模拟中我强烈建议优先用null。delete会删除网格本体，一旦删掉，后面要在开挖临空面生成shell单元，surface信息已经不存在了，支护根本挂不上。我在项目里就吃过一次亏，上台阶区域删得干干净净，初支怎么都生成不了，最后只能重新导入模型再来一遍。null本构的好处是网格还在，区域可以随时改回其他本构，适合工序模拟。

5.2 地应力位移清零的时机别手滑

初始地应力弹性平衡后，必须要做位移清零。如果没有清零，后续监测到的“开挖沉降”会叠加上自重平衡阶段的沉降，数据整体偏大，很容易误导决策。清零操作需要在开挖前完成，而且清完之后不要再动重力或边界条件。如果中途改了约束，就需要重新平衡再清零，否则历史曲线里面有一段“假位移”，对账时很头疼。

5.3 计算时间和收敛控制

FLAC3D是显式求解，台阶法这种分步开挖计算量不小。我的建议是正式跑批量工况前，先做一个循环的试算，看看每一步大概多少步数能收敛、耗时多少，再决定全模型怎么算。收敛目标设在1e-5不平衡比基本够用，再往下压只是在烧时间。另外隧道开挖这类问题一般用小变形模式就行，如果开了大变形，计算时间成倍增加，而位移结果往往变化不大。

5.4 6.0和5.0、以及老版本命令的差异

还是那句，FLAC3D 6.0延续了5.0的新命令体系，和4.0完全不是一回事。网上一搜flac3d台阶法命令流，很多是旧版gen zone brick那一套，直接复制到6.0基本跑不了。新版要看zone create、zone cmodel assign这种形式。另外6.0命令支持Tab补全，我写长命令流时会频繁按Tab看关键字，比翻手册快得多。

5.5 遇到断层破碎带怎么处理

如果隧道段落里刚好有断层穿越，别把它和普通围岩混在一起建模。我的做法是单独分一个group，给相对较低的粘聚力和摩擦角，比如完整围岩强度较高，断层带就按弱化参数处理。这样开挖到断层带时，塑性区和支护受力会明显异常，也更容易判断断层位置对施工的影响。真要模拟断层接触面的滑移，那要上interface单元，属于另一个复杂场景，不在本文讨论范围内。

这段命令流实践整体走下来，我最深的体会是：命令流不是写代码，而是把你脑子里的施工工序一步步翻译成FLAC3D的计算指令。台阶法模拟比全断面麻烦的地方，不是命令本身，而是你得先想清楚先挖哪、后挖哪、支护之前留多少应力释放步。如果每次改工况都需要从头梳理工序，那大概率是分组设计没做到位。把分组和循环封装做好，后面再换参数、换工序，也就是十分钟的事。