*HEADING
**
** Contact between aligned blocks in 3D.
**
*NODE,NSET=NB1
*ELEMENT,ELSET=EB1,TYPE=SC8RT
*ELGEN,ELSET=EB1
*NSET,NSET=FIX_N
*BOUNDARY
FIX_N,2,2
1,1,3
5,2,3
**
*NODE,NSET=NB2
*ELEMENT,ELSET=EB2,TYPE=SC8RT
*ELGEN,ELSET=EB2
**
*ELSET,ELSET=EALL
EB1,EB2
**
*SHELL SECTION,ELSET=EALL,MATERIAL=MAT1
1.0
*MATERIAL,NAME=MAT1
*ELASTIC
1.E6,0.2
*DENSITY
1.,
*CONDUCTIVITY
10.,
*SPECIFIC HEAT
0.001,
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16,STEFAN=1.E-6
*SURFACE,TYPE=ELEMENT,NAME=SB1
EB1,
*SURFACE,TYPE=ELEMENT,NAME=SB2
EB2,
**
*AMPLITUDE,NAME=SMOOTH,DEFINITION=SMOOTH STEP
0.,0.,1.E-3,-0.1,5.E-3,1.
**
*NSET, NSET=QA_TEST_ALLNODES
NB1, NB2
**
*INITIAL CONDITIONS,TYPE=TEMP
NB1, 100.
NB2, 0.
**
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT,EXPLICIT
**
*CONTACT PAIR,INTERACTION=LINEAR, MECHANICAL CONSTRAINT=PENALTY, WEIGHT=1.0
SB1,SB2
*SURFACE INTERACTION,NAME=LINEAR,PADTHICKNESS=0.1
*GAP CONDUCTANCE
**Gap conductance = k_interface / t_interface
*GAP HEAT GENERATION
*CONTACT DAMPING,DEFINITION=CRITICAL DAMPING FRACTION
*SURFACE BEHAVIOR,NO SEPARATION,PRESSURE-OVERCLOSURE=LINEAR
1.E6,
**Pressure-overclosure slope = E_interface / t_interface
*FRICTION,ROUGH,SHEAR TRACTION SLOPE=416666.7
**Shear traction slope = G_interface / t_interface
**                     = E_int / t_interface / (2*(1+nu_interface))
**
**
*OUTPUT, FIELD, NUMBER INTERVAL=1
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
U, NT
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=5e-4
*END STEP