*HEADING 
: CONTACT PAIR WITH LOAD CONTROL FOR ELASTIC-PLASTIC MATERIAL 
          HIGH INITIAL STIFFNESS [GK3D6LN]
** DEPENDENCY ON FIELD VARIABLE AND TEMPERATURE FOR GASKET ELEMENTS
*NODE,NSET=NALL
1, 3.0,  0.0
2, 3.5,  0.0
3, 4.0,  0.0
4, 3.0,  1.0
5, 3.5,  1.0
6, 4.0,  1.0
*ELEMENT,TYPE=GK3D6LN,ELSET=EINT
1,1,2,3,4,5,6 
*GASKET SECTION,ELSET=EINT,BEHAVIOR=IS
*GASKET BEHAVIOR,NAME=IS
*GASKET ELASTICITY,COMPONENT=MEMBRANE,DEP=1
2000.,0.2,0.0,50.0 
2000.,0.2,100.0,50.0 
2000.,0.2,0.0,150.0 
2000.,0.2,100.0,150.0 
*GASKET ELASTICITY,COMPONENT=TRANSVERSE SHEAR,VAR=FORCE,DEP=1
200.,0.0,50.0 
200.,100.0,50.0 
200.,0.0,150.0 
200.,100.0,150.0 
*GASKET THICKNESS  BEHAVIOR,DIRECTION=LOADING,TYPE=ELASTIC-PLASTIC,VAR=FORCE,
YIELD=0.0050,DEP=1
0. ,   0.   ,0.0,50.0 
20.,   0.005,0.0,50.0 
45.,   0.05 ,0.0,50.0 
70.,   0.135,0.0,50.0 
82.,   0.220,0.0,50.0 
83.,   0.287,0.0,50.0 
0. ,   0.   ,100.0,50.0 
20.,   0.005,100.0,50.0 
45.,   0.05 ,100.0,50.0 
70.,   0.135,100.0,50.0 
82.,   0.220,100.0,50.0 
83.,   0.287,100.0,50.0 
0. ,   0.   ,0.0,150.0 
20.,   0.005,0.0,150.0 
45.,   0.05 ,0.0,150.0 
70.,   0.135,0.0,150.0 
82.,   0.220,0.0,150.0 
83.,   0.287,0.0,150.0 
0. ,   0.   ,100.0,150.0 
20.,   0.005,100.0,150.0 
45.,   0.05 ,100.0,150.0 
70.,   0.135,100.0,150.0 
82.,   0.220,100.0,150.0 
83.,   0.287,100.0,150.0 
*GASKET THICKNESS BEHAVIOR,DIRECTION=UNLOADING,TYPE=ELASTIC-PLASTIC,VAR=FORCE,
DEP=1
0.00,   0.08,  0.08,0.0,50.0 
14.0,   0.12,  0.08,0.0,50.0 
70.0,   0.135, 0.08,0.0,50.0 
0.00,   0.170, 0.17,0.0,50.0 
25.0,   0.200, 0.17,0.0,50.0 
82.0,   0.220, 0.17,0.0,50.0 
0.00,   0.230, 0.23,0.0,50.0 
36.0,   0.265, 0.23,0.0,50.0 
83.0,   0.287, 0.23,0.0,50.0 
0.00,   0.08,  0.08,100.0,50.0 
14.0,   0.12,  0.08,100.0,50.0 
70.0,   0.135, 0.08,100.0,50.0 
0.00,   0.170, 0.17,100.0,50.0 
25.0,   0.200, 0.17,100.0,50.0 
82.0,   0.220, 0.17,100.0,50.0 
0.00,   0.230, 0.23,100.0,50.0 
36.0,   0.265, 0.23,100.0,50.0 
83.0,   0.287, 0.23,100.0,50.0 
0.00,   0.08,  0.08,0.0,150.0 
14.0,   0.12,  0.08,0.0,150.0 
70.0,   0.135, 0.08,0.0,150.0 
0.00,   0.170, 0.17,0.0,150.0 
25.0,   0.200, 0.17,0.0,150.0 
82.0,   0.220, 0.17,0.0,150.0 
0.00,   0.230, 0.23,0.0,150.0 
36.0,   0.265, 0.23,0.0,150.0 
83.0,   0.287, 0.23,0.0,150.0 
0.00,   0.08,  0.08,100.0,150.0 
14.0,   0.12,  0.08,100.0,150.0 
70.0,   0.135, 0.08,100.0,150.0 
0.00,   0.170, 0.17,100.0,150.0 
25.0,   0.200, 0.17,100.0,150.0 
82.0,   0.220, 0.17,100.0,150.0 
0.00,   0.230, 0.23,100.0,150.0 
36.0,   0.265, 0.23,100.0,150.0 
83.0,   0.287, 0.23,100.0,150.0 
*EXPANSION
1.e-2,
*INITIAL CONDITIONS, TYPE=FIELD, VARIABLE=1
NALL, 112.0
*INITIAL CONDITIONS, TYPE=TEMPERATURE
NALL,  0.0
*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,-20., 2.0,-35., 3.0,-0.1, 
4.0,-35., 5.0,-50., 6.0,-0.1, 7.0,-50., 
8.0,-65., 9.0,-0.1, 10.0,-65., 11.0,-80., 
12.0,-0.1, 13.0,-80., 14.0,-95., 15.0,-0.1, 
16.0,-95., 17.0,-110., 18.0,-0.1, 19.0,-110.,
20.0,-125., 21.0,-0.1, 22.0,-125., 23.0,-140., 
24.0,-0.1, 25.0,-140., 26.0,-155., 27.0,-0.1, 
28.0,-155., 29.0,-170., 30.0,-0.1
*BOUNDARY
1,1
2,1
3,1
10001,1,1
10001,3,6
*******************************************
*STEP,INC=10000
Test Normal Behavior of Gasket
*STATIC
0.01,20.,,0.01
*CONTROLS,PARAMETER=TIME INCREMENT

,0.25
*CLOAD,AMP=HIST
10001,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
*CONTROLS,RESET
*BOUNDARY,OP=NEW
1,1
2,1
3,1
10001,1,6
*CLOAD,OP=NEW
*TEMPERATURE,OP=NEW
NALL,100.
*END STEP
*******************************************