*HEADING
: CONTACT PAIR DISPLACEMENT CONTROL FOR ELASTIC-PLASTIC MATERIAL,
          HIGH INITIAL STIFFNESS [GK3D6N]
*NODE,NSET=NALL
1,0.0,0.0,0.0
2,1.0,0.0,1.0
3,2.0,0.0,1.0
4,0.0,1.0,0.0
5,1.0,1.0,1.0
6,2.0,1.0,0.0
7,0.0,1.0,0.0
8,1.0,1.0,1.0
9,2.0,1.0,0.0
10,1.0,2.0,0.5
*ELEMENT, TYPE=R3D3, ELSET=RIGID
2,7,8,9
*RIGID BODY, ELSET=RIGID, REF NODE=10
1.0, 
*ELEMENT,TYPE=GK3D6N,ELSET=EINT
1,1,2,3,4,5,6
*GASKET SECTION,ELSET=EINT,BEHAVIOR=IS
*GASKET BEHAVIOR,NAME=IS
*GASKET ELASTICITY,COMPONENT=MEMBRANE
2000.,0.2
*GASKET ELASTICITY,COMPONENT=TRANSVERSE SHEAR
200.,
**
** loading type : EPH - Elastic-Plastic with High Initial Stiffness
**
*GASKET THICKNESS BEHAVIOR,DIRECTION=LOADING
0.,0.
20.,0.005
45.,0.05
70.,0.135
82.,0.22
83.,0.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,-0., 1.0,-0.089, 2.0,-0.098, 3.0,-0.0, 
4.0,-0.098, 5.0,-0.005, 6.0,-0.0, 7.0,-0.105, 
8.0,-0.11, 9.0,-0.0, 10.0,-0.11, 11.0,-0.12, 
12.0,-0.0, 13.0,-0.12, 14.0,-0.13, 15.0,-0.0, 
16.0,-0.13, 17.0,-0.139, 18.0,-0.0, 19.0,-0.139,
20.0,-0.217, 21.0,-0.0, 22.0,-0.217, 23.0,-0.240, 
24.0,-0.0, 25.0,-0.240, 26.0,-0.265, 27.0,-0.0, 
28.0,-0.27, 29.0,-0.295, 30.0,-0.0
*BOUNDARY
1,1
2,1
3,1
10,1,1
10,3,6
*******************************************
*surface, NAME=SURFA
2,SNEG
*surface, NAME=SURFB
1,SPOS
*CONTACT PAIR, INTERACTION=SMOOTH,TIED,ADJUST=0.01
SURFB,SURFA
*SURFACE INTERACTION, NAME=SMOOTH
*STEP,INC=10000
Test Normal Behavior of Gasket
*STATIC
0.01,30.,,0.01
*BOUNDARY,AMP=HIST
10,2,2,1.0
*NODE PRINT,FREQ=0
U, 
RF, 
*NODE FILE,FREQ=148
U, 
RF, 
*EL PRINT,FREQ=0, ELSET=EINT
S, 
E, 
NE, 
EE, 
PE, 
EVOL, 
THE, 
TEMP, 
*EL FILE,FREQ=148, ELSET=EINT
E, 
EE, 
NE, 
PE, 
S, 
EVOL, 
THE, 
*OUTPUT,FIELD,FREQ=148,VAR=PRE
*END STEP
***************************************
*STEP
Free Thermal Expansion of Element
*STATIC
*BOUNDARY,OP=NEW
1,1
2,1
3,1
*TEMPERATURE,OP=NEW
NALL,100.
*END STEP
***************************************