Stress-Strain Relations for

Elastic Isotropic Materials

Isotropic materials have identical properties in all direc-
tions, such that the physical properties of the material
are the same regardless of the direction in which they
are measured.
• The constitutive equation for isotropic materials that
relates a general state of stress to strain is based on
a generalized form of Hooke’s law.

For an isotropic material, the stress is linearly propor-
tional to the strain and is described by two material
constants. For example,
I Young’s modulus E and Poisson’s ratio ν.
I Alternatively, the stress-strain relationship can be
parameterized by the Lamé constants λ and µ (where
µ is also known as the shear modulus).
• The stress within a solid represents the internal forces
per unit area that develop within a material in response
to external loads.
• These stress components are crucial in understanding
how a material will deform and potentially fail under
various loading conditions.
LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 1 of 45
Physical Interpretation of Stress Components
1. Normal Stresses: (σ11, σ22, σ33) represent forces act-
ing perpendicular to the respective planes within the
material. Positive values indicate tension, while neg-
ative values indicate compression.
2. Shear Stresses: (σ12, σ13, σ23) represent forces acting
parallel to the respective planes within the material.
Shear stresses tend to deform the material by sliding
layers over one another.

By understanding these components, engineers can an-

alyze how materials respond to different types of load-
ing conditions, such as tension, compression, shear, bend-
ing, and torsion, and design structures to withstand these
loads safely.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 2 of 45

Stress Component Indexes i and j
y
x2

σyy
σ22

σ21 σyx
σ12 σyz σxy
σ23
σ32 σzy

σ11 σxx
σ13 σxz
σ31 σzx
σ33 x1 σzz x

x3 z

Stress components on volume element – in the limit, this element

shrinks to a point.

To grasp the physical meaning of the surface and direc-

tion of these stress components, let’s break down the
notation σij :
1. Subscript i: First subscript i
w Direction of the normal to the perpendicular sur-
face on which the stress acts.
2. Subscript j: Second subscript j
w Direction of the stress component acting on that
plane. This indicates the direction in which the inter-
nal force is applied.
LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 3 of 45
 Interpretation of Stress Components σij
In a three-dimensional space, stress is represented by a
second-order tensor (or matrix) with nine components:
 
σ11 σ12 σ13


σ = σ21 σ22 σ23
 



σ31 σ32 σ33
 

42 2. FORCES AND BALANCE PRINCIPLES

Index convention for σij components: Symbol σij represents the

F IGURE 2.5.1. Index convention for ij components: Symbol ij represents the jth com-
jth component offorce
ponent of the surface theacting
surface forcewhose
on the surface acting
normalon the tosurface
is parallel the xi -[Link]
All components
normal is parallelshown
toherein
the are positive. All components shown are positive.
xi-axis.
are positive. In the context of the normal components, our sign convention implies that tension is positive
while compression is negative.6
Each component index (ij) can be interpreted as in
Note that, if we wished, we could reverse our sign convention without affecting Newton’s third law of
the following Table.
motion. In other words,

(2.5.3) T(n) = T( n)

is always true; herein, T(n) is due to the force acting on the plane whose outward normal is n, and T( n)

isLinkedIn
due to the
� -force acting
Dr. Lonny on the other
Thompson, side of
Clemson this plane.
University, Sept. 1, 2024 4 of 45

As formulated herein, the entries of matrix (2.5.2) determine the stress state within a continuum at a given
 Stress Component Description
