*HEADING
 SIMPLE CONTACT TEST, COUPLED TEMPERATURE-DISPLACEMENT ()
 SURFACE-BASED TIE CONSTRAINT IS USED
 CPEG8T ELEMENT
*NODE
 1,1.,5.
 3,2.,5.
 9,2.,10.
 7,1.,10.
 11,1.,10.0000001
 15,2.,10.0000001
 21,1.,12.5000001
 25,2.,12.5000001
 31,1.,12.5000001
 35,2.,12.5000001
 41,1.,15.0000001
 45,2.,15.0000001
*NODE,NSET=XTRA
 998,0.,0.
*NGEN,NSET=NALL
 1,3
 3,9,3
 7,9
 1,7,3
 11,15
 21,25
 11,21,5
 15,25,5
 16,20,2
 31,35,2
 41,45,2
 31,41,5
 35,45,5
*NSET,NSET=BOTTOM
 1,2,3
*NSET,NSET=TOP
 41,43,45
*NSET,NSET=LEFT
 1,4,7,11,16,31,36,41
*ELEMENT,TYPE=CPEG8T,ELSET=SOLID
 1,1,3,9,7,2,6,8,4
 2,11,13,23,21,12,18,22,16
 3,13,15,25,23,14,20,24,18
 4,31,35,45,41,33,40,43,36
*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, 15,2
 31,35,4
 41,45,4
**
** SURFACES
**
*SURFACE, NAME=ASURF
 1, S3
*SURFACE, NAME=BSURF
 2, S1
 3, S1
*CONTACT PAIR, SMALL SLIDING, INTERACTION=GAP
 BSURF, ASURF
*SURFACE, NAME=ATIED1
 2, S3
*SURFACE, NAME=ATIED2
 3, S3
*SURFACE, NAME=BTIED
 4, S1
*TIE,NAME=TIE,POSITION TOLERANCE=0.01
 ATIED1, BTIED
 ATIED2, BTIED
**
*** APPLY DIFFERENT TEMPERATURES TO THE NODES ABOVE AND
*** BELOW THE CONTACT SURFACE
**
*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
 BOTTOM,2,2
 TOP,2,2
 NELEM1, 11, 11, 700.0
 NELEM23, 11, 11, 500.0
 1,1,1,0.0
 41,1,1,0.0
 998,3,5,0.
*CONTACT PRINT,MASTER=ASURF,SLAVE=BSURF
 CSTRESS,CDISP
 HFL,HFLA,HTL
 SFDR,SFDRA
*CONTACT FILE,MASTER=ASURF,SLAVE=BSURF
 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
*BOUNDARY
 11, 1, 1, 0.01
 16, 1, 1, 0.01
 31, 1, 1, 0.01
 36, 1, 1, 0.01
 41, 1, 1, 0.01
 NELEM1, 11, 11, 700
 NELEM23, 11, 11, 500
 TOP, 2, 2
 BOTTOM, 2, 2
 998,3,5,0.
*BOUNDARY,FIXED
 1, 1, 1
 4, 1, 1
 7, 1, 1
*EL FILE
 S,
 E,
*NODE FILE
 RF,
 NT, RFL
*CONTACT PRINT,MASTER=ASURF,SLAVE=BSURF
 CSTRESS,CDISP
 HFL,HFLA,HTL
 SFDR, SFDRA
*CONTACT FILE,MASTER=ASURF,SLAVE=BSURF
 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,0.1,0.1
*BOUNDARY, OP=NEW
 NELEM1, 11, 11, 700
 NELEM23, 11, 11, 500
 TOP, 2, 2
 BOTTOM, 2, 2
 1,  1, 1, 0.0
 4,  1, 1, 0.0 
 7,  1, 1, 0.0 
 11, 1, 1, 0.0
 16, 1, 1, 0.0
 31, 1, 1, 0.0
 36, 1, 1, 0.0
 41, 1, 1, 0.0
 998,3,5,0.
*EL FILE
 S,
 E,
*NODE FILE
 RF,
 NT, RFL
*CONTACT PRINT,MASTER=ASURF,SLAVE=BSURF
 CSTRESS,CDISP
 HFL,HFLA,HTL
 SFDR,SFDRA
*CONTACT FILE,MASTER=ASURF,SLAVE=BSURF
 CSTRESS,CDISP
 HFL,HFLA,HTL
 SFDR,SFDRA
*END STEP
**
**** TEST FOR LINEAR PERTURBATION
***
*STEP,NLGEOM,PERTURBATION
*STATIC
*BOUNDARY,OP=NEW
 BOTTOM,2,2
 TOP ,2,2 , 0.05
 1,  1, 1, 0.0
 41, 1, 1, 0.0
 998,3,5,0.
*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,MASTER=ASURF,SLAVE=BSURF
*CONTACT FILE,MASTER=ASURF,SLAVE=BSURF
*END STEP