## Energy: Expressions for Linear Elastic Strain Energy Functions

### Strain Energy in Linear Elastic Materials:

In the previous section we ended up showing that the increment in the energy per unit volume of the undeformed configuration is equal to:

(1)

As described in the section on hyperelastic materials, elasticity of a material implies that there exists a strain energy function from which the stresses can be derived. For a function to exist, then Equation 1 has to be an exact differential implying that can be obtained by differentiating a function with respect to (See hyperelastic materials).

In the derivation of 1 we used the following definition:

For linear elastic materials in small deformations, we can assume that . Additionally, we can assume that the derivatives with respect to are equal to the derivatives with respect to and thus, we can use the engineering strain tensor as defined in the section on strain measures so that:

Additionally, we can assume that the stress components are linear functions of the strain components and are obtained by differentiating with respect to the strain components as shown in the linear elastic constitutive laws section. In this case, the strain energy stored during deformation can be calculated as follows:

(2)

where the factor is due to the linear relationship between the stress and the strain.

#### Strain Energy in Linear Elastic Isotropic Materials:

If we adopt the constitutive relationship for isotropic linear elastic materials, then we can replace the stresses as functions of the strains or replace the strains as functions of the stresses. We can then write the following two equivalent forms of Equation 2 as applied to isotropic linear elastic materials:

(3)

(4)

##### Volumetric and Deviatoric Strain Energies in Linear Elastic Isotropic Materials:

The strain energy of a linear elastic isotropic material per unit volume can be divided into two additive components: a volumetric strain energy function that accompanies volumetric changes and a deviatoric strain energy function that accompanies shape changes at constant volume. Such decomposition is valid for small strains when the volumetric strain can be measured as the trace of the infinitesimal strain matrix. In this case, the relationship between the volumetric strain and the applied stresses using the linear elastic constitutive laws is:

The above relationship shows that nonzero volumetric strains exist only when the hydrostatic stress is nonzero. To separate the volumetric from the deviatoric strain energy functions, the deviatoric stress will be used. By replacing with and in Equation 2 the following is obtained:

(5)

The right hand side can be separated into two terms:

The first term is called the deviatoric strain energy term while the second is called the volumetric strain energy term. Using the linear elastic constitutive laws, these terms can be further simplified to have the following form:

with

The reader should attempt to manipulate Equation 5 to obtain the above simplified version in which the deviatoric strain energy is dependent solely on the von Mises stress , while the volumetric strain energy is dependent solely on the hydrostatic stress .

#### Linear Elastic Strain Energy in bars under Axial Loading with Small Deformations:

Consider a linear elastic bar of length such that its neutral axis before deformation is aligned with the axis. Then, using the expressions for the stresses and strains in an a plane bar under axial loading, the total strain energy stored in the bar can be evaluated using the following integral:

#### Linear Elastic Strain Energy in Euler Bernoulli Beams with Small Deformations:

Consider a linear elastic Euler Bernoulli beam with a constant cross-sectional area and with a length such that its neutral axis before deformation is aligned with the axis. Then, using the expressions for the stresses and strains in an Eulber Bernoulli beam, the total strain energy stored in the beam can be evaluated using one of the forms of the following integral:

#### Linear Elastic Strain Energy in Timoshenko Beams with Small Deformations:

Consider a linear elastic Euler Tiomoshenko beam with a constant cross-sectional area and with a length such that its neutral axis before deformation is aligned with the axis. Then, using the expressions for the stresses and strains in a Timoshenko beam, the total strain energy stored in the beam can be evaluated using one of the forms of the following integral:

### Examples and Problems:

#### Example 1:

The strains in a linear elastic orthotropic material are given by the strain matrix:

If the material axes of orthotropy are aligned with the orthonormal basis set used, then find the stresses and the strain energy per unit volume if , , , and

##### Solution:

The stress strain relationships for orthotropic materials shown can be used to find the stresses. Therefore, the stress matrix has the following form:

The strain energy per unit volume for a general linear elastic material is given by Equation 2:

View Mathematica Code

