skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Phase field theory for fracture at large strains including surface stresses
Phase field theory for fracture is developed at large strains with an emphasis on a correct introduction of surface stresses. This is achieved by multiplying the cohesion and gradient energies by the local ratio of the crack surface areas in the deformed and undeformed configurations and with the gradient energy in terms of the gradient of the order parameter in the reference configuration. This results in an expression for the surface stresses which is consistent with the sharp surface approach. Namely, the structural part of the Cauchy surface stress represents an isotropic biaxial tension, with the magnitude of a force per unit length equal to the surface energy. The surface stresses are a result of the geometric nonlinearities, even when strains are infinitesimal. They make multiple contributions to the Ginzburg-Landau equation for damage evolution, both in the deformed and undeformed configurations. Important connections between material parameters are obtained using an analytical solution for two separating surfaces, as well as an analysis of the stress-strain curves for homogeneous tension for different degradation and interpolation functions. A complete system of equations is presented in the undeformed and deformed configurations. All the phase field parameters are obtained utilizing the existing first principle simulations for the uniaxial tension of Si crystal in the [100] and [111] directions.  more » « less
Award ID(s):
1904830 1943710
PAR ID:
10220744
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
ArXivorg
ISSN:
2331-8422
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Digital image correlation (DIC) is an increasingly popular and effective non-contact method for measuring full-field displacements and strains of deformable bodies under load. Current DIC methods applied to bodies undergoing large displacements and rotations require a large measurement area for both the reference (i.e., undeformed) image and the deformed images. This can limit the resulting resolution of the displacement and strain fields. To address this issue, we propose a two-stage dynamic DIC method capable of measuring displacements and strains under material convection with high resolution. During the first stage, the reference image is assembled from smaller, high-resolution images of the undeformed object obtained using a spatially-fixed or moving frame. Following capture, each sub-image is rigidly translated and rotated into its appropriate place, thereby producing a full, high-resolution image of the reference body. In stage two, images of the loaded and deformed body, again obtained using a small camera frame with high resolution, are aligned with matching regions of the undeformed composite image using BRISK feature detection before performing DIC.We demonstrate the method on a contact problem whereby an elastomeric roller travels along a rigid surface. In doing so, we obtain fine resolution measurements of the state of strain of the region of the roller sidewall in contact with the substrate, even as new material convects through the region of interest. We present these measurements as a series of images and videos capturing strain evolution as the roller transitions from static loads to a fully dynamic steady-state, documenting the effectiveness of the method. 
    more » « less
  2. Abstract Additively manufactured (AM) metallic materials often comprise as-printed dislocation cells inside grains. These dislocation cells can give rise to substantial microscale internal stresses in both initial undeformed and plastically deformed samples, thereby affecting the mechanical properties of AM metallic materials. Here we develop models of microscale internal stresses in AM stainless steel by focusing on their back stress components. Three sources of microscale back stresses are considered, including the printing and deformation-induced back stresses associated with as-printed dislocation cells as well as the deformation-induced back stresses associated with grain boundaries. We use a three-dimensional discrete dislocation dynamics model to demonstrate the manifestation of printing-induced back stresses. We adopt a dislocation pile-up model to evaluate the deformation-induced back stresses associated with as-printed dislocation cells. The extracted back stress relation from the pile-up model is incorporated into a crystal plasticity model that accounts for the other two sources of back stresses as well. The crystal plasticity finite element simulation results agree with the experimentally measured tension-compression asymmetry and macroscopic back stress, the latter of which represents the effective resultant of microscale back stresses of different origins. Our results provide an in-depth understanding of the origins and evolution of microscale internal stresses in AM metallic materials. 
    more » « less
  3. A thermodynamically consistent multiphase phase-field approach for stress and temperature-induced martensitic phase transformation at the nanoscale and under large strains is developed. A total of N independent order parameters are considered for materials with N variants, where one of the order parameters describes A ↔ M transformations and the remaining N − 1 independent order parameters describe the transformations between the variants. A non-contradictory gradient energy is used within the free energy of the system to account for the energies of the interfaces. In addition, a non-contradictory kinetic relationships for the rate of the order parameters versus thermodynamic driving forces is suggested. As a result, a system of consistent coupled Ginzburg-Landau equations for the order parameters are derived. The crystallographic solution for twins within twins is presented for the cubic to tetragonal transformations. A 3D complex twins within twins microstructure is simulated using the developed phase-field approach and a large-strain-based nonlinear finite element method. A comparative study between the crystallographic solution and the simulation result is presented. 
    more » « less
  4. In this work a transformation strain gradient enhancement is introduced into a phenomenological constitutive model for the pseudoelastic behavior of shape memory alloys. The constitutive model is able to capture several unique features of the constitutive response of these materials during the transformation between austenite and martensite during the pseudoelastic response. These features include the asymmetry in the initial transformation stresses in tension versus compression, the asymmetry in the transformation strains in tension and compression, and finally the asymmetry in the hardening behavior in tension and compression. In fact, experiments have shown that untrained NiTi exhibits hardening during its transformation in compression, but softening for tensile loading. It is this softening behavior that motivates the need for the introduction of the transformation strain gradient into the constitutive modeling. Transformation strain gradient effects are introduced via a phase variable that describes the extent of transformation. The free energy of the material then depends on gradients of the phase variable, which introduces a material length scale into the theory. The governing equation for the phase variable is developed from a microforce balance and continuum thermodynamics analysis. The model is implemented in the commercial finite element software Abaqus through user defined subroutines and several numerical simulations are performed to illustrate the model response and lack of numerical mesh-dependency of the results. 
    more » « less
  5. Recent works have shown that in contrast to classical linear elastic fracture mechanics, endowing crack fronts in a brittle Green-elastic solid with Steigmann-Ogden surface elasticity yields a model that predicts bounded stresses and strains at the crack tips for plane-strain problems. However, singularities persist for anti-plane shear (mode-III fracture) under far field loading, even when Steigmann-Ogden surface elasticity is incorporated. This work is motivated by obtaining a model of brittle fracture capable of predicting bounded stresses and strains for all modes of loading. We formulate an exact general theory of a three-dimensional solid containing a boundary surface with strain-gradient surface elasticity. For planar reference surfaces parameterized by flat coordinates, the form of surface elasticity reduces to that introduced by Hilgers and Pipkin, and when the surface energy is independent of the surface covariant derivative of the stretching, the theory reduces to that of Steigmann and Ogden. We discuss material symmetry using Murdoch and Cohen’s extension of Noll’s theory. We present a model small-strain surface energy that incorporates resistance to geodesic distortion, satisfies strong ellipticity, and requires the same material constants found in the Steigmann-Ogden theory. Finally, we derive and apply the linearized theory to mode-III fracture in an infinite plate under far-field loading. We prove that there always exists a unique classical solution to the governing integro-differential equation, and in contrast to using Steigmann-Ogden surface elasticity, our model is consistent with the linearization assumption in predicting finite stresses and strains at the crack tips. 
    more » « less