Invariants of tensors

In mathematics, in the fields of multilinear algebra and representation theory, invariants of tensors are coefficients of the characteristic polynomial of the tensor A:[1]

\ p(\lambda):=\det (\mathbf{A}-\lambda \mathbf{E}),

where \mathbf E is the identity tensor and \lambda\in\mathbb C is the polynomial's indeterminate (it is important to bear in mind that a polynomial's indeterminate \lambda may also be a non-scalar as long as power, scaling and adding make sense for it, e.g., p(\mathbf A) is legitimate, and in fact, quite useful).

The first invariant of an n×n tensor A (I_A) is the coefficient for \lambda^{n-1} (because the coefficient for  \lambda^n is always 1), the second invariant (II_A) is the coefficient for \lambda^{n-2}, etc., the nth invariant is the free term.

The definition of the invariants of tensors and specific notations used throughout the article were introduced into the field of rheology by Ronald Rivlin and became extremely popular there. In fact even the trace of a tensor A is usually denoted as I_A in the textbooks on rheology.

Properties

The first invariant (trace) is always the sum of the diagonal components:

\ I_A=A_{11}+A_{22}+ \cdots + A_{nn}=\mathrm{tr}(\mathbf{A}) \,

The nth invariant is just \pm \det \mathbf{A}, the determinant of \mathbf{A} (up to sign).

The invariants do not change with rotation of the coordinate system (they are objective). Obviously, any function of the invariants only is also objective.

Calculation of the invariants of symmetric 3×3 tensors

Most tensors used in engineering are symmetric 3×3. For this case the invariants can be calculated as:


  \mathrm{I}_A= \mathrm{tr}(\mathbf{A}) = A_{11}+A_{22}+A_{33} = A_1+A_2+A_3 \,

\mathrm{II}_A  = \frac{1}{2} \left( ( \mathrm{tr}\mathbf{A})^2 - \mathrm{tr}(\mathbf{A}^2 ) \right) = A_{11}A_{22}+A_{22}A_{33}+A_{11}A_{33}-A_{12}A_{21}-A_{23}A_{32}-A_{13}A_{31}

(the sum of principal minors)


\begin{align}
& = A_1A_2+A_2A_3+A_1A_3 \\
\mathrm{III}_A & = \det (\mathbf{A})= A_1 A_2 A_3
\end{align}

where A_1, A_2, A_3 are the eigenvalues of tensor A.

Because of the Cayley–Hamilton theorem the following equation is always true:

\ \mathbf{A}^3 - \mathrm{I}_A \mathbf{A}^2 +\mathrm{II}_A \mathbf{A} -\mathrm{III}_A \mathbf{E}= 0

where E is the second-order identity tensor.

A similar equation holds for tensors of higher order.

Engineering application

A scalar valued tensor function f that depends merely on the three invariants of a symmetric 3×3 tensor \mathbf{A} is objective, i.e., independent from rotations of the coordinate system. Moreover, every objective tensor function depends only on the tensor's invariants. Thus, objectivity of a tensor function is fulfilled if, and only if, for some function \ g:R^3\to R we have

\ f(\mathbf A)=g(I_A,II_A,III_A). \,

A common application to this is the evaluation of the potential energy as function of the strain tensor, within the framework of linear elasticity. Exhausting the above theorem the free energy of the system reduces to a function of 3 scalar parameters rather than 6. Within linear elasticity the free energy has to be quadratic in the tensor's elements, which eliminates an additional scalar. Thus, for an isotropic material only two independent parameters are needed to describe the elastic properties, known as Lamé coefficients. Consequently, experimental fits and computational efforts may be eased significantly.

Other invariants

Since the invariants are constant in any reference, functions of invariants are also constant. Some sources define the three invariants of the 3×3 tensors as

I_1 = tr(\mathbf{A}) = A_{ii}
I_2 = tr(\mathbf{A}^2) = A_{ij}A_{ji}
I_3 = tr(\mathbf{A}^3) = A_{ij}A_{jk}A_{ki}

Therefore,

I_A = I_1
II_A = \frac{1}{2}(I_1^2 - I_2)
III_A = \frac{1}{6}(I_1^3 - 3 I_1 I_2 + 2 I_3)

See also

References

  1. SPENCER, A. J. M. Continuum Mechanics. Longman, 1980.
This article is issued from Wikipedia - version of the Friday, November 06, 2015. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.