*HEADING
THERMAL STRESS IN A CYLINDRICAL SHELL [C3D4T]
**
** 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
*ELEMENT,TYPE=C3D4T,ELSET=ALLE
1, 1,9,11,13
2, 1,5,13,8
3, 11,15,8,13
4, 1,11,4,8
5, 1,11,8,13
6, 9,2,3,6
7, 9,13,6,15
8, 3,7,15,6
9, 9,3,11,15
10, 9,3,15,6
*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
*NSET,NSET=Z1
1,2,5,6,9,13,17,18
*NSET,NSET=Z2
3,4,7,8,11,15,19,20
*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
*STEP
STRESS ANALYSIS USING TEMP DISTRIBUTION
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 ,0.4
*FIXED MASS SCALING, ELSET=ALLE, FACTOR=1.E4
*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