*HEADING
: DIRECT ENFORCEMENT WITH LOAD CONTROL FOR ELASTIC-PLASTIC MATERIAL, 
          HIGH INITIAL STIFFNESS, INITIAL GAP [GKPE4]
*NODE,NSET=NALL
1,
2,1.
3,0.,1.
4,1.,1.
*ELEMENT,TYPE=GKPE4,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
0.,0.
20.,0.005
45.,0.05
70.,0.135
82,.22
83,.287
*GASKET THICKNESS BEHAVIOR,DIRECTION=UNLOADING
0.,0.08,0.08
14.,0.12,0.08
70.,0.135,0.08
0.,0.17,0.17
25.,0.20,0.17
82.,0.22,0.17
0.,0.23,0.23
36.,0.265,0.23
83.,0.287,0.23
*EXPANSION
1.e-2,
*******************************************
*AMPLITUDE,NAME=HIST
 0.,  0.000,  1.,-15.000,  2.,-25.000,  3., -0.100,
 4.,-25.000,  5.,-40.000,  6., -0.100,  7.,-45.000,
 8.,-45.000,  9.,-60.000, 10.,-75.000, 11., -0.100, 
12.,-75.000, 13.,-82.500, 14., -0.100, 15.,-82.500,
16.,-95.000 
*BOUNDARY
1,1,2
2,1,2
*******************************************
*STEP,INC=10000
Test Thickness Behavior of Gasket
*STATIC
0.01,16.,,0.01
*CONTROLS,PARAMETER=TIME INCREMENT

,0.25
*CLOAD,AMP=HIST
3,2,.5
4,2,.5
*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,1,2,0.0
4,2,2,0.0
*CLOAD,OP=NEW
2,1,0.5
4,1,0.5
*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,.5
4,1,.5
*END STEP
*******************************************
*STEP
Thermal Expansion of Element
*STATIC
*CLOAD,OP=NEW
*BOUNDARY,OP=NEW
1,1,2
2,2
*TEMPERATURE,OP=NEW
NALL,100.
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
*NODE OUTPUT
RF,
*END STEP
*******************************************