## Displacement and Strain: Strain Measures

In this section, we give various definitions and interpretations of strain measures. We will start the section with the definition of uniaxial strains:

### Uniaxial Strain Measures:

#### Definitions:

Assume that a bar of length is uniformly stretched to a length . The quantity is termed the stretch ratio. There are various uniaxial strain measures that we can use to define how the length of the bar changes:

The Engineering Strain is defined as:

The True (Logarithmic) Strain is defined as:

The Green (Lagrangian) Strain is defined as:

The stretch ratio has the limits , which limits the strain measures as follows:

The above three strain measures are related as follows:

There are advantages to each of the above strain measures. The engineering strain is the most intuitive measure as it is simply the ratio of the change in length over the original length. The true strain is additive which led to the terminology: “True”! For example, if a bar is stretch from to and then to . Then, the true strain measured from to is also equal to the sum of the true strain from to plus the true strain from to .

The Green strain, on the other hand, can easily be computed in a three dimensional object with nonlinear deformation.

There is a very subtle difference between the engineering and the true strain measures when viewed as integrals. The true strain is the integral of the function while the engineering strain is the integral of the function :

The following tool draws the variation of the above three strain measures when the stretch ratio varies between 0 and 2. By changing the slider you can also pick a value for the stretch ratio between 0 and 2 and the tool will display the corresponding values of the strain measures. Notice that the three strain measures agree when the stretch ratio is close to 1, i.e., for small strains, the three measures agree.

#### Nonuniform Axial Strain:

In the above discussion, we assumed that the bar deforms uniformly. I.e., the extension or contraction is the same at every point on the bar. However, a more general case is when the bar deforms non-uniformly. In this case, we can assume that the original position of the bar is given by a function where varies between and . The new position is given by . We can think of the local stretch as a the limit:

In this case, the three strain measures would be given by:

##### Example:

Assume a bar of length with the ends situated at and . Assuming that the new position , then, the stretch is different at every point and is given by:

I.e., the stretch is zero at the end , while is the highest at the end . In this case, the strain is not uniformly distributed, but rather is given by the functions:

The following tool draws the sketches of the reference and deformed configurations. The vertical lines in the reference configuration are equidistant but they are not after deformation. Change the density of the vertical lines and notice that after deformation, the distances between the vertical lines on the right side of the bar are more than the distances on the left side.

The following tool draws the distribution of the strain as a function of the original position . As predicted, the strain on the left side of the bar after deformation, is higher than the strain on the right side. You can change the slider to pick a specific value for for which the tool will calculate and display the value of the different measures of strain.

### Three Dimensional Strain Measures:

In general, strain measures give a measure of how the lengths or angles change when we compare vectors in the undeformed configuration with vectors in the deformed configuration. If a tangent vector in the reference configuration is denoted , while its deformed image is denoted , then we can calculate the various strain measures along the direction by the simple computations:

where we can use the relationship between the dot product and the norm when convenient:

However, a more convenient and traditional way to measure strain is to define strain tensors by which we can calculate strain along certain directions. We will discuss the two most common strain tensors in the literature:

#### Small (Infinitesimal) Strain Tensor:

The small of infinitesimal strain tensor is defined as the symmetric part of the displacement gradient :

Which has the following component form:

which can be written in a simple form as follows :

In the case of **small deformations**, the small strain tensor can be used to compute the engineering longitudinal and shear strains as shown below.

##### Calculating Longitudinal Strains along General Vectors:

In the case of **small deformations**, the engineering longitudinal strain along a tangent vector in the reference configuration can be calculated using the relationship:

To show the above relationship, we are going to use the relationship between the deformed tangent vector and its original vector as shown in the displacement gradient tensor section:

The length of can be estimated as follows:

Since is a skewsymmetric tensor, we have: . For small deformations, the following terms can be neglected , and which simplifies the above relationship to:

Since the term , the above relationship can be further simplified as follows:

Thus, for ** small deformations**, the engineering strain along the reference configuration tangent vector can be calculated as follows:

##### Calculating Longitudinal Strains Along the Basis Vectors:

The strain along the basis vectors , and in the reference configuration can be calculated using the above relationship as follows :

Therefore,

Therefore, the diagonal components of the strain matrix give the value of the longitudinal strains along the basis vectors of the reference configuration.

##### Calculating Angle Change between General Vectors:

Given two vectors and separated by an angle in the reference configuration that deform into the two vectors and separated by an angle in the deformed configuration, then half the change in the dot produce between the two vectors before and after deformation can be calculated as:

Using the relationship between the deformed tangent vectors , and their original vectors , as shown in the displacement gradient tensor section we can write and as functions of and :

Since is symmetric and is skewsymmetric and for **small deformations**, the dot product can be approximated as follows:

Then half the change in the dot produce between the two vectors before and after deformation can be calculated using as follows:

When the deformations are small such that and then half the difference in the cosines of the angles before and after deformation can be calculated as follows:

If the original vectors and are orthogonal to each other, then . In that case, the above relationship gives half the engineering shear strain of planes parallel to and perpendicular to .

##### Calculating Engineering Shear Strain in the Planes of the Basis Vectors:

Using the above relationship, the engineering shear strains in the planes of the basis vectors and with can be calculated as follows:

i.e., the off diagonal components give the engineering shear strains in the planes –, – and –

##### Example:

In the following example, a two dimensional square centred around the origin is deformed using the displacement relationship:

The following tool lets you vary the values of and between -0.8 and 1.2 and the values of and between -0.2 and 0.2. The tool then draws the square before deformation and then after deformation. You can also select two unit vectors and by varying their respective angles of inclination with the horizontal axis and . The two vectors are drawn before and after deformation. is drawn in blue and is drawn in red. The tool calculates the deformation gradient , the gradient of displacement tensor , the small strain matrix and the infinitesimal rotation matrix . The longitudinal strains along and are calculated underneath the figures in addition to where is the angle between and which are the respective images of and . Vary the values of to see their effect on the deformation and the components of the strain matrix . Also, check your strain computations against the results shown for various choices of and . Note that the computation of is approximate because it assumes small deformation and that

##### Symmetry of the Strain Tensor:

The symmetry of the infinitesimal strain matrix implies that there is a coordinate system in which the strain matrix is diagonal. The eigenvalues of the strain matrix are called the principal strains. If an orthonormal coordinate system aligned with the eigenvectors is chosen for the basis set, then the strain matrix is diagonal (no shear strains). In other words, in this coordinate system, no shear strains exist on the planes perpendicular to the basis vectors.

##### Example:

The example shown above is repeated. In the first row, the tool draws the original shape and the deformed shape showing the deformation of the original coordinate system. In the second row, the tool draws the original shape of a small square aligned with the eigenvectors of the strain matrix which are drawn in blue and red. Notice that the square deforms into a rectangle (i.e., no shear strain!). The new rectangle might be slightly rotated as the deformation in this case is decomposed into an infinitesimal stretch described by and an infinitesimal rotation described by . You can vary the components to investigate their effect on the deformed shape. Which components control the stretch and which components control the rotation?

##### Calculating Volumetric Strains:

For **small deformations**, the diagonal components of the strain matrix can be used to calculate the volumetric strain as follows:

where and are the deformed and original volumes of a small element inside the deformed object.

##### Behaviour under Pure Rotations:

A major disadvantage of the infinitesimal strain matrix is that it predicts strains for bodies that undergo large rotations even when the actual strain is zero. For example, consider a two dimensional displacement described by the rotation matrix:

i.e.,

Then, the displacement vector at every point is described by:

The gradient of displacement matrix is given by:

The infinitesimal strain matrix is then given by:

For small values of , . However, at larger values of rotation, the strain matrix will predict high strains even though the object rotating is not being stretched. The value for the compressive strain in the direction of the basis vectors and should not come as a surprise. If a square of unit length rotates with an angle , is actually the length of the projections of the square’s sides onto the coordinate axis!

#### Green Strain Tensor:

The Green strain tensor is is defined as follows:

By utilizing the relationship , we can replace to obtain the relationship:

In component form, the Green strain tensor components can be written as:

Similar to the small strain tensor, the Green strain tensor is symmetric, i.e., there is a coordinate system in which the strain matrix has diagonal components.

##### Calculating Longitudinal Green Strains along General Vectors:

The Green strain tensor can be used to calculate the uniaxial (longitudinal) Green strain along general vectors. Let be a vector in the reference configuration and be its deformed image, then the Green strain along can be calculated as follows:

Utilizing the relationship we get:

Therefore:

##### Calculating Change of the Dot Product between General Vectors:

Similarly, given two vectors and before deformation and their respective images and after deformation, **half** the difference between the dot product in the deformed configuration and the dot product in the reference configuration can be calculated using the Green strain tensor as follows:

##### Behaviour under Pure Rotations:

One of the major advantages of the Green strain tensor is that it predicts zero strain in the cases of pure rotation. Assuming that a deformation is described by a rotation matrix , i.e.: . Then, the Green strain tensor is calculated as follows:

i.e., under pure rotations, the Green strain matrix has zero components!

##### Eigenvectors:

The eigenvectors of the Green strain matrix stay perpendicular to each other after deformation. This can be shown as follows. Let and be two distinct eigenvectors of . Since is symmetric, then . After deformation we have:

However, since is an eigenvector of , then there exists a real number such that . Therefore:

Since we have:

This is not necessarily true for the eigenvectors of the small strain matrix. You can find a counter example for the small strain matrix using the following tool!

##### Example:

The above example is repeated but with higher ranges of the components to account for larger deformations. In the first row, the square and its deformed shape are shown. In the second row, a square aligned with the eigenvectors of the small strain matrix is drawn along with its image. In the third row, a square aligned with the eigenvectors of the Green strain matrix is drawn along with its image. Vary the components and try to answer the following:

- At what values of do the Green strain matrix and the small strain matrix have similar or different values?
- What happens to the eigenvectors of the small strain matrix at higher rotations? Are they still perpendicular to each other after deformation?
- What happens to the eigenvectors of the Green strain matrix at higher rotations? Are they still perpendicular to each other after deformation?

### Examples and Problems:

#### Example 1:

Calculate the infinitesimal and Green strain matrices for the following position function:

Also, find the displacement function, the uniaxial small and Green strains along the direction of the vector .

##### Solution:

The displacement function is:

The small strain tensor is:

The Green strain tensor is:

The deformation is very large as shown by applying this deformation to a unit cube (see figure below), so the strain measures are different.

The uniaxial small and Green strain along the vector can be obtained as follows:

View Mathematica Code

F={{1.2,0.2,0.2},{0.2,1.3,0.1},{0.9,0.5,1}}; X={X1,X2,X3}; x=F.X; u=x-X; Gradu=Table[D[u[[i]],X[[j]]],{i,1,3},{j,1,3}] einfinitesimal=1/2*(Gradu+Transpose[Gradu]); %//MatrixForm egreen=1/2*(Gradu+Transpose[Gradu]+Transpose[Gradu].Gradu); %//MatrixForm u={1,1,1}; n=u/Norm[u]; n.einfinitesimal.n n.egreen.n Box1=Cuboid[{0,0,0},{1,1,1}]; Box2=GeometricTransformation[Box1,F]; Graphics3D[{EdgeForm[Thickness[0.01]],Box1,EdgeForm[Thickness[0.005]],Box2},Axes->True,AxesOrigin->{0,0,0},BaseStyle->Directive[Bold,15],AxesLabel->{e1,e2,e3}]

#### Example 2:

A cube of unit dimensions undergoes a rigid body rotation in of , then , and then around the basis vectors , , and respectively. Use Mathematica to visualize the box. Compare the two strain measures and .

##### Solution:

The deformation gradient is given by:

Note that is on the right since it is applied first.

The Green strain is identically zero while the small strain tensor predicts the following strains:

View Mathematica Code

