*HEADING
SIMPLE CONTACT TEST, COUPLED TEMPERATURE-DIPLACEMENT []
*RESTART,WRITE
*NODE
       1,          1.,          5.
       2,         1.5,          5.
       3,          2.,          5.
       4,          1.,         7.5
       5,          2.,         7.5
       6,          1.,         10.
       7,         1.5,         10.
       8,          2.,         10.
       9,          1.,          5.,         0.5
      10,          2.,          5.,         0.5
      11,          1.,         10.,         0.5
      12,          2.,         10.,         0.5
      13,          1.,          5.,          1.
      14,         1.5,          5.,          1.
      15,          2.,          5.,          1.
      16,          1.,         7.5,          1.
      17,          2.,         7.5,          1.
      18,          1.,         10.,          1.
      19,         1.5,         10.,          1.
      20,          2.,         10.,          1.
      21,          1.,        12.5
      22,          1.,        11.25
      23,          1.,        10.,
      24,        1.25,        12.5,
      25,        1.25,        10.,
      26,         1.5,        12.5
      27,         1.5,        11.25
      28,         1.5,        10.
      29,          1.,        12.5,         0.5
      30,          1.,        10.,          0.5
      31,         1.5,        12.5,         0.5
      32,         1.5,        10.,          0.5
      33,          1.,        12.5,          1.
      34,          1.,        11.25,         1.
      35,          1.,        10.,           1.
      36,        1.25,        12.5,          1.
      37,        1.25,        10.,           1.
      38,         1.5,        12.5,          1.
      39,         1.5,        11.25,         1.
      40,         1.5,        10.,           1.,
      44,        1.75,        12.5,
      45,        1.75,        10.
      46,          2.,        12.5
      47,          2.,        11.25
      48,          2.,        10.
      51,          2.,        12.5,          0.5
      52,          2.,        10.,           0.5
      56,        1.75,        12.5,          1.
      57,        1.75,        10.,           1.
      58,          2.,        12.5,          1.
      59,          2.,        11.25,         1.
      60,          2.,        10.,           1.
      61,          1.,        12.5
      62,         1.5,        12.5
      63,          2.,        12.5
      64,          1.,        13.75
      65,          2.,        13.75
      66,          1.,         15.
      67,         1.5,         15.
      68,          2.,         15.
      69,          1.,        12.5,         0.5
      70,          2.,        12.5,         0.5
      71,          1.,         15.,         0.5
      72,          2.,         15.,         0.5
      73,          1.,        12.5,          1.
      74,         1.5,        12.5,          1.
      75,          2.,        12.5,          1.
      76,          1.,        13.75,         1.
      77,          2.,        13.75,         1.
      78,          1.,         15.,          1.
      79,         1.5,         15.,          1.
      80,          2.,         15.,          1.
*NSET,NSET=BOTTOM
1,2,3,9,10,13,14,15,
*NSET,NSET=TOP
66,67,68,71,72,
78,79,80,
*NSET,NSET=BACK
1,2,3
66,67,68,
*ELEMENT, TYPE=C3D20T, ELSET=SOLID
       1,       1,       3,       8,       6,      13,      15,      20,
      18,       2,       5,       7,       4,      14,      17,      19,
      16,       9,      10,      12,      11
       2,      21,      23,      28,      26,      33,      35,      40,
      38,      22,      25,      27,      24,      34,      37,      39,
      36,      29,      30,      32,      31
       3,      26,      28,      48,      46,      38,      40,      60,
      58,      27,      45,      47,      44,      39,      57,      59,
      56,      31,      32,      52,      51
       4,      61,      63,      68,      66,      73,      75,      80,
      78,      62,      65,      67,      64,      74,      77,      79,
      76,      69,      70,      72,      71
*SOLID SECTION, ELSET=SOLID, MATERIAL=SOLMAT
*MATERIAL, NAME=SOLMAT
*ELASTIC
30.E6,0.3
*EXPANSION
10.E-6,
*CONDUCTIVITY
43.0, 
*DENSITY
7700., 
*SPECIFIC HEAT
600., 
*SURFACE INTERACTION, NAME=GAP
1.0,
*SURFACE BEHAVIOR, AUGMENTED LAGRANGE
*GAP CONDUCTANCE
1.,0.
1.,1.
*GAP RADIATION
3.4664e-2, 3.4664e-2
1.,0.
1.,1.
** converted from old data: 1.E-9, 1.E-9
*GAP HEAT GENERATION
0.5,0.5
**
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16,STEFANB=5.6697E-8
**
*** nset definitions for prescribing temperatures.
**
*NSET, NSET=NELEM1
1,3,13,15,
18,20,6,
8,
*NSET, NSET=NELEM23
35,40,60,
23,28,48,
73,75,63
61,
78,80,66,
68,
**
** surfaces
**
*surface, NAME=ASURF
1, S5
*surface, NAME=BSURF
2, S4
3, S4
*CONTACT PAIR, SMALL SLIDING, INTERACTION=GAP
BSURF, ASURF
*surface, NAME=ATIED
2, S6
3, S6
*surface, NAME=BTIED
4, S3
*CONTACT PAIR, SMALL SLIDING, TIED, ADJUST=0.01, INTERACTION=GAP
ATIED, BTIED
**
*** Apply different temperatures to the nodes above and
*** below the contact surface
**
*AMPLITUDE,name=STEP
0,1,100,1
*STEP
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
*ENERGY FILE
*EL FILE
 S,
 E,
 EE,
