*HEADING 
: CONTACT PAIR WITH DISPLACEMENT CONTROL FOR ELASTIC-PLASTIC MATERIAL 
          HIGH INITIAL STIFFNESS, INITIAL GAP [GK3D4LN]
** SLOPE DROP OPTION FOR GASKET
*NODE,NSET=NALL
1,    3.0,    0.0
2,    4.0,    0.0
3,    3.0,    1.0
4,    4.0,    1.0
*ELEMENT,TYPE=GK3D4LN,ELSET=EINT
1,1,2,3,4
*GASKET SECTION,ELSET=EINT,BEHAVIOR=IS
,0.1
*GASKET BEHAVIOR,NAME=IS
*GASKET ELASTICITY,COMPONENT=MEMBRANE
2000.,0.2
*GASKET ELASTICITY,COMPONENT=TRANSVERSE SHEAR,VAR=FORCE
200.,
*GASKET THICKNESS  BEHAVIOR,DIRECTION=LOADING,TYPE=ELASTIC-PLASTIC,VAR=FORCE,
SLOPE DROP=0.861111
0.,0.
20.0, 0.005
45.0, 0.05
70.0, 0.135
82.0, 0.22
83.0, 0.287
*GASKET THICKNESS BEHAVIOR,DIRECTION=UNLOADING,TYPE=ELASTIC-PLASTIC,VAR=FORCE
0.,   0.08,  0.08
14.,  0.12,  0.08
70.,  0.135, 0.08
0.0,  0.17,  0.17
25.0, 0.20,  0.17
82.0, 0.22,  0.17
0.0,  0.23,  0.23
36.0, 0.265, 0.23
83.0, 0.287, 0.23
*EXPANSION
1.e-2,
*rigid body,analytical surface=RIGID, REFNODE=10001 
*surface,NAME=RIGID,TYPE=CYLINDER
2.0, 1.0, +1., 5.0, 1.0, +1.0
2.0, 1.0,  0.0
START, 3.0, 0.0
LINE,  0.0, 0.0
*surface, NAME=DEFORM
EINT, SPOS
*CONTACT PAIR, INTERACTION=INTER, TIED,ADJUST=0.1
DEFORM, RIGID
*SURFACE INTERACTION, NAME=INTER
*******************************************
*AMPLITUDE,NAME=HIST
0.0,0.0, 1.0,-0.189, 2.0,-0.198, 3.0,-0., 
4.0,-0.198, 5.0,-0.205, 6.0,-0., 7.0,-0.205,
8.0,-0.31, 9.0,-0.0, 10.0,-0.31, 11.0,-0.32, 
12.0,-0.0, 13.0,-0.32, 14.0,-0.33, 15.0,-0.0, 
16.0,-0.33, 17.0,-0.339, 18.0,-0., 19.0,-0.339,
20.0,-0.417, 21.0,-0.0, 22.0,-0.417, 23.0,-0.440, 
24.0,-0., 25.0,-0.440, 26.0,-0.465, 27.0,-0.0, 
28.0,-0.47, 29.0,-0.495, 30.0,-0.
*BOUNDARY
1,1
2,1
10001,1,1
10001,3,6
*******************************************
*STEP,INC=10000
Test Normal Behavior of Gasket
*STATIC
0.01,20.,,0.01
*BOUNDARY,AMP=HIST
10001, 2, 2, 1.0
*NODE PRINT,FREQ=0
U, 
RF, 
*EL PRINT,FREQ=0
S, 
E, 
NE, 
EE, 
PE, 
EVOL, 
THE, 
TEMP, 
*OUTPUT, FIELD, FREQ=148
*NODE OUTPUT
U, 
RF, 
*ELEMENT OUTPUT
E, 
EE, 
NE, 
PE, 
S, 
EVOL, 
THE, 
*END STEP
******************************************
*STEP
Thermal Expansion of Element
*STATIC
*BOUNDARY,OP=MOD
10001,2,2,0.0
*TEMPERATURE,OP=NEW
NALL,100.
*END STEP
*******************************************