最近在折腾UDEC模拟岩石单轴压缩试验，发现要画出完整的应力应变曲线还真得注意几个关键操作。直接上干货，分享个实测可用的案例代码，顺便聊聊实现细节
UDEC 7.0单轴压缩案例代码，含全应力应变曲线
先看核心代码骨架：
model new model large-strain on block create polygon 0,0 5,0 5,10 0,10 block zone generate quad 10 block zone cmodel elastic block zone property density 2500 young 5e9 poisson 0.25 block contact jmodel assign mohr block contact property stiffness-normal 2e9 stiffness-shear 1e9 friction 35 block fix x range y -0.1 0.1 block fix x range y 9.9 10.1 block history velocity-y range y 9.9 10.1 program call 'strain.fis'这个模型做了个10m高的长方体试件。重点在接触面设置——用Mohr-Coulomb准则模拟潜在破裂面。材料参数这里用的花岗岩典型值，实际操作记得根据岩样调整。
监测点定义是出曲线的关键，用FISH函数抓数据：
def strain_monitor global ydisp = block.head.y(2) - block.head.y(1) global strain = ydisp / 10.0 # 试件高度10m global stress = -block.gp.force.unbal.y(2) / 5.0 # 截面面积5m² end这里有个坑：UDEC的节点力需要取不平衡力才能反映真实应力状态。很多新手直接用接触力会导致曲线异常。
UDEC 7.0单轴压缩案例代码，含全应力应变曲线
加载循环控制是另一个要点：
block apply velocity-y -1e-5 range y 9.9 10.1 model cycle 5000 block apply velocity-y 0 model solve ratio 1e-3分阶段加载能更好捕捉峰后行为。前5000步匀速加载，之后改用自动收敛控制，这样既能保证加载速度，又能在试件破坏后继续计算。
数据处理脚本示例（Python）：
import matplotlib.pyplot as plt strain = [line.split()[1] for line in open('strain.dat')] stress = [line.split()[3] for line in open('stress.dat')] plt.plot(strain, stress, 'r-', lw=2) plt.xlabel('Axial Strain (%)') plt.ylabel('Axial Stress (MPa)') plt.show()典型输出曲线会经历弹性段-屈服-峰后跌落三个阶段。如果曲线出现锯齿状抖动，可能是接触刚度设置不合理，建议调整stiffness-normal参数到1e10量级试试。
最后给个调参小技巧：当模型总是不收敛时，先试试把材料强度临时调低，运行几步后再恢复真实参数，相当于给模型"热启动"，这招对处理复杂破坏形态特管用。