*HEADING
 VERIFICATION TEST FOR FINITE SLIDING CONTACT PATCH ALGORITHM
 USING C3D10M ELEMENTS
*RESTART, WRITE, FREQ=1
*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=C3D10M, ELSET=EBASE
 1,    1,  203,  201, 1201,  102,  202,  101,  601,  702,  701 
 2,    1,    3,  203, 1003,    2,  103,  102,  502,  503,  603
 3,    1, 1201, 1001, 1003,  601, 1101,  501,  502, 1102, 1002
 4,    1, 1201, 1003,  203,  601, 1102,  502,  102,  702,  603
 5, 1003, 1201, 1203,  203, 1102, 1202, 1103,  603,  702,  703
*ELGEN, ELSET=E1
 1, 10, 2, 5, 1,,, 6, 1000, 200
*ELGEN, ELSET=E2
 2, 10, 2, 5, 1,,, 6, 1000, 200
*ELGEN, ELSET=E3
 3, 10, 2, 5, 1,,, 6, 1000, 200
*ELGEN, ELSET=E4
 4, 10, 2, 5, 1,,, 6, 1000, 200
*ELGEN, ELSET=E5
 5, 10, 2, 5, 1,,, 6, 1000, 200
*ELSET, ELSET=EBASE
 E1, E2, E3, E4, E5
*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=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
*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, TYPE=SURFACE
 SBLOCK1, SBASE				      
*CONTACT PAIR, INTERACTION=NOSEP, TYPE=SURFACE
 SBLOCK2, SBASE
*SURFACE INTERACTION, NAME=NOSEP
 1.0,
**
** STEP 1: Establish contact
**
*STEP, NLGEOM, INC=10, UNSYMM=YES
 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.1
 NMOVE2, 2, 2, -0.1
 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