σ11 Surface: Acts on a surface where the normal is aligned with the x1
(or x) axis.
Direction: Acts in the x1 direction (along the x-axis).
Meaning: Normal stress acting perpendicular to a plane normal to
the x1 -axis.
σ22 Surface: Acts on a surface where the normal is aligned with the x2
(or y) axis.
Direction: Acts in the x2 direction (along the y-axis).
Meaning: Normal stress acting perpendicular to a plane normal to
the x2 -axis.
σ33 Surface: Acts on a surface where the normal is aligned with the x3
(or z) axis.
Direction: Acts in the x3 direction (along the z-axis).
Meaning: Normal stress acting perpendicular to a plane normal to
the x3 -axis.
σ12 , σ21 Surface: σ12 acts on a surface normal to the x1 -axis.
Direction: σ12 acts in the x2 direction (parallel to the y-axis).
Surface: σ21 acts on a surface normal to x2 -axis.
Direction: σ21 acts in the x1 direction (parallel to the x-axis).
Meaning: Shear stress tries to slide one layer of the material relative
to another along the x2 -direction while acting on a plane
normal to the x1 -direction, or vice versa.
σ13 , σ31 Surface: σ13 acts on a surface normal to the x1 -axis.
Direction: σ13 acts in the x3 direction (along the z-axis).
Surface: σ31 acts on a surface normal to the x3 -axis.
Direction: σ31 acts in the x1 direction (along the x-axis).
Meaning: Shear stress tries to slide material layers along the z-
direction while acting on a plane normal to the x1 -
direction, or vice versa.
σ23 , σ32 Surface: σ23 acts on a surface normal to the x2 -axis.
Direction: σ23 acts in the x3 direction (along the z-axis).
Surface: σ32 acts on a surface normal to the x3 -axis.
Direction: σ32 acts in the x2 direction (along the y-axis).
Meaning: Shear stress tries to slide material layers along the z-
direction while acting on a plane normal to the x2 -
LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 5 of 45
direction, or vice versa.
Symmetry of the Stress Tensor
In most real-world materials (especially in isotropic and
linearly elastic materials), the stress tensor is symmetric,
meaning that σij = σji.
x2

σ22

σ21
σ23 σ12
σ32
σ11
σ13
σ31
σ33 x1

x3


This symmetry reflects the fact that shear stresses on
perpendicular planes must be equal in magnitude but
act in opposite directions to maintain rotational equi-
librium. Therefore:

σ12 = σ21, σ13 = σ31, σ23 = σ32

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 6 of 45

Constitutive Equation for Isotropic Materials
The constitutive equations for isotropic linear elastic ma-
terials describe the relationship between stress and strain
tensors in terms of material properties.

The material is isotropic, meaning its mechanical prop-
erties are the same in all directions. As a result, the
relationship between stress and strain is independent
of the coordinate system.

When expressed in terms of Lamé parameters λ and
µ, these equations are derived from Hooke’s law:
σij = λδij tr() + 2µij
- σij is the stress tensor component in the ith direc-
tion and jth direction.
- ij is the strain tensor component in the ith direc-
tion and jth direction.
- δij is the Kronecker delta, which is 1 if i = j and 0
otherwise.
- tr() = P3k=1 kk = 11 + 22 + 33 is the trace of the
strain tensor, representing the volumetric strain.
- λ and µ are the Lamé parameters, with λ being the
first Lamé parameter and µ (also denoted as G )
being the second, equivalent to shear modulus.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 7 of 45

Interpretation of terms in Constitutive Equation
I Both material parameters λ and µ represent the stiff-
ness or elasticity of the material to load and deforma-
tion.

Expanding the volumetric strain, the relationship be-

tween stress and strain for isotropic materials can be ex-
pressed:

σij = λ (11 + 22 + 33) δij + 2µij

1. The first term λ (11 + 22 + 33) δij accounts for the vol-
umetric change, which contributes to the normal stresses
but not to shear stresses.
2. The second term 2µij represents the deviatoric (shear)
stress component, which accounts for changes in shape
without changing the volume.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 8 of 45

Relationship with Other Elastic Constants
The Lamé parameters λ and µ are related to other elas-
tic constants like Young’s modulus E, Poisson’s ratio ν
and the bulk modulus K by:
Eν
λ=
(1 + ν)(1 − 2ν)
E
µ=G=
2(1 + ν)
2µ E
K =λ+ =
3 3(1 − 2ν)
• Young’s Modulus E is a measure of the stiffness of a
material in tension or compression.
• Poisson’s Ratio ν describes the ratio of lateral strain to
axial strain in a material subjected to uniaxial stress.
• Shear Modulus G (which is the same as µ ) describes
the material’s response to shear stress.
• Bulk Modulus K is a measure of a material’s incom-
pressibility: the higher the bulk modulus, the less com-
pressible the material is under uniform pressure.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 9 of 45

