Variational damage model: A new paradigm for fractures

Abstract

The computational modeling of fracture in solids using damage mechanics faces challenges with complex crack topology. This can be addressed by using a variational framework to reformulate the damage mechanics. In this paper, we propose several mathematically elegant variational damage models (VDMs) for fracture mechanics without explicitly using damage variables. Based on the energy density ϕ, the fracture energy density is formulated as \tilde{\varphi } = ϕ/(1 + ℓϕ/G_c) and the damage variable is expressed as s = ϕ/(ϕ + G_c/ℓ), which satisfy \tilde{\varphi }|_ϕ→∞ = G_c/ℓ and s|_ϕ→∞ = 1 as ϕ approaches infinity. These limits demonstrate that the new energy density converges to the Griffith energy release rate at full-damage state. The VDM profoundly modified the energy functional, implicitly incorporating the damage field. As a generalization of previous model, we propose a family of VDMs of varying orders. Additionally, we develop a multi-damage model to account for different types of energy densities, such as elastic thin plate and gradient elasticity. Using this functional, it is straightforward to deduce the governing equation for automatically evolving fractures. These formulations can be employed in conventional finite element method or other numerical methods with minimal modifications. Compared to the phase field method with the same mesh density, a sharper crack interface can be achieved. We demonstrate the capabilities of the proposed variational damage formulations using representative numerical examples.

Graphical abstract

Keywords

fracture / damage mechanics / multi-damage / variational principle

References

[1] Belytschko T, Möes N, Usui S, Parimi C. Arbitrary discontinuities in finite elements. International Journal for Numerical Methods in Engineering, 2001, 50(4): 993–1013

[2] Belytschko T, Chen H, Xu J, Zi G. Dynamic crack propagation based on loss of hyperbolicity and a new discontinuous enrichment. International Journal for Numerical Methods in Engineering, 2003, 58(12): 1873–1905

[3] Shen Y X, Lew A. An optimally convergent discontinuous Galerkin-based extended finite element method for fracture mechanics. International Journal for Numerical Methods in Engineering, 2010, 82(6): 716–755

[4] Liu W K, Jun S, Zhang Y F. Reproducing kernel particle methods. International Journal for Numerical Methods in Fluids, 1995, 20(8–9): 1081–1106

[5] Belytschko T, Lu Y Y, Gu L. Element-free Galerkin methods. International Journal for Numerical Methods in Engineering, 1994, 37(2): 229–256

[6] Rabczuk T, Areias P M A, Belytschko T. A simplified mesh-free method for shear bands with cohesive surfaces. International Journal for Numerical Methods in Engineering, 2007, 69(5): 993–1021

[7] Rabczuk T, Belytschko T. A three dimensional large deformation meshfree method for arbitrary evolving cracks. Computer Methods in Applied Mechanics and Engineering, 2007, 196(29–30): 2777–2799

[8] Zhang Y M, Zhuang X Y. Cracking elements method for dynamic brittle fracture. Theoretical and Applied Fracture Mechanics, 2019, 102: 1–9

[9] Areias P, Msekh M A, Rabczuk T. Damage and fracture algorithm using the screened Poisson equation and local remeshing. Engineering Fracture Mechanics, 2016, 158: 116–143

[10] Zhang Y M, Lackner R, Zeiml M, Mang H A. Strong discontinuity embedded approach with standard SOS formulation: Element formulation, energy-based crack-tracking strategy, and validations. Computer Methods in Applied Mechanics and Engineering, 2015, 287: 335–366

[11] Ruiz G, Pandolfi A, Ortiz M. Three-dimensional cohesive modeling of dynamic mixed-mode fracture. International Journal for Numerical Methods in Engineering, 2001, 52(1–2): 97–120

[12] Radovitzky R, Seagraves A, Tupek M, Noels L. A scalable 3D fracture and fragmentation algorithm based on a hybrid, discontinuous galerkin, cohesive element method. Computer Methods in Applied Mechanics and Engineering, 2011, 200(1–4): 326–344

[13] Miehe C, Welschinger F, Hofacker M. Thermodynamically consistent phase-field models of fracture: Variational principles and multi-field FE implementations. International Journal for Numerical Methods in Engineering, 2010, 83(10): 1273–1311

