*HEADING
: SIMPLE CONTACT TEST, 3-D SMALL SLIDING, DEFORMABLE, C3D10M
** Contact Features Tested:
**	*CONTACT INTERFERENCE w/SHRINK, *FRICTION,
**      *CHANGE FRICTION (mu=0 and mu=0.3),
*RESTART,WRITE, FREQ=4
*NODE, NSET=BOTTOM
*NODE,NSET=BSLIDE
*NODE
*NODE,NSET=TSLID
*NODE
*NODE,NSET=TOP
*NODE
*NGEN, NSET=BOTTOM
*NGEN, NSET=BSLIDE
*NFILL
BOTTOM,BSLIDE, 2, 200
**
*NGEN, NSET=TOP
*NGEN ,NSET=TSLID
*NFILL
TSLID, TOP, 2, 200
*NSET,NSET=QA_TEST,GEN
**
*ELEMENT,TYPE=C3D10M,ELSET=SOLID
**
*ELSET,ELSET=ETOP
*ELSET,ELSET=EBOT
*surface,NAME=TSURF
2,S1
3,S1
*surface,NAME=BSURF
9,S1
6,S1
*CONTACT PAIR,INTERACTION=SLIDE,SMALL SLIDING, TYPE=SURFACE TO SURFACE
TSURF,BSURF
*SURFACE INTERACTION,NAME=SLIDE
*FRICTION
0.,
*SOLID SECTION,ELSET=SOLID,MATERIAL=SOLMAT
*MATERIAL,NAME=SOLMAT
*ELASTIC
30.E6,.0
*BOUNDARY
TOP,1,3
BOTTOM,1,3
**
** step 1: Shrink fit with *CONTACT INTERFERENCE.
** 	   Establish base state.
**
*STEP ,NLGEOM
*STATIC
0.25,1.
*CONTACT INTERFERENCE, SHRINK
TSURF,BSURF
*ENERGY FILE,FREQ=4
*CONTACT PRINT,SLAVE=TSURF,MASTER=BSURF
*CONTACT FILE,NSET=QA_TEST,F=4
*EL FILE ,FREQ=4
S,
E,
*NODE FILE,FREQ=4,NSET=QA_TEST
RF,
*ENERGY PRINT
*EL PRINT ,SUM=NO
S,
E,
*NODE PRINT,SUM=NO
RF,
*END STEP
**
** step 2: slide top cube to the right w/o friction
**
*STEP,PERT
*STATIC
*BOUNDARY
TOP,1,,0.1
*ENERGY FILE,FREQ=6
*EL FILE ,FREQ=6
S,
E,
*NODE FILE,FREQ=6,NSET=QA_TEST
U,
RF,
*CONTACT FILE,NSET=QA_TEST,F=6
*ENERGY PRINT
*EL PRINT ,SUM=NO
S,
E,
*NODE PRINT,SUM=NO
RF,
*END STEP
**
** step 3: change surface friction
**
*STEP,NLGEOM
*STATIC
1.0,1.0
*CHANGE FRICTION,INTERACTION=SLIDE
*FRICTION
0.3,
*ENDSTEP
**
** step 4: slide to the left w/ friction
**
*STEP,PERT
*STATIC
*BOUNDARY
TOP,1,,-0.1
*EL FILE, FREQ=6
S,
E,
*NODE FILE, FREQ=6
U,
RF,
*CONTACT FILE ,NSET=QA_TEST, FREQ=6
*CONTACT PRINT,NSET=QA_TEST, FREQ=6
*PRINT,CONTACT=YES
*ENDSTEP