*HEADING
 SIMPLE CONTACT TEST, COUPLED TEMPERATURE-DISPLACEMENT ()
 CPEG4RHT ELEMENT
*RESTART, WRITE
*NODE
 1,1.,5.
 3,2.,5.
 9,2.,10.
 7,1.,10.
 11,1.,10.0000001
 13,2.,10.0000001
 19,2.,15.0000001
 17,1.,15.0000001
 114,1.,12.5
 116,2.,12.5
*NODE,NSET=XTRA
 998,0.,0.
*NGEN,NSET=NALL
 1,3
 3,9,3
 7,9
 1,7,3
 2,8,3
 11,13
 13,19,3
 17,19
 11,17,3
 12,18,3
*NSET,NSET=BOTTOM
 1,2,3
*NSET,NSET=TOP
 17,18,19
*NSET,NSET=SIDE
 1,4,7,11,17
*ELEMENT,TYPE=CPEG4RHT,ELSET=SOLID
 1, 1,3,9,7
 2, 11,12,15,14
 3, 12,13,16,15
 4, 114,116,19,17
*SOLID SECTION,ELSET=SOLID,MATERIAL=SOLMAT,REFNODE=XTRA
*MATERIAL,NAME=SOLMAT
*ELASTIC
 30.E6,.3
*EXPANSION
 10.E-6, 
*CONDUCTIVITY
 43.0, 
*DENSITY
 7700., 
*SPECIFIC HEAT
 600., 
*SURFACE INTERACTION, NAME=GAP
 1.0, 
*FRICTION
 0.01, 
*GAP CONDUCTANCE
 1.,0.
 1.,1.
*GAP RADIATION
 3.4664e-2,3.4664e-2
 1.,0.
 1.,1.
*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, 7, 9
*NSET, NSET=NELEM23, GEN
 11, 13,1
 17,19,1
 114,114
 116,116
*NSET,NSET=NALL
 114,116
**
** SURFACES
**
*SURFACE, NAME=ASURF
 1, S3
*SURFACE, NAME=BSURF
 2, S1
 3, S1
*CONTACT PAIR, SMALL SLIDING, INTERACTION=GAP, NOTHICK
 BSURF, ASURF
*SURFACE, NAME=ATIED
 2, S3
 3, S3
*SURFACE, NAME=BTIED
 4, S1
*CONTACT PAIR, SMALL SLIDING, TIED, ADJUST=0.01, INTERACTION=GAP
 ATIED, BTIED
**
*** APPLY DIFFERENT TEMPERATURES TO THE NODES ABOVE AND
*** BELOW THE CONTACT SURFACE
**
*STEP
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
*BOUNDARY
 BOTTOM,2,2
 TOP,2,2
 1, 1, 1
 17, 1, 1
 NELEM1, 11, 11, 700.0
 NELEM23, 11, 11, 500.0
 998,3,5,0.
*OUTPUT, FIELD
*ELEMENT OUTPUT
 S,
*NODE OUTPUT
 U, 
*EL FILE
 S,
*NODE FILE,NSET=NALL
 U,
*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
*SOLUTION TECHNIQUE,TYPE=SEPARATED
*BOUNDARY, OP=NEW
 11, 1, 1, 0.01
 17, 1, 1, 0.01
 NELEM1, 11, 11, 700
 NELEM23, 11, 11, 500
 TOP, 2, 2
 BOTTOM, 2, 2
 1, 1, 1
 7, 1, 1
 114,1,1,0.01
 998,3,5,0.
*END STEP
**
*** STEP 3. PUSH TOP BLOCK BACK INTO POSITION.
**
*STEP, NLGEOM, INC=10
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
 0.1,1.0,0.1,0.1
*BOUNDARY, OP=NEW
 11, 1, 1
 17, 1, 1
 7, 1, 1
 NELEM1, 11, 11, 700
 NELEM23, 11, 11, 500
 TOP, 2, 2
 BOTTOM, 2, 2
 1, 1, 1
 114,1,1
 998,3,5,0.
*END STEP
**
**** TEST FOR LINEAR PERTURBATION
***
*STEP,NLGEOM,PERTURBATION
*STATIC
*BOUNDARY,OP=NEW
 BOTTOM,2,2
 TOP ,2,2 , 0.05
 1, 1, 1
 17, 1, 1
 998,3,5,0.
*OUTPUT, FIELD
*ELEMENT OUTPUT
 S,
*NODE OUTPUT
 U,
*EL FILE
 S,
*NODE FILE,NSET=NALL
 U,
*END STEP