*HEADING 
: DIRECT ENFORCEMENT WITH LOAD CONTROL FOR ELASTIC-PLASTIC MATERIAL, 
          LOW INITIAL STIFFNESS, INITIAL GAP [GKAX4]
** DEPENDENCY ON FIELD VARIABLE AND TEMPERATURE FOR GASKET ELEMENTS
*NODE,NSET=NALL
1,3.0, 0.0
2,4.0, 0.0
3,3.0, 1.0
4,4.0, 1.0
*ELEMENT,TYPE=GKAX4,ELSET=EINT
1,1,2,3,4
*GASKET SECTION,ELSET=EINT,BEHAVIOR=IS
,0.1
*GASKET BEHAVIOR,NAME=IS
*GASKET ELASTICITY,COMPONENT=MEMBRANE,DEP=1
2000.,0.2,0.0,50.0
2000.,0.2,100.0,50.0
2000.,0.2,0.0,150.0
2000.,0.2,100.0,150.0
*GASKET ELASTICITY,COMPONENT=TRANSVERSE SHEAR,DEP=1
200.,0.0,50.0
200.,100.0,50.0
200.,0.0,150.0
200.,100.0,150.0
*GASKET THICKNESS BEHAVIOR,DIRECTION=LOADING,TYPE=ELASTIC-PLASTIC,
YIELD=0.08700,TENSILE STIFFNESS FACTOR=1.0E-03, DEP=1
0.,0.      ,0.0,50.0
10.71,0.052,0.0,50.0
28.57,0.087,0.0,50.0
46.43,0.139,0.0,50.0
82.14,0.217,0.0,50.0
164.3,0.287,0.0,50.0
0.,0.      ,100.0,50.0
10.71,0.052,100.0,50.0
28.57,0.087,100.0,50.0
46.43,0.139,100.0,50.0
82.14,0.217,100.0,50.0
164.3,0.287,100.0,50.0
0.,0.      ,0.0,150.0
10.71,0.052,0.0,150.0
28.57,0.087,0.0,150.0
46.43,0.139,0.0,150.0
82.14,0.217,0.0,150.0
164.3,0.287,0.0,150.0
0.,0.      ,100.0,150.0
10.71,0.052,100.0,150.0
28.57,0.087,100.0,150.0
46.43,0.139,100.0,150.0
82.14,0.217,100.0,150.0
164.3,0.287,100.0,150.0
*GASKET THICKNESS BEHAVIOR,DIRECTION=UNLOADING,TYPE=ELASTIC-PLASTIC,
DEP=1
0.,0.087,0.087   ,0.0,50.0
14.28,0.122,0.087,0.0,50.0
46.43,0.139,0.087,0.0,50.0 
0.,0.169,0.169   ,0.0,50.0
25.,0.20,0.169   ,0.0,50.0
82.14,0.217,0.169,0.0,50.0
0.,0.23,0.23     ,0.0,50.0
35.7,0.265,0.23  ,0.0,50.0
164.3,0.287,0.23 ,0.0,50.0
0.,0.087,0.087   ,100.0,50.0
14.28,0.122,0.087,100.0,50.0
46.43,0.139,0.087,100.0,50.0 
0.,0.169,0.169   ,100.0,50.0
25.,0.20,0.169   ,100.0,50.0
82.14,0.217,0.169,100.0,50.0
0.,0.23,0.23     ,100.0,50.0
35.7,0.265,0.23  ,100.0,50.0
164.3,0.287,0.23 ,100.0,50.0
0.,0.087,0.087   ,0.0,150.0
14.28,0.122,0.087,0.0,150.0
46.43,0.139,0.087,0.0,150.0 
0.,0.169,0.169   ,0.0,150.0
25.,0.20,0.169   ,0.0,150.0
82.14,0.217,0.169,0.0,150.0
0.,0.23,0.23     ,0.0,150.0
35.7,0.265,0.23  ,0.0,150.0
164.3,0.287,0.23 ,0.0,150.0
0.,0.087,0.087   ,100.0,150.0
14.28,0.122,0.087,100.0,150.0
46.43,0.139,0.087,100.0,150.0 
0.,0.169,0.169   ,100.0,150.0
25.,0.20,0.169   ,100.0,150.0
82.14,0.217,0.169,100.0,150.0
0.,0.23,0.23     ,100.0,150.0
35.7,0.265,0.23  ,100.0,150.0
164.3,0.287,0.23 ,100.0,150.0
*EXPANSION
1.e-2,
*INITIAL CONDITIONS,TYPE=FIELD,VARIABLE=1
NALL,120.0
*INITIAL CONDITIONS, TYPE=TEMPERATURE
NALL, 0.0
*******************************************
*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.1
*BOUNDARY
1,1,2
2,1,2
*******************************************
*STEP,INC=10000,NLGEOM
Test Normal Behavior of Gasket
*STATIC
.01,20.,,0.01
*CLOAD,AMP=HIST
3,2,10.47236
4,2,11.51936
*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, 
*END STEP
*******************************************
*STEP
Test Membrane Behavior of Element
*STATIC
*BOUNDARY,OP=NEW
1,1,2,0.0
2,2,2,0.0
3,1,2,0.0
4,2,2,0.0
*CLOAD,OP=NEW
2,1,-50.0
4,1,-50.0
*END STEP
*******************************************
*STEP
Test Transverse Shear Stiffness of Element
*STATIC
*BOUNDARY,OP=NEW
1,1,2,0.0
3,2,2,0.0
2,1,2,0.0
4,2,2,0.0
*CLOAD,OP=NEW
3,1,-4.0608
4,1,-7.2936
*END STEP
*******************************************
*STEP
Thermal Expansion of Element
*STATIC
*BOUNDARY,OP=NEW
1,2
2,2
*CLOAD,OP=NEW
*TEMPERATURE,OP=NEW
NALL,100.
*END STEP
*******************************************