*HEADING
 VERIFICATION TEST FOR FINITE SLIDING CONTACT PATCH ALGORITHM
 USING C3D10M ELEMENTS
*RESTART, WRITE, FREQ=1
*NODE, NSET=NB1
*NGEN, NSET=NB1
*NCOPY, CHANGE NUMBER=500,  OLDSET=NB1, SHIFT, NEW SET=NB2
*NCOPY, CHANGE NUMBER=1000, OLDSET=NB1, SHIFT, NEW SET=NB3
*NCOPY, CHANGE NUMBER=1500, OLDSET=NB1, SHIFT, NEW SET=NB4
*NCOPY, CHANGE NUMBER=2000, OLDSET=NB1, SHIFT, NEW SET=NB5
*NCOPY, CHANGE NUMBER=2500, OLDSET=NB1, SHIFT, NEW SET=NB5
*NCOPY, CHANGE NUMBER=3000, OLDSET=NB1, SHIFT, NEW SET=NB6
*NCOPY, CHANGE NUMBER=3500, OLDSET=NB1, SHIFT, NEW SET=NB7
*NCOPY, CHANGE NUMBER=4000, OLDSET=NB1, SHIFT, NEW SET=NB8
*NCOPY, CHANGE NUMBER=4500, OLDSET=NB1, SHIFT, NEW SET=NB9
*NCOPY, CHANGE NUMBER=5000, OLDSET=NB1, SHIFT, NEW SET=NB10
*NCOPY, CHANGE NUMBER=5500, OLDSET=NB1, SHIFT, NEW SET=NB11
*NCOPY, CHANGE NUMBER=6000, OLDSET=NB1, SHIFT, NEW SET=NB12
*NSET, NSET=NBASE
 NB1, NB2, NB3, NB4, NB5, NB6, NB8, NB9, NB10, NB11, NB12
*SYSTEM
*NODE, NSET=NBLOCK1
*SYSTEM
*NODE, NSET=NBLOCK2
*NSET, NSET=NMOVE1
*NSET, NSET=NMOVE2
*NSET, NSET=NBASE, GEN
*NSET, NSET=NALL
 NBASE,
 NBLOCK1,
 NBLOCK2,
*ELEMENT, TYPE=C3D10M, ELSET=EBASE
*ELGEN, ELSET=E1
*ELGEN, ELSET=E2
*ELGEN, ELSET=E3
*ELGEN, ELSET=E4
*ELGEN, ELSET=E5
*ELSET, ELSET=EBASE
 E1, E2, E3, E4, E5
*ELEMENT, TYPE=C3D8R, ELSET=EBLOCK1
*ELEMENT, TYPE=C3D8R, ELSET=EBLOCK2
*ELSET, ELSET=EB1
*ELSET, ELSET=EB2
*ELSET, ELSET=ET1
*ELSET, ELSET=ET2
*ELSET, ELSET=EALL
 EBLOCK1,
 EBLOCK2,
*ELSET, ELSET=ETOP
 E1, E5
*SOLID SECTION, ELSET=EBASE, MAT=MAT1
*MATERIAL,NAME=MAT1
*ELASTIC
 3e6, 0.0
*SOLID SECTION, ELSET=EALL, MAT=MAT2
*MATERIAL,NAME=MAT2
*ELASTIC
 3e6, 0.0
*SURFACE, NAME=SBASE, TRIM=NO
 ETOP, S3,
*SURFACE, NAME=SBLOCK1, TRIM=YES
 EB1, S1
*SURFACE, NAME=SBLOCK2, TRIM=YES
 EB2, S1
*SURFACE, NAME=SECT1, TRIM=YES
 ET1, S1
*SURFACE, NAME=SECT2, TRIM=YES
 ET2, S1
*SURFACE, NAME=LOAD1, TRIM=YES
 ET1, S2
*SURFACE, NAME=LOAD2, TRIM=YES
 ET2, S2
*CONTACT PAIR, INTERACTION=NOSEP, EXTENSION ZONE=0.10
 SBLOCK1, SBASE
*CONTACT PAIR, INTERACTION=NOSEP, EXTENSION ZONE=0.10
 SBLOCK2, SBASE
*SURFACE INTERACTION, NAME=NOSEP
 1.0,
*SURFACE BEHAVIOR, NO SEPARATION
**
** STEP 1: Establish contact
**
*STEP, NLGEOM, INC=10
 STEP 1 -- Establish Contact
