*HEADING
: CONTACT PAIR WITH LOAD CONTROL FOR ELASTIC MATERIAL 
          WITH DAMAGE, NORMAL ONLY BEHAVIOR [GK3D2N]
*NODE,NSET=NALL
1,0.,0.,1.
2,0.,1.,1.
*ELEMENT,TYPE=GK3D2N,ELSET=EINT
1,1,2
*GASKET SECTION,ELSET=EINT,BEHAVIOR=IS
,0.1,,1.1,0.,1.,0.
*GASKET BEHAVIOR,NAME=IS
*GASKET THICKNESS BEHAVIOR,TYPE=DAMAGE,DIRECTION=LOADING
0.,0.
10.71,0.052
28.57,0.087
46.43,0.139
82.14,0.217
164.3,0.287
*GASKET THICKNESS BEHAVIOR,TYPE=DAMAGE,DIRECTION=UNLOADING
0.,0.,0.139
14.28,0.122,0.139
46.43,0.139,0.139
0.,0.,0.217
12.,0.17,0.217
40.,0.21,0.217
82.14,0.217,0.217
0.,0.,0.287
11.0,0.24,0.287
164.3,0.287,0.287
*EXPANSION
1.e-2,
*******************************************
*AMPLITUDE,NAME=HIST
0.0,0.0, 1.0,-30., 2.0,-35., 3.0,-0., 
4.0,-35., 5.0,-10., 6.0,-0., 7.0,-37., 
8.0,-40., 9.0,-0.0, 10.0,-40., 11.0,-42., 
12.0,-0.0, 13.0,-42., 14.0,-44., 15.0,-0.0, 
16.0,-44., 17.0,-46.43, 18.0,-0., 19.0,-46.43,
20.0,-82.14, 21.0,-0.0, 22.0,-82.14, 23.0,-85., 
24.0,-0., 25.0,-85., 26.0,-100., 27.0,-0.0, 
28.0,-130., 29.0,-170., 30.0,-0.
*BOUNDARY
1,1,1
*******************************************
*NODE
10000,0.,5.,1.
1001,-1.,1.,0.
1002,-1.,1.,2.
1003,1.,1.,2.
1004,1.,1.,0.
*ELEMENT,TYPE=R3D4,ELSET=RBODY
100,1001,1002,1003,1004
*RIGID BODY,ELSET=RBODY,REFNODE=10000
*surface,type=node,NAME=CNSET
2,1.0
*surface,NAME=RSURF
RBODY,SNEG
*CONTACT PAIR,INTERACTION=INT1,TIED,ADJUST=0.01
CNSET,RSURF
*SURFACE INTERACTION,NAME=INT1
1.,
*******************************************
*STEP,INC=10000
Test Thickness Behavior of Gasket
*STATIC
0.01,30.,,0.01
*PRINT, FREQUENCY=100
*OUTPUT, HISTORY, FREQUENCY=100
*BOUNDARY
10000,1,1
10000,3,6
*CLOAD,AMP=HIST
10000,2,1.
*NODE PRINT,FREQ=0
U, 
RF, 
*NODE FILE,FREQ=148
U, 
RF, 
*EL PRINT,FREQ=0,ELSET=EINT
S, 
E, 
NE, 
EE, 
EVOL, 
THE, 
TEMP, 
*EL FILE,FREQ=148,ELSET=EINT
E, 
EE, 
NE, 
S, 
EVOL, 
THE, 
*OUTPUT,FIELD,FREQ=148,VAR=PRE
*CONTACT PRINT,MASTER=RSURF,SLAVE=CNSET,FREQ=0
*CONTACT FILE,MASTER=RSURF,SLAVE=CNSET,FREQ=148
*END STEP
*******************************************
**No Test of Membrane Behavior of Element
*******************************************
**No Test Transverse Shear Stiffness of Element
*******************************************
*STEP
Thermal Expansion of Element
*STATIC
*BOUNDARY,OP=NEW
10000,1,6
1,1
*TEMPERATURE,OP=NEW
NALL,100.
*END STEP
*******************************************