flac3d隧道开挖面失稳模拟。
凌晨三点的办公室键盘声格外清脆，我盯着屏幕上跳动的应力云图，红色区域像癌细胞扩散般在隧道掌子面蔓延。这已经是第三次调整参数了，混凝土衬砌的位移量还是超过警戒值。FLAC3D的黑色命令行窗口不断吐出计算日志，仿佛在嘲笑我的无知。
新建模型时总得先处理网格这个磨人的小妖精。用zone create brick命令拉出个长方体区域，尺寸设置得比实际隧道大3倍——吃过边界效应的大亏，上个月就是因为模型范围太小，模拟出来的塑性区直接怼到了边墙上。材料参数表摊在左手边，摩尔-库仑模型里的内摩擦角数值被我反复圈画，这个值每调整1度，掌子面稳定性就会发生翻天覆地的变化。
; 创建隧道主体网格 zone create brick size 12 8 12 zone cmodel assign mohr-coulomb zone property density 2500 young 5e9 poisson 0.3 ... friction 35 cohesion 1e5 dilation 5这段FISH代码里的dilation参数特别有意思，它控制着材料屈服后的体积变化。有次我把这个值设成0，结果掌子面像融化的冰淇淋一样塌下来，吓得甲方连夜打电话问是不是设计有问题。其实这正说明岩体的剪胀效应在支护设计中不容忽视。
flac3d隧道开挖面失稳模拟。
开挖步的设置是门艺术，用zone relax命令释放应力时，每次开挖进尺控制在0.5米比较稳妥。但上次为了赶进度，我直接设置2米一步，结果计算到第5步就出现负体积警告，模型直接崩了。现在学乖了，老老实实写循环：
; 分步开挖实现 def excavate loop n (1,10) command zone relax range cylinder end1 0 0 0 end2 0 0 10 radius 3.5 solve save ['excavate_' + string(n) + '.sav'] endcommand ; 实时监测拱顶位移 gp = gp.near(0,0,5) disp = gp.disp.z if disp > 0.1 then exit endif endloop end excavate监测点的位移判断是个巧妙的逃生阀，当z方向位移超过10cm就自动终止计算，避免无意义的迭代。有次忘记加这个判断，模型跑了整晚，早上发现掌子面位移都到两米了——现实里早塌方八百回了。
后处理阶段盯着位移矢量图，那些红色箭头像无数把匕首指向隧道中心。突然发现靠近拱腰的位置出现了反向位移，这可不是好兆头。调出土体塑性区分布图，果然在拱脚位置看到了连成片的剪切屈服区。这时候就得祭出土体加固大法，在关键区域插入锚杆单元：
; 锚杆加固实现 struct cable create by-line 0 0 5 0 3 5 segments 10 struct cable property young 2e11 yield-tension 3e5 ... cross-sectional-area 0.005 struct link attach xyz 0 3 5 range cylinder end1 0 0 5 end2 0 3 5 radius 0.2锚杆的预紧力设置得像重庆小面的辣度，少了没劲，过了会脆断。有回把yield-tension设成5e5，结果计算到第三步锚杆就开始大面积失效，跟放鞭炮似的在模型里噼里啪啦报错。
现在看收敛曲线就像看心电图，当不平衡力在1e3量级反复震荡时，就知道该调整求解器参数了。set mech ratio这个值调到0.01能提高计算精度，但耗时直接翻倍。有时候在深夜听着主机风扇的嘶吼，会恍惚觉得是隧道里的风钻声——或许这就是数值模拟的魅力，在代码与方程构建的虚拟世界里，提前预演着现实工程的血与火。
