An official website of the United States government
In this work, we present and prove results underlying a method which uses functionals derived from the interaction integral to approximate the stress intensity factors along a three-dimensional crack front. We first prove that the functionals possess a pair of important properties. The functionals are well-defined and continuous for square-integrable tensor fields, such as the gradient of a finite element solution. Furthermore, the stress intensity factors are representatives of such functionals in a space of functions over the crack front. Our second result is an error estimate for the numerical stress intensity factors computed via our method. The latter property of the functionals provides a recipe for numerical stress intensity factors; we apply the functionals to the gradient of a finite element approximation for a specific set of crack front variations, and we calculate the stress intensity factors by inverting the mass matrix for those variations.
more »
« less
A method to compute mixed‐mode stress intensity factors for nonplanar cracks in three dimensions
Summary Methods to compute the stress intensity factors along a three‐dimensional (3D) crack front often display a tenuous rate of convergence under mesh refinement or, worse, do not converge, particularly when applied on unstructured meshes. In this work, we propose an alternative formulation of the interaction integral functional and a method to compute stress intensity factors along the crack front which can be shown to converge. The novelty of our method is the decoupling of the two discretizations: the bulk mesh for the finite element solution and the mesh along the crack front for the numerical stress intensity factors, and hence we term it the multiple mesh interaction integral (MMII) method. Through analysis of the convergence of the functional and method, we find scalings of these two mesh sizes to guarantee convergence of the computed stress intensity factors in a variety of norms, including maximum pointwise error and total variation. We demonstrate the MMII on four examples: a semiinfinite straight crack with the asymptotic displacement fields, the same geometry with a nonuniform stress intensity factor along the crack front, a spherical cap crack in a cylinder under tension, and the elliptical crack under far‐field tension and shear.
more »
« less
- Award ID(s):
- 1662452
- PAR ID:
- 10376882
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- International Journal for Numerical Methods in Engineering
- Volume:
- 121
- Issue:
- 19
- ISSN:
- 0029-5981
- Page Range / eLocation ID:
- p. 4292-4328
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Coupon specimens of poled and depoled lead zirconate titanate (PZT) are examined under combined stress wave and electric loading conditions. Mode‐I crack initiation and fracture behavior is examined using ultrahigh‐speed imaging and two‐dimensional digital image correlation. The dynamic critical stress intensity factor () is extracted using measured displacement fields ahead of the impulsively loaded crack tip, and compared between poled and depoled plates that were either under no electric field, positive 0.46 kV/mm electric field, or negative 0.46 kV/mm electric field. Poled specimens had a poling direction and applied electric field direction normal to the crack front. The addition of an electric field resulted in a crack‐enhancing effect, where the dynamic fracture toughness of poled specimens under0.46 kV/mm was almost half that of samples with no electric field. Depoled samples experienced almost no change in dynamic fracture toughness with the addition of an electric field.more » « less
-
Abstract In traditional modeling approaches, earthquakes are often depicted as displacement discontinuities across zero‐thickness surfaces embedded within a linear elastodynamic continuum. This simplification, however, overlooks the intricate nature of natural fault zones and may fail to capture key physical phenomena integral to fault processes. Here, we propose a diffuse interface description for dynamic earthquake rupture modeling to address these limitations and gain deeper insight into fault zones' multifaceted volumetric failure patterns, mechanics, and seismicity. Our model leverages a steady‐state phase‐field, implying time‐independent fault zone geometry, which is defined by the contours of a signed distance function relative to a virtual fault plane. Our approach extends the classical stress glut method, adept at approximating fault‐jump conditions through inelastic alterations to stress components. We remove the sharp discontinuities typically introduced by the stress glut approach via our spatially smooth, mesh‐independent fault representation while maintaining the method's inherent logical simplicity within the well‐established spectral element method framework. We verify our approach using 2D numerical experiments in an open‐source spectral element implementation, examining both a kinematically driven Kostrov‐like crack and spontaneous dynamic rupture in diffuse fault zones. The capabilities of our methodology are showcased through mesh‐independent planar and curved fault zone geometries. Moreover, we highlight that our phase‐field‐based diffuse rupture dynamics models contain fundamental variations within the fault zone. Dynamic stresses intertwined with a volumetrically applied friction law give rise to oblique plastic shear and fault reactivation, markedly impacting rupture front dynamics and seismic wave radiation. Our results encourage future applications of phase‐field‐based earthquake modeling.more » « less
-
Abstract Cracking during sintering is a common problem in powder processing and is usually caused by constraint that prevents the sintering material from shrinking in one or more directions. Different factors influence sintering‐induced cracking, including temperature schedule, packing density, and specimen geometry. Here we use the discrete element method to directly observe the stress distribution and sinter‐cracking behavior in edge notched panels sintered under a uniaxial restraint. This geometry allows an easy comparison with traditional fracture mechanics parameters, facilitating analysis of sinter‐cracking behavior. We find that cracking caused by self‐stress during sintering resembles the growth of creep cracks in fully dense materials. By deriving the constrained densification rate from the appropriate constitutive equations, we discover that linear shrinkage transverse to the loading axis is accelerated by a contribution from the effective Poisson's ratio of a sintering solid. Simulation of different notch geometries and initial relative densities reveals conditions that favor densification and minimize crack growth, alluding to design methods for avoiding cracking in actual sintering processes. We combine the far‐field stress and crack length to compute the net section stress, finding that it characterizes the stress profile between the notches and correlates with the sinter‐crack growth rate, demonstrating its potential to quantitatively describe sinter‐cracking.more » « less
-
Living tissues and some engineering materials contain water. When a wet material loses water, high triaxial tensile stress may build up and cause instability. The mechanism of instability under triaxial tension has attracted great attention, but quantitative study remains an ongoing chal- lenge. Here we develop an experimental method to apply well-controlled triaxial tensile stress and observe osmotic instability in situ. We synthesize a hydrogel in an elastomer tube with strong adhesion between them. The elastomer dissolves minute amount of water, but allows water to diffuse out and places the hydrogel under homogeneous, equal-triaxial, tensile stress. We develop a method to determine the stress as a function of time. The transparent setup enables observation of various types of osmotic instabilities, including cavity nucleation, crack propagation, and surface undulation. Notably, our method enables the measurement of crack speed from ~10−5 m/ s to a limit comparable to the Rayleigh wave speed ~1 m/s. We observe a large jump in crack speed at a critical energy release rate. This work opens opportunities to study the physics of soft materials under high triaxial tension.more » « less
