*HEADING
THERMAL STRESS IN A CYLINDRICAL SHELL [C3D8T]
**
** See ch. 1.3.16 of ABAQUS/Verification Manual
**
*NODE,NSET=ALLN,SYSTEM=C
 1, 0.0995,0.0,0.0
 2, 0.0995,11.25,0.0
 3, 0.0995,11.25,0.02
 4, 0.0995,0.0,0.02
 5, 0.1005,0.0,0.0
 6, 0.1005,11.25,0.0
 7, 0.1005,11.25,0.02
 8, 0.1005,0.0,0.02
 9, 0.0995,5.625,0.0
10, 0.0995,11.25,0.01
11, 0.0995,5.625,0.02
12, 0.0995,0.0,0.01
13, 0.1005,5.625,0.0
14, 0.1005,11.25,0.01                                   
15, 0.1005,5.625,0.02
16, 0.1005,0.0,0.01
17, 0.1,0.0,0.0
18, 0.1,11.25,0.0
19, 0.1,11.25,0.02
20, 0.1,0.0,0.02
21, 0.1,5.625,0.0
22, 0.1,5.625,0.02
*ELEMENT,TYPE=C3D8T,ELSET=ALLE
1,1,9,11,4,17,21,22,20
2,17,21,22,20,5,13,15,8
3,9,2,3,11,21,18,19,22
4,21,18,19,22,13,6,7,15
*ELSET,ELSET=QA_TEST
ALLE
*SOLID SECTION,ELSET=ALLE,MATERIAL=STEEL,ORIE=OR1
*ORIENTATION,SYS=CYL,NAME=OR1
 0.0,0.0,0.0, 0.0, 0.0, 1.0
 1, 0.0
*MATERIAL,NAME=STEEL
*DENSITY
7800., 
*CONDUCTIVITY
52., 
*SPECIFIC HEAT
586., 
*ELASTIC
200.E9, 0.3                          
*EXPANSION
1.2E-5, 
*NSET,NSET=ROTATE
2,3,6,7,9,11,13,15,18,19,21,22
*NSET,NSET=Z1
1,2,5,6,9,13,17,18,21
*NSET,NSET=Z2
3,4,7,8,11,15,19,20,22
*NSET,NSET=RIN
 1,2,3,4,9,11
*NSET,NSET=ROUT
 5,6,7,8,13,15
*TRANSFORM,NSET=ROTATE,TYPE=C
0.,0.,0.,0.,0.,1.
*AMPLITUDE,NAME=RAMP
0.0,0.0,0.4,1.0
*EQUATION
** RADIAL
2, 
RIN,1,1.0,ROUT,1,-1.0
**AXIAL
2    ,      ,         
RIN,3,1.0,ROUT,3,-1.0
**RADIAL MID NODES
2    ,      
17,1,1.0,5,1,-1.0
2    ,      
18,1,1.0,6,1,-1.0
2    ,      
19,1,1.0,7,1,-1.0
2    ,      
20,1,1.0,8,1,-1.0
**AXIAL MID NODES
2    ,      
17,3,1.0,20,3,1.0
2    ,      
18,3,1.0,19,3,1.0
2    ,      
19,3,1.0,20,3,-1.0
2, 
21,3,1.0,22,3,1.0
2, 
22,3,1.0,20,3,-1.0
*STEP
STRESS ANALYSIS USING TEMP DISTRIBUTION
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 ,0.4
*FIXED MASS SCALING, ELSET=ALLE, FACTOR=1.E4
*MONITOR, NODE=1, DOF=11
*BOUNDARY,AMP=RAMP
** THERMAL B.C.
 RIN, 11,11,200.0
 ROUT,11,11,100.0
*BOUNDARY,TYPE=VELOCITY
**DISPLACEMENT B.C.
ALLN,2
***************************
*FILE OUTPUT, NUM=1
*EL FILE, ELSET=ALLE
Temp,S
***************************
*OUTPUT,FIELD,NUMBER INTERVAL=1
*ELEMENT OUTPUT, ELSET=QA_TEST
Temp,S
**************************
***OUTPUT,HISTORY,TIME INTERVAL=4.E-2
***ELEMENT OUTPUT, ELSET=ALLE
**Temp,S,
***********
*OUTPUT,HISTORY,TIME INTERVAL=0.4
*ENERGY OUTPUT,VAR=ALL
************************
*ENDSTEP