Recognition: no theorem link
Deep learning-based phase-field modelling of brittle fracture in anisotropic media
Pith reviewed 2026-05-15 07:07 UTC · model grok-4.3
The pith
A variational deep learning framework solves higher-order anisotropic phase-field fracture models by minimising total energy with B-spline enriched trial spaces.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
This work presents a variational physics-informed deep learning framework for phase-field modelling of brittle crack propagation in anisotropic media. For the first time in a variational deep learning setting it introduces a family of higher-order anisotropic phase-field models through a generalised crack density functional, solves the fracture problem by minimising the total energy with the Deep Ritz Method, and enriches the trial space with higher-order B-spline basis functions to represent higher-order gradients accurately and stably without conventional automatic differentiation. The methodology is assessed for isotropic, cubic and orthotropic fracture surface energy densities, and the 1
What carries the argument
Generalised crack density functional solved via Deep Ritz Method on a trial space enriched with higher-order B-spline basis functions
If this is right
- Direction-dependent crack growth is reproduced in cubic and orthotropic media without mesh-dependent artifacts.
- The same energy-minimisation procedure works for any fracture surface energy density that can be written as a generalised crack density functional.
- Higher-order gradients are handled stably, allowing phase-field models of arbitrary order to be treated within the Deep Ritz framework.
- The approach removes the requirement for automatic differentiation when computing variational derivatives of the energy.
Where Pith is reading between the lines
- The B-spline enrichment could be transferred to other physics-informed neural network problems that involve fourth- or higher-order differential operators.
- Because the method is mesh-free in spirit, it may scale more readily to three-dimensional anisotropic fracture problems than traditional finite-element phase-field codes.
- The generalised crack density functional provides a natural route to include additional physics such as thermal or electrical coupling inside the same variational deep learning setting.
Load-bearing premise
Higher-order B-spline basis functions in the trial space represent higher-order gradients accurately and stably, eliminating the need for automatic differentiation.
What would settle it
Run the method on a standard orthotropic plate with a central crack under uniaxial tension and compare the predicted crack path angle against the analytically expected direction; mismatch in the angle would falsify the claim that the enriched trial space captures anisotropic behaviour correctly.
Figures
read the original abstract
This work presents a variational physics-informed deep learning framework for phase-field modelling of brittle crack propagation in anisotropic media. Previous Deep Ritz Method (DRM) approaches have focused on second-order, isotropic phase-field fracture formulations. In contrast, the present work introduces, for the first time within a variational deep learning setting, a family of higher-order anisotropic phase-field models through a generalised crack density functional. The resulting fracture problem is solved by minimising the total energy using the DRM. The trial space is enriched with higher-order B-spline basis functions to represent higher-order gradients accurately and stably, thereby eliminating the need for conventional automatic differentiation. The methodology is assessed for isotropic, cubic, and orthotropic fracture surface energy densities. Numerical examples demonstrate direction-dependent crack growth in anisotropic cases, highlighting the capability of the method to accurately capture this behaviour.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a variational physics-informed deep learning framework based on the Deep Ritz Method for phase-field modeling of brittle fracture in anisotropic media. It introduces a generalized higher-order anisotropic crack density functional, solves the energy minimization problem via DRM, and enriches the trial space with higher-order B-spline basis functions to represent higher-order gradients without conventional automatic differentiation. The approach is assessed on isotropic, cubic, and orthotropic fracture surface energy densities, with numerical examples claimed to demonstrate direction-dependent crack growth.
Significance. If the stability and accuracy of the B-spline enrichment are rigorously established, the work would usefully extend DRM-based variational methods from isotropic second-order phase-field models to anisotropic higher-order cases. The avoidance of automatic differentiation via spline enrichment is a technically interesting choice that could improve robustness in energy minimization for fracture problems.
major comments (2)
- [Numerical examples] Numerical examples section: the central claim that the method 'accurately captures' direction-dependent crack growth in anisotropic cases is not supported by any reported error norms, convergence rates, or quantitative comparisons against reference solutions (e.g., established FEM implementations of the same anisotropic phase-field models). Only qualitative demonstration is described.
- [Method (B-spline enrichment)] Section describing the B-spline-enriched trial space: the assertion that higher-order B-splines 'represent higher-order gradients accurately and stably' and thereby eliminate the need for automatic differentiation lacks a priori stability estimates, dispersion analysis, or numerical convergence studies for the cubic/orthotropic crack-density functionals. This is load-bearing for the claim that the enriched space is suitable for the generalized anisotropic model.
minor comments (1)
- [Abstract] Abstract: the phrase 'for the first time within a variational deep learning setting' should be accompanied by a brief citation or clarification of the exact prior DRM works that are being extended.
Simulated Author's Rebuttal
We thank the referee for the constructive comments and positive assessment of the work's potential. We address each major comment below and will revise the manuscript to strengthen the quantitative support for our claims.
read point-by-point responses
-
Referee: [Numerical examples] Numerical examples section: the central claim that the method 'accurately captures' direction-dependent crack growth in anisotropic cases is not supported by any reported error norms, convergence rates, or quantitative comparisons against reference solutions (e.g., established FEM implementations of the same anisotropic phase-field models). Only qualitative demonstration is described.
Authors: We agree that quantitative validation is needed to support the accuracy claims. In the revised manuscript, we will add L2 error norms for the phase-field and displacement fields, crack path deviation metrics, and convergence rates with respect to network parameters and spline order, using established FEM solutions as references for the isotropic, cubic, and orthotropic cases. revision: yes
-
Referee: [Method (B-spline enrichment)] Section describing the B-spline-enriched trial space: the assertion that higher-order B-splines 'represent higher-order gradients accurately and stably' and thereby eliminate the need for automatic differentiation lacks a priori stability estimates, dispersion analysis, or numerical convergence studies for the cubic/orthotropic crack-density functionals. This is load-bearing for the claim that the enriched space is suitable for the generalized anisotropic model.
Authors: We acknowledge the value of additional analysis. While a full a priori theoretical stability proof is beyond the scope of this computational paper, the revised version will include numerical convergence studies, error tables, and dispersion analysis specifically for the B-spline enrichment on the cubic and orthotropic crack-density functionals to empirically demonstrate stability and accuracy. revision: partial
Circularity Check
No significant circularity; extends DRM with B-spline enrichment on independent numerical validation
full rationale
The derivation introduces a generalised anisotropic crack density functional and solves the energy minimisation via DRM with higher-order B-spline enrichment of the trial space. This enrichment is presented as a technical choice to avoid automatic differentiation, building on established DRM literature without reducing the central claim (accurate capture of direction-dependent crack growth in cubic/orthotropic cases) to a fitted parameter, self-definition, or self-citation chain. Numerical examples are used to demonstrate the behaviour rather than deriving it tautologically from the inputs. Any prior DRM citations are not load-bearing for the anisotropic extension or the stability claim, satisfying the criteria for a low circularity score with self-contained content against external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (1)
- B-spline order
axioms (1)
- domain assumption Minimization of total energy via DRM solves the anisotropic fracture problem
Reference graph
Works this paper leans on
- [1]
- [2]
-
[3]
Heon Lee and Sungho Jeon. An experimental and numerical study of fracture coalescence in pre- cracked specimens under compression. International Journal of Solids and Structures , 48(6):979– 999, 2011
work page 2011
- [4]
-
[5]
B. Bourdin, G. A. Francfort, and J.-J. Marigo. The Variational Approach to Fracture. Journal of Elasticity, 91(1-3):5–148, 2008
work page 2008
- [6]
-
[7]
Christian Miehe, Martina Hofacker, and Fabian Welschinger. A phase field model for rate- independent crack propagation: Robust algorithmic implementation based on operator splits. Computer Methods in Applied Mechanics and Engineering , 199(45–48):2765–2778, 2010
work page 2010
-
[8]
Michael J. Borden, Clemens V. Verhoosel, Michael A. Scott, Thomas J.R. Hughes, and Chad M. Landis. A phase-field description of dynamic brittle fracture. Computer Methods in Applied Mechanics and Engineering, 217–220:77–95, 2012
work page 2012
-
[9]
Tushar Kanti Mandal, Vinh Phu Nguyen, Jian-Ying Wu, Chi Nguyen-Thanh, and Alban de Vaucorbeil. Fracture of thermo-elastic solids: Phase-field modeling and new results with an efficient monolithic solver. Computer Methods in Applied Mechanics and Engineering, 376:113648, 2021
work page 2021
-
[10]
Jike Han, Seishiro Matsubara, Shuji Moriguchi, and Kenjiro Terada. Variational crack phase-field model for ductile fracture with elastic and plastic damage variables.Computer Methods in Applied Mechanics and Engineering, 400:115577, 2022
work page 2022
-
[11]
Chenyi Luo, Lorenzo Sanavia, and Laura De Lorenzis. Phase-field modeling of drying-induced cracks: Choice of coupling and study of homogeneous and localized damage. Computer Methods in Applied Mechanics and Engineering , 410:115962, 2023
work page 2023
-
[12]
Fan Fei, Andrea Costa, John E. Dolbow, Randolph R. Settgast, and Marco Cusini. A phase-field model for hydraulic fracture nucleation and propagation in porous media. International Journal for Numerical and Analytical Methods in Geomechanics , 47(16):3065–3089, 2023
work page 2023
-
[13]
Michael J. Borden, Thomas J.R. Hughes, Chad M. Landis, Amin Anvari, and Isaac J. Lee. A phase-field formulation for fracture in ductile materials: Finite deformation balance law deriva- tion, plastic degradation, and stress triaxiality effects. Computer Methods in Applied Mechanics and Engineering, 312:130–166, 2016. 30
work page 2016
-
[14]
Plohr, Gianmarco Manzini, and Hashem M
Lampros Svolos, JeeYeon N. Plohr, Gianmarco Manzini, and Hashem M. Mourad. On the convex- ity of phase-field fracture formulations: Analytical study and comparison of various degradation functions. International Journal of Non-Linear Mechanics , 154:104359, 2023
work page 2023
-
[15]
Nima Sarmadi, Mohaddeseh Mousavi Nezhad, and Quentin J. Fisher. On the numerical and mesh-dependent parameters in a computationally enhanced phase-field fracture model coupled with a novel mesh refinement strategy. Engineering with Computers , 39(6):3973–3992, 2023
work page 2023
-
[16]
Maziar Raissi, Paris Perdikaris, and George Em Karniadakis. Physics-informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations. Journal of Computational Physics , 378:686–707, 2019
work page 2019
-
[17]
A physics- informed deep learning framework for inversion and surrogate modelling in solid mechanics
Ehsan Haghighat, Maziar Raissi, Adrian Moure, Hector Gomez, and Ruben Juanes. A physics- informed deep learning framework for inversion and surrogate modelling in solid mechanics. Com- puter Methods in Applied Mechanics and Engineering , 379:113741, 2021
work page 2021
-
[18]
Physics-informed deep learning for computational elastodynamics without labeled data
Chengping Rao, Hao Sun, and Yang Liu. Physics-informed deep learning for computational elastodynamics without labeled data. Journal of Engineering Mechanics , 147(8):04021043, 2021
work page 2021
-
[19]
Critical investigation of failure modes in physics-informed neural networks
Shamsulhaq Basir and Inanc Senocak. Critical investigation of failure modes in physics-informed neural networks. AIAA SCITECH 2022 Forum , 2022
work page 2022
-
[20]
Transfer learn- ing enhanced physics informed neural network for phase-field modeling of fracture
Somdatta Goswami, Cosmin Anitescu, Souvik Chakraborty, and Timon Rabczuk. Transfer learn- ing enhanced physics informed neural network for phase-field modeling of fracture. Theoretical and Applied Fracture Mechanics, 106:102447, 2020
work page 2020
-
[21]
Somdatta Goswami, Cosmin Anitescu, and Timon Rabczuk. Adaptive fourth-order phase field analysis using deep energy minimization.Theoretical and Applied Fracture Mechanics, 107:102527, 2020
work page 2020
-
[22]
Hamdia, Xiaoying Zhuang, and Timon Rabczuk
Esteban Samaniego, Cosmin Anitescu, Somdatta Goswami, Vinh Phu Nguyen-Thanh, Haolin Guo, Karam S. Hamdia, Xiaoying Zhuang, and Timon Rabczuk. An energy approach to the solu- tion of partial differential equations in computational mechanics via machine learning: Concepts, implementation and applications. Computer Methods in Applied Mechanics and Engineerin...
work page 2020
-
[23]
The deep ritz method: A deep learning-based numerical algorithm for solving variational problems
Weinan E and Bing Yu. The deep ritz method: A deep learning-based numerical algorithm for solving variational problems. Communications in Mathematics and Statistics , 6(1):1–12, 2018
work page 2018
- [24]
-
[25]
T. C. T. Ting. Anisotropic Elasticity: Theory and Applications . Oxford University Press, New York, 1996
work page 1996
-
[26]
Hao Wang, Tian Ma, Yang Liu, Kai Zhang, and Xin Li. Experimental investigation on the 3D anisotropic fracture behavior of layered shales under Mode-I loading. Rock Mechanics and Rock Engineering, 57:3459–3482, 2024
work page 2024
-
[27]
Gang Bao, S. Ho, Zhigang Suo, and Bing Fan. The role of material orthotropy in fracture specimens for composites. International Journal of Solids and Structures , 29(9):1105–1116, 1992. 31
work page 1992
-
[28]
Gurtin and Paolo Podio-Guidugli
Morton E. Gurtin and Paolo Podio-Guidugli. Configurational forces and a constitutive theory for crack propagation that allows for kinking and curving. Journal of the Mechanics and Physics of Solids, 46(8):1343–1378, 1998
work page 1998
- [29]
-
[30]
Bin Li, Christian Peco, Daniel Mill´ an, Irene Arias, and Marino Arroyo. Phase-field modeling and simulation of fracture in brittle materials with strongly anisotropic surface energy. International Journal for Numerical Methods in Engineering , 102(3–4):711–727, 2015
work page 2015
-
[31]
E. G. Kakouris and S. P. Triantafyllou. Material Point Method for Crack Propagation in Anisotropic Media: A Phase Field Approach. Archive of Applied Mechanics, 88:287–316, 2018
work page 2018
-
[32]
Matthew Tancik, Pratul P. Srinivasan, Ben Mildenhall, Sara Fridovich-Keil, Nithin Ragha- van, Utkarsh Singhal, Ravi Ramamoorthi, Jonathan T. Barron, and Ren Ng. Fourier fea- tures let networks learn high frequency functions in low dimensional domains. arXiv preprint arXiv:2006.10739, 2020
-
[33]
Thomas Bachlechner, Bodhisattwa Prasad Majumder, Huanru Henry Mao, Garrison W. Cottrell, and Julian McAuley. Rezero is all you need: Fast convergence at large depth. arXiv preprint arXiv:2003.04887, 2020
-
[34]
Deep Residual Learning for Image Recognition
Kaiming He, Xiangyu Zhang, Shaoqing Ren, and Jian Sun. Deep residual learning for image recognition. arXiv preprint arXiv:1512.03385 , 2015
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[35]
Emmanouil G. Kakouris and Savvas P. Triantafyllou. Phase-Field Material Point Method for Dynamic Brittle Fracture with Isotropic and Anisotropic Surface Energy. Computer Methods in Applied Mechanics and Engineering , 357:112503, 2019
work page 2019
-
[36]
Javier Bonet and Richard D. Wood. Nonlinear Continuum Mechanics for Finite Element Analysis. Cambridge University Press, Cambridge, 2nd edition, 2008
work page 2008
- [37]
-
[38]
A review on phase-field models of brittle fracture and a new fast hybrid formulation
Marreddy Ambati, Tymofiy Gerasimov, and Laura De Lorenzis. A review on phase-field models of brittle fracture and a new fast hybrid formulation. Computational Mechanics, 55(2):383–405, 2015
work page 2015
-
[39]
On penalization in variational phase-field models of brittle fracture
Timur Gerasimov and Laura De Lorenzis. On penalization in variational phase-field models of brittle fracture. Computer Methods in Applied Mechanics and Engineering , 354:990–1026, 2019
work page 2019
-
[40]
Wheeler, Thomas Wick, and Wolfgang Wollner
Matthew F. Wheeler, Thomas Wick, and Wolfgang Wollner. An augmented-lagrangian method for the phase-field approach for pressurized fractures. Computer Methods in Applied Mechanics and Engineering, 271:69–85, 2014
work page 2014
-
[41]
Yousef Ghaffari Motlagh, Peter K. Jimack, and Ren´ e de Borst. Deep learning phase-field model for brittle fractures. International Journal for Numerical Methods in Engineering, 124(3):620–638, 2023
work page 2023
-
[42]
A Practical Guide to Splines, volume 27 ofApplied Mathematical Sciences
Carl de Boor. A Practical Guide to Splines, volume 27 ofApplied Mathematical Sciences. Springer, New York, revised edition edition, 2001. 32
work page 2001
-
[43]
Thomas J. R. Hughes, John A. Cottrell, and Yuri Bazilevs. Isogeometric analysis: Cad, finite elements, nurbs, exact geometry and mesh refinement. Computer Methods in Applied Mechanics and Engineering, 194(39–41):4135–4195, 2005
work page 2005
-
[44]
Ritukesh Bharali, Somdatta Goswami, Cosmin Anitescu, and Timon Rabczuk. A robust mono- lithic solver for phase-field fracture integrated with fracture energy based arc-length method and under-relaxation. Computer Methods in Applied Mechanics and Engineering , 394:114927, 2022
work page 2022
-
[45]
Tao Jin, Zhao Li, and Kuiying Chen. A novel phase-field monolithic scheme for brittle crack prop- agation based on the limited-memory bfgs method with adaptive mesh refinement. International Journal for Numerical Methods in Engineering , 125(22):e7572, 2024
work page 2024
-
[46]
On the Spectral Bias of Neural Networks
Nasim Rahaman, Aristide Baratin, Devansh Arpit, Felix Draxler, Min Lin, Fred A. Hamprecht, Yoshua Bengio, and Aaron Courville. On the spectral bias of neural networks. arXiv preprint arXiv:1806.08734, 2019
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[47]
Understanding the difficulty of training deep feedforward neural networks
Xavier Glorot and Yoshua Bengio. Understanding the difficulty of training deep feedforward neural networks. In Proceedings of the 13th International Conference on Artificial Intelligence and Statistics (AISTATS), volume 9 of Proceedings of Machine Learning Research, pages 249–256. PMLR, 2010
work page 2010
-
[48]
Johannes M¨ uller and Marius Zeinhofer. Deep ritz revisited. arXiv preprint arXiv:1912.03937 , 2020
-
[49]
Dying relu and initialization: Theory and numerical examples
Lu Lu, Yeonjong Shin, Yanhui Su, and George Em Karniadakis. Dying relu and initialization: Theory and numerical examples. Communications in Computational Physics , 28(5):1671–1706, 2020
work page 2020
-
[50]
N. Sukumar and Ankit Srivastava. Exact imposition of boundary conditions with distance func- tions in physics-informed deep neural networks. Computer Methods in Applied Mechanics and Engineering, 389:114333, 2022
work page 2022
-
[51]
Urb´ an, J´ erˆ ome Darbon, and George Em Karniadakis
Elham Kiyani, Khemraj Shukla, Jorge F. Urb´ an, J´ erˆ ome Darbon, and George Em Karniadakis. Optimizing the optimizer for physics-informed neural networks and kolmogorov-arnold networks. Computer Methods in Applied Mechanics and Engineering , 446:118308, 2025
work page 2025
-
[52]
Byrd, Jorge Nocedal, and Robert B
Richard H. Byrd, Jorge Nocedal, and Robert B. Schnabel. Representations of quasi-newton matrices and their use in limited memory methods. Mathematical Programming, 63(1–3):129–156, 1994
work page 1994
-
[53]
Adam: A method for stochastic optimization
Diederik P Kingma and Jimmy Ba. Adam: A method for stochastic optimization. arXiv, 2014
work page 2014
-
[54]
Nojus Plung˙ e, P. Brommer, R. S. Edwards, and E. G. Kakouris. Drm fracture jax, April 2026. https://doi.org/10.5281/zenodo.19553764. A Constitutive matrices for anisotropic elasticity The constitutive matrices are expressed in Voigt notation and are given in the material principal axes. The in-plane formulation corresponds to a two-dimensional setting un...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.