pith. sign in

arxiv: 2606.27177 · v1 · pith:LHKCTHJTnew · submitted 2026-06-25 · 💻 cs.CE

A hybrid IFENN solver for generalizable modeling of phase-field fracture initiation and propagation

Pith reviewed 2026-06-26 01:47 UTC · model grok-4.3

classification 💻 cs.CE
keywords phase-field fracturehybrid solverfinite element neural networkcrack initiationcrack propagationphysics-informed traininggeneralizable geometries
0
0 comments X

The pith

The IFENN hybrid couples a FEM solver to neural networks that approximate the phase-field equation, enabling fracture predictions on unseen geometries after one training run on a benchmark shape.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that a hybrid solver can capture the complete evolution of phase-field fracture, from the start of a crack through its growth, on shapes the model has never encountered before. It achieves this by training two networks once on a single benchmark geometry with a physics-informed loss that uses only the maximum strain energy and the phase-field variable at a small number of strategically chosen points. One network handles crack initiation while the other manages propagation, and artificial boundary conditions keep far-field values near zero during inference. A sympathetic reader would care because the method promises to cut the cost of repeated full simulations when the same trained networks are reused across different component shapes.

Core claim

The IFENN framework tightly couples a standard finite element solver for mechanical equilibrium with a pre-trained neural network that approximates the phase-field diffusion equation; a DeepOKAN network is used for the initiation stage and a CNN for the propagation stage, both trained physics-informed on one benchmark geometry with few increments and limited Gauss points sampled from the fracture process zone, then applied to both the training geometry and arbitrary unseen geometries.

What carries the argument

The Integrated Finite Element Neural Network (IFENN) hybrid scheme, which pairs a conventional FEM solver for equilibrium with a neural network surrogate for the phase-field equation.

If this is right

  • The method models both crack initiation and subsequent propagation within a single framework.
  • Training cost drops sharply because only a small number of increments and Gauss points from the fracture zone are required.
  • The same trained networks produce usable results on both the original benchmark and previously unseen geometries.
  • Artificial boundary conditions allow the networks to extrapolate near-zero phase-field values far from the crack tip.

Where Pith is reading between the lines

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

  • Engineering workflows that repeatedly analyze fracture in families of similar parts could replace many full simulations with a single training step plus fast inference.
  • The same hybrid coupling pattern could be tested on other coupled-field problems such as diffusion-reaction systems or thermo-mechanical problems.
  • Systematic checks on geometries whose topology differs markedly from the training shape would reveal the practical limits of the claimed generalizability.

Load-bearing premise

A neural network trained exclusively on one benchmark geometry with artificial boundary conditions will produce accurate phase-field predictions on arbitrary unseen geometries without retraining or accumulating large errors.

What would settle it

A side-by-side comparison on a new geometry with a different topology, measuring whether the hybrid solver's predicted phase-field values and crack paths deviate substantially from a full finite-element reference solution.

Figures

Figures reproduced from arXiv: 2606.27177 by Diab Abueidda, Fouad Amin, Mostafa E. Mobasher, Panos Pantidis.

