*HEADING
 : Adjust contact surface normal at symmetry plane
 Small sliding example with C3D8 elements
 Global coordinate system
*NODE, NSET=NROT
       1,         1.5,          0.,
       2,          2.,          0.,
       3,         2.5,          0.,
       4,     1.38582,    0.574025,
       5,     1.84776,    0.765367,
       6,      2.3097,    0.956709,
       7,     1.06066,     1.06066,
       8,     1.41421,     1.41421,
       9,     1.76777,     1.76777,
      10,    0.574025,     1.38582,
      11,    0.765366,     1.84776,
      12,    0.956708,      2.3097,
      13, -6.55671E-8,         1.5,
      14, -8.74228E-8,          2.,
      15, -1.09278E-7,         2.5,
      16,         1.5,          0.,          1.
      17,          2.,          0.,          1.
      18,         2.5,          0.,          1.
      19,     1.38582,    0.574025,          1.
      20,     1.84776,    0.765367,          1.
      21,      2.3097,    0.956709,          1.
      22,     1.06066,     1.06066,          1.
      23,     1.41421,     1.41421,          1.
      24,     1.76777,     1.76777,          1.
      25,    0.574025,     1.38582,          1.
      26,    0.765367,     1.84776,          1.
      27,    0.956708,      2.3097,          1.
      28, -6.55671E-8,         1.5,          1.
      29, -8.74228E-8,          2.,          1.
      30, -1.09278E-7,         2.5,          1.
      31,         2.6,          0.,
      32,        3.05,          0.,
      33,         3.5,          0.,
      34,     2.40209,    0.994976,
      35,     2.81783,     1.16718,
      36,     3.23358,     1.33939,
      37,     1.83848,     1.83848,
      38,     2.15668,     2.15668,
      39,     2.47487,     2.47487,
      40,    0.994978,     2.40209,
      41,     1.16719,     2.81783,
      42,     1.33939,     3.23358,
      43,  -1.1365E-7,         2.6,
      44,  -1.3332E-7,        3.05,
      45,  -1.5299E-7,         3.5,
      46,         2.6,          0.,          1.
      47,        3.05,          0.,          1.
      48,         3.5,          0.,          1.
      49,     2.40209,    0.994977,          1.
      50,     2.81783,     1.16718,          1.
      51,     3.23358,     1.33939,          1.
      52,     1.83848,     1.83848,          1.
      53,     2.15668,     2.15668,          1.
      54,     2.47487,     2.47487,          1.
      55,    0.994977,     2.40209,          1.
      56,     1.16718,     2.81783,          1.
      57,     1.33939,     3.23358,          1.
      58,  -1.1365E-7,         2.6,          1.
      59,  -1.3332E-7,        3.05,          1.
      60,  -1.5299E-7,         3.5,          1.
*ELEMENT, TYPE=C3D8, ELSET=EALL
       1,       1,       2,       5,       4,      16,      17,
      20,      19
       2,       2,       3,       6,       5,      17,      18,
      21,      20
       3,       4,       5,       8,       7,      19,      20,
      23,      22
       4,       5,       6,       9,       8,      20,      21,
      24,      23
       5,       7,       8,      11,      10,      22,      23,
      26,      25
       6,       8,       9,      12,      11,      23,      24,
      27,      26
       7,      10,      11,      14,      13,      25,      26,
      29,      28
       8,      11,      12,      15,      14,      26,      27,
      30,      29
       9,      31,      32,      35,      34,      46,      47,
      50,      49
      10,      32,      33,      36,      35,      47,      48,
      51,      50
      11,      34,      35,      38,      37,      49,      50,
      53,      52
      12,      35,      36,      39,      38,      50,      51,
      54,      53
      13,      37,      38,      41,      40,      52,      53,
      56,      55
      14,      38,      39,      42,      41,      53,      54,
      57,      56
      15,      40,      41,      44,      43,      55,      56,
      59,      58
      16,      41,      42,      45,      44,      56,      57,
      60,      59
*SOLID SECTION, ELSET=EALL, MATERIAL=MAT1
 1.0, 
*MATERIAL, NAME=MAT1
*ELASTIC
 3000.0,0.2
*ELSET, ELSET=EIN
       2,       4,       6,        8
*ELSET, ELSET=EOUT
       9,      11,      13,       15
*NSET, NSET=NXSYMM
      13,      14,      15,      28,      29,      30,      43,      44,
      45,      58,      59,       60
*NSET, NSET=NYSYMM
       1,       2,       3,      16,      17,      18,      31,      32,
      33,      46,      47,       48
*ELSET, ELSET=EPRESS
      10,      12,      14,       16
*surface, NAME=MASTER
 EIN,
*surface, NAME=SLAVE
 EOUT,
*CONTACT PAIR,INTERACTION=NONE,SMALL SLIDING
 SLAVE,MASTER
*SURFACE INTERACTION, NAME=NONE
*BOUNDARY
 NXSYMM, XSYMM
 NYSYMM, YSYMM
 NXSYMM, 3,3
 NYSYMM, 3,3
**
** STEP 1: Apply pressure load
**
*STEP, NLGEOM, INC=50
*STATIC
 ,1.0
*DLOAD
 EPRESS,P4,100.0
*PRINT, CONTACT=YES
*NODE PRINT, F=0
*EL PRINT, F=0
*CONTACT PRINT, FREQ=100
*CONTACT FILE, MASTER=MASTER, SLAVE=SLAVE, FREQ=100
 CSTRESS, CDISP
*END STEP
*STEP, NLGEOM, INC=50
**
** STEP 2 : Remove pressure load
**
*STATIC
 ,1.0
*DLOAD, OP=NEW
 EPRESS,P4,0.0
*PRINT, CONTACT=YES
*NODE PRINT, F=0
*EL PRINT, F=0
*CONTACT PRINT, FREQ=100
*CONTACT FILE, MASTER=MASTER, SLAVE=SLAVE, FREQ=100
 CSTRESS, CDISP
*END STEP