混凝土塑性损伤模型表格解读by自习菌

ABAQUS塑性损伤模型计算表格解读by自习菌（wx公众号）

受压本构：

fc,r：砼单轴抗压强度标准值，可根据需要取多种值，此处取fck轴心抗压强度标准值

fck：C30,20.1MPa;C35,23.4;C40,26.8;C45,29.6;C50,32.4;C55,35.5;C60,38.5

εcr：与单轴抗压强度fc,r相应的峰值压应变，规范附录公式

αc：单轴受压应力-应变曲线下降段的形状参数，规范附录公式

εcu：应力-应变曲线下降段0.5 fc,r对应的压应变

εcu/εcr：规范附录公式

可适当修正抗压强度代表值fcr，峰值压应变εcr，以及曲线形状参数αc，砼规C.2.4附录。

Ec：弹性模量，只是辅助计算的一个临时取值。C30,3e4MPa;C40,3.25e4;C50,3.45e4

ρc：规范公式

n：规范公式

x：穷举数列，按规范公式与ε、εcr相关

dc：单轴受压损伤演化参数，以x=1为界限，规范为分段公式

ε：由x计算出，规范公式

σ：规范公式

σ修正：在应力-应变曲线上选定弹性阶段与弹塑性阶段的分界点，按Susoo88取0.4 fc,r，

或取1/3~1/2 fc,r，可见这也是一个可调整的值。通过这个选定的点的应力应变，计算弹性阶段的斜率，即E0弹性模量，这个弹性模量就是所采用本构的弹性模量，用E0和ε再重新反算弹性阶段的σ，即得“σ修正”。

对于C30砼，fc,r 取fck=20.1MPa，0.4*20.1=8.04MPa，在表格中插入一行，定义一个ε值，使σ无限逼近8.04（此时尚需重新定义表格这一行x列公式，使之由ε列导出）。根据这个应力应变值，求出E0，再由E0修正弹性阶段的应力值（即插入行之上的部分）。

【Susoo88：受压曲线与受拉曲线弹性临界点不一样，会产生两个弹性模量，需要在输入时选较大值，不然在后面导算等效塑性应变时会出错…】

σtrue，εtrue：之前得到的应力应变是“名义”应力应变，需要在此转换成真实应力应变。

表格中的红色部分为abaqus中的输入数据。

屈服应力，采用真实应力

非弹性应变，公式见表格及Susoo88，第一行要手动取0，否则软件报错

bc、塑性应变、损伤因子，公式见表格及Susoo88文章

最后一列的塑性应变检查，塑性应变值应随着非弹性应变的增加而增大，且不为负值，否则软件报错。表格中Subopions是啥？

受拉本构：

按规范公式，受拉应力应变曲线在从0到ftr并非线性变化，Susoo88文章有误。

但表格中E0的计算却也认为是线性变化，与规范不一致？

是否也应仿照受压本构选取一个点指定为弹性阶段的临界点？

ρt通过Ec定义，为什么不直接用E0定义？

此处E0是否应按Susoo88所说跟受压时的E0取大值？

按照受拉时E0的计算方法，临界点即取为峰值点，故由此计算的弹性模量几乎恒小于受压时的E0。据此，受拉本构在abaqus中只用输入峰值之后的下降段。

表格中的红色部分为abaqus中的输入数据。

开裂应变，定义方式同“非弹性应变”。

bt、塑性应变、损伤因子，公式见表格及Susoo88文章

最后一列的塑性应变检查，塑性应变值应随着非弹性应变的增加而增大，且不为负值，否则软件报错。

by自习菌