Abstract

Phase-field fracture models provide a powerful approach to modeling fracture, potentially enabling the unguided prediction of the crack growth in complex patterns. To ensure that only tensile stresses and not compressive stresses drive crack growth, several models have been proposed that aim to distinguish between the compressive and tensile loads. However, these models have a critical shortcoming: they do not account for the crack direction, and hence they cannot distinguish between crack-normal tensile stresses that drive crack growth and crack-parallel stresses that do not. In this study, we apply a phase-field fracture model, developed in our earlier work, that uses the crack direction to distinguish crack-parallel stresses from crack-normal stresses. This provides a transparent energetic formulation that drives cracks to grow in when crack faces open or slide past each other, while the cracks respond like the intact solid when the crack faces contact under normal compressive loads. We compare our approach against two widely used approaches, Spectral splitting and the volumetric-deviatoric splitting, and find that these predict unphysical crack growth and unphysical stress concentrations under loading conditions in which these should not occur. Specifically, we show that the splitting models predict spurious crack growth and stress concentration under pure crack-parallel normal stresses. However, our formulation, which resolves the crack-parallel stresses from the crack-normal stresses, predicts these correctly.

1 Introduction

Phase-field fracture models introduce a new field, in addition to the deformation, that indicates if a material element is intact or fractured, and then regularize this phase-field using higher-order derivatives. The higher-order derivatives endow the phase-field with sufficient smoothness to enable the use of standard numerical methods to predict the crack growth without specifying crack paths beforehand. In other words, rather than sharp singular cracks as in the classical approach to fracture, cracks are now regularized to be nonsingular variations of the phase-field, making them amenable to standard numerical approaches. The elastic response is coupled to the phase-field by setting that the regions that are fractured cannot carry any load, i.e., the elastic energy density goes to zero in fractured regions.

However, the combination of crack regularization to a finite volume and setting the fractured elastic energy to zero causes unphysical behavior, such as crack growth under compressive loads. To address these issues, Miehe et al. [1] proposed a partition of the elastic energy into tensile and compressive parts (denoted the “spectral” split), and allowed the damaged region to sustain only the compressive part. An alternative class of approaches, following Amor et al. [2], proposed a partition of the energy into compressive hydrostatic, tensile hydrostatic, and deviatoric parts (denoted the “volumetric-deviatoric” or “VolDev” split), and allowed the damaged region to sustain only the compressive hydrostatic part. This has been extended to the finite deformation setting in Refs. [36]; however, none of these works consider the orientation of the crack. Other modifications of the energy decompositions include Ref. [7] that proposed the star-convex model as a modification of the volumetric-deviatoric decomposition; Ref. [8] that extended the spectral decomposition to account for mixed mode fracture; and Ref. [9] to account for anisotropy. These approaches are reviewed and compared in Refs. [1013] and elsewhere.

However, all of the spectral and VolDev splitting approaches have a key shortcoming: they do not account for the orientation of the crack. Compressive stresses along the crack face and across the crack face are both treated in exactly the same way. For instance, tension across the crack faces that drives the crack growth is not distinguished from tension along the length of the crack, that does not directly drive crack growth but it is important in other ways. For instance, Bazant et al. [14] highlighted that phase-field models do not accurately capture the influence of crack-parallel stresses, which significantly impact the fracture energy and the fracture process zone width; they argue that the inability of phase-field models to incorporate the tensorial nature of stress interactions leads to discrepancies when compared to experimental results from the gap test. Similarly, Nguyen et al. [15] emphasized that phase-field models lack the ability to explicitly account for the history-dependent effects of T-stress, which are crucial for realistic simulations of mixed-mode fracture in concrete and quasi-brittle materials.

An important departure from the idea of splitting approaches is the work by Steinke and Kaliske [16], introduced the idea that the crack orientation plays an essential role in determining the response. Based on this idea, Refs. [1619] proposed modifications of the stress tensor to account for the difference in crack-normal and crack-parallel tractions. While our approach in Refs. [20,21] builds on the ideas of Ref. [16], an important distinction is that we modify the energy rather than the stress, and hence our model is hyperelastic. In contrast, as shown in Ref. [20], the stress-based approach is not generally hyperelastic. Our approach enables us to distinguish between the crack-parallel stresses and crack-closing stresses as well as allow for crack faces sliding past each other. Our approach applies a QR (upper triangular) decomposition of the deformation gradient in the basis associated with the local crack normal, enabling us to distinguish between the crack faces sliding past each other—which is formulated to cost no energy; crack opening—which is also formulated to cost no energy; and crack-closing and crack face contact—which are formulated to have an elastic response that is identical to the intact material. This is shown schematically in Fig. 1.

