*HEADING
HT PATCH TEST FOR COUPLED TEMPERATURE-DISPLACEMENT EULERIAN ELEMENTS  [EC3D8RT]
T=C0+C1x+C2y+C3z;C0=0.0,C1=200,C2=100,C3=200;HFL1=-kC1,HFL2=-kC2,HFL3=-kC3
Two materials with identical properties and adding up to
1.0 for volume fraction (fully filled eulerian element)
**
** See ch. 1.11.6 of ABAQUS/Verification Manual
**
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE, NSET=REFNODE
9999, 0.,0.,0.
*NODE, NSET=NALL
10,0.,0.,1.
20,1.,0.,1.
120,1.,0.,0.
110,0.,0.,0.
100,0.32,0.186,0.643
202,0.677,0.305,0.683
302,.826,.288,.288
400,.249,.342,.192
102,.165,.745,.702
200,.788,.693,.644
300,.850,.649,.263
402,.273,.750,.230
12,0.,1.,1.
22,1.,1.,1.
320,1.,1.,0.
420,0.,1.,0.
*NSET,NSET=SIDES
10, 
20, 
110, 
120, 
12, 
22, 
320, 
420, 
*ELEMENT,TYPE=EC3D8RT, ELSET=EALL
1,10,20,120,110,100,202,302,400
2,10,100,400,110,12,102,402,420
3,100,202,302,400,102,200,300,402
4,202,20,120,302,200,22,320,300
5,102,200,300,402,12,22,320,420
6,302,120,110,400,300,320,420,402
7,20,202,100,10,22,200,102,12
*ELSET,ELSET=EOUT
3
*ELSET,ELSET=QA_TEST
EOUT
*EULERIAN SECTION,ELSET=EALL
A1,emat1
A2,emat2
*MATERIAL,NAME=A1
*ELASTIC
1.0,0.0
*CONDUCTIVITY
4.85E-4, 
*DENSITY
0.283, 
*SPECIFIC HEAT
0.116, 
*EXPANSION
0.0, 
*MATERIAL,NAME=A2
*ELASTIC
1.0,0.0
*CONDUCTIVITY
4.85E-4, 
*DENSITY
0.283, 
*SPECIFIC HEAT
0.116, 
*EXPANSION
0.0, 
*INITIAL CONDITIONS, TYPE=TEMP
NALL,0.10
*initial conditions, type=volume fraction
eall, emat1, 0.75
eall, emat2, 0.25
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
,20.0
*BOUNDARY
 10,11, ,200.
 20,11, ,400.
120,11, ,200.
110,11, , 0.0
 12,11, ,300.
 22,11, ,500.
320,11, ,300.
420,11, ,100.
******************************
*OUTPUT,FIELD,NUM=20
*ELEMENT OUTPUT
TEMP,
*NODE OUTPUT, NSET=NALL
NT, RFL
***********
*OUTPUT,FIELD,NUM=1
*ELEMENT OUTPUT,ELSET=QA_TEST
TEMP,HFL
*OUTPUT, HISTORY,TIME INTERVAL=1.0
*ENERGY OUTPUT,VAR=ALL
******************************
*ENDSTEP