*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=Cgax4rt,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, Stefan Boltzman=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, Nlgeom
*Coupled Temperature-Displacement,Steady State
*Output,History,Frequency=1
*Output,Field
*Element Output
 S,
 Le,
 Ee,
 Temp,
*Output,Field
*Node Output
 Rf,
 Nt,
 Rfl,
*Energy Print
*El Print ,Sum=No
 S, 
 Le, 
 Ee,
 Temp, 
*Node Print,Sum=No
 Rf, 
 Nt,Rfl
 U, 
*Boundary
 Bottom,2,2
 Bottom,5,5
 Top,2,2
 Top,5,5
 Nelem1, 11, 11, 700.0
 Nelem23, 11, 11, 500.0
*Contact Print
 Cstress,Cdisp
 Hfl,Hfla,Htl
 Sfdr,Sfdra
*Output,Field
*Contact Output
 Cstress,Cdisp
 Hfl,Hfla,Htl
 Sfdr,Sfdra
*El File
 S,
 Le,
 Ee,
 Temp, 
*Node File
 Rf,
 Nt, 
 Rfl, 
*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,
 Le,
*Output,Field
*Node Output
 Rf,
 Nt,Rfl
*Contact Print
 Cstress,Cdisp
 Hfl,Hfla,Htl
 Sfdr, Sfdra
*Output,Field
*Contact Output
 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
 Top, 5, 5
 Bottom, 2, 2
 Bottom, 5, 5
*Output,Field
*Element Output
 S,
 Le,
*Output,Field
*Node Output
 Rf,
 Nt,Rfl
*Contact Print
 Cstress,Cdisp
 Hfl,Hfla,Htl
 Sfdr,Sfdra
*Output,Field
*Contact Output
 Cstress,Cdisp
 Hfl,Hfla,Htl
 Sfdr,Sfdra
*End Step
**
**** Test For Linear Perturbation
***
*Step,Nlgeom,Perturbation
*Static
*Boundary,Op=New
 Bottom, 2, 2
 Bottom, 5, 5
 Top ,2,2 , 0.05
*Output,History,Frequency=1
*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
 Cstress,Cdisp
 Hfl,Hfla,Htl
 Sfdr,Sfdra
*End Step