In hypoelasticity, we took the Cauchy stress, , and developed a constitutive relationship with the Left-Cauchy Green Strain Tensor, . This relationship was quite complicated (eq. 9 in Section 6: Hyperelasticity). Recall also that these tensors were both defined in spatial coordinates, by definition. Then, we derived an even more complicated “hypoelastic” formula (eq. 7 in Section 7: Hypoelasticity) that involved the rate of deformation tensor, . Lastly, a simple three-line pseudo-algorithm was given to illustrate how actual FEA (“solid” element) codes handle such “rate-form” constitutive expressions to solve, say, transient dynamic problems involving one or more bodies. We saw that this algorithm works equally nicely for linear infinitesimal elasticity as well.
In terms of both the non-rate form of the constitutive relationship and the rate-form, linear infinitesimal elasticity simplifies things considerably. Recall that the Second Piola Kirchhoff Stress Tensor, , can be used for infinitesimal elasticity (review the example at the end of Section 4: Alternative Measures of Stress). Since is defined in material coordinates, we will develop a constitutive relationship with a strain tensor that is similarly invariant to rigid body rotation, viz., .
Some authors instead choose to use the Cauchy stress and the Eulerian strain for the formulation of their linear infinitesimal constitutive expression. We are free to choose either stress-strain pair ( or ) in linear infinitesimal elasticity. We saw previously that our “time-stepping algorithm” was independent of this choice.
In this chapter we’ll take a look only at linear infinitesimal elasticity. If a material is subjected to large strain (even if it is linear) or if a material is nonlinear (even if subjected to small strain), then it should be treated using the approaches developed in the previous chapters on hyperelasticity and hypoelasticity. Just as in our chapter on hyperelasticity, the focus here will be on developing a non-rate constitutive expression that relates stress and strain. In the elastic regime, a wide variety of materials, including most metals, behave linearly.
The theory of linear infinitesimal elasticity was developed independent of hyperelastic theory, with the former being mostly developed, originally, from intuition and experience [Malvern]. Nevertheless, the linear theory can be derived as well, using the same starting assumptions that we used at in the previous hyperelasticity derivation.
Infinitesimal elasticity ,
note: The linear infinitesimal stress will be denoted by for the remainder of this chapter. So long as and are work conjugate, we can proceed in deriving our constitutive relationship.
We know the following three expressions to be true as well:
In index notation,
For linear materials, we can define , using a Taylor expansion, as follows:
note: This is a taylor series with higher terms neglected
note: In eq. 1,
If there is no initial stress:
Note that the order or doesn’t matter
has potential terms ( x matrix) to begin with.
Symmetry: and 36 terms (6 x 6 matrix)
Symmetry in (i.e. eq. 3 above)
21 constants (general anisotropic)
A plane of symmetry is defined as follows: for every pair of coordinate systems that are mirror images of each other about the plane of symmetry, the elastic constants are the same. A single plane of symmetry reduces the number of independent constants to 13. Orthotropic symmetry (three orthogonal planes of symmetry) reduces the number of independent constants to 9. The proof can be found in [Malvern].
Isotropic only 2 independent constants
Recall (valid outside of isotropy)
note: Since is a function of strain only, path dependence is not considered (i.e. we are still only considering elasticity)
Also note: is quadratic
recall that for hyperelastic isotropy,
Let’s assume , using the following reasoning:
We know that is quadratic. Thus, doesn’t depend on , since is cubic (). is linear and so is, accordingly, squared, in order to obtain a quadratic. Our resulting term is also multiplied by a coefficient, . is, similarly, multiplied by a constant, .
Recall that and is a function of both and .
Also recall that , since is a symmetric tensor.
* – note that “” must equal “” in the expression that is multiplied by “” in order to be non-zero.
In eq. 4, and are elastic constants. We can see that we have two independent constants, as expected. These constants are determined experimentally, as will be discussed shortly.
The following alternative expression is commonly found in literature on FEA:
We will see later what is, as well as why we use the “” operator.
note: is symmetric in matrix form
Also note: The shear terms are multiplied by “2” because we used symmetry of to reduce the 9 x 9 tensor, , to a 6 x 6.
The matrix “” can be written in “Voigt” notation, “” as well:
Using Voigt notation, we can construct the following table:
The above table can be simplified in this way due to the symmetries of the tensor apparent upon inspection.
e.g. the way you read the row depends on what information you need:
when I = 1, i = 1, j = 1
when J = 1, k = 1, l =1
This is easily confirmed by comparing, say, the “1-1” term in each of the above two matrices.
note: Eq. 6 is how stiffness matrices are formed in FEA
Sometimes we like to write the constitutive equation in the manner shown in eq. 5. We can prove that this expression is true by noting that the “:” operator is essentially two dot products (it can be easily shown that our previous definition of : for second order tensors still holds true as well).
, which is as expected.
note: symmetry of was invoked in the last step of the derivation