A Brief Introduction to Tensors and their properties

1. BASIC PROPERTIES OF TENSORS

1.1 Examples of Tensors

The gradient of a vector field is a good example of a second-order tensor.  Visualize a vector field: at every point in space, the field has a vector value $u\left({x}_{1},{x}_{2},{x}_{3}\right)$.  Let $G=\nabla u\text{\hspace{0.17em}}$ represent the gradient of u.  By definition, G enables you to calculate the change in u when you move from a point x in space to a nearby point at $x+dx$:

$du=G\cdot dx$

G is a second order tensor.  From this example, we see that when you multiply a vector by a tensor, the result is another vector.

This is a general property of all second order tensors.  A tensor is a linear mapping of a vector onto another vector.  Two examples, together with the vectors they operate on, are:

The stress tensor

$t=\text{\hspace{0.17em}}n\cdot \sigma$

where n is a unit vector normal to a surface, $\sigma$ is the stress tensor and t is the traction vector acting on the surface.

$dw=F\cdot dx$

where dx is an infinitesimal line element in an undeformed solid, and dw is the vector representing the deformed line element.

1.2 Matrix representation of a tensor

To evaluate and manipulate tensors, we express them as components in a basis, just as for vectors.  We can use the displacement gradient to illustrate how this is done.  Let $u\left({x}_{1},{x}_{2},{x}_{3}\right)$ be a vector field, and let $G=\nabla u$ represent the gradient of u.  Recall the definition of G

$du=G\cdot dx$

Now, let $\left\{{e}_{1},{e}_{2},{e}_{3}\right\}$ be a Cartesian basis, and express both du and dx as components.  Then, calculate the components of du in terms of dx using the usual rules of calculus

$\begin{array}{l}d{u}_{1}=\frac{\partial {u}_{1}}{\partial {x}_{1}}d{x}_{1}+\frac{\partial {u}_{1}}{\partial {x}_{2}}d{x}_{2}+\frac{\partial {u}_{1}}{\partial {x}_{3}}d{x}_{3}\\ d{u}_{2}=\frac{\partial {u}_{2}}{\partial {x}_{1}}d{x}_{1}+\frac{\partial {u}_{2}}{\partial {x}_{2}}d{x}_{2}+\frac{\partial {u}_{2}}{\partial {x}_{3}}d{x}_{3}\\ d{u}_{3}=\frac{\partial {u}_{3}}{\partial {x}_{1}}d{x}_{1}+\frac{\partial {u}_{3}}{\partial {x}_{2}}d{x}_{2}+\frac{\partial {u}_{3}}{\partial {x}_{3}}d{x}_{3}\end{array}$

We could represent this as a matrix product

$\left[\begin{array}{c}d{u}_{1}\\ d{u}_{2}\\ d{u}_{3}\end{array}\right]=\left[\begin{array}{ccc}\frac{\partial {u}_{1}}{\partial {x}_{1}}& \frac{\partial {u}_{1}}{\partial {x}_{2}}& \frac{\partial {u}_{1}}{\partial {x}_{3}}\\ \frac{\partial {u}_{2}}{\partial {x}_{1}}& \frac{\partial {u}_{2}}{\partial {x}_{2}}& \frac{\partial {u}_{2}}{\partial {x}_{3}}\\ \frac{\partial {u}_{3}}{\partial {x}_{1}}& \frac{\partial {u}_{3}}{\partial {x}_{2}}& \frac{\partial {u}_{3}}{\partial {x}_{3}}\end{array}\right]\left[\begin{array}{c}d{x}_{1}\\ d{x}_{2}\\ d{x}_{3}\end{array}\right]$

Alternatively, using index notation

$d{u}_{i}=\frac{\partial {u}_{i}}{\partial {x}_{j}}d{x}_{j}$

From this example we see that G can be represented as a $3×3$ matrix.  The elements of the matrix are known as the components of G in the basis $\left\{{e}_{1},{e}_{2},{e}_{3}\right\}$.  All second order tensors can be represented in this form.  For example, a general second order tensor S could be written as

$S\equiv \left[\begin{array}{ccc}{S}_{11}& {S}_{12}& {S}_{13}\\ {S}_{21}& {S}_{22}& {S}_{23}\\ {S}_{31}& {S}_{32}& {S}_{33}\end{array}\right]$

You have probably already seen the matrix representation of stress and strain components in introductory courses.

Since S can be represented as a matrix, all operations that can be performed on a $3×3$ matrix can also be performed on S.  Examples include sums and products, the transpose, inverse, and determinant.  One can also compute eigenvalues and eigenvectors for tensors, and thus define the log of a tensor, the square root of a tensor, etc.  These tensor operations are summarized below.

Note that the numbers ${S}_{11}$, ${S}_{12}$, … ${S}_{33}$ depend on the basis $\left\{{e}_{1},{e}_{2},{e}_{3}\right\}$, just as the components of a vector depend on the basis used to represent the vector.  However, just as the magnitude and direction of a vector are independent of the basis, so the properties of a tensor are independent of the basis.  That is to say, if S is a tensor and u is a vector, then the vector

$v=S\cdot u$

has the same magnitude and direction, irrespective of the basis used to represent u, v, and S.

1.3 The difference between a matrix and a tensor

If a tensor is a matrix, why is a matrix not the same thing as a tensor?  Well, although you can multiply the three components of a vector u by any $3×3$ matrix,

$\left[\begin{array}{c}{b}_{1}\\ {b}_{2}\\ {b}_{3}\end{array}\right]=\left[\begin{array}{ccc}{a}_{11}& {a}_{12}& {a}_{13}\\ {a}_{21}& {a}_{22}& {a}_{23}\\ {a}_{31}& {a}_{32}& {a}_{33}\end{array}\right]\left[\begin{array}{c}{u}_{1}\\ {u}_{2}\\ {u}_{3}\end{array}\right]$

