2.1.1 Torsion of a hollow cylinder

Product: ABAQUS/Standard  

This problem verifies and illustrates the use of axisymmetric solid elements with twist in ABAQUS. An Airy stress function provided by Fung (1977) is used to obtain the stress components in the cylindrical coordinate system for this problem.

Problem description

The physical problem consists of a hollow cylinder fixed at one end in the translational degrees of freedom. Opposing twists, equal in magnitude, are applied to the inner and outer diameters of the hollow cylinder. Figure 2.1.1–1 shows the model geometry used in this analysis. Most of the meshes used for this problem are uniform; however, since the stresses are a quadratic function of the radius, mesh refinement in regions of high strain gradients is suggested. Two input files have refined meshes near the inner diameter of the hollow cylinder with appropriate mesh refinement *MPCs. Some input files use the *COUPLING and *KINEMATIC COUPLING options to couple the boundary surfaces to reference nodes; the twists are applied to the reference nodes. All axisymmetric solid elements with twist are tested in this problem; coupled temperature-displacement elements are tested in a *COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE step with arbitrarily applied boundary temperatures and a coefficient of thermal expansion, 0, to prevent coupling with the mechanical solution. Mesh convergence studies have not been performed. To test the use of substructures in problems involving axisymmetric elements with twist, the problem is also solved using substructuring with the entire model treated as a single substructure.

Results and discussion

The derived analytical solution for this problem, based on an Airy stress function provided by Fung (1977), is

where is the twist angle in radians, 1 = 2.05128 × 104, and 2 = –1.6667 × 10–2.

Figure 2.1.1–2 shows the variation of the shear stress, , with respect to the radius of the hollow cylinder for the CGAX8 element model (where S13) compared to the analytical solution. The results for all elements agree well with the analytical solution. The results from the substructure analysis match the results that are obtained when substructures are not used.

Input files

torsholcyl_cgax3_couplingk.inp

CGAX3 model using the *COUPLING and *KINEMATIC options to impose the twist.

torsholcyl_cgax3_kincoupl.inp

CGAX3 model using the *KINEMATIC COUPLING option to impose the twist.

torsholcyl_cgax4r_meshrefine.inp

CGAX4R model with mesh refinement.

torsholcyl_cgax4rh_eh.inp

CGAX4RH model with enhanced hourglass control.

torsholcyl_cgax8_substruct.inp

CGAX8 model using substructuring.

torsholcyl_cgax8_substruct_gen1.inp

Substructure generation referenced in the analysis torsholcyl_cgax8_substruct.inp.

torsholcyl_cgax8h_neohook.inp

CGAX8H model with a neo-Hookean incompressible hyperelastic material.

torsholcyl_cgax8r_meshrefine.inp

CGAX8R model with mesh refinement.

Reference

Figures

Figure 2.1.1–1 Torsion of a hollow cylinder.

Figure 2.1.1–2 Variation of shear stress with radius.