*HEADING
: CONTACT PAIR WITH DISPLACEMENT CONTROL FOR ELASTIC MATERIAL 
          WITH DAMAGE, THROUGH THICKNESS ONLY BEHAVIOR, INITIAL GAP, 
          NOT ALIGNED W/ GLOBAL AXES [GK2D2N]
*SYSTEM
0.0,0.0,0.0,1.0,-1.0,0.0
1.0,1.0,0.0
*NODE,NSET=NALL
1,0.,0.
2,0.,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.
******************************************************
** Add weak 1-node springs to remove R.B.M. at REFNODE 9999
** The R.B.M. occurs in Step 2 and can be a problem on some
** machines if not constrained properly. DOF 6 at node 9999
** is also constrained to prevent rotational R.B.M.
**
*element,type=spring1,elset=ground1
9991,9999
*element,type=spring1,elset=ground2
9992,9999
*spring,elset=ground1
1, 
1.e-6, 
*spring,elset=ground2
2, 
1.e-6, 
*elset,elset=eout
1,1001,1002
******************************************************
*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,type=node,NAME=TOP
2,1.0
*surface,type=node,NAME=BOTTOM
1,1.0
*CONTACT PAIR,INTERACTION=ROUGH,ADJUST=0.001,TIED
BOTTOM,MLOWER
TOP,MUPPER
*SURFACE INTERACTION,NAME=ROUGH
*******************************************************
*ELEMENT,TYPE=GK2D2N,ELSET=EINT
1,1,2
*GASKET SECTION,ELSET=EINT,BEHAVIOR=IS
,.1
*GASKET BEHAVIOR,NAME=IS
*GASKET THICKNESS BEHAVIOR,TYPE=DAMAGE,DIRECTION=LOADING
0.,0.
10.71,0.052
28.57,0.087
46.43,0.139
82.14,0.217
164.3,0.287
*GASKET THICKNESS BEHAVIOR,TYPE=DAMAGE,DIRECTION=UNLOADING
0.,0.,0.139
14.28,0.122,0.139
46.43,0.139,0.139
0.,0.,0.217
12.,0.17,0.217
40.,0.21,0.217
82.14,0.217,0.217
0.,0.0,0.287
11.0,0.24,0.287
164.3,0.287,0.287
*EXPANSION
1.e-2,
*******************************************
*AMPLITUDE,NAME=HIST
 0.,-.100,    1.,-.160,   2.,-.100,   3.,-.187,
 4.,-.200,    5.,-.100,   6.,-.238,   7.,-.240,
 8.,-.180,    9.,-.100,  10.,-.268,  11.,-.270,
12.,-.200,   13.,-.100,  14.,-.299,  15.,-.301,
16.,-.200,   17.,-.100,  18.,-.316,  19.,-.318,
20.,-.200,   21.,-.100,  22.,-.364,  23.,-.366,
24.,-.200,   25.,-.100,  26.,-.386,  27.,-.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,27.,,0.01
*PRINT, FREQUENCY=100
*OUTPUT, HISTORY, FREQUENCY=100
*BOUNDARY,AMP=HIST
9999,2,,1.0
*NODE PRINT,FREQ=0
U,RF
*NODE FILE,FREQ=148
U, 
RF, 
*OUTPUT,FIELD,FREQ=148
*NODE OUTPUT
U,
RF,
*EL PRINT,FREQ=0
S, 
E, 
NE, 
EE, 
EVOL, 
THE, 
TEMP, 
*EL FILE,FREQ=148,ELSET=EOUT
E, 
EE, 
NE, 
S, 
EVOL, 
THE, 
*OUTPUT,FIELD,FREQ=148
*ELEMENT OUTPUT,ELSET=EOUT
E,
EE,
NE,
S,
EVOL,
THE,
*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
*******************************************