E11=100; E22=100; E33=200; G12=G13=G23=35; eps={{0.01,0.015,0.001},{0.015,-0.02,0.002},{0.001,0.002,0.001}}; stress=Table[0,{i,1,3},{j,1,3}]; stress[[1,1]]=eps[[1,1]]*E11; stress[[2,2]]=eps[[2,2]]*E22; stress[[3,3]]=eps[[3,3]]*E33; stress[[1,2]]=stress[[2,1]]=2eps[[1,2]]*G12; stress[[1,3]]=stress[[3,1]]=2eps[[1,3]]*G13; stress[[2,3]]=stress[[3,2]]=2eps[[2,3]]*G23; stress//MatrixForm eps//MatrixForm Energy=Sum[stress[[i,j]]*eps[[i,j]]/2,{i,1,3},{j,1,3}]

#### Example 2:

If the state of stress in a linear elastic isotropic material is given by the stress matrix:

Find the strain energy density and its deviatoric and volumetric components stored inside the material if and .

##### Solution:

Using Equation 3, the strain energy density is:

The deviatoric and volumetric components are:

View Mathematica Code

s = {{20, 30, 0}, {30, -10, 0}, {0, 0, 25}}; Nu = 0.1; Ee = 10000; G = Ee/2/(1 + Nu); U = (1 + Nu)/2/Ee*(Sum[s[[i, j]]^2, {i, 1, 3}, {j, 1, 3}]) - Nu/2/Ee*(Sum[s[[i, i]], {i, 1, 3}])^2 p = (s[[1, 1]] + s[[2, 2]] + s[[3, 3]] )/ 3 Udeviatoric = 1/12/G*((s[[1, 1]] - s[[2, 2]])^2 + (s[[3, 3]] - s[[2, 2]])^2 + (s[[1, 1]] - s[[3, 3]])^2 + 6*(s[[1, 2]]^2 + s[[1, 3]]^2 + s[[2, 3]]^2)) Uvolumetric = 3 (1 - 2 Nu)/2/Ee*p^2

#### Example 3:

The position function in dimensions of m. of a plate follows the following form:

The plate has dimensions of 2m in the direction of and 1m in the direction of with the bottom left corner coinciding with the origin of the coordinate system. Assuming that the material is linear elastic isotropic with Young’s modulus of 210 GPa and Poisson’s ratio of 0.3, then:

- Find an expression for the strain energy density as a function of the position inside the plate, and draw the contour plot of the strain energy density function on the plate.
- If the thickness of the plate is 10mm, find the total strain energy stored inside the plate.
- Find an expression for the deviatoric and volumetric strain energy density functions, and draw the contour plot of each on the plate.

##### Solution:

The displacement function can be calculated as follows:

The gradient of the displacement is given by the following matrix:

The engineering small strain is given by:

The vector representation of the strain is given by:

Using the constitutive equation for linear elastic isotropic materials, the vector representation of the stress is given by:

The strain energy density is given by:

The total energy is given by:

The deviatoric and volumetric strain energy densities are given by:

The required contour plots are:

