VM38

VM38
Internal Pressure Loading of a Thick-Walled Cylinder

Overview

Reference: S. Timoshenko, Strength of Material, Part II, Elementary Theory and Problems, 3rd Edition, D. Van Nostrand Co., Inc., New York, NY, 1956, pg. 388, article 70.
Analysis Type(s): Static Analysis (ANTYPE = 0)
Element Type(s):
2D Structural Solid Elements (PLANE182)
2D Structural Surface Effect Elements (SURF153)
3D Structural Solid Elements (SOLID185)
3D Structural Surface Effect Elements (SURF154)
Input Listing: vm38.dat

Test Case

A long thick-walled cylinder is subjected to an internal pressure p (with no end cap load). Determine the radial stress, σr, and the tangential (hoop) stress, σt, at locations near the inner and outer surfaces of the cylinder for a pressure, pel, just below the yield strength of the material, a fully elastic material condition.

Figure 51: Thick-Walled Cylinder Problem Sketch

Thick-Walled Cylinder Problem Sketch

Material PropertiesGeometric PropertiesLoading
E = 30 x 106 psi
σyp = 30,000 psi
υ = 0.3
a = 4 in
b = 8 in
pel = 12,990 psi

Analysis Assumptions and Modeling Notes

The theory available for this problem is based on the Tresca (maximum shear) yield criterion while Mechanical APDL uses the von Mises yield criterion. The applied pult pressure is calculated from the Tresca theory by using . This procedure is sufficient to calculate approximate loads but the resulting nonlinear stress components should not be compared directly.

The problem is solved first using axisymmetric solid elements (PLANE182) and then using 3D solid elements (SOLID185). Since the problem is axisymmetric, only a small sector (5°) is modeled with SOLID185. In order to ensure constant axial strain (implied by the "long" cylinder definition), nodal coupling is used with PLANE182 and SOLID185. To illustrate the use of surface effect elements, the internal pressure P is applied using 2D structural surface effect elements (SURF153) in the first analysis, whereas 3D structural surface effect elements (SURF154) are used in the second analysis. Results are obtained from the solution phase and from the element centroid data.

Results Comparison

 TargetMechanical APDLRatio
PLANE182: Fully ElasticStressr, psi (X=4.4 in)-9,984.-9,986.1.000
Stresst, psi (X=4.4 in)18,645.18,809.1.009
Stressr, psi (X=7.6 in)-468.-467.0.999
Stresst, psi (X=7.6 in)9,128.9,120.0.999
SOLID185: Fully ElasticStressr, psi (X=4.4 in)-9,984.-9,986.1.000
Stresst, psi (X=4.4 in)18,645.18,809.1.009
Stressr, psi (X=7.6 in)-468.-467.0.999
Stresst, psi (X=7.6 in)9,128.9,120.0.999

  1. Output quantity SEQV