Fig. 1
The top row shows different loadings, and the middle and lower rows show the idealized deformation for intact and cracked specimens, respectively. Based on this idealization, we assign zero energy to modes (a) and (b). The dotted lines in the second and third rows show the undeformed configuration. This decomposition follows Refs. [16,21].
Fig. 1
The top row shows different loadings, and the middle and lower rows show the idealized deformation for intact and cracked specimens, respectively. Based on this idealization, we assign zero energy to modes (a) and (b). The dotted lines in the second and third rows show the undeformed configuration. This decomposition follows Refs. [16,21].
Close modal

In this article, we compare the phase-field fracture model developed by Hakimzadeh et al. [20,21], who distinguished between the crack-parallel stresses and crack-closing stresses, against the spectral [1] and VolDev [2] splitting approaches. Specifically, in Sec. 3, we examine the effect of tensile and compressive crack-parallel stresses and show that the crack growth and stresses predicted by the spectral and VolDev approaches are unphysical, in contrast to the model from Ref. [20].

2 Formulation of the Model

2.1 Classical Phase-Field Fracture Model.

The variational fracture model by Francfort and Marigo [22] describes brittle fracture through energy minimization
(2.1)
where y is the displacement, Γ is the crack surface, Ω is the reference domain, W is the elastic energy density, Gc is the work to fracture, and Hn1(Γ) is the Hausdorff area of Γ. The model asserts that the deformation and the crack surface are obtained by minimizing E subject to irreversibility, i.e., the crack cannot heal.
Due to the computational challenges in working with this model directly—the solution lives in function spaces that are very challenging to approximate numerically—a regularized phase-field formulation was introduced by Ambrosio and Tortorelli [23,24], inspired by Modica [25,26]
(2.2)
where ϕ(x) is a scalar phase-field that quantifies the level of damage at every referential point x, and is sufficiently smooth to allow standard finite element method (FEM) approximation. Specifically, ϕ, bounded by 0ϕ1, indicates that the material is intact at ϕ(x)1 and completely fractured at ϕ(x)0. The small parameter ϵ controls the crack width. The method has been used very widely for fracture calculations and we do not review this extensive literature.

2.2 Accounting for Crack Face Contact.

A key limitation of the standard phase-field fracture model is its inability to distinguish between tensile, shear, and compressive loadings across or along the crack face, often leading to unphysical crack propagation under compression. Various modifications have been proposed to address this issue, but they fail to incorporate the crack orientation explicitly [1,2,610,27]. As noted in Ref. [16] and further demonstrated in Refs. [20,21] and elsewhere, correctly capturing the crack behavior under mixed-mode loading requires a formulation that accounts for the crack face normal.

To overcome this limitation, we adopt the phase-field formulation introduced in Ref. [20], which incorporates an effective crack energy density Wd that depends on the crack orientation. Using the QR (or Gram–Schmidt) decomposition of the deformation gradient y, e.g., Ref. [28], in the basis of the crack enables us to distinguish the various modes shown in Fig. 1. We then assigned zero energy to modes in which the crack faces slide past each other or open up, and assign the intact response to modes in which the crack faces contact. This ensures a physically consistent response, particularly under compressive loading, by preventing unphysical interpenetration of crack faces, as demonstrated through several examples in Refs. [20,21].

Our approach uses a free energy of the form
(2.3)

The key feature of this model is the introduction of the effective crack energy density Wd that is active in the fractured volumes. It imparts the properties of an idealized sharp crack to the regularized (finite volume) phase-field crack. The function Wd takes as argument the deformation gradient, as well as a vector field d that corresponds to the normal to the crack face which is essential to be able to decompose the deformation into the modes shown in Fig. 1. However, d is not an independent field; rather it is a function of ϕ as described in Ref. [21]. For brevity, we do not present the detailed expressions for Wd but refer instead to Refs. [20,21]. In this article, we use an isotropic linear elastic response, with Lamé moduli λ and μ.

2.3 Energy Splitting Approaches.

A simplistic approach to modeling the energy density of damaged volumes is to assume that the damaged region carries zero elastic energy. However, this assumption often leads to unphysical behaviors, such as the interpenetration of crack faces, incorrect mechanical response when the crack closes, and crack growth under compression. To address these shortcomings, alternative formulations have been proposed, notably by the spectral split by Miehe et al. [1] and the VolDev split by Amor et al. [2], which decompose the energy and selectively associate specific energy components with the damaged region. While these methods offer improvements over the zero-energy assumption, they do not account for crack orientation in the energetic decomposition, leading to unphysical predictions under complex stress states.

2.3.1 Spectral Splitting.

Following Ref. [1], the spectral splitting approach defines the compressive (ψ0) and tensile (ψ0+) energies as
(2.4)
where ε1,ε2,ε3 are the principal strains with the corresponding principal directions n1,n2,n3. Here x+:=max{0,x} and x:=min{0,x}, and only the compressive energy ψ0 is taken to contribute to the energy in the cracked region. This results in the stress
(2.5)
which inherently prevents tensile stresses in all directions.

