*HEADING
THERMAL STRESS IN A CYLINDRICAL SHELL [CGAX4RHT]
**
** See ch. 1.3.16 of ABAQUS/Verification Manual
**
*NODE,NSET=ALLN
1, 0.0995,0.0
2, 0.1005,0.0
3, 0.1005,0.02
4, 0.0995,0.02
5, 0.099833,0.0
6, 0.1005,0.01
7, 0.099833,0.02
8, 0.0995,0.01
9, 0.1001667,0.0
10,0.1001667,0.02
**
*NSET,NSET=INSIDE
1,4
*NSET,NSET=OUTSIDE
2,3
*NSET,NSET=MID1
5,9
*NSET,NSET=MID2
7,10
*ELEMENT,TYPE=CGAX4RHT,ELSET=ALLE
1,1,5,7,4
2,5,9,10,7
3,9,2,3,10
*ELSET,ELSET=QA_TEST
ALLE,
*SOLID SECTION,ELSET=ALLE,MATERIAL=STEEL
*MATERIAL,NAME=STEEL
*DENSITY
7800., 
*CONDUCTIVITY
52., 
*SPECIFIC HEAT
586., 
*ELASTIC
200.E9, 0.3                          
*EXPANSION
1.2E-5, 
*EQUATION
**RADIAL
2, 
INSIDE,1,1.0,OUTSIDE,1,-1.0
**AXIAL
2, 
INSIDE,2,1.0,OUTSIDE,2,-1.0
**RADIAL MID NODES
2    ,      
MID1,1,1.0,2,1,-1.0
2    ,      
MID2,1,1.0,3,1,-1.0         
**AXIAL MID NODES
2    ,     
5,2,1.0,7,2,1.0
2    ,     
9,2,1.0,10,2,1.0
*STEP
STRESS ANALYSIS USING TEMP DISTRIBUTION
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
*BOUNDARY
ALLN,5,,
OUTSIDE, 11,, 100.
INSIDE,  11,, 200.
*EL PRINT,ELSET=ALLE
S,
*NODE PRINT, FREQ=0
*EL PRINT
 HFL,
*EL PRINT,POSITION=NODES
 HFL,
*EL PRINT,POSITION=AVERAGED AT NODES
 HFL,
*EL PRINT,POSITION=CENTROID
 HFL,
*ENDSTEP