Lode coordinates

Lode coordinates ${\displaystyle (z,r,\theta )}$ or Haigh–Westergaard coordinates ${\displaystyle (\xi ,\rho ,\theta )}$.^[1] are a set of tensor invariants that span the space of real, symmetric, second-order, 3-dimensional tensors and are isomorphic with respect to principal stress space. This right-handed orthogonal coordinate system is named in honor of the German scientist Dr. Walter Lode because of his seminal paper written in 1926 describing the effect of the middle principal stress on metal plasticity.^[2] Other examples of sets of tensor invariants are the set of principal stresses ${\displaystyle ({\sigma }_1,{\sigma }_2,{\sigma }_3)}$ or the set of kinematic invariants ${\displaystyle (I_1,J_2,J_3)}$. The Lode coordinate system can be described as a cylindrical coordinate system within principal stress space with a coincident origin and the z-axis parallel to the vector ${\displaystyle ({\sigma }_1,{\sigma }_2,{\sigma }_3)=(1,1,1)}$.

The Lode coordinates are most easily computed using the mechanics invariants. These invariants are a mixture of the invariants of the Cauchy stress tensor, ${\displaystyle \mathit{\sigma }}$, and the stress deviator, ${\displaystyle \mathit{s}}$, and are given by^[3]

${\displaystyle I_1=\mathrm{t}\mathrm{r}(\mathit{\sigma })}$

${\displaystyle J_2=\frac{1}{2}[\text{tr}({\mathit{\sigma }}^2)-\frac{1}{3}\text{tr}(\mathit{\sigma }{)}^2]=\frac{1}{2}\mathrm{t}\mathrm{r}(\mathit{s}\cdot \mathit{s})=\frac{1}{2}\Vert \mathit{s}{\Vert }^2}$

${\displaystyle J_3=\mathrm{d}\mathrm{e}\mathrm{t}(\mathit{s})=\frac{1}{3}\mathrm{t}\mathrm{r}(\mathit{s}\cdot \mathit{s}\cdot \mathit{s})}$

which can be written equivalently in Einstein notation

${\displaystyle I_1={\sigma }_{kk}}$

${\displaystyle J_2=\frac{1}{2}[\text{tr}({\mathit{\sigma }}^2)-\frac{1}{3}\text{tr}(\mathit{\sigma }{)}^2]=\frac{1}{2}{s}_{ij}{s}_{ji}=\frac{1}{2}{s}_{ij}{s}_{ij}}$

${\displaystyle J_3=\frac{1}{6}{ϵ}_{ijk}{ϵ}_{pqr}{\sigma }_{ip}{\sigma }_{jq}{\sigma }_{kr}=\frac{1}{3}{s}_{ij}{s}_{jk}{s}_{ki}}$

where ${\displaystyle ϵ}$ is the Levi-Civita symbol (or permutation symbol) and the last two forms for ${\displaystyle J_2}$ are equivalent because ${\displaystyle \mathit{s}}$ is symmetric (${\displaystyle s_{ij}=s_{ji}}$).

The gradients of these invariants^[4] can be calculated by

${\displaystyle \frac{\partial I_1}{\partial \mathit{\sigma }}=\mathit{I}}$

${\displaystyle \frac{\partial J_2}{\partial \mathit{\sigma }}=\mathit{s}=\mathit{\sigma }-\frac{\mathrm{t}\mathrm{r}\left(\mathit{\sigma }\right)}{3}\mathit{I}}$

${\displaystyle \frac{\partial J_3}{\partial \mathit{\sigma }}=\mathit{T}=\mathit{s}\cdot \mathit{s}-\frac{2J_2}{3}\mathit{I}}$

where ${\displaystyle \mathit{I}}$ is the second-order identity tensor and ${\displaystyle \mathit{T}}$ is called the Hill tensor.

The ${\displaystyle z}$-coordinate is found by calculating the magnitude of the orthogonal projection of the stress state onto the hydrostatic axis.

${\displaystyle z={\mathit{E}}_{\mathit{z}}:\mathit{\sigma }=\frac{\mathrm{t}\mathrm{r}(\mathit{\sigma })}{\sqrt{3}}=\frac{I_1}{\sqrt{3}}}$

where

${\displaystyle {\mathit{E}}_{\mathit{z}}=\frac{\mathit{I}}{\Vert \mathit{I}\Vert }=\frac{\mathit{I}}{\sqrt{3}}}$