the resulting three numbers $\left({b}_{1},{b}_{2},{b}_{3}\right)$ may or may not represent the components of a vector.  If they are the components of a vector, then the matrix represents the components of a tensor A, if not, then the matrix is just an ordinary old matrix.

To check whether $\left({b}_{1},{b}_{2},{b}_{3}\right)$ are the components of a vector, you need to check how $\left({b}_{1},{b}_{2},{b}_{3}\right)$ change due to a change of basis.  That is to say, choose a new basis, calculate the new components of u in this basis, and calculate the new matrix in this basis (the new elements of the matrix will depend on how the matrix was defined.  The elements may or may not change $–$ if they don’t, then the matrix cannot be the components of a tensor).  Then, evaluate the matrix product to find a new left hand side, say $\left({\beta }_{1},{\beta }_{2},{\beta }_{3}\right)$.  If  $\left({\beta }_{1},{\beta }_{2},{\beta }_{3}\right)$ are related to $\left({b}_{1},{b}_{2},{b}_{3}\right)$ by the same transformation that was used to calculate the new components of u, then $\left({b}_{1},{b}_{2},{b}_{3}\right)$ are the components of a vector, and, therefore, the matrix represents the components of a tensor.

1.4 Formal definition

Tensors are rather more general objects than the preceding discussion suggests.   There are various ways to define a tensor formally.  One way is the following:

A tensor is a linear vector valued function defined on the set of all vectors

More specifically, let $S\left(v\right)$ denote a tensor operating on a vector.  Linearity then requires that, for all vectors $v,w$ and scalars $\alpha$

·         $S\left(v+w\right)=S\left(v\right)+S\left(w\right)$

·         $S\left(\alpha v\right)=\alpha S\left(v\right)$

Alternatively, one can define tensors as sets of numbers that transform in a particular way under a change of coordinate system.  In this case we suppose that n dimensional space can be parameterized by a set of n  real numbers ${x}_{i}$.   We could change coordinate system by introducing a second set of real numbers ${{x}^{\prime }}_{i}\left({x}_{k}\right)$ which are invertible functions of ${x}_{i}$.   Tensors can then be defined as sets of real numbers that transform in a particular way under this change in coordinate system.  For example

·         A tensor of zeroth rank is a scalar that is independent of the coordinate system.

·         A covariant tensor of rank 1 is a vector that transforms as ${{v}^{\prime }}_{i}=\frac{\partial {x}_{j}}{\partial {{x}^{\prime }}_{i}}{v}_{j}$

·         A contravariant tensor of rank 1 is a vector that transforms as ${{v}^{\prime }}^{i}=\frac{\partial {{x}^{\prime }}_{i}}{\partial {x}_{j}}{v}^{j}$

·         A covariant tensor of rank 2 transforms as ${{S}^{\prime }}_{ij}=\frac{\partial {x}_{i}}{\partial {{x}^{\prime }}_{k}}\frac{\partial {x}_{j}}{\partial {{x}^{\prime }}_{l}}{S}_{kl}$

·         A contravariant tensor of rank 2 transforms as ${{S}^{\prime }}^{kl}=\frac{\partial {{x}^{\prime }}_{i}}{\partial {x}_{k}}\frac{\partial {{x}^{\prime }}_{i}}{\partial {x}_{l}}{S}^{ij}$

·         A mixed tensor of rank 2 transforms as ${{S}^{\prime }}_{j}^{i}=\frac{\partial {{x}^{\prime }}_{k}}{\partial {x}_{i}}\frac{\partial {x}_{j}}{\partial {{x}^{\prime }}_{l}}{S}_{l}^{k}$

Higher rank tensors can be defined in similar ways.  In solid and fluid mechanics we nearly always use Cartesian tensors, (i.e. we work with the components of tensors in a Cartesian coordinate system) and this level of generality is not needed (and is rather mysterious).  We might occasionally use a curvilinear coordinate system, in which we do express tensors in terms of covariant or contravariant components $–$ this gives some sense of what these quantities mean.   But since solid and fluid mechanics live in Euclidean space we don’t see some of the subtleties that arise, e.g. in the theory of general relativity.

1.5 Creating a tensor using a dyadic product of two vectors.

Let a and b be two vectors.  The dyadic product of a and b is a second order tensor S denoted by

$S=a\otimes b\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{ij}={a}_{i}{b}_{j}$.

with the property

