- Award ID(s):
- 1903308
- NSF-PAR ID:
- 10167757
- Date Published:
- Journal Name:
- Soft Matter
- ISSN:
- 1744-683X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
null (Ed.)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
-
Energy dissipation around a propagating crack is the primary mechanism for the enhanced fracture toughness in viscoelastic solids. Such dissipation is spatially non-uniform and is highly coupled to the crack propagation process due to the history-dependent nature of viscoelasticity. We present an experimental approach to map the dissipation field during crack propagation in soft viscoelastic solid. Specifically, we track randomly distributed tracer particles to measure the evolving deformation field. The measured deformation field is then put into a nonlinear constitutive model to determine the dissipation field. Our methodology was used to investigate the deformation and dissipation fields around a propagating crack in a Polyampholyte (PA) hydrogel. The deformation field measurements allowed us to assess whether the commonly assumed translational invariance in viscoelastic fracture theories holds true in practical experiments. Furthermore, by combining the obtained deformation fields with a nonlinear viscoelastic model, we captured the complete history of the dissipation field during crack propagation. We found that dissipation occurred even at material points that are a few millimeters away from the crack tip. The mapped dissipation field also enabled the separate determination of the intrinsic and dissipative components of fracture toughness for the viscoelastic hydrogel.more » « less
-
null (Ed.)Dynamic networks contain crosslinks that re-associate after disconnecting, imparting them with viscoelastic properties. While continuum approaches have been developed to analyze their mechanical response, these approaches can only describe their evolution in an average sense, omitting local, stochastic mechanisms that are critical to damage initiation or strain localization. To address these limitations, we introduce a discrete numerical model that mesoscopically coarse-grains the individual constituents of a dynamic network to predict its mechanical and topological evolution. Each constituent consists of a set of flexible chains that are permanently cross-linked at one end and contain reversible binding sites at their free ends. We incorporate nonlinear force–extension of individual chains via a Langevin model, slip-bond dissociation through Eyring's model, and spatiotemporally-dependent bond attachment based on scaling theory. Applying incompressible, uniaxial tension to representative volume elements at a range of constant strain rates and network connectivities, we then compare the mechanical response of these networks to that predicted by the transient network theory. Ultimately, we find that the idealized continuum approach remains suitable for networks with high chain concentrations when deformed at low strain rates, yet the mesoscale model proves necessary for the exploration of localized stochastic events, such as variability of the bond kinetics, or the nucleation of micro-cavities that likely conceive damage and fracture.more » « less
-
Abstract The atomic‐scale cracking mechanism in clay is vital in discovering the cracking mechanism of clay at the continuum scale in that clay is a nanomaterial. In this article, we investigate mechanisms of modes I and II crack propagations in pyrophyllite and Ca‐montmorillonite with water adsorption through reactive molecular dynamics (MD) with a bond‐order force field. Clay water adsorption is considered by adding water molecules to the clay surface. During the equilibration stage, water adsorption could cause bending deformation of the predefined edge crack region. The relatively small orientating angle of water molecules indicates the formation of hydrogen bonds in the crack propagation process. The peak number density of adsorbed water decreases with the increasing strains. The atomistic structure evolution of the crack tip under loading is analyzed to interpret the nanoscale crack propagation mechanism. The numerical results show that the crack tip first gets blunted with a significant increase in the radius of the curvature of the crack tip and a slight change in crack length. The crack tip blunting is studied by tracking the crack tip opening distance and O–Si–O angle in the tetrahedral Si–O cell in modes I and II cracks. We compare bond‐breaking behaviors between Al–O and Si–O. It is found that Si–O bond breaking is primarily responsible for crack propagation. The critical stress intensity factor and critical energy release rate are determined from MD simulation results.
-
Abstract The recently conceived gap test and its simulation revealed that the fracture energy Gf (or Kc, Jcr) of concrete, plastic-hardening metals, composites, and probably most materials can change by ±100%, depending on the crack-parallel stresses σxx, σzz, and their history. Therefore, one must consider not only a finite length but also a finite width of the fracture process zone, along with its tensorial damage behavior. The data from this test, along with ten other classical tests important for fracture problems (nine on concrete, one on sandstone), are optimally fitted to evaluate the performance of the state-of-art phase-field, peridynamic, and crack band models. Thanks to its realistic boundary and crack-face conditions as well as its tensorial nature, the crack band model, combined with the microplane damage constitutive law in its latest version M7, is found to fit all data well. On the contrary, the phase-field models perform poorly. Peridynamic models (both bond based and state based) perform even worse. The recent correction in the bond-associated deformation gradient helps to improve the predictions in some experiments, but not all. This confirms the previous strictly theoretical critique (JAM 2016), which showed that peridynamics of all kinds suffers from several conceptual faults: (1) It implies a lattice microstructure; (2) its particle–skipping interactions are a fiction; (4) it ignores shear-resisted particle rotations (which are what lends the lattice discrete particle model (LDPM) its superior performance); (3) its representation of the boundaries, especially the crack and fracture process zone faces, is physically unrealistic; and (5) it cannot reproduce the transitional size effect—a quintessential characteristic of quasibrittleness. The misleading practice of “verifying” a model with only one or two simple tests matchable by many different models, or showcasing an ad hoc improvement for one type of test while ignoring misfits of others, is pointed out. In closing, the ubiquity of crack-parallel stresses in practical problems of concrete, shale, fiber composites, plastic-hardening metals, and materials on submicrometer scale is emphasized.more » « less