*HEADING
SIMPLE CONTACT TEST, COUPLED TEMPERATURE-DISPLACEMENT ()
*RESTART, WRITE
*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
*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,3
*NSET,NSET=TOP
41,45
*NSET,NSET=LEFT
1,11,31,41
*ELEMENT,TYPE=CAX4T,ELSET=SOLID
1,1,3,9,7
2,11,13,23,21
3,13,15,25,23
4,31,35,45,41
*SOLID SECTION,ELSET=SOLID,MATERIAL=SOLMAT
*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.
** 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, 7, 9
*NSET, NSET=NELEM23, GEN
11, 15,2
31,35,2
41,45,2
**
** 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
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
EE,
TEMP,
*OUTPUT,FIELD
*NODE OUTPUT
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
*CONTACT PRINT
CSTRESS,CDISP
HFL,HFLA,HTL
SFDR,SFDRA
*OUTPUT,FIELD
*CONTACT OUTPUT,VARIABLE=PRESELECT
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.03
31, 1, 1, 0.03
41, 1, 1, 0.03
NELEM1, 11, 11, 700
NELEM23, 11, 11, 500
TOP, 2, 2
BOTTOM, 2, 2
*BOUNDARY,FIXED
1, 1, 1
7, 1, 1
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*OUTPUT,FIELD
*NODE OUTPUT
RF,
NT,RFL
*CONTACT PRINT
CSTRESS,CDISP
HFL,HFLA,HTL
SFDR, SFDRA
*OUTPUT,FIELD
*CONTACT OUTPUT,VARIABLE=PRESELECT
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
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*OUTPUT,FIELD
*NODE OUTPUT
RF,
NT,RFL
*CONTACT PRINT
CSTRESS,CDISP
HFL,HFLA,HTL
SFDR,SFDRA
*OUTPUT,FIELD
*CONTACT OUTPUT,VARIABLE=PRESELECT
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
*OUTPUT,HISTORY,FREQUENCY=1
*ENERGY OUTPUT,VARIABLE=PRESELECT
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*OUTPUT,FIELD
*NODE OUTPUT
RF,
*ENERGY PRINT
*EL PRINT ,SUM=NO
 S,E
TEMP, 
*NODE PRINT,SUM=NO
RF, 
NT, 
*CONTACT PRINT
*OUTPUT,FIELD
*CONTACT OUTPUT,VARIABLE=PRESELECT
*END STEP