$S\cdot u=\left(a\otimes b\right)\cdot u=a\left(b\cdot u\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{ij}{u}_{j}=\left({a}_{i}{b}_{k}\right){u}_{k}={a}_{i}\left({b}_{k}{u}_{k}\right)$

for all vectors u.  (Clearly, this maps u onto a vector parallel to a with magnitude $|a|\left(b\cdot u\right)$ )

The components of $a\otimes b$ in a basis $\left\{{e}_{1},{e}_{2},{e}_{3}\right\}$ are

$\left[\begin{array}{ccc}{a}_{1}{b}_{1}& {a}_{1}{b}_{2}& {a}_{1}{b}_{3}\\ {a}_{2}{b}_{1}& {a}_{2}{b}_{2}& {a}_{2}{b}_{3}\\ {a}_{3}{b}_{1}& {a}_{3}{b}_{2}& {a}_{3}{b}_{3}\end{array}\right]$

Note that not all tensors can be constructed using a dyadic product of only two vectors (this is because $\left(a\otimes b\right)\cdot u$ always has to be parallel to a, and therefore the representation cannot map a vector onto an arbitrary vector).  However, if a, b, and c are three independent vectors (i.e. no two of them are parallel) then all tensors can be constructed as a sum of scalar multiples of the nine possible dyadic products of these vectors.

2. OPERATIONS ON SECOND ORDER TENSORS

Tensor components

Let $\left\{{e}_{1},{e}_{2},{e}_{3}\right\}$ be a Cartesian basis, and let S be a second order tensor.  The components of S in $\left\{{e}_{1},{e}_{2},{e}_{3}\right\}$ may be represented as a matrix

$\left[\begin{array}{ccc}{S}_{11}& {S}_{12}& {S}_{13}\\ {S}_{21}& {S}_{22}& {S}_{23}\\ {S}_{31}& {S}_{32}& {S}_{33}\end{array}\right]$

where

$\begin{array}{l}{S}_{11}={e}_{1}\cdot \left(S\cdot {e}_{1}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{12}={e}_{1}\cdot \left(S\cdot {e}_{2}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{11}={e}_{1}\cdot \left(S\cdot {e}_{3}\right),\\ {S}_{21}={e}_{2}\cdot \left(S\cdot {e}_{1}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{22}={e}_{2}\cdot \left(S\cdot {e}_{2}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{21}={e}_{2}\cdot \left(S\cdot {e}_{3}\right),\\ {S}_{31}={e}_{3}\cdot \left(S\cdot {e}_{1}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{32}={e}_{3}\cdot \left(S\cdot {e}_{2}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{31}={e}_{3}\cdot \left(S\cdot {e}_{3}\right),\end{array}$

The representation of a tensor in terms of its components can also be expressed in dyadic form as

$S=\sum _{j=1}^{3}\sum _{i=1}^{3}{S}_{ij}{e}_{i}\otimes {e}_{j}$

This representation is particularly convenient when using polar coordinates, or when using a general non-orthogonal coordinate system.

#### Addition Let S and T be two tensors.  Then $U=S+T$ is also a tensor.

Denote the Cartesian components of U, S and T by matrices as defined above.  The components of U are then related to the components of S and T by

$\left[\begin{array}{ccc}{U}_{11}& {U}_{12}& {U}_{13}\\ {U}_{21}& {U}_{22}& {U}_{23}\\ {U}_{31}& {U}_{32}& {U}_{33}\end{array}\right]=\left[\begin{array}{ccc}{S}_{11}+{T}_{11}& {S}_{12}+{T}_{12}& {S}_{13}+{T}_{13}\\ {S}_{21}+{T}_{21}& {S}_{22}+{T}_{22}& {S}_{23}+{T}_{23}\\ {S}_{31}+{T}_{31}& {S}_{32}+{T}_{32}& {S}_{33}+{T}_{33}\end{array}\right]$

#### In index notation we would write

${U}_{ij}={S}_{ij}+{T}_{ij}$

#### Product of a tensor and a vector

Let u be a vector and S a second order tensor.  Then

$v=S\cdot u$

is a vector.

Let $\left({u}_{1},{u}_{2},{u}_{3}\right)$ and $\left({v}_{1},{v}_{2},{v}_{3}\right)$ denote the components of vectors u and v in a Cartesian basis $\left\{{e}_{1},{e}_{2},{e}_{3}\right\}$, and denote the Cartesian components of S as described above.  Then

$\left[\begin{array}{c}{v}_{1}\\ {v}_{2}\\ {v}_{3}\end{array}\right]=\left[\begin{array}{ccc}{S}_{11}& {S}_{12}& {S}_{13}\\ {S}_{21}& {S}_{22}& {S}_{23}\\ {S}_{31}& {S}_{32}& {S}_{33}\end{array}\right]\left[\begin{array}{c}{u}_{1}\\ {u}_{2}\\ {u}_{3}\end{array}\right]=\left[\begin{array}{c}{S}_{11}{u}_{1}+{S}_{12}{u}_{2}+{S}_{13}{u}_{3}\\ {S}_{21}{u}_{1}+{S}_{22}{u}_{2}+{S}_{23}{u}_{3}\\ {S}_{31}{u}_{1}+{S}_{32}{u}_{2}+{S}_{33}{u}_{3}\end{array}\right]$

Alternatively, using index notation

${v}_{i}={S}_{ij}{u}_{j}$

The product

$v=u\cdot S$

is also a vector.  In component form

#### or

${v}_{i}={u}_{j}{S}_{ji}$

#### Product of two tensors

Let T and S be two second order tensors.  Then $U=T\cdot S$ is also a tensor.

Denote the components of U, S and T by $3×3$ matrices.  Then,

$\begin{array}{l}\left[\begin{array}{ccc}{U}_{11}& {U}_{12}& {U}_{13}\\ {U}_{21}& {U}_{22}& {U}_{23}\\ {U}_{31}& {U}_{32}& {U}_{33}\end{array}\right]=\left[\begin{array}{ccc}{T}_{11}& {T}_{12}& {T}_{13}\\ {T}_{21}& {T}_{22}& {T}_{23}\\ {T}_{31}& {T}_{32}& {T}_{33}\end{array}\right]\left[\begin{array}{ccc}{S}_{11}& {S}_{12}& {S}_{13}\\ {S}_{21}& {S}_{22}& {S}_{23}\\ {S}_{31}& {S}_{32}& {S}_{33}\end{array}\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}=\left[\begin{array}{ccc}{T}_{11}{S}_{11}+{T}_{12}{S}_{21}+{T}_{13}{S}_{31}& {T}_{11}{S}_{12}+{T}_{12}{S}_{22}+{T}_{13}{S}_{32}& {T}_{11}{S}_{13}+{T}_{12}{S}_{23}+{T}_{13}{S}_{33}\\ {T}_{21}{S}_{11}+{T}_{22}{S}_{21}+{T}_{23}{S}_{31}& {T}_{21}{S}_{12}+{T}_{22}{S}_{22}+{T}_{23}{S}_{32}& {T}_{21}{S}_{12}+{T}_{22}{S}_{22}+{T}_{23}{S}_{32}\\ {T}_{31}{S}_{11}+{T}_{32}{S}_{21}+{T}_{33}{S}_{31}& {T}_{31}{S}_{12}+{T}_{32}{S}_{22}+{T}_{33}{S}_{32}& {T}_{31}{S}_{13}+{T}_{32}{S}_{23}+{T}_{33}{S}_{33}\end{array}\right]\end{array}$

Alternatively, using index notation

${U}_{ij}={T}_{ik}{S}_{kj}$

Note that tensor products, like matrix products, are not commutative; i.e. $T\cdot S\ne S\cdot T$

#### Transpose

Let S be a tensor.  The transpose of S is denoted by ${S}^{T}$ and is defined so that

$u\cdot {S}^{T}=S\cdot u$

Denote the components of S by a 3x3 matrix.  The components of  ${S}^{T}$ are then

${S}^{T}\equiv \left[\begin{array}{ccc}{S}_{11}& {S}_{21}& {S}_{31}\\ {S}_{12}& {S}_{22}& {S}_{32}\\ {S}_{13}& {S}_{23}& {S}_{33}\end{array}\right]$

i.e. the rows and columns of the matrix are switched.

Note that, if A and B are two tensors, then

${\left(A\cdot B\right)}^{T}={B}^{T}\cdot {A}^{T}$

Trace

Let S be a tensor, and denote the components of S by a $3×3$ matrix.  The trace of S is denoted by tr(S) or trace(S), and can be computed by summing the diagonals of the matrix of components

$\text{trace}\left(S\right)={S}_{11}+{S}_{22}+{S}_{33}$

More formally, let $\left\{{e}_{1},{e}_{2},{e}_{3}\right\}$ be any Cartesian basis.  Then

$\text{trace}\left(S\right)={e}_{1}\cdot S\cdot {e}_{1}+{e}_{2}\cdot S\cdot {e}_{2}+{e}_{3}\cdot S\cdot {e}_{3}$

The trace of a tensor is an example of an invariant of the tensor $–$ you get the same value for trace(S) whatever basis you use to define the matrix of components of S.

In index notation, the trace is written ${S}_{kk}$

Contraction.

Inner Product: Let S and T be two second order tensors.  The inner product of S and T is a scalar, denoted by $S:T$.  Represent S and T by their components in a basis.  Then

$\begin{array}{l}S:T={S}_{11}{T}_{11}+{S}_{12}{T}_{12}+{S}_{13}{T}_{13}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+{S}_{21}{T}_{21}+{S}_{22}{T}_{22}+{S}_{23}{T}_{23}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+{S}_{31}{T}_{31}+{S}_{32}{T}_{32}+{S}_{33}{T}_{33}\end{array}$

In index notation $S:T\equiv {S}_{ij}{T}_{ij}$

Observe that $S:T=T:S$, and also that $S:I=\text{trace(}S\text{)}$, where I is the identity tensor.

Outer product: Let S and T be two second order tensors.  The outer product of S and T is a scalar, denoted by $S\cdot \cdot T$.  Represent S and T by their components in a basis.  Then

$\begin{array}{l}S\cdot \cdot T={S}_{11}{T}_{11}+{S}_{21}{T}_{12}+{S}_{31}{T}_{13}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+{S}_{12}{T}_{21}+{S}_{22}{T}_{22}+{S}_{32}{T}_{23}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+{S}_{13}{T}_{31}+{S}_{23}{T}_{32}+{S}_{33}{T}_{33}\end{array}$

In index notation $S\cdot \cdot T\equiv {S}_{ij}{T}_{ji}$

Observe that $S\cdot \cdot T={S}^{T}:T$

Determinant

The determinant of a tensor is defined as the determinant of the matrix of its components in a basis.  For a second order tensor

$\begin{array}{l}\mathrm{det}S=\mathrm{det}\left[\begin{array}{l}{S}_{11}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{12}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{13}\\ {S}_{21}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{22}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{23}\text{\hspace{0.17em}}\\ {S}_{31}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{32}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{33}\end{array}\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}={S}_{11}\left({S}_{22}{S}_{33}-{S}_{23}{S}_{32}\right)+{S}_{12}\left({S}_{23}{S}_{31}-{S}_{21}{S}_{33}\right)+{S}_{13}\left({S}_{21}{S}_{32}-{S}_{31}{S}_{22}\right)\end{array}$

In index notation this would read

$\mathrm{det}\left(S\right)=\frac{1}{6}{\in }_{ijk}{\in }_{lmn}{S}_{li}{S}_{mj}{S}_{nk}\text{\hspace{0.17em}}$

Note that if S and T are two tensors, then

$\mathrm{det}\left(S\right)=\mathrm{det}\left({S}^{T}\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{det}\left(S\cdot T\right)=\mathrm{det}\left(S\right)\mathrm{det}\left(T\right)$

Inverse

Let S be a second order tensor.  The inverse of S exists if and only if $\mathrm{det}\left(S\right)\ne 0$, and is defined by

${S}^{-1}\cdot S=I$

where ${S}^{-1}$ denotes the inverse of S and I is the identity tensor.

The inverse of a tensor may be computed by calculating the inverse of the matrix of its components.  Formally, the inverse of a second order tensor can be written in a simple form using index notation as

${S}_{ji}^{-1}=\frac{1}{2\mathrm{det}\left(S\right)}{\in }_{ipq}{\in }_{jkl}{S}_{pk}{S}_{ql}$

In practice it is usually faster to compute the inverse using methods such as Gaussian elimination.

Change of Basis.

Let S be a tensor, and let $\left\{{e}_{1},{e}_{2},{e}_{3}\right\}$ be a Cartesian basis.  Suppose that the components of S in the basis $\left\{{e}_{1},{e}_{2},{e}_{3}\right\}$ are known to be

$\left[{S}^{\left(e\right)}\right]=\left[\begin{array}{ccc}{S}_{11}^{\left(e\right)}& {S}_{12}^{\left(e\right)}& {S}_{13}^{\left(e\right)}\\ {S}_{21}^{\left(e\right)}& {S}_{22}^{\left(e\right)}& {S}_{23}^{\left(e\right)}\\ {S}_{31}^{\left(e\right)}& {S}_{32}^{\left(e\right)}& {S}_{33}^{\left(e\right)}\end{array}\right]$

Now, suppose that we wish to compute the components of  S in a second Cartesian basis, $\left\{{m}_{1},{m}_{2},{m}_{3}\right\}$.  Denote these components by

$\left[{S}^{\left(m\right)}\right]=\left[\begin{array}{ccc}{S}_{11}^{\left(m\right)}& {S}_{12}^{\left(m\right)}& {S}_{13}^{\left(m\right)}\\ {S}_{21}^{\left(m\right)}& {S}_{22}^{\left(m\right)}& {S}_{23}^{\left(m\right)}\\ {S}_{31}^{\left(m\right)}& {S}_{32}^{\left(m\right)}& {S}_{33}^{\left(m\right)}\end{array}\right]$

To do so, first compute the components of the transformation matrix [Q]

$\left[Q\right]=\left[\begin{array}{l}{m}_{1}\cdot {e}_{1}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}_{1}\cdot {e}_{2}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}_{1}\cdot {e}_{3}\\ {m}_{2}\cdot {e}_{1}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}_{2}\cdot {e}_{2}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}_{2}\cdot {e}_{3}\\ {m}_{3}\cdot {e}_{1}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}_{3}\cdot {e}_{2}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}_{3}\cdot {e}_{3}\end{array}\right]$

(this is the same matrix you would use to transform vector components from $\left\{{e}_{1},{e}_{2},{e}_{3}\right\}$ to $\left\{{m}_{1},{m}_{2},{m}_{3}\right\}$ ).  Then,

$\left[{S}^{\left(m\right)}\right]=\left[Q\right]\left[{S}^{\left(e\right)}\right]{\left[Q\right]}^{T}$

or, written out in full

$\left[\begin{array}{ccc}{S}_{11}^{\left(m\right)}& {S}_{12}^{\left(m\right)}& {S}_{13}^{\left(m\right)}\\ {S}_{21}^{\left(m\right)}& {S}_{22}^{\left(m\right)}& {S}_{23}^{\left(m\right)}\\ {S}_{31}^{\left(m\right)}& {S}_{32}^{\left(m\right)}& {S}_{33}^{\left(m\right)}\end{array}\right]=\left[\begin{array}{l}{m}_{1}\cdot {e}_{1}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}_{1}\cdot {e}_{2}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}_{1}\cdot {e}_{3}\\ {m}_{2}\cdot {e}_{1}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}_{2}\cdot {e}_{2}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}_{2}\cdot {e}_{3}\\ {m}_{3}\cdot {e}_{1}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}_{3}\cdot {e}_{2}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}_{3}\cdot {e}_{3}\end{array}\right]\left[\begin{array}{ccc}{S}_{11}^{\left(e\right)}& {S}_{12}^{\left(e\right)}& {S}_{13}^{\left(e\right)}\\ {S}_{21}^{\left(e\right)}& {S}_{22}^{\left(e\right)}& {S}_{23}^{\left(e\right)}\\ {S}_{31}^{\left(e\right)}& {S}_{32}^{\left(e\right)}& {S}_{33}^{\left(e\right)}\end{array}\right]\left[\begin{array}{l}{m}_{1}\cdot {e}_{1}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}_{2}\cdot {e}_{1}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}_{3}\cdot {e}_{1}\\ {m}_{1}\cdot {e}_{2}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}_{2}\cdot {e}_{2}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}_{3}\cdot {e}_{2}\\ {m}_{1}\cdot {e}_{3}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}_{2}\cdot {e}_{3}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}_{3}\cdot {e}_{3}\end{array}\right]$

