*HEADING
SIMPLE CONTACT TEST, COUPLED THERMAL-ELECTRICAL-STRUCTURAL
*RESTART,WRITE
*NODE
*NSET,NSET=BOTTOM
*NSET,NSET=TOP
*NSET,NSET=BACK
*ELEMENT, TYPE=Q3D4, ELSET=SOLID
*SOLID SECTION, ELSET=SOLID, MATERIAL=SOLMAT
*MATERIAL, NAME=SOLMAT
*ELASTIC
 30.E6,0.3
*EXPANSION
 10.E-6,
*CONDUCTIVITY
 43.0,
*ELECTRICAL CONDUCTIVITY
*JOULE HEAT FRACTION
*DENSITY
 7700.,
*SPECIFIC HEAT
 600.,
*SURFACE INTERACTION, NAME=GAP
 1.0,
*GAP CONDUCTANCE
*GAP RADIATION
*GAP HEAT GENERATION
*GAP ELECTRICAL CONDUCTANCE
**
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16,STEFANB=5.6697E-8
**
*** nset definitions for prescribing temperatures.
**
*NSET,NSET=NELEM15,GEN
*NSET, NSET=NELEM620
*NSET, NSET=NELEM620, GEN
**
** surfaces
**
*surface, NAME=ASURF
 3, S1
 4, S1
*surface, NAME=BSURF
 6, S1
 7, S1
 11,S1
 12,S1
*CONTACT PAIR, SMALL SLIDING, INTERACTION=GAP, TYPE=SURFACE TO SURFACE
 BSURF, ASURF
*surface, NAME=ATIED
 8, S1
 9, S1
 13,S1
 14,S1
*surface, NAME=BTIED
 16,S1
 17,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
**
*AMPLITUDE,name=STEP
 0,1,100,1
*STEP
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE, ELECTRICAL
 1.,1.
*ENERGY FILE
*EL FILE
 S,
 E,
 EE,
 TEMP,
*NODE FILE
 RF,
 EPOT,NT,RFL
*ENERGY PRINT
*EL PRINT ,SUM=NO
 S,
 E,
 EE,
TEMP,
*NODE PRINT,SUM=NO
 RF,
 EPOT,NT,RFL
 U,
*BOUNDARY
 BACK,3,3
 BOTTOM,2,2
 TOP,2,2
 1, 1, 1
 13,1, 1
 266,1, 1
 278,1, 1
 NELEM15, 9, 9, 700.0
 NELEM620, 9, 9, 500.0
 NELEM15, 11, 11, 700.0
 NELEM620, 11, 11, 500.0
*CONTACT PRINT
 CSTRESS,CDISP
 ECD,HFL,HFLA,HTL
 SFDR,SFDRA
*CONTACT FILE
 CSTRESS,CDISP
 ECD,HFL,HFLA,HTL,
 SFDR,SFDRA,
*OUTPUT, FIELD, VARIABLE=ALL
*END STEP
**
**** STEP 2. Push top block to the right.
**
*STEP, NLGEOM ,INC=10
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE, ELECTRICAL
 0.1,1.0,0.1,0.1
*CHANGE FRICTION,INTERACTION=GAP,amplitude=step
*FRICTION
 0.01,
*BOUNDARY, OP=NEW
 BACK,3,3
 123, 1, 1, 0.01
 135, 1, 1, 0.01
 261, 1, 1, 0.01
 273, 1, 1, 0.01
 266, 1, 1, 0.01
 278, 1, 1, 0.01
 NELEM15, 9, 9, 700.0
 NELEM620, 9, 9, 500.0
 NELEM15, 11, 11, 700
 NELEM620, 11, 11, 500
 TOP, 2, 2
 BOTTOM, 2, 2
 1, 1, 1
 13,1, 1
*EL FILE
 S,
 E,
*NODE FILE
 RF,
 EPOT,NT,RFL,
*CONTACT PRINT
 CSTRESS,CDISP
 ECD,HFL,HFLA,HTL,
 SFDR,SFDRA,
*CONTACT FILE
 CSTRESS,CDISP
 ECD,HFL,HFLA,HTL,
 SFDR,SFDRA,
*OUTPUT, FIELD, VARIABLE=ALL
*END STEP
**
*** STEP 3. Push top block back into position.
**
*STEP, NLGEOM, INC=10
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE,ELECTRICAL
 0.1,1.,0.1,0.1
*BOUNDARY, OP=NEW
 BACK,3,3
 123, 1, 1,
 135, 1, 1,
 261, 1, 1,
 273, 1, 1,
 266, 1, 1,
 278, 1, 1,
 NELEM15, 9, 9, 700.0
 NELEM620, 9, 9, 500.0
 NELEM15, 11, 11, 700
 NELEM620, 11, 11, 500
 TOP, 2, 2
 BOTTOM, 2, 2
 1, 1, 1
 13,1, 1
*EL FILE
 S,
 E,
*NODE FILE
 RF,
 EPOT,NT, RFL,
*CONTACT PRINT
 CSTRESS,CDISP
 ECD,HFL,HFLA,HTL,
 SFDR,SFDRA,
*CONTACT FILE
 CSTRESS,CDISP
 ECD,HFL,HFLA,HTL,
 SFDR,SFDRA,
*OUTPUT, FIELD, VARIABLE=ALL
*END STEP
**
**** test for linear perturbation
***
*STEP,NLGEOM,PERTURBATION
*STATIC
*BOUNDARY,OP=NEW
 BACK,3,3
 BOTTOM,2,2
 TOP ,2,2 , 0.05
 1, 1, 1
 13,1, 1
 266,1, 1
 278,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,
 EPOT,NT,
*CONTACT PRINT
*CONTACT FILE
*OUTPUT, FIELD, VARIABLE=ALL
*END STEP