skip to main content

Title: 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) more » 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. « less
; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Journal of the mechanics and physics of solids
Sponsoring Org:
National Science Foundation
More Like this
  1. Double-network gels are a class of tough soft materials comprising two elastic networks with contrasting structures. The formation of a large internal damage zone ahead of the crack tip by the rupturing of the brittle network accounts for the large crack resistance of the materials. Understanding what determines the damage zone is the central question of the fracture mechanics of double-network gels. In this work, we found that at the onset of crack propagation, the size of necking zone, in which the brittle network breaks into fragments and the stretchable network is highly stretched, distinctly decreases with the increase of the solvent viscosity, resulting in a reduction in the fracture toughness of the material. This is in sharp contrast to the tensile behavior of the material that does not change with the solvent viscosity. This result suggests that the dynamics of stretchable network strands, triggered by the rupture of the brittle network, plays a role. To account for this solvent viscosity effect on the crack initiation, a delayed blunting mechanism regarding the polymer dynamics effect is proposed. The discovery on the role of the polymer dynamic adds an important missing piece to the fracture mechanism of this unique material.
  2. 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.
  3. Toughness of soft materials such as elastomers and gels depends on their ability to dissipate energy and to reduce stress concentration at the crack tip. The primary energy dissipation mechanism is viscoelasticity. Most analyses and models of fracture are based on linear viscoelastic theory (LVT) where strains are assumed to be small and the relaxation mechanisms are independent of stress or strain history. A well-known paradox is that the size of the dissipative zone predicted by LVT is unrealistically small. Here we use a physically based nonlinear viscoelastic model to illustrate why the linear theory breaks down. Using this nonlinear model and analogs of crack problems, we give a plausible resolution to this paradox. In our model, viscoelasticity arises from the breaking and healing of physical cross-links in the polymer network. When the deformation is small, the kinetics of bond breaking and healing are independent of the strain/stress history and the model reduces to the standard linear theory. For large deformations, localized bond breaking damages the material near the crack tip, reducing stress concentration and dissipating energy at the same time. The damage zone size is a new length scale which depends on the strain required to accelerate bond breakingmore »kinetics. These effects are illustrated by considering two cases with stress concentrations: the evolution of spherical damage in a viscoelastic body subjected to internal pressure, and a zero degree peel test.« less
  4. 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 wheremore »the propagating cracks slowed down, deviated, and/or stopped, improving the resistance of the shape memory alloy specimen to fracture.

    « less
  5. Fracturing in brittle rocks with an existing crack results in the development of a significant nonlinear region surrounding the crack tip called the fracture process zone. Various experimental and numerical studies have shown that the crack tip parameters such as the crack tip opening displacement (CTOD) and the fracture energy are critically important in characterizing the fracture process zone. In this study, numerical simulations of rock specimens with a center notch subjected to three-point bending were conducted using the extended finite element method (XFEM) along with the cohesive zone model (CZM) to account for fracture process zone. The input parameters of CZM such as the elastic and critical crack opening displacements were first estimated based on the results of three-point bending tests on the center notched Barre granite specimens. Displacements were measured using the two dimensional digital image correlation technique and used to characterize the evolution of the fracture process zone and estimate the parameters of the cohesive zone model. The results from the numerical simulations showed that CZM provided a good agreement with experimental data as it predicted all three stages of cracking from fracture process initiation to macro-crack growth.