*HEADING
SIMPLE CONTACT TEST, COUPLED TEMPERATURE-DISPLACEMENT ()
CPEG3HT 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=CPEG3HT,ELSET=SOLID
1,  1,  3,  7
2,  3,  9,  7
3, 11, 12, 14
4, 12, 15, 14
5, 12, 13, 15
6, 13, 16, 15
7,114,116, 17
8,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
2, S2
*SURFACE, NAME=BSURF
3, S1
5, S1
*CONTACT PAIR, SMALL SLIDING, INTERACTION=GAP, TYPE=SURFACE TO SURFACE
BSURF, ASURF
*SURFACE, NAME=ATIED
4, S2
6, S2
*SURFACE, NAME=BTIED
7, S1
*CONTACT PAIR, SMALL SLIDING, TIED, ADJUST=0.01, INTERACTION=GAP, TYPE=SURFACE TO SURFACE
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 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