An official website of the United States government
Abstract 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.
more »
« less
Oscillatory and tip-splitting instabilities in 2D dynamic fracture: The roles of intrinsic material length and time scales
Recent theoretical and computational progress has led to unprecedented understanding of symmetry-breaking instabilities in 2D dynamic fracture. At the heart of this progress resides the identification of two intrinsic, near crack tip length scales — a nonlinear elastic length scale ℓ and a dissipation length scale ξ — that do not exist in Linear Elastic Fracture Mechanics (LEFM), the classical theory of cracks. In particular, it has been shown that at a propagation velocity v of about 90% of the shear wave-speed, cracks in 2D brittle materials undergo an oscillatory instability whose wavelength varies linearly with ℓ, and at larger loading levels (corresponding to yet higher propagation velocities), a tip-splitting instability emerges, both in agreements with experiments. In this paper, using phase-field models of brittle fracture, we demonstrate the following properties of the oscillatory instability: (i) It exists also in the absence of near-tip elastic nonlinearity, i.e. in the limit ℓ→0, with a wavelength determined by the dissipation length scale ξ. This result shows that the instability crucially depends on the existence of an intrinsic length scale associated with the breakdown of linear elasticity near crack tips, independently of whether the latter is related to nonlinear elasticity or to dissipation. (ii) It is a supercritical Hopf bifurcation, featuring a vanishing oscillations amplitude at onset. (iii) It is largely independent of the phenomenological forms of the degradation functions assumed in the phase-field framework to describe the cohesive zone, and of the velocity-dependence of the fracture energy Γ(v) that is controlled by the dissipation time scale in the Ginzburg-Landau-type evolution equation for the phase-field. These results substantiate the universal nature of the oscillatory instability in 2D. In addition, we provide evidence indicating that the tip-splitting instability is controlled by the limiting rate of elastic energy transport inside the crack tip region. The latter is sensitive to the wave-speed inside the dissipation zone, which can be systematically varied within the phase-field approach. Finally, we describe in detail the numerical implementation scheme of the employed phase-field fracture approach, allowing its application in a broad range of materials failure problems.
more »
« less
- 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
-
-
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
-
ABSTRACT 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
-
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
-
NA (Ed.)We derive sharp-interface models for one-dimensional brittle fracture via the inverse-deformation approach. Methods of Γ -convergence are employed to obtain the singular limits of previously proposed models. The latter feature a local, non-convex stored energy of inverse strain, augmented by small interfacial energy, formulated in terms of the inverse-strain gradient. They predict spontaneous fracture with exact crack-opening discontinuities, without the use of damage (phase) fields or pre-existing cracks; crack faces are endowed with a thin layer of surface energy. The models obtained herewith inherit the same properties, except that surface energy is now concentrated at the crack faces in the Γ -limit. Accordingly, we construct energy-minimizing configurations. For a composite bar with a breakable layer, our results predict a pattern of equally spaced cracks whose number is given as an increasing function of applied load.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/104372
