*HEADING
: SIMPLE CONTACT TEST, 3-D SMALL SLIDING, DEFORMABLE, C3D20
** Contact Features Tested:
**	*CONTACT INTERFERENCE w/SHRINK, *FRICTION,
**      *CHANGE FRICTION (mu=0 and mu=0.3),
*RESTART,WRITE, FREQ=4
*NODE, NSET=BOTTOM
1,0.,0.,8.,
5,1.,0.,8.,
25,1.,1.,8.
21,0.,1.,8.
*NODE,NSET=BSLIDE
401,0.,0.,10.
405,1.,0.,10.
421,0.,1.,10.
425,1.,1.,10.
*NODE,NSET=TSLID
501,-0.2,-0.2,9.9
505,1.2,-0.2,9.9
521,-0.2,1.2,9.9
525,1.2,1.2,9.9
*NODE,NSET=TOP
901,-0.2,-0.2,11.9
905,1.2,-0.2,11.9
921,-0.2,1.2,11.9
925,1.2,1.2,11.9
*NGEN, NSET=BOTTOM
1,5,2
21,25,2
1,21,10
5,25,10
*NGEN, NSET=BSLIDE
401,405,2
421,425,2
401,421,10
405,425,10
*NFILL
BOTTOM,BSLIDE, 2, 200
**
*NGEN, NSET=TOP
901,905,2
921,925,2
901,921,10
905,925,10
*NGEN ,NSET=TSLID
501,505,2
521,525,2
501,521,10
505,525,10
*NFILL
TSLID, TOP, 2, 200
**
*ELEMENT,TYPE=C3D20,ELSET=SOLID
1, 1,5,25,21,401,405,425,421,3,15,23,11,403,415,423,
411,201,205,225,221
3, 501,505,525,521,901,905,925,921,503,515,523,511,903,915,923,
911,701,705,725,721
**
*ELSET,ELSET=ETOP
3,
*ELSET,ELSET=EBOT
1,
*surface,NAME=TSURF
EBOT,S2
*surface,NAME=BSURF
ETOP,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
CSTRESS,CDISP
CFN,
CFS,
CAREA,
*CONTACT FILE,SLAVE=TSURF,MASTER=BSURF,F=4
CSTRESS,CDISP
CFN,
CFS,
CAREA,
*EL FILE ,FREQ=4
S,
E,
*NODE FILE,FREQ=4
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
U,
RF,
*CONTACT FILE,SLAVE=TSURF,MASTER=BSURF,F=6
**CSTRESS,CDISP
CFN,
CFS,
CAREA,
*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
*ENERGY FILE,FREQ=6
*EL FILE, FREQ=6
S,
E,
*NODE FILE, FREQ=6
U,
RF,
*CONTACT FILE,SLAVE=TSURF,MASTER=BSURF, FREQ=6
CSTRESS,CDISP
CFN,
CFS,
CAREA,
*CONTACT PRINT,SLAVE=TSURF,MASTER=BSURF, FREQ=6
CSTRESS,CDISP
CFN,
CFS,
CAREA,
*PRINT,CONTACT=YES
*ENDSTEP