导读：为简便钢筋混凝土构件或者结构的本构模型设置，本期给大家推荐一款Abaqus混凝土CDP模型插件，供大家应用参考。这个插件无需繁琐的Excel操作，仅需选择混凝土等级即可在Abaqus前处理界面一键生成混凝土CDP本构曲线，且可任意调整本构曲线长度，并可对极限强度进行修正，且适用于不同的力、位移单位，可用于各类混凝土构件及结构的精细化分析。

一、写在文前

对于钢筋混凝土构件或者结构而言，正确合理的本构模型是对构件或结构进行非线性分析的关键。ABAQUS提供三种混凝土本构关系模型，分别为脆性开裂模型、弥散开裂模型及损伤塑性模型，其中，混凝土损伤塑性(Concrete Damaged Plasticity，CDP)模型是通过将各向同性下损伤弹性与拉伸和压缩塑性相结合的方式来对混凝土的非弹性行为进行描述的，适用于Standard和Explicit两大求解模块，可用于模拟混凝土在任意荷载作用下的受力情况，同时考虑了由于拉、压塑性应变导致的弹性刚度的退化以及循环荷载作用下刚度的恢复，具有较好的收敛性。有关CDP模型的介绍及应用可见推文：

【JY】浅谈混凝土损伤模型及Abaqus中CDP的应用

二、程序可解决的问题

采用ABAQUS模拟梁柱节点时，ABAQUS中CDP模型损伤系数计算到0.9和损伤系数计算到0.99所得的滞回曲线相差甚大，笔者建立了现浇梁柱节点模型对此进行了验证。

CDP模型本构曲线末尾段的选取，对滞回曲线下降段的影响较大。相信大部分人是根据《混凝土结构设计规范》编制Excel表格来计算混凝土应力应变关系和损伤系数-塑性应变的关系，采用计算时往往会遇到一个问题，表格中所取的点太少，导致曲线不够精确，对于塑性应变较大时的曲线没有进行计算，混凝土过早地退出工作，与实际试验所得曲线不符。如果自己人为地进行计算的话，计算工作量较大，十分麻烦。

因此为了方便大家的应用，笔者开发了ABAQUS混凝土CDP模型插件，省去了繁琐的Excel计算，只需选择混凝土等级并选择性地修改参数，便可直接在Abaqus前处理界面直接生成混凝土本构，与大家一起分享。

三、操作界面

操作界面各参数意义：

四、本构生成

五、混凝土CDP模型插件正确性验证

1、有限元模型建立

为了验证所编子程序的合理性与正确性，选用清华大学陆新征教授钢筋混凝土框架结构拟静力倒塌试验中的边柱构件。

分别建立了三个有限元模型，有限元模型如图2。编号分别为Column-1、Column-2、Column-3。其中Colunm-1和Column-2区别在于：Column-1混凝土本构模型采用编写的JYCDP子程序，Column-2采用的是GB 50010-2010推荐的混凝土本构关系，其余参数均一致。Column-1和Column-3区别在于，Column-1采用子程序单位为N,mm，Column-3采用的子程序单位为kN,m，其余参数均一致。

滞回曲线的捏缩现象通常是由于钢筋混凝土界面粘结滑移和混凝土保护层剥落等引起的。往复荷载作用下，钢筋材料需要考虑包辛格效应，实际中包辛格效应的影响因素非常复杂。为了便于分析计算，Clough较早提出了带有再加载刚度退化的双折线滞回模型。

清华大学曲哲将Clough提出的双折线滞回模型进行了修正，提出了改进的 Clough钢筋滞回本构模型，该模型在反向再加载时，指向按卸载刚度加载至历史最大点对应的应力的 0.2 倍，即 0.2fmax，再指向历史最大点，能较好地考虑钢筋加载-卸载-反向加载过程产生的包辛格效应。

为了反映出滞回曲线的捏缩效果，采用了考虑承载力退化的Clough模型，直接在ABAQUS中通过用户子程序调用。

子程序钢筋属性：

提交作业界面：

2. 本构正确性验证

Column-1所得结果与试验结果对比：

Column-1所得结果和Column-2所得结果对比：

3. 力、位移单位正确性验证：

为了符合大家的使用习惯，所编CDP插件的力、位移单位可选择为kN和m，将两种单位下计算结果相对比：

六、写在文后

1、公众号后台回复“JYCDP”，可获得相关下载链接方法。

2、点击“阅读原文”也可获得相应链接。

【插件使用方法】

• 电脑路径下输入 %homepath%\abaqus_plugins并回车；

• 将下载的文件夹解压至当前目录下(注意文件夹名称为：JYCDP)；

• 打开Abaqus，菜单栏中点击plug-ins，里面找出JYCDP并点击；

• 选择相关参数并点击确定。

（完）

欢迎关注，下期更加精彩~

往期精彩

#性能分析

【JY】基于性能的抗震设计浅析（一）

【JY】基于性能的抗震设计浅析（二）

【JY】浅析消能附加阻尼比

【JY】近断层结构设计策略分析与讨论

【JY】浅析各动力求解算法及其算法数值阻尼(人工阻尼)

理念

【JY|体系】结构概念设计之(结构体系概念)

【JY|理念】结构概念设计之(设计理念进展)

【JY】有限单元分析的常见问题及单元选择

【JY】结构动力学之显隐式

【JY】浅谈结构设计

【JY】浅谈混凝土损伤模型及Abaqus中CDP的应用

【JY】结构概念设计之（隔震概念设计）

【JY】结构概念设计之（消能减震黏滞阻尼器）

#概念机理

【JY】基于Ramberg-Osgood本构模型的双线性计算分析

【JY】结构动力学初步-单质点结构的瞬态动力学分析

【JY】从一根悬臂梁说起

【JY】反应谱的详解与介绍

【JY】结构瑞利阻尼与经济订货模型

【JY】主成分分析与振型分解

【JY】浅谈结构多点激励之概念机理（上）

【JY】浅谈结构多点激励之分析方法（下）

【JY】板壳单元的分析详解

【JY】橡胶支座的简述和其力学性能计算

【JY】振型求解之子空间迭代

【JY】橡胶支座精细化模拟与有限元分析注意要点

【JY】推开土木工程振型求解之兰索斯法(Lanczos法)的大门

【JY】基于OpenSees和Sap2000静力动力计算案例分析

【JY】建筑结构施加地震波的方法与理论机理

【JY】消能减震黏滞阻尼器的力学原理与应用

#软件讨论

【JY】复合材料分析利器—内聚力单元

【JY】SDOF计算教学软件开发应用分享

【JY】Abaqus案例—天然橡胶隔震支座竖(轴)向力学性能

【JY】Abaqus6.14-4如何关联fortran？

【JY】如何利用python来编写GUI？

【JY】如何解决MATLAB GUI编程软件移植运行问题？

【JY】浅谈结构分析与设计软件

【JY|STR】求解器之三维结构振型分析

【JY】SignalData软件开发应用分享

【JY】基于Matlab的双线性滞回代码编写教程

#其他

【JY】位移角还是有害位移角？

【JY】如何利用python来编写GUI？

【JY】今日科普之BIM

~关注未来更精彩~