*HEADING
: DIRECT ENFORCEMENT LOAD CONTROL FOR ELASTIC MATERIAL WITH DAMAGE, 
          INITIAL GAP, INITIAL VOID W/ ORTHO THERMAL EXP. COEFF.,
          ORIENTATION [GK3D8]
*NODE,NSET=NALL
1,0.0,0.0,1.0
2,1.0,0.0,1.0
3,1.0,0.0,0.0
4,0.0,0.0,0.0
5,0.0,1.0,1.0
6,1.0,1.0,1.0
7,1.0,1.0,0.0
8,0.0,1.0,0.0
*ORIENTATION,NAME=MATAXES
1.0,0.0,0.0,0.0,1.0,0.0
2,0.0
*ELEMENT,TYPE=GK3D8,ELSET=EINT
1,1,2,3,4,5,6,7,8
*GASKET SECTION,ELSET=EINT,BEHAVIOR=IS,ORIENTATION=MATAXES
,0.10,0.10,
*GASKET BEHAVIOR,NAME=IS
*GASKET ELASTICITY,COMPONENT=MEMBRANE
2000.,0.2
*GASKET ELASTICITY,COMPONENT=TRANSVERSE SHEAR
200.,
**
** loading type : DE - Damaged Elasticity
**
*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.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, TYPE=ORTHO
1.e-2,0.5e-2,0.5e-02
*******************************************
*AMPLITUDE,NAME=HIST
0.0,0.0, 1.0,-20., 2.0,-35., 3.0,-0.1,
4.0,-35., 5.0,-50., 6.0,-0.1, 7.0,-50.,
8.0,-65., 9.0,-0.1, 10.0,-65., 11.0,-80.,
12.0,-0.1, 13.0,-80., 14.0,-95., 15.0,-0.1,
16.0,-95., 17.0,-110., 18.0,-0.1, 19.0,-110.,
20.0,-125., 21.0,-0.1, 22.0,-125., 23.0,-140.,
24.0,-0.1, 25.0,-140., 26.0,-155., 27.0,-0.1,
28.0,-155., 29.0,-170., 30.0,-0.0
*BOUNDARY
1,1,3
2,1,3
3,1,3
4,1,3
*******************************************
*STEP,INC=10000
Test Normal Behavior of Gasket
*STATIC
0.1,30.,,0.1
*CONTROLS,PARAMETER=TIME INCREMENT

,0.25
*PRINT, FREQ=100
*CLOAD,AMP=HIST
5,2,0.25
6,2,0.25
7,2,0.25
8,2,0.25
*NODE PRINT,FREQ=0
U, 
RF, 
*NODE FILE,FREQ=148
U, 
RF, 
*EL PRINT,FREQ=0, ELSET=EINT
S, 
E, 
NE, 
EE, 
EVOL, 
THE, 
TEMP, 
*EL FILE,FREQ=148, ELSET=EINT
E, 
EE, 
NE, 
S, 
EVOL, 
THE, 
*OUTPUT,FIELD,FREQ=148,VAR=PRE
*END STEP
***************************************
*STEP
Test Membrane Behavior of Element
*STATIC
*CONTROLS,RESET
*BOUNDARY,OP=NEW
1,1,2
2,2,2
3,2,3
4,1,3
5,1,1
7,3,3
8,1,1
8,3,3
*CLOAD,OP=NEW
2,1,0.25
3,1,0.25
7,1,0.25
6,1,0.25
*END STEP
*******************************************
*STEP
Test Transverse Shear Stiffness of Element
*STATIC
*BOUNDARY,OP=NEW
1,1,3
2,1,3
3,1,3
4,1,3
*CLOAD,AMP=HIST,OP=NEW
5,1,0.25
5,3,0.25
6,1,0.25
6,3,0.25
7,1,0.25
7,3,0.25
8,1,0.25
8,3,0.25
*END STEP
*******************************************
*STEP
Free Thermal Expansion of Element
*STATIC
*BOUNDARY,OP=NEW
1,1,2
2,2,2
3,2,2
4,1,3
*TEMPERATURE,OP=NEW
NALL,100.
*CLOAD,OP=NEW
*END STEP
***************************************