is the unit normal in the direction of the hydrostatic axis.

The ${\displaystyle r}$-coordinate is found by calculating the magnitude of the stress deviator (the orthogonal projection of the stress state into the deviatoric plane).

${\displaystyle r={\mathit{E}}_{\mathit{r}}:\mathit{\sigma }=\Vert \mathit{s}\Vert =\sqrt{2J_2}}$

where

${\displaystyle {\mathit{E}}_{\mathit{r}}=\frac{\mathit{s}}{\Vert \mathit{s}\Vert }}$

Derivation

The relation that ${\displaystyle r=\sqrt{2J_2}}$ can be found by expanding the relation ${\displaystyle r=\frac{\mathit{s}}{\Vert \mathit{s}\Vert }:\mathit{\sigma }}$ and writing ${\displaystyle \sigma }$ in terms of the isotropic and deviatoric parts while expanding the magnitude of ${\displaystyle \mathit{s}}$ ${\displaystyle r=\frac{\mathit{s}}{\sqrt{\mathit{s}:\mathit{s}}}:(\mathit{s}+z{\mathit{E}}_{\mathit{z}})=\frac{1}{\sqrt{\mathit{s}:\mathit{s}}}(\mathit{s}:\mathit{s}+z\mathit{s}:{\mathit{E}}_{\mathit{z}})}$. Because ${\displaystyle {\mathit{E}}_{\mathit{z}}}$ is isotropic and ${\displaystyle \mathit{s}}$ is deviatoric, their product is zero. Which leaves us with ${\displaystyle r=\frac{\mathit{s}:\mathit{s}}{\sqrt{\mathit{s}:\mathit{s}}}=\sqrt{\mathit{s}:\mathit{s}}}$ Applying the identity ${\displaystyle \mathit{A}:\mathit{B}=\mathrm{t}\mathrm{r}({\mathit{A}}^T\cdot \mathit{B})}$ and using the definition of ${\displaystyle J_2=\frac{1}{2}\mathrm{t}\mathrm{r}(\mathit{s}\cdot \mathit{s})}$ ${\displaystyle r=\sqrt{\mathrm{t}\mathrm{r}(\mathit{s}\cdot \mathit{s})}=\sqrt{2J_2}}$

is a unit tensor in the direction of the radial component.