*STATIC
 0.0001,0.0001
*BOUNDARY, OP=NEW
 NBASE,  1, 3,  0.0
 NMOVE1, 1, 1,  0.0
 NMOVE2, 1, 1,  0.0
 NMOVE1, 2, 2, -0.001
 NMOVE2, 2, 2, -0.001
 NMOVE1, 3, 3,  0.0
 NMOVE2, 3, 3,  0.0
*PRINT, CONTACT=NO
*NODE PRINT, F=0
*EL PRINT, F=0
*CONTACT PRINT, FREQ=0
*CONTACT FILE, FREQ=0
 CSTRESS, CDISP,
*END STEP
**
** STEP 2: Apply pressure on blocks
**
*STEP, NLGEOM, INC=10
 STEP 2 -- Apply pressure on blocks
*STATIC
 0.0001,0.0001
*BOUNDARY, OP=NEW
 NBASE,  1, 3, 0.0
 NMOVE1, 1, 1, 0.0
 NMOVE2, 1, 1, 0.0
 NMOVE1, 3, 3, 0.0
 NMOVE2, 3, 3, 0.0
*DSLOAD, OP=NEW
 LOAD1, P, 100.0
 LOAD2, P, 200.0
*PRINT, CONTACT=YES
*NODE PRINT, F=10, NSET=NBLOCK1
 U,
*NODE PRINT, F=10, NSET=NBLOCK2
 U,
*NODE FILE, F=10, NSET=NBLOCK1
 U,
*NODE FILE, F=10, NSET=NBLOCK2
 U,
*EL PRINT, ELSET=EBLOCK1, F=10
 SINV,
*EL PRINT, ELSET=EBLOCK2, F=10
 SINV,
*EL FILE, ELSET=EBLOCK1, F=10
 SINV,
*EL FILE, ELSET=EBLOCK2, F=10
 SINV,
*CONTACT PRINT, SLAVE=SBLOCK1, MASTER=SBASE, FREQ=1
*CONTACT PRINT, SLAVE=SBLOCK2, MASTER=SBASE, FREQ=1
*CONTACT FILE,  SLAVE=SBLOCK1, MASTER=SBASE, FREQ=1
 CSTRESS, CDISP,
*CONTACT FILE,  SLAVE=SBLOCK2, MASTER=SBASE, FREQ=1
 CSTRESS, CDISP,
*OUTPUT, FIELD, FREQUENCY=1
*NODE OUTPUT, NSET=NBLOCK1
 U,
*NODE OUTPUT, NSET=NBLOCK2
 U,
*CONTACT OUTPUT, SLAVE=SBLOCK1, MASTER=SBASE, VARIABLE=PRESELECT
*CONTACT OUTPUT, SLAVE=SBLOCK2, MASTER=SBASE, VARIABLE=PRESELECT
*END STEP
**
** STEP 3: Slide block 1 in pos 1-dir by 8 units
**         Slide block 2 in neg 1-dir by 8 units
**
*STEP, NLGEOM, INC=20
 STEP 3 -- Slide blocks
*STATIC
 1.0,10.0
*BOUNDARY, TYPE=VELOCITY, OP=MOD
 NMOVE1, 1, 1,  0.8
 NMOVE2, 1, 1, -0.8
*END STEP
**
** STEP 4: Slide block 1 in pos 3-dir by 4 units
**         Slide blcok 2 diagonally
**
*STEP, NLGEOM, INC=20
 STEP 4 -- Slide blocks
*STATIC
 1.0,10.0
*BOUNDARY, TYPE=VELOCITY, OP=MOD
 NMOVE1, 1, 1, 0.0
 NMOVE1, 3, 3, 0.4
 NMOVE2, 1, 1, 0.8
 NMOVE2, 3, 3,-0.4
*END STEP
**
** STEP 5: Slide block 1 diagonally
**         Slide block 2 inin pos 3-dir by 4 units
**
*STEP, NLGEOM, INC=20
 STEP 5 -- Slide blocks diagonally to original position
*STATIC
 1.0, 10.0
*BOUNDARY, TYPE=VELOCITY, OP=MOD
 NMOVE1, 1, 1, -0.8
 NMOVE1, 3, 3, -0.4
 NMOVE2, 1, 1,  0.0
 NMOVE2, 3, 3,  0.4
*END STEP