Unlike micromechanics failure models that have a well-defined crack path, phase-field fracture models are capable of predicting the crack path in arbitrary geometries and dimensions by utilizing a diffuse representation of cracks. However, such models rely on the calibration of a fracture energy (Gc) and a regularization length-scale (lc) parameter, which do not have a strong micromechanical basis. Here, we construct the equivalent crack-tip cohesive zone laws representing a phase-field fracture model, to elucidate the effects of Gc and lc on the fracture resistance and crack growth mechanics under mode I K-field loading. Our results show that the cohesive zone law scales with increasing Gc while maintaining the same functional form. In contrast, increasing lc broadens the process zone and results in a flattened traction-separation profile with a decreased but sustained peak cohesive traction over longer separation distances. While Gc quantitatively captures the fracture initiation toughness, increasing Gc coupled with decreasing lc contributes to a rising fracture resistance curve and a higher steady-state toughness—both these effects cumulate in an evolving cohesive zone law with crack progression. We discuss the relationship between these phase-field parameters and process zone characteristics in the material.
- Award ID(s):
- 1827343
- PAR ID:
- 10216227
- Date Published:
- Journal Name:
- Journal of the mechanics and physics of solids
- ISSN:
- 0022-5096
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract -
Abstract There are several issues to be solved in the fracture mechanics of shape memory alloys, one of them being the resistance to crack growth and therefore to fracture. This paper discusses the crack growth in a single crystal CoNiAl shape memory alloy under cyclic loading and the effect of micro-structural barriers. To observe the crack growth in detail, tests are conducted on edge-notched specimens. The displacement field is obtained using digital image correlation (DIC), and the fracture parameters are calculated by fitting anisotropic crack tip displacement equations to DIC data. Similar crack growth behaviors are observed in both superelastic and shape memory specimens, with a comparatively higher crack growth rate in the superelastic case: first a crack initiates at the notch and grows, then new cracks are observed to form near the tip of the main crack, or on the notch when the growth slows down. Then, further cyclic loading leads to the growth of the main crack and the new crack simultaneously with the two cracks merging at the end. Test specimens are examined post-failure with optical microscopy to better understand this complicated behavior. Results showed the presence of a non-transforming secondary (
γ ) phase around the regions where the propagating cracks slowed down, deviated, and/or stopped, improving the resistance of the shape memory alloy specimen to fracture. -
Strongly anisotropic geomaterials, such as layered shales, have been observed to undergo fracture under compressive loading. This paper applies a phase‐field fracture model to study this fracture process. While phase‐field fracture models have several advantages—primarily that the fracture path is not predetermined but arises naturally from the evolution of a smooth non‐singular damage field—they provide unphysical predictions when the stress state is complex and includes compression that can cause crack faces to contact. Building on a recently developed phase‐field model that accounts for compressive traction across the crack face, this paper extends the model to the setting of anisotropic fracture. The key features of the model include the following: (1) a homogenized anisotropic elastic response and strongly anisotropic model for the work to fracture; (2) an effective damage response that accounts consistently for compressive traction across the crack face, that is derived from the anisotropic elastic response; (3) a regularized crack normal field that overcomes the shortcomings of the isotropic setting, and enables the correct crack response, both across and transverse to the crack face. To test the model, we first compare the predictions to phase‐field fracture evolution calculations in a fully resolved layered specimen with spatial inhomogeneity, and show that it captures the overall patterns of crack growth. We then apply the model to previously reported experimental observations of fracture evolution in laboratory specimens of shales under compression with confinement, and find that it predicts well the observed crack patterns in a broad range of loading conditions. We further apply the model to predict the growth of wing cracks under compression and confinement. Prior approaches to simulate wing cracks have treated the initial cracks as an external boundary, which makes them difficult to apply to general settings. Here, the effective crack response model enables us to treat the initial crack simply as a nonsingular damaged zone within the computational domain, thereby allowing for easy and general computations.more » « less
-
Recent works have shown that in contrast to classical linear elastic fracture mechanics, endowing crack fronts in a brittle Green-elastic solid with Steigmann-Ogden surface elasticity yields a model that predicts bounded stresses and strains at the crack tips for plane-strain problems. However, singularities persist for anti-plane shear (mode-III fracture) under far field loading, even when Steigmann-Ogden surface elasticity is incorporated. This work is motivated by obtaining a model of brittle fracture capable of predicting bounded stresses and strains for all modes of loading. We formulate an exact general theory of a three-dimensional solid containing a boundary surface with strain-gradient surface elasticity. For planar reference surfaces parameterized by flat coordinates, the form of surface elasticity reduces to that introduced by Hilgers and Pipkin, and when the surface energy is independent of the surface covariant derivative of the stretching, the theory reduces to that of Steigmann and Ogden. We discuss material symmetry using Murdoch and Cohen’s extension of Noll’s theory. We present a model small-strain surface energy that incorporates resistance to geodesic distortion, satisfies strong ellipticity, and requires the same material constants found in the Steigmann-Ogden theory. Finally, we derive and apply the linearized theory to mode-III fracture in an infinite plate under far-field loading. We prove that there always exists a unique classical solution to the governing integro-differential equation, and in contrast to using Steigmann-Ogden surface elasticity, our model is consistent with the linearization assumption in predicting finite stresses and strains at the crack tips.more » « less
-
Crack initiation emerges due to a combination of elasticity, plasticity, and disorder, and it displays strong dependence on the material’s microstructural details. The characterization of the structural uncertainty in the original microstructure is typically empirical and systematic characterization protocols are lacking. In this paper, we propose an investigational tool in the form of the curvature an ellipsoidal notch: As the radius of curvature at the notch increases, there is a dynamic phase transition from notch-induced crack initiation to bulk-disorder crack nucleation. The notch length scale associated with this transition may provide an additional characteristic of the original material microstructure. We investigate brittle but elastoplastic metals with coarse-grained, microstructural disorder that could originate in a material’s manufacturing process, such as alloying. We perform extensive and realistic simulations using a phase-field approach coupled to crystal plasticity. The microstructural disorder and notch width are systematically varied. We identify this transition for various disorder strengths in terms of the damage evolution. We identify detectable precursors to crack initiation that we quantify in terms of the expected stress drops during mode I fracture loading. Finally, we discuss ways to observe and analyze this brittle to quasi-brittle transition in experiments.more » « less