*HEADING
Stress analysis: 9 x 9 mesh with c3d8h. read and interpolate temperature from pgc38ths.inp
*NODE,NSET=ALL
 1, 3.
91, 8.
100, 8.,5.
10, 3.,5.
101, 3.,, .5
191, 8.,, .5
200, 8.,5., .5
110, 3.,5., .5
201, 3.,, 1.
291, 8.,, 1.
300, 8.,5., 1.
210, 3.,5., 1.
*NGEN,NSET=NL1
1,10,1
*NGEN,NSET=NR1
91,100
*NFILL,NSET=ALL1
NL1,NR1,9,10
*NGEN,NSET=NL2
101,110,1
*NGEN,NSET=NR2
191,200
*NFILL,NSET=ALL2
NL2,NR2,9,10
*NGEN,NSET=NL3
201,210,1
*NGEN,NSET=NR3
291,300,1
*NFILL,NSET=ALL3
NL3,NR3,9,10
*NSET,NSET=FIX1,GEN
 11,91,10
** 111,191,10
 211,291,10
*NSET,NSET=FIX2,GEN
91,100
** 191,200
291,300
*NSET,NSET=FIX3,GEN
10,100,10
** 110,200,10
210,300,10
*NSET,NSET=FIX4,GEN
1,10
** 101,110
201,210
*NSET,NSET=ALL
ALL1,ALL2,ALL3
*NSET,NSET=NALL
FIX1,FIX2,FIX3,FIX4
*ELEMENT, TYPE=C3D8H, ELSET=EALL
1,1,11,12,2, 201,211,212,202
*ELGEN,ELSET=EALL
1,9,1,1,9,10,10
*ELSET,ELSET=EFLUX2,GEN
81,89
*SOLID SECTION,MATERIAL=A1, ELSET=EALL
*MATERIAL,NAME=A1
*ELASTIC
 30.E6, .3
*EXPANSION
 1.E-5, 
*CONDUCTIVITY
 3.77E-5, 
*SPECIFIC HEAT
 .39, 
*DENSITY
 82.9, 
*RESTART,WRITE
*STEP
*STATIC
 1.0,1.0
*BOUNDARY
 ALL1,  3
 ALL3,  3
 FIX1, 1, 2
*DLOAD,OP=NEW
 EFLUX2, P2, -10000.
*TEMPERATURE, file=pgc38ths.odb, interpolate
*ENERGY FILE
*EL FILE
 LOADS
*NODE FILE
 U, 
 RF, 
 NT, 
 RFL, 
*EL PRINT
 LOADS
*NODE PRINT
 U, 
 RF, 
 NT, 
 RFL, 
*END STEP