*HEADING
: DIRECT ENFORCEMENT WITH DISPLACEMENT CONTROL FOR ELASTIC-PLASTIC MATERIAL, 
          LOW INITIAL STIFFNESS, STACKED, NOT ALIGNED W/ GLOBAL AXES [GKPS6]
*SYSTEM
0.0,0.0,0.0,1.0,-1.0,0.0
1.0,1.0,0.0
*TRANSFORM,NSET=NALL
1.0,-1.0,0.0,1.0,1.0,0.0
*NODE,NSET=NALL
1,0.,0.
2,.5,0.
3,1.,0.
1001,0.,0.5
1002,.5,0.5
1003,1.,0.5
4,0.,1.
5,.5,1.
6,1.,1.
*ELEMENT,TYPE=GKPS6,ELSET=EINT
1,1,2,3,1001,1002,1003
2,1001,1002,1003,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.,
*GASKET THICKNESS BEHAVIOR,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,DIRECTION=UNLOADING
0.,0.087,0.087
14.28,0.122,0.087
46.43,0.139,0.087
0.,0.169,0.169
25.,0.20,0.169
82.14,0.217,0.169
0.,0.23,0.23
35.7,0.265,0.23
164.3,0.287,0.23
*EXPANSION
1.e-2,
*******************************************
*AMPLITUDE,NAME=HIST
 0.,-.000,    1.,-.060,   2.,-.000,   3.,-.087,
 4.,-.100,    5.,-.000,   6.,-.138,   7.,-.140,
 8.,-.080,    9.,-.000,  10.,-.168,  11.,-.170,
12.,-.100,   13.,-.000,  14.,-.199,  15.,-.201,
16.,-.100,   17.,-.000,  18.,-.216,  19.,-.218,
20.,-.100,   21.,-.000,  22.,-.264,  23.,-.266,
24.,-.100,   25.,-.000,  26.,-.286,  27.,-.300   
*BOUNDARY
1,1,2
2,1,2
3,1,2
*******************************************
*STEP,INC=10000
Test Thickness Behavior of Gasket
*STATIC
0.01,27.,,0.01
*BOUNDARY,AMP=HIST
4,2,,1.0
5,2,,1.0
6,2,,1.0
*NODE PRINT,FREQ=0
U,RF
*EL PRINT,FREQ=0
S, 
E, 
NE, 
EE, 
PE, 
EVOL, 
THE, 
TEMP, 
*OUTPUT, FIELD, FREQ=148
*NODE OUTPUT
U, 
RF, 
*ELEMENT OUTPUT
E, 
EE, 
NE, 
PE, 
S, 
EVOL, 
THE, 
*END STEP
*******************************************
*STEP
Test Membrane Behavior of Element
*STATIC
*BOUNDARY,OP=NEW
1,1,2,0.0
1001,1,2,0.0
4,1,2,0.0
2,2,,0.0
5,2,,0.0
3,2,,0.0
6,2,,0.0
3,1,,.025
1003,1,,.025
6,1,,.025
*END STEP
*******************************************
*STEP
Test Transverse Shear Stiffness of Element
*STATIC
*BOUNDARY,OP=NEW
1,1,2,0.0
2,1,2,0.0
3,1,2,0.0
4,2,2,0.0
5,2,2,0.0
6,2,2,0.0
4,1,1,0.025
5,1,1,0.025
6,1,1,0.025
*END STEP
*******************************************
*STEP
Free Thermal Expansion of Element
*STATIC
*BOUNDARY,OP=NEW
1,1,2,0.0
2,2,,0.0
3,2,,0.0
4,1,,0.0
*TEMPERATURE,OP=NEW
NALL,100.
*END STEP
*******************************************