Adaptive fourth-order phase field analysis using deep energy minimization

Abstract Phase field modeling of fracture is computationally expensive as it demands a very fine mesh to resolve the damage region. Hence, the practical application of such models are severely limited. Local refinement techniques are often necessary. In our recent work (Goswami et al., 2019), for solving brittle fracture problems using physics informed neural network (PINN), the crack path is resolved by minimizing the variational energy of the system. However, in Goswami et al. (2019) we used a pre-refined domain based on prior information of the failure path, which is not always available. In this work, we propose an adaptive h -refinement scheme to locally refine the domain along the path of the growth of the crack. The phase field parameter, ϕ and a residual-based posteriori error estimator are the proposed convenient measures to determine the need for refinement. For ϕ , a critical threshold is chosen such that it is lower than the value at which crack nucleation occurs and the fracture region is easily identified. On the other hand, for the residual-based error estimator, elements contributing to the highest error are marked for refinement. The proposed algorithm takes as input the geometry described using NURBS patches. For the evaluation of the basis functions, we develop a procedure based on the Bezier representation and integrate it with the adaptive refinement formulation. The results obtained using the adaptive refinement integrated variational energy based PINN approach is validated with the available analytical solution for several examples from the literature. The proposed approach is implemented on several two and three-dimensional examples to illustrate the effectiveness of the formulation. Code and data necessary for replicating the results of the examples in the article will be made available through a GitHub repository.

[1]  John A. Evans,et al.  Isogeometric finite element data structures based on Bézier extraction of NURBS , 2011 .

[2]  Ted Belytschko,et al.  Cracking particles: a simplified meshfree method for arbitrary evolving cracks , 2004 .

[3]  Gilles A. Francfort,et al.  Revisiting brittle fracture as an energy minimization problem , 1998 .

[4]  E Weinan,et al.  The Deep Ritz Method: A Deep Learning-Based Numerical Algorithm for Solving Variational Problems , 2017, Communications in Mathematics and Statistics.

[5]  Fernando Fraternali,et al.  Eigenfracture: An Eigendeformation Approach to Variational Fracture , 2009, Multiscale Model. Simul..

[6]  Dominik Schillinger,et al.  Isogeometric collocation for phase-field fracture models , 2015 .

[7]  Yoshua Bengio,et al.  Understanding the difficulty of training deep feedforward neural networks , 2010, AISTATS.

[8]  T. Rabczuk,et al.  A Nonlocal Operator Method for Partial Differential Equations with Application to Electromagnetic Waveguide Problem , 2019, Computers, Materials & Continua.

[9]  Ratna Kumar Annabattula,et al.  A FEniCS implementation of the phase field method for quasi-static brittle fracture , 2018, Frontiers of Structural and Civil Engineering.

[10]  Timon Rabczuk,et al.  Adaptive phase field analysis with dual hierarchical meshes for brittle fracture , 2019, Engineering Fracture Mechanics.

[11]  Maziar Raissi,et al.  Deep Hidden Physics Models: Deep Learning of Nonlinear Partial Differential Equations , 2018, J. Mach. Learn. Res..

[12]  B. Bourdin,et al.  The Variational Approach to Fracture , 2008 .

[13]  S. Timoshenko,et al.  Theory of elasticity , 1975 .

[14]  Paris Perdikaris,et al.  Physics-informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations , 2019, J. Comput. Phys..

[15]  Anthony Gravouil,et al.  2D and 3D Abaqus implementation of a robust staggered phase-field solution for modeling brittle fracture , 2017 .

[16]  D. S. Dugdale.  Yielding of steel sheets containing slits , 1960 .

[17]  Jitesh H. Panchal,et al.  Simulator-free solution of high-dimensional stochastic elliptic partial differential equations using deep neural networks , 2019, J. Comput. Phys..

[18]  Y. Yoon,et al.  Multiscale failure analysis with coarse-grained micro cracks and damage , 2014 .

[19]  Paris Perdikaris,et al.  Physics-Constrained Deep Learning for High-dimensional Surrogate Modeling and Uncertainty Quantification without Labeled Data , 2019, J. Comput. Phys..

[20]  Ronald Krueger,et al.  The Virtual Crack Closure Technique : History , Approach and Applications , 2002 .

[21]  Timon Rabczuk,et al.  Dual‐horizon peridynamics , 2015, 1506.05146.

[22]  Timon Rabczuk,et al.  Transfer learning enhanced physics informed neural network for phase-field modeling of fracture , 2019, Theoretical and Applied Fracture Mechanics.

[23]  Tae-Yeon Kim,et al.  Strong form-based meshfree collocation method for wind-driven ocean circulation , 2019, Computer Methods in Applied Mechanics and Engineering.

[24]  J. Michopoulos,et al.  Strong-Form Collocation Method for Solidification and Mechanical Analysis of Polycrystalline Materials , 2019, Journal of Engineering Mechanics.

[25]  Cv Clemens Verhoosel,et al.  A phase-field description of dynamic brittle fracture , 2012 .

[26]  Thomas J. R. Hughes,et al.  Isogeometric Analysis: Toward Integration of CAD and FEA , 2009 .

[27]  Les A. Piegl,et al.  The NURBS Book , 1995, Monographs in Visual Communication.

[28]  John A. Evans,et al.  Isogeometric boundary element analysis using unstructured T-splines , 2013 .

[29]  Thomas J. R. Hughes,et al.  A higher-order phase-field model for brittle fracture: Formulation and analysis within the isogeometric analysis framework , 2014 .

[30]  T. Rabczuk,et al.  A Deep Collocation Method for the Bending Analysis of Kirchhoff Plate , 2021, Computers, Materials & Continua.

[31]  Young-Cheol Yoon,et al.  Extended particle difference method for weak and strong discontinuity problems: part I. Derivation of the extended particle derivative approximation for the representation of weak and strong discontinuities , 2013, Computational Mechanics.

[32]  T. Rabczuk,et al.  T-spline based XIGA for fracture analysis of orthotropic media , 2015 .

[33]  Timon Rabczuk,et al.  An adaptive isogeometric analysis collocation method with a recovery-based error estimator , 2019, Computer Methods in Applied Mechanics and Engineering.

[34]  Timon Rabczuk,et al.  Adaptive fourth-order phase field analysis for brittle fracture , 2020 .

[35]  Ted Belytschko,et al.  Multiscale aggregating discontinuities: A method for circumventing loss of material stability , 2008 .

[36]  Y. Yoon,et al.  Extended particle difference method for moving boundary problems , 2014, Computational Mechanics.

[37]  Christian Miehe,et al.  A phase field model for rate-independent crack propagation: Robust algorithmic implementation based on operator splits , 2010 .

[38]  T. Rabczuk,et al.  Phase field simulations of coupled microstructure solidification problems via the strong form particle difference method , 2017, International Journal of Mechanics and Materials in Design.

[39]  Arnulf Jentzen,et al.  Solving high-dimensional partial differential equations using deep learning , 2017, Proceedings of the National Academy of Sciences.

[40]  B. Bourdin,et al.  Numerical experiments in revisited brittle fracture , 2000 .

[41]  A. A. Griffith.  The Phenomena of Rupture and Flow in Solids , 1921 .

[42]  Christian Miehe,et al.  Thermodynamically consistent phase‐field models of fracture: Variational principles and multi‐field FE implementations , 2010 .