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
Phase field approach for nanoscale interactions between crack propagation and phase transformation
The phase field approach (PFA) for the interaction of fracture and martensitic phase transformation (PT) is developed, which includes the change in surface energy during PT and the effect of unexplored scale parameters proportional to the ratio of the widths of the crack surface and the phase interface, both at the nanometer scale. The variation of these two parameters causes unexpected qualitative and quantitative effects: shift of PT away from the crack tip, “wetting” of the crack surface by martensite, change in the structure and geometry of the transformed region, crack trajectory, and process of interfacial damage evolution, as well as transformation toughening. The results suggest additional parameters controlling coupled fracture and PTs.
more »
« less
- Award ID(s):
- 1904830
- PAR ID:
- 10146395
- Date Published:
- Journal Name:
- Nanoscale
- Volume:
- 11
- Issue:
- 46
- ISSN:
- 2040-3364
- Page Range / eLocation ID:
- 22243 to 22247
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
A phase-field model for thermomechanically-induced fracture in NiTi at the single crystal level, i.e., fracture under loading paths that may take advantage of either of the functional properties of NiTi–superelasticity or shape memory effect–, is presented, formulated within the kinematically linear regime. The model accounts for reversible phase transformation from austenite to martensite habit plane variants and plastic deformation in the austenite phase. Transformation-induced plastic deformation is viewed as a mechanism for accommodation of the local deformation incompatibility at the austenite–martensite interfaces and is accounted for by introducing an interaction term in the free energy derived based on the Mori–Tanaka and Kröner micromechanical assumptions and the hypothesis of martensite instantaneous growth within austenite. Based on experimental observations suggesting that NiTi fractures in a stress-controlled manner, damage is assumed to be driven by the elastic energy, i.e., phase transformation and plastic deformation are assumed to contribute in crack formation and growth indirectly through stress redistribution. The model is restricted to quasistatic mechanical loading (no latent heat effects), thermal loading sufficiently slow with respect to the time rate of heat transfer by conduction (no thermal gradients), and a temperature range below 𝑀𝑑, which is the temperature above which the austenite phase is stable, i.e., stress-induced martensitic transformation is suppressed. The numerical implementation of the model is based on an efficient scheme of viscous regularization in both phase transformation and plastic deformation, an explicit numerical integration via a tangent modulus method, and a staggered scheme for the coupling of the unknown fields. The model is shown able to capture transformation-induced toughening, i.e., stable crack advance attributed to the shielding effect of inelastic deformation left in the wake of the growing crack under nominal isothermal loading, actuation-induced fracture under a constant bias load, and crystallographic dependence on crack pattern.more » « less
-
Athermal resistance to the motion of a phase interface due to a precipitate is investigated. The coupled phase field and elasticity equations are solved for the phase transformation (PT). The volumetric misfit strain due to the precipitate is included using the error and rectangular functions. Due to the presence of precipitates, the critical thermal driving forces remarkably differ between the direct and reverse PTs, resulting in a hysteresis behavior. For the precipitate radius small compared to the interface width, the misfit strain does not practically show any effect on the critical thermal driving force. Also, the critical thermal driving force value nonlinearly increases vs. the precipitate concentration for both the direct and reverse PTs. Change in the precipitate surface energy significantly changes the PT morphology and the critical thermal driving forces. The critical thermal driving force shows dependence on the misfit strain for large precipitate sizes compared to the interface width. For both the constant surface energy (CSE) and variable surface energy (VSE) boundary conditions (BCs) at the precipitate surface, the critical thermal driving force linearly increases vs. the misfit strain coefficient for the direct PT while it is almost independent of it for the reverse PT. For larger precipitates, the critical thermal driving force nonlinearly increases vs. the precipitate concentration for the direct PT. For the reverse PT, its value for the CSE BCs linearly increases vs. the precipitate concentration while it is independent of the precipitate concentration for the VSE BCs. Also, for any concentration, the VSE BCs result in higher thermal critical driving forces, a smaller hysteresis range, and a larger transformation rate. The critical microstructure and thermal driving forces are validated using the thermodynamic phase equilibrium condition for stationary interfaces.more » « less
-
Athermal resistance to the motion of a phase interface due to a precipitate is investigated. The coupled phase field and elasticity equations are solved for the phase transformation (PT). The volumetric misfit strain due to the precipitate is included using the error and rectangular functions. Due to the presence of precipitates, the critical thermal driving forces remarkably differ between the direct and reverse PTs, resulting in a hysteresis behavior. For the precipitate radius small compared to the interface width, the misfit strain does not practically show any effect on the critical thermal driving force. Also, the critical thermal driving force value nonlinearly increases vs. the precipitate concentration for both the direct and reverse PTs. Change in the precipitate surface energy significantly changes the PT morphology and the critical thermal driving forces. The critical thermal driving force shows dependence on the misfit strain for large precipitate sizes compared to the interface width. For both the constant surface energy (CSE) and variable surface energy (VSE) boundary conditions (BCs) at the precipitate surface, the critical thermal driving force linearly increases vs. the misfit strain coefficient for the direct PT while it is almost independent of it for the reverse PT. For larger precipitates, the critical thermal driving force nonlinearly increases vs. the precipitate concentration for the direct PT. For the reverse PT, its value for the CSE BCs linearly increases vs. the precipitate concentration while it is independent of the precipitate concentration for the VSE BCs. Also, for any concentration, the VSE BCs result in higher thermal critical driving forces, a smaller hysteresis range, and a larger transformation rate. The critical microstructure and thermal driving forces are validated using the thermodynamic phase equilibrium condition for stationary interfaces.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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/C9NR05960A
