abaqus、MIDAS GTS NX基坑隧道开挖模拟

基坑开挖模拟这事儿吧，搞岩土工程的都懂有多头秃。两年前我接手个地铁隧道项目，甲方非要同时看Abaqus和MIDAS GTS NX的对比报告，结果连着三周都在和单元生死技术死磕。今天就唠唠这两个软件在开挖模拟里的骚操作。

Abaqus的倔强

这货的单元生死逻辑就像个强迫症患者，非得用场变量（Field Variable）来标记激活状态。下面这段Python脚本估计不少人都改吐过：

mdb.models['Model-1'].fieldOutputRequests['F-Output-1'].setValues(
    variables=('S','PE','PEEQ','PEMAG','UVARM'))
mdb.models['Model-1'].historyOutputRequests['H-Output-1'].setValues(
    variables=PRESELECT)

看起来平平无奇？重点在UVARM这个用户自定义场变量。当你要干掉某层土体时，得在材料定义里搞个类似这样的判断：

if UVARM > 0.5:
    material.stiffness = 0
else:
    material.stiffness = original_value

这招其实是在材料层面玩狸猫换太子，把刚度矩阵直接归零。但新手常踩的坑是没同步处理接触关系，结果开挖后的结构体还在和空气摩擦较劲...

MIDAS的温柔

转战MIDAS GTS NX时发现，这货直接把开挖做成了可视化操作。右键点开施工阶段管理器，新增激活/钝化组的时候，简直像在玩模拟城市。不过它的命令流也有自己的脾气：

!施工阶段定义
CONST, NAME=Stage1
DEACTIVATE, ELEMENT, GROUP=soil_layer1
ACTIVATE, ELEMENT, GROUP=lining1

看起来清爽？但遇到复杂地层组合时，分组命名能让人疯掉。有次把"soillayer1"误写成"soillayar1"，程序居然不报错直接跳过，白算了六个小时才发现是拼写错误...

abaqus、MIDAS GTS NX基坑隧道开挖模拟

网格的玄学

两种软件都逃不过网格质量的毒打。Abaqus做隧道周边网格加密时，这个参数调节能把人逼疯：

p.setElementType(elemTypes=(ElemType(elemCode=C3D8R, elemLibrary=STANDARD), ), regions=(
    p.sets['tunnel_zone'], ))

C3D8R减缩积分单元确实节省算力，但塑性区发展大了就容易出现沙漏模式。反观MIDAS的自动网格划分，虽然号称智能，但生成的全六面体网格经常在支护结构转角处出现畸变单元。

后处理的黑魔法

算完看着云图发懵是常态。Abaqus里要抓支护结构轴力，得这么折腾ODB：

session.viewports['Viewport:1'].odbDisplay.setPrimaryVariable(
    variableLabel='S', outputPosition=INTEGRATION_POINT, )

而MIDAS直接右键就能出内力图，但导出数据时发现数值震荡得比股票K线还刺激——后来才发现是没在施工阶段设置结果平滑参数。

结个尾吧：用Abaqus就像开手动挡跑车，每个参数都得亲自调教；MIDAS则是自动挡SUV，上手快但想玩花活就得拆中控。最近发现把Abaqus的计算内核接到MIDAS前处理上用，居然能避开两者的致命缺陷...算了，这骚操作下次再唠。