Clear[X1,X2,X3,x1,x2,x3,Ee,Nu] X={X1,X2,X3}; x={x1,x2,x3}; x1=X1+0.001X2+0.0002X1^2; x2=1.001X2-0.0002X2^2; x3=X3; u=x-X; gradu=Table[D[u[[i]],X[[j]]],{i,1,3},{j,1,3}]; esmall=1/2(gradu+Transpose[gradu]); esmall//MatrixForm Strainvector={esmall[[1,1]],esmall[[2,2]],esmall[[3,3]],2*esmall[[1,2]],2*esmall[[1,3]],2*esmall[[2,3]]}; Strainvector//MatrixForm Ee=210000; Nu=0.3; G=Ee/2/(1+Nu); Cc={{1/Ee,-Nu/Ee,-Nu/Ee,0,0,0},{-Nu/Ee,1/Ee,-Nu/Ee,0,0,0},{-Nu/Ee,-Nu/Ee,1/Ee,0,0,0},{0,0,0,1/G,0,0},{0,0,0,0,1/G,0},{0,0,0,0,0,1/G}}; Dd=FullSimplify[Inverse[Cc]]; stressvector=FullSimplify[Chop[Dd.Strainvector]]; stressvector//MatrixForm S={{stressvector[[1]],stressvector[[4]],stressvector[[5]]},{stressvector[[4]],stressvector[[2]],stressvector[[6]]},{stressvector[[5]],stressvector[[6]],stressvector[[3]]}}; S=Chop[FullSimplify[S]]; S//MatrixForm StrainEnergy=FullSimplify[Sum[1/2*S[[i,j]]*esmall[[i,j]],{i,1,3},{j,1,3}]] TotalEnergy=Integrate[StrainEnergy,{X1,0,2},{X2,0,1},{X3,0,0.01}] Udeviatoric=FullSimplify[1/12/G*((S[[1,1]]-S[[2,2]])^2+(S[[3,3]]-S[[2,2]])^2+(S[[1,1]]-S[[3,3]])^2+6*(S[[1,2]]^2+S[[1,3]]^2+S[[2,3]]^2))] Uvolumetric=FullSimplify[(1-2Nu)/6/Ee*(S[[1,1]]+S[[2,2]]+S[[3,3]])^2] (*Check that the above expressions for the energy add up*) Chop[FullSimplify[StrainEnergy-Udeviatoric-Uvolumetric]] ContourPlot[StrainEnergy,{X1,0,2},{X2,0,1},AspectRatio->Automatic,ContourLabels->All,PlotLabel->"StrainEnergy MN m/m^3"] ContourPlot[Udeviatoric,{X1,0,2},{X2,0,1},AspectRatio->Automatic,ContourLabels->All,PlotLabel->"Deviatoric StrainEnergy MN m/m^3"] ContourPlot[Uvolumetric,{X1,0,2},{X2,0,1},AspectRatio->Automatic,ContourLabels->All,PlotLabel->"Volumetric StrainEnergy MN m/m^3"]

#### Example 4:

Consider a linear elastic small deformations cantilever beam. If a load of value is applied at the free end, find an expression for the strain energy per unit volume as a function of the position, , and the beam length . Also find the total strain energy stored. Assume Young’s modulus to be and that the beam has a rectangular cross section with a moment of inertia .

##### Solution:

The bending moment in a cantilever beam loaded by a concentrated load at the free end is given by:

The stress component and the strain component are given by:

Therefore, the strain energy per unit volume at any point is given by:

The total strain energy stored in the beam is given by:

Note that the displacement at the free end is equal to:

By assuming a gradual increase of the external load , the work done during the application of the load is equal to:

Therefore, the work done on the beam is equal to the total strain energy stored in the beam.

#### Problems:

- A linear elastic isotropic material with a shear modulus and a bulk modulus is stretched such that the small strain matrix is:
Find the strain energy density stored during the application of the deformation, and find its deviatoric and volumetric components.(Answer: 7.58, 3.38, 4.21 kN.m/m^3)

- Two states of stress are given by the stress matrices:
where , , and is the identity matrix. Find the deviatoric and the volumetric strain energy densities associated with both matrices if the material is linear elastic isotropic with Young’s modulus and Poisson’s ratio .

- Compare the total strain energy stored in a simple versus fixed ends beam of length , Young’s modulus , moment of inertia loaded with a constant distributed load . (Answer: Simple beam has 6 times more energy stored for the same distributed loading!) .
- For a linear isotropic elastic material with a Young’s modulus and Poisson’s ratio , express the strain energy density function under the condition of plane stress in terms of:
- , and .
- , and .

- For a linear isotropic elastic material with a Young’s modulus and Poisson’s ratio , express the strain energy density function under the condition of plane strain in terms of:
- , and .
- , and .

Hi! I was just wondering that in example 4, the total strain energy of the beam is P^2*L^3 / 6EI and the displacement in the free end is PL^3 / 3EI and therefore Energy = P*Delta / 2. If we consider that the load P is a mass hanging from the end of the beam, shouldn’t the potential energy change of the mass equal the work done? The potential energy is calculated as E_pot = mgh which is with the symbols used above E_pot = Delta * P. This is two times larger than the work.

In statics applications, the load is assumed to increase gradually when evaluating the internal energy of the beam. If a load is applied suddenly when the displacement is zero, then the beam would oscillate around the equilibrium position and would have both potential and kinetic energy at every point in time. In the “dynamics” case where P*D is used (i.e., mgh) for the potential energy, the beam would never be in a state of static equilibrium.