*HEADING
: SIMPLE CONTACT TEST, 3-D SMALL SLIDING, DEFORMABLE
** C3D20, C3D15, S8R5
** Contact Features Tested:
**	*CONTACT INTERFERENCE w/SHRINK, *FRICTION,
**      *CHANGE FRICTION (mu=0 and mu=0.3),
**	*surface,type=node
**    NOTE: BOTH FRICTION STEPS ARE PERTURBATION STEPS
*RESTART,WRITE, FREQ=6
*NODE
*NGEN,NSET=BOTBR
*NCOPY,OLDSET=BOTBR,NEWSET=N2,SHIFT,CHANGE=10
*NCOPY,OLDSET=N2,NEWSET=N3,SHIFT,CHANGE=10
*NSET,NSET=NBRICK
BOTBR,N2,N3
*NCOPY,OLDSET=NBRICK,NEWSET=NWEDGE,SHIFT,CHANGE=100
*NCOPY,OLDSET=NWEDGE,NEWSET=NSHELL,SHIFT,CHANGE=100
**
*NGEN, NSET=T1
*NCOPY,OLDSET=T1,NEWSET=T2,SHIFT,CHANGE=10
*NCOPY,OLDSET=T2,NEWSET=T3,SHIFT,CHANGE=10
*NSET, NSET=TBRK
T1,T2,T3
**
*NCOPY,OLDSET=TBRK,NEWSET=TWEDGE,SHIFT,CHANGE=100
*NCOPY,OLDSET=TWEDGE,NEWSET=TSHELL,SHIFT,CHANGE=100
**
**  BRICK ELEMENTS:
*ELEMENT,TYPE=C3D20
*ELGEN,ELSET=EBRICK
**
**  WEDGE ELEMENTS:
*ELEMENT,TYPE=C3D15
*ELGEN,ELSET=EWEDGE
**
**  SHELL ELEMENTS:
*ELEMENT,TYPE=S8R5,ELSET=ESHELL
**
*SOLID SECTION,ELSET=EBRICK,MATERIAL=MAT
*SOLID SECTION,ELSET=EWEDGE,MATERIAL=MAT
*SHELL SECTION,ELSET=ESHELL,MATERIAL=MAT,TEMP=3
0.1,
**
*MATERIAL,NAME=MAT
*ELASTIC
3.0E6,0.3
*EXPANSION
1.0E-6,
**
*surface,NAME=SURF1A
1,S2
*surface,NAME=SURF1B
11,S1
*surface,NAME=SURF2A
3,S4
*surface,NAME=SURF2B
12,S3
*surface,NAME=SURF3B
14,SNEG
*surface,NAME=SURF3A
4,SPOS
*CONTACT PAIR,INTERACTION=CNTCT,SMALL, TYPE=SURFACE TO SURFACE
SURF1A,SURF1B
SURF2A,SURF2B
SURF3A,SURF3B
*SURFACE INTERACTION,NAME=CNTCT
*FRICTION
0.,
**
*NSET,NSET=BOTTOM,GEN
*NSET,NSET=TOP,GEN
*NSET,NSET=TOPSHL,GEN
*NSET, NSET=BOTSHL, GEN
**
*ELEMENT,TYPE=SPRING1,ELSET=SPR
*SPRING,ELSET=SPR
3,
1.0E-2,
**
*BOUNDARY
TOP, 1,3
BOTTOM, 1,3
**221,1,2
221,4,6
**227,1,2
TOPSHL,1,2
BOTSHL,1,2
**
** step 1: Shrink fit with *CONTACT INTERFERENCE.
** 	   Establish base state.
**
*STEP ,NLGEOM
*STATIC
0.25,1.
*CONTROLS, PARAMETER=FIELD
,2.E-02
*CONTACT INTERFERENCE, SHRINK
SURF1A, SURF1B
*CONTACT INTERFERENCE, SHRINK
SURF2A, SURF2B
*CONTACT INTERFERENCE, SHRINK
SURF3A, SURF3B
*ENERGY FILE,FREQ=4
*CONTACT PRINT,SLAVE=SURF1A,MASTER=SURF1B
*CONTACT PRINT,SLAVE=SURF2A,MASTER=SURF2B
*CONTACT PRINT,SLAVE=SURF3A,MASTER=SURF3B
*CONTACT FILE,SLAVE=SURF3A,MASTER=SURF3B,F=4
*CONTACT FILE,SLAVE=SURF2A,MASTER=SURF2B,F=4
*CONTACT FILE,SLAVE=SURF1A,MASTER=SURF1B,F=4
*EL FILE ,FREQ=4
S,
E,
*NODE FILE,FREQ=4
RF,
*ENERGY PRINT, FREQ=4
*EL PRINT ,SUM=NO, FREQ=4
S,
E,
*NODE PRINT,SUM=NO, FREQ=4
RF,
*END STEP
**
** step 2: slide top cube to the right w/o friction
**
*STEP,PERT
*STATIC
*CONTROLS, PARAMETER=FIELD, RESET
*BOUNDARY
TOP,1,,0.1
TOPSHL,1,, 0.1
*ENERGY FILE,FREQ=6
*EL FILE ,FREQ=6
S,
E,
*NODE FILE,FREQ=6
RF,
*CONTACT FILE,SLAVE=SURF3A,MASTER=SURF3B,F=6
*CONTACT FILE,SLAVE=SURF2A,MASTER=SURF2B,F=6
*CONTACT FILE,SLAVE=SURF1A,MASTER=SURF1B,F=6
*ENERGY PRINT, FREQ=6
*EL PRINT ,SUM=NO, FREQ=6
S,
E,
*NODE PRINT,SUM=NO, FREQ=6
RF,
*END STEP
**
** step 3: change surface friction
**
*STEP,NLGEOM
*STATIC
1.0,1.0
*CHANGE FRICTION,INTERACTION=CNTCT
*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,
*CONTACT PRINT,SLAVE=SURF1A,MASTER=SURF1B, FREQ=6
*CONTACT PRINT,SLAVE=SURF2A,MASTER=SURF2B, FREQ=6
*CONTACT PRINT,SLAVE=SURF3A,MASTER=SURF3B, FREQ=6
*CONTACT FILE,SLAVE=SURF3A,MASTER=SURF3B,F=6
*CONTACT FILE,SLAVE=SURF2A,MASTER=SURF2B,F=6
*CONTACT FILE,SLAVE=SURF1A,MASTER=SURF1B,F=6
*ENDSTEP