COMSOL多场耦合地应力平衡开挖与衬砌支护案例：带衬砌与钢衬支护的实践研究

COMSOL 地应力平衡后开挖及衬砌支护案例（带衬砌、钢衬）

隧道开挖模拟最头疼的就是初始地应力场的平衡问题。前些天用COMSOL折腾了个带衬砌支护的案例，今天把关键步骤拆开说说。咱们直接从地应力平衡开始，到开挖后钢衬安装一气呵成。

地应力平衡阶段最容易被忽视的是重力补偿。这里用固体力学接口设置初始应力场时，记得勾选"包含体载荷"选项。下面这段代码是平衡阶段的核心参数设置：

model.param.set('rho_rock', '2700[kg/m^3]') # 岩体密度 model.component("comp1").material("mat1").propertyGroup("def").set("youngs", "10E9[Pa]") model.component("comp1").material("mat1").propertyGroup("def").set("poissons", "0.25")

这里Young's modulus别设太高，否则后续开挖位移会小到看不见。重点来了，在初始研究步骤里要禁用自重载荷，用预定义位移场代替。相当于先给模型做个"全身按摩"让应力自然分布。

COMSOL 地应力平衡后开挖及衬砌支护案例（带衬砌、钢衬）

开挖模拟其实是个"欺骗"过程。通过材料切换实现，在开挖区域新建选择集：

model.component("comp1").geom("geom1").create("sel1", "Box"); model.component("comp1").geom("geom1").feature("sel1").set("entitydim", "3"); model.component("comp1").geom("geom1").feature("sel1").set("size", new String[]{"5", "3", "2"}); // 开挖尺寸

重点是把原岩体材料换成真空材料，同时保留位移场。这时候之前的平衡状态就被打破，围岩开始变形。注意要开启几何非线性选项，不然大变形计算会翻车。

衬砌支护的骚操作来了。新建壳接口，材料参数要跟围岩匹配：

model.component("comp1").material("concrete").propertyGroup("def").set("density", "2400[kg/m^3]"); model.component("comp1").material("concrete").propertyGroup("def").set("youngs", "30E9[Pa]");

这里有个坑——壳的局部坐标系方向必须正确。建议用边界选择工具指定法线方向，不然衬砌可能"长"在岩体外面。钢衬更麻烦，得用接触对处理界面滑移：

model.component("comp1").physics('cont').feature('con1').set('type', 'frictionless'); model.component("comp1").physics('cont').feature('con1').set('gap', '0.01[m]'); // 预留安装间隙

最后的重头戏是分步求解：先平衡，冻结初始条件；再开挖，记录位移场；最后激活衬砌约束。这时候应该看到典型的"猫眼"形变分布，最大位移出现在拱腰位置。要是出现应力奇点，八成是网格在开挖角点处太粗糙，局部加密就能解决。

整个模型跑下来最直观的发现是：二次衬砌施加时机对最终位移影响巨大。早支护虽然位移小，但可能造成过大围岩压力；晚支护又容易引发塌方。建议用参数化扫描找最佳支护时机，这个留着下次再细说。