Relationship to Engineering Constants
Understanding Lamé constants and their relationship to
engineering constants (e.g., E, G, ν, and K) allows en-
gineers to convert between different sets of material
properties which can be necessary when dealing with
different types of analyses or when working with data
from various sources.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 10 of 45

Applications of the Constitutive Equations
• Stress Analysis: The constitutive equations are fun-
damental in stress analysis and design of structures,
particularly in predicting how materials will respond
under various loading conditions.
• FEA Modeling: These equations are used in finite
element analysis (FEA) software to simulate material
behavior in engineering applications, allowing for ac-
curate predictions of stress and strain distributions.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 11 of 45

Let’s dive deeper into the constitutive equa-
tion relating stress to strain for isotropic
materials
Constitutive Equation in Tensor Form:
The general form of the constitutive equation for isotropic
materials in terms of the stress tensor σij components
and the strain tensor components ij is: for i and j rang-
ing from 1 to 3 dimensions:
3
σij = λδij kk + 2µij
X

k=1

- σij : symmetric stress tensor.

- ij : symmetric strain tensor.
- kk = tr() = 11 + 22 + 33 is the trace of the strain
tensor, representing the volumetric strain.
- δij is the Kronecker delta (identity tensor), which is 1
when i = j and 0 otherwise.
- λ: first Lamé constant.
- µ: second Lamé constant, also known as the shear
modulus.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 12 of 45

Lamé Constants Terminology
• Lamé Constants are perhaps the most commonly used
term, emphasizing values are fixed for a given mate-
rial.
• Lamé Parameters is another widely accepted term. It
is especially common in mathematical contexts where
these values are treated as parameters within equa-
tions or models.
• The use of the term constants underscores the fixed
nature of these values for a given material, while pa-
rameters highlight their role within the material’s con-
stitutive equations.
• In some contexts, λ and µ might be grouped under
the broader term ‘elastic constants,’ which includes
other material properties like Young’s modulus E, Pois-
son’s ratio ν, shear modulus G, and bulk modulus K.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 13 of 45

Matrix Representation of the Tensors
For a 3D isotropic material, the components of the second-
order stress and strain tensors, as well as the identity
tensor, can be represented as 3 × 3 matrices.

σ = λ tr() I + 2µ 

Stress σ , Strain  and Identify Tensors:

     
σ11 σ12 σ13  11 12 13   1 0 0 
     
σ = σ21 σ22 σ23 ,  = 21 22 23 , I =  0 1 0 
     
  

σ31 σ32 σ33 31 32 33 0 0 1

     

x2

σ22

σ21
σ23 σ12
σ32
σ11
σ13
σ31
σ33 x1

x3

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 14 of 45

Constitutive Equation in Matrix Form:
Using the above matrices, the constitutive equation can
be expressed:
     


σ11 σ12 σ13 1 0 0
 
11 12 13 
 
σ21 σ22 σ23 = λ (11 +22 +33) 0 1 0 +2µ 21 22 23
     
    


σ31 σ32 σ33 0 0 1 31 32 33
     

This can be further expanded to:

 
λ (11 + 22 + 33 ) + 2µ 11 2µ 12 2µ 13
σ= 2µ 21 λ (11 + 22 + 33 ) + 2µ 22 2µ 23
 
 

2µ 31 2µ 32 λ (11 + 22 + 33 ) + 2µ 33

Component Form of the Constitutive Equation:

In index notation, the constitutive equation can be ex-
panded for each component of the stress tensor:
σ11 = λ (11 + 22 + 33) + 2µ11
σ22 = λ (11 + 22 + 33) + 2µ22
σ33 = λ (11 + 22 + 33) + 2µ33
σ12 = 2µ12, σ23 = 2µ23, σ31 = 2µ31

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 15 of 45

Importance of the Constitutive Equation
The constitutive equation for isotropic materials expresses
the linear relationship between stress and strain tensors.

It shows that stress in a material is proportional to
both the volumetric strain (through the Lamé constant
λ ) and the shear strain (through the shear modulus µ).

