*HEADING 
: DIRECT ENFORCEMENT WITH LOAD CONTROL FOR ELASTIC MATERIAL 
          WITH DAMAGE, INITIAL GAP [GKAX4]
*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=GKAX4,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
200.,
*GASKET THICKNESS  BEHAVIOR,DIRECTION=LOADING,TYPE=DAMAGE
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,DIRECTION=UNLOADING,TYPE=DAMAGE
0.0, 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.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
2,1,2
*******************************************
*STEP,INC=10000,NLGEOM
Test Normal Behavior of Gasket
*STATIC
.01,20.,,0.01
*CLOAD,AMP=HIST
3,2,10.47236
4,2,11.51936
*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, 
*OUTPUT,FIELD,FREQ=148,VAR=PRESELECT
*END STEP
*******************************************
*STEP
Test Membrane Behavior of Element
*STATIC
*BOUNDARY,OP=NEW
1,1,2,0.0
2,2,2,0.0
3,1,2,0.0
4,2,2,0.0
*CLOAD,OP=NEW
2,1,-50.0
4,1,-50.0
*END STEP
*******************************************
*STEP
Test Transverse Shear Stiffness of Element
*STATIC
*BOUNDARY,OP=NEW
1,1,2,0.0
3,2,2,0.0
2,1,2,0.0
4,2,2,0.0
*CLOAD,OP=NEW
3,1,-4.0608
4,1,-7.2936
*END STEP
*******************************************
*STEP
Thermal Expansion of Element
*STATIC
*BOUNDARY,OP=NEW
1,2
2,2
*CLOAD,OP=NEW
*TEMPERATURE,OP=NEW
NALL,100.
*EL FILE,FREQ=148
S, 
EVOL, 
THE, 
*NODE FILE
RF
*END STEP
*******************************************