Stress: Cauchy Sensor Tensor
- Identify that in a continuum, stress at a point is defined using a symmetric matrix, and not a number. Describe what the components of the matrix mean.
- Identify that the stress matrix is symmetric. Recognize the implications as having a coordinate system in which no shear stresses appear in the matrix. Employ the process of diagonalization of a symmetric matrix to find this coordinate system.
- Compute the normal and shear stresses on different planes using a stress matrix.
The topics presented here were first introduced by Augustin-Louis Cauchy in the nineteenth century. Cauchy introduced the idea of stress or traction vectors which are force vectors per unit area. He then introduced the stress as a symmetric matrix. The stress matrix is perhaps one of the early concepts that promoted the study of vectors and matrices. The word that is used to describe a matrix is perhaps due to the fact, as will be shown here, that the stress tensor acts on area vectors to produce traction vectors.
Before defining the “Stress Tensor/Matrix”, we first define the force vector per unit area. The force vector per unit area is termed the “Traction Vector”. The traction vector at a point inside a body can be defined after removing the adjacent material and considering that a force vector is exerted by the removed material on an area surrounding that point. The traction vector at this point considering the surface with the normal vector is then (Figure 1 ):
By invoking Newton’s third law of motion (every action has an equal and opposite reaction), the traction vector acting on the corresponding point on the removed material whose normal vector is -n is equal to (Figure 1 ):
The Cauchy Stress Tetrahedron
In this section we present the proof according to the French Mathematician Augustin-Louis Cauchy that shows that the state of stress at a particular point inside a continuum is well defined using a symmetric matrix, which is called the stress matrix or stress tensor. Knowing this matrix allows the calculation of any traction vector on any plane passing through that point. First, we consider an infinitesimal cube (Figure 2) at a point inside the material and we consider the three faces perpendicular to the three basis vectors , , and . The force per unit area vector on each of these three faces has three components and therefore, in total, there are nine components. We then set as the stress component on the face with normal i and acting in the direction . Notice that in some textbooks, this definition could be reversed. A tetrahedron can be formed by considering four faces perpendicular to , , , and an arbitrary vector (Figure 3). We are going to use the stress components to find the traction vector acting on the arbitrary area whose normal vector is . If is the area of the surface with normal , then:
where , , and are the areas of the triangular faces perpendicular to , , and respectively (Can you derive the above formula using the properties of the cross product?). Setting with :
Invoking Newton’s equations of equilibrium, the sum of the forces acting on the tetrahedron (including the gravity force where , , and are the components of the gravity acceleration) are equal to its mass multiplied by its acceleration vector . Therefore:
The mass can be related to the height of the tetrahedron perpendicular to using the mass density as follows:
Substituting for in the equilibrium equations, rearranging, and taking the limit as the volume of the tetrahedron goes to zero, (i.e., as ) lead to the linear relationship between the traction vector and the unit normal vector :
I.e., the Cauchy stress tensor is a linear operator that acts as a linear function from such that where is a unit vector, the result is the traction vector (force vector per unit area) acting on the surface with normal .
Symmetry of the Stress Matrix
Because of moment equilibrium whether the body is in static or dynamic equilibrium, it will be shown that in common materials, the Cauchy Stress Tensor is a symmetric tensor, i.e., . According to Euler’s laws of motion, the rate of change of the angular momentum in the material is equal to the moment of the external forces. Euler’s laws can be applied to an infinitesimal cuboid with external forces acting on its faces as shown in Figure 4. The moment around the axis perpendicular to the plane shown and passing through the bottom left corner (black dot), neglecting smaller quantities, can be evaluated as follows:
Assuming the cuboid rotates around the bottom left corner with angular velocity , then has the following form:
where is the density, is the angular acceleration and I is the polar moment of inertia. Equating with and taking the limit as the volume goes to zero(,):
Similarly, and . Therefore:
Normal and Shear Stress
Given a state of stress described by and given a plane with perpendicular , the normal stress is a real number representing the component of the force acting perpendicular to that plane and is given by (positive for tension and negative for compression):
While the shear stress is the magnitude of the component of the traction vector acting parallel to the surface of the plane:
Two Dimensional Illustrative Example
In this example a two dimensional state of stress at a point is described by the Cauchy stress tensor . Consider a plane with normal . The black arrow in the figure below indicates the direction of the plane, the traction vector is indicated by the blue arrow, the green arrow (outwards indicate tension and inwards indicate compression) shows the normal stress while the red arrow indicates the direction of the shear stress vector. The values are given next to the table. Try it out by changing the values in the boxes below and hitting enter (Notice that the vector n is normalized before any computations are performed):
Try to answer the following questions:
Find a state of stress and a direction when the green arrow is pointing opposite to the black arrow. What happens to the directions of the other vectors?
Find a state of stress and a direction such that the normal stress is zero while the shear stress is not and vice versa. What happens to the stress state when you change in each case?
What happens when you choose n to be horizontal or vertical?
What happens when you choose values for the stress components such that and ?
Three Dimensional Illustrative Example
Now try to answer the same questions on the three dimensional example below by changing the values in the box and hitting enter:
Principal Stresses and Principal Directions
Since is symmetric, there exists a coordinate system in which the component form of is diagonal, i.e., there is a coordinate transformation such that the planes perpendicular to the coordinate system have no shear stresses! This coordinate system is the one aligned with the eigenvectors of the stress tensor. The eigenvectors are called “The Principal Directions” while the eigenvalues are called “The Principal Stresses”.
The principal stresses are denoted , and . If , and are the corresponding normalized eigenvectors, then
is an orthogonal matrix of transformation into the coordinate system described by the eigenvectors. The order of the eigenvectors should be chosen to preserve the right hand orientation of the new coordinate system and in this case is a rotation. The stress matrix in this particular coordinate system will be diagonal:
Two Dimensional Illustrative Example
Consider a two dimensional state of stress at a point described by the Cauchy stress tensor . The square on the left in the figure below is aligned with the original coordinate system. The red and blue arrows represent the traction vectors acting on the horizontal and vertical planes, respectively. The example below finds the eigenvalues and eigenvectors of the stress matrix. Then, a coordinate transformation described by the matrix whose rows represent the eigenvectors of is used to rotate everything into that new coordinate system. In the new coordinate system, the force vectors are perpendicular to the faces of the square!
Three Dimensional Illustrative Example
Similarly, the example below illustrates the concept of rotating a cube into the coordinate system described by the eigenvectors of . The sizes of the arrows showing the traction or principal stress vectors qualitatively represent the magnitudes of the vectors. I.e they are not proportionate to the magnitude but only represent which vector has the larger or smaller magnitude relative to others.
Mohr’s Circle for Plane Stress
Visit Wikipedia’s entry on Mohr Circle to learn about the history, the construction and the applications of Mohr’s circle!
In the following example, you can change the values of the stress matrix entries from -10 to 10 units to reconstruct Mohr’s circle. The black arrow indicates the point on the circle corresponding to a plane oriented at an angle (theta) with respect to the dotted line which corresponds to the plane perpendicular to the horizontal direction. The angle on the circle between the dotted and solid arrows is (2theta).