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)这里Youngs 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]); // 预留安装间隙最后的重头戏是分步求解先平衡冻结初始条件再开挖记录位移场最后激活衬砌约束。这时候应该看到典型的猫眼形变分布最大位移出现在拱腰位置。要是出现应力奇点八成是网格在开挖角点处太粗糙局部加密就能解决。整个模型跑下来最直观的发现是二次衬砌施加时机对最终位移影响巨大。早支护虽然位移小但可能造成过大围岩压力晚支护又容易引发塌方。建议用参数化扫描找最佳支护时机这个留着下次再细说。