10.2. Thermoelasticity

The capability to do a thermoelastic analysis exists in the following elements:

PLANE222 - 2D 4-Node Coupled-Field Solid
PLANE223 - 2D 8-Node Coupled-Field Solid
SOLID225 - 3D 8-Node Coupled-Field Solid
SOLID226 - 3D 20-Node Coupled-Field Solid
SOLID227 - 3D 10-Node Coupled-Field Solid
LINK228 - 3D Coupled-Field Link

These elements support both the thermal expansion and piezocaloric effects, and use the strong (matrix) coupling method.

In addition to the above, the following elements support the thermal expansion effect only in the form of a thermal strain load vector, i.e. use weak coupling method:

SOLID5 - 3D 8-Node Coupled-Field Solid
PLANE13 - 2D 4-Node Coupled-Field Solid
SOLID98 - 3D 10-Node Coupled-Field Solid

Constitutive Equations of Thermoelasticity

The coupled thermoelastic constitutive equations (Nye([356])) are:

(10–37)

(10–38)

where:

{ε} = total strain vector = [εx εy εz εxy εyz εxz]T
S = entropy density
{σ} = stress vector = [σx σy σz σxy σyz σxz]T
ΔT = temperature change = T - Tref
T = current temperature
T0 = absolute reference temperature = Tref + Toff
Tref = reference temperature (input on TREF command or as REFT on MP command)
Toff = offset temperature from absolute zero to zero (input on TOFFST command)
[D] = elastic stiffness matrix (inverse defined in Equation 2–4 or input using TB,ANEL command)
{α} = vector of coefficients of thermal expansion = [α x α y α z 0 0 0]T (input using, for example, ALPX, ALPY, ALPZ on MP command)
ρ = density (input as DENS on MP command)
Cp = specific heat at constant stress or pressure (input as C on MP command)

Using {ε} and ΔT as independent variables, and replacing the entropy density S in Equation 10–38 by heat density Q using the second law of thermodynamics for a reversible change

(10–39)

we obtain

(10–40)

(10–41)

where:

{β} = vector of thermoelastic coefficients = [D] {α}
= specific heat at constant strain or volume (see Note)

Substituting Q from Equation 10–41 into the heat flow equation Equation 6–1 produces:

(10–42)

where:

= thermal conductivity matrix
, , = thermal conductivities (input as KXX, KYY, KZZ on MP command)

Derivation of Thermoelastic Matrices

Applying the variational principle to stress equation of motion and the heat flow conservation equation coupled by the thermoelastic constitutive equations, produces the following finite element matrix equation:

(10–43)

where:

[M] = element mass matrix (defined by Equation 2–58)
[C] = element structural damping matrix (discussed in Damping Matrices)
[K] = element stiffness matrix (defined by Equation 2–58)
{u} = displacement vector
{F} = sum of the element nodal force (defined by Equation 2–56) and element pressure (defined by Equation 2–58) vectors
[Ct] = element specific heat matrix (defined by Equation 6–28)
[Kt] = element thermal conductivity matrix (defined by Equation 6–28)
{T} = temperature vector
{Q} = sum of the element heat generation load and element convection surface heat flow vectors (defined by Equation 6–28)
= element thermoelastic stiffness matrix
[B] = strain-displacement matrix (see Equation 2–44)
{N} = element shape functions
= element thermoelastic damping matrix (see Note)

Energy Calculation

In static and transient thermoelastic analyses, the element instantaneous total strain energy is calculated as:

(10–44)

where:

Ut = total strain energy (output as an NMISC element item UT).

Note that Equation 10–44 uses the total strain, whereas the standard strain energy (output as SENE) uses the elastic strain.

In a harmonic thermoelastic analysis, the time-averaged element total strain energy is given by:

(10–45)

where:

{ε}* = complex conjugate of the total strain

The real part of Equation 10–45 represents the average stored strain energy, while its imaginary part - the average energy loss due to thermoelastic damping.

The thermoelastic damping can be quantified by the quality factor Q derived from the total strain energy Equation 10–45 using the real and imaginary solution sets:

(10–46)

where:

Ne = number of thermoelastic elements


Note:  When KEYOPT(9) is set to 2 in a transient analysis, the program uses the actual absolute temperature instead of the reference temperature for improved accuracy.