在COMSOL中运用水平集法和蠕动流模块模拟裂隙注浆过程考虑浆液—岩体的耦合作用。 一般而言裂隙开度越大浆液所需注入压力越小。 本算例从结果来看可以验证此定律。 裂隙变形的本构取之于已发表的文献。 本算例中初始时刻裂隙内部只存在水1mPas然后对裂隙进行长达40秒的注浆浆液粘度为10mPas。 保持注浆速率不变注浆压力逐渐增大符合一般注浆规律。 诚然由于comsol的收敛较难在初始时刻注浆压力可能出现数值震荡。打开COMSOL的模型树先往材料库里塞两个流体——水相和浆液相。这里有个细节容易被新手忽略给材料赋值时得注意单位制是否统一。特别是粘度参数1mPas的水和10mPas的浆液差了整整一个数量级这种量级差直接影响压力场分布。% 材料参数设置伪代码 material(water) .set(dynamic_viscosity, 1e-3); % 1mPas换算为Pa·s material(grout) .set(dynamic_viscosity, 10e-3);水平集函数像染色剂一样标记两种流体界面这里建议开启重新初始化功能。最近帮学弟调模型时发现不勾选这个选项的话计算到第20秒左右界面会开始模糊就像老式电视机信号不好时的雪花噪点。蠕动流模块的处理需要点小技巧。当裂隙宽度变化超过初始值的5%时建议启动几何非线性选项。遇到过不收敛的情况吗试着把相对容差从默认的0.01调到0.005虽然计算时间会翻倍但压力曲线的毛刺明显减少。在COMSOL中运用水平集法和蠕动流模块模拟裂隙注浆过程考虑浆液—岩体的耦合作用。 一般而言裂隙开度越大浆液所需注入压力越小。 本算例从结果来看可以验证此定律。 裂隙变形的本构取之于已发表的文献。 本算例中初始时刻裂隙内部只存在水1mPas然后对裂隙进行长达40秒的注浆浆液粘度为10mPas。 保持注浆速率不变注浆压力逐渐增大符合一般注浆规律。 诚然由于comsol的收敛较难在初始时刻注浆压力可能出现数值震荡。观察压力监测点的数据时前5秒的震荡曲线特别有意思。这其实是数值计算中的启动效应就像踩油门时发动机的短暂抖动。可以通过设置初始步长0.1秒来缓解但别指望完全消除——真实注浆现场的传感器数据也会有类似的抖动。% 求解器配置示例 solver.create(st1, StudyStep); solver.feature(st1).set(tlist, range(0,0.5,40)); solver.feature(tols).set(rtol, 0.005);在结果分析阶段双击压力云图时注意颜色标尺范围。有个反直觉的现象最大压力点并不在注浆口而是出现在裂隙中段。这其实是粘性流体在狭窄通道里的瓶颈效应和挤牙膏时中间鼓起的现象异曲同工。验证开度-压力关系时推荐用参数化扫描同时跑三个不同开度的模型。数据导出后做对数坐标系下的线性拟合斜率绝对值越大说明耦合作用越明显。记得检查网格质量——粗糙的网格会让斜率失真就像用马赛克图片做边缘检测。最后提醒当看到注浆压力曲线呈现类指数增长时别慌这恰恰验证了流固耦合的作用机制。保存结果前务必做网格无关性验证特别是裂隙边缘的边界层网格至少要保证三层单元才能捕捉到流速梯度变化。