Voigt notation

In mathematics, Voigt notation or Voigt form in multilinear algebra is a way to represent a symmetric tensor by reducing its order.^[1] There are a few variants and associated names for this idea: Mandel notation, Mandel–Voigt notation and Nye notation are others found. Kelvin notation is a revival by Helbig^[2] of old ideas of Lord Kelvin. The differences here lie in certain weights attached to the selected entries of the tensor. Nomenclature may vary according to what is traditional in the field of application. The notation is named after physicists Woldemar Voigt^[1] & John Nye (scientist).

For example, a 2×2 symmetric tensor X has only three distinct elements, the two on the diagonal and the other being off-diagonal. Thus its rank can be reduced by expressressing it as a vector without loss of information:

$${\displaystyle X={\begin{bmatrix}x_{11}&x_{12}\\x_{12}&x_{22}\end{bmatrix}}={\begin{bmatrix}x_{11}\\x_{22}\\x_{12}\end{bmatrix}}.}$$

Voigt notation is used in materials science to simplify the representation of the rank-2 stress and strain tensors, and fourth-rank stiffness and compliance tensors.

The 3×3 stress and strain tensors in their full forms can be written as:

${\displaystyle {\boldsymbol {\sigma }}={\begin{bmatrix}\sigma _{11}&\sigma _{12}&\sigma _{13}\\\sigma _{21}&\sigma _{22}&\sigma _{23}\\\sigma _{31}&\sigma _{32}&\sigma _{33}\end{bmatrix}}\quad }$ and ${\displaystyle \quad {\boldsymbol {\varepsilon }}={\begin{bmatrix}\varepsilon _{11}&\varepsilon _{12}&\varepsilon _{13}\\\varepsilon _{21}&\varepsilon _{22}&\varepsilon _{23}\\\varepsilon _{31}&\varepsilon _{32}&\varepsilon _{33}\end{bmatrix}}}$.

Voigt notation then utilises the symmetry of these matrices (${\displaystyle \sigma _{12}=\sigma _{21}}$ and so on) to express them instead as a 6×1 vector:

${\displaystyle {\underline {\sigma }}={\begin{bmatrix}\sigma _{1}\\\sigma _{2}\\\sigma _{3}\\\sigma _{4}\\\sigma _{5}\\\sigma _{6}\end{bmatrix}}:={\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{33}\\\sigma _{23}\\\sigma _{13}\\\sigma _{12}\end{bmatrix}}\quad }$ and ${\displaystyle \quad {\underline {\varepsilon }}={\begin{bmatrix}\varepsilon _{1}\\\varepsilon _{2}\\\varepsilon _{3}\\\varepsilon _{4}\\\varepsilon _{5}\\\varepsilon _{6}\end{bmatrix}}:={\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{33}\\\gamma _{23}\\\gamma _{13}\\\gamma _{12}\end{bmatrix}}}$

where ${\displaystyle \gamma _{12}=2\varepsilon _{12}}$, ${\displaystyle \gamma _{23}=2\varepsilon _{23}}$, and ${\displaystyle \gamma _{13}=2\varepsilon _{13}}$ are the engineering shear strains.

The benefit of using different representations for stress and strain is that the scalar invariance $${\displaystyle {\boldsymbol {\sigma }}\cdot {\boldsymbol {\varepsilon }}=\sigma _{ij}\varepsilon _{ij}={\underline {\sigma }}\cdot {\underline {\varepsilon }}}$$ is preserved.

This notation now allows the three-dimensional symmetric fourth-order stiffness, ${\displaystyle C}$, and compliance, ${\displaystyle S}$, tensors to be reduced to 6×6 matrices:

$${\displaystyle C_{ijkl}\Rightarrow C_{\alpha \beta }={\begin{bmatrix}C_{11}&C_{12}&C_{13}&C_{14}&C_{15}&C_{16}\\C_{12}&C_{22}&C_{23}&C_{24}&C_{25}&C_{26}\\C_{13}&C_{23}&C_{33}&C_{34}&C_{35}&C_{36}\\C_{14}&C_{24}&C_{34}&C_{44}&C_{45}&C_{46}\\C_{15}&C_{25}&C_{35}&C_{45}&C_{55}&C_{56}\\C_{16}&C_{26}&C_{36}&C_{46}&C_{56}&C_{66}\end{bmatrix}}.}$$