This equation is fundamental in elasticity theory and
forms the basis for analyzing stress and deformation
in isotropic materials under various loading conditions.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 16 of 45

Relationship with Young’s Modulus and Poisson’s Ra-
tio:
The Lamé constants λ and µ can be related to the com-
monly used material properties: Young’s modulus E and
Poisson’s ratio ν.
Eν E
λ= , µ=
(1 + ν)(1 − 2ν) 2(1 + ν)
• Young’s modulus E represents a measure of stiffness
or elasticity.
• Poisson’s ratio ν is the ratio of lateral strain to axial
strain under uniaxial stress.
Using these relationships, the constitutive equation can
also be expressed in terms of E and ν.

Eν E
σij = δij kk + ij
(1 + ν)(1 − 2ν) 1+ν
We can factor out the common terms in the normal
stress components to simplify the expressions:

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 17 of 45

Normal Stress Components:
E
σ11 = [ν (11 + 22 + 33) + (1 − ν)11]
(1 + ν)(1 − 2ν)
E
σ22 = [ν (11 + 22 + 33) + (1 − ν)22]
(1 + ν)(1 − 2ν)
E
σ33 = [ν (11 + 22 + 33) + (1 − ν)33]
(1 + ν)(1 − 2ν)
Shear Stress Components: The shear stress compo-
nents remain as follows:
E
σ12 = 12
1+ν
E
σ23 = 23
1+ν
E
σ31 = 31
1+ν

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 18 of 45

 In summary, the constitutive relations for each com-
ponent of the stress tensor in terms of the strain tensor,
Young’s modulus E, and Poisson’s ratio ν are:
E
σ11 = [ν (11 + 22 + 33) + (1 − ν)11]
(1 + ν)(1 − 2ν)
E
σ22 = [ν (11 + 22 + 33) + (1 − ν)22]
(1 + ν)(1 − 2ν)
E
σ33 = [ν (11 + 22 + 33) + (1 − ν)33]
(1 + ν)(1 − 2ν)
E E E
σ12 = 12, σ23 = 23, σ31 = 31
1+ν 1+ν 1+ν

x2

σ22

σ12
σ23 σ12
σ23
σ11
σ13
σ13
σ33 x1

x3

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 19 of 45

Inverse Relationship
(Strain as a Function of Stress)
To find the inverse relationship, we express the strain
components ij in terms of the stress components σij .
For an isotropic material, the inverse constitutive re-
lationship can be derived as:

1+ν ν
ij = σij − δij σkk
E E
where:
ij is the strain tensor.
σij is the stress tensor.
σkk is the trace of the stress tensor, representing the
sum of the normal stresses: σkk = σ11 + σ22 + σ33.
E is Young’s modulus.
ν is Poisson’s ratio.

Expanded Form of the Inverse Constitutive

Equations:
Expanding the above equation for each component of
the strain tensor:
LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 20 of 45
 Normal Strain Components:
For the normal strain components 11, 22, and 33 :

1
11 = [σ11 − ν (σ22 + σ33)]
E
1
22 = [σ22 − ν (σ11 + σ33)]
E
1
33 = [σ33 − ν (σ11 + σ22)]
E
Shear Strain Components:
For the shear strain components 12, 23, and 13 :

1+ν
12 = σ12
E
1+ν
23 = σ23
E
1+ν
13 = σ13
E

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 21 of 45

Summary of Inverse Constitutive Equations:
The inverse constitutive relations for strain in terms of
stress for an isotropic material are:
1
11 = [σ11 − ν (σ22 + σ33)]
E
1
22 = [σ22 − ν (σ11 + σ33)]
E
1
33 = [σ33 − ν (σ11 + σ22)]
E
1+ν 1+ν 1+ν
12 = σ12, 23 = σ23, 31 = σ31
E E E
These Generalized Hooke’s law equations describe how
the strain in each direction is related to the stresses in
all three directions, accounting for the material’s elas-
ticity properties E and ν. The normal strains are influ-
enced by both the corresponding normal stresses and
the stresses in the perpendicular directions (due to Pois-
son’s effect). In contrast, the shear strains are directly
proportional to the corresponding shear stresses.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 22 of 45