TEMP, 
*NODE FILE
RF,
NT,RFL
*ENERGY PRINT
*EL PRINT ,SUM=NO
 S,
 E,
 EE,
TEMP, 
*NODE PRINT,SUM=NO
RF, 
NT,RFL
U, 
*BOUNDARY
BACK,3,3
BOTTOM,2,2
TOP,2,2
1, 1, 1
9, 1, 1
13,1, 1
66,1, 1
71,1, 1
78,1, 1
NELEM1, 11, 11, 700.0
NELEM23, 11, 11, 500.0
*CONTACT PRINT
CSTRESS,CDISP
HFL,HFLA,HTL
SFDR,SFDRA
*CONTACT FILE
CSTRESS,CDISP
HFL,HFLA,HTL,
SFDR,SFDRA,
*END STEP
**
**** STEP 2. Push top block to the right.
**
*STEP, NLGEOM ,INC=10
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
0.1,1.0,0.1,0.1
*CHANGE FRICTION,INTERACTION=GAP,amplitude=step
*FRICTION
0.01, 
*BOUNDARY, OP=NEW
BACK,3,3
23, 1, 1, 0.01
30, 1, 1, 0.01
35, 1, 1, 0.01
22, 1, 1, 0.01
34, 1, 1, 0.01
61, 1, 1, 0.01
69, 1, 1, 0.01
73, 1, 1, 0.01
64, 1, 1, 0.01
76, 1, 1, 0.01
66, 1, 1, 0.01
71, 1, 1, 0.01
78, 1, 1, 0.01
NELEM1, 11, 11, 700
NELEM23, 11, 11, 500
TOP, 2, 2
BOTTOM, 2, 2
1, 1, 1
9, 1, 1
13,1, 1
4,1,1
16,1,1
6,1,1
11,1,1
18,1,1
*EL FILE
S,
E,
*NODE FILE
RF,
NT,RFL,
*CONTACT PRINT
CSTRESS,CDISP
HFL,HFLA,HTL,
SFDR,SFDRA,
*CONTACT FILE
CSTRESS,CDISP
HFL,HFLA,HTL,
SFDR,SFDRA,
*END STEP
**
*** STEP 3. Push top block back into position.
**
*STEP, NLGEOM, INC=10
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
0.1,1.,0.1,0.1
*BOUNDARY, OP=NEW
BACK,3,3
23, 1, 1,
30, 1, 1, 
35, 1, 1, 
22, 1, 1,
34, 1, 1, 
61, 1, 1, 
69, 1, 1, 
73, 1, 1, 
64, 1, 1,
76, 1, 1, 
66, 1, 1, 
71, 1, 1,
78, 1, 1,
NELEM1, 11, 11, 700
NELEM23, 11, 11, 500
TOP, 2, 2
BOTTOM, 2, 2
1, 1, 1
9, 1, 1
13,1, 1
4,1,1
16,1,1
6,1,1
11,1,1
18,1,1
*EL FILE
S,
E,
*NODE FILE
RF,
NT, RFL,
*CONTACT PRINT
CSTRESS,CDISP
HFL,HFLA,HTL,
SFDR,SFDRA,
*CONTACT FILE
CSTRESS,CDISP
HFL,HFLA,HTL,
SFDR,SFDRA,
*END STEP
**
**** test for linear perturbation
***
*STEP,NLGEOM,PERTURBATION
*STATIC
*BOUNDARY,OP=NEW
BACK,3,3
BOTTOM,2,2
TOP ,2,2 , 0.05
1, 1, 1
9, 1, 1
13,1, 1
66,1, 1
71,1, 1
78,1, 1
*ENERGY FILE
*EL FILE
 S,
 E,
*NODE FILE
 RF,
*ENERGY PRINT
*EL PRINT ,SUM=NO
 S,
 E,
TEMP, 
*NODE PRINT,SUM=NO
RF, 
NT, 
*CONTACT PRINT
*CONTACT FILE
*END STEP