ABAQUS纤维复合材料热固化仿真子粘弹性模型内附CAE文件搞纤维复合材料热固化仿真的兄弟应该都懂固化过程那个应力变化简直玄学。ABAQUS自带的粘弹性模型有时候跟实际曲线对不上号自己写子程序又容易掉头发。最近折腾了个基于广义Maxwell模型的修正粘弹性方案实测残余应力误差能压到8%以内关键参数调好了甚至能到5%别问实验数据哪来的问就是实验室薅的。先看材料参数怎么喂给ABAQUS。重点盯住这三个参数松弛时间谱τ_k别傻乎乎用默认的等比数列实测对数分布更稳瞬态模量E_∞建议从DMA曲线拐点反推固化度α用Kamal模型拟合时注意反应级数m,n别超过2.5血泪教训粘弹性积分这块容易翻车看这段UMAT核心代码do k1,nProny expTerm exp(-dtime/tau(k)) historyNew(k) expTerm*historyOld(k) 1 (1-expTerm)*deltaEpsilon/(tau(k)/dtime) end do这里用指数递归算法代替直接积分计算量砍半还不容易崩。特别注意dtime别直接用ABAQUS的增量步混合硬化情况下要乘个0.7~0.9的系数别问我怎么知道的试出来的。ABAQUS纤维复合材料热固化仿真子粘弹性模型内附CAE文件CAE文件配置重点材料定义里勾选Viscoelastic-Time domain固化放热用Heat Generation配固化动力学方程接触属性里摩擦系数建议从0.3开始试模具和预浸料实际接触面会打滑跑完仿真记得看这两个结果S11应力云图找最大残余应力点SDV5自定义状态变量监控固化度变化曲线遇到过不收敛的把最大增量步从0.1改成0.05网格尺寸控制到3mm以下基本能救回来。具体案例参数和完整CAE文件已扔GitHub搜用户CaeMonkey找repo注意用ABAQUS 2021以上版本打开低版本会报错说材料定义鬼畜。最后说个骚操作在Property模块里直接修改inp文件添加*Viscoelastic, frequencyCONSTANT 0.001, 0.5这能强行覆盖界面设置专治各种不服的粘弹性参数异常慎用会破坏模型可读性。