The Lode angle can be considered, rather loosely, a measure of loading type. The Lode angle varies with respect to the middle eigenvalue of the stress. There are many definitions of Lode angle that each utilize different trigonometric functions: the positive sine^[5] (It is incorrect (false) information that Chakrabarty defines shear stress mode angle with sin(Teta). Chakrabarty adopted Lode parameter with opposite sign than it has in the original Lode paper from 1926, see Page 59, formula (3) in the book «Chakrabarty, Jagabanduhu; 2006, Theory of Plasticity: Third edition, Elsevier, Amsterdam.». Effectively Chakrabarty shear stress mode angle is defined with -sin(Teta) just like it did for the first time Novozhilov in his work from 1951, Novozhilov 183-194, «О СВЯЗИ МЕЖДУ НАПРЯЖЕНИЯМИ И ДЕФОРМАЦИЯМИ В НЕЛИНЕЙНО УПРУГОЙ СРЕДЕ» (On relations between stresses and strains in non-linear elastic media) from 1951. The first work in which arguments and discussion is delivered why it is wise and useful to define shear stress mode angle with sin(Teta) is the Ziolkowski’s work «Parametrization of Cauchy Stress Tensor Treated as Autonomous Object Using Isotropy Angle and Skewness Angle» from 2022, https://et.ippt.pan.pl/index.php/et/article/view/2210. Accepting definition with sin(3Teta) in fact means that pure shears are accepted as comparison reference states, what has serious consequences.), negative sine,^[6] and positive cosine^[7] (here denoted ${\displaystyle {\theta }_s}$, ${\displaystyle {\overline{\theta }}_s}$, and ${\displaystyle {\theta }_c}$, respectively)

${\displaystyle \sin (3{\theta }_s)=-\sin (3{\overline{\theta }}_s)=\cos (3{\theta }_c)=\frac{J_3}{2}{\left(\frac{3}{J_2}\right)}^{3/2}}$

and are related by

${\displaystyle {\theta }_s=\frac{\pi }{6}-{\theta }_c{\theta }_s=-{\overline{\theta }}_s}$

Derivation

The relation between ${\displaystyle {\theta }_c}$ and ${\displaystyle {\theta }_s}$ can be shown by applying a trigonometric identity relating sine and cosine by a shift ${\displaystyle \cos (3{\theta }_c)=\sin (3{\theta }_s)}$ ${\displaystyle \cos (3{\theta }_c)=\sin ((3{\theta }_s-\frac{\pi }{2})+\frac{\pi }{2})}$ ${\displaystyle \cos (3{\theta }_c)=\cos (3{\theta }_s-\frac{\pi }{2})}$. Because cosine is an even function and the range of the inverse cosine is usually ${\displaystyle 0\le {\cos }^{-1}(\theta )\le 1}$ we take the negative possible value for the ${\displaystyle {\theta }_s}$ term, thus ensuring that ${\displaystyle {\theta }_c}$ is positive. ${\displaystyle 3{\theta }_c=\pm (3{\theta }_s-\frac{\pi }{2})}$ ${\displaystyle {\theta }_c=\frac{\pi }{6}-{\theta }_s}$

These definitions are all defined for a range of ${\displaystyle \pi /3}$.

The unit normal in the angular direction which completes the orthonormal basis can be calculated for ${\displaystyle {\theta }_s}$^[8] and ${\displaystyle {\theta }_c}$^[9] using

${\displaystyle {\mathit{E}}_{\mathit{\theta }\mathit{s}}=\frac{\mathit{T}/\Vert \mathit{T}\Vert -\sin (3{\theta }_s){\mathit{E}}_{\mathit{r}}}{\cos (3{\theta }_s)}{\mathit{E}}_{\mathit{\theta }\mathit{c}}=\frac{\mathit{T}/\Vert \mathit{T}\Vert -\cos (3{\theta }_c){\mathit{E}}_{\mathit{r}}}{\sin (3{\theta }_c)}}$.

The meridional profile is a 2D plot of ${\displaystyle (z,r)}$ holding ${\displaystyle \theta }$ constant and is sometimes plotted using scalar multiples of ${\displaystyle (z,r)}$. It is commonly used to demonstrate the pressure dependence of a yield surface or the pressure-shear trajectory of a stress path. Because ${\displaystyle r}$ is non-negative the plot usually omits the negative portion of the ${\displaystyle r}$-axis, but can be included to illustrate effects at opposing Lode angles (usually triaxial extension and triaxial compression).

One of the benefits of plotting the meridional profile with ${\displaystyle (z,r)}$ is that it is a geometrically accurate depiction of the yield surface.^[8] If a non-isomorphic pair is used for the meridional profile then the normal to the yield surface will not appear normal in the meridional profile. Any pair of coordinates that differ from ${\displaystyle (z,r)}$ by constant multiples of equal absolute value are also isomorphic with respect to principal stress space. As an example, pressure ${\displaystyle p=-I1/3}$ and the Von Mises stress ${\displaystyle {\sigma }_v=\sqrt{3J_2}}$ are not an isomorphic coordinate pair and, therefore, distort the yield surface because

${\displaystyle p=-\frac{1}{\sqrt{3}}z}$

${\displaystyle {\sigma }_v=\sqrt{\frac{3}{2}}r}$

and, finally, ${\displaystyle |-1/\sqrt{3}|\ne |\sqrt{3/2}|}$.

The octahedral profile is a 2D plot of ${\displaystyle (r,\theta )}$ holding ${\displaystyle z}$ constant. Plotting the yield surface in the octahedral plane demonstrates the level of Lode angle dependence. The octahedral plane is sometimes referred to as the 'pi plane'^[10] or 'deviatoric plane'.^[11]

The octahedral profile is not necessarily constant for different values of pressure with the notable exceptions of the von Mises yield criterion and the Tresca yield criterion which are constant for all values of pressure.

The term Haigh-Westergaard space is ambiguously used in the literature to mean both the Cartesian principal stress space^[12][13] and the cylindrical Lode coordinate space^[14][15]

• Yield (engineering)
• Plasticity
• Stress
• Henri Tresca
• von Mises stress
• Mohr–Coulomb theory
• Strain
• Strain tensor
• Stress–energy tensor
• Stress concentration
• 3-D elasticity

0.00 (0 votes) Original source: https://en.wikipedia.org/wiki/Lode coordinates. Read more