相场断裂模型分步算法在ABAQUS中的实现  

Implementation of a staggered algorithm for a phase field model in ABAQUS

作　　者：刘国威[1] 李庆斌[1] 左正[1,2] 

机构地区：[1]清华大学水沙科学与水利水电工程国家重点实验室,北京100084 [2]中国电能成套设备有限公司,北京100080

出　　处：《岩石力学与工程学报》2016年第5期1019-1030,共12页Chinese Journal of Rock Mechanics and Engineering

基　　金：国家自然科学基金资助项目(51339003,51579134);; 国家重点基础研究发展计划(973)项目(2013CB035902)~~

摘　　要：首次基于ABAQUS平台实现相场断裂模型的分步算法,利用UMAT/VUMAT子程序接口对非线性平衡方程进行隐式/显式求解,提出位移场和相场的分步求解算法。准静态加载下的I和II型裂缝扩展问题、冲击荷载作用下的I型动态裂缝分叉和II型动态裂缝扩展问题,计算所得结果与已有文献和试验成果基本一致,说明了方法的可靠性;同时,对岩石力学领域中的双平行翼型裂缝相交和动态裂缝三维曲面扩展等问题,也利用该方法进行模拟。计算实例表明,结构储存的较高弹性应变能是脆性材料动态裂缝分叉的主要原因,而相场断裂模型是一种能够计算包括裂缝起裂、相交、分叉和三维曲面扩展的有效断裂计算方法,借助通用有限元平台可以方便地实现其分步算法。A staggered updated method for a phase field model was implemented in the commercial finite element software ABAQUS through UMAT and VUMAT subroutines. In order to verify the reliability of the algorithm,crack propagation in modes I and II under quasi-static and dynamic loads was calculated. All the results are generally consistent with the testing results in the existed references. In addition,simulations for wing cracks and curved surface cracks were also carried out. The results show that the main reason of dynamic crack branching is the high elastic strain energy stored in solids. The algorithm of phase field model is effective to simulate crack initiation,intersection,bifurcation and propagation in three-dimensional space,and can be executed conveniently in commercial FEM software.

关 键 词：岩石力学 ABAQUS UMAT/VUMAT 相场模型 分步算法 裂缝分叉 

分 类 号：TU45[建筑科学—岩土工程]