1.3.34 Cylinder subjected to an asymmetric temperature field: CAXA elements

Product: ABAQUS/Standard  

Elements tested

CAXA4n    CAXA4Rn    CAXA8n    CAXA8Rn   

(n = 1, 2, 3, 4)


Problem description

A hollow cylinder of circular cross-section, inner radius , outer radius , and length , is subjected to an asymmetric temperature distribution that is a linear function of the spatial coordinates:

where is the constant temperature at the outside surface of the cylinder at 0° and r, , and z (see displacement solution, below) are the cylindrical coordinates. For a linear elastic material of Young's modulus E, Poisson's ratio , and thermal expansion coefficient , the solution for a structure subjected to such a temperature distribution is stress-free, with displacements as follows:

Only one-half of the structure is considered, with a symmetry plane at 0. The form of the displacement solution, which is a quadratic function in both r and z, indicates that a single second-order element can model the structure adequately and yield accurate results. This problem is also solved with an 8 × 12 mesh of fully integrated first-order elements and a 16 × 24 mesh of reduced integration first-order elements.

Material:

Linear elastic, Young's modulus = 30 × 106, Poisson's ratio = 0.33, coefficient of thermal expansion = 1 × 10–4.

Boundary conditions:

0 on the 0 plane; 0.06 is applied at and 0 to eliminate the rigid body motion in the global x-direction. This value of is obtained from the equation for above.

Loading:

A temperature field of the form is applied. This is accomplished by calculating the temperature at each node and defining the temperature value using the *TEMPERATURE option.

Results and discussion

The analytical solution and the ABAQUS results for the CAXA8n, CAXA8Rn, CAXA4n, and CAXA4Rn (n = 1, 2, 3 or 4) elements are tabulated below for a structure with these parameters: 6, 2, 6, and 300. The output locations are at points , , , and on the 0° plane, as shown in the figure on the previous page, and at points , and H, which are at the corresponding locations on the 180° plane. While both the CAXA8n and CAXA8Rn elements match the exact solution precisely with a zero state of stress, the models using the CAXA4n and CAXA4Rn elements fail to predict a stress-free state, even though the displacement solutions predicted are quite reasonable. However, the CAXA4Rn models give much more accurate results than the CAXA4n models. This example demonstrates that the fully integrated first-order elements do not handle bending problems very well.

VariableExactCAXA8nCAXA8RnCAXA4nCAXA4Rn
at A000–140710.0168
at A6 × 10–26 × 10–26 × 10–26 × 10–26 × 10–2
at A00000
at B00011664–3.2186
at B–3 × 10–2–3 × 10–2–3 × 10–2–2.9644 × 10–2–2.9999 × 10–2
at B6 × 10–26 × 10–26 × 10–26.0312 × 10–26.0001 × 10–2
at C000–140760.0162
at C1.4 × 10–21.4 × 10–21.4 × 10–21.3993 × 10–21.4 × 10–2
at C00000
at D00011108–3.5190
at D5 × 10–25 × 10–25 × 10–25.0306 × 10–25.0001 × 10–2
at D18 × 10–218 × 10–218 × 10–217.95 × 10–218 × 10–2
at E000–14071–0.0168
at E–6 × 10–2–6 × 10–2–6 × 10–2–6 × 10–2–6 × 10–2
at E00000
at F000–116643.2186
at F3 × 10–23 × 10–23 × 10–22.9644 × 10–22.9999 × 10–2
at F–6 × 10–2–6 × 10–2–6 × 10–2–6.0312 × 10–2–6.0001 × 10–2
at G000140763.5100
at G–1.4 × 10–2–1.4 × 10–2–1.4 × 10–2–1.3993 × 10–2–1.4 × 10–2
at G00000
at H000111083.5100
at H–5 × 10–2–5 × 10–2–5 × 10–2–5.0306 × 10–2–5.0001 × 10–2
at H–18 × 10–2–18 × 10–2–18 × 10–2–17.95 × 10–2–18 × 10–2

Note:  The results are independent of n, the number of Fourier modes.

Figure 1.3.34–1 through Figure 1.3.34–4 show plots of the undeformed and deformed meshes, the applied asymmetric temperature field, the contours of , and the contours of , respectively, for the CAXA84 model.

Input files

ecnssfsl.inp

CAXA41 elements.

ecntsfsl.inp

CAXA42 elements.

ecnusfsl.inp

CAXA43 elements.

ecnvsfsl.inp

CAXA44 elements.

ecnssrsl.inp

CAXA4R1 elements.

ecntsrsl.inp

CAXA4R2 elements.

ecnusrsl.inp

CAXA4R3 elements.

ecnvsrsl.inp

CAXA4R4 elements.

ecnwsfsl.inp

CAXA81 elements.

ecnxsfsl.inp

CAXA82 elements.

ecnysfsl.inp

CAXA83 elements.

ecnzsfsl.inp

CAXA84 elements.

ecnwsrsl.inp

CAXA8R1 elements.

ecnxsrsl.inp

CAXA8R2 elements.

ecnysrsl.inp

CAXA8R3 elements.

ecnzsrsl.inp

CAXA8R4 elements.

Figures

Figure 1.3.34–1 Deformed mesh.

Figure 1.3.34–2 Applied temperature field.

Figure 1.3.34–3 Contours of r-displacement.

Figure 1.3.34–4 Contours of z-displacement.