Figure 1
Figure 1. Figure 1: Schematic of a cracked domain with a sample phase-field damage contour [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Architecture of the adopted DeepOKAN, with Radial Basis Functions (RBFs) as the activation functions for all the [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: a. First sampling method (SM1), where colored Gauss points indicate sensor locations across zones A–D. b. Second sampling method (SM2), featuring uniformly and randomly selected Gauss points. c. Representative profiles of H and ϕ, illustrating the information contained within each zone. 8 [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Schematic representation of the weight-function adopted for filtering the far-field [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: d maps the selected training sensors against their aligned testing counterparts. The spatial proximity between each matched pair is evident. As it will be shown in the numerical examples, this approach enables the DeepOKAN to readily predict the localized phase-field initiation on completely unseen geometries, without requiring re-training or architectural modifications [PITH_FULL_IMAGE:figures/full_fig_p… view at source ↗
Figure 6
Figure 6. Figure 6: Illustration of the physics-informed CNN setup and training procedure. [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Schematic representations of a. train geometry (Model A) and b. test geometries (Model B). Once the FEM analysis is complete, we proceed with the one-time training of the networks. For the initiation stage we develop two variants of the DeepOKAN, one that is trained on TS1 and one on TS2. These are henceforth termed as DK-TS1 and DK-TS2. Both networks follow the general architecture of [PITH_FULL_IMAGE:fi… view at source ↗
Figure 8
Figure 8. Figure 8: Training performance of DK-TS1. a. True ϕ values. b. Predicted ϕ values. c. Absolute error ϕ values. d. Evolution of training loss. e. Metrics of ϕ, including its L2 norm and maximum value. The results of the DK-TS1 and DK-TS2 training are depicted in Figs. 8 and 9 respectively. In both figures, the top row displays the true ϕ profiles, the middle row presents the corresponding network predictions, and the… view at source ↗
Figure 9
Figure 9. Figure 9: Training performance of DK-TS2. a. True ϕ values. b. Predicted ϕ values. c. Absolute error ϕ values. d. Evolution of training loss. e. Metrics of ϕ, including its L2 norm and maximum value. Next, we develop the CNN model to capture the propagation stage of phase-field. The network consists of four consecutive convolution layers, with subsequent activation functions as presented in Section 3.2. The H profil… view at source ↗
Figure 10
Figure 10. Figure 10: Training performance of the developed CNN. The true and predicted [PITH_FULL_IMAGE:figures/full_fig_p018_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Numerical results: Case 1. a. Reaction force-displacement curves. b. Evolution of phase-field second norm. c. 3D views of phase-field contours at various initiation increments from FEM (top) and IFENN (bottom). Next, we demonstrate the capability of the hybrid analysis to operate within a multi-pass staggered scheme, again on Model A. The numerical tolerance is set to tol = 10−3 , defining the threshold b… view at source ↗
Figure 12
Figure 12. Figure 12: Numerical results: Case 2. a. Reaction force-displacement curves. b. Evolution of phase-field second norm. c. 3D views of phase-field contours at various initiation increments from FEM (top) and IFENN (bottom). Having demonstrated the performance of IFENN under the weight-function approach, we now evaluate its implementation utilizing a DeepOKAN that is trained via the artificial boundary condition method… view at source ↗
Figure 13
Figure 13. Figure 13: Numerical results: Case 3. a. Reaction force-displacement curves. b. Evolution of phase-field second norm. c. 3D views of phase-field contours at various initiation increments from FEM (top) and IFENN (bottom). So far, all IFENN analyses have been conducted on Model A, which is the same geometry that was used for the training of both networks. In this final example, we extend the applicability of these pr… view at source ↗
Figure 14
Figure 14. Figure 14: Numerical results: Case 4. a. Reaction force-displacement curves. b. Evolution of phase-field second norm. c. 3D views of phase-field contours at various initiation increments from FEM (top) and IFENN (bottom). 6. Summary, conclusions and outlook In this work, we present for the first time a complete IFENN framework for capturing the entire evolution of phase-field fracture in quasi-brittle materials. Bui… view at source ↗
read the original abstract

In this paper we demonstrate how the Integrated Finite Element Neural Network (IFENN) framework can effectively model the entire evolution of phase-field fracture, including the initiation and propagation stage, across generalizable geometries. IFENN is a hybrid scheme for coupled computational mechanics problems, tightly coupling a standard FEM solver (mechanical equilibrium) with a pre-trained neural network (coupled field). In this work, the phase-field diffusion equation is approximated with: i) a DeepONet architecture with Kolmogorov-Arnold networks in the trunk and branch (DeepOKAN) for the initiation stage, and ii) a Convolution Neural Network (CNN) for the propagation stage. Both networks are trained only once, on a benchmark geometry, using a purely physics-informed approach based on the maximum strain energy and the phase-field variable. The training process utilizes an extremely small number of training increments and only a limited number of Gauss points that are strategically sampled from the fracture process zone. These features enable a substantial decrease of the offline training cost. To address the extrapolation of the DeepOKAN predictions in regions away from the crack tip during the inference stage, we implement a set of artificial boundary conditions to enforce the near-zero values in the far-field predictions. We showcase the flexibility and numerical accuracy of the proposed methodology across both the training and unseen geometries.

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

3 major / 2 minor

Summary. The manuscript introduces a hybrid IFENN framework that couples a standard FEM solver for mechanical equilibrium with pre-trained neural networks to solve the phase-field diffusion equation for fracture. DeepOKAN (Kolmogorov-Arnold networks in trunk and branch) handles the initiation stage and a CNN handles propagation; both are trained once on a single benchmark geometry via a physics-informed loss based on maximum strain energy evaluated at a small number of strategically sampled Gauss points inside the process zone. Artificial far-field boundary conditions are added during inference to enforce near-zero predictions away from the crack tip. The central claim is that this yields accurate modeling of the full fracture evolution (initiation and propagation) on both the training geometry and arbitrary unseen geometries.

Significance. If the generalizability claim is substantiated with quantitative validation, the work would demonstrate a practical route to low-cost, one-time training of hybrid physics-informed networks for phase-field fracture on arbitrary domains. The emphasis on extremely limited training increments and Gauss points, together with the explicit handling of extrapolation via artificial BCs, addresses a known cost bottleneck in data-driven fracture modeling and could be extended to other coupled mechanics problems.

major comments (3)
  1. [Abstract] Abstract and numerical-results section: the central claim of 'numerical accuracy' and 'flexibility' across unseen geometries is unsupported by any reported quantitative metrics (L2 errors on the phase field, crack-path deviation, global energy balance, or convergence under mesh refinement). Without these, the generalizability assertion cannot be evaluated.
  2. [Method] Method and inference description: the artificial boundary conditions introduced to correct DeepOKAN extrapolation away from the crack tip are presented as a fix, yet no analysis quantifies their effect on local accuracy near the process zone or on geometries whose far-field conditions differ from the artificial prescription.
  3. [Numerical examples] Training protocol: the networks are trained exclusively on one benchmark geometry with a fixed, small set of Gauss points; the manuscript provides no description of the number, geometric diversity, or loading conditions of the 'unseen' test cases used to demonstrate generalizability, which is load-bearing for the primary claim.
minor comments (2)
  1. [Introduction] Notation for the phase-field variable and the maximum-strain-energy loss should be defined explicitly at first use rather than assumed from prior literature.
  2. [Figures] Figure captions should state the number of Gauss points used and the precise locations of the artificial boundary conditions so that the limited-data regime is reproducible.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive feedback on our manuscript. We address each major comment below with explanations and indicate the revisions we will make to strengthen the quantitative support for our claims and the clarity of the presentation.

read point-by-point responses
  1. Referee: [Abstract] Abstract and numerical-results section: the central claim of 'numerical accuracy' and 'flexibility' across unseen geometries is unsupported by any reported quantitative metrics (L2 errors on the phase field, crack-path deviation, global energy balance, or convergence under mesh refinement). Without these, the generalizability assertion cannot be evaluated.

    Authors: We acknowledge that the current manuscript primarily demonstrates accuracy and flexibility through visual comparisons of phase-field contours and crack paths against reference FEM solutions. While these comparisons show close qualitative agreement on both the training geometry and the unseen cases, we agree that explicit quantitative metrics are needed to rigorously support the central claims. In the revised manuscript we will add L2 errors on the phase-field variable, crack-path deviation metrics, global energy balance checks, and a brief mesh-convergence study for the key examples. revision: yes

  2. Referee: [Method] Method and inference description: the artificial boundary conditions introduced to correct DeepOKAN extrapolation away from the crack tip are presented as a fix, yet no analysis quantifies their effect on local accuracy near the process zone or on geometries whose far-field conditions differ from the artificial prescription.

    Authors: The artificial far-field boundary conditions are introduced to enforce the expected near-zero phase-field values outside the process zone during inference. We recognize that the manuscript does not include a dedicated quantification of their influence. We will add a short analysis in the revised method section that compares DeepOKAN predictions with and without these boundary conditions on the benchmark geometry (reporting local errors inside the process zone) and will verify robustness on one additional unseen geometry whose far-field setup differs from the artificial prescription. revision: yes

  3. Referee: [Numerical examples] Training protocol: the networks are trained exclusively on one benchmark geometry with a fixed, small set of Gauss points; the manuscript provides no description of the number, geometric diversity, or loading conditions of the 'unseen' test cases used to demonstrate generalizability, which is load-bearing for the primary claim.

    Authors: The unseen geometries are shown in the numerical examples, yet we agree that a concise, systematic description of their number, geometric variations, and loading conditions is currently insufficient. In the revised manuscript we will expand Section 4 with a table listing each unseen case, its geometric parameters (e.g., notch location, specimen aspect ratio), and the applied boundary/loading conditions, thereby making the generalizability evaluation fully reproducible. revision: yes

Circularity Check

0 steps flagged

No circularity; physics-informed training derives from governing PDEs independent of model outputs

full rationale

The paper trains DeepOKAN and CNN components via physics-informed losses on the phase-field diffusion equation residuals using maximum strain energy, applied to one benchmark geometry with limited Gauss points. This setup follows standard PINN methodology and does not reduce any prediction to a fitted parameter or self-referential definition. Artificial boundary conditions are introduced explicitly as an engineering patch for extrapolation, not as a definitional closure. No load-bearing self-citations, uniqueness theorems, or ansatzes imported from prior author work are described in the provided text. The generalization claim to unseen geometries is presented as an empirical outcome rather than a tautological reduction of inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on the standard phase-field fracture formulation and finite-element discretization; the neural-network approximations are the primary addition and are trained against the physics residual.

axioms (2)
  • domain assumption The phase-field model with maximum strain energy criterion accurately captures fracture initiation and propagation
    Invoked as the training target for both networks.
  • domain assumption Neural networks can approximate solutions to the phase-field diffusion equation to sufficient accuracy for engineering use
    Basis for replacing the FEM solve of the diffusion equation with DeepOKAN and CNN.

pith-pipeline@v0.9.1-grok · 5778 in / 1244 out tokens · 32954 ms · 2026-06-26T01:47:20.913162+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

45 extracted references

  1. [1]

    J. N. Reddy, An introduction to the finite element method, in: Dynamics of Earth’s Fluid System, CRC Press, 2026, pp. 199–226

  2. [2]

    T. J. Hughes, The finite element method: linear static and dynamic finite element analysis, Courier Corporation, 2003

  3. [3]

    G. J. Kennedy, J. R. Martins, A parallel finite-element framework for large-scale gradient-based design optimization of high-performance structures, Finite Elements in Analysis and Design 87 (2014) 56–73

  4. [4]

    P. K. Kristensen, E. Martínez-Pañeda, Phase field fracture modelling using quasi-newton methods and a new adaptive step scheme, Theoretical and Applied Fracture Mechanics 107 (2020) 102446

  5. [5]

    Neiva, S

    E. Neiva, S. Badia, A. F. Martín, M. Chiumenti, A scalable parallel finite element framework for growing geometries. application to metal additive manufacturing, International Journal for Numerical Methods in Engineering 119 (11) (2019) 1098–1125

  6. [6]

    Farhat, F.-X

    C. Farhat, F.-X. Roux, A method of finite element tearing and interconnecting and its parallel solution algorithm, International journal for numerical methods in engineering 32 (6) (1991) 1205–1227

  7. [7]

    Cornejo, V

    A. Cornejo, V. Mataix, F. Zárate, E. Oñate, Combination of an adaptive remeshing technique with a coupled fem–dem approach for analysis of crack propagation problems, Computational particle mechanics 7 (4) (2020) 735–752

  8. [8]

    Ambati, T

    M. Ambati, T. Gerasimov, L. De Lorenzis, A review on phase-field models of brittle fracture and a new fast hybrid formulation, Computational Mechanics 55 (2015) 383–405

  9. [9]

    Marigo, Revisiting brittle fracture as an energy minimization problem, Journal of the Mechanics and Physics of Solids 46 (8) (1998) 1319–1342

    J.-J. Marigo, Revisiting brittle fracture as an energy minimization problem, Journal of the Mechanics and Physics of Solids 46 (8) (1998) 1319–1342

  10. [10]

    Zhou, C.-C

    S.-W. Zhou, C.-C. Xia, Propagation and coalescence of quasi-static cracks in brazilian disks: an insight from a phase field model, Acta Geotechnica 14 (4) (2019) 1195–1214

  11. [11]

    Spatschek, E

    R. Spatschek, E. Brener, A. Karma, Phase field modeling of crack propagation, Philosophical Magazine 91 (1) (2011) 75–95

  12. [12]

    Svolos, H

    L. Svolos, H. M. Mourad, G. Manzini, K. Garikipati, A fourth-order phase-field fracture model: Formulation and numerical solution using a continuous/discontinuous galerkin method, Journal of the Mechanics and Physics of Solids 165 (2022) 104910

  13. [13]

    Storvik, J

    E. Storvik, J. W. Both, J. M. Sargado, J. M. Nordbotten, F. A. Radu, An accelerated staggered scheme for variational phase-field models of brittle fracture, Computer Methods in Applied Mechanics and Engineering 381 (2021) 113822. 23

  14. [14]

    S. Yang, Y. Shen, An acceleration scheme for the phase field fatigue fracture simulation with a concurrent temporal homogenization method, Computer Methods in Applied Mechanics and Engineering 416 (2023) 116294

  15. [15]

    Zhang, H

    W. Zhang, H. A. Alkhazaleh, M. Samavatian, V. Samavatian, Machine learning-assisted investigation of anisotropic elasticity in metallic alloys, Materials Today Communications 40 (2024) 109950

  16. [16]

    N. N. Vlassis, R. Ma, W. Sun, Geometric deep learning for computational mechanics part i: Anisotropic hyperelasticity, Computer Methods in Applied Mechanics and Engineering 371 (2020) 113299

  17. [17]

    Huang, J

    D. Huang, J. N. Fuhg, C. Weißenfels, P. Wriggers, A machine learning based plasticity model using proper orthogonal decomposition, Computer Methods in Applied Mechanics and Engineering 365 (2020) 113008

  18. [18]

    Ghaffari Motlagh, P

    Y. Ghaffari Motlagh, P. K. Jimack, R. de Borst, Deep learning phase-field model for brittle fractures, International Journal for Numerical Methods in Engineering 124 (3) (2023) 620–638

  19. [19]

    Kiyani, M

    E. Kiyani, M. Manav, N. Kadivar, L. De Lorenzis, G. E. Karniadakis, Predicting crack nucleation and propagation in brittle materials using deep operator networks with diverse trunk architectures, Computer Methods in Applied Mechanics and Engineering 441 (2025) 117984

  20. [20]

    Goswami, C

    S. Goswami, C. Anitescu, S. Chakraborty, T. Rabczuk, Transfer learning enhanced physics informed neural network for phase-field modeling of fracture, Theoretical and Applied Fracture Mechanics 106 (2020) 102447

  21. [21]

    Goswami, M

    S. Goswami, M. Yin, Y. Yu, G. E. Karniadakis, A physics-informed variational deeponet for predicting crack path in quasi-brittle materials, Computer Methods in Applied Mechanics and Engineering 391 (2022) 114587

  22. [22]

    Manav, R

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

  23. [23]

    Zheng, T

    B. Zheng, T. Li, H. Qi, L. Gao, X. Liu, L. Yuan, Physics-informed machine learning model for computational fracture of quasi-brittle materials without labelled data, International Journal of Mechanical Sciences 223 (2022) 107282

  24. [24]

    Dammaß, K

    F. Dammaß, K. A. Kalina, M. Kästner, Neural networks meet phase-field: A hybrid fracture model, Computer Methods in Applied Mechanics and Engineering 440 (2025) 117937

  25. [25]

    Aldakheel, E

    F. Aldakheel, E. S. Elsayed, Y. Heider, O. Weeger, Physics-based machine learning for computational fracture mechanics, Machine Learning for Computational Science and Engineering 1 (1) (2025) 18

  26. [26]

    F. M. Amin, D. W. Abueidda, P. Pantidis, M. E. Mobasher, I-fenn with deeponets: accelerating simulations in coupled multiphysics problems, Computer Methods in Applied Mechanics and Engineering 451 (2026) 118645. 24

  27. [27]

    D. W. Abueidda, M. E. Mobasher, Variational temporal convolutional networks for i-fenn thermoelasticity, Computer Methods in Applied Mechanics and Engineering 429 (2024) 117122

  28. [28]

    D. W. Abueidda, M. E. Mobasher, I-fenn for thermoelasticity based on physics-informed temporal convolutional network (pi-tcn), Computational Mechanics 74 (6) (2024) 1229–1259

  29. [29]

    Pantidis, M

    P. Pantidis, M. E. Mobasher, Integrated finite element neural network (i-fenn) for non-local continuum damage mechanics, Computer Methods in Applied Mechanics and Engineering 404 (2023) 115766

  30. [30]

    Pantidis, H

    P. Pantidis, H. Eldababy, C. M. Tagle, M. E. Mobasher, Error convergence and engineering-guided hyperparameter search of pinns: Towards optimized i-fenn performance, Computer Methods in Applied Mechanics and Engineering 414 (2023) 116160

  31. [31]

    Pantidis, H

    P. Pantidis, H. Eldababy, D. Abueidda, M. E. Mobasher, I-FENN with temporal convolutional networks: Expediting the load-history analysis of non-local gradient damage propagation, Computer Methods in Applied Mechanics and Engineering 425 (2024) 116940

  32. [32]

    Pantidis, L

    P. Pantidis, L. Svolos, D. Abueidda, M. E. Mobasher, Integrated finite element neural network (ifenn) for phase-field fracture with minimal input and generalized geometry-load handling, Computer Methods in Applied Mechanics and Engineering 448 (2026) 118485

  33. [33]

    Bourdin, G

    B. Bourdin, G. A. Francfort, J.-J. Marigo, Numerical experiments in revisited brittle fracture, Journal of the Mechanics and Physics of Solids 48 (4) (2000) 797–826

  34. [34]

    Ambrosio, V

    L. Ambrosio, V. M. Tortorelli, Approximation of functional depending on jumps by elliptic functional via t-convergence, Communications on Pure and Applied Mathematics 43 (8) (1990) 999–1036

  35. [35]

    L. Lu, P. Jin, G. Pang, Z. Zhang, G. E. Karniadakis, Learning nonlinear operators via deeponet based on the universal approximation theorem of operators, Nature machine intelligence 3 (3) (2021) 218–229

  36. [36]

    J. He, S. Kushwaha, J. Park, S. Koric, D. Abueidda, I. Jasiuk, Sequential deep operator networks (s- deeponet) for predicting full-field solutions under time-dependent loads, Engineering Applications of Artificial Intelligence 127 (2024) 107258

  37. [37]

    Huang, Y

    P. Huang, Y. Leng, C. Lian, H. Liu, Porous-deeponet: Learning the solution operators of parametric reactive transport equations in porous media, Engineering 39 (2024) 94–103

  38. [38]

    D. W. Abueidda, P. Pantidis, M. E. Mobasher, Deepokan: Deep operator network based on kolmogorov arnold networks for mechanics problems, Computer Methods in Applied Mechanics and Engineering 436 (2025) 117699

  39. [39]

    Z. Liu, Y. Wang, S. Vaidya, F. Ruehle, J. Halverson, M. Soljacic, T. Hou, M. Tegmark, Kan: Kolmogorov– arnold networks, in: International conference on learning representations, Vol. 2025, 2025, pp. 70367– 70413. 25

  40. [40]

    Somvanshi, S

    S. Somvanshi, S. A. Javed, M. M. Islam, D. Pandit, S. Das, A survey on kolmogorov-arnold network, ACM Computing Surveys 58 (2) (2025) 1–35

  41. [41]

    R. H. Peerlings, R. de Borst, W. Brekelmans, M. G. Geers, Gradient-enhanced damage modelling of concrete fracture, Mechanics of Cohesive-frictional Materials: An International Journal on Experiments, Modelling and Computation of Materials and Structures 3 (4) (1998) 323–342

  42. [42]

    Z. Lai, L. Zhao, Q. Shao, Locally enhanced neural networks for discontinuities in solid mechanics, International Journal of Mechanical Sciences (2026) 111660

  43. [43]

    Arndt, W

    D. Arndt, W. Bangerth, M. Bergbauer, M. Feder, M. Fehling, J. Heinz, T. Heister, L. Heltai, M. Kronbichler, M. Maier, et al., The deal. ii library, version 9.5, Journal of Numerical Mathematics 31 (3) (2023) 231–246

  44. [44]

    Balay, S

    S. Balay, S. Abhyankar, M. Adams, J. Brown, P. Brune, K. Buschelman, L. Dalcin, A. Dener, V. Eijkhout, W. Gropp, et al., Petsc users manual (2019)

  45. [45]

    Miehe, M

    C. Miehe, M. Hofacker, F. 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) (2010) 2765–2778. 26