The Kinematic Equations

David Roylance
Department of Materials Science and Engineering
Massachusetts Institute of Technology
Cambridge, MA 02139

September 19, 2000

Introduction
The kinematic or strain-displacement equations describe how the strains – the stretching and
distortion – within a loaded body relate to the body’s displacements. The displacement com-
ponents in the x, y, and z directions are denoted by the vector u ≡ ui ≡ (u, v, w), and are
functions of position within the body: u = u(x, y, z). If all points within the material experi-
ence the same displacement (u = constant), the structure moves as a rigid body, but does not
stretch or deform internally. For stretching to occur, points within the body must experience
different displacements.

Infinitesimal strain

Figure 1: Incremental deformation.

Consider two points A and B separated initially by a small distance dx as shown in Fig. 1, and
experiencing motion in the x direction. If the displacement at point A is uA , the displacement
at B can be expressed by a Taylor’s series expansion of u(x) around the point x = A:

∂u
uB = uA + du = uA + dx
∂x
where here the expansion has been truncated after the second term. The differential motion δ
between the two points is then
 
∂u ∂u
δ = uB − uA = uA + dx − uA = dx
∂x ∂x
In our concept of stretching as being the differential displacement per unit length, the x com-
ponent of strain is then

1
 δ ∂u
x = = (1)
dx ∂x
Hence the strain is a displacement gradient. Applying similar reasoning to differential motion
in the y direction, the y-component of strain is the gradient of the vertical displacement v with
respect to y:

∂v
y = (2)
∂y

Figure 2: Shearing distortion.

The distortion of the material, which can be described as the change in originally right
angles, is the sum of the tilts imparted to vertical and horizontal lines. As shown in Fig. 2,
the tilt of an originally vertical line is the relative horizontal displacement of two nearby points
along the line:
 
∂u ∂u
δ = uB − uA = uA + dy − uA = dy
∂y ∂y
The change in angle is then

δ ∂u
γ1 ≈ tan γ1 = =
dy ∂y
Similarly (see Fig. 3), the tilt γ2 of an originally horizontal line is the gradient of v with respect
to x. The shear strain in the xy plane is then

∂v ∂u
γxy = γ1 + γ2 = + (3)
∂x ∂y
This notation, using  for normal strain and γ for shearing strain, is sometimes known as the
“classical” description of strain.

Matrix Formulation
The “indicial notation” described in the Module on Matrix and Index Notation provides a concise
method of writing out all the components of three-dimensional states of strain:
!
1 ∂ui ∂uj 1
ij = + ≡ (ui,j + uj,i ) (4)
2 ∂xj ∂xi 2

2
 Figure 3: Shearing strain.

where the comma denotes differentiation with respect to the following spatial variable. This
double-subscript index notation leads naturally to a matrix arrangement of the strain compo-
nents, in which the i-j component of the strain becomes the matrix element in the ith row and
the j th column:
     
∂u 1 ∂u ∂v 1 ∂u ∂w
  ∂x  2 ∂y + ∂x 2  ∂z
+ ∂x  
 1 ∂u ∂v ∂v 1 ∂v ∂w 
ij =  2  ∂y + ∂x  ∂z +  (5)
  ∂y  2 ∂y 
1 ∂w ∂u 1 ∂v ∂w ∂w
2 ∂x + ∂z 2 ∂z + ∂y ∂z

Note that the strain matrix is symmetric, i.e. ij = ji . This symmetry means that there are six
rather than nine independent strains, as might be expected in a 3 × 3 matrix. Also note that the
indicial description of strain yields the same result for the normal components as in the classical
description: 11 = x . However, the indicial components of shear strain are half their classical
counterparts: 12 = γxy /2.
In still another useful notational scheme, the classical strain-displacement equations can be
written out in a vertical list, similar to a vector:
   

 x 
 
 ∂u/∂x 


 y 
 
 


 
 
 ∂v /∂ y 


 
 
 

 z ∂w/∂ z
=

 γyz 
 
 ∂v/∂z + ∂w/∂ y 


 
 
 


 γxz 
 
 ∂u/∂z + ∂w/∂x 


 
 
 

