An official website of the United States government
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
Critical Comparison of Phase-Field, Peridynamics, and Crack Band Model M7 in Light of Gap Test and Classical Fracture Tests
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
- Award ID(s):
- 2029641
- PAR ID:
- 10346006
- Date Published:
- Journal Name:
- Journal of Applied Mechanics
- Volume:
- 89
- Issue:
- 6
- ISSN:
- 0021-8936
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
The 2023 smooth Lagrangian Crack-Band Model (slCBM), inspired by the 2020 invention of the gap test, prevented spurious damage localization during fracture growth by introducing the second gradient of the displacement field vector, named the “sprain,” as the localization limiter. The key idea was that, in the finite element implementation, the displacement vector and its gradient should be treated as independent fields with the lowest ( ) continuity, constrained by a second-order Lagrange multiplier tensor. Coupled with a realistic constitutive law for triaxial softening damage, such as microplane model M7, the known limitations of the classical Crack Band Model were eliminated. Here, we show that the slCBM closely reproduces the size effect revealed by the gap test at various crack-parallel stresses. To describe it, we present an approximate corrective formula, although a strong loading-path dependence limits its applicability. Except for the rare case of zero crack-parallel stresses, the fracture predictions of the line crack models (linear elastic fracture mechanics, phase-field, extended finite element method (XFEM), cohesive crack models) can be as much as 100% in error. We argue that the localization limiter concept must be extended by including the resistance to material rotation gradients. We also show that, without this resistance, the existing strain-gradient damage theories may predict a wrong fracture pattern and have, for Mode II and III fractures, a load capacity error as much as 55%. Finally, we argue that the crack-parallel stress effect must occur in all materials, ranging from concrete to atomistically sharp cracks in crystals.more » « less
-
Abstract Whereas various simplistic microplane models of limited applicability, defined by stress–strain curves on the microplane, can function as either explicit or implicit, the explicit‐to‐implicit conversion of realistic versatile microplane models for plain or fiber‐reinforced concrete, shale and composites has remained a challenge for quarter century. The reason is that these realistic models use microplane stress–strain boundaries defined by inequalities. Here, we show how the conversion can be easily achieved on the microplane level and then transferred to a tangent stiffness tensor or an inelastic stiffness tensor to be used in Newton–Raphson iterations within a loading step. To ensure convergence, a minor adjustment in the M7 algorithm is introduced to achieve continuity. Power‐law convergence, almost quadratic in most cases, is also demonstrated. Seven examples of crack‐band finite element simulations of challenging laboratory tests document nearly identical implicit and explicit results, as well as good match of test data. Three of them, including the vertex effect in compression‐torsion tests, pure Mode II shear fracture, and the “gap test” of the crack‐parallel compression effect on Mode I load‐deflection curve, have not been reproduced by other models before. The coding of implicit M7 subroutine, usable in, for example, UMAT of ABAQUS, is posted for a free download.more » « less
-
Abstract The crack band model, which was shown to provide a superior computational representation of fracture of quasibrittle materials (in this journal, May 2022), still suffers from three limitations: (1) The material damage is forced to be uniform across a one-element wide band because of unrestricted strain localization instability; (2) the width of the fracture process zone is fixed as the width of a single element; and (3) cracks inclined to rectangular mesh lines are represented by a rough zig-zag damage band. Presented is a generalization that overcomes all three, by enforcing a variable multi-element width of the crack band front controlled by a material characteristic length l0. This is achieved by introducing a homogenized localization energy density that increases, after a certain threshold, as a function of an invariant of the third-order tensor of second gradient of the displacement vector, called the sprain tensorη, representing (in isotropic materials) the magnitude of its Laplacian (not expressible as a strain-gradient tensor). The continuum free energy density must be augmented by additional sprain energy Φ(l0η), which affects only the postpeak softening damage. In finite element discretization, the localization resistance is effected by applying triplets of self-equilibrated in-plane nodal forces, which follow as partial derivatives of Φ(l0η). The force triplets enforce a variable multi-element crack band width. The damage distribution across the fracture process zone is non-uniform but smoothed. The standard boundary conditions of the finite element method apply. Numerical simulations document that the crack band propagates through regular rectangular meshes with virtually no directional bias.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1115/1.4054221
