*HEADING
: CONTACT PAIR WITH DISPLACEMENT CONTROL FOR ELASTIC-PLASTIC MATERIAL,
          HIGH INITIAL STIFFNESS, INITIAL GAP, INITIAL VOID W/ ISO THERMAL EXP. COEFF.,
          STACKED [GKPS6N]
*NODE,NSET=NALL
*NODE,NSET=REFNODES
******************************************************
*ELEMENT,TYPE=R2D2,ELSET=RIGID1
*ELEMENT,TYPE=R2D2,ELSET=RIGID2
*RIGID BODY,ELSET=RIGID1,REFNODE=9998
*RIGID BODY,ELSET=RIGID2,REFNODE=9999
*surface,NAME=MLOWER
RIGID1,SPOS
*surface,NAME=MUPPER
RIGID2,SNEG
*surface,type=node,NAME=TOP
*surface,type=node,NAME=BOTTOM
*CONTACT PAIR,INTERACTION=ROUGH,ADJUST=0.001,TIED
BOTTOM,MLOWER
TOP,MUPPER
*SURFACE INTERACTION,NAME=ROUGH
*******************************************************
*ELEMENT,TYPE=GKPS6N,ELSET=EINT
*GASKET SECTION,ELSET=EINT,BEHAVIOR=IS
*GASKET BEHAVIOR,NAME=IS
*GASKET THICKNESS BEHAVIOR,DIRECTION=LOADING
*GASKET THICKNESS BEHAVIOR,DIRECTION=UNLOADING
*EXPANSION
1.e-2,
*******************************************
*AMPLITUDE,NAME=HIST
 0.,-.104,    1.,-.100,   2.,-.106,   3.,-.100,
 4.,-.140,    5.,-.106,   6.,-.100,   7.,-.234,
 8.,-.236,    9.,-.100,  10.,-.319,  11.,-.321,
12.,-.386,   13.,-.100,  14.,-.387,  15.,-.400
*BOUNDARY
9998,1,2
9998,6,6
9999,1,1
9999,6,6
*******************************************
*STEP,INC=10000
Test Thickness Behavior of Gasket
*STATIC
0.01,15.,,0.01
*PRINT, FREQ=100
*BOUNDARY,AMP=HIST
9999,2,,1.0
*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,FREQ=148,VAR=PRE
*END STEP
*******************************************
**No Test of Membrane Behavior of Element
*******************************************
**No Test of Transverse Shear Stiffness of Element
*******************************************
*STEP
Free Thermal Expansion of Element
*STATIC
*BOUNDARY,OP=NEW
9998,1,2,0.0
9998,6,6,0.0
*TEMPERATURE,OP=NEW
NALL,100.
*END STEP
*******************************************