*HEADING
: NAFEMS TEST LE11: SOLID CYLINDER/TAPER/SPHERE - TEMPERATURE
COARSE MESH, C3D20 ELEMENTS
*NODE
   1,1.    , 0., 0.
  13,0.7071, 0., 0.7071
  21,0.7071, 0., 1.79
 401,1.4   , 0., 0.
 409,1.2124, 0., 0.7
 413,1.1062, 0., 1.045
 417,1.    , 0., 1.39
 421,1.    , 0., 1.79
*NGEN,LINE=C  
1,13,1,,0.,0.,0.
*NGEN
13,21
*NGEN,LINE=C
401,409,1,,0.,0.,0.
*NGEN
409,413
413,417
417,421
*NSET,NSET=IN,GEN
1,21,1
*NSET,NSET=OUT,GEN
401,421,1
*NFILL,NSET=XZPLANE
IN,OUT,4,100
*NCOPY,SHIFT,CHANGE NUMBER=1000,MULTIPLE=12,OLD SET=XZPLANE,
NEWSET=REST
0.,0.,0.
0.,0.,0.,0.,0.,1.,7.5
*NSET,NSET=YZPLANE,GEN
12001,12421,1
*NSET,NSET=XYPLANE,GEN
   1,12001,1000
 101,12101,1000
 201,12201,1000
 301,12301,1000
 401,12401,1000
*NSET,NSET=HI,GEN
  21,12021,1000
 121,12121,1000
 221,12221,1000
 321,12321,1000
 421,12421,1000
*NSET,NSET=ALLN
XZPLANE,REST
**
*ELEMENT,TYPE=C3D20
1, 1,5,405,401,4001,4005,4405,4401,3,205,403,201,4003,4205,4403,
4201,2001,2005,2405,2401
*ELGEN,ELSET=ALLE
1, 5,4,1, 3,4000,100
*ELSET,ELSET=EOUT
201, 
*SOLID SECTION,ELSET=ALLE,MATERIAL=MAT
*MATERIAL,NAME=MAT
*ELASTIC
210.E9,.3
*EXPANSION
2.3E-4,
*BOUNDARY
XZPLANE,2
YZPLANE,1
XYPLANE,3
HI,3
*STEP
*STATIC
*TEMPERATURE,USER
ALLN, 
*NODE PRINT,F=0
*EL PRINT,ELSET=EOUT,POSITION=AVERAGED AT NODES
S, 
*EL FILE,ELSET=EOUT,POSITION=AVERAGED AT NODES
S, 
*ENDSTEP