*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
13,2.,10.0000001
19,2.,15.0000001
17,1.,15.0000001
114,1.,12.5
116,2.,12.5
*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=CPE4RHT,ELSET=SOLID
1, 1,3,9,7
2, 11,12,15,14
3, 12,13,16,15
4, 114,116,19,17
*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,LAGRANGE
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, 13,1
17,19,1
114,114
116,116
**
** surfaces
**
*surface, NAME=ASURF
1, S3
*surface, NAME=BSURF
2, S1
3, S1
*CONTACT PAIR, SMALL SLIDING, INTERACTION=GAP, TYPE=SURFACE TO SURFACE
BSURF, ASURF
*surface, NAME=ATIED
2, S3
3, S3
*surface, NAME=BTIED
4, 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
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
*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
*EL FILE
S,
E,
*NODE FILE
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
*EL FILE
S,
E,
*NODE FILE
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
*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
*CONTACT FILE
*END STEP