Lode coordinates explained
Lode coordinates
or
Haigh–Westergaard coordinates
.
[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
(\sigma1,\sigma2,\sigma3)
or the set of kinematic invariants
. 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
(\sigma1,\sigma2,\sigma3)=(1,1,1)
.
Mechanics invariants
The Lode coordinates are most easily computed using the mechanics invariants. These invariants are a mixture of the invariants of the Cauchy stress tensor,
, and the stress deviator,
, and are given by
[3] I1=tr(\boldsymbol{\sigma})
J2=
\left[tr(\boldsymbol{\sigma}2)-
tr(\boldsymbol{\sigma})2\right]=
tr\left(\boldsymbol{s} ⋅ \boldsymbol{s}\right)=
\lVert\boldsymbol{s}\rVert2
J3=det(\boldsymbol{s})=
tr\left(\boldsymbol{s} ⋅ \boldsymbol{s} ⋅ \boldsymbol{s}\right)
which can be written equivalently in Einstein notation
J2=
\left[tr(\boldsymbol{\sigma}2)-
tr(\boldsymbol{\sigma})2\right]=
sijsji=
sijsij
J3=
\epsilonijk\epsilonpqr\sigmaip\sigmajq\sigmakr=
sijsjkski
where
is the
Levi-Civita symbol (or permutation symbol) and the last two forms for
are equivalent because
is symmetric (
).
The gradients of these invariants[4] can be calculated by
| \partialI1 |
\partial\boldsymbol{\sigma |
} = \boldsymbol
| \partialJ2 |
\partial\boldsymbol{\sigma |
} = \boldsymbol = \boldsymbol - \frac\boldsymbol
| \partialJ3 |
\partial\boldsymbol{\sigma |
} = \boldsymbol = \boldsymbol\cdot\boldsymbol - \frac\boldsymbol
where
is the second-order identity tensor and
is called the Hill tensor.
Axial coordinate
The
-coordinate is found by calculating the magnitude of the orthogonal projection of the stress state onto the hydrostatic axis.
z=\boldsymbol{Ez}\colon\boldsymbol{\sigma}=
| tr(\boldsymbol{\sigma |
)}{\sqrt{3}} |
=
}
where
} = \frac
is the unit normal in the direction of the hydrostatic axis.
Radial coordinate
The
-coordinate is found by calculating the magnitude of the stress deviator (the orthogonal projection of the stress state into the deviatoric plane).
r=\boldsymbol{Er}\colon\boldsymbol{\sigma}=\lVert\boldsymbol{s}\rVert=\sqrt{2J2}
where
}
is a unit tensor in the direction of the radial component.
Lode angle – angular coordinate
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] negative sine,[6] and positive cosine[7] (here denoted
,
, and
, respectively)
\sin(3\thetas)=-\sin(3\bar{\theta}s)=\cos(3\thetac)=
\right)3/2
and are related by
\thetas=
-\thetac
\thetas=-\bar{\theta}s
These definitions are all defined for a range of
.
| Stress State \sigma1\geq\sigma2\geq\sigma3
|
|
|
|
---|
range | |
|
|
|
Triaxial Compression (TXC) | \sigma1\geq\sigma2=\sigma3
|
|
|
|
Shear (SHR) | \sigma2=(\sigma1+\sigma3)/2
|
|
|
|
Triaxial Extension (TXE) | \sigma1=\sigma2\geq\sigma3
|
|
|
|
|
The unit normal in the angular direction which completes the orthonormal basis can be calculated for
[8] and
[9] using
} = \frac \qquad \boldsymbol = \frac.
Meridional profile
The meridional profile is a 2D plot of
holding
constant and is sometimes plotted using scalar multiples of
. It is commonly used to demonstrate the pressure dependence of a
yield surface or the pressure-shear trajectory of a stress path. Because
is non-negative the plot usually omits the negative portion of the
-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
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
by constant multiples of equal absolute value are also isomorphic with respect to principal stress space. As an example, pressure
and the
Von Mises stress
are not an isomorphic coordinate pair and, therefore, distort the yield surface because
}z
} r
and, finally,
|-1/\sqrt{3}| ≠ |\sqrt{3/2}|
.
Octahedral profile
The octahedral profile is a 2D plot of
holding
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.
A note on terminology
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]
See also
References
- Menetrey, P.H., Willam, K.J., 1995, Triaxial Failure Criterion for Concrete and Its Generalization, ACI Structural Journal
- Lode, W. (1926). Versuche über den Einfuss der mittleren Hauptspannung auf das Fliessen der Metalle Eisen Kupfer und Nickel. Zeitung Phys., vol. 36, pp. 913–939.
- Asaro, R.J., Lubarda, V.A., 2006, Mechanics of Solids and Materials, Cambridge University Press
- Brannon, R.M., 2009, KAYENTA: Theory and User's Guide, Sandia National Laboratories, Albuquerque, New Mexico.
- Chakrabarty, J., 2006, Theory of Plasticity: Third edition, Elsevier, Amsterdam.
- de Souza Neto, E.A., Peric, D., Owen, D.R.J., 2008, Computational Methods for Plasticity, Wiley
- Han, D.J., Chen, W.F., 1985, A Nonuniform Hardening Plasticity Model for Concrete Materials, Mechanics of Materials
- Brannon, R.M., 2007, Elements of Phenomenological Plasticity: Geometrical Insight, Computational Algorithms, and Topics in Shock Physics, Shock Wave Science and Technology Reference Library: Solids I, Springer-New York
- Bigoni, D., Piccolroaz, A., 2004, Yield criteria for quasibrittle and frictional materials, Int. J. Solids Struct.
- Lubliner, J., 1990, Plasticity Theory, Pearson Education
- Chaboche, J.L., 2008, A review of some plasticity and viscoplasticity theories, Int. J. Plasticity
- Mouazen, A.M., Nemenyi, M., 1998, A review of the finite element modelling techniques of soil tillage, Mathematics and Computers in Simulation
- Keryvin, V., 2008, Indentation as a probe for pressure sensitivity of metallic glasses, J. Phys.: Condens. Matter
- Cervenka, J., Papanikolaou, V.K., 2008, Three dimensional combined fracture-plastic material model for concrete, Int. J. of Plasticity
- Piccolroaz, A., Bigoni, D., 2009, Yield criteria for quasibrittle and frictional materials: A generalization to surfaces with corners, Int. J. of Solids and Struc.