Data Driven Prognosis of Fracture Dynamics Using Tensor Train and Gaussian Process Regression

Accurate prediction of the mechanical failure of structural components plays an important role in the design of engineering structures. However, the fracture process is challenging to model numerically due to the existence of an elastic singularity at the crack tip. The phase field model (PFM) is one promising approach for modeling brittle fracture using an auxiliary field variable to regularize discontinuities associated with sharp cracks. Unfortunately, it is generally computationally expensive due to the need to solve a fourth-order partial differential equation. In order to reduce the computational burden for parameterized problems, different reduced order modeling approaches such as proper orthogonal decomposition (POD) and discrete empirical interpolation method (DEIM), have been proposed. However, these frameworks are model-based and intrusive approaches and need in-depth efforts in modifying existing complex simulation codes. In this work, a tensor train (TT) approach is proposed in combination with Gaussian process regression (GPR) for modeling and predicting the dynamics of fracture in composite materials. As opposed to POD and DEIM, the TT-GPR approach is a fully data-driven and non-intrusive approach. In this work, we study the fracture of a brittle elastic 2D rectangular slab with a pre-existing crack under Mode I (crack opening) loading conditions. The high dimensional training data for the TT model was generated using PFM with the finite difference method. The predictions by the TT-GPR model are compared with the results from the finite difference method. The TT-GPR is robust enough to predict the catastrophic evolution of the crack growth trend with an accuracy of up to 95% with computational load reduced by two orders of magnitude.

[1]  Xian-Da Zhang,et al.  Matrix Analysis and Applications , 2017 .

[2]  Stefan Volkwein,et al.  Galerkin Proper Orthogonal Decomposition Methods for a General Equation in Fluid Dynamics , 2002, SIAM J. Numer. Anal..

[3]  S. B. Biner Programming Phase-Field Modeling , 2017 .

[4]  Daniel Kressner,et al.  Exploiting Hidden Structure in Matrix Computations: Algorithms and Applications , 2016 .

[5]  R. Murray,et al.  Model reduction for compressible flows using POD and Galerkin projection , 2004 .

[6]  M. Hinze,et al.  Proper Orthogonal Decomposition Surrogate Models for Nonlinear Dynamical Systems: Error Estimates and Suboptimal Control , 2005 .

[7]  Thi-Hong-Hieu Le,et al.  Mesh-Based and Meshfree Reduced Order Phase-Field Models for Brittle Fracture: One Dimensional Problems , 2019, Materials.

[8]  Boris N. Khoromskij,et al.  Tensor Numerical Methods in Quantum Chemistry , 2018 .

[9]  Boris N. Khoromskij,et al.  Tensor-Structured Galerkin Approximation of Parametric and Stochastic Elliptic PDEs , 2011, SIAM J. Sci. Comput..

[10]  Carl E. Rasmussen,et al.  Gaussian processes for machine learning , 2005, Adaptive computation and machine learning.

[11]  T. Anderson,et al.  Fracture mechanics - Fundamentals and applications , 2017 .

[12]  Wr Graham,et al.  OPTIMAL CONTROL OF VORTEX SHEDDING USING LOW-ORDER MODELS. PART I-OPEN-LOOP MODEL DEVELOPMENT , 1999 .

[13]  S. Ravindran A reduced-order approach for optimal control of fluids using proper orthogonal decomposition , 2000 .

[14]  I.B.C.M. Rocha,et al.  Accelerating crack growth simulations through adaptive model order reduction , 2020, International Journal for Numerical Methods in Engineering.

[15]  Ann E. Jeffers,et al.  Feature-preserving rational Bézier triangles for isogeometric analysis of higher-order gradient damage models , 2019 .

[16]  Pham Luu Trung Duong,et al.  Application of multi-output Gaussian process regression for remaining useful life prediction of light emitting diodes , 2018, Microelectron. Reliab..

[17]  Ngai Wong,et al.  Tensor Network alternating linear scheme for MIMO Volterra system identification , 2016, Autom..

[18]  K. Kunisch,et al.  Control of the Burgers Equation by a Reduced-Order Approach Using Proper Orthogonal Decomposition , 1999 .

[19]  J. Lua,et al.  A simplified continuum damage mechanics based modeling strategy for cumulative fatigue damage assessment of metallic bolted joints , 2020 .

[21]  Alain Karma,et al.  Helical crack-front instability in mixed-mode fracture , 2010, Nature.

[22]  Ted Belytschko,et al.  Elastic crack growth in finite elements with minimal remeshing , 1999 .

[23]  Pham Luu Trung Duong,et al.  Tensor Train Decomposition for Data-Driven Prognosis of Fracture Dynamics in Composite Materials , 2020, 2020 IEEE Aerospace Conference.

[24]  Juš Kocijan,et al.  Modelling and Control of Dynamic Systems Using Gaussian Process Models , 2015 .

[25]  A. Karma,et al.  Phase-field model of mode III dynamic fracture. , 2001, Physical review letters.

[26]  Ivan Oseledets,et al.  Tensor-Train Decomposition , 2011, SIAM J. Sci. Comput..