*HEADING
:DIRECT ENFORCEMENT WITH DISP. CONTROL FOR ELASTIC MATERIAL WITH DAMAGE [GK3D12M]
** SOLID ELEMENT NUMBERING OPTION ON GASKET ELEMENTS 
*NODE,NSET=NALL
1,0.,0.,0.
2,0.,0.,1.
3,1.,0.,1.
4,1.,0.,0.
5,0.,0.,0.5
6,0.5,0.,1.
7,1.,0.,0.5
8,0.5,0.,0.
9,0.5,0.,0.5
10,0.,1.,0.
11,0.,1.,1.
12,1.,1.,1.
13,1.,1.,0.
14,0.,1.,0.5
15,0.5,1.,1.
16,1.,1.,0.5
17,0.5,1.,0.
18,0.5,1.,0.5
*ELEMENT,TYPE=GK3D12M,ELSET=EINT,SOLID ELEMENT NUMBERING
1,1,2,4,10,11,13,5,9,8,14,18,17
2,2,3,4,11,12,13,6,7,9,15,16,18
**1,1,2,4,5,9,8,10,11,13,14,18,17
**2,2,3,4,6,7,9,11,12,13,15,16,18
*GASKET SECTION,ELSET=EINT,BEHAVIOR=IS
,0.1
*GASKET BEHAVIOR,NAME=IS
*GASKET ELASTICITY,COMPONENT=MEMBRANE
2000.,0.2
*GASKET ELASTICITY,COMPONENT=TRANSVERSE SHEAR
200.,
*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.287
11.0,0.24,0.287
164.3,0.287,0.287
*EXPANSION
1.e-2,
*******************************************
*AMPLITUDE,NAME=HIST
0.0,-0., 1.0,-0.189, 2.0,-0.198, 3.0,-0.1, 
4.0,-0.198, 5.0,-0.105, 6.0,-0.1, 7.0,-0.205, 
8.0,-0.21, 9.0,-0.1, 10.0,-0.21, 11.0,-0.22, 
12.0,-0.1, 13.0,-0.22, 14.0,-0.23, 15.0,-0.1, 
16.0,-0.23, 17.0,-0.239, 18.0,-0.1, 19.0,-0.239,
20.0,-0.317, 21.0,-0.1, 22.0,-0.317, 23.0,-0.340, 
24.0,-0.1, 25.0,-0.340, 26.0,-0.365, 27.0,-0.1, 
28.0,-0.37, 29.0,-0.395, 30.0,-0.1
*NSET,NSET=LNODES,GENERATE
2,9,1
*NSET,NSET=TNODES,GENERATE
10,18,1
*NSET,NSET=LEFT
1,5,2,10,14,11
*NSET,NSET=RIGHT
4,7,3,13,16,12
*NSET,NSET=MIDDLE
8,9,6,17,18,15
*BOUNDARY
LNODES,2,2
1,1,3
2,1,1
*******************************************
*STEP,INC=10000
Test Thickness Behavior of Gasket
*STATIC
0.01,30.,,0.01
*BOUNDARY,AMP=HIST
TNODES,2,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, 
*OUTPUT,FIELD,FREQ=148,VAR=PRE
*EL FILE,FREQ=148
E, 
EE, 
NE, 
S, 
EVOL, 
THE, 
*END STEP
*******************************************
*STEP
Test of Membrane Behavior of Element
*STATIC
*BOUNDARY,OP=NEW
LEFT,1,2,0.0
1,3,3,0.0
RIGHT,1,1,0.025
RIGHT,2,2,0.0
MIDDLE,2,2,0.0
*END STEP
*******************************************
*STEP
Test Transverse Shear Stiffness of Element
*STATIC
*BOUNDARY,OP=NEW
LNODES,1,3,0.0
1,1,3,0.0
TNODES,1,1,0.025
TNODES,3,3,0.025
*END STEP
*******************************************
*STEP
Free Thermal Expansion of Element
*STATIC
*BOUNDARY,OP=NEW
1,1,3
2,1
LNODES,2
*TEMPERATURE,OP=NEW
NALL,100.
*END STEP
*******************************************