*HEADING
: DIRECT ENFORCEMENT WITH LOAD CONTROL FOR ELASTIC MATERIAL
          WITH DAMAGE, INITIAL GAP [GKPS4]
*NODE,NSET=NALL
*ELEMENT,TYPE=GKPS4,ELSET=EINT
*GASKET SECTION,ELSET=EINT,BEHAVIOR=IS
*GASKET BEHAVIOR,NAME=IS
*GASKET ELASTICITY,COMPONENT=MEMBRANE
*GASKET ELASTICITY,COMPONENT=TRANSVERSE SHEAR
*GASKET THICKNESS BEHAVIOR,TYPE=DAMAGE,DIRECTION=LOADING
*GASKET THICKNESS BEHAVIOR,TYPE=DAMAGE,DIRECTION=UNLOADING
*EXPANSION
1.e-2,
*******************************************
*AMPLITUDE,NAME=HIST
 0.,  0.000,  1., -5.000,  2.,  -9.000,  3., -0.100,
 4., -9.000,  5.,-12.000,  6., -0.100,  7.,-12.000,
 8.,-45.000,  9.,-48.000, 10.,-30.000, 11., -0.100,
12.,-48.000, 13.,-82.500, 14., -0.100, 15.,-82.500,
16.,-160.000,17., -0.001, 18.,-160.0, 19.,-200.
*BOUNDARY
1,1,2
2,1,2
*******************************************
*STEP,INC=10000
Test Thickness Behavior of Gasket
*STATIC
0.01,19.,,0.01
*CLOAD,AMP=HIST
3,2,.5
4,2,.5
*NODE PRINT,FREQ=0
U,RF
*EL PRINT,FREQ=0
S,
E,
NE,
EE,
EVOL,
THE,
TEMP,
*OUTPUT, FIELD, FREQ=148
*NODE OUTPUT
U,
RF,
*ELEMENT OUTPUT
E,
EE,
NE,
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,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
Free Thermal Expansion of Element
*STATIC
*CLOAD,OP=NEW
*BOUNDARY,OP=NEW
1,1,2
2,2
*TEMPERATURE,OP=NEW
NALL,100.
*END STEP
*******************************************