Simplified Notations
• The notations σnxσny and σnz are sometimes used by
engineers to emphasize that these stress components
act normal to the surface,
• In the context of elasticity, these stresses are typically
referred to as normal stresses. The subscript n is used
to explicitly indicate that these stresses are normal to
the respective surfaces.

Using this notation:

1
x = [σnx − ν (σny + σnz )]
E
1
y = [σny − ν (σnx + σnz )]
E
1
z = [σnz − ν (σnx + σny )]
E
τxy τyz τxz
γxy = ; γyz = ; γxz =
G G G
where

γxy = 2xy , γyz = 2yz , γxz = 2xz

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 23 of 45

Alternative Notations for Stress Components
• Normal Stress Notation: In traditional notation, nor-
mal stresses are denoted as σxx, σyy and σzz , where
the first subscript indicates the plane on which the
stress acts, and the second subscript indicates the di-
rection in which the stress is applied. For example,
σxx represents the normal stress acting on the plane
perpendicular to the x-axis in the x-direction.
• Alternative Notation: σnx, σny and σnz :

– Used to emphasize the normal nature of these stress

components. Here, n stands for ‘normal,’ making
it clear that these stresses are perpendicular to the
surface or plane of interest.
– Used where it’s important to distinguish between
normal stresses (which act perpendicular to a sur-
face) and shear stresses (which act parallel to a sur-
face). Using this notation helps reinforce the con-
cept that these are specifically normal stresses.

• Shear stress Notation: Often the shear stress com-

ponents σ12 = τxy , σ23 = τyz , and σ13 = τxz are used
with ‘tau’ (τ ) notation to distinguish shear from nor-
mal stress components.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 24 of 45

General Inverse Constitutive Equation:
The inverse relationship, where the strain tensor  is ex-
pressed as a function of the stress tensor σ, is given by:

1
 = ((1 + ν) σ − ν tr(σ) I)
E
σ is the stress tensor.
 is the strain tensor.
tr() is the trace of the strain tensor, which is a scalar.
I is the identity tensor.
λ and µ are the Lamé constants.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 25 of 45

Inverse Constitutive Equation in Matrix Form:
Using the matrix representations of the stress tensor σ,
strain tensor , and identity matrix I, the inverse consti-
tutive equation relating strain to stress is:
1
= ((1 + ν) σ − ν tr(σ) I)
E
In matrix form:
    
σ11 σ12 σ13 σ11 + σ22 + σ33 0 0
1 
 = (1 + ν) σ21 σ22 σ23  − ν  0 σ11 + σ22 + σ33 0
   
   
E σ31 σ32 σ33 0 0 σ11 + σ22 + σ33


Simplifying,
 
σ − ν(σ22 + σ33) (1 + ν)σ12 (1 + ν)σ13
1  11 

 =  (1 + ν)σ21 σ22 − ν(σ11 + σ33) (1 + ν)σ23


E


(1 + ν)σ31 (1 + ν)σ32 σ33 − ν(σ11 + σ22)
 

• This matrix equation shows how each component of

the strain tensor  is related to the components of the
stress tensor σ.
• The normal strains are influenced by all the normal
stresses due to Poisson’s effect, while the shear strains
are directly proportional to the corresponding shear
stresses.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 26 of 45

Summary in Complete Tensor Notation:
Stress as a Function of Strain:

σ = λ tr() I + 2µ 

or equivalently,
Eν E
σ= tr() I + 
(1 + ν)(1 − 2ν) 1+ν
Strain as a Function of Stress:
1
 = ((1 + ν) σ − ν tr(σ) I)
E
These equations in complete tensor notation are com-
pact and express the relationship between stress and
strain in a form that is independent of the coordinate
system, making them particularly useful for analyzing
stress and strain in complex three-dimensional prob-
lems involving isotropic materials.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 27 of 45

