Imposing equilibrium on measured 3-D stress fields using Hodge decomposition and FFT-based optimization

We present a methodology to impose micromechanical constraints, i.e. stress equilibrium at grain and sub-grain scale, to an arbitrary (non-equilibrated) voxelized stress field obtained, for example, by means of synchrotron X-ray diffraction techniques. The method consists in finding the equilibrated stress field closest (in L-norm sense) to the measured non-equilibrated stress field, via the solution of an optimization problem. The extraction of the divergence-free (equilibrated) part of a general (non-equilibrated) field is performed using the Hodge decomposition of a symmetric matrix field, which is the generalization of the Helmholtz decomposition of a vector field into the sum of an irrotational field and a solenoidal field. The combination of: a) the Euler-Lagrange equations that solve the optimization problem, and b) the Hodge decomposition, gives a differential expression that contains the bi-harmonic operator and two times the curl operator acting on the measured stress field. These high-order derivatives can be efficiently performed in Fourier space. The method is applied to filter the non-equilibrated parts of a synthetic piecewise constant stress fields with a known ground truth, and stress fields in Gum Metal, a beta-Ti-based alloy measured in-situ using Diffraction Contrast Tomography (DCT). In both cases, the largest corrections were obtained near grain boundaries.

[1]  W. Ludwig,et al.  Three-dimensional reconstruction of intragranular strain and orientation in polycrystals by near-field X-ray diffraction , 2020, Current Opinion in Solid State and Materials Science.

[2]  R. Suter,et al.  Voxel-based strain tensors from near-field High Energy Diffraction Microscopy , 2020 .

[3]  Jonathan P. Wright,et al.  Reconstructing intragranular strain fields in polycrystalline materials from scanning 3DXRD data , 2020, Journal of applied crystallography.

[4]  Y. Hirose,et al.  Intragranular three-dimensional stress tensor fields in plastically deformed polycrystals , 2019, Science.

[5]  A. Beaudoin,et al.  Utilizing a novel lattice orientation based stress characterization method to study stress fields of shear bands , 2019, Journal of the Mechanics and Physics of Solids.

[6]  P. Dawson,et al.  Non-Destructive Characterization of Subsurface Residual Stress Fields and Correlation with Microstructural Conditions in a Shot-Peened Inconel Component , 2018, Experimental Mechanics.

[7]  Hugh T. Philipp,et al.  Study of residual stresses in Ti-7Al using theory and experiments , 2017 .

[8]  D. Kochmann,et al.  Predicting the effective response of bulk polycrystalline ferroelectric ceramics via improved spectral phase field methods , 2017 .

[9]  R. Lebensohn,et al.  Instantiation of crystal plasticity simulations for micromechanical modelling with direct input from microstructural data collected at light sources , 2017 .

[10]  H. Proudhon,et al.  Nanox: a miniature mechanical stress rig designed for near-field X-ray diffraction imaging techniques , 2016, Journal of synchrotron radiation.

[11]  W. Ludwig,et al.  Three-dimensional full-field X-ray orientation microscopy , 2016, Scientific Reports.

[12]  Claude Fressengeas,et al.  A numerical spectral approach for solving elasto-static field dislocation and g-disclination mechanics , 2014 .

[13]  P. Dawson,et al.  A computational framework for evaluating residual stress distributions from diffraction-based lattice strain data , 2013 .

[14]  P. Dawson,et al.  Quantifying Three-Dimensional Residual Stress Distributions Using Spatially-Resolved Diffraction Measurements and Finite Element Based Data Reduction , 2013 .

[15]  Robert M. Suter,et al.  Adaptive reconstruction method for three-dimensional orientation imaging , 2013 .

[16]  W. J. Palenstijn,et al.  Advances in X-ray diffraction contrast tomography: Flexibility in the setup geometry and application to multiphase materials , 2013 .

[17]  Matthew P. Miller,et al.  A two-scale methodology for determining the residual stresses in polycrystalline solids using high energy X-ray diffraction data , 2013 .

[18]  Philip Eisenlohr,et al.  An elasto-viscoplastic formulation based on fast Fourier transforms for the prediction of micromechanical fields in polycrystalline materials , 2012 .

[19]  M. Miller,et al.  Far-field high-energy diffraction microscopy: a tool for intergranular orientation and strain analysis , 2011 .

[20]  M. Miller,et al.  High-energy diffraction microscopy at the advanced photon source , 2011 .

[21]  Jonathan P. Wright,et al.  Determining grain resolved stresses in polycrystalline materials using three-dimensional X-ray diffraction , 2010 .

[22]  J. Buffière,et al.  Three-dimensional grain mapping by x-ray diffraction contrast tomography and the use of Friedel pairs in diffraction data analysis. , 2009, The Review of scientific instruments.

[23]  Giuseppe Geymonat,et al.  Hodge decomposition for symmetric matrix fields and the elasticity complex in Lipschitz domains , 2008 .

[24]  Robert M. Suter,et al.  Forward modeling method for microstructure reconstruction using x-ray diffraction microscopy: Single-crystal verification , 2006 .

[25]  B. Ralph.  Three-dimensional X-ray diffraction microscopy (mapping polycrystals and their dynamics) , 2005 .

[26]  Dennis DeTurck,et al.  Vector Calculus and the Topology of Domains in 3-Space , 2002, Am. Math. Mon..

[27]  Weimin Han,et al.  Some mixed finite element methods for biharmonic equation , 2000 .

[28]  Murli M. Gupta,et al.  Direct solution of the biharmonic equation using noncoupled approach , 1979 .

[29]  Steven G. Johnson.  Notes on FFT-based differentiation , 2011 .

[30]  D. Juul Jensen,et al.  Electronic Reprint Applied Crystallography Three-dimensional Maps of Grain Boundaries and the Stress State of Individual Grains in Polycrystals and Powders Applied Crystallography Three-dimensional Maps of Grain Boundaries and the Stress State of Individual Grains in Polycrystals and Powders , 2001 .