VM33

VM33
Transient Thermal Stress in a Cylinder

Overview

Reference:R. J. Roark, W. C. Young, Formulas for Stress and Strain, 5th Edition, McGraw-Hill Book Co., Inc., New York, NY, 1975, pg. 585.
Analysis Type(s):Coupled-field Analysis (ANTYPE = 4)
Element Type(s):
3D Coupled-field Solid Element (SOLID5)
3D 20-Node Coupled-Field Solid (SOLID226)
Input Listing:vm33.dat

Test Case

A long thick-walled cylinder, initially at a uniform temperature To, has its outer radius temperature raised at a constant rate of 1.0°/sec to temperature Tf. After a steady state of heat flow has been reached, determine the tangential stress at the inner and outer surfaces. Display the outer-to-inner surface temperature difference and the tangential stress as a function of time.

Figure 43: Cylinder Problem Sketch

Cylinder Problem Sketch

Material PropertiesGeometric PropertiesLoading
E = 30 x 106 psi
υ = 0.3
α = 8.4 x 10-6 in/in-°F
k = .000625 BTU/sec in-°F
ρ = 0.284 lb-sec2/in4
c = .10 BTU/lb-°F
a = 1.0 in
b = 3.0 in
h = .20 in
Tb = 500°F
To = 70°F

Analysis Assumptions and Modeling Notes

Due to symmetry, only a wedge of arbitrary height is required for modeling. A 5° wedge is selected to minimize curved geometry effects when using a lower order element. The thermal steady state condition is satisfied when the inner and outer wall temperature difference is constant. A transient thermal-stress analysis is required with a sufficient time period to allow the steady-state condition to be obtained. A time period of t = 430 sec is selected. The temperature Tf is assigned a value of T = 500°F such that, for a ramped load condition, the constant temperature rise of: Tf-To/Δt = 500-70/430 = 1°F/sec is obtained. Since the structural dynamic effects are not of concern, inertial and damping structural effects can be ignored, by specifying time integration for the temperature degree of freedom only. A sufficient number of elements (15) is modeled through the thickness such that an accurate thermal transient and nodal stress results are obtained.

Symmetric structural boundary conditions are used at the radial and bottom planes. Since the cylinder being modeled is long, nodes at z = h are coupled in UZ to enforce a constant axial strain condition. The reported values at t = 430 sec. should be fairly accurate since thermal steady state is achieved.

Results Comparison

Tangential StressTargetMechanical APDLRatio
SOLID5
Stressypsi(r=b)-13396-130630.975
Stressypsi(r=a)10342103921.005
SOLID226
Stressypsi(r=b)-13396-128980.963
Stressypsi(r=a)1034299700.964

Figure 44: Outer-to-inner Surface Temperature Difference

Outer-to-inner Surface Temperature Difference

Figure 45: Tangential Stress as a Function of Time

Tangential Stress as a Function of Time