Section 1.
10
1.10 Special Second Order Tensors & Properties of
Second Order Tensors
In this section will be examined a number of special second order tensors, and special
properties of second order tensors, which play important roles in tensor analysis. Many
of the concepts will be familiar from Linear Algebra and Matrices. The following will be
discussed:
The Identity tensor
Transpose of a tensor
Trace of a tensor
Norm of a tensor
Determinant of a tensor
Inverse of a tensor
Orthogonal tensors
Rotation Tensors
Change of Basis Tensors
Symmetric and Skew-symmetric tensors
Axial vectors
Spherical and Deviatoric tensors
Positive Definite tensors
1.10.1 The Identity Tensor
The linear transformation which transforms every tensor into itself is called the identity
tensor. This special tensor is denoted by I so that, for example,
Ia a for any vector a
In particular, Ie1 e1 , Ie 2 e 2 , Ie 3 e 3 , from which it follows that, for a Cartesian
coordinate system, I ij ij . In matrix form,
1 0 0
I 0 1 0 (1.10.1)
0 0 1
1.10.2 The Transpose of a Tensor
The transpose of a second order tensor A with components Aij is the tensor A T with
components A ji ; so the transpose swaps the indices,
A Aij e i e j , A T A ji e i e j Transpose of a Second-Order Tensor (1.10.2)
Solid Mechanics Part III 84 Kelly
Section 1.10
In matrix notation,
A11 A12 A13 A11 A21 A31
A A21 A22 A23 , A T
A12 A22 A32
A31 A32 A33 A13 A23 A33
Some useful properties and relations involving the transpose are {▲Problem 2}:
A T T
A
A B A T B T
T
u v T v u
Tu uT T , uT T T u
(1.10.3)
AB T
B A T T
A : B AT : BT
(u v ) A u ( A T v )
A : BC B T A : C AC T : B
A formal definition of the transpose which does not rely on any particular coordinate
system is as follows: the transpose of a second-order tensor is that tensor which satisfies
the identity1
u Av v A T u (1.10.4)
for all vectors u and v. To see that Eqn. 1.10.4 implies 1.10.2, first note that, for the
present purposes, a convenient way of writing the components Aij of the second-order
tensor A is A ij . From Eqn. 1.9.4, A ij e i Ae j and the components of the transpose
can be written as A T ij e i A T e j . Then, from 1.10.4,
A T
ij e i A T e j e j Ae i A ji A ji .
1.10.3 The Trace of a Tensor
The trace of a second order tensor A, denoted by trA , is a scalar equal to the sum of the
diagonal elements of its matrix representation. Thus (see Eqn. 1.4.3)
trA Aii Trace (1.10.5)
A more formal definition, again not relying on any particular coordinate system, is
trA I : A Trace (1.10.6)
1
as mentioned in §1.9, from the linearity of tensors, uA v u Av and, for this reason, this expression is
usually written simply as uAv
Solid Mechanics Part III 85 Kelly
Section 1.10
and Eqn. 1.10.5 follows from 1.10.6 {▲Problem 4}. For the dyad u v {▲Problem 5},
tr u v u v (1.10.7)
Another example is
tr(E 2 ) I : E 2
ij e i e j : E pq E qr e p e r (1.10.8)
Eiq E qi
This and other important traces, and functions of the trace are listed here {▲Problem 6}:
trA Aii
trA 2 Aij A ji
trA 3 Aij A jk Aki (1.10.9)
trA 2 Aii A jj
trA 3 Aii A jj Akk
Some useful properties and relations involving the trace are {▲Problem 7}:
trA T trA
tr AB tr BA
tr A B trA trB (1.10.10)
tr A trA
A : B tr A T B tr AB T tr B T A tr BA T
The double contraction of two tensors was earlier defined with respect to Cartesian
coordinates, Eqn. 1.9.16. This last expression allows one to re-define the double
contraction in terms of the trace, independent of any coordinate system.
Consider again the real vector space of second order tensors V 2 introduced in §1.8.5.
The double contraction of two tensors as defined by 1.10.10e clearly satisfies the
requirements of an inner product listed in §1.2.2. Thus this scalar quantity serves as an
inner product for the space V 2 :
A, B A : B tr A T B (1.10.11)
and generates an inner product space.
Just as the base vectors e i form an orthonormal set in the inner product (vector dot
product) of the space of vectors V, so the base dyads e i e j form an orthonormal set in
the inner product 1.10.11 of the space of second order tensors V 2 . For example,
Solid Mechanics Part III 86 Kelly
Section 1.10
e1 e1 , e1 e1 e1 e1 : e1 e1 1 (1.10.12)
Similarly, just as the dot product is zero for orthogonal vectors, when the double
contraction of two tensors A and B is zero, one says that the tensors are orthogonal,
A : B tr A T B 0 , A, B orthogonal (1.10.13)
1.10.4 The Norm of a Tensor
Using 1.2.8 and 1.10.11, the norm of a second order tensor A, denoted by A (or A ), is
defined by
A A:A (1.10.14)
This is analogous to the norm a of a vector a, aa .
1.10.5 The Determinant of a Tensor
The determinant of a second order tensor A is defined to be the determinant of the
matrix A of components of the tensor:
A11 A12 A13
det A det A21 A22 A23
A31 A32 A33
ijk Ai1 A j 2 Ak 3 (1.10.15)
ijk A1i A2 j A3k
Some useful properties of the determinant are {▲Problem 8}
det( AB) det A det B
det A T det A
det(A ) 3 det A
(1.10.16)
det u v 0
pqr det A ijk Aip A jq Akr
Ta Tb Tc det Ta b c
Note that det A , like trA , is independent of the choice of coordinate system / basis.
Solid Mechanics Part III 87 Kelly
Section 1.10
1.10.6 The Inverse of a Tensor
The inverse of a second order tensor A, denoted by A 1 , is defined by
AA 1 I A 1 A (1.10.17)
The inverse of a tensor exists only if it is non-singular (a singular tensor is one for
which det A 0 ), in which case it is said to be invertible.
Some useful properties and relations involving the inverse are:
( A 1 ) 1 A
(A) 1 (1 / ) A 1
(1.10.18)
( AB) 1 B 1 A 1
det( A 1 ) (det A) 1
Since the inverse of the transpose is equivalent to the transpose of the inverse, the
following notation is used:
A T ( A 1 ) T ( A T ) 1 (1.10.19)
1.10.7 Orthogonal Tensors
An orthogonal tensor Q is a linear vector transformation satisfying the condition
Qu Qv u v (1.10.20)
for all vectors u and v. Thus u is transformed to Qu , v is transformed to Qv and the dot
product u v is invariant under the transformation. Thus the magnitude of the vectors
and the angle between the vectors is preserved, Fig. 1.10.1.
Q Qu
v
Qv
Figure 1.10.1: An orthogonal tensor
Since
Qu Qv uQT Qv u QTQ v (1.10.21)
Solid Mechanics Part III 88 Kelly
Section 1.10
it follows that for u v to be preserved under the transformation, Q T Q I , which is also
used as the definition of an orthogonal tensor. Some useful properties of orthogonal
tensors are {▲Problem 10}:
QQ T I Q T Q, Qik Q jk ij Qki Qkj
Q 1 Q T (1.10.22)
det Q 1
1.10.8 Rotation Tensors
If for an orthogonal tensor, det Q 1 , Q is said to be a proper orthogonal tensor,
corresponding to a rotation. If det Q 1 , Q is said to be an improper orthogonal
tensor, corresponding to a reflection. Proper orthogonal tensors are also called rotation
tensors.
1.10.9 Change of Basis Tensors
Consider a rotation tensor Q which rotates the base vectors e1 , e 2 , e 3 into a second set,
e1 , e2 , e3 , Fig. 1.10.2.
ei Qe i i 1, 2, 3 (1.10.23)
Such a tensor can be termed a change of basis tensor from e i to ei . The transpose
QT rotates the base vectors ei back to e i and is thus change of basis tensor from ei to
e i . The components of Q in the e i coordinate system are, from 1.9.4, Qij e i Qe j and
so, from 1.10.23,
Q Qij e i e j , Qij e i ej , (1.10.24)
which are the direction cosines between the axes (see Fig. 1.5.5).
e3 e3
e2
e2 Q
e1
e1
Figure 1.10.2: Rotation of a set of base vectors
Solid Mechanics Part III 89 Kelly
Section 1.10
The change of basis tensor can also be expressed in terms of the base vectors from both
bases:
Q ei e i , (1.10.25)
from which the above relations can easily be derived, for example ei Qe i , QQ T I ,
etc.
Consider now the operation of the change of basis tensor on a vector:
Qv vi Qe i vi ei (1.10.26)
Thus Q transforms v into a second vector v , but this new vector has the same
components with respect to the basis ei , as v has with respect to the basis e i , vi vi .
Note the difference between this and the coordinate transformations of §1.5: here there
are two different vectors, v and v .
Example
Consider the two-dimensional rotation tensor
0 1
Q e i e j ei e i
1 0
which corresponds to a rotation of the base vectors through / 2 . The vector v 1 1
T
then transforms into (see Fig. 1.10.3)
1 1
Qv e i ei
1 1
e1 e2
Qv v
e2 e1
Figure 1.10.3: a rotated vector
Similarly, for a second order tensor A, the operation
QAQ T QAij e i e j Q T Aij Qe i e j Q T Aij Qe i Qe j Aij ei ej
(1.10.27)
Solid Mechanics Part III 90 Kelly
Section 1.10
results in a new tensor which has the same components with respect to the ei , as A has
with respect to the e i , Aij Aij .
1.10.10 Symmetric and Skew Tensors
A tensor T is said to be symmetric if it is identical to the transposed tensor, T T T , and
skew (antisymmetric) if T T T .
Any tensor A can be (uniquely) decomposed into a symmetric tensor S and a skew tensor
W, where
symA S
1
2
A AT
(1.10.28)
1
skewA W A A T
2
and
S ST , W W T (1.10.29)
In matrix notation one has
S11 S12 S13 0 W12 W13
S S12 S 22 S 23 ,
W W12 0 W23 (1.10.30)
S13 S 23 S 33 W13 W23 0
Some useful properties of symmetric and skew tensors are {▲Problem 13}:
S : B S : B T S : 12 B B T
W : B W : B W : 12 B T
T
S:W 0
(1.10.31)
tr SW 0
v Wv 0
det W 0 has no inverse
where v and B denote any arbitrary vector and second-order tensor respectively.
Note that symmetry and skew-symmetry are tensor properties, independent of coordinate
system.
Solid Mechanics Part III 91 Kelly
Section 1.10
1.10.11 Axial Vectors
A skew tensor W has only three independent coefficients, so it behaves “like a vector”
with three components. Indeed, a skew tensor can always be written in the form
Wu ω u (1.10.32)
where u is any vector and ω characterises the axial (or dual) vector of the skew tensor
W. The components of W can be obtained from the components of ω through
Wij e i We j e i ω e j e i k e k e j
e i k kjp e p kji k (1.10.33)
ijk k
If one knows the components of W, one can find the components of ω by inverting this
equation, whence {▲Problem 14}
ω W23e1 W13e 2 W12 e 3 (1.10.34)
Example (of an Axial Vector)
Decompose the tensor
1 2 3
T Tij 4 2 1
1 1 1
into its symmetric and skew parts. Also find the axial vector for the skew part. Verify
that Wa ω a for a e1 e3 .
Solution
One has
1 2 3 1 4 1 1 3 2
1
1
S T T 4 2 1 2 2 1 3
2
T
2
2 1
1 1 1 3 1 1 2 1 1
1 2 3 1 4 1 0 1 1
1 1
W T T 4 2 1 2 2 1 1
T
0 0
2 2
1 1 1 3 1 1 1 0 0
The axial vector is
ω W23e1 W13e 2 W12 e 3 e 2 e 3
Solid Mechanics Part III 92 Kelly
Section 1.10
and it can be seen that
0 1 1 1 1 Wa Wij (e i e j )(e1 e 3 ) Wij ( j1 j 3 )e i (Wi1 Wi 3 )e i
1 0 0 0 1 or (W11 W13 )e1 (W21 W23 )e 2 (W31 W33 )e 3
1 0 0 1 1 e1 e 2 e 3
and
e1 e 2 e3
ωa 0 1 1 e1 e 2 e 3
1 0 1
■
The Spin Tensor
The velocity of a particle rotating in a rigid body motion is given by v ω x , where ω
is the angular velocity vector and x is the position vector relative to the origin on the axis
of rotation (see Problem 9, §1.1). If the velocity can be written in terms of a skew-
symmetric second order tensor w , such that wx v , then it follows from wx ω x
that the angular velocity vector ω is the axial vector of w . In this context, w is called
the spin tensor.
1.10.12 Spherical and Deviatoric Tensors
Every tensor A can be decomposed into its so-called spherical part and its deviatoric
part, i.e.
A sphA devA (1.10.35)
where
sphA 1
3
trA I
13 A11 A22 A33 0 0
0 1
3
A11 A22 A33 0
0 0 1
3
A11 A22 A
33
devA A sphA
A11 13 A11 A22 A33 A12 A13
A21 A22 13 A11 A22 A33 A23
A31 A32 A33 13 A11 A22 A33
(1.10.36)
Solid Mechanics Part III 93 Kelly
Section 1.10
Any tensor of the form I is known as a spherical tensor, while devA is known as a
deviator of A, or a deviatoric tensor.
Some important properties of the spherical and deviatoric tensors are
tr (devA ) 0
sph (devA ) 0 (1.10.37)
devA : sphB 0
1.10.13 Positive Definite Tensors
A positive definite tensor A is one which satisfies the relation
vAv 0 , v o (1.10.38)
The tensor is called positive semi-definite if vAv 0 .
In component form,
vi Aij v j A11v12 A12 v1v 2 A13 v1v3 A21v 2 v1 A22 v 22 (1.10.39)
and so the diagonal elements of the matrix representation of a positive definite tensor
must always be positive.
It can be shown that the following conditions are necessary for a tensor A to be positive
definite (although they are not sufficient):
(i) the diagonal elements of A are positive
(ii) the largest element of A lies along the diagonal
(iii) det A 0
(iv) Aii A jj 2 Aij for i j (no sum over i, j )
These conditions are seen to hold for the following matrix representation of an example
positive definite tensor:
2 2 0
A 1 4 0
0 0 1
A necessary and sufficient condition for a tensor to be positive definite is given in the
next section, during the discussion of the eigenvalue problem.
One of the key properties of a positive definite tensor is that, since det A 0 , positive
definite tensors are always invertible.
An alternative definition of positive definiteness is the equivalent expression
Solid Mechanics Part III 94 Kelly
Section 1.10
A:vv 0 (1.10.40)
1.10.14 Problems
1. Show that the components of the (second-order) identity tensor are given by I ij ij .
2. Show that
(a) (u v ) A u ( A T v )
(b) A : BC B T A : C AC T : B
3. Use (1.10.4) to show that I T I .
4. Show that (1.10.6) implies (1.10.5) for the trace of a tensor.
5. Show that tr u v u v .
6. Formally derive the index notation for the functions
trA 2 , trA 3 , ( trA ) 2 , ( trA) 3
7. Show that A : B tr ( A T B) .
8. Prove (1.10.16f), Ta Tb Tc det T a b c .
9. Show that ( A 1 ) T : A 3 . [Hint: one way of doing this is using the result from
Problem 7.]
10. Use 1.10.16b and 1.10.18d to prove 1.10.22c, det Q 1 .
11. Use the explicit dyadic representation of the rotation tensor, Q ei e i , to show that
the components of Q in the “second”, ox1 x 2 x3 , coordinate system are the same as
those in the first system [hint: use the rule Qij ei Qej ]
12. Consider the tensor D with components (in a certain coordinate system)
1 / 2 1 / 2 1 / 2
0 1/ 2 1/ 2
1 / 2 1 / 2 1 / 2
Show that D is a rotation tensor (just show that D is proper orthogonal).
13. Show that tr SW 0 .
14. Multiply across (1.10.32), Wij ijk k , by ijp to show that ω 12 ijk Wij e k . [Hint:
use the relation 1.3.19b, ijp ijk 2 pk .]
15. Show that 12 a b b a is a skew tensor W. Show that its axial vector is
ω 12 b a . [Hint: first prove that b u a a u b u a b b a u .]
16. Find the spherical and deviatoric parts of the tensor A for which Aij 1 .
Solid Mechanics Part III 95 Kelly