*HEADING
: PERTURBATION ANALYSIS/MULTIPLE LOAD CASES FOR ELASTIC_PLASTIC MATERIAL WITH LOW
          INITIAL STIFFNESS, THICKNESS DIRECTION BEHAVIOR ONLY [GK2D2N]
*NODE,NSET=NALL
1,
2,,1.
*ELEMENT,TYPE=GK2D2N,ELSET=EINT
1,1,2
*GASKET SECTION,ELSET=EINT,BEHAVIOR=IS
,0.1
*GASKET BEHAVIOR,NAME=IS
*GASKET THICKNESS BEHAVIOR,DIRECTION=LOADING,VARIABLE=FORCE
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,VARIABLE=FORCE
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
*GASKET CONTACT AREA
1.,0.
1.5,0.3
*BOUNDARY
1,1
*STEP,NLGEOM,PERTURBATION
*STATIC
*CLOAD
2,1,-0.010298
*EL PRINT,FREQ=0
S, 
CS11, 
E, 
EE, 
*EL FILE,FREQ=18
S, 
CS11, 
E, 
EE, 
*OUTPUT,FIELD,FREQ=18,VAR=PRE
*output,history
*element output,elset=eint
ee,e,s
*END STEP
*STEP,NLGEOM
*STATIC
*CLOAD
2,1,-46.43
*NODE PRINT
U, 
RF, 
*EL PRINT,FREQ=0
S, 
CS11, 
E, 
EE, 
PE, 
EVOL, 
*EL FILE,FREQ=18
S, 
CS11, 
E, 
PE, 
EVOL, 
*output,history
*element output,elset=eint
ee,e,s
*END STEP
*STEP,NLGEOM,PERTURBATION
*STATIC
*CLOAD
2,1,-94.
*EL PRINT,FREQ=0
S, 
CS11, 
E, 
EE, 
*EL FILE,FREQ=18
S, 
CS11, 
E, 
EE, 
*output,history
*element output,elset=eint
ee,e,s
*END STEP
*STEP,NLGEOM
*STATIC
*CLOAD
2,1,-14.
*NODE PRINT,FREQ=0
U, 
RF, 
*EL PRINT,FREQ=0
S, 
CS11, 
E, 
EE, 
PE, 
EVOL, 
*NODE FILE,FREQ=18
U, 
RF, 
*EL FILE,FREQ=18
S, 
CS11, 
E, 
PE, 
EVOL, 
*output,history
*element output,elset=eint
ee,e,s
*END STEP
*STEP,NLGEOM,PERTURBATION
*STATIC
*CLOAD
2,1,-94.
*EL PRINT,FREQ=0
S, 
CS11, 
E, 
EE, 
*EL FILE,FREQ=18
S, 
CS11, 
E, 
EE, 
*output,history
*element output,elset=eint
ee,e,s
*END STEP

*STEP,NLGEOM,PERTURBATION
*STATIC
*LOAD CASE,NAME=PRESSURE-1
*CLOAD
2,1,-94.
*END LOAD CASE
*LOAD CASE,NAME=PRESSURE-2
*CLOAD
2,1,-150
*END LOAD CASE
*EL PRINT,FREQ=0
S, 
CS11, 
E, 
EE, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
EE
*output,history
*element output,elset=eint
ee,e,s
*END STEP