flac3d耦合pfc3d隧道开挖模拟。 位移连续性良好地表沉降规律合理。隧道施工总让人头大尤其是遇到软弱围岩的时候。上次帮设计院做地铁暗挖段模拟传统连续体方法死活算不出颗粒破碎后的应力重分布。灵机一动把FLAC3D和PFC3D这对冤家凑成了CP结果竟然比相亲角还和谐。搞耦合模拟就像给挖掘机装机械臂——得找准连接点。FLAC负责大范围连续介质PFC管颗粒间的微观运动。关键在数据交换的频率这里偷懒直接用了zone-node耦合模式; FLAC3D脚本片段 zone cmodel assign elastic zone node initialize velocity-x 0.0 zone couple pfc id1 ; 指定耦合区域编号对应的PFC端得设置镜像边界; PFC脚本片段 wall generate id 1 type segment ... ball attribute displacement multiply 0.0 ; 冻结边界颗粒这里有个坑——FLAC的网格节点和PFC的颗粒必须严格对齐。我直接导出FLAC网格顶点坐标生成匹配的PFC墙体就像用3D打印机复刻模具似的。开挖步设置最考验节奏感。每推进0.5米就触发一次数据同步program call excavate.dat zone delete range group tunnel ; FLAC开挖命令 ball activate range geometry face ; 激活PFC休眠颗粒这时候FLAC的应力释放会自动传递给PFC颗粒系统。有次手滑把同步间隔设成2米结果地表沉降曲线直接变成心电图——典型的步长过大引发数值震荡。flac3d耦合pfc3d隧道开挖模拟。 位移连续性良好地表沉降规律合理。位移连续性的验证得玩大家来找茬。截取FLAC网格节点和邻近PFC颗粒的位移量用Python画个散点图import matplotlib.pyplot as plt flac_disp [d[:,0] for d in flac_data] pfc_disp [p[:,0] for p in pfc_data] plt.scatter(flac_disp, pfc_disp, alpha0.5) plt.plot([-5,5],[-5,5],r--) # 理想拟合线发现90%的数据点都落在±2mm误差带内比甲方要求的5mm标准还严苛。特别是拱顶下沉的吻合度连老师傅看了都嘀咕这数值模型怕不是装了陀螺仪地表沉降槽的形态更有意思。传统方法算出的沉降曲线像平滑的抛物线耦合模型却在隧道中线两侧10米处多了两个小鼓包。翻出现场监测报告一对比——好家伙上周刚出现的局部隆起被这模型蒙对了后来发现是PFC颗粒在应力释放时发生了侧向弹跳。这种跨尺度模拟最吃硬件。原本单用FLAC只要半小时的计算加上PFC后得挂着服务器跑通宵。有次半夜收到报警短信以为模型崩了结果只是咖啡机漏电跳闸——虚惊一场也算模拟之外的乐趣吧。