考虑对称性,案例采用1/4圆进行模拟,单层单元。建模命令流如附件。
采用拉格朗日流固耦合算法,建模命令流中装药不耦合系数1.5,空气单独建模。在实际计算过程中,未考虑空气受热膨胀对裂隙尖端产生的拉力,希望有相同研究兴趣的可以一块交流一下。
采用ALE算法(设置空单元)还是拉格朗日算法(流固耦合)取决于*SECTION_SOLID_ALE关键字中ELFORM的设置(11or12),设置11为流固耦合,需要控制*ALE_MULTI-MATERIAL_GROUP关键字,设置12为空单元,需要控制*INITIAL_VOID_PART和*CONSTRAINED_LAGRANGE_IN_SOLID关键字。
裂隙的产生扩展通过添加*MAT_ADD_EROSION控制材料失效标准。本案例设置的是mnpres失效,即拉应力失效,具体数值根据材料动载破坏强度和期望效果调整。