To prove this result, let u and v be vectors satisfying

$v=S\cdot u$

Denote the components of u and v in the two bases by  $\underset{_}{{u}^{\left(e\right)}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\underset{_}{{u}^{\left(m\right)}}$ and $\underset{_}{{v}^{\left(e\right)}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\underset{_}{{v}^{\left(m\right)}}$, respectively.  Recall that the vector components are related by

$\begin{array}{l}\underset{_}{{u}^{\left(m\right)}}=\left[Q\right]\underset{_}{{u}^{\left(e\right)}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\underset{_}{{u}^{\left(e\right)}}={\left[}^{Q}\underset{_}{{u}^{\left(m\right)}}\\ \underset{_}{{v}^{\left(m\right)}}=\left[Q\right]\underset{_}{{v}^{\left(e\right)}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\underset{_}{{v}^{\left(e\right)}}={\left[}^{Q}\underset{_}{{v}^{\left(m\right)}}\end{array}$

Now, we could express the tensor-vector product in either basis

$\underset{_}{{v}^{\left(m\right)}}=\left[{S}^{\left(m\right)}\right]\underset{_}{{u}^{\left(m\right)}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\underset{_}{{v}^{\left(e\right)}}=\left[{S}^{\left(e\right)}\right]\underset{_}{{u}^{\left(e\right)}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}$

Substitute for $\underset{_}{{u}^{\left(e\right)}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\underset{_}{{v}^{\left(e\right)}}$ from above into the second of these two relations, we see that

${\left[Q\right]}^{T}\underset{_}{{v}^{\left(m\right)}}=\left[{S}^{\left(e\right)}\right]{\left[Q\right]}^{T}\underset{_}{{u}^{\left(m\right)}}\text{\hspace{0.17em}}$

Recall that

$\left[Q\right]{\left[Q\right]}^{T}=\left[I\right]\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left[I\right]\underset{_}{{v}^{\left(m\right)}}=\underset{_}{{v}^{\left(m\right)}}$

so multiplying both sides by [Q] shows that

$\underset{_}{{v}^{\left(m\right)}}=\left[Q\right]\left[{S}^{\left(e\right)}\right]{\left[Q\right]}^{T}\underset{_}{{u}^{\left(m\right)}}\text{\hspace{0.17em}}$

so, comparing with the first of equation (1)

$\left[{S}^{\left(m\right)}\right]=\left[Q\right]\left[{S}^{\left(e\right)}\right]{\left[Q\right]}^{T}$

as stated.

In index notation, we would write

${S}_{ij}^{\left(m\right)}={Q}_{ik}{S}_{kl}^{\left(e\right)}{Q}_{jl}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{Q}_{ij}={m}_{i}\cdot {e}_{j}$

Another, perhaps cleaner, way to derive this result is to expand the two tensors as the appropriate dyadic products of the basis vectors

$\begin{array}{l}{S}_{kl}^{\left(m\right)}{m}_{k}\otimes {m}_{l}={S}_{kl}^{\left(e\right)}{e}_{k}\otimes {e}_{l}\\ ⇒{m}_{i}\cdot \left[{S}_{kl}^{\left(m\right)}{m}_{k}\otimes {m}_{l}\right]\cdot {m}_{j}={m}_{i}\cdot {S}_{kl}^{\left(e\right)}{e}_{k}\otimes {e}_{l}\cdot {m}_{j}\\ ⇒{S}_{ij}^{\left(m\right)}=\left({m}_{i}\cdot {e}_{k}\right){S}_{kl}^{\left(e\right)}\left({e}_{l}\cdot {m}_{j}\right)\end{array}$

Invariants

Invariants of a tensor are scalar functions of the tensor components which remain constant under a basis change.  That is to say, the invariant has the same value when computed in two arbitrary bases $\left\{{e}_{1},{e}_{2},{e}_{3}\right\}$ and $\left\{{m}_{1},{m}_{2},{m}_{3}\right\}$.  A symmetric second order tensor always has three independent invariants.

Examples of invariants are

1.      The three eigenvalues

2.      The determinant

3.      The trace

4.      The inner and outer products

These are not all independent $–$ for example any of 2-4 can be calculated in terms of 1.

In practice, the most commonly used invariants are:

$\begin{array}{l}{I}_{1}=\mathrm{trace}\left(S\right)={S}_{kk}\\ {I}_{2}=\frac{1}{2}\left(\mathrm{trace}{\left(}^{S}-S\cdot \cdot S\right)=\frac{1}{2}\left({S}_{ii}{S}_{jj}-{S}_{ij}{S}_{ji}\right)\\ {I}_{3}=\mathrm{det}\left(S\right)=\frac{1}{6}{\in }_{ijk}{\in }_{pqr}{S}_{ip}{S}_{jq}{S}_{kr}={\in }_{ijk}{S}_{i1}{S}_{j2}{S}_{k3}\end{array}$

Eigenvalues and Eigenvectors (Principal values and direction)

Let S be a second order tensor.  The scalars $\lambda$ and unit vectors m which satisfy

$S\cdot m=\lambda m$

are known as the eigenvalues and eigenvectors of S, or the principal values and principal directions of S. Note that $\lambda$ may be complex.  For a second order tensor in three dimensions, there are generally three values of $\lambda$ and three unique unit vectors m which satisfy this equation.  Occasionally, there may be only two or one value of $\lambda$.  If this is the case, there are infinitely many possible vectors m that satisfy the equation.  The eigenvalues of a tensor, and the components of the eigenvectors, may be computed by finding the eigenvalues and eigenvectors of the matrix of components.

The eigenvalues of a symmetric tensor are always real, and its eigenvectors are mutually perpendicular (these two results are important and are proved below).  The eigenvalues of a skew tensor are always pure imaginary or zero.

The eigenvalues of a second order tensor are computed using the condition $\mathrm{det}\left(S-\lambda I\right)=0$.  This yields a cubic equation, which can be expressed as

${\lambda }^{3}-{I}_{1}{\lambda }^{2}+{I}_{2}\lambda -{I}_{3}=0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{I}_{1}=trace\left(S\right),\text{​}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{I}_{2}=\left({I}_{1}^{2}-S\cdot \cdot S\right)/2\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{I}_{3}=\mathrm{det}\left(S\right)$

There are various ways to solve the resulting cubic equation explicitly $–$ a solution for symmetric S is given below, but the results for a general tensor are too messy to be given here.   The eigenvectors are then computed from the condition $\left(S-\lambda I\right)m=0$.

The Cayley-Hamilton Theorem

Let S be a second order tensor and let ${I}_{1}=trace\left(S\right),\text{​}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{I}_{2}=\left({I}_{1}^{2}-S\cdot \cdot S\right)/2\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{I}_{3}=\mathrm{det}\left(S\right)$ be the three invariants.   Then

${S}^{3}-{I}_{1}{S}^{2}+{I}_{2}S-{I}_{3}=0$

(i.e. a tensor satisfies its characteristic equation).   There is an obscure trick to show this… Consider the tensor $S-\alpha I$ (where $\alpha$ is an arbitrary scalar), and let T be the adjoint of $S-\alpha I$, (the adjoint is just the inverse multiplied by the determinant) which satisfies

$T\left(S-\alpha I\right)=\mathrm{det}\left(S-\alpha I\right)I=\left(-{\alpha }^{3}+{I}_{1}{\alpha }^{2}-{I}_{2}\alpha +{I}_{3}\right)I$

Assume that T= ${\alpha }^{2}{T}_{1}+\alpha {T}_{2}+{T}_{3}$.   Substituting in the preceding equation shows that

${T}_{1}=I\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{T}_{1}S-{T}_{2}={I}_{1}I\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{T}_{3}-{T}_{2}S={I}_{2}I\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{T}_{3}S={I}_{3}I$

Use these to substitute for ${I}_{1},{I}_{2},{I}_{3}$ into

${S}^{3}-{I}_{1}{S}^{2}+{I}_{2}S-{I}_{3}={S}^{3}-\left(S-{T}_{2}\right){S}^{2}+\left({T}_{3}-{T}_{2}S\right)S-{T}_{3}S=0$

3 SPECIAL TENSORS

Identity tensor  The identity tensor I is the tensor such that, for any tensor S or vector v

$\begin{array}{l}I\cdot v=v\cdot I=v\\ S\cdot I=I\cdot S=S\end{array}$

In any basis, the identity tensor has components

$\left[\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 0& 1\end{array}\right]$

Symmetric Tensor A symmetric tensor S has the property

$S={S}^{T}$

The components of a symmetric tensor have the form

$\left[\begin{array}{ccc}{S}_{11}& {S}_{12}& {S}_{13}\\ {S}_{12}& {S}_{22}& {S}_{23}\\ {S}_{13}& {S}_{23}& {S}_{33}\end{array}\right]$

so that there are only six independent components of the tensor, instead of nine.   Symmetric tensors have some nice properties:

·         The eigenvectors of a symmetric tensor with distinct eigenvalues are orthogonal.   To see this, let $u,v$ be two eigenvectors, with corresponding eigenvalues ${\lambda }_{u},{\lambda }_{v}$. Then $v\cdot \left[S\cdot u\right]=u\cdot \left[{S}^{T}\cdot v\right]=u\cdot \left[S\cdot v\right]⇒v\cdot {\lambda }_{u}u=u\cdot {\lambda }_{v}v⇒\left({\lambda }_{u}-{\lambda }_{v}\right)u\cdot v=0⇒u\cdot v=0$.

·         The eigenvalues of a symmetric tensor are realTo see this, suppose that $\lambda ,u$ are a complex eigenvalue/eigenvector pair, and let $\overline{\lambda },\overline{u}$ denote their complex conjugates.   Then, by definition $S\cdot u=\lambda u⇒S\overline{u}=\overline{\lambda }\overline{u}$.  And hence $\overline{u}\cdot \left[S\cdot u\right]=\lambda \overline{u}\cdot u,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}u\cdot \left[S\cdot \overline{u}\right]=\overline{\lambda }u\cdot \overline{u}$.  But note that for a symmetric tensor $\overline{u}\cdot \left[S\cdot u\right]=u\cdot \left[{S}^{T}\cdot \overline{u}\right]=u\cdot \left[S\cdot \overline{u}\right]$.  Thus $\lambda \overline{u}\cdot u=\overline{\lambda }u\cdot \overline{u}⇒\lambda =\overline{\lambda }$.

The eigenvalues of a symmetric tensor can be computed as

$\begin{array}{l}{\lambda }_{k}=\frac{{I}_{1}}{3}+2\sqrt{\frac{-p}{3}}\mathrm{cos}\left\{\frac{1}{3}{\mathrm{cos}}^{-1}\left(\frac{3q}{2p}\sqrt{\frac{-3}{p}}\right)-\frac{2\left(k-1\right)\pi }{3}\right\}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}k=1,2,3\\ p={I}_{2}-\frac{1}{3}{I}_{1}^{2}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}q=-\frac{2{I}_{1}^{3}-9{I}_{1}{I}_{2}+27{I}_{3}}{27}\\ {I}_{1}=trace\left(S\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{I}_{2}=\frac{1}{2}\left({I}_{1}^{2}-S:S\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{I}_{3}=\mathrm{det}\left(S\right)\end{array}$

The eigenvectors can then be found by back-substitution into $\left[S-\lambda I\right]\cdot m=0$. To do this, note that the matrix equation can be written as

$\left[\begin{array}{ccc}{S}_{11}-\lambda & {S}_{12}& {S}_{13}\\ {S}_{12}& {S}_{22}-\lambda & {S}_{23}\\ {S}_{13}& {S}_{23}& {S}_{33}-\lambda \end{array}\right]\left[\begin{array}{c}{m}_{1}\\ {m}_{2}\\ {m}_{3}\end{array}\right]=\left[\begin{array}{c}0\\ 0\\ 0\end{array}\right]$

Since the determinant of the matrix is zero, we can discard any row in the equation system and take any column over to the right hand side.   For example, if the tensor has at least one eigenvector with ${m}_{3}\ne 0$ then the values of ${m}_{1},{m}_{2}$ for this eigenvector can be found by discarding the third row, and writing

$\left[\begin{array}{cc}{S}_{11}-\lambda & {S}_{12}\\ {S}_{12}& {S}_{22}-\lambda \end{array}\right]\left[\begin{array}{c}{m}_{1}\\ {m}_{2}\end{array}\right]=-{m}_{3}\left[\begin{array}{c}{S}_{13}\\ {S}_{23}\end{array}\right]$

·         Spectral decomposition of a symmetric tensor  Let S be a symmetric second order tensor, and let $\left\{{\lambda }_{i},{e}_{i}\right\}$ be the three eigenvalues and eigenvectors of S.  Then S can be expressed as

$S=\sum _{i=1}^{3}{\lambda }_{i}{e}_{i}\otimes {e}_{i}$

To see this, note that S can always be expanded as a sum of 9 dyadic products of an orthogonal basis.  $S={S}_{ij}{e}_{i}\otimes {e}_{j}$But since ${e}_{i}$ are eigenvectors it follows that ${e}_{m}\cdot \left({S}_{ij}{e}_{i}\otimes {e}_{j}\right)\cdot {e}_{k}={S}_{mk}=\left\{\begin{array}{c}{\lambda }_{k}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}m=k\\ 0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}m\ne k\end{array}$

Skew Tensor.  A skew tensor S has the property

${S}^{T}=-S$

The components of a skew tensor have the form

$\left[\begin{array}{ccc}0& {S}_{12}& {S}_{13}\\ -{S}_{12}& 0& {S}_{23}\\ -{S}_{13}& -{S}_{23}& 0\end{array}\right]$

Every second-order skew tensor has a dual vector w that satisfies

$S\cdot u=w×u$

for all vectors u.   You can see this by noting that ${S}_{12}=-{w}_{3}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{13}={w}_{2}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{S}_{23}=-{w}_{1}$ and expanding out the tensor and cross products explicitly.  In index notation, we can also write

${S}_{ij}=-{\in }_{ijk}{w}_{k}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{w}_{i}=-\frac{1}{2}{\in }_{ijk}{S}_{jk}$.

Orthogonal Tensors An orthogonal tensor R has the property

$\begin{array}{l}R\cdot {R}^{T}={R}^{T}\cdot R=I\\ {R}^{-1}={R}^{T}\end{array}$

An orthogonal tensor must have $\mathrm{det}\left(S\right)=±1$; a tensor with $\mathrm{det}\left(S\right)=+1$ is known as a proper orthogonal tensor.  Orthogonal tensors also have some interesting and useful properties:

·         Orthogonal tensors map a vector onto another vector with the same length.  To see this, let u be an arbitrary vector.  Then, note that ${|R\cdot u|}^{2}=\left[R\cdot u\right]\cdot \left[R\cdot u\right]=u\cdot {R}^{T}R\cdot u=u\cdot u={|u|}^{2}$

·         The eigenvalues of an orthogonal tensor are $1,{e}^{±i\theta }$ for some value of $\theta$.  To see this, let $u$ be an eigenvector, with corresponding eigenvalue $\lambda$.  By definition, $R\cdot u=\lambda u$.  Hence, $\left[R\cdot u\right]\cdot \left[R\cdot u\right]=\lambda u\cdot \lambda u⇒u\cdot u={\lambda }^{2}u\cdot u⇒{\lambda }^{2}=1$. Similarly, $\lambda \overline{\lambda }=1$.  Since the characteristic equation is cubic, there must be at most three eigenvalues, and at least one eigenvalue must be real.

Proper orthogonal tensors can be visualized physically as rotations.  A rotation can also be represented in several other forms besides a proper orthogonal tensor.   For example

·         The Rodriguez representation quantifies a rotation as an angle of rotation $\theta$ (in radians) about some axis n (specified by a unit vector).  Given R, there are various ways to compute n  and $\theta$.  For example, one way would be find the eigenvalues and the real eigenvector.  The real eigenvector (suitably normalized) must correspond to n; the complex eigenvalues give ${e}^{i\theta }$.  A faster method is to note that

$\text{trace}\left(R\right)=1+2\mathrm{cos}\theta \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}2\mathrm{sin}\theta n=\text{dual}\left(R-{R}^{T}\right)$

·         Alternatively, given n and $\theta$, R can be computed from

$R=\mathrm{cos}\theta I+W\cdot W\left(1-\mathrm{cos}\theta \right)+W\mathrm{sin}\theta$

where W is the skew tensor that has n as its dual vector, i.e. ${W}_{ij}=-{\in }_{ijk}{n}_{k}$.  In index notation, this formula is

${R}_{ij}=\mathrm{cos}\theta {\delta }_{ij}+{n}_{i}{n}_{j}\left(1-\mathrm{cos}\theta \right)-\mathrm{sin}\theta {\in }_{ijk}{n}_{k}$

Another useful result is the Polar Decomposition Theorem, which states that invertible second order tensors can be expressed as a product of a symmetric tensor with an orthogonal tensor:

$A=R\cdot U=V\cdot R\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}R{R}^{T}=I\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}U={U}^{T}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}V={V}^{T}$

Moreover, the tensors $R,U,V$ are unique.  To see this, note that

·         ${A}^{T}A$ is symmetric and has positive eigenvalues (to see that it’s symmetric, simply take the transpose, and to see that the eigenvalues are positive, note that $dx\cdot \left({A}^{T}\cdot A\right)\cdot dx>0$ for all vectors dx).

·         Let  ${\lambda }_{k}^{2}$ and ${m}_{k}$ be the three eigenvalues and eigenvectors of ${A}^{T}A$.  Since the eigenvectors are orthogonal, we can write ${A}^{T}A=\sum _{k=1}^{3}{\lambda }_{k}^{2}{m}_{k}\otimes {m}_{k}$

·         We can then set $U=\sum _{k=1}^{3}{\lambda }_{k}{m}_{k}\otimes {m}_{k}$ and define $R=A{U}^{-1}$U is clearly symmetric, and also ${U}^{2}={A}^{T}A$. To see that R is orthogonal note that: ${R}^{T}R={U}^{-T}{A}^{T}A{U}^{-1}={U}^{-T}{U}^{2}{U}^{-1}=I$.

·         Given that U and R exist we can write $R\cdot U=\left[R\cdot U\cdot {R}^{T}\right]\cdot R$ so if we define $V=RU{R}^{T}$ then $A=VR$.   It is easy to show that V is symmetric.

·         To see that the decomposition is unique, suppose that $A=\stackrel{^}{R}\stackrel{^}{U}$ for some other tensors $\stackrel{^}{R},\stackrel{^}{U}$.  Then ${A}^{T}A={\stackrel{^}{U}}^{2}$.  But ${A}^{T}A$ has a unique square root so $\stackrel{^}{U}=U$.  The uniqueness of R follows immediately.