Engineering Matrix Form
To express the constitutive equations for isotropic ma-
terials in matrix stiffness form using Voigt notation, we
represent the stress and strain tensors as column vec-
tors and the stiffness matrix (also known as the elastic
matrix) as a 6 × 6 matrix. In Voigt notation, the six in-
dependent components of the stress and strain tensors
are arranged as follows:
   
σ11 σxx


σ22 σyy 
  



σ33  σzz 
  


σ= 
 =  
σ23  τyz 




σ13  τxz 




σ12  τxy 


Similarly, for the Strain column vector:

   
11   xx 


22   yy 
   



33   zz 
   


= 
  =  
223 γyz 









213 γxz 





212 γxy 


LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 28 of 45

Constitutive Equation in Voigt Matrix Form:
The constitutive equation in matrix form relates the stress
vector σ to the strain vector  via the stiffness matrix C:

σ = C

For an isotropic material, the stiffness matrix C (also

called the elastic stiffness matrix or material stiffness
matrix) is given by:
 
C11 

C12 C12 0 0 0 
C12 C11 C12 0 0 0 
 



C12 C12 C11 0 0 0 
 


C=


 
0 



0 0 C44 0 0 
0 



0 0 0 C44 0 
0 
0 0 0 0 C44 
The coefficients C11, C12, and C44 are related to the ma-
terial properties: Young’s modulus E and Poisson’s ratio
ν as follows:
E(1 − ν)
C11 =
(1 + ν)(1 − 2ν)
Eν
C12 =
(1 + ν)(1 − 2ν)
E
C44 = =G
2(1 + ν)

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 29 of 45

Expanded Constitutive Equations:
The explicit form of the elastic matrix equation σ = C
is:
    


σ11   C11 C12 C12 0 0 0   11 
σ22   C12 C11 C12 0 0 0   22 
    



σ33   C12 C12 C11 0 0 0   33 
    
  
   
 = 

  




σ23   0 0 0 C44 0 0   223 

 




σ13   0 0 0 0 C44 0   213 

σ12   0 0 0 0 0 C44   212 
This expands to the following six equations:

σ11 = C1111 + C1222 + C1233

σ22 = C1211 + C1122 + C1233
σ33 = C1211 + C1222 + C1133
σ23 = C44 · 223
σ13 = C44 · 213
σ12 = C44 · 212

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 30 of 45

 Alternatively, using notation to better distinguish be-
tween normal and shear stress and engineering strain
angles:
E
σxx = [(1 − ν)xx + ν(yy + zz )]
(1 + ν)(1 − 2ν)
E
σyy = [(1 − ν)yy + ν(xx + zz )]
(1 + ν)(1 − 2ν)
E
σzz = [(1 − ν)zz + ν(xx + yy )]
(1 + ν)(1 − 2ν)
τyz = G γyz
τxz = G γxz
τxy = G γxy

σyy

τxy
τyz τxy
τyz

σxx
τxz
τxz
σzz x

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 31 of 45

 In Voigt notation, the constitutive equations for an
isotropic material can be expressed as:

σ = C
(1−ν)
 
ν ν

 (1+ν)(1−2ν) (1+ν)(1−2ν) (1+ν)(1−2ν) 0 0 0 

 
ν (1−ν) ν
0 0 0
 
 
(1+ν)(1−2ν) (1+ν)(1−2ν) (1+ν)(1−2ν)
 
 
 

ν ν (1−ν) 

(1+ν)(1−2ν) (1+ν)(1−2ν) (1+ν)(1−2ν) 0 0 0 

C=E
 
 
 
1
0 0 0 0 0
 
 
2(1+ν)
 
 
1
 
0 0 0 0 0
 
 

 2(1+ν) 


1 

0 0 0 0 0 2(1+ν)


• This matrix relates the six independent stress com-

ponents to the six independent strain components,
using the material properties E and ν for an isotropic
material.
• The stiffness matrix C for an isotropic material can
also be expressed more compactly in terms of the
Lamé coefficients λ and µ.
• The Lamé coefficients are related to the material’s elas-
tic properties, and using them allows for a more fun-
damental representation of the material behavior.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 32 of 45