In the Appendix of Ref. [20], we identified significant challenges associated with spectral splitting. Specifically, this approach does not accurately capture the crack response under complex loading conditions, such as uniaxial tension parallel to the crack. Despite these limitations, many studies on the fracture modeling continue to rely on this approach or variants like those proposed by Lo et al. [29], which builds on the work of Ref. [1]. Recent studies in studying the crack behavior under complex loadings, including Ref. [30], utilize this extended method by Lo et al. [29]. Other examples of fracture modeling include recent work by Clayton et al. [31], which adopts the formulation introduced by Miehe et al. [32] and Stocek et al. [33], which is built on the frameworks of Refs. [34] and [1].

2.3.2 Volumetric-Deviatoric Splitting.

The VolDev split proposed by Amor et al. [2] aims to allow cracks to resist compressive hydrostatic stresses while excluding tensile hydrostatic and deviatoric stresses. The energy is expressed as
(2.6)
where κ=λ+2μ/3 is the bulk modulus, and εD:=ε13(trε)I is the deviatoric strain component. Here only the compressive hydrostatic term 12κtrε2 is retained for the cracked region’s energy, giving the stress response
producing a purely compressive hydrostatic response.

As discussed in the Appendix of our previous work [20], the regularized approach proposed by Amor et al. [2] encounters limitations in accurately capturing the crack response under complex stress states. These include uniaxial tension parallel to the crack plane and uniaxial compression normal to the crack. Despite these limitations, this approach remains widely used in several studies that model crack behavior in complex loadings.

2.4 Numerical Implementation.

We employ the FEM for numerical simulations, using the open-source FEniCS library, which has been extensively utilized in mechanics applications [35,36]. The details are provided in Ref. [37] for the spectral split approach, Ref. [38] for the VolDev split approach, and Ref. [20] for our approach.

3 Crack Growth and Stress Distribution Under Crack-Parallel Stresses

Here we analyze the crack growth and the stress distribution when the crack is subject to crack-parallel (tensile and compressive) stresses using spectral splitting, VolDev splitting, and our approach.

Consider a specimen containing a crack with normal e2, subjected to a far-field stress of σ0e1e1, with σ0>0. Assuming isotropic elasticity, this corresponds equivalently to the following far-field strain:

(3.1)
The resulting stress response within the crack using the spectral splitting model is
(3.2)

We highlight two key unphysical aspects of this result. First, when the Poisson ratio ν is positive, we find an unphysical compressive stress in the e2 direction. Second, the stress along the e1 direction is zero, whereas under a far-field stress parallel to the crack, the material is expected to behave as an intact material, sustaining the far-field stress as it is. This incorrect stress distribution leads to spurious crack growth. The artificial lack of stiffness in tension parallel to the crack creates spurious stress concentrations near the crack tip, as the intact material ahead of the tip sustains the far-field stress in the e1 direction, while the cracked region does not. Moreover, this discrepancy can significantly affect the calculation of the T-stress [14].

In the VolDev model, we obtain σsplit=0. However, the material should fully sustain the far-field tensile stress parallel to the crack. This unphysical response leads to spurious stress concentrations at the crack tip, which influence both the crack driving force and the calculation of the T-stress.

3.1 Numerical Evaluation of the Crack Growth and Stress Distributions.

We next numerically evaluate the crack path and stress distribution in a finite domain. For a clear comparison, we use the same domain, mesh, and material properties across all three approaches, and first set ν=0 to eliminate crack faces separating or contacting due to the crack-parallel loading. We then examine having ν>0 and ν<0.

In our calculations, we set the vertical displacement and the horizontal tractions to zero on the top and bottom faces; the horizontal displacement is controlled on the left and right faces; and the vertical tractions are set to zero on the left and right faces.

3.1.1 Crack-Parallel Tension.

We first examine the effect of tensile loading in Fig. 2, which shows the initial crack, the subsequent crack growth, and the stress distribution at the first load step, for all three models. For simplicity of interpretation, we set ν=0. We observe that our approach predicts essentially a uniform stress distribution, that corresponds to the cracked specimen having the same response as an intact specimen, and there is no crack growth. Both of these features correspond to physical expectations.

Fig. 2
Comparison of crack growth paths and stress fields under tensile crack-parallel stress: (a) initial crack, shown by ϕ, (b) no crack growth observed in our model, (c) uniform stress distribution in our model, (d) unphysical crack growth with spectral splitting, (e) spurious stress concentrations with spectral splitting, (f) unphysical crack growth with VolDev splitting, and (g) spurious stress concentrations with VolDev splitting
Fig. 2
Comparison of crack growth paths and stress fields under tensile crack-parallel stress: (a) initial crack, shown by ϕ, (b) no crack growth observed in our model, (c) uniform stress distribution in our model, (d) unphysical crack growth with spectral splitting, (e) spurious stress concentrations with spectral splitting, (f) unphysical crack growth with VolDev splitting, and (g) spurious stress concentrations with VolDev splitting
Close modal

