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
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
- NSF-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
-
-
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
Summary The discrete crack mechanics (DCM) method is a dislocation‐based crack modeling technique where cracks are constructed using Volterra dislocation loops. The method allows for the natural introduction of displacement discontinuities, avoiding numerically expensive techniques. Mesh dependence in existing computational modeling of crack growth is eliminated by utilizing a superposition procedure. The elastic field of cracks in finite bodies is separated into two parts: the infinite‐medium solution of discrete dislocations and an finite element method solution of a correction problem that satisfies external boundary conditions. In the DCM, a crack is represented by a dislocation array with a fixed outer loop determining the crack tip position encompassing additional concentric loops free to expand or contract. Solving for the equilibrium positions of the inner loops gives the crack shape and stress field. The equation of motion governing the crack tip is developed for quasi‐static growth problems. Convergence and accuracy of the DCM method are verified with two‐ and three‐dimensional problems with well‐known solutions. Crack growth is simulated under load and displacement (rotation) control. In the latter case, a semicircular surface crack in a bent prismatic beam is shown to change shape as it propagates inward, stopping as the imposed rotation is accommodated.
-
more » « less
Abstract Strain localization and resulting plasticity and failure play an important role in the evolution of the lithosphere. These phenomena are commonly modeled by Stokes flows with viscoplastic rheologies. The nonlinearities of these rheologies make the numerical solution of the resulting systems challenging, and iterative methods often converge slowly or not at all. Yet accurate solutions are critical for representing the physics. Moreover, for some rheology laws, aspects of solvability are still unknown. We study a basic but representative viscoplastic rheology law. The law involves a yield stress that is independent of the dynamic pressure, referred to as von Mises yield criterion. Two commonly used variants, perfect/ideal and composite viscoplasticity, are compared. We derive both variants from energy minimization principles, and we use this perspective to argue when solutions are unique. We propose a new stress‐velocity Newton solution algorithm that treats the stress as an independent variable during the Newton linearization but requires solution only of Stokes systems that are of the usual velocity‐pressure form. To study different solution algorithms, we implement 2‐D and 3‐D finite element discretizations, and we generate Stokes problems with up to 7 orders of magnitude viscosity contrasts, in which compression or tension results in significant nonlinear localization effects. Comparing the performance of the proposed Newton method with the standard Newton method and the Picard fixed‐point method, we observe a significant reduction in the number of iterations and improved stability with respect to problem nonlinearity, mesh refinement, and the polynomial order of the discretization.
-
more » « less
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 under 0.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.