Stiffness Matrix in Terms of
Lamé Coefficients
The stiffness (elastic) matrix C in Voigt notation using
the Lamé coefficients λ and µ is compactly expressed
by:
 
λ + 2µ

λ λ 0 0 0 
λ λ + 2µ λ 0 0 0 
 



λ λ λ + 2µ 0 0 0 
 


C=


 
0 



0 0 µ 0 0 
0 



0 0 0 µ 0 
0 
0 0 0 0 µ
• The block in the first three rows and columns indicates
the coupling of normal stress with each normal strain
due to direct strain in that direction and the influence
of the other two directions (due to Poisson’s effect).
• The µ terms on the diagonal in the lower part of the
matrix represent the material’s response to shear stresses.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 33 of 45

Expanded Form of the Matrix Equations
Expanding the matrix equation σ = C, we get the fol-
lowing relationships:

σ11 = (λ + 2µ)11 + λ22 + λ33

σ22 = λ11 + (λ + 2µ)22 + λ33
σ33 = λ11 + λ22 + (λ + 2µ)33
σ23 = µ · 223
σ13 = µ · 213
σ12 = µ · 212

which, after collecting, are the same relations as ob-

tained earlier with the tensor expressions:
σ11 = λ (11 + 22 + 33) + 2µ11
σ22 = λ (11 + 22 + 33) + 2µ22
σ33 = λ (11 + 22 + 33) + 2µ33
σ12 = 2µ12, σ23 = 2µ23, σ31 = 2µ31

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 34 of 45

Compliance Matrix
The compliance matrix S is the inverse of the stiffness
matrix C, and relates the strain vector  to the stress
vector σ. When expressed in Voigt notation, the com-
pliance matrix S is given by:

 = Sσ

For an isotropic material, the compliance matrix S can

be expressed in terms of E and ν as:
 
1 

−ν −ν 0 0 0 

−ν 1 −ν 0 0 0
 
 
 
 
1 −ν −ν 1 0 0 0
 
 

S = C−1
 
= 

1+ν


E 0 0 0 0 0
 
2
 
 
1+ν
 
0 

 0 0 0 2 0 


1+ν
 
0 
0 0 0 0 2


This compliance matrix S can also be rewritten in terms

of the Lamé coefficients λ and µ by using the relation-
ships between the material properties:
µ(3λ + 2µ) λ
E= , ν=
λ+µ 2(λ + µ)

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 35 of 45

Misconception: Misinterpreting the Stiffness
and Compliance Matrices
We might be tempted to think of the stiffness (elasticity)
matrix C or the compliance matrix S as simply scaling
factors that relate stress and strain without appreciat-
ing their role as operators that consider the material’s
directional response to loading.
• The stiffness matrix C and the compliance matrix S
are not simple scalar factors. They represent fourth-
order tensors (here represented in matrix form) that
fully describe the material’s response to multi-axial
stress states.

C: Relates the entire stress state to the entire strain

state, considering how each component of strain
influences each component of stress.
S: Similarly, S describes how each component of
stress influences the overall strain state.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 36 of 45

Which pair of material constants should I
use for Isotropic Materials?
The choice between expressing the constitutive equa-
tions in terms of Lamé constants λ and µ or terms of
Young’s modulus E and Poisson’s ratio ν depends on
the context of the problem and the specific advantages
and disadvantages associated with each set of param-
eters.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 37 of 45

Advantages of Using Lamé Constants λ and µ
1. Simplicity in Mathematical Formulation:

• Compact Expression: The constitutive equations

for isotropic materials are often more compact when
expressed using λ and µ. This is particularly true
in the context of continuum mechanics and elas-
ticity theory, where the equations involve the di-
vergence and gradient of displacement fields, and
relationships between these vector operators can
be established.
• Natural Fit in Elasticity Theory: λ and µ arise nat-
urally from the general theory of elasticity, where
they are directly linked to the volumetric and shear
deformation responses of the material. λ relates
to the material’s response to volumetric changes
(dilatation), while µ (the shear modulus) directly re-
lates to the material’s resistance to shear.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 38 of 45

2. Theoretical Insights:

• Connection to Stress-Strain Relationships: The use

