Linear Elastic Material

For a linear elastic material, Hooke’s law relates the stress tensor to the elastic strain tensor:

is the 4th order elasticity tensor, “:” stands for the double-dot tensor product (or double contraction). The elastic strain εel is the difference between the total strain ε and all inelastic strains εinel. There may also be an extra stress contribution σex with contributions from initial stresses and viscoelastic stresses. In case of geometric nonlinearity, the second Piola-Kirchhoff stress tensor and the Green-Lagrange strain tensor are used.

The elastic strain energy density is

This expression assumes that the initial stress contribution is constant during the straining of the material.

Because of the symmetry, the strain tensor can be written as the following matrix:

Due to the symmetry, the elasticity tensor can be completely represented by a symmetric 6-by-6 matrix as:

which is the elasticity matrix.

In this case, the elasticity matrix becomes

Different pairs of elastic moduli can be used, and as long as two moduli are defined, the others can be computed according to Table 3-1.

Table 3-1: Expressions for the elastic moduli.
Description	Variable	D(E,ν)	D(E,G)	D(K,G)	D(λ,μ)
Young’s modulus	E =	E	E
Poisson’s ratio	ν =	ν
Bulk modulus	K =			K
Shear modulus	G =		G	G	μ
Lamé parameter λ	λ =				λ
Lamé parameter μ	μ =		G	G	μ
Pressure- wave speed	cp =
Shear-wave speed	cs =

According to Table 3-1, the elasticity matrix D for isotropic materials is written in terms of Lamé parameters λ and μ,

or in terms of the bulk modulus K and shear modulus G:

There are two different ways to represent orthotropic or anisotropic data. The converts the indices as:

thus, Hooke’s law is presented in the form involving the elasticity matrix D and the following vectors:

COMSOL Multiphysics uses the complete tensor representation internally to perform the coordinate system transformations correctly.

Beside the , the elasticity coefficients can be entered following the notation. In the , the sorting of indices is:

The last three rows and columns in the elasticity matrix D are thus swapped.

The elasticity matrix for orthotropic material in the has the following structure:

where the components are as follows:

where

The values of Ex, Ey, Ez, νxy, νyz, νxz, Gxy, Gyz, and Gxz are supplied in designated fields in the physics interface. COMSOL Multiphysics deduces the remaining components—νyx, νzx, and νzy—using the fact that the matrices D and D−1 are symmetric. The compliance matrix has the following form:

	The values of νxy and νyx are different for an orthotropic material. For a certain set of given material data, you must make sure that the definition of the indices is consistent with the definition used in COMSOL Multiphysics.

The elasticity matrix in the changes the sorting of the last three elements in the elasticity matrix:

In the general case of fully anisotropic material, you provide explicitly all 21 components of the symmetric elasticity matrix D, in either or .

For the linear elastic material, the stress components in coordinate system are

For anisotropic and orthotropic materials, the 4th-order elasticity tensor is defined from D matrix according to:

The user input D matrix always contains the physical components of the elasticity tensor

and the corresponding tensor components are computed internally according to:

For an isotropic material:

where λ and μ are the first and second Lamé elastic parameters and g is the metric tensor.

For a hyperelastic material, the second Piola-Kirchhoff stress tensor is computed as

which is computed as the contravariant components of the stress in the local coordinate system:

The energy variation is computed as

which can be also written as

The free energy for the linear thermoelastic material can be written as

where the strain energy density Ws(ε, T) is given by Equation 3-12. Hence, the stress can be found as

and the entropy per unit volume can be calculated as

where T0 is a reference temperature, the volumetric heat capacity ρCp can be assumed to be independent of the temperature (Dulong-Petit law), and the elastic entropy is

where α is the thermal expansion coefficient tensor. For an isotropic material, it simplifies into

The heat balance equation can be written as

where k is the thermal conductivity matrix, and the heat source caused by the dissipation is

where

is the strain-rate tensor and the tensor τ represents all possible inelastic stresses (for example, a viscous stress).

Using the tensor components, the heat balance can be rewritten as:

In many cases, the second term can be neglected in the left-hand side of Equation 3-13 because all Tαmn are small. The resulting approximation is often called uncoupled thermoelasticity.