对于纤维增强复合材料的模拟，在Abaqus中，集成了二维Hashin失效准则与多种损伤演化准则，但缺少三维的复合材料本构模型。参考一篇已发表的SCI文章，基于Abaqus使用Fortran语言建立三维平纹织物复合材料弹塑性、渐进损伤本构模型-Vumat子程序。平纹织物复合材料在1方向和2方向丝束性能近似相同。
首先介绍该子程序的使用方法与效果
1. 在Abaqus中建立三维复合材料模型，这里建立一个简单的方块。赋给材料方向，1,2方向分别表示丝束的方向，3方向表示垂直于1，2的方向，也就是面外方向。
2. 建立材料属性
3. 建立显示Explicit计算时间步，时间0.005，在场输出中勾选输出 SDV和 STATUS.
4. 划分网格，赋给Explicit 3D stress单元类型，边界条件根据需要设定即可，此处不再赘述。此处划分为一个单元，使用12方向往复加载卸载。建立Job，提交模型前在Job中选择该子程序，设置双精度计算。
5. 查看结果，等效塑性应变在卸载时没有变化，再次加载时剪切应力按照原来的路径返回，剪切损伤在卸载时也保持不变。
6. 将该子程序应用在低速冲击模型中，可以顺利运行。
接下来简要介绍该子程序的相关理论，子程序、参考的论文名称以及输入材料参数的对应含义打包在附件中。
弹性阶段总应力与总弹性应变之间的关系为
式中，是柯西应力，S0是柔度矩阵，
;是弹性应变。
考虑损伤，有效应力与实际的应力之间的关系为：
损伤演化准则为：
上述R大于1时，损伤起始，计算w：
然后计算损伤系数D：
塑性判断
等效应力：
a66=0.146为实验获取的值。
流动应力：
σ0为初始屈服应力，n=3.58和A=7.18e-8为实验值，屈服时有
使用关联流动法则：
等效塑性应变增量根据下式计算：
☛附件内容：1. inp算例模型(低速冲击工况，1/4模型，层间使用cohesive element) 2. 子程序 3 .使用方法 4.参考论文名称
▼ 附件内容扫描下方二维码获取 ▼
或点击文末阅读原文了解详情
往期推荐
基于Abaqus中Cohesive单元的钢筋混凝土粘接性能仿真分析（含详细操作过程及inp文件）
Abaqus中Cohesive单元传热仿真分析案例讲解-COH2D4T（含模型文件）
基于Abaqus的旋压成形分析案例详解（含仿真源文件）
基于Abaqus Dload子程序在圆柱面上施加旋转载荷仿真分析案例讲解（附子程序及CAE文件）
