*HEADING
 SIMPLE CONTACT TEST, COUPLED TEMPERATURE-DISPLACEMENT ()
 CPEG4RHT ELEMENT
*RESTART, WRITE
*NODE
*NODE,NSET=XTRA
*NGEN,NSET=NALL
*NSET,NSET=BOTTOM
*NSET,NSET=TOP
*NSET,NSET=SIDE
*ELEMENT,TYPE=CPEG4RHT,ELSET=SOLID
*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
*GAP RADIATION
*GAP HEAT GENERATION
**
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16,STEFANB=5.6697E-8
**
*** NSET DEFINITIONS FOR PRESCRIBING TEMPERATURES.
**
*NSET, NSET=NELEM1
*NSET, NSET=NELEM23, GEN
*NSET,NSET=NALL
**
** SURFACES
**
*SURFACE, NAME=ASURF
 1, S3
*SURFACE, NAME=BSURF
 2, S1
 3, S1
*CONTACT PAIR, SMALL SLIDING, INTERACTION=GAP
 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
*ENERGY FILE
*EL FILE
 S,
 E,
 EE,
 TEMP,
*NODE FILE,NSET=NALL
 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
 1, 1, 1
 17, 1, 1
 NELEM1, 11, 11, 700.0
 NELEM23, 11, 11, 500.0
 998,3,5,0.
*CONTACT CONTROLS, FRICTION ONSET=DELAYED
*CONTACT PRINT
 CSTRESS,CDISP
 HFL,HFLA,HTL
 SFDR,SFDRA,
 CAREA,CFN,
*CONTACT FILE
 CSTRESS,CDISP
 HFL,HFLA,HTL
 SFDR,SFDRA,
 CAREA,CFN,
*SECTION PRINT, NAME=BTIED, SURFACE=BTIED
*SECTION FILE, NAME=BTIED, SURFACE=BTIED
*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.
*EL FILE
 S,
 E,
*NODE FILE,NSET=NALL
 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,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.
*EL FILE
 S,
 E,
*NODE FILE,NSET=NALL
 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
 BOTTOM,2,2
 TOP ,2,2 , 0.05
 1, 1, 1
 17, 1, 1
 998,3,5,0.
*ENERGY FILE
*EL FILE
 S,
 E,
*NODE FILE,NSET=NALL
 RF,
*ENERGY PRINT
*EL PRINT ,SUM=NO
 S,E
 TEMP,
*NODE PRINT,SUM=NO
 RF,
 NT,
*CONTACT PRINT
*CONTACT FILE
*END STEP