*HEADING
: DIRECT ENFORCEMENT WITH LOAD CONTROL FOR ELASTIC-PLASTIC MATERIAL,
          HIGH INITIAL STIFFNESS, INITIAL GAP [GKPE4]
*NODE,NSET=NALL
*ELEMENT,TYPE=GKPE4,ELSET=EINT
*GASKET SECTION,ELSET=EINT,BEHAVIOR=IS
*GASKET BEHAVIOR,NAME=IS
*GASKET ELASTICITY,COMPONENT=MEMBRANE
*GASKET ELASTICITY,COMPONENT=TRANSVERSE SHEAR
*GASKET THICKNESS BEHAVIOR,DIRECTION=LOADING
*GASKET THICKNESS BEHAVIOR,DIRECTION=UNLOADING
*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.
*END STEP
*******************************************