*HEADING
: CONTACT PAIR WITH LOAD CONTROL FOR ELASTIC-PLASTIC MATERIAL, 
          HIGH INITIAL STIFFNESS, NOT ALIGNED W/ GLOBAL AXES [GKPS4N]
*PREPRINT,MODEL=YES
*SYSTEM
0.0,0.0,0.0,1.0,-1.0,0.0
1.0,1.0,0.0
*NODE,NSET=NALL
1,
2,1.
3,0.,1.
4,1.,1.
1001,-2.,0.
1002, 2.,0.
1003,-2.,1.
1004, 2.,1.
*TRANSFORM,NSET=REFNODES
1.0,-1.0,0.0,1.0,1.0,0.0
*NODE,NSET=REFNODES
9998,0.,0.
9999,0.,1.
******************************************************
*ELEMENT,TYPE=R2D2,ELSET=RIGID1
1001,1001,1002
*ELEMENT,TYPE=R2D2,ELSET=RIGID2
1002,1003,1004
*RIGID BODY,ELSET=RIGID1,REFNODE=9998
*RIGID BODY,ELSET=RIGID2,REFNODE=9999
*surface,NAME=MLOWER
RIGID1,SPOS
*surface,NAME=MUPPER
RIGID2,SNEG
*surface,NAME=TOP
EINT,SPOS
*surface,NAME=BOTTOM
EINT,SNEG
*CONTACT PAIR,INTERACTION=ROUGH,ADJUST=0.001,TIED
BOTTOM,MLOWER
TOP,MUPPER
*SURFACE INTERACTION,NAME=ROUGH
*******************************************************
*ELEMENT,TYPE=GKPS4N,ELSET=EINT
1,1,2,3,4
*GASKET SECTION,ELSET=EINT,BEHAVIOR=IS
*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
9998,1,2
9998,6,6
9999,1,1
9999,6,6
*******************************************
*STEP,INC=10000
Test Thickness Behavior of Gasket
*STATIC
0.01,16.,,0.01
*CONTROLS,PARAMETER=TIME INCREMENT

,0.25
*PRINT, FREQUENCY=100
*OUTPUT, HISTORY, FREQUENCY=100
*CLOAD,AMP=HIST
9999,2,1.0
*NODE PRINT,FREQ=148,GLOBAL=YES
U,RF
*NODE PRINT,FREQ=148,GLOBAL=NO
U,RF
*NODE FILE,FREQ=148,GLOBAL=YES
U, 
RF, 
*NODE FILE,FREQ=148,GLOBAL=NO
U, 
RF, 
*EL PRINT,ELSET=EINT,FREQ=0
S, 
E, 
NE, 
EE, 
PE, 
EVOL, 
THE, 
TEMP, 
*EL FILE,ELSET=EINT,FREQ=148
E, 
EE, 
NE, 
PE, 
S, 
EVOL, 
THE, 
*OUTPUT,FIELD,FREQ=148,VAR=PRE
*OUTPUT, HISTORY, FREQ=148
*NODE OUTPUT, NSET=REFNODES, GLOBAL=NO
 U,RF
*END STEP
*******************************************
** No Test Membrane Behavior of Element
** No Test Transverse Shear Stiffness Behavior of Element
** No Test Free Thermal Expansion of Element
*******************************************