*HEADING
 VERIFICATION TEST FOR FINITE SLIDING CONTACT PATCH ALGORITHM
 USING C3D8 ELEMENTS
*RESTART, WRITE, OVERLAY
*NODE, NSET=NB1
   1,   0.0,  0.0, 0.0
  21,  10.0,  0.0, 0.0
 101,   0.0,  0.5, 0.0
 121,  10.0,  0.5, 0.0
 201,   0.0,  1.0, 0.0
 221,  10.0,  1.0, 0.0
*NGEN, NSET=NB1
   1,  21, 1
 101, 121, 1
 201, 221, 1
*NCOPY, CHANGE NUMBER=500,  OLDSET=NB1, SHIFT, NEW SET=NB2
  0.0, 0.0, 0.5

*NCOPY, CHANGE NUMBER=1000, OLDSET=NB1, SHIFT, NEW SET=NB3
  0.0, 0.0, 1.0

*NCOPY, CHANGE NUMBER=1500, OLDSET=NB1, SHIFT, NEW SET=NB4
  0.0, 0.0, 1.5

*NCOPY, CHANGE NUMBER=2000, OLDSET=NB1, SHIFT, NEW SET=NB5
  0.0, 0.0, 2.0

*NCOPY, CHANGE NUMBER=2500, OLDSET=NB1, SHIFT, NEW SET=NB5
  0.0, 0.0, 2.5

*NCOPY, CHANGE NUMBER=3000, OLDSET=NB1, SHIFT, NEW SET=NB6
  0.0, 0.0, 3.0

*NCOPY, CHANGE NUMBER=3500, OLDSET=NB1, SHIFT, NEW SET=NB7
  0.0, 0.0, 3.5

*NCOPY, CHANGE NUMBER=4000, OLDSET=NB1, SHIFT, NEW SET=NB8
  0.0, 0.0, 4.0

*NCOPY, CHANGE NUMBER=4500, OLDSET=NB1, SHIFT, NEW SET=NB9
  0.0, 0.0, 4.5

*NCOPY, CHANGE NUMBER=5000, OLDSET=NB1, SHIFT, NEW SET=NB10
  0.0, 0.0, 5.0

*NCOPY, CHANGE NUMBER=5500, OLDSET=NB1, SHIFT, NEW SET=NB11
  0.0, 0.0, 5.5

*NCOPY, CHANGE NUMBER=6000, OLDSET=NB1, SHIFT, NEW SET=NB12
  0.0, 0.0, 6.0

*NSET, NSET=NBASE
 NB1, NB2, NB3, NB4, NB5, NB6, NB8, NB9, NB10, NB11, NB12
*SYSTEM
 0.5, 1.0, 0.5
*NODE, NSET=NBLOCK1
 251, 0.0, 0.0, 1.0
 252, 0.5, 0.0, 1.0
 253, 1.0, 0.0, 1.0
 254, 0.0, 0.0, 0.5
 255, 0.5, 0.0, 0.5
 256, 1.0, 0.0, 0.5
 257, 0.0, 0.0, 0.0
 258, 0.5, 0.0, 0.0
 259, 1.0, 0.0, 0.0
 261, 0.0, 0.5, 1.0
 262, 0.5, 0.5, 1.0
 263, 1.0, 0.5, 1.0
 264, 0.0, 0.5, 0.5
 265, 0.5, 0.5, 0.5
 266, 1.0, 0.5, 0.5
 267, 0.0, 0.5, 0.0
 268, 0.5, 0.5, 0.0
 269, 1.0, 0.5, 0.0
 271, 0.0, 1.0, 1.0
 272, 0.5, 1.0, 1.0
 273, 1.0, 1.0, 1.0
 274, 0.0, 1.0, 0.5
 275, 0.5, 1.0, 0.5
 276, 1.0, 1.0, 0.5
 277, 0.0, 1.0, 0.0
 278, 0.5, 1.0, 0.0
 279, 1.0, 1.0, 0.0
*SYSTEM
 8.5, 1.0, 4.5
