*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
2,8,3
11,15
21,25
11,21,5
15,25,5
16,20
31,35,2
41,45,2
31,41,5
33,43,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=CGAX6MT,ELSET=SOLID
1, 1, 9, 7, 5, 8,4
2, 9, 1, 3, 5, 2,6
3,11,23,21,17,22,16
4,23,11,13,17,12,18
5,13,25,23,19,24,18
6,25,13,15,19,14,20
7,31,45,41,38,43,36
8,45,31,35,38,33,40
*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=NELEM12, GEN
1,9
*NSET, NSET=NELEM36, GEN
11,15
31,35,2
41,45,2
**
** surfaces
**
*surface, NAME=ASURF
1, S2
*surface, NAME=BSURF
4, S2
6, S2
*CONTACT PAIR, SMALL SLIDING, INTERACTION=GAP
BSURF, ASURF
*surface, NAME=ATIED
3, S2
5, S2
*surface, NAME=BTIED
8, S2
*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,
*EL FILE,FREQ=999
S,
E,
*NODE FILE,FREQ=999
RF, 
NT,RFL
U,
*BOUNDARY
BOTTOM,2,2
TOP,2,2
NELEM12, 11, 11, 700.0
NELEM36, 11, 11, 500.0
NALL,5,5
*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
16, 1, 1, 0.03
31, 1, 1, 0.03
36, 1, 1, 0.03
41, 1, 1, 0.03
NELEM12, 11, 11, 700
NELEM36, 11, 11, 500
TOP, 2, 2
BOTTOM, 2, 2
NALL,5,5
*BOUNDARY,FIXED
1, 1, 1
4, 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
NELEM12, 11, 11, 700
NELEM36, 11, 11, 500
TOP, 2, 2
BOTTOM, 2, 2
NALL,5,5
*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
NALL,5,5
*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