γxy ∂u/∂y + ∂v/∂x
This vector-like arrangement of the strain components is for convenience only, and is sometimes
called a pseudovector. Strain is actually a second-rank tensor, like stress or moment of inertia,
and has mathematical properties very different than those of vectors. The ordering of the
elements in the pseudovector form is arbitrary, but it is conventional to list them as we have
here by moving down the diagonal of the strain matrix of Eqn. 5 from upper left to lower right,
then move up the third column, and finally move one column to the left on the first row; this
gives the ordering 1,1; 2,2; 3,3; 2,3; 1,3; 1,2.
Following the rules of matrix multiplication, the strain pseudovector can also be written in

3
terms of the displacement vector as
   

 x 
 ∂/∂x 0 0

 
   

 y 
  0 ∂/∂y 0  u 

 
   
z  0 0 ∂/∂z 
=  v (6)

 γyz 
  0 ∂/∂z ∂/∂y  

 
   w 

 γxz  
 ∂/∂z 0 ∂/∂x 

 

γxy ∂/∂y ∂/∂x 0

The matrix in brackets above, whose elements are differential operators, can be abbreviated as
L:
 
∂/∂x 0 0
 0 ∂/∂y 0 
 
 
 0 0 ∂/∂y 
L=  (7)
 0 ∂/∂z ∂/∂y 
 
 ∂/∂z 0 ∂/∂x 
∂/∂y ∂/∂x 0
The strain-displacement equations can then be written in the concise “pseudovector-matrix”
form:

 = Lu (8)
Equations such as this must be used in a well-defined context, as they apply only when the
somewhat arbitrary pseudovector listing of the strain components is used.

Volumetric strain
Since the normal strain is just the change in length per unit of original length, the new length
L0 after straining is found as

L0 − L0
= ⇒ L0 = (1 + )L0 (9)
L0
If a cubical volume element, originally of dimension abc, is subjected to normal strains in all
three directions, the change in the element’s volume is

∆V a0 b0 c0 − abc a(1 + x ) b(1 + y ) c(1 + z ) − abc

= =
V abc abc

= (1 + x ) (1 + y ) (1 + z ) − 1 ≈ x + y + z (10)
where products of strains are neglected in comparison with individual values. The volumetric
strain is therefore the sum of the normal strains, i.e. the sum of the diagonal elements in the
strain matrix (this is also called the trace of the matrix, or Tr[]). In index notation, this can
be written simply

∆V
= kk
V
This is known as the volumetric, or “dilatational” component of the strain.

4
 Example 1
To illustrate how volumetric strain is calculated, consider a thin sheet of steel subjected to strains in its
plane given by x = 3, y = −4, and γxy = 6 (all in µin/in). The sheet is not in plane strain, since it can
undergo a Poisson strain in the z direction given by z = −ν(x + y ) = −0.3(3 − 4) = 0.3. The total
state of strain can therefore be written as the matrix
 
3 6 0
[] =  6 −4 0  × 10−6
0 0 0.3
where the brackets on the [] symbol emphasize that the matrix rather than pseudovector form of the
strain is being used. The volumetric strain is:
∆V
= (3 − 4 + 0.3) × 10−6 = −0.7 × 10−6
V
Engineers often refer to “microinches” of strain; they really mean microinches per inch. In the case of
volumetric strain, the corresponding (but awkward) unit would be micro-cubic-inches per cubic inch.

Finite strain
The infinitesimal strain-displacement relations given by Eqns. 3.1–3.3 are used in the vast major-
ity of mechanical analyses, but they do not describe stretching accurately when the displacement
gradients become large. This often occurs when polymers (especially elastomers) are being con-
sidered. Large strains also occur during deformation processing operations, such as stamping of
steel automotive body panels. The kinematics of large displacement or strain can be complicated
and subtle, but the following section will outline a simple description of Lagrangian finite strain
to illustrate some of the concepts involved.
Consider two orthogonal lines OB and OA as shown in Fig. 4, originally of length dx and
dy, along the x-y axes, where for convenience we set dx = dy = 1. After strain, the endpoints of
these lines move to new positions A1 O1 B1 as shown. We will describe these new positions using
the coordinate scheme of the original x-y axes, although we could also allow the new positions
to define a new set of axes. In following the motion of the lines with respect to the original
positions, we are using the so-called Lagrangian viewpoint. We could alternately have used the
final positions as our reference; this is the Eulerian view often used in fluid mechanics.
After straining, the distance dx becomes
 
