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

doi:10.13722/j.cnki.jrme.2015.1264

摘要
图/表
参考文献
相关文章 (15)

全文: PDF (24927 KB) HTML (1 KB)
输出: BibTeX | EndNote (RIS)

摘要首次基于ABAQUS平台实现相场断裂模型的分步算法，利用UMAT/VUMAT子程序接口对非线性平衡方程进行隐式/显式求解，提出位移场和相场的分步求解算法。准静态加载下的I和II型裂缝扩展问题、冲击荷载作用下的I型动态裂缝分叉和II型动态裂缝扩展问题，计算所得结果与已有文献和试验成果基本一致，说明了方法的可靠性；同时，对岩石力学领域中的双平行翼型裂缝相交和动态裂缝三维曲面扩展等问题，也利用该方法进行模拟。计算实例表明，结构储存的较高弹性应变能是脆性材料动态裂缝分叉的主要原因，而相场断裂模型是一种能够计算包括裂缝起裂、相交、分叉和三维曲面扩展的有效断裂计算方法，借助通用有限元平台可以方便地实现其分步算法。

	服务

	把本文推荐给朋友
	加入我的书架
	加入引用管理器
	E-mail Alert
	RSS
	作者相关文章
	刘国威1
	李庆斌1
	左正1
	2

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

Abstract：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.

Key words： rock mechanics ABAQUS UMAT/VUMAT phase field model staggered algorithm crack branching

引用本文:

刘国威1，李庆斌1，左正1，2. 相场断裂模型分步算法在ABAQUS中的实现[J]. 岩石力学与工程学报, 2016, 35(5): 1019-1030.
LIU Guowei1，LI Qingbin1，ZUO Zheng1，2. Implementation of a staggered algorithm for a phase field model in ABAQUS. , 2016, 35(5): 1019-1030.

链接本文:

http://www.rockmech.org/CN/10.13722/j.cnki.jrme.2015.1264 或 http://www.rockmech.org/CN/Y2016/V35/I5/1019

