*HEADING
: DIRECT ENFORCEMENT WITH DISPLACEMENT CONTROL FOR ELASTIC-PLASTIC MATERIAL, 
          LOW INITIAL STIFFNESS [GK2D2], DEPENDENCY FOR CONTACT AREA
*NODE,NSET=NALL
1,0.,0.,1.
2,0.,1.,1.
*ELEMENT,TYPE=GK2D2,ELSET=EINT
1,1,2
*GASKET SECTION,ELSET=EINT,BEHAVIOR=IS
,,,,0.,1.,0.
*GASKET BEHAVIOR,NAME=IS
*GASKET ELASTICITY,COMPONENT=TRANSVERSE SHEAR,VARIABLE=FORCE,
DEPENDENCIES=1
200.,0.0,0.0
200.,0.0,100.0
*GASKET THICKNESS BEHAVIOR,DIRECTION=LOADING,VARIABLE=FORCE,
DEPENDENCIES=1
0.,0., 0.0, 0.0
10.71,0.052,0.0, 0.0
28.57,0.087,0.0, 0.0
46.43,0.139,0.0, 0.0
82.14,0.217,0.0, 0.0
164.3,0.287,0.0, 0.0
0.,0., 0.0, 100.0
10.71,0.052,0.0, 100.0
28.57,0.087,0.0, 100.0
46.43,0.139,0.0, 100.0
82.14,0.217,0.0, 100.0
164.3,0.287,0.0, 100.0
*GASKET THICKNESS BEHAVIOR,DIRECTION=UNLOADING,VARIABLE=FORCE,
DEPENDENCIES=1
0.,0.087,0.087,0.0,0.0
14.28,0.122,0.087,0.0,0.0
46.43,0.139,0.087,0.0,0.0
0.,0.169,0.169,0.0,0.0
25.,0.20,0.169,0.0,0.0
82.14,0.217,0.169,0.0,0.0
0.,0.23,0.23,0.0,0.0
35.7,0.265,0.23,0.0,0.0
164.3,0.287,0.23,0.0,0.0
0.,0.087,0.087,0.0,100.0
14.28,0.122,0.087,0.0,100.0
46.43,0.139,0.087,0.0,100.0
0.,0.169,0.169,0.0,100.0
25.,0.20,0.169,0.0,100.0
82.14,0.217,0.169,0.0,100.0
0.,0.23,0.23,0.0,100.0
35.7,0.265,0.23,0.0,100.0
164.3,0.287,0.23,0.0,100.0
*GASKET CONTACT AREA,DEPENDENCIES=1
1.0,0.0,0.0,0.0
1.5,0.3,0.0,0.0
3.0,0.8,0.0,0.0
1.0,0.0,0.0,100.0
1.5,0.3,0.0,100.0
3.0,0.8,0.0,100.0
*EXPANSION
1.e-2,
*******************************************
*AMPLITUDE,NAME=HIST
 0.,-.000,    1.,-.060,   2.,-.000,   3.,-.087,
 4.,-.100,    5.,-.000,   6.,-.138,   7.,-.140,
 8.,-.080,    9.,-.000,  10.,-.168,  11.,-.170,
12.,-.100,   13.,-.000,  14.,-.199,  15.,-.201,
16.,-.100,   17.,-.000,  18.,-.216,  19.,-.218,
20.,-.100,   21.,-.000,  22.,-.264,  23.,-.266,
24.,-.100,   25.,-.000,  26.,-.286,  27.,-.300   
*BOUNDARY
1,1,2
*******************************************
*INITIAL CONDITIONS, TYPE=FIELD, VARIABLE=1
NALL,50.0
*STEP,INC=10000
Test Thickness Behavior of Gasket
*STATIC
0.01,27.,,0.01
*PRINT, FREQUENCY=100
*OUTPUT, FIELD, FREQUENCY=100
*OUTPUT, HISTORY, FREQUENCY=100
*BOUNDARY,AMP=HIST
2,2,,1.0
*NODE PRINT,FREQ=0
U,RF
*NODE FILE,FREQ=148
U, 
RF, 
*EL PRINT,FREQ=0
S, 
E, 
NE, 
EE, 
PE, 
EVOL, 
THE, 
TEMP, 
CS11, 
*EL FILE,FREQ=148
E, 
EE, 
NE, 
PE, 
S, 
CS11, 
EVOL, 
THE, 
*END STEP
*******************************************
**No Test of Membrane Behavior of Element
*******************************************
*STEP
Test Transverse Shear Stiffness of Element
*STATIC
*BOUNDARY,OP=NEW
1,1,2,0.0
2,2,2,0.0
2,1,1,0.025
*END STEP
*******************************************
*STEP
Free Thermal Expansion of Element
*STATIC
*BOUNDARY,OP=NEW
1,1,2,0.0
*TEMPERATURE,OP=NEW
NALL,100.
*END STEP
*******************************************