*NODE, NSET=NBLOCK2
 351, 0.0, 0.0, 1.0
 352, 0.5, 0.0, 1.0
 353, 1.0, 0.0, 1.0
 354, 0.0, 0.0, 0.5
 355, 0.5, 0.0, 0.5
 356, 1.0, 0.0, 0.5
 357, 0.0, 0.0, 0.0
 358, 0.5, 0.0, 0.0
 359, 1.0, 0.0, 0.0
 361, 0.0, 0.5, 1.0
 362, 0.5, 0.5, 1.0
 363, 1.0, 0.5, 1.0
 364, 0.0, 0.5, 0.5
 365, 0.5, 0.5, 0.5
 366, 1.0, 0.5, 0.5
 367, 0.0, 0.5, 0.0
 368, 0.5, 0.5, 0.0
 369, 1.0, 0.5, 0.0
 371, 0.0, 1.0, 1.0
 372, 0.5, 1.0, 1.0
 373, 1.0, 1.0, 1.0
 374, 0.0, 1.0, 0.5
 375, 0.5, 1.0, 0.5
 376, 1.0, 1.0, 0.5
 377, 0.0, 1.0, 0.0
 378, 0.5, 1.0, 0.0
 379, 1.0, 1.0, 0.0
*NSET, NSET=NMOVE1
 261, 262, 263, 264, 265, 266, 267, 268, 269 
 271, 272, 273, 274, 275, 276, 277, 278, 279 
*NSET, NSET=NMOVE2
 361, 362, 363, 364, 365, 366, 367, 368, 369 
 371, 372, 373, 374, 375, 376, 377, 378, 379 
*NSET, NSET=NBASE, GEN
    1,   21, 1
  501,  521, 1
 1001, 1021, 1
 1501, 1521, 1
 2001, 2021, 1
 2501, 2521, 1
 3001, 3021, 1
 3501, 3521, 1
 4001, 4021, 1
 4501, 4521, 1
 5001, 5021, 1
 5501, 5521, 1
 6001, 6021, 1
 6501, 6521, 1
*NSET, NSET=NALL
 NBASE,
 NBLOCK1,
 NBLOCK2,
*ELEMENT, TYPE=C3D8, ELSET=EBASE
 1, 1001, 1003, 3, 1, 1201, 1203, 203, 201
*ELGEN, ELSET=EBASE
 1, 10, 2, 1,1,,,6,1000,10
*ELEMENT, TYPE=C3D8R, ELSET=EBLOCK1
 100, 251, 252, 255, 254, 261, 262, 265, 264
 101, 252, 253, 256, 255, 262, 263, 266, 265
 102, 254, 255, 258, 257, 264, 265, 268, 267
 103, 255, 256, 259, 258, 265, 266, 269, 268
 104, 261, 262, 265, 264, 271, 272, 275, 274
 105, 262, 263, 266, 265, 272, 273, 276, 275
 106, 265, 266, 269, 268, 275, 276, 279, 278
 107, 264, 265, 268, 267, 274, 275, 278, 277
*ELEMENT, TYPE=C3D8R, ELSET=EBLOCK2
 108, 351, 352, 355, 354, 361, 362, 365, 364
 109, 352, 353, 356, 355, 362, 363, 366, 365
 110, 354, 355, 358, 357, 364, 365, 368, 367
 111, 355, 356, 359, 358, 365, 366, 369, 368
 112, 361, 362, 365, 364, 371, 372, 375, 374
 113, 362, 363, 366, 365, 372, 373, 376, 375
 114, 365, 366, 369, 368, 375, 376, 379, 378
 115, 364, 365, 368, 367, 374, 375, 378, 377
*ELSET, ELSET=EALL
 EBLOCK1,
 EBLOCK2,
*ELSET, ELSET=EB1
 100, 101, 102, 103
*ELSET, ELSET=EB2
 108, 109, 110, 111
*ELSET, ELSET=ET1
 104, 105, 106, 107
*ELSET, ELSET=ET2
 112, 113, 114, 115
*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
 EBASE, S2,
*SURFACE, NAME=SBLOCK1, TRIM=YES
 EB1, S1
*SURFACE, NAME=SBLOCK2, TRIM=YES
 EB2, S1
*SURFACE, NAME=LOAD1, TRIM=YES
 ET1, S2
*SURFACE, NAME=LOAD2, TRIM=YES
 ET2, S2
*SURFACE, NAME=SECT1, TRIM=YES
 ET1, S1
*SURFACE, NAME=SECT2, TRIM=YES
 ET2, S1
*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
*SECTION FILE, NAME=SECT1, SURFACE=SECT1
 SOAREA, 
 SOF,
*SECTION PRINT, NAME=SECT1, SURFACE=SECT1
 SOAREA, 
 SOF,
*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