Qx=RotationMatrix[thx,{1,0,0}]; Qy=RotationMatrix[thy,{0,1,0}]; Qz=RotationMatrix[thz,{0,0,1}]; thx=30Degree thy=20Degree thz=20Degree F=Qx.Qy.Qz; X={X1,X2,X3}; x=F.X; u=x-X; Gradu=Table[D[u[[i]],X[[j]]],{i,1,3},{j,1,3}]; einfinitesimal=1/2*(Gradu+Transpose[Gradu]); einfinitesimal=FullSimplify[einfinitesimal]; %//MatrixForm egreen=1/2*(Gradu+Transpose[Gradu]+Transpose[Gradu].Gradu); egreen=FullSimplify[egreen]; %//MatrixForm Box1=Cuboid[{0,0,0},{1,1,1}]; Box2=GeometricTransformation[Box1,F]; Graphics3D[{Box1,Box2},Axes->True,AxesOrigin->{0,0,},AxesLabel->{e1,e2,e3}]

#### Example 3:

A cube of unit length undergoes a simple shearing motion of an angle in the direction of and perpendicular to . Find the infinitesimal and Green strain measures that describe this motion. Comment on the difference between the two strain measures.

##### Solution:

The position function is given by:

The deformation gradient is given by:

The small strain and Green strain tensors are given by:

Both measures predict the same shear strains. For small deformations, both measures predict zero normal strains. As the increases, the vertical lines start developing large axial strains which the can predict while cannot. The undeformed and deformed shapes of the cube are shown below for .

F={{1,Tan[theta],0},{0,1,0},{0,0,1}}; X={X1,X2,X3}; x=F.X; u=x-X; Gradu=Table[D[u[[i]],X[[j]]],{i,1,3},{j,1,3}] einfinitesimal=1/2*(Gradu+Transpose[Gradu]); einfinitesimal=FullSimplify[einfinitesimal]; %//MatrixForm egreen= 1/2*(Gradu+Transpose[Gradu]+Transpose[Gradu].Gradu); egreen=FullSimplify[egreen]; %//MatrixForm Box1=Cuboid[{0,0,0},{1,1,1}]; Box2=GeometricTransformation[Box1,F]; theta=30Degree; Graphics3D[{Box1,Box2},Axes->True,AxesOrigin->{0,0,0},AxesLabel->{e1,e2,e3}]

#### Example 4:

The following small strain matrix represents a two dimensional state of strain

Find the principal strains and their directions. If a coordinate system is aligned with the principal directions, find the components of the strain matrix in this new coordinate system.

##### Solution:

The principal strains and and their directions , and be obtained using Mathematica by finding the eigenvalues and eigenvectors of . The following are the eiganvalues:

The following are the normalized eigenvectors.

The transformation matrix between the original coordinate system and the new coordinate system is:

The strain matrix in the new coordinate system has the form:

View Mathematica Code

eps={{0.01,0.012},{0.012,0}}; s=Eigensystem[eps] Q={-s[[2,1]],-s[[2,2]]} epsdash=Chop[Q.eps.Transpose[Q]]

#### Example 5:

A state of uniform small strain is described by the following small strain matrix:

Determine the following:

- The longitudinal strain along the direction of the vector .
- The change in the cosines of the angles between the vectors and and comment on whether the angle between the vectors decreased or increased after deformation.
- The shear strains of planes parallel to the vector and perpendicular to the vector .

##### Solution:

The strains along can be calculated as follows:

The change in the cosines of the angles between the vectors and can be calculated as follows:

The angle between the vectors and before deformation can be calculated using the dot product:

Therefore, the angle after deformation between the vectors is larger because:

Since the vectors and are perpendicular, the shear strains of planes parallel to the vector and perpendicular to the vector can be obtained as follows:

View Mathematica Code

