A Brief Introduction To Tensors and Their Properties
A Brief Introduction To Tensors and Their Properties
A Brief Introduction To Tensors and Their Properties
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
. Let
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
:
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:
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
be a vector field, and let
the definition of G
Now, let
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
You have probably already seen the matrix representation of stress and strain components in
introductory courses.
, 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
has the same magnitude and direction, irrespective of the basis used to represent u, v, and S.
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
matrix,
To check whether
check how
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 dont, then the matrix cannot be the components of a tensor). Then,
evaluate the matrix product to find a new left hand side, say
. If
are related to
by the
same transformation that was used to calculate the new components of u, then
are the components of a vector, and, therefore, the matrix represents the
components of a tensor.
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
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
of real numbers
which are invertible functions of
. Tensors can then be
defined as sets of real numbers that transform in a particular way under this change in coordinate
system. For example
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 dont see some of the subtleties that arise, e.g. in the theory of
general relativity.
Let a and b be two vectors. The dyadic product of a and b is a second order tensor S denoted by
.
with the property
for all vectors u. (Clearly, this maps u onto a vector parallel to a with magnitude
)
The components of
in a basis
are
Note that not all tensors can be constructed using a dyadic product of only two vectors (this is because
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.
Tensor components.
Let
The components of S in
where
The representation of a tensor in terms of its components can also be expressed in dyadic form as
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
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
is a vector.
Let
and
The product
denote
the
, and denote
or
Observe that
(unless S is symmetric).
is also a tensor.
matrices. Then,
Note that tensor products, like matrix products, are not commutative; i.e.
Transpose
are then
Trace
Contraction.
Inner Product: Let S and T be two second order tensors. The inner product of S and T is a scalar,
denoted by
. Represent S and T by their components in a basis. Then
In index notation
Observe that
Outer product: Let S and T be two second order tensors. The outer product of S and T is a scalar,
denoted by
. Represent S and T by their components in a basis. Then
In index notation
Observe that
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
Inverse
where
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
In practice it is usually faster to compute the inverse using methods such as Gaussian elimination.
Change of Basis.
Now, suppose that we wish to compute the components of S in a second Cartesian basis,
. Denote these components by
(this is the same matrix you would use to transform vector components from
to
). Then,
and
Substitute for
we see that
Recall that
as stated.
Another, perhaps cleaner, way to derive this result is to expand the two tensors as the appropriate
dyadic products of the basis vectors
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
and
symmetric second order tensor always has three independent invariants.
are known as the eigenvalues and eigenvectors of S, or the principal values and principal directions
of S. Note that
may be complex. For a second order tensor in three dimensions, there are
generally three values of
and three unique unit vectors m which satisfy this equation.
Occasionally, there may be only two or one value of
. 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
. This yields a cubic equation, which can be expressed as
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
be the
three invariants. Then
(i.e. a tensor satisfies its characteristic equation). There is an obscure trick to show this Consider
the tensor
(where
is an arbitrary scalar), and let T be the adjoint of
, (the adjoint is just the inverse multiplied by the determinant) which satisfies
Assume that T=
preceding equation shows that
. Substituting in the
into
3 SPECIAL TENSORS
Identity tensor The identity tensor I is the tensor such that, for any tensor S or vector v
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
. Then
The eigenvalues of a symmetric tensor are real. To see this, suppose that
complex eigenvalue/eigenvector pair, and let
Then,
by
are a
definition
.
And hence
. But note that for a
symmetric
tensor
.
Thus
.
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
then the values of
discarding the third row, and writing
To see this, note that S can always be expanded as a sum of 9 dyadic products of an
orthogonal basis.
. But since
are
; a tensor with
Orthogonal tensors map a vector onto another vector with the same length. To see this, let u
be an arbitrary vector. Then, note that
. Similarly,
. 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
where W is the skew tensor that has n as its dual vector, i.e.
. In index notation, this formula is
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:
Let
for all
and
and
define
. Then
for some
. But
. The uniqueness of R follows