*HEADING
 HEAT, GAP CONDUCTANCE, MF DEP USER SUBROUTINE - UTFVMHM
*NODE,NSET=PLANEA
*NCOPY,CHANGE NUMBER=100,OLD SET=PLANEA,SHIFT,NEW SET=PLANEB
*NCOPY,CHANGE NUMBER=200,OLD SET=PLANEA,SHIFT,NEW SET=PLANEC
*NCOPY,CHANGE NUMBER=300,OLD SET=PLANEA,SHIFT,NEW SET=PLANED
*NSET,NSET=BCPLANE
 PLANEB, PLANEC
*ELEMENT,TYPE=DCC3D8,ELSET=ONE
*ELEMENT,TYPE=DCC3D8,ELSET=TWO
*surface, NAME=SURF1
ONE, S1
*surface, NAME=SURF2
TWO, S2
*CONTACT PAIR, INTERACTION=GAP
SURF1, SURF2
*ELSET,ELSET=BOTH
ONE,TWO
*SOLID SECTION,ELSET=BOTH,MATERIAL=BLOCK
*MATERIAL,NAME=BLOCK
*CONDUCTIVITY
150.,
*DENSITY
1.,
*SPECIFIC HEAT
1.,
*SURFACE INTERACTION, NAME=GAP
*GAP CONDUCTANCE
*NSET,NSET=EDGE1,GEN
*NSET,NSET=EDGE3,GEN
*NSET,NSET=EDGE5,GEN
*NSET,NSET=EDGE7,GEN
*STEP,AMP=STEP
*HEAT TRANSFER,STEADY STATE
.25,1.,.25,.25
*BOUNDARY
PLANEA,11,11,1000.
PLANED,11,11
*MASS FLOW RATE,USER
EDGE1,
EDGE3,
EDGE5,
EDGE7,
*EL PRINT,F=0
*NODE PRINT
NT,
*NODE FILE
NT,
*OUTPUT,FIELD
*NODE OUTPUT
NT,
*CONTACT PRINT
*CONTACT FILE
*OUTPUT,FIELD
*CONTACT OUTPUT,VARIABLE=PRESELECT
*END STEP