Eps={{-0.01,0.002,-0.023},{0.002,0.05,0},{-0.023,0,0.02}}; dX={1,2,1}; Norm[dX] Eps.dX dX.Eps.dX/Norm[dX]^2 dX={1,0,1}; dY={1,1,1}; anglechange=dX.Eps.dY/Norm[dX]/Norm[dY] THETA=N[ArcCos[dX.dY/Norm[dX]/Norm[dY]]] THETA/Degree theta=ArcCos[2*anglechange+Cos[THETA]] theta/Degree dX={-1,1,0}; dY={1,1,0}; ShearSTrain=dX.Eps.dY/Norm[dX]/Norm[dY]

#### Example 6:

The longitudinal engineering strains on the surface of a test specimen were measured using a strain gauge rosette to be 0.005, 0.002 and –0.001 along the three directions: , , and respectively, where , , and . If the material is assumed to be in a small strain state, find the principal strains and their directions on the surface of the test specimen.

##### Solution:

The state of strain on the surface of the material can be represented by the symmetric small strain matrix:

Since the longitudinal strain is know along three directions, three equations can be written to find the three unknown components of the strain matrix as follows:

The norm of each of , , and is 1. Therefore, the three equations are:

Solving the above three equations yields:

The principal strains and their directions are the eigenvalues and eigenvectors of the small strain matrix. The principal strains are:

The principal directions are:

View Mathematica Code

eps=Table[Subscript[e,i,j],{i,1,2},{j,1,2}] a={Cos[30Degree],Sin[30Degree]}; b={0,1}; c={-Cos[30Degree],Sin[30Degree]}; FullSimplify[a.eps.a/Norm[a]^2] FullSimplify[b.eps.b/Norm[b]^2] FullSimplify[c.eps.c/Norm[c]^2] a=Solve[{a.eps.a/Norm[a]^2==0.005,b.eps.b/Norm[b]^2==0.002,c.eps.c/Norm[c]^2==-0.001,eps[[1,2]]==eps[[2,1]]},Flatten[eps]] eps=eps/.a[[1]]; Eigensystem[eps]//MatrixForm

#### Example 7:

In the above examples, the deformation was described by a uniform strain matrix (i.e., a strain matrix that is not function in position). In this example, a nonuniform strain, i.e., the strain matrix is function of position is explored. Assume that the two dimensional displacement function of a 2units by 2units plate has the following form:

- Find the deformed configuration function .
- Find and .
- Draw the vector plot of the displacement function on the plate.
- Draw the contour plots of , , , and on the plate.

##### Solution:

The position function of the deformed configuration is given by:

The deformation gradient is given by:

The displacement gradient is given by:

Therefore, the small strain matrix is given by:

The Green strain matrix is given by:

The required plots are shown below.

View Mathematica Code

Clear[x1,x2,X1,X2] X={X1,X2} ; u={0.2X1,0.09X2*X1} ; x=u+X; Gradu=Table[D[u[[i]],X[[j]]],{i,1,2},{j,1,2}] einfinitesimal=1/2*(Gradu+Transpose[Gradu]); egreen=1/2*(Gradu+Transpose[Gradu]+Transpose[Gradu].Gradu); einfinitesimal//MatrixForm FullSimplify[egreen]//MatrixForm VectorPlot[u,{X1,0,2},{X2,0,2},BaseStyle->Directive[Bold,15],AspectRatio->Automatic,PlotLabel->"u"] ContourPlot[einfinitesimal[[1,1]],{X1,0,2},{X2,0,2} ,BaseStyle->Directive[Bold,15],AspectRatio->Automatic,ContourLabels->All,PlotLabel->"einf11"] ContourPlot[egreen[[1,1]],{X1,0,2},{X2,0,2} ,BaseStyle->Directive[Bold,15],AspectRatio->Automatic,ContourLabels->All,PlotLabel->"egreen11"] ContourPlot[einfinitesimal[[1,2]],{X1,0,2},{X2,0,2} ,BaseStyle->Directive[Bold,15],AspectRatio->Automatic,ContourLabels->All,PlotLabel->"einf12"] ContourPlot[egreen[[1,2]],{X1,0,2},{X2,0,2} ,BaseStyle->Directive[Bold,15],AspectRatio->Automatic,ContourLabels->All,PlotLabel->"egreen12"] ContourPlot[einfinitesimal[[2,2]],{X1,0,2},{X2,0,2} ,BaseStyle->Directive[Bold,15],AspectRatio->Automatic,ContourLabels->All,PlotLabel->"einf22"] ContourPlot[egreen[[2,2]],{X1,0,2},{X2,0,2} ,BaseStyle->Directive[Bold,15],AspectRatio->Automatic,ContourLabels->All,PlotLabel->"egreen22"]

