*HEADING
: DIRECT ENFORCEMENT WITH LOAD CONTROL FOR ELASTIC-PLASTIC MATERIAL,
          HIGH INITIAL STIFFNESS [GK3D18]
*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=GK3D18,ELSET=EINT
1,1,2,3,4,5,6,7,8,10,11,12,13,14,15,16,17,9,18
*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.,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.0, 1.0,-30., 2.0,-35., 3.0,-0., 
4.0,-35., 5.0,-10., 6.0,-0., 7.0,-37., 
8.0,-40., 9.0,-0.0, 10.0,-40., 11.0,-42., 
12.0,-0.0, 13.0,-42., 14.0,-44., 15.0,-0.0, 
16.0,-44., 17.0,-46.43, 18.0,-0., 19.0,-46.43,
20.0,-82.14, 21.0,-0.0, 22.0,-82.14, 23.0,-85., 
24.0,-0., 25.0,-85., 26.0,-100., 27.0,-0.0, 
28.0,-130., 29.0,-170., 30.0,-0.
*NSET,NSET=LNODES,GENERATE
2,9,1
*NSET,NSET=TNODES,GENERATE
10,18,1
*NSET,NSET=TCRNDS,GENERATE
10,13,1
*NSET,NSET=TMDNDS,GENERATE
14,17,1
*NSET,NSET=LEFT
1,5,2,10,14,11
*NSET,NSET=RIGHT
4,7,3,12,16,13
*NSET,NSET=RCRN
4,13,3,12
*NSET,NSET=RMID
7,16
*NSET,NSET=FRCRN
2,11,3,12
*NSET,NSET=FRMID
6,15
*NSET,NSET=MIDDLE
8,9,6,17,18,15
*BOUNDARY
1,1,3
LNODES,2,2
2,1,1
*******************************************
*STEP,INC=10000
Test Thickness Behavior of Gasket
*STATIC
0.01,30.,,0.01
*CONTROLS,PARAMETER=TIME INCREMENT

,0.25
*PRINT, FREQ=100
*CLOAD,AMP=HIST
TCRNDS,2,0.027778
TMDNDS,2,0.11111
18,2,0.44444
*NODE PRINT,FREQ=0
U,RF
*NODE FILE,FREQ=148
U,RF
*EL PRINT,FREQ=0
S, 
E, 
NE, 
EE, 
PE, 
EVOL, 
THE, 
TEMP, 
*EL FILE,FREQ=148
E, 
EE, 
NE, 
PE, 
S, 
EVOL, 
THE, 
*OUTPUT,FIELD,FREQ=148,VAR=PRESELECT
*END STEP
*******************************************
*STEP
Test of Membrane Behavior of Element
*STATIC
*CONTROLS,RESET
*BOUNDARY,OP=NEW
LEFT,1,2,0.0
1,3,3,0.0
RIGHT,2,2,0.0
MIDDLE,2,2,0.0
*CLOAD,OP=NEW
RCRN,1,0.25
RMID,1,1.0
*END STEP
*******************************************
*STEP
Test Transverse Shear Stiffness of Element
*STATIC
*BOUNDARY,OP=NEW
LNODES,1,3,0.0
1,1,3,0.0
*CLOAD,OP=NEW
TCRNDS,1,0.125
TMDNDS,1,0.5
18,1,2.0
TCRNDS,3,0.125
TMDNDS,3,0.5
18,3,2.0
*END STEP
*******************************************
*STEP
Free Thermal Expansion of Element
*STATIC
*BOUNDARY,OP=NEW
1,1,3
2,1
LNODES,2
*CLOAD,OP=NEW
*TEMPERATURE,OP=NEW
NALL,100.
*END STEP
*******************************************