Generation of a large network of hydraulic cracks is of key importance not only for the success of fracking of shale but also for the recent scheme of sequestration of CO2 in deep formations of basalt and peridotite, which are mafic and ultramafic rocks that combine chemically with CO2. In numerical simulation of the creation of a fracture network in porous rock, an important goal is to enhance the rock permeability. The objective of this article is to calculate the effect of osmotic pressure gradients caused by gradients of concentration of the ions of Ca, Mg, Na, etc. on the effective permeability of the rock. The basic differential equations are formulated, and their explicit solutions for appropriate initial and boundary conditions are obtained under certain plausible simplifications. The main result is explicit approximate formulas for the critical time before which no water permeation through a test specimen can be observed. Depending on various parameters, this time can be unacceptably long, which is manifested as a zero water outflow. The solution may also explain the unreasonably small permeability values reported for some shales.
Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher.
Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?
Some links on this page may take you to nonfederal websites. Their policies may differ from this site.

Abstract Free, publiclyaccessible full text available December 1, 2024 
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 workconjugate to “sprain.” The partial derivatives of the associated sprain energy density yielded in the preceeding study, sets of curvature resisting selfequilibrated 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 thirdorder 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 strainplussprain 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 degreesoffreedom yield the set of equilibrium equations of the structure for each loading step. One and twodimensional examples of crack growth in fracture specimens are given. It is demonstrated that the simulation results of the threepoint 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.
Free, publiclyaccessible full text available March 1, 2025 
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 oneelement 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 zigzag damage band. Presented is a generalization that overcomes all three, by enforcing a variable multielement 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 thirdorder 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 straingradient 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 selfequilibrated inplane nodal forces, which follow as partial derivatives of Φ(l0η). The force triplets enforce a variable multielement crack band width. The damage distribution across the fracture process zone is nonuniform 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.

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.

Free, publiclyaccessible full text available April 1, 2025

Motivated by the extraordinary strength of nacre, which exceeds the strength of its fragile constituents by an order of magnitude, the fishnet statistics became in 2017 the only analytically solvable probabilistic model of structural strength other than the weakestlink and fiberbundle models. These two models lead, respectively, to the Weibull and Gaussian (or normal) distributions at the largesize limit, which are hardly distinguishable in the central range of failure probability. But they differ enormously at the failure probability level of 10−6 , considered as the maximum tolerable for engineering structures. Under the assumption that no more than three fishnet links fail prior to the peak load, the preceding studies led to exact solutions intermediate between Weibull and Gaussian distributions. Here massive Monte Carlo simulations are used to show that these exact solutions do not apply for fishnets with more than about 500 links. The simulations show that, as the number of links becomes larger, the likelihood of having more than three failed links up to the peak load is no longer negligible and becomes large for fishnets with many thousands of links. A differential equation is derived for the probability distribution of nottoolarge fishnets, characterized by the size effect, the mean and the coefficient of variation. Although the largesize asymptotic distribution is beyond the reach of the Monte Carlo simulations, it can by illuminated by approximating the largescale fishnet as a continuum with a crack or a circular hole. For the former, instability is proven via complex variables, and for the latter via a known elasticity solution for a hole in a continuum under antiplane shear. The fact that rows or enclaves of link failures acting as cracks or holes can form in the largescale continuum at many random locations necessarily leads to the Weibull distribution of the large fishnet, given that these cracks or holes become unstable as soon they reach a certain critical size. The Weibull modulus of this continuum is estimated to be more than triple that of the central range of small fishnets. The new model is expected to allow spinoffs for printed materials with octet architecture maximizing the strength–weight ratio.more » « lessFree, publiclyaccessible full text available January 1, 2025

Free, publiclyaccessible full text available December 1, 2024

Free, publiclyaccessible full text available December 1, 2024

Free, publiclyaccessible full text available September 1, 2024

Free, publiclyaccessible full text available August 1, 2024