[28]	HESCH C，WEINBERG K. Thermodynamically consistent algorithms for a finite-deformation phase-field approach to fracture[J]. International Journal for Numerical Methods in Engineering，2014，99(12)：906-924. " target="_blank">
[9]	BELYTSCHKO T，BLACK T. Elastic crack growth in finite elements with minimal remeshing[J]. International Journal for Numerical Methods in Engineering，1999，45(5)：601-620. " target="_blank">
[2]	GORDON R. J，STRYK R. A. Eroding interface and improved tetrahedral element algorithms for high-velocity impact computations in three dimensions[J]. International Journal of Impact Engineering，1987，5(1/4)：411-421. " target="_blank">
[6]	XU X P，NEEDLEMAN A. Numerical simulations of fast crack growth in brittle solids[J]. Journal of the Mechanics and Physics of Solids，1994，42(9)：1 397-1 434. " target="_blank">
[12]	MENOUILLARD T，BELYTSCHKO T. Dynamic fracture with meshfree enriched XFEM[J]. Acta Mechanica，2010，213(1/2)：53-69. " target="_blank">
[17]	HENRY H，LEVINE H. Dynamic instabilities of fracture under biaxial strain using a phase field model[J]. Physical Review Letters，2004，93(10). DOI:http://dx.doi.org/10.1103/PhysRevLett.93.105504.
[20]	BOURDIN B，FRANCFORT G A，MARIGO J J. Numerical experiments in revisited brittle fracture[J]. Journal of the Mechanics and Physics of Solids，2000，48(4)：797-826. " target="_blank">
[22]	KUHN C，MULLER R. A continuum phase field model for fracture[J]. Engineering Fracture Mechanics，2010，77(18)：3 625-3 634. " target="_blank">
[23]	MIEHE C，HOFACKER M，WELSCHINGER F. A phase field model for rate-independent crack propagation：robust algorithmic implementation based on operator splits[J]. Computer Methods in Applied Mechanics and Engineering，2010，199(45/48)：2 765-2 778. " target="_blank">
[26]	SCHLÜTER A， WILLENBÜCHER A， KUHN C，et al. Phase field approximation of dynamic brittle fracture[J]. Computational Mechanics，2014，54(5)：1 141-1 161. " target="_blank">
[27]	MIEHE C，LAMBRECHT M. Algorithms for computation of stresses and elasticity moduli in terms of Seth-Hill family of generalized strain tensors[J]. Communications in Numerical Methods in Engineering，2001，17(5)：337-353. " target="_blank">
[29]	BORDEN M J，VERHOOSEL C V，SCOTT M A，et al. A phase-field description of dynamic brittle fracture[J]. Computer Methods in Applied Mechanics and Engineering，2012，217：77-95. " target="_blank">
[30]	LINDER C，ARMERO F. Finite elements with embedded branching[J]. Finite Elements in Analysis and Design，2009，45(4)：280-293.
[3]	FAN R，FISH J. The rs-method for material failure simulations[J]. International Journal for Numerical Methods in Engineering，2008，73(11)：1 607-1 623. " target="_blank">
[13]	XU D D，LIU Z L，LIU X M，et al. Modeling of dynamic crack branching by enhanced extended finite element method[J]. Computational Mechanics，2014，54(2)：489-502. " target="_blank">
[8]	SAM C H，PAPOULIA K D，VAVASIS S A. Obtaining initially rigid cohesive finite element models that are temporally convergent[J]. Engineering Fracture Mechanics，2005，72(14)：2 247-2 267. " target="_blank">
[18]	PONS J，KARMA A. Helical crack-front instability in mixed-mode fracture[J]. Nature，2010，464(7285)：85-89. " target="_blank">
[5]	SONG J H，WANG H W，BELYTSCHKO T. A comparative study on finite element methods for dynamic fracture[J]. Computational Mechanics，2008，42(2)：239-250.
[10]	MOES N，DOLBOW J，BELYTSCHKO T. A finite element method for crack growth without remeshing[J]. International Journal for Numerical Methods in Engineering，1999，46(1)：131-150. " target="_blank">
[16]	KARMA，KESSLER D A，Levine H. Phase-field model of mode III dynamic fracture[J]. Physical Review Letters，2001，87(4). DOI:http: //dx.doi.org/10.1103/PhysRevLett.87.045501
[24]	HOFACKER M，MIEHE C. A phase field model of dynamic fracture：Robust field updates for the analysis of complex crack patterns[J]. International Journal for Numerical Methods in Engineering，2013，93(3)：276-301.
[31]	FINEBERG J，MARDER M. Instability in dynamic fracture[J]. Physics Reports-Review Section of Physics Letters，1999，313(1/2)：1-108. " target="_blank">
[4]	LIU Y，FILONOVA V，HU N，et al. A regularized phenomenological multiscale damage model[J]. International Journal for Numerical Methods in Engineering，2014，99(12)：867-887. " target="_blank">
[14]	RASHID Y R. Ultimate strength analysis of prestressed concrete pressure vessels[J]. Nuclear Engineering and Design，1968，7(4)：334-344. " target="_blank">
[1]	BELYTSCHKO T，LIN J I. A 3-dimensional impact penertation algorithm with erosion[J]. Computers and Structures，1987，25(1)：95-104. " target="_blank">
[11]	BELYTSCHKO T，CHEN H，XU J X，et al. Dynamic crack propagation based on loss of hyperbolicity and a new discontinuous enrichment[J]. International Journal for Numerical Methods in Engineering，2003，58(12)：1 873-1 905. " target="_blank">
[15]	BAZANT Z P，CEDOLIN L. Finite-element modeling of crack band propagation[J]. Journal of Structural Engineering，ASCE，1983，109(1)：69-92. " target="_blank">
[7]	CAMACHO G T，ORTIZ M. Computational modelling of impact damage in brittle materials[J]. International Journal of Solids and Structures，1996，33(20/22)：2 899-2 938. " target="_blank">
[21]	MIEHE C，WELSCHINGER F，HOFACKER M. Thermodynamically consistent phase-field models of fracture：variational principles and multi-field FE implementations[J]. International Journal for Numerical Methods in Engineering，2010，83(10)：1 273-1 311. " target="_blank">
[25]	MSEKH M A，SARGADO J M，JAMSHIDIAN M，et al. Abaqus implementation of phase-field model for brittle fracture[J]. Computational Materials Science，2015，96(B)：472-484. " target="_blank">
[19]	FRANCFORT G A，MARIGO J J. Revisiting brittle fracture as an energy minimization problem[J]. Journal of the Mechanics and Physics of Solids，1998，46(8)：1 319-1 342. " target="_blank">