*HEADING
: SIMPLE CONTACT TEST, 3-D SMALL SLIDING, DEFORMABLE, C3D10MH
** 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
426,.5,.5,10.
26,.5,.5,8.
*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
526,.5,.5,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
*NODE
926,.5,.5,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
*NSET,NSET=QA_TEST,GEN
1,400
450,1000
**
*ELEMENT,TYPE=C3D10M,ELSET=SOLID
1, 1,5,21,401,3,26,11,201,203,211
2, 401,421,425,21,411,423,426,211,221,223
3, 405,401,425,5,403,426,415,205,203,215
4, 21,5,25,425,26,15,23,223,215,225,
5, 401,425,5,21,426,215,203,211,223,26
6, 501,505,521,901,503,526,511,701,703,711
7, 901,921,925,521,911,923,926,711,721,723
8, 905,901,925,505,903,926,915,705,703,715
9, 521,505,525,925,526,515,523,723,715,725
10,901,925,505,521,926,715,703,711,723,526
**
*ELSET,ELSET=ETOP
5,6
*ELSET,ELSET=EBOT
1,2
*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
*CONTACT PRINT,SLAVE=TSURF,MASTER=BSURF
*ENERGY PRINT
*EL PRINT ,SUM=NO
S,
E,
*NODE PRINT,SUM=NO
RF,
*OUTPUT,FIELD,FREQ=0
*OUTPUT,HISTORY
*CONTACTOUTPUT, SLAVE=TSURF,MASTER=BSURF
CSTRESS
*END STEP
**
** step 2: slide top cube to the right w/o friction
**
*STEP,PERT
*STATIC
*BOUNDARY
TOP,1,,0.1
*ENERGY PRINT
*EL PRINT ,SUM=NO
S,
E,
*NODE PRINT,SUM=NO
RF,
*OUTPUT,FIELD,FREQ=0
*OUTPUT,HISTORY
*CONTACTOUTPUT, SLAVE=TSURF,MASTER=BSURF
CSTRESS
*END STEP
**
** step 3: change surface friction
**
*STEP,NLGEOM
*STATIC
1.0,1.0
*CHANGE FRICTION,INTERACTION=SLIDE
*FRICTION
0.3,
*OUTPUT,FIELD,FREQ=0
*OUTPUT,HISTORY
*CONTACTOUTPUT, SLAVE=TSURF,MASTER=BSURF
CSTRESS
*ENDSTEP
**
** step 4: slide to the left w/ friction
**
*STEP,PERT
*STATIC
*BOUNDARY
TOP,1,,-0.1
*OUTPUT,FIELD,FREQ=2
*ELEMENT OUTPUT,VAR=PRESELECT
*NODE OUTPUT,VAR=PRESELECT
*OUTPUT,HISTORY
*CONTACTOUTPUT, SLAVE=TSURF,MASTER=BSURF
CSTRESS
*ENDSTEP