*HEADING 
: DIRECT ENFORCEMENT WITH LOAD CONTROL FOR ELASTIC-PLASTIC MATERIAL, 
          LOW INITIAL STIFFNESS, INITIAL GAP [GKAX4]
** UVARM FOR GASKET 
*RESTART,WRITE,OVERLAY
*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=ELASTIC-PLASTIC,
YIELD=0.08700,TENSILE STIFFNESS FACTOR=1.0E-03
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=ELASTIC-PLASTIC
0.,0.087,0.087
14.28,0.122,0.087
46.43,0.139,0.087 
0.,0.169,0.169
25.,0.20,0.169
82.14,0.217,0.169
0.,0.23,0.23
35.7,0.265,0.23
164.3,0.287,0.23
*EXPANSION
1.e-2,
*USER OUTPUT VARIABLES
1, 
*******************************************
*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=148
S, 
E, 
NE, 
EE, 
PE, 
EVOL, 
THE, 
TEMP, 
UVARM1,
*EL FILE,FREQ=148
E, 
EE, 
NE, 
PE, 
S, 
EVOL, 
THE, 
*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.
**END STEP
*******************************************