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戴个增高垫这招往往能让数值稳定性起死回生。这个模型的妙处在于能同时捕捉瞬时变形和长期蠕变。下次做边坡长期稳定性分析时不妨把博格斯摩尔模型拉出来遛遛说不定能挖出些常规分析漏掉的时间效应彩蛋。