abaqus复合材料、钢筋混凝土、编织结构仿真合集

Too many attempts made for this increment，复合材料渐进损伤分析里Standard又卡住了。一加Hashin准则和损伤演化，收敛就成玄学，三个方向的模型轮着崩，记一下各自的关键点。

复合材料层合板，Property → Material → Elastic → Engineering Constants，下面挂Hashin Damage和Damage Evolution，再往下Cohesive Behavior和Damage Stabilization。E1取120-180 GPa，E2取8-12 GPa，G12取4-7 GPa，ν12取0.2-0.3。纤维拉伸强度XT取2200 MPa，压缩强度XC取1500 MPa，横向拉伸强度YT取50 MPa。Hashin准则不能用于实体单元，只能用壳或连续壳。Damage Evolution用Energy类型，输入断裂能Gf，黏性系数μ取1e-6到1e-5 Pa·s这个量级。

单元选择上，常规壳S4R适合薄壁层合板，连续壳SC8R必须用扫掠网格划分，扫掠方向就是铺层堆叠方向，厚度方向只允许一层单元。薄板结构用实体单元长宽比会很恶劣，实体单元加复合材料截面的组合基本不要碰。铺层定义的时候铺叠方向从顶面还是底面开始必须明确，0°/90°铺层顺序搞反了拉弯耦合刚度矩阵全错。曲面结构用离散方向，不能用笛卡尔坐标系，方向矢量要随曲面连续变化。Cohesive单元模拟层间分层，双线性本构加BK失效准则，损伤起始用二次应力准则，损伤演化用BK准则η取1.5左右，GIc、GIIc从材料手册查。

钢筋混凝土这边，混凝土用C3D8R，钢筋用T3D2，钢筋通过Embedded Region嵌入混凝土，把钢筋设为embedded region，混凝土设为主区域host region。CDP模型的膨胀角取36°，黏性系数取0.001，fb0/fc0按规范取1.16，Kc取0.667。拉伸硬化那个曲线如果不设，算出来的挠度和裂缝宽度会严重失真。嵌入约束不约束转动自由度，节点区精细分析要留意这一点。钢筋建模的时候箍筋弯钩部位别用直线段近似，会导致钢筋长度和实际不符。CDP模型中最常见的收敛困难源于拉伸裂纹局部化，黏性正则化参数是解决收敛的关键，取值通常在1e-4到1e-7之间。

编织结构跟上面两类都不一样。ABAQUS本身没有针对编织复合材料的原生支持，需要自己写子程序。介观尺度RVE用TexGen生成，施加周期性边界条件来保证变形协调，通过ABAQUS的equation约束实现主从节点相对位移的绑定。ABD矩阵通过六次单独的有限元模拟计算，每次施加一个方向的面内或面外变形。编织复合材料的纤维束在1方向和2方向性能近似相同，本构模型可以用UMAT实现弹塑性耦合损伤，纤维用弹性损伤模型，基体用耦合弹塑性损伤模型，抛物线屈服准则表征拉压差异。跟Ls-Dyna的MAT_162比，ABAQUS做编织冲击需要自己搭Cohesive加本构，配置复杂但可控性高。

有一次做一个碳纤维层合板开孔压缩，铺层顺序[0/45/-45/90]s，孔边应力集中明显。提交之后跑了六个小时，msg里全是“local-1 direction has zero length”，检查发现曲面过渡区的材料方向定义失效了。当时刚过完年回来上班，办公室暖气烧得太足，整个人昏昏沉沉的。以为是网格问题，把孔边网格从2mm加密到0.5mm，算得更慢了，报错一样。后来把曲面分区，每块单独定义离散方向，问题解决，时间已经耗掉了两天。同一个模型里如果还要挂UMAT和DFLUX，所有子程序必须写在同一个.for文件里，CMNAME分支区分材料类型，SDV编号提前规划好别重叠。

收敛调参上，复合材料Damage Stabilization的黏性系数取1e-6到1e-5之间，大于1e-4损伤演化会被过度平滑，峰值载荷偏高曲线偏缓。*CONTROLS, PARAMETERS=TIME INCREMENTATION里把平衡迭代次数从默认5次提到15次，给求解器更多机会。Standard实在收不住的时候用*DYNAMIC, EXPLICIT做准静态，加载时间拉到一阶模态周期的10倍以上，ALLKE/ALLIE控制在5%以内。后处理输出单层结果，Field Output Request里展开Stresses和Strains，输出位置选All Section Points in All Plies，损伤分析把SDV勾上，不然看不到损伤变量演化。

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删