混凝土这玩意儿在工程界真是让人又爱又恨。今天咱们聊聊用CDP模型做双轴压-压的细观模拟，特别是骨料、砂浆、界面过渡区这三兄弟全用上损伤塑性模型的骚操作
混凝土双轴细观数值模型，有压-压，三相材料均采用cdp模型，
先说个有意思的现象——混凝土在双轴受压时强度比单轴还能提高10%-20%，但你要是直接拿宏观模型硬怼，结果可能比甲方改需求还离谱。这时候就得靠细观建模把骨料、砂浆、界面过渡区拆开搞事情。
ABAQUS里头的CDP模型参数设置就像调火锅蘸料，比例不对直接翻车。拿界面过渡区举个栗子：
mdb.models['Model-1'].materials['ITZ'].ConcreteDamagedPlasticity( dilationAngle=30.0, eccentricity=0.1, fb0=1.16, viscosity=0.0005) mdb.models['Model-1'].materials['ITZ'].concreteCompressionHardening( table=((1.35e-5, 0.0), (2.5e-3, 26.8), (5.0e-3, 30.0)))这个黏性系数viscosity=0.0005是个暗坑，新手容易直接照搬砂浆参数。实测发现界面区黏度得比砂浆低20%左右，否则计算收敛速度比蜗牛爬还慢。压缩硬化曲线里那个5.0e-3的极限应变值，实验室数据要是超过这个数，模型分分钟表演数值发散。
骨料的刚度可不是越大越好，见过有人把弹性模量怼到80GPa，结果应力云图直接出现俄罗斯方块效应。建议试试分段刚度：
*Elastic, dependencies=1 50.0e3, 0.2, 0.0 *Plastic 0.0, 0.0 0.005, 70.0e3第一行初始模量给50GPa，等应变超过0.005后跳涨到70GPa。这招能模拟骨料在高压下的晶体结构调整，比死磕固定刚度靠谱多了。
混凝土双轴细观数值模型，有压-压，三相材料均采用cdp模型，
压-压工况最骚的是边界条件设定。见过有人直接锁死两侧位移，结果破坏模式变成抽象派艺术。推荐用动态荷载比例控制：
loadProportions = [ (0.5, 0.5), (0.7, 0.3), (1.0, 1.0) ] for ratio in loadProportions: applyBiaxialLoad(model, ratio[0], ratio[1]) solveIncrement(0.01)这个分步加载策略比一次性加载稳定三个段位。特别是最后那个(1.0,1.0)的全荷载步，必须把时间增量压到0.01以下，不然混凝土还没破坏，你的CPU先烧了。
最后说个血泪经验：CDP模型里的损伤变量输出千万别只看最后结果。建议每个荷载步都输出拉伸损伤和压缩损伤，你会惊奇地发现界面区的损伤演化路线比网红打卡路线还曲折。有次模拟结果和试验对不上，后来发现是压缩损伤权重设成了0.5，实际应该用0.45，就这0.05的差距让整个破坏模式从剪切破坏变成劈裂破坏。
搞细观建模就像拼乐高，参数调不好分分钟塌房。不过当你看到应力云图里骨料像跳棋一样把荷载传递出去，界面区损伤带像血管网一样蔓延时，那种快感比通关黑魂还爽。