#### Problems:

- Which of the following functions describing motion is/are not physically possible and why?
- .
- .
- .
- .

- Let a cube of unit length be represented by the set of vectors such that three of the cube’s sides are aligned with the three orthonormal basis set vectors , , and and one of the cube vertices lies on the origin of the coordinate system. Let be the deformed configuration such that is defined as :
Find the following:

- The displacement function .
- The new position of the 8 vertices of the cube.
- The deformed curves of the three edges of the cube that are aligned with the basis vectors , , and .
- The infinitesimal strain matrix and the volumetric strain as a function of the position inside the cube.
- The Green strain matrix as a function of the position inside the cube.
- Evaluate the infinitesimal strain matrix at two of the eight cube vertices.
- Evaluate the Green strain matrix at two of the eight cube vertices.

- A cube undergoes a deformation such that the infinitesimal strain is described by the matrix:
Find the following:

- The principal strains and their directions.
- The longitudinal strains along the direction of the vector:
- The approximate change in the angles between the vectors and
- The coordinate system transformation in which the strain matrix is diagonal. Express the strain matrix in the new coordinate system.

- The longitudinal engineering strains on the surface of a test specimen were measured using a strain rosette to be 0.005, 0.002, and –0.001 along the three directions: , , and where is along the direction of the basis vector , is oriented anticlockwise from , and is oriented anticlockwise from . If the material is assumed to be in a small strain state, find the principal strains and their directions on the surface of the test specimen.
- The measured strains for the three axes of a strain gauge rosette are: , , and . If coincides with the basis vector while and coincide with vectors that are anticlockwise and respectively from . Find the components of the strain tensor.
- Assume that a two dimensional rectangular plate with dimensions 5 units along the basis vector and 1 unit along the basis vector is situated such that the origin of the coordinate system is at the midpoint of the left side of the plate. If the displacement function of the plate is described by:
find the following:

- The position function .
- The infinitesimal strain matrix as a function of the position.
- The Green strain matrix as a function of the position.
- Draw the vector plot of the displacement function.
- Draw the contour plots of , , , and .

- The general equation for the coordinates of a point on an ellipsoidal surface whose axes are oriented along the Cartesian coordinate system is:
where , , and are the three semi-diameters of the ellipsoidal surface. When , the surface is a sphere with radius . Assuming that the motion of an object is described by the homogeneous deformation gradient such that the diagonal components of are the only non-zero components.

- Show that a sphere of radius deforms into an ellipsoid surface.
- Find the relationship between the semi-diameters of the ellipsoid and the components of the deformation gradient.
- Find the relationship between the semi-diameters of the ellipsoid and the components of the infinitesimal strain tensor.

- Three unit square elements in deform as shown in the figure below.

For each case:- Write expressions for the displacement components and as functions of , , and the variables shown for each square.
- Determine the components of the infinitesimal strain tensor and the infinitesimal rotation tensor.
- Determine the components of the Green strain tensor.

hey Samer, for Problem 2 should x belong to omega and X belongs to omega0?

Yes. You are absolutely correct. Thank you! I have corrected it!

hey Samer, for calculating angles between vectors section, first formula, right side of = sign, should it be ||dX||||dY||cos, instead of ||dX||||dX||?

Hey Samer, in 4.5.2 where you ffirst introduced green strain you used “epsilon Green along dx” direction, laterin green strain tensor section you used “epsilon Green along dX” direction. Are these two supposed to be interchangeable?

Thank you for noticing this.

For consistency, I changed the strain measures to read “The strain along the direction dX”. dX is a direction vector in the reference configuration whose image is dx. The strain is a measure of how the length of dx is different from the length of dX.