0 ∂u
(dx) = 1 + dx
∂x
Using our earlier “small” thinking, the x-direction strain would be just ∂u/∂x. But when the
strains become larger, we must also consider that the upward motion of point B1 relative to O1 ,
that is ∂v/∂x, also helps stretch the line OB. Considering both these effects, the Pythagorean
theorem gives the new length O1 B1 as
s
 2  2
∂u ∂v
O1 B1 = 1+ +
∂x ∂x
We now define our Lagrangian strain as

5
 Figure 4: Finite displacements.

O1 B1 − OB
x = = O1 B1 − 1
OB
s
 2  2
∂u ∂u ∂v
= 1+2 + + −1
∂x ∂x ∂x
√
Using the series expansion 1 + x = 1 + x/2 + x2 /8 + · · · and neglecting terms beyond first
order, this becomes
( "  2  2 #)
1 ∂u ∂u ∂v
x ≈ 1+ 2 + + −1
2 ∂x ∂x ∂x
" 2  2 #
∂u 1 ∂u ∂v
= + + (11)
∂x 2 ∂x ∂x
Similarly, we can show
" 2  2 #
∂v 1 ∂v ∂u
y = + + (12)
∂y 2 ∂y ∂y
∂u ∂v ∂u ∂u ∂v ∂v
γxy = + + + (13)
∂y ∂x ∂y ∂x ∂y ∂x
When the strains are sufficiently small that the quadratic terms are negligible compared with
the linear ones, these reduce to the infinitesimal-strain expressions shown earlier.

Example 2
The displacement function u(x) for a tensile specimen of uniform cross section and length L, fixed at
one end and subjected to a displacement δ at the other, is just the linear relation
x
u(x) = δ
L
The Lagrangian strain is then given by Eqn. 11 as

6
  2
δ 1 δ
x = +
L 2 L
The first term is the familiar small-strain expression, with the second nonlinear term becoming more
important as δ becomes larger. When δ = L, i.e. the conventional strain is 100%, there is a 50%
difference between the conventional and Lagrangian strain measures.

The Lagrangian strain components can be generalized using index notation as

1
ij =(ui,j + uj,i + ur,i ur,j ).
2
A pseudovector form is also convenient occasionally:
 
       u,x 
  
 x 
 
 u,x 
 u v
1  ,x ,x
0 0 
 v 

 ,x

y = v,y +  0 0 u,y v,y 

 γ 
   u +v   2 u  u 
xy ,y ,x ,y v,y u,x v,x  

,y 


v,y
  
    ∂/∂x 0  
 ∂/∂x 0 u,x v,x 0 0  0 ∂/∂x  
=  0 ∂/∂y + 1 0 0 u,y v,y    u
 2  ∂/∂y 0  v
∂/∂y ∂/∂x u,y v,y u,x v,x
0 ∂/∂y
which can be abbreviated

 = [L + A(u)] u (14)
The matrix A(u) contains the nonlinear effect of large strain, and becomes negligible when
strains are small.

Problems
1. Write out the abbreviated strain-displacement equation  = Lu (Eqn. 8) for two dimen-
sions.

2. Write out the components of the Lagrangian strain tensor in three dimensions:

1
ij = (ui,j + uj,i + ur,i ur,j )
2

3. Show that for small strains the fractional volume change is the trace of the infinitesimal
strain tensor:

∆V
≡ kk = x + y + z
V

4. When the material is incompressible, show the extension ratios are related by

λ x λy λz = 1

7
5. Show that the kinematic (strain-displacement) relations in for polar coordinates can be
written

∂ur
r =
∂r

1 ∂uθ ur
θ = +
r ∂θ r

1 ∂ur ∂uθ uθ
γrθ = + −
r ∂θ ∂r r

8
MIT OpenCourseWare
http://ocw.mit.edu

3.11 Mechanics of Materials

Fall 1999

For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.