*HEADING 
: DIRECT ENFORCEMENT LOAD CONTROL FOR ELASTIC-PLASTIC MATERIAL 
          HIGH INITIAL STIFFNESS, INITIAL GAP [GK3D6L]
*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=GK3D6L,ELSET=EINT
1,1,2,3,4,5,6 
*GASKET SECTION,ELSET=EINT,BEHAVIOR=IS
,0.1
*GASKET BEHAVIOR,NAME=IS
*GASKET ELASTICITY,COMPONENT=MEMBRANE
2000.,0.2
*GASKET ELASTICITY,COMPONENT=TRANSVERSE SHEAR,VAR=FORCE
200.,
*GASKET THICKNESS  BEHAVIOR,DIRECTION=LOADING,TYPE=ElASTIC-PLASTIC,VAR=FORCE
0.,0.
20.0, 0.005
45.0, 0.05
70.0, 0.135
82.0, 0.22
83.0, 0.287
*GASKET THICKNESS  BEHAVIOR,DIRECTION=UNLOADING,TYPE=ElASTIC-PLASTIC,VAR=FORCE
0.,   0.08,  0.08
14.,  0.12,  0.08
70.,  0.135, 0.08
0.0,  0.17,  0.17
25.0, 0.20,  0.17
82.0, 0.22,  0.17
0.0,  0.23,  0.23
36.0, 0.265, 0.23
83.0, 0.287, 0.23
*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,3
2,1,3
3,1,3
*******************************************
*STEP,INC=10000,NLGEOM
Test Normal Behavior of Gasket
*STATIC
0.01,20.,,0.01
*CONTROLS,PARAMETER=TIME INCREMENT

,0.25
*CLOAD,AMP=HIST
4, 2,0.1667
5, 2,0.6667
6, 2,0.1667
*NODE PRINT,FREQ=0
U,RF
*NODE FILE,FREQ=148
U, 
RF, 
*EL PRINT,FREQ=0
S, 
E, 
NE, 
EE, 
PE, 
EVOL, 
THE, 
TEMP, 
*EL FILE,FREQ=148
E, 
EE, 
NE, 
PE, 
S, 
EVOL, 
THE, 
*OUTPUT,FIELD,VAR=PRESELECT,FREQ=148
*END STEP
*******************************************
*STEP
Test Membrane Behavior of Element
*STATIC
*CONTROLS,RESET
*BOUNDARY,OP=NEW
1, 1, 2, 0.0
2, 2, 2, 0.0
3, 2, 2, 0.0
4, 1, 2, 0.0
5, 2, 2, 0.0
6, 2, 2, 0.0
*CLOAD,OP=NEW
3, 1, 50.0
6, 1, 50.0
*END STEP
*******************************************
*STEP
Test Transverse Shear Stiffness of Element
*STATIC
*BOUNDARY,OP=NEW
1,  1,  2,  0.0
2,  1,  2,  0.0
3,  1,  2,  0.0
4,  2,  2,  0.0
5,  2,  2,  0.0
6,  2,  2,  0.0
*CLOAD,OP=NEW
4,  1,  0.8335
5,  1,  3.3335
6,  1,  0.8335
*END STEP
*******************************************
*STEP
Free Thermal Expansion of Element
*STATIC
*BOUNDARY,OP=NEW
1,1,2
2,2
3,2
*CLOAD,OP=NEW
*TEMPERATURE,OP=NEW
NALL,100.
*OUTPUT,FIELD
*node output
RF
*element output
S
*END STEP
*******************************************