*HEADING
: DIRECT ENFORCEMENT WITH DISPLACEMENT CONTROL FOR ELASTIC MATERIAL 
          WITH DAMAGE [GK2D2]
*NODE,NSET=NALL
1,0.,0.,1.
2,0.,1.,1.
*ELEMENT,TYPE=GK2D2,ELSET=ELEMENTS-AT-INTERFACE
1,1,2
*ELSET,ELSET=ALL_ELEMENTS
ELEMENTS-AT-INTERFACE,
*GASKET SECTION,ELSET=ELEMENTS-AT-INTERFACE,BEHAVIOR=INTERFACE-BEHAVIOR
,,,,0.,1.,0.
*GASKET BEHAVIOR,NAME=INTERFACE-BEHAVIOR
*GASKET ELASTICITY,COMPONENT=TRANSVERSE SHEAR,VARIABLE=FORCE
200.,
*GASKET THICKNESS BEHAVIOR,TYPE=DAMAGE,DIRECTION=LOADING,VARIABLE=FORCE
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,VARIABLE=FORCE
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,0.287
11.0,0.24,0.287
164.3,0.287,0.287
*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
*******************************************
*STEP,INC=10000
Test Thickness Behavior of Gasket
*STATIC
0.01,27.,,0.01
*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, 
EVOL, 
THE, 
TEMP, 
*EL FILE,FREQ=148
E, 
EE, 
NE, 
S, 
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
*******************************************