In both the spectral and VolDev splits, the stresses are significantly non-uniform with spurious stress concentrations at the crack tips that were anticipated by the simple analysis above. The spurious stress distribution, in turn, drives unphysical crack growth.

3.1.2 Crack-Parallel Compression.

We next examine the effect of compressive loading in Fig. 3, which shows the initial crack, the subsequent crack growth, and the stress distribution at the first load step, for all three models; we set ν=0. We observe again that our approach predicts essentially a uniform stress distribution, that corresponds to the cracked specimen having the same response as an intact specimen, and there is no crack growth. Both of these features correspond to physical expectations.

Fig. 3
Comparison of crack growth paths and stress fields under compressive crack-parallel stress: (a) initial crack, shown by ϕ, (b) no crack growth observed in our model, (c) uniform stress distribution in our model, (d) no crack growth with spectral splitting, (e) spurious stress concentrations with spectral splitting, (f) unphysical crack growth with VolDev splitting, and (g) spurious stress concentrations with VolDev splitting
Fig. 3
Comparison of crack growth paths and stress fields under compressive crack-parallel stress: (a) initial crack, shown by ϕ, (b) no crack growth observed in our model, (c) uniform stress distribution in our model, (d) no crack growth with spectral splitting, (e) spurious stress concentrations with spectral splitting, (f) unphysical crack growth with VolDev splitting, and (g) spurious stress concentrations with VolDev splitting
Close modal

Further, in both the spectral and VolDev splits, the stresses are significantly non-uniform with spurious stress concentrations at the crack tips that were anticipated by the simple analysis above. The spurious stress distribution, in turn, drives unphysical crack growth that should simply not occur given that the specimen is under compression.

3.1.3 Effect of Poisson’s Ratio.

In the calculations above, we have set ν=0 to aid the interpretation of the stress fields. Here we compare the models with crack-parallel tension and compression, with positive and negative ν. We examine a horizontal crack in a large rectangular specimen, with the top and bottom traction-free, and the displacement fixed—to be zero in the vertical direction and nonzero in the horizontal direction—on the sides.

Figure 4 compares the stresses and deformations under crack-parallel compression with ν>0. In our model, the stress field is uniform and the deformation corresponds to uniform barreling due to the lateral expansion, as expected physically. However, in both the spectral splitting and VolDev models, the stress is clearly inhomogeneous when it should not be, and the deformation around the crack shows unphysical features.

Fig. 4
Comparison of stress fields and deformation without crack growth under compressive crack-parallel stress with ν>0. Deformations are magnified 150× for visualization: (a) stress field in our model, (b) deformation with our model, (c) stress field with spectral splitting, (d) deformation with spectral splitting, (e) stress field with VolDev splitting, and (f) deformation with VolDev splitting
Fig. 4
Comparison of stress fields and deformation without crack growth under compressive crack-parallel stress with ν>0. Deformations are magnified 150× for visualization: (a) stress field in our model, (b) deformation with our model, (c) stress field with spectral splitting, (d) deformation with spectral splitting, (e) stress field with VolDev splitting, and (f) deformation with VolDev splitting
Close modal

Figure 5 compares the models under crack-parallel compression with ν<0. In our model, as expected, there is lateral contraction that causes the crack faces to lose contact and move apart; the stress field is therefore inhomogeneous with stress concentrations at the crack tips, as is typical. However, in both the spectral splitting and VolDev models, while they appear to capture the lateral contraction, the stress distribution is not at all as expected around crack tips, and the deformation shows unphysical features near the crack tips.

Fig. 5
Comparison of stress fields and deformation without crack growth under compressive crack-parallel stress with ν<0. Deformations are magnified 150× for visualization: (a) stress field in our model, (b) deformation with our model, (c) stress field with spectral splitting, (d) deformation with spectral splitting, (e) stress field with VolDev splitting, and (f) deformation with VolDev splitting
Fig. 5
Comparison of stress fields and deformation without crack growth under compressive crack-parallel stress with ν<0. Deformations are magnified 150× for visualization: (a) stress field in our model, (b) deformation with our model, (c) stress field with spectral splitting, (d) deformation with spectral splitting, (e) stress field with VolDev splitting, and (f) deformation with VolDev splitting
Close modal

Figure 6 compares the models under crack-parallel tension with ν>0. In our model, as expected, there is lateral contraction that causes the crack faces to lose contact and move apart, and the stress field is therefore inhomogeneous with stress concentrations at the crack tips. However, in both the spectral splitting and VolDev models, while they capture the lateral contraction, the stress distribution is not at all as expected around crack tips, and the deformation shows unphysical features near the crack tips.

