*HEADING
HT PATCH TEST FOR COUPLED TEMPERATURE-DISPLACEMENT ELEMENTS [C3D4T]
T=C0+C1x+C2y+C3z;C0=0.0,C1=200,C2=100,C3=200;HFL1=-kC1,HFL2=-kC2,HFL3=-kC3
**
** 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=OUT
 10,  0.,0.,1.
 30,  1.,0.,1.
 70,  0.,0.,0.
 90,  1.,0.,0.
 210, 0.,1.,1.
 230, 1.,1.,1.
 270, 0.,1.,0.
 290, 1.,1.,0.
*NODE,NSET=IN
 1010, .32,.286,.643
 1030, .677,.305,.683
 1070, .249,.342,.192
 1090, .826,.288,.288
 1210, .165,.745,.702
 1230, .788,.693,.644
 1270, .273,.750,.230
 1290, .850,.649,.263
*ELEMENT,TYPE=C3D4T,ELSET=E1
   1,1070,1090,1290,1030
   2,1030,1290,1210,1230
   3,1070,1290,1270,1210
   4,1070,1030,1210,1010
   5,1070,1030,1290,1210
   6,1070,1270,1290, 270
   7, 270,1290,  90, 290
   8,1070,1290,1090,  90
   9,1070, 270,  90,  70
  10,1070, 270,1290,  90
  11,1030,1230,1210, 230
  12, 230,1210,  10, 210
  13,1030,1210,1010,  10
  14,1030, 230,  10,  30
  15,1030, 230,1210,  10
  16,1070,1090,1030,  90
  17,  90,1030,  10,  30
  18,1070,1030,1010,  10
  19,1070,  90,  10,  70
  20,1070,  90,1030,  10
  21,1030,1090,1290,  90
  22,  90,1290, 230, 290
  23,1030,1290,1230, 230
  24,1030,  90, 230,  30
  25,1030,  90,1290, 230
  26,1210,1230,1290, 230
  27, 230,1290, 270, 290
  28,1210,1290,1270, 270
  29,1210, 230, 270, 210
  30,1210, 230,1290, 270
  31,1210,1270,1070, 270
  32, 270,1070,  10,  70
  33,1210,1070,1010,  10
  34,1210, 270,  10, 210
  35,1210, 270,1070,  10
*NSET,NSET=NALL
IN,OUT,
*ELSET,ELSET=EOUT
1,2,3,4,5,6,8,10,11,13,15,16,18,
20,21,23,25,26,28,30,31,33,35
*ELSET,ELSET=QA_TEST
EOUT,
*SOLID SECTION,MATERIAL=A1, ELSET=E1
*MATERIAL,NAME=A1
*ELASTIC
1.0,0.0
*CONDUCTIVITY
4.85E-4, 
*DENSITY
0.283, 
*SPECIFIC HEAT
0.116, 
*EXPANSION
0.0, 
*INITIAL CONDITION,TYPE=TEMP
OUT,0.0
IN,0.0
*RIGID BODY, ELSET = E1, REF NODE = REFNODE, ISOTHERMAL = NO
*STEP
 STEADY STATE
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
,20.0
*BOUNDARY
 10,11, , 200.
 30,11, , 400.
 70,11, ,  00.0
 90,11, , 200.
210,11, , 300.
230,11, , 500.
270,11, , 100.
290,11, , 300.
******************************
*FILE OUTPUT, NUM=1
*EL FILE,ELSET=EOUT
TEMP,HFL,
******************************
*OUTPUT,FIELD,NUM=1
*ELEMENT OUTPUT, ELSET=QA_TEST
TEMP,HFL,
*NODE OUTPUT, NSET=NALL
NT,
******************************
***OUTPUT, HISTORY,TIME INTERVAL=2
***NODE OUTPUT, NSET=NALL
**NT,
**************
*OUTPUT, HISTORY,TIME INTERVAL=20
*ENERGY OUTPUT,VAR=ALL
******************************
*ENDSTEP