FLAC3D蠕变三轴压缩试验：博格斯摩尔本构，应变时间曲线
在岩土工程数值模拟里，蠕变试验就像给材料做"慢动作回放"。今天咱们拿FLAC3D折腾个博格斯摩尔（Burgers-Malvern）模型的蠕变三轴压缩试验，重点看看怎么让应变曲线在时间轴上"跳起华尔兹"。
先给模型搭个台子。用六面体网格搞个标准圆柱试样，高度2m直径1m。这里有个细节：网格密度直接影响蠕变计算的稳定性。我习惯在轴向多分几层，毕竟蠕变变形主要沿着加载方向走。
; 生成圆柱试样 gen zone cyl p0 0 0 0 p1 1 0 0 p2 0 0 2 size 8 16 16材料参数设置是重头戏。博格斯摩尔模型把粘弹性和塑性揉在一起，参数多得能凑两桌麻将。特别注意粘性系数η1和η2——这俩兄弟控制着蠕变的快慢节奏。
; 材料参数设置 model config creep model burger-malvern property bulk 1e8 shear 5e7 property coh 1e5 ten 1e5 fric 35 property creep-mu 0.5 creep-lambda 0.3 property eta1 1e15 eta2 1e18边界条件要模拟真实三轴仪。底部固定，侧向用伺服控制保持围压。加载策略讲究"先快后慢"：初期快速加载到预定应力，之后保持恒定让蠕变自由发挥。
; 加载步骤设置 fix z range z 0 0.1 ; 底部固定 apply szz -1e6 hist gp 0 0 2 ; 顶部加载 apply sxx -0.5e6 syy -0.5e6 ; 侧向围压 ; 蠕变计算阶段 set creep dt auto solve time 1e5 ; 10^5秒约等于27小时当应变曲线开始画图时，注意看它的"三级跳"：先是弹性段垂直拉升，接着进入衰减蠕变阶段（曲线变缓），最后是稳态蠕变的线性爬升。要是看到曲线突然翘尾巴，可能是参数设飘了——这时候该检查粘性系数是不是太小。
FLAC3D蠕变三轴压缩试验：博格斯摩尔本构，应变时间曲线
数据提取时别只会用FLAC自带的绘图工具。我习惯把关键点的应变数据倒出来用Python处理：
hist gp zdisp 0 0 2 ; 记录顶部位移 hist write 1 vs 2 file strain_time.txt ; 时间-应变数据最后来个实战技巧：蠕变计算特别吃时间步长。看到计算不收敛时，别急着摔键盘，试试把set creep dtmin调小两个数量级。实在不行就给粘性系数η2戴个"增高垫"，这招往往能让数值稳定性起死回生。
这个模型的妙处在于能同时捕捉瞬时变形和长期蠕变。下次做边坡长期稳定性分析时，不妨把博格斯摩尔模型拉出来遛遛，说不定能挖出些常规分析漏掉的时间效应彩蛋。