[14] Wu J Y. A unified phase-field theory for the mechanics of damage and quasi-brittle failure. Journal of the Mechanics and Physics of Solids, 2017, 103: 72–99

[15] Zhou S W, Zhuang X Y, Zhu H H, Rabczuk T. Phase field modelling of crack propagation, branching and coalescence in rocks. Theoretical and Applied Fracture Mechanics, 2018, 96: 174–192

[16] Bie Y H, Ren H L, Yan H H, Chen J Y. The unified nonlocal peridynamics-based phase-field damage theory. Theoretical and Applied Fracture Mechanics, 2023, 126: 103980

[17] Silling S A, Epton M, Weckner O, Xu J, Askari E. Peridynamic states and constitutive modeling. Journal of Elasticity, 2007, 88(2): 151–184

[18] Ren H L, Zhuang X Y, Cai Y C, Rabczuk T. Dual-horizon peridynamics. International Journal for Numerical Methods in Engineering, 2016, 108(12): 1451–1476

[19] Ren H L, Zhuang X Y, Rabczuk T. A nonlocal operator method for solving partial differential equations. Computer Methods in Applied Mechanics and Engineering, 2020, 358: 112621

[20] Ren H L, Zhuang X Y, Rabczuk T. A higher order nonlocal operator method for solving partial differential equations. Computer Methods in Applied Mechanics and Engineering, 2020, 367: 113132

[21] Francfort G A, Marigo J J. Revisiting brittle fracture as an energy minimization problem. Journal of the Mechanics and Physics of Solids, 1998, 46(8): 1319–1342

[22] Karma A, Kessler D A, Levine H. Phase-field model of mode III dynamic fracture. Physical Review Letters, 2001, 87(4): 045501

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

[24] Ulmer H, Hofacker M, Miehe C. Phase field modeling of brittle and ductile fracture. Proceedings in Applied Mathematics and Mechanics, 2013, 13(1): 533–536

[25] Wheeler M F, Wick T, Wollner W. An augmented-Lagrangian method for the phase-field approach for pressurized fractures. Computer Methods in Applied Mechanics and Engineering, 2014, 271: 69–85

[26] Aldakheel F, Hudobivnik B, Hussein A, Wriggers P. Phase-field modeling of brittle fracture using an efficient virtual element scheme. Computer Methods in Applied Mechanics and Engineering, 2018, 341: 443–466

[27] de Borst R, Verhoosel C V. Gradient damage vs phase-field approaches for fracture: Similarities and differences. Computer Methods in Applied Mechanics and Engineering, 2016, 312: 78–94

[28] SimoJ CHughes T J R. Computational Inelasticity. New York: Springer New York, 2006

[29] Hackl K, Fischer F D. On the relation between the principle of maximum dissipation and inelastic evolution given by dissipation potentials. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2008, 464(2089): 117–132

[30] KachanovL. Introduction to Continuum Damage Mechanics. Dordrecht: Springer Dordrecht, 1986

[31] OchsnerA. Continuum Damage Mechanics. Berlin: Springer, 2016

[32] ZhangW HCai Y Q. Continuum Damage Mechanics and Numerical Applications. Berlin: Springer, 2010

[33] Feng X Q, Yu S W. Damage micromechanics for constitutive relations and failure of microcracked quasi-brittle materials. International Journal of Damage Mechanics, 2010, 19(8): 911–948

[34] Ren X D, Chen J S, Li J, Slawson T R, Roth M J. Micro-cracks informed damage models for brittle solids. International Journal of Solids and Structures, 2011, 48(10): 1560–1571

[35] MurakamiS. Continuum Damage Mechanics: A Continuum Mechanics Approach to the Analysis of Damage and Fracture. Dordrecht: Springer Dordrecht, 2012

[36] Jirásek M, Zeman J. Localization study of a regularized variational damage model. International Journal of Solids and Structures, 2015, 69: 131–151

[37] Griffith A A. The phenomena of rupture and flow in solids. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 1921, 221: 163–198

[38] Borden M J, Hughes T J R, Landis C M, Anvari A, Lee I J. A phase-field formulation for fracture in ductile materials: Finite deformation balance law derivation, plastic degradation, and stress triaxiality effects. Computer Methods in Applied Mechanics and Engineering, 2016, 312: 130–166

[39] Wolfram S. The mathematica book. Assembly Automation, 1999, 19(1): 77

[40] Silling S A. Reformulation of elasticity theory for discontinuities and long-range forces. Journal of the Mechanics and Physics of Solids, 2000, 48(1): 175–209

[41] Ren H L, Zhuang X Y, Fu X L, Li Z Y, Rabczuk T. Bond-based nonlocal models by nonlocal operator method in symmetric support domain. Computer Methods in Applied Mechanics and Engineering, 2024, 418: 116230

[42] Varma S S, Mangipudi K R, Budarapu P R. A coupled quantum-molecular mechanics approach for performance analysis of defective silicon based photovoltaic solar cells. Physica Scripta, 2023, 98(3): 035007

[43] BatheK J. Finite element method. In: Wah B W, ed. Wiley Encyclopedia of Computer Science and Engineering. Hoboken, NJ: John Wiley & Sons, Inc., 2007, 1–12

[44] ReddyJ N. Introduction to the Finite Element Method. New York: McGraw-Hill Education, 2019

[45] Wu Y C, Xiao J Z. Implementation of the multiscale stochastic finite element method on elliptic PDE problems. International Journal of Computational Methods, 2017, 14(1): 1750003

[46] Hughes T J R, Cottrell J A, Bazilevs Y. Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement. Computer Methods in Applied Mechanics and Engineering, 2005, 194(39–41): 4135–4195

[47] Nguyen V P, Anitescu C, Bordas S P A, Rabczuk T. Isogeometric analysis: An overview and computer implementation aspects. Mathematics and Computers in Simulation, 2015, 117: 89–116

[48] LiuG RGu Y T. An Introduction to Meshfree Methods and Their Programming. Dordrecht: Springer Dordrecht, 2005

[49] Ren H L, Zhuang X Y, Rabczuk T, Zhu H H. Dual-support smoothed particle hydrodynamics in solid: Variational principle and implicit formulation. Engineering Analysis with Boundary Elements, 2019, 108: 15–29

[50] Bie Y H, Ding K J, Zhao Z F, Wei Y G. The adaptive coupling of dual-horizon peridynamic element and finite element for the progressive failure of materials. International Journal of Fracture, 2024, 245: 89–114

[51] Hudobivnik B, Aldakheel F, Wriggers P. A low order 3D virtual element formulation for finite elasto–plastic deformations. Computational Mechanics, 2019, 63(2): 253–269

[52] Raissi M, Perdikaris P, Karniadakis G E. Physics-informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations. Journal of Computational Physics, 2019, 378: 686–707

[53] Samaniego E, Anitescu C, Goswami S, Nguyen-Thanh V M, Guo H W, Hamdia K, Zhuang X Y, Rabczuk T. An energy approach to the solution of partial differential equations in computational mechanics via machine learning: Concepts, implementation and applications. Computer Methods in Applied Mechanics and Engineering, 2020, 362: 112790

[54] Siruvuri S V, Budarapu P R, Paggi M. Influence of cracks on fracture strength and electric power losses in silicon solar cells at high temperatures: Deep machine learning and molecular dynamics approach. Applied Physics A: Materials Science & Processing, 2023, 129(6): 408

[55] Miehe C, Lambrecht M. Algorithms for computation of stresses and elasticity moduli in terms of Seth–Hill’s family of generalized strain tensors. Communications in Numerical Methods in Engineering, 2001, 17(5): 337–353

[56] WinklerB J. Traglastuntersuchungen von unbewehrten und bewehrten Betonstrukturen auf der Grundlage eines objektiven Werkstoffgesetzes für Beton. Dissertation for the Doctoral Degree. Innsbruck: University of Innsbruck, 2001

[57] Gravouil A, Möes N, Belytschko T. Non-planar 3D crack growth by the extended finite element and level sets––Part II: Level set update. International Journal for Numerical Methods in Engineering, 2002, 53(11): 2569–2586

Acknowledgements

H. Ren and X. Zhuang gratefully acknowledge the support from ERC Grant COTOFLEXI (No. 802205). We would like to thank Ms. Cancan Su for the useful discussion during the revision of the manuscript.

Funding note

Open Access funding enabled and organized by DEAL agreement.

Open Access

This article is licensed under a Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/), which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Competing interests

The authors declare that they have no competing interests.

RIGHTS & PERMISSIONS

2025 The Author(s). This article is published with open access at link.springer.com and journal.hep.com.cn