*HEADING
**
** Analytical contact with linear softening in 3D.
**
*NODE,NSET=REFNOD
1,2.,0.,0.5
*RIGID BODY,REF NODE=1,ANALYTICAL SURFACE=SB1,ISOTHERMAL=YES
*HEATCAP,ELSET=PTCAPAC
0.004,
*ELEMENT,TYPE=HEATCAP,ELSET=PTCAPAC
1,1
*BOUNDARY
1,1,6
**
*NODE,NSET=NB2
101,0.,1.1,0.
102,1.,1.1,0.
103,2.,1.1,0.
104,3.,1.1,0.
105,4.,1.1,0.
111,0.,2.1,0.
112,1.,2.1,0.
113,2.,2.1,0.
114,3.,2.1,0.
115,4.,2.1,0.
121,0.,1.1,1.
122,1.,1.1,1.
123,2.,1.1,1.
124,3.,1.1,1.
125,4.,1.1,1.
131,0.,2.1,1.
132,1.,2.1,1.
133,2.,2.1,1.
134,3.,2.1,1.
135,4.,2.1,1.
*ELEMENT,ELSET=EB2,TYPE=C3D8RT
101,101,102,112,111,121,122,132,131
*ELGEN,ELSET=EB2
101,4,1,1
*NSET,NSET=LOAD_N
111,112,113,114,115,131,132,133,134,135
**
*ELSET,ELSET=EALL
EB2,
**
*SOLID SECTION,ELSET=EALL,MATERIAL=MAT1
*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=CYLINDER,NAME=SB1
-1.,1.,2.,0.,1.,2.
-1.,1.,1.
START,0.,0.
LINE,6.,0.
*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
NB2, REFNOD
**
*INITIAL CONDITIONS,TYPE=TEMP
REFNOD, 100.
NB2, 0.
**
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT,EXPLICIT
,5.E-3
**
*BOUNDARY,AMPLITUDE=SMOOTH
LOAD_N,2,2,0.001
**
*CONTACT PAIR,INTERACTION=LINEAR, MECHANICAL CONSTRAINT=PENALTY
SB1,SB2
*SURFACE INTERACTION,NAME=LINEAR,PADTHICKNESS=0.1
*GAP CONDUCTANCE
1.,0.
1.,0.1
**Gap conductance = k_interface / t_interface
*GAP HEAT GENERATION
0.,
*CONTACT DAMPING,DEFINITION=CRITICAL DAMPING FRACTION
0.,
*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))
**
*FILE OUTPUT,NUM=1
*NODE FILE
U, NT
**
*OUTPUT, FIELD, NUMBER INTERVAL=1
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
U, NT
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.005
*END STEP