pith. machine review for the scientific record. sign in

arxiv: 2603.20120 · v2 · submitted 2026-03-20 · ⚛️ physics.comp-ph

Recognition: no theorem link

Deep learning-based phase-field modelling of brittle fracture in anisotropic media

Authors on Pith no claims yet

Pith reviewed 2026-05-15 07:07 UTC · model grok-4.3

classification ⚛️ physics.comp-ph
keywords phase-field modellingbrittle fractureanisotropic mediaDeep Ritz MethodB-spline basis functionsvariational deep learningcrack propagationphysics-informed neural networks
0
0 comments X

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.

The paper introduces a physics-informed deep learning approach to phase-field modelling of brittle fracture that extends prior Deep Ritz Method work to anisotropic media. It defines a generalised crack density functional that incorporates higher-order anisotropic fracture surface energies and minimises the resulting total energy functional. The trial space is built from higher-order B-spline basis functions so that higher-order gradients are represented accurately and stably without automatic differentiation. Numerical tests on isotropic, cubic and orthotropic cases show that the method reproduces the expected direction-dependent crack paths.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2603.20120 by E. G. Kakouris, N. Plung\.e, P. Brommer, R. S. Edwards.

Figure 1
Figure 1. Figure 1: (a) Solid body Ω with crack path Γ and (b) phase field approximation of the crack path Γ. [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: (a) Polar representation of the fracture surface energy density [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Schematic illustration of the DRM employed in this work. Spatial coordinates are input [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: (a) Effect of the scaling coefficient r on the tanh activation function. Smaller values extend the quasi-linear regime. (b) Non-smooth piecewise constraint used to enforce admissible phase-field values. In classical finite element formulations, these effects may be alleviated through mesh refinement. In NN-based approximations, however, refinement of the discretisation is not directly applicable. Instead, … view at source ↗
Figure 5
Figure 5. Figure 5: Displacement field enforcement via distance functions. The shaded regions represent non [PITH_FULL_IMAGE:figures/full_fig_p015_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Square plate under pure tension: geometry and boundary conditions: (a) homogeneous plate [PITH_FULL_IMAGE:figures/full_fig_p020_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: (a) Polar representation of the fracture surface energy density [PITH_FULL_IMAGE:figures/full_fig_p020_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Results for the isotropic symmetry case. (a) Evolution of elastic and fracture energies [PITH_FULL_IMAGE:figures/full_fig_p021_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Cubic material symmetry for two material orientations. Material orientation [PITH_FULL_IMAGE:figures/full_fig_p022_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Orthotropic material symmetry for two material orientations. Material orientation [PITH_FULL_IMAGE:figures/full_fig_p024_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Orthotropic material symmetry of −30◦ tested with network architecture introduced by Manav et al. for second-order fracture problems: (a) variation of elastic and fracture energies with applied displacement, (b) phase-field distribution obtained using the finite element method (FEM), and (c) phase-field distribution obtained using the proposed DRM at a displacement of 0.015 m. The fracture trajectory is g… view at source ↗
Figure 12
Figure 12. Figure 12: Results for the layered orthotropic material case. (a) Evolution of elastic and fracture [PITH_FULL_IMAGE:figures/full_fig_p026_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Displacement step convergence study. (a) Mean absolute error and standard deviation [PITH_FULL_IMAGE:figures/full_fig_p027_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Mesh convergence study. (a) Mean absolute error and standard deviation of elastic energy, [PITH_FULL_IMAGE:figures/full_fig_p028_14.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

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)
  1. [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.
  2. [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)
  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

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

1 free parameters · 1 axioms · 0 invented entities

Abstract-only review; specific functional forms, hyperparameters, and training details are not provided, so ledger entries are limited to explicitly invoked assumptions.

free parameters (1)
  • B-spline order
    Chosen to represent higher-order gradients accurately and stably
axioms (1)
  • domain assumption Minimization of total energy via DRM solves the anisotropic fracture problem
    Central variational principle assumed to hold for the generalised crack density functional

pith-pipeline@v0.9.0 · 5450 in / 1068 out tokens · 35038 ms · 2026-05-15T07:07:57.203432+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

54 extracted references · 54 canonical work pages · 2 internal anchors

  1. [1]

    Sharon, S

    E. Sharon, S. P. Gross, and J. Fineberg. Local crack branching as a mechanism for instability in dynamic fracture. Physical Review Letters, 74(25):5096–5099, 1995

  2. [2]

    Sharon, S

    E. Sharon, S. P. Gross, and J. Fineberg. Microbranching instability and the dynamic fracture of brittle materials. Physical Review B, 54(10):7128–7139, 1996

  3. [3]

    An experimental and numerical study of fracture coalescence in pre- cracked specimens under compression

    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

  4. [4]

    Anderson

    Ted L. Anderson. Fracture Mechanics: Fundamentals and Applications. CRC Press, Boca Raton, FL, 4th edition, 2017

  5. [5]

    Bourdin, G

    B. Bourdin, G. A. Francfort, and J.-J. Marigo. The Variational Approach to Fracture. Journal of Elasticity, 91(1-3):5–148, 2008

  6. [6]

    Griffith

    Alan A. Griffith. The Phenomena of Rupture and Flow in Solids. Philosophical Transactions of the Royal Society A , 221:163–198, 1921

  7. [7]

    A phase field model for rate- independent crack propagation: Robust algorithmic implementation based on operator splits

    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

  8. [8]

    Borden, Clemens V

    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

  9. [9]

    Fracture of thermo-elastic solids: Phase-field modeling and new results with an efficient monolithic solver

    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

  10. [10]

    Variational crack phase-field model for ductile fracture with elastic and plastic damage variables.Computer Methods in Applied Mechanics and Engineering, 400:115577, 2022

    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

  11. [11]

    Phase-field modeling of drying-induced cracks: Choice of coupling and study of homogeneous and localized damage

    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

  12. [12]

    Dolbow, Randolph R

    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

  13. [13]

    Borden, Thomas J.R

    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

  14. [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

  15. [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

  16. [16]

    Physics-informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations

    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

  17. [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

  18. [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

  19. [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

  20. [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

  21. [21]

    Adaptive fourth-order phase field analysis using deep energy minimization.Theoretical and Applied Fracture Mechanics, 107:102527, 2020

    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

  22. [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...

  23. [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

  24. [24]

    Manav, R

    M. Manav, R. Molinaro, S. Mishra, and L. De Lorenzis. Phase-field modeling of fracture with physics-informed deep learning. Computer Methods in Applied Mechanics and Engineer- ing, 429:117104, 2024

  25. [25]

    T. C. T. Ting. Anisotropic Elasticity: Theory and Applications . Oxford University Press, New York, 1996

  26. [26]

    Experimental investigation on the 3D anisotropic fracture behavior of layered shales under Mode-I loading

    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

  27. [27]

    Ho, Zhigang Suo, and Bing Fan

    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

  28. [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

  29. [29]

    Takei, B

    A. Takei, B. Roman, J. Bico, E. Hamm, and F. Melo. Forbidden directions for the fracture of thin anisotropic sheets: An analogy with the Wulff plot. Physical Review Letters, 110(14):144301, 2013

  30. [30]

    Phase-field modeling and simulation of fracture in brittle materials with strongly anisotropic surface energy

    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

  31. [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

  32. [32]

    Srinivasan, Ben Mildenhall, Sara Fridovich-Keil, Nithin Ragha- van, Utkarsh Singhal, Ravi Ramamoorthi, Jonathan T

    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. [33]

    Cottrell, and Julian McAuley

    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. [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

  35. [35]

    Kakouris and Savvas P

    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

  36. [36]

    Javier Bonet and Richard D. Wood. Nonlinear Continuum Mechanics for Finite Element Analysis. Cambridge University Press, Cambridge, 2nd edition, 2008

  37. [37]

    Kakouris

    Emmanouil G. Kakouris. Linear Elastic Fracture Mechanics via the Material Point Method: A Phase Field Approach. Phd thesis, University of Nottingham, Nottingham, UK, 2019

  38. [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

  39. [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

  40. [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

  41. [41]

    Jimack, and Ren´ e de Borst

    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

  42. [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

  43. [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

  44. [44]

    A robust mono- lithic solver for phase-field fracture integrated with fracture energy based arc-length method and under-relaxation

    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

  45. [45]

    A novel phase-field monolithic scheme for brittle crack prop- agation based on the limited-memory bfgs method with adaptive mesh refinement

    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

  46. [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

  47. [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

  48. [48]

    Deep ritz revisited

    Johannes M¨ uller and Marius Zeinhofer. Deep ritz revisited. arXiv preprint arXiv:1912.03937 , 2020

  49. [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

  50. [50]

    Sukumar and Ankit Srivastava

    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

  51. [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

  52. [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

  53. [53]

    Adam: A method for stochastic optimization

    Diederik P Kingma and Jimmy Ba. Adam: A method for stochastic optimization. arXiv, 2014

  54. [54]

    Brommer, R

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