ABAQUS岩土工程案例集，边坡隧道基坑实战解析

边坡强度折减算到折减系数1.35的时候Standard崩了，msg里全是塑性区贯通导致的收敛失败。岩土这三类模型，边坡、隧道、基坑，看着不一样，地应力平衡和材料本构这两关过不了哪个都跑不动。

边坡用Mohr-Coulomb，内摩擦角25°，黏聚力15kPa，弹性模量30MPa，泊松比0.3，重度18kN/m³。坡高10m，坡比1:1.5，模型范围取坡脚到边界2倍坡高，坡顶到边界2倍坡高，底部全固定，两侧滚轴。强度折减法用FIELD, VARIABLE=1，把摩擦角和黏聚力都跟场变量挂钩，INITIAL CONDITIONS, TYPE=FIELD, VARIABLE=1，初始值1.0。STEP, NLGEOM=ON，STATIC, DIRECT，*CONTROLS, PARAMETERS=FIELD，收敛失败时的折减系数就是安全系数。判定标准用塑性区贯通或者坡顶位移突变，两个一起看。

隧道用Drucker-Prager或者修正剑桥模型，衬砌用弹性。开挖用MODEL CHANGE, REMOVE，衬砌激活用MODEL CHANGE, ADD=STRAIN FREE。地应力平衡先算一步，GEOSTATIC，侧压力系数K0取0.5。开挖释放荷载用MODEL CHANGE配合*BOUNDARY, OP=NEW，释放系数0.7左右。隧道直径6m，埋深20m，衬砌厚0.3m，C30混凝土，弹性模量30GPa。

基坑用Mohr-Coulomb或者HS模型，地连墙用壳单元或者实体单元，内支撑用梁单元。开挖分步，每步挖2m，MODEL CHANGE, REMOVE，支撑架设用MODEL CHANGE, ADD=STRAIN FREE。围护结构跟土体的接触面用Surface-to-Surface，摩擦系数0.3。基坑深15m，宽30m，长50m，地连墙厚0.8m，深25m。

有一次做基坑开挖模拟，分五步开挖，每步挖3m。前四步跑得好好的，第五步刚开始就崩，msg里全是“excessive distortion”和“zero pivot”。那天下午刚下过雨，办公室窗户没关严，地上湿了一片，隔壁工位同事在拆快递，胶带撕得响。第一反应是接触面问题，把地连墙跟土体的摩擦系数从0.3降到0.2，加了*CONTACT CONTROLS, STABILIZE，重新提交，还是崩。怀疑是网格不匹配，土体靠近地连墙的网格从1m加密到0.5m，跑了一晚上，还是不行。翻dat文件发现地连墙底部的接触面定义反了，主面设成了土体，从面设成了墙，硬接触下刚度大的做主面，反了之后穿透控制失效，墙底穿进土里了。改回来之后好了。后来又发现第五步开挖释放系数设成了0.7，前四步都是0.5，释放太快导致土体回弹过大，改成0.5之后位移跟实测数据对上了。

地应力平衡是岩土模拟的第一道坎，*INITIAL CONDITIONS, TYPE=STRESS, GEOSTATIC，给土体一个初始应力场，K0按经验取值，砂土0.4-0.5，黏土0.5-0.7。平衡完第一步的位移应该接近零，如果位移量级跟后续开挖位移差不多，说明地应力没平衡好。PLAXIS和FLAC3D在地应力平衡这块比ABAQUS省事，有专门的K0 procedure和initial stress generation，但ABAQUS胜在UMAT接口灵活，土体本构可以自己写，边界面模型、亚塑性模型都能塞进去。

后处理边坡提取坡顶位移和折减系数的关系曲线，用FIELD OUTPUT, VARIABLE=1, FV1输出场变量，沿坡顶路径提取U2。隧道提取拱顶沉降和水平收敛，OUTPUT, HISTORY，U2和U1。基坑提取墙顶水平位移和地表沉降，沿垂直基坑方向路径提取U1和U2。

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删