DDD-XRD Matlab implementation and validation

DDD-XRD approach

DDD-XRD: Matlab implementation and validation

Nicolas Bertin and Wei Cai

January 2018

---

This page provides Matlab functions to compute the displacement gradient associated with discrete dislocation segments based on the non-singular formulation presented in Bertin and Cai, CMS, 2018.

Analytical non-singular displacement gradient formulation

It can be shown that the deformation gradient Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle G_{ml} = u_{m,l}} produced by a dislocation loop Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle C} with Burgers vector Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \vec{b}} can be written as (see Bertin and Cai, CMS, 2018):

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} u_{m,l}(\vec{x}) = &-\frac{1}{8\pi} \oint_C b_m \epsilon_{jlk} R_{,ppj} dx'_k -\frac{1}{8\pi} \oint_C b_i \epsilon_{mik} R_{,ppl} \, dx'_k \\ &-\frac{1}{8\pi(1-\nu)} \oint_C b_i \epsilon_{ijk} R_{,mjl} \, dx'_k \end{align} }

where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle R = \|\vec{x'}-\vec{x}\|} is the norm of the distance vector linking the field point Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \vec{x}} and the coordinate ${\vec {x'}}$ spanning the dislocation line, $\epsilon _{ijk}$ is the permutation tensor, and $\nu$ is Poisson's ratio of the medium. $R_{,ijk}=\partial ^{3}R/\partial x_{i}\partial x_{j}\partial x_{k}$ denotes the third derivative of the radius vector wrt. the field coordinate.

Following the isotropic Burgers distribution proposed in Cai et al., JMPS, 2006, the singularity in the above expression can be eliminated by employing the modified radius vector $R_{a}$ defined as:

$R_{a}={\sqrt {R^{2}+a^{2}}}={\sqrt {R_{i}R_{i}+a^{2}}}={\sqrt {(x'_{i}-x_{i})(x'_{i}-x_{i})+a^{2}}}$

where $a$ denotes the dislocation core radius. With this, the displacement gradient produced by a straight dislocation segment with end points ${\vec {x}}_{a}$ and ${\vec {x}}_{b}$ can be analytically expressed as:

${\begin{aligned}u_{m,l}({\vec {x}})=&-{\frac {1}{8\pi }}b_{m}\epsilon _{jlk}A_{jk}({\vec {x}})-{\frac {1}{8\pi }}b_{i}\epsilon _{mik}A_{lk}({\vec {x}})\\&-{\frac {1}{8\pi (1-\nu )}}b_{i}\epsilon _{ijk}B_{jklm}({\vec {x}})\end{aligned}}$

where line integrals $A_{jk}({\vec {x}})$ and $B_{jklm}({\vec {x}})$ are given by:

$A_{jk}({\vec {x}})=t_{k}\left[3a^{2}d_{j}J_{05}+2d_{j}J_{03}+3a^{2}t_{j}J_{15}+2t_{j}J_{13}\right]$

${\begin{aligned}B_{jklm}({\vec {x}})=&t_{k}\left[(\delta _{mj}d_{l}+\delta _{jl}d_{m}+\delta _{lm}d_{j})J_{03}+(\delta _{mj}t_{l}+\delta _{jl}t_{m}+\delta _{lm}t_{j})J_{13}-3(d_{m}d_{j}d_{l})J_{05}\right.\\&\left.-3(d_{m}d_{j}t_{l}+d_{m}t_{j}d_{l}+t_{m}d_{j}d_{l})J_{15}-3(d_{m}t_{j}t_{l}+t_{m}d_{j}t_{l}+t_{m}t_{j}d_{l})J_{25}-3(t_{m}t_{j}t_{l})J_{35}\right]\end{aligned}}$

where ${\vec {t}}$ denotes the unit dislocation line tangent, and ${\vec {d}}={\vec {x}}_{0}-{\vec {x}}$ is the vector linking field point ${\vec {x}}$ to its orthogonal projection ${\vec {x}}_{0}$ on the dislocation line.

When adopting the following segment parametric representation

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \vec{x'} = \vec{x}_0 + s\vec{t}, \; s \in (s_1,s_2) }

an analytical expression for the line integrals Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle J_{ij}} is obtained as follows:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \begin{align} &J_{03} = \int_{s_1}^{s_2} \frac{1}{R_a^3} ds = \left. \frac{s}{( \vec{d} \cdot \vec{d} + a^2 ) R_a } \right|_{s_1}^{s_2} \\ &J_{13} = \int_{s_1}^{s_2} \frac{s}{R_a^3} ds = \left. -\frac{1}{R_a} \right|_{s_1}^{s_2} \\ &J_{05} = \int_{s_1}^{s_2} \frac{1}{R_a^5} ds = \left. \frac{2s^3}{3( \vec{d} \cdot \vec{d} + a^2 )^2 R_a^3 } \right|_{s_1}^{s_2} + \left. \frac{s}{( \vec{d} \cdot \vec{d} + a^2 ) R_a^3 } \right|_{s_1}^{s_2} \\ &J_{15} = \int_{s_1}^{s_2} \frac{s}{R_a^5} ds = \left. -\frac{1}{3 R_a^3} \right|_{s_1}^{s_2} \\ &J_{25} = \int_{s_1}^{s_2} \frac{s^2}{R_a^5} ds = \left. \frac{s^3}{3( \vec{d} \cdot \vec{d} + a^2 ) R_a^3 } \right|_{s_1}^{s_2} \\ &J_{35} = \int_{s_1}^{s_2} \frac{s^3}{R_a^5} ds = \left. -\frac{2(\vec{d} \cdot \vec{d} + a^2)}{3 R_a^3} \right|_{s_1}^{s_2} -\left. \frac{s^2}{R_a^3 } \right|_{s_1}^{s_2} \end{align} }

Matlab implementation

Code

The computation of the non-singular displacement gradient based on the expressions provided in the above (Bertin and Cai, CMS, 2018) is implemented in a Matlab code available here:

ddd_xrd_matlab.tar

Function displacement_gradient_seg.m readily implements the displacement gradient for a straight dislocation segment using the expressions presented in the above. An optimized version (faster) is provided in displacement_gradient_seg_opt.m. The calculation of the displacement gradient field for a set of dislocation segments can be performed using function displacement_gradient_structure.m. Note that this implementation uses the DDLab data structure to represent the dislocations.

Validation: triangular loop test case

The code includes a test case (test_triangular_loop.m) that (i) provides an example of how to use the calculation functions, and (ii) performs calculations that validate the non-singular displacement gradient formulation.

This test case does the following:

• generates a random triangular dislocation loop
• computes the displacement gradient field along a line using the non-singular expression provided in Bertin and Cai, CMS, 2018.
• compares it with the displacement gradient field obtained by numerically differentiating the displacement field
• compares the stress obtained from the displacement gradient field with the non-singular stress expression provided in Cai et al., JMPS, 2006