Fig. 6
Comparison of stress fields and deformation without crack growth under tensile crack-parallel stress with ν>0. Deformations are magnified 150× for visualization: (a) stress field in our model, (b) deformation with our model, (c) stress field with spectral splitting, (d) deformation with spectral splitting, (e) stress field with VolDev splitting, and (f) deformation with VolDev splitting
Fig. 6
Comparison of stress fields and deformation without crack growth under tensile crack-parallel stress with ν>0. Deformations are magnified 150× for visualization: (a) stress field in our model, (b) deformation with our model, (c) stress field with spectral splitting, (d) deformation with spectral splitting, (e) stress field with VolDev splitting, and (f) deformation with VolDev splitting
Close modal

Finally, Fig. 7 compares the models under crack-parallel tension with ν<0. In our model, as expected, there is lateral expansion and the overall deformation is uniform as the crack faces remain in contact; the stress field is similarly homogeneous around the crack. However, in both the spectral splitting and VolDev models, the deformation is completely unphysical: the material within the damaged region shrinks laterally, corresponding to the crack faces inter-penetrating. Further, the stress distribution is not at all as expected around crack tips.

Fig. 7
Comparison of stress fields and deformation without crack growth under tensile crack-parallel stress with ν<0. Deformations are magnified 150× for visualization: (a) stress field in our model, (b) deformation with our model, (c) stress field with spectral splitting, (d) deformation with spectral splitting, (e) stress field with VolDev splitting, and (f) deformation with VolDev splitting
Fig. 7
Comparison of stress fields and deformation without crack growth under tensile crack-parallel stress with ν<0. Deformations are magnified 150× for visualization: (a) stress field in our model, (b) deformation with our model, (c) stress field with spectral splitting, (d) deformation with spectral splitting, (e) stress field with VolDev splitting, and (f) deformation with VolDev splitting
Close modal

4 Concluding Remarks

We have shown that the widely-used spectral and VolDev splitting approaches lead to unphysical crack growth and spurious stress concentrations in the important settings of crack-parallel stresses. We have further shown that accounting for the crack face contact through the crack-response model proposed in Ref. [20] provides the correct predictions in these settings, suggesting the use of this model in situations where crack-parallel stresses play an important role [14].

We have idealized the crack as having no friction between the crack faces for simplicity. For future studies aiming to infuse more realistic physics into the model, it is important to further incorporate friction in the model [17,18]. Another important area of application is to couple the model to poromechanics, e.g., by adapting Refs. [28,30,3942]. A further important area is to combine the crack-response model with crack nucleation models, following Refs. [4347]. For all of these future directions, it is essential to predict the correct stress fields that in turn provide the correct crack growth driving force.

Our focus in this article has been on the small-strain regime for simplicity. However, the crack-response model of Ref. [20] is valid in the large deformation setting and can be readily incorporated in such models, e.g., Refs. [13,4851].

Acknowledgment

We thank Meenu Krishnan for useful discussions. We acknowledge financial support from NSF (2108784, 2012259, and 2342349), ARO (MURI W911NF-24-2-0184), AEI (PID2021-124195NB-C32, PCI2024-155023-2, CEX2023-001347-S funded by MICIU/AEI/10.13039/501100011033), and ERC (834728); and NSF for XSEDE computing resources provided by Pittsburgh Supercomputing Center.

Conflict of Interest

There are no conflicts of interest.

Data Availability Statement

The code implementing all three approaches for the calculations in this article is available online.2

Footnote

References

