•
|
A position x in an element is identified, where the value is to be computed.
|
•
|
The total strain ε is computed from the derivatives of the shape functions at x. The strain depends on the values of the displacement degrees of freedom at the nodes and on the shape functions in the element.
|
•
|
The expression giving the thermal strain εth is evaluated at x. This expression will depend on a temperature, which may be prescribed or computed in another physics interface.
|
•
|
The creep strain εcr is a state variable, which is stored at the Gauss points in the element. The value is picked from the Gauss point closest to x.
|
•
|
The inelastic strain εinel (in this example, the creep strains is now subtracted from the total strains, to form the elastic strain εel.
|
•
|
The elastic strain εel is multiplied by the elasticity tensor , evaluated at x, to give the stress σ. The material properties may depend on the location, either explicitly, or for example by a temperature dependency.
|
•
|