skip to main content


Title: New perspective of fracture mechanics inspired by gap test with crack-parallel compression

The line crack models, including linear elastic fracture mechanics (LEFM), cohesive crack model (CCM), and extended finite element method (XFEM), rest on the century-old hypothesis of constancy of materials’ fracture energy. However, the type of fracture test presented here, named the gap test, reveals that, in concrete and probably all quasibrittle materials, including coarse-grained ceramics, rocks, stiff foams, fiber composites, wood, and sea ice, the effective mode I fracture energy depends strongly on the crack-parallel normal stress, in-plane or out-of-plane. This stress can double the fracture energy or reduce it to zero. Why hasn’t this been detected earlier? Because the crack-parallel stress in all standard fracture specimens is negligible, and is, anyway, unaccountable by line crack models. To simulate this phenomenon by finite elements (FE), the fracture process zone must have a finite width, and must be characterized by a realistic tensorial softening damage model whose vectorial constitutive law captures oriented mesoscale frictional slip, microcrack opening, and splitting with microbuckling. This is best accomplished by the FE crack band model which, when coupled with microplane model M7, fits the test results satisfactorily. The lattice discrete particle model also works. However, the scalar stress–displacement softening law of CCM and tensorial models with a single-parameter damage law are inadequate. The experiment is proposed as a standard. It represents a simple modification of the three-point-bend test in which both the bending and crack-parallel compression are statically determinate. Finally, a perspective of various far-reaching consequences and limitations of CCM, LEFM, and XFEM is discussed.

 
more » « less
Award ID(s):
2029641
NSF-PAR ID:
10160406
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Proceedings of the National Academy of Sciences
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Volume:
117
Issue:
25
ISSN:
0027-8424
Page Range / eLocation ID:
p. 14015-14020
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract In the standard fracture test specimens, the crack-parallel normal stress is negligible. However, its effect can be strong, as revealed by a new type of experiment, briefly named the gap test. It consists of a simple modification of the standard three-point-bend test whose main idea is to use plastic pads with a near-perfect yield plateau to generate a constant crack-parallel compression and install the end supports with a gap that closes only when the pads yield. This way, the test beam transits from one statically determinate loading configuration to another, making evaluation unambiguous. For concrete, the gap test showed that moderate crack-parallel compressive stress can increase up to 1.8 times the Mode I (opening) fracture energy of concrete, and reduce it to almost zero on approach to the compressive stress limit. To model it, the fracture process zone must be characterized tensorially. We use computer simulations with crack-band microplane model, considering both in-plane and out-of-plane crack-parallel stresses for plain and fiber-reinforced concretes, and anisotropic shale. The results have broad implications for all quasibrittle materials, including shale, fiber composites, coarse ceramics, sea ice, foams, and fone. Except for negligible crack-parallel stress, the line crack models are shown to be inapplicable. Nevertheless, as an approximation ignoring stress tensor history, the crack-parallel stress effect may be introduced parametrically, by a formula. Finally we show that the standard tensorial strength models such as Drucker–Prager cannot reproduce these effects realistically. 
    more » « less
  2. 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
  3. Abstract

    This work presents a stabilized formulation for phase‐field fracture of hyperelastic materials near the limit of incompressibility. At this limit, traditional mixed displacement and pressure formulations must satisfy the inf‐sup condition for solution stability. The mixed formulation coupled with the damage field can lead to an inhibition of crack opening as volumetric changes are severely penalized effectively creating a pressure‐bubble. To overcome this bottleneck, we utilize a mixed formulation with a perturbed Lagrangian formulation which enforces the incompressibility constraint in the undamaged material and reduces the pressure effect in the damaged material. A mesh‐dependent stabilization technique based on the residuals of the Euler–Lagrange equations multiplied with a differential operator acting on the weight space is used, allowing for linear interpolation of all field variables of the elastic subproblem. This formulation was validated with three examples at finite deformations: a plane‐stress pure‐shear test, a two‐dimensional geometry in plane‐stress, and a three‐dimensional notched sample. In the last example, we incorporate a hybrid formulation with an additive strain energy decomposition to account for different behaviors in tension and compression. The results show close agreement with analytical solutions for crack tip opening displacements and performs well at the limit of incompressibility.

     
    more » « less
  4. 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 breaking 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. 
    more » « less
  5. ABSTRACT We investigate the appropriateness of calving or crevasse models from the literature using linear elastic fracture mechanics (LEFM). To this end, we compare LEFM model-predicted stress intensity factors (SIFs) against numerically computed SIFs using the displacement correlation method in conjunction with the finite element method. We present several benchmark simulations wherein we calculate the SIF at the tips of water-filled surface and basal crevasses penetrating through rectangular ice slabs under different boundary conditions, including grounded and floating conditions. Our simulation results indicate that the basal boundary condition significantly influences the SIF at the crevasse tips. We find that the existing calving models using LEFM are not generally accurate for evaluating SIFs in grounded glaciers or floating ice shelves. We also illustrate that using the ‘single edge crack’ weight function in the LEFM formulations may be appropriate for predicting calving from floating ice shelves, owing to the low fracture toughness of ice; whereas, using the ‘double edge crack’ or ‘central through crack’ weight functions is more appropriate for predicting calving from grounded glaciers. To conclude, we recommend using the displacement correlation method for SIF evaluation in real glaciers and ice shelves with complex geometries and boundary conditions. 
    more » « less