*HEADING
UINTERCT: SIMILAR TO UFRICXXY, BUT USES
USER SUBROUTINE UINTER INSTEAD OF USER
SUBROUTINE FRIC
*RESTART,WRITE,F=5
*NODE
1,
2,2.,
3,2.,1.
4,,1.
5,,,2.
6,2.,,2.
7,2.,1.,2.
8,,1.,2.
9,0.5,1.,0.5
10,1.5,1.,0.5
11,1.5,2.,0.5
12,0.5,2.,0.5
13,0.5,1.,1.5
14,1.5,1.,1.5
15,1.5,2.,1.5
16,0.5,2.,1.5
*NSET,NSET=TOPE,GENERATE
9,16
*NSET,NSET=BOTE,GENERATE
1,8
*NSET,NSET=BOT
1,2,5,6
*NSET,NSET=TOP
11,12,15,16
*NSET,NSET=MAST
3,4,7,8
*ELEMENT,TYPE=C3D8T,ELSET=EALL
1,1,2,3,4,5,6,7,8
2,9,10,11,12,13,14,15,16
*SOLID SECTION,MATERIAL=SOLMAT,ELSET=EALL
*MATERIAL,NAME=SOLMAT
*ELASTIC
200.E3,0.
*EXPANSION
0.0, 
*CONDUCTIVITY
43.0, 
*SURFACE INTERACTION,NAME=FRI,USER,PROPERTIES=3

200.e7,0.001,0.00001
*surface,NAME=ASURF
1,S5
*surface,NAME=BSURF
2,S3
*CONTACT PAIR,INTERACTION=FRI,TYPE=SURFACE
BSURF,ASURF
*BOUNDARY
BOT,1,3
TOP,1
TOP,3
*AMPLITUDE,NAME=SPEC
0.,0.,0.1,0.1,0.2,0.2,0.3,0.3,
0.4,0.4,0.5,0.5,0.6,0.6,0.7,0.7,
0.8,0.8,0.9,0.9,1.0,1.0
*STEP,NLGEOM, UNSYMM=YES
BRING BLOCKS INTO CONTACT
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
*BOUNDARY
TOP,2,,-9.123E-2
TOPE,11,,0.0
BOTE,11,,0.0
*EL FILE
S, 
E, 
*NODE FILE
U, 
RF, 
NT, 
*PRINT,CONTACT=YES
*CONTACT FILE
CSTRESS,CDISP
HFL,HFLA,HTL
SFDR,SFDRA
*CONTACT PRINT
CSTRESS,CDISP
HFL,HFLA,HTL
SFDR,SFDRA
*OUTPUT,HISTORY,FREQ=0
*OUTPUT,FIELD,FREQ=0
*END STEP
*STEP,NLGEOM, UNSYMM=YES
MAINTAIN COMPRESSIVE STRESS
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
*BOUNDARY,OP=NEW
TOPE,11,,0.0
BOTE,11,,0.0
BOT,1,3
*CLOAD
TOP,2,-4000.
*PRINT,CONTACT=YES
*OUTPUT,HISTORY,FREQ=0
*OUTPUT,FIELD,FREQ=0
*END STEP
*STEP,NLGEOM,AMP=STEP, UNSYMM=YES
APPLY SLIDING MOTION
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
0.1,1.0,0.1,0.1
*BOUNDARY,FIXED
MAST,1,3
*BOUNDARY,AMP=SPEC
TOPE,11,,200.
BOTE,11,,0.0
*CLOAD
TOP,1,25.
*PRINT,CONTACT=YES
*EL PRINT,F=0
*NODE PRINT
NT, 
RF,
*OUTPUT,HISTORY,FREQ=0
*OUTPUT,FIELD,FREQ=0
*END STEP