A simple mnemonic rule for memorizing Voigt notation is as follows:

• Write down the second order tensor in matrix form (in the example, the stress tensor)
• Strike out the diagonal
• Continue on the third column
• Go back to the first element along the first row.

Voigt indexes are numbered consecutively from the starting point to the end (in the example, the numbers in blue).

The diagram below also shows the order of the indices: $${\displaystyle {\begin{matrix}ij&=\\\Downarrow &\\\alpha &=\end{matrix}}{\begin{matrix}11&22&33&23,32&13,31&12,21\\\Downarrow &\Downarrow &\Downarrow &\Downarrow &\Downarrow &\Downarrow &\\1&2&3&4&5&6\end{matrix}}}$$

For a symmetric tensor of second rank $${\displaystyle {\boldsymbol {\sigma }}={\begin{bmatrix}\sigma _{11}&\sigma _{12}&\sigma _{13}\\\sigma _{21}&\sigma _{22}&\sigma _{23}\\\sigma _{31}&\sigma _{32}&\sigma _{33}\end{bmatrix}}}$$ only six components are distinct, the three on the diagonal and the others being off-diagonal. Thus it can be expressed, in Mandel notation,^[3] as the vector $${\displaystyle {\tilde {\sigma }}^{M}=\langle \sigma _{11},\sigma _{22},\sigma _{33},{\sqrt {2}}\sigma _{23},{\sqrt {2}}\sigma _{13},{\sqrt {2}}\sigma _{12}\rangle .}$$

The main advantage of Mandel notation is to allow the use of the same conventional operations used with vectors, for example: $${\displaystyle {\tilde {\sigma }}:{\tilde {\sigma }}={\tilde {\sigma }}^{M}\cdot {\tilde {\sigma }}^{M}=\sigma _{11}^{2}+\sigma _{22}^{2}+\sigma _{33}^{2}+2\sigma _{23}^{2}+2\sigma _{13}^{2}+2\sigma _{12}^{2}.}$$

A symmetric tensor of rank four satisfying ${\displaystyle D_{ijkl}=D_{jikl}}$ and ${\displaystyle D_{ijkl}=D_{ijlk}}$ has 81 components in three-dimensional space, but only 36 components are distinct. Thus, in Mandel notation, it can be expressed as $${\displaystyle {\tilde {D}}^{M}={\begin{pmatrix}D_{1111}&D_{1122}&D_{1133}&{\sqrt {2}}D_{1123}&{\sqrt {2}}D_{1113}&{\sqrt {2}}D_{1112}\\D_{2211}&D_{2222}&D_{2233}&{\sqrt {2}}D_{2223}&{\sqrt {2}}D_{2213}&{\sqrt {2}}D_{2212}\\D_{3311}&D_{3322}&D_{3333}&{\sqrt {2}}D_{3323}&{\sqrt {2}}D_{3313}&{\sqrt {2}}D_{3312}\\{\sqrt {2}}D_{2311}&{\sqrt {2}}D_{2322}&{\sqrt {2}}D_{2333}&2D_{2323}&2D_{2313}&2D_{2312}\\{\sqrt {2}}D_{1311}&{\sqrt {2}}D_{1322}&{\sqrt {2}}D_{1333}&2D_{1323}&2D_{1313}&2D_{1312}\\{\sqrt {2}}D_{1211}&{\sqrt {2}}D_{1222}&{\sqrt {2}}D_{1233}&2D_{1223}&2D_{1213}&2D_{1212}\\\end{pmatrix}}.}$$

It is useful, for example, in calculations involving constitutive models to simulate materials, such as the generalized Hooke's law, as well as finite element analysis,^[4] and Diffusion MRI.^[5]

Hooke's law has a symmetric fourth-order stiffness tensor with 81 components (3×3×3×3), but because the application of such a rank-4 tensor to a symmetric rank-2 tensor must yield another symmetric rank-2 tensor, not all of the 81 elements are independent. Voigt notation enables such a rank-4 tensor to be represented by a 6×6 matrix. However, Voigt's form does not preserve the sum of the squares, which in the case of Hooke's law has geometric significance. This explains why weights are introduced (to make the mapping an isometry).

A discussion of invariance of Voigt's notation and Mandel's notation can be found in Helnwein (2001).^[6]

• Vectorization (mathematics)
• Hooke's law
• Linear_elasticity#Anisotropic_homogeneous_media