of Lamé constants provides direct insight into how
a material responds to different types of deforma-
tion. For instance, µ is directly related to the ma-
terial’s shear modulus, which is crucial for under-
standing shear stress and strain. λ is associated
with the material’s ability to resist changes in vol-
ume under pressure.
• Utility in Advanced Theories: In advanced elastic-
ity theories and the study of seismic waves or ge-
omechanics, λ and µ are often more directly ap-
plicable than E and ν1, making them the preferred
choice in these contexts.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 39 of 45

Disadvantages of using Lamé Constants are they are
Less Intuitive for Practical Applications
1. Engineering Usage: E (Young’s modulus) and ν (Pois-
son’s ratio) are more commonly used and understood
in engineering practice. Engineers typically design
and specify materials using E and ν, so using λ and
µ may require additional conversion steps or might
be less intuitive.
2. Material Data Availability: Most material property database
and handbooks list materials in terms of E and ν. Us-
ing λ and µ requires conversion, adding an extra step
in practical engineering calculations.
Eν E
λ= , µ=
(1 + ν)(1 − 2ν) 2(1 + ν)

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 40 of 45

Advantages of using Young’s Modulus E and Pois-
son’s Ratio ν
1. Intuitive Understanding:

• Direct Physical Interpretation: E and ν have clear,

intuitive physical meanings. E describes how a
material stretches under tension, and ν describes
how it contracts laterally when stretched longitu-
dinally. This makes them more accessible to engi-
neers and more directly applicable to design and
analysis.
• Common Usage in Engineering: In engineering
applications, E and ν are the standard parameters
used to describe the elastic properties of materi-
als. They are often directly available in material
specifications, simplifying their use in calculations
and simulations.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 41 of 45

2. Standardization and Data Availability:

• Widely Available Data: Material properties are typ-

ically provided in terms of E and ν in engineer-
ing standards and databases. This makes it easier
for engineers to apply these values directly in their
work.
• Compatibility with FEA Software: Most finite el-
ement analysis (FEA) software packages are de-
signed to take E and ν as input parameters. This
compatibility simplifies the setup and execution of
simulations.

3. Ease of Application in Design:

• Direct Application to Design Formulas: Many en-

gineering design formulas, particularly those re-
lated to beams, shafts, and structural elements,
are derived using E and ν.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 42 of 45

Disadvantages of Using Young’s Modulus Poisson’s
Ratio
1. Less Fundamental in Theoretical Contexts.
2. Increased Complexity in Nonlinear and Anisotropic
Materials: Handling Complex Material Behaviors:
For materials that do not behave isotropically or lin-
early, the use of E and ν might be less convenient. In
such cases, λ and µ can provide a more flexible and
general framework for describing material behavior.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 43 of 45

Which to use in my application?
The choice between using Lamé constants λ and µ ver-
sus Young’s modulus E and Poisson’s ratio ν depends on
the context. E and ν are more intuitive and widely used
in engineering practice, making them ideal for design
and analysis in typical engineering applications. How-
ever, λ and µ offer advantages in theoretical elasticity,
complex material modeling, and certain advanced ap-
plications. Engineers and scientists should choose the
parameters that best fit the problem at hand, consid-
ering both the practical and theoretical aspects of the
material behavior they are analyzing. Either way, (λ, µ)
and (E, ν) can easily be converted back and forth.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 44 of 45

References
1. Hibbeler, R. C. (Any Edition). Mechanics of Materials,
Pearson.
2. Beer, F. P., Johnston, E. R., [Link]. (Any Edition). Me-
chanics of materials, McGraw-Hill.
3. Any textbook or chapter with the keyword Elasticity.
For example:

 Sadd, M. H. (2009). Elasticity: Theory, Applica-

tions, and Numerics. Academic Press.
 Timoshenko, S. P., & Goodier, J. N. (1970). The-
ory of Elasticity. McGraw-Hill.
 Boresi, A. P., & Schmidt, R. J. (2003). Advanced
Mechanics of Materials. Wiley.

LinkedIn � - Dr. Lonny Thompson, Clemson University, Sept. 1, 2024 45 of 45