1.
Miehe
,
C.
,
Welschinger
,
F.
, and
Hofacker
,
M.
,
2010
, “
Thermodynamically Consistent Phase-Field Models of Fracture: Variational Principles and Multi-field FE Implementations
,”
Int. J. Numer. Methods Eng.
,
83
(
10
), pp.
1273
1311
.
2.
Amor
,
H.
,
Marigo
,
J.-J.
, and
Maurini
,
C.
,
2009
, “
Regularized Formulation of the Variational Brittle Fracture With Unilateral Contact: Numerical Experiments
,”
J. Mech. Phys. Solids
,
57
(
8
), pp.
1209
1229
.
3.
Tian
,
F.
,
Zeng
,
J.
,
Zhang
,
M.
, and
Li
,
L.
,
2022
, “
Mixed Displacement–Pressure-Phase Field Framework for Finite Strain Fracture of Nearly Incompressible Hyperelastic Materials
,”
Comput. Methods Appl. Mech. Eng.
,
394
, p.
114933
.
4.
Xing
,
C.
,
Yu
,
T.
,
Sun
,
Y.
, and
Wang
,
Y.
,
2023
, “
An Adaptive Phase-Field Model With Variable-Node Elements for Fracture of Hyperelastic Materials at Large Deformations
,”
Eng. Fract. Mech.
,
281
, p.
109115
.
5.
Najmeddine
,
A.
, and
Shakiba
,
M.
,
2024
, “
Efficient BFGS Quasi-Newton Method for Large Deformation Phase-Field Modeling of Fracture in Hyperelastic Materials
,”
Eng. Fract. Mech.
,
310
, p.
110463
.
6.
Tang
,
S.
,
Zhang
,
G.
,
Guo
,
T. F.
,
Guo
,
X.
, and
Liu
,
W. K.
,
2019
, “
Phase Field Modeling of Fracture in Nonlinearly Elastic Solids Via Energy Decomposition
,”
Comput. Methods Appl. Mech. Eng.
,
347
, pp.
477
494
.
7.
Vicentini
,
F.
,
Zolesi
,
C.
,
Carrara
,
P.
,
Maurini
,
C.
, and
de Lorenzis
,
L.
,
2024
, “
On the Energy Decomposition in Variational Phase-Field Models for Brittle Fracture Under Multi-axial Stress States
,”
Int. J. Fract.
,
247
, p.
291
.
8.
Wang
,
F.
,
Shao
,
J.
, and
Huang
,
H.
,
2021
, “
A Phase-Field Modeling Method for the Mixed-Mode Fracture of Brittle Materials Based on Spectral Decomposition
,”
Eng. Fract. Mech.
,
242
, p.
107473
.
9.
van Dijk
,
N. P.
,
Espadas-Escalante
,
J. J.
, and
Isaksson
,
P.
,
2020
, “
Strain Energy Density Decompositions in Phase-Field Fracture Theories for Orthotropy and Anisotropy
,”
Int. J. Solids Struct.
,
196
, pp.
140
153
.
10.
Zhang
,
S.
,
Jiang
,
W.
, and
Tonks
,
M. R.
,
2022
, “
Assessment of Four Strain Energy Decomposition Methods for Phase Field Fracture Models Using Quasi-Static and Dynamic Benchmark Cases
,”
Mater. Theory
,
6
(
1
), p.
6
.
11.
Gupta
,
A.
,
Nguyen
,
D. T.
,
Hirshikesh
, and
Duddu
,
R.
,
2024
, “
Damage Mechanics Challenge: Predictions From an Adaptive Finite Element Implementation of the Stress-Based Phase-Field Fracture Model
,”
Eng. Fract. Mech.
,
306
, p.
110252
.
12.
Rahaman
,
M. M.
,
2022
, “
An Open-Source Implementation of a Phase-Field Model for Brittle Fracture Using Gridap in Julia
,”
Math. Mech. Solids
,
27
(
11
), pp.
2404
2427
.
13.
Clayton
,
J. D.
,
2021
, “
Nonlinear Thermodynamic Phase Field Theory With Application to Fracture and Dynamic Inelastic Phenomena in Ceramic Polycrystals
,”
J. Mech. Phys. Solids
,
157
, p.
104633
.
14.
Bažant
,
Z. P.
,
Nguyen
,
H. T.
, and
Dönmez
,
A. A.
,
2022
, “
Critical Comparison of Phase-Field, Peridynamics, and Crack Band Model M7 in Light of Gap Test and Classical Fracture Tests
,”
ASME J. Appl. Mech.
,
89
(
6
), p.
061008
.
15.
Nguyen
,
H. T.
,
Pathirage
,
M.
,
Cusatis
,
G.
, and
Bažant
,
Z. P.
,
2020
, “
Gap Test of Crack-Parallel Stress Effect on Quasibrittle Fracture and Its Consequences
,”
ASME J. Appl. Mech.
,
87
(
7
), p.
071012
.
16.
Steinke
,
C.
, and
Kaliske
,
M.
,
2019
, “
A Phase-Field Crack Model Based on Directional Stress Decomposition
,”
Comput. Mech.
,
63
(
5
), pp.
1019
1046
.
17.
Fei
,
F.
, and
Choo
,
J.
,
2020
, “
A Phase-Field Model of Frictional Shear Fracture in Geologic Materials
,”
Comput. Methods Appl. Mech. Eng.
,
369
, p.
113265
.
18.
Fei
,
F.
, and
Choo
,
J.
,
2020
, “
A Phase-Field Method for Modeling Cracks With Frictional Contact
,”
Int. J. Numer. Methods Eng.
,
121
(
4
), pp.
740
762
.
19.
Agrawal
,
V.
,
2016
, “
Multiscale Phase-Field Model for Phase Transformation and Fracture
,” Ph.D. thesis,
Carnegie Mellon University
,
Pittsburgh, PA
.
20.
Hakimzadeh
,
M.
,
Agrawal
,
V.
,
Dayal
,
K.
, and
Mora-Corral
,
C.
,
2022
, “
Phase-Field Finite Deformation Fracture With an Effective Energy for Regularized Crack Face Contact
,”
J. Mech. Phys. Solids
,
167
, p.
104994
.
21.
Hakimzadeh
,
M.
,
Mora-Corral
,
C.
,
Walkington
,
N.
,
Buscarnera
,
G.
, and
Dayal
,
K.
,
2025
, “
Phase-Field Modeling of Fracture Under Compression and Confinement in Anisotropic Geomaterials
,”
Int. J. Numer. Anal. Methods Geomech.
,
49
(
4
), pp.
1319
1335
.
22.
Francfort
,
G. A.
, and
Marigo
,
J.-J.
,
1998
, “
Revisiting Brittle Fracture as an Energy Minimization Problem
,”
J. Mech. Phys. Solids
,
46
(
8
), pp.
1319
1342
.
23.
Ambrosio
,
L.
, and
Tortorelli
,
V. M.
,
1990
, “
Approximation of Functionals Depending on Jumps by Elliptic Functionals Via Γ-Convergence
,”
Comm. Pure Appl. Math.
,
43
(
8
), pp.
999
1036
.
24.
Ambrosio
,
L.
, and
Tortorelli
,
V. M.
,
1992
, “
On the Approximation of Free Discontinuity Problems
,”
Boll. Un. Mat. Ital. B (7)
,
6
(
1
), pp.
105
123
.
25.
Modica
,
L.
, and
Mortola
,
S.
,
1977
, “
Un Esempio Di Γ-Convergenza
,”
Boll. Un. Mat. Ital. B (5)
,
14
(
1
), pp.
285
299
.
26.
Modica
,
L.
,
1987
, “
The Gradient Theory of Phase Transitions and the Minimal Interface Criterion
,”
Arch. Rational Mech. Anal.
,
98
(
2
), pp.
123
142
.
27.
De Lorenzis
,
L.
, and
Maurini
,
C.
,
2022
, “
Nucleation Under Multi-axial Loading in Variational Phase-Field Models of Brittle Fracture
,”
Int. J. Fract.
,
237
(
1-2
), pp.
61
81
.
28.
Clayton
,
J. D.
, and
Freed
,
A. D.
,
2020
, “
A Constitutive Framework for Finite Viscoelasticity and Damage Based on the Gram–Schmidt Decomposition
,”
Acta Mechanica
,
231
(
8
), pp.
3319
3362
.
29.
Lo
,
Y.-S.
,
Borden
,
M. J.
,
Ravi-Chandar
,
K.
, and
Landis
,
C. M.
,
2019
, “
A Phase-Field Model for Fatigue Crack Growth
,”
J. Mech. Phys. Solids
,
132
, p.
103684
.
30.
Sun
,
X.
, and
Duddu
,
R.
,
2021
, “
A Poro-Damage Phase Field Model for Hydrofracturing of Glacier Crevasses
,”
Extreme Mech. Lett.
,
45
, p.
101277
.
31.
Clayton
,
T.
,
Duddu
,
R.
,
Siegert
,
M.
, and
Martinez-Paneda
,
E.
,
2022
, “
A Stress-Based Poro-Damage Phase Field Model for Hydrofracturing of Creeping Glaciers and Ice Shelves
,”
Eng. Fract. Mech.
,
272
, p.
108693
.
32.
Miehe
,
C.
,
Schaenzel
,
L.-M.
, and
Ulmer
,
H.
,
2015
, “
Phase Field Modeling of Fracture in Multi-physics Problems. Part I. Balance of Crack Surface and Failure Criteria for Brittle Crack Propagation in Thermo-Elastic Solids
,”
Comput. Methods Appl. Mech. Eng.
,
294
, pp.
449
485
.
33.
Stocek
,
J.
,
Arthern
,
R. J
, and
Marsh
,
O. J.
,
2023
, “A Viscoelastic Phase-Field Model for Iceberg Calving,” arXiv preprint arXiv:2306.08405.
34.
Miehe
,
C.
,
Aldakheel
,
F.
, and
Teichtmeister
,
S.
,
2017
, “
Phase-Ffield Modeling of Ductile Fracture at Finite Strains: A Robust Variational-Based Numerical Implementation of a Gradient-Extended Theory by Micromorphic Regularization
,”
Int. J. Numer. Methods Eng.
,
111
(
9
), pp.
816
863
.
35.
Logg
,
A.
,
Mardal
,
K.-A.
, and
Wells
,
G.
,
2012
,
Automated Solution of Differential Equations by the Finite Element Method: The FEniCS Book
, Vol.
84
,
Springer Science & Business Media
.
36.
Barchiesi
,
E.
,
Yang
,
H.
,
Tran
,
C. A.
,
Placidi
,
L.
, and
Müller
,
W. H.
,
2021
, “
Computation of Brittle Fracture Propagation in Strain Gradient Materials by the FEniCS Library
,”
Math. Mech. Solids
,
26
(
3
), pp.
325
340
.
37.
Hirshikesh
,
Natarajan
,
S.
, and
Annabattula
,
R. K.
,
2019
, “
A Fenics Implementation of the Phase Field Method for Quasi-Static Brittle Fracture
,”
Front. Struct. Civ. Eng.
,
13
, pp.
380
396
.
38.
Hrishikesh
,
Natarajan
,
S.
,
Annabattula
,
R. K.
, and
Martínez-Pañeda
,
E.
,
2019
, “
Phase Field Modelling of Crack Propagation in Functionally Graded Materials
,”
Compos. B Eng.
,
169
, pp.
239
248
.
39.
Karimi
,
M.
,
Massoudi
,
M.
,
Walkington
,
N.
,
Pozzi
,
M.
, and
Dayal
,
K.
,
2022
, “
Energetic Formulation of Large-Deformation Poroelasticity
,”
Int. J. Numer. Anal. Methods Geomech.
,
46
(
5
), pp.
910
932
.
40.
Karimi
,
M.
,
Massoudi
,
M.
,
Dayal
,
K.
, and
Pozzi
,
M.
,
2023
, “
High-Dimensional Nonlinear Bayesian Inference of Poroelastic Fields From Pressure Data
,”
Math. Mech. Solids
,
28
(
9
), pp.
2108
2131
.
41.
Chua
,
J.
,
Karimi
,
M.
,
Kozlowski
,
P.
,
Massoudi
,
M.
,
Narasimhachary
,
S.
,
Kadau
,
K.
,
Gazonas
,
G.
, and
Dayal
,
K.
,
2024
, “
Deformation Decomposition Versus Energy Decomposition for Chemo-and Poro-Mechanics
,”
ASME J. Appl. Mech.
,
91
(
1
), p.
014501
.
42.
Duddu
,
R.
,
Jiménez
,
S.
, and
Bassis
,
J.
,
2020
, “
A Non-Local Continuum Poro-Damage Mechanics Model for Hydrofracturing of Surface Crevasses in Grounded Glaciers
,”
J. Glaciol.
,
66
(
257
), pp.
415
429
.
43.
Agrawal
,
V.
, and
Dayal
,
K.
,
2015
, “
A Dynamic Phase-Field Model for Structural Transformations and Twinning: Regularized Interfaces With Transparent Prescription of Complex Kinetics and Nucleation. Part I: Formulation and One-Dimensional Characterization
,”
J. Mech. Phys. Solids
,
85
, p.
270
.
44.
Agrawal
,
V.
, and
Dayal
,
K.
,
2015
, “
A Dynamic Phase-Field Model for Structural Transformations and Twinning: Regularized Interfaces With Transparent Prescription of Complex Kinetics and Nucleation. Part II: Two-Dimensional Characterization and Boundary Kinetics
,”
J. Mech. Phys. Solids
,
85
, p.
291
.
45.
Chua
,
J.
,
Agrawal
,
V.
,
Breitzman
,
T.
,
Gazonas
,
G.
, and
Dayal
,
K.
,
2022
, “
Phase-Field Modeling and Peridynamics for Defect Dynamics, and an Augmented Phase-Field Model With Viscous Stresses
,”
J. Mech. Phys. Solids
,
159
, p.
104716
.
46.
Chua
,
J.
,
Agrawal
,
V.
,
Walkington
,
N.
,
Gazonas
,
G.
, and
Dayal
,
K.
,
2024
, “
Interplay Between Nucleation and Kinetics in Dynamic Twinning
,”
ASME J. Appl. Mech.
,
91
(
12
), p.
121002
.
47.
Agrawal
,
V.
, and
Dayal
,
K.
,
2017
, “
Dependence of Equilibrium Griffith Surface Energy on Crack Speed in Phase-Field Models for Fracture Coupled to Elastodynamics
,”
Int. J. Fract.
,
207
(
2
), pp.
243
249
.
48.
Clayton
,
J. D.
, and
Knap
,
J.
,
2014
, “
A Geometrically Nonlinear Phase Field Theory of Brittle Fracture
,”
Int. J. Fract.
,
189
(
2
), pp.
139
148
.
49.
Clayton
,
J. D.
, and
Knap
,
J.
,
2015
, “
Nonlinear Phase Field Theory for Fracture and Twinning With Analysis of Simple Shear
,”
Philos. Mag.
,
95
(
24
), pp.
2661
2696
.
50.
Mohammad Mousavi
,
S.
,
Ang
,
I.
,
Mulderrig
,
J.
, and
Bouklas
,
N.
,
2024
, “
Evaluating Fracture Energy Predictions Using Phase-Field and Gradient-Enhanced Damage Models for Elastomers
,”
ASME J. Appl. Mech.
,
91
(
12
), p.
121009
.
51.
Mousavi
,
S. M.
,
Mulderrig
,
J.
,
Talamini
,
B.
, and
Bouklas
,
N.
,
2025
, “A Chain Stretch-Based Gradient-Enhanced Model for Damage and Fracture in Elastomers,”
arXiv preprint
https://arxiv.org/abs/2502.02822.