*HEADING
EC3KTRSJ: THERMAL STRESS IN A CYLINDRICAL SHELL  [C3D20RT]
*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=C3D20RT,ELSET=ALLE
1, 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,
16,17,18,19,20
*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,10,11,13,14,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=SYMM,GENERATE
1,8,1
10,16,2
17,20,1            
*NSET,NSET=RIN
 1,2,3,4,9,10,11,12
*NSET,NSET=ROUT
 5,6,7,8,13,14,15,16
*TRANSFORM,NSET=ROTATE,TYPE=C
0.,0.,0.,0.,0.,1.
*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
*STEP
STRESS ANALYSIS USING TEMP DISTRIBUTION
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
1.,1.
*BOUNDARY
** THERMAL B.C.
 1, 11,11, 200.
 2, 11,11, 200.
 3, 11,11, 200.
 4, 11,11, 200.
 9, 11,11, 200.
10, 11,11, 200.
11, 11,11, 200.
12, 11,11, 200.
 5, 11,11, 100.
 6, 11,11, 100.
 7, 11,11, 100.
 8, 11,11, 100.
13, 11,11, 100.
14, 11,11, 100.
15, 11,11, 100.
16, 11,11, 100.
**DISPLACEMENT B.C.
ALLN, 2
*EL PRINT,ELSET=ALLE
S, 
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=ALLE
S, 
*NODE PRINT
U, 
NT, 
*ENDSTEP