An official website of the United States government
Abstract A preceding 2023 study argued that the resistance of a heterogeneous material to the curvature of the displacement field is the most physically realistic localization limiter for softening damage. The curvature was characterized by the second gradient of the displacement vector field, which includes the material rotation gradient, and was named the “sprain” tensor, while the term “spress” is here proposed as the force variable work-conjugate to “sprain.” The partial derivatives of the associated sprain energy density yielded in the preceeding study, sets of curvature resisting self-equilibrated nodal sprain forces. However, the fact that the sprain forces had to be applied on the adjacent nodes of a finite element greatly complicated the programming and extended the simulation time in a commercial code such as abaqus by almost two orders of magnitude. In the present model, Smooth Lagrangian Crack Band Model (slCBM), these computational obstacles are here overcome by using finite elements with linear shape functions for both the displacement vector and for an approximate displacement gradient tensor. A crucial feature is that the nodal values of the approximate gradient tensor are shared by adjacent finite elements. The actual displacement gradient tensor calculated from the nodal displacement vectors is constrained to the approximate displacement gradient tensor by means of a Lagrange multiplier tensor, either one for each element or one for each node. The gradient tensor of the approximate gradient tensor then represents the approximate third-order displacement curvature tensor, or Hessian of the displacement field. Importantly, the Lagrange multiplier behaves as an externally applied generalized moment density that, similar to gravity, does not affect the total strain-plus-sprain energy density of material. The Helmholtz free energy of the finite element and its associated stiffness matrix are formulated and implemented in a user’s element of abaqus. The conditions of stationary values of the total free energy of the structure with respect to the nodal degrees-of-freedom yield the set of equilibrium equations of the structure for each loading step. One- and two-dimensional examples of crack growth in fracture specimens are given. It is demonstrated that the simulation results of the three-point bend test are independent of the orientation of a regular square mesh, capture the width variation of the crack band, the damage strain profile across the band, and converge as the finite element mesh is refined.
more »
« less
Sprain energy consequences for damage localization and fracture mechanics
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
- Award ID(s):
- 2029641
- PAR ID:
- 10626892
- Publisher / Repository:
- Proceedings of the National Academy of Sciences
- Date Published:
- Journal Name:
- Proceedings of the National Academy of Sciences
- Volume:
- 121
- Issue:
- 40
- ISSN:
- 0027-8424
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
A theoretical analysis on crack formation and propagation was performed based on the coupling between the electrochemical process, classical elasticity, and fracture mechanics. The chemical potential of oxygen, thus oxygen partial pressure, at the oxygen electrode-electrolyte interface ( ) was investigated as a function of transport properties, electrolyte thickness and operating conditions (e.g., steam concentration, constant current, and constant voltage). Our analysis shows that: a lower ionic area specific resistance (ASR), and a higher electronic ASR ( ) of the oxygen electrode/electrolyte interface are in favor of suppressing crack formation. The thus local pO2, are sensitive towards the operating parameters under galvanostatic or potentiostatic electrolysis. Constant current density electrolysis provides better robustness, especially at a high current density with a high steam content. While constant voltage electrolysis leads to greater variations of Constant current electrolysis, however, is not suitable for an unstable oxygen electrode because can reach a very high value with a gradually increased A crack may only occur under certain conditions whenmore » « less
-
Abstract This paper proposes and investigates the two-grid stabilized lowest equal-order finite element method for the time-independent dual-permeability-Stokes model with the Beavers-Joseph-Saffman-Jones interface conditions. This method is mainly based on the idea of combining the two-grid and the two local Gauss integrals for the dual-permeability-Stokes system. In this technique, we use a difference between a consistent mass matrix and an under-integrated mass matrix for the pressure variable of the dual-permeability-Stokes model using the lowest equal-order finite element quadruples. In the two-grid scheme, the global problem is solved using the standard$$ P_1-P_1-P_1-P_1 $$ finite element approximations only on a coarse grid with grid sizeH. Then, a coarse grid solution is applied on a fine grid of sizehto decouple the interface terms and the mass exchange terms for solving the three independent subproblems such as the Stokes equations, microfracture equations, and the matrix equations on the fine grid. On the other hand, microfracture and matrix equations are decoupled through the mass exchange terms. The weak formulation is reported, and the optimal error estimate is derived for the two-grid schemes. Furthermore, the numerical results validate that the two-grid stabilized lowest equal-order finite element method is effective and has the same accuracy as the coupling scheme when we choose$$ h=H^2 $$ .more » « less
-
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1073/pnas.2410668121
