Torsion of a hollow cylinder

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.

This page discusses:

ProductsAbaqus/Standard

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 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 multi-point constraints. Some input files use a kinematic coupling 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 steady-state coupled temperature-displacement 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

ϕ=(1+ν)C1/Er2-C2,
σrθ=-C1/r2,

where ϕ is the twist angle in radians, C1 = 2.05128 × 104, and C2 = −1.6667 × 10−2.

Figure 2 shows the variation of the shear stress, σrθ, with respect to the radius of the hollow cylinder for the CGAX8 element model (where σrθ=-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_cgax3h.inp

CGAX3H model.

torsholcyl_cgax4.inp

CGAX4 model.

torsholcyl_cgax4h.inp

CGAX4H model.

torsholcyl_cgax4ht.inp

CGAX4HT model.

torsholcyl_cgax4r_meshrefine.inp

CGAX4R model with mesh refinement.

torsholcyl_cgax4rh.inp

CGAX4RH model.

torsholcyl_cgax4rh_eh.inp

CGAX4RH model with enhanced hourglass control.

torsholcyl_cgax4t.inp

CGAX4T model.

torsholcyl_cgax6.inp

CGAX6 model.

torsholcyl_cgax6h.inp

CGAX6H model.

torsholcyl_cgax6m.inp

CGAX6M model.

torsholcyl_cgax6mh.inp

CGAX6MH model.

torsholcyl_cgax8.inp

CGAX8 model.

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.inp

CGAX8H model.

torsholcyl_cgax8h_neohook.inp

CGAX8H model with a neo-Hookean incompressible hyperelastic material.

torsholcyl_cgax8ht.inp

CGAX8HT model.

torsholcyl_cgax8r_meshrefine.inp

CGAX8R model with mesh refinement.

torsholcyl_cgax8rh.inp

CGAX8RH model.

torsholcyl_cgax8rht.inp

CGAX8RHT model.

torsholcyl_cgax8rt.inp

CGAX8RT model.

torsholcyl_cgax8t.inp

CGAX8T model.

References

  1. Fung Y. C.Foundations of Solid Mechanics, Prentice-Hall Inc., New Jersey, 1977.

Figures

Figure 1. Torsion of a hollow cylinder.

Figure 2. Variation of shear stress with radius.