*HEADING
: DISPLACEMENT CONTROL WITH VISCOELASTIC MATERIAL 
          WITH DAMAGE [GK3D2, GK3D2N]
*NODE,NSET=NALL
1,0.,0.,1.
2,0.,1.,1.
3,0.,0.,2.
4,0.,1.,2.
5,0.,1.,2.
*ELEMENT,TYPE=GK3D2,ELSET=EINT
1,1,2
*ELEMENT,TYPE=GK3D2N,ELSET=EINT1
2,3,4
*GASKET SECTION,ELSET=EINT,BEHAVIOR=IS,ORIENTATION=LOCAL
,,,1.1,0.,1.,0.
*GASKET BEHAVIOR,NAME=IS
*GASKET ELASTICITY,COMPONENT=TRANSVERSE SHEAR,VARIABLE=FORCE
200.,
*GASKET THICKNESS BEHAVIOR,TYPE=DAMAGE,DIRECTION=LOADING,VARIABLE=FORCE
0.,0.
200.,0.5
300.,1.
*GASKET THICKNESS BEHAVIOR,TYPE=DAMAGE,DIRECTION=UNLOADING,VARIABLE=FORCE
0.,0.,0.75
50.,0.5,0.75
250.,0.75,0.75
*VISCOELASTIC,FREQUENCY=TABULAR,TYPE=TRACTION,PRELOAD=UNIAXIAL
4.644,412.84,1.,0.
3.1396,408.888,15.8,0.
2.1416,406.132,25.1,0.
1.4556,404.248,39.8,0.
1.0172,402.9528,63.1,0.
0.71,402.0464,100.,0.
**
4.644,412.84,1.,0.45
3.1396,408.888,15.8,0.45
2.1416,406.132,25.1,0.45
1.4556,404.248,39.8,0.45
1.0172,402.9528,63.1,0.45
0.71,402.0464,100.,0.45
**
2.322,206.42,1.,0.55
1.5698,204.444,15.8,0.55
1.0708,203.066,25.1,0.55
0.7278,202.124,39.8,0.55
0.5086,201.4764,63.1,0.55
0.355,201.0232,100.,0.55
**
2.322,206.42,1.,1.
1.5698,204.444,15.8,1.
1.0708,203.066,25.1,1.
0.7278,202.124,39.8,1.
0.5086,201.4764,63.1,1.
0.355,201.0232,100.,1.
*GASKET CONTACT AREA
1.0,0.0
1.5,1.0
*ORIENTATION,NAME=LOCAL,SYSTEM=RECT
1.0,0.0,0.0,0.0,0.0,-1.0
3,90.
*GASKET SECTION,ELSET=EINT1,BEHAVIOR=IS1
,,,1.1,0.,1.,0.
*GASKET BEHAVIOR,NAME=IS1
*GASKET ELASTICITY,COMPONENT=TRANSVERSE SHEAR,VARIABLE=FORCE
200.,
*GASKET THICKNESS BEHAVIOR,TYPE=DAMAGE,DIRECTION=LOADING,VARIABLE=FORCE
0.,0.
200.,0.5
300.,1.
*GASKET THICKNESS BEHAVIOR,TYPE=DAMAGE,DIRECTION=UNLOADING,VARIABLE=FORCE
0.,0.,0.75
50.,0.5,0.75
250.,0.75,0.75
*VISCOELASTIC,FREQUENCY=TABULAR,TYPE=TRACTION
1.161E-2,-3.21E-2,1.
7.849E-3,-2.222E-2,15.8
5.354E-3,-1.533E-2,25.1
3.639E-3,-1.062E-2,39.8
2.543E-3,-7.382E-3,63.1
1.775E-3,-5.116E-3,100.0
*GASKET CONTACT AREA
1.0,0.0
1.5,1.0
*ELEMENT,TYPE=MASS,ELSET=MASS
3,2
4,5
*MASS,ELSET=MASS
0.5
*BOUNDARY
1,1,3
3,1,1
*MPC
TIE,5,4
*******************************************
*STEP,NLGEOM
*STATIC
1.,1.
*BOUNDARY
2,2,2,-0.25
4,1,1,-0.25
*EL PRINT
S, 
E,  
CS11, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
8.4
*BOUNDARY
2,2,2,-0.05
4,1,1,-0.05
*EL PRINT
S, 
E,  
CS11, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
8.4
*BOUNDARY,LOAD CASE=2
2,2,2,-0.05
4,1,1,-0.05
*EL PRINT
S, 
E,  
CS11, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*END STEP
*STEP,NLGEOM
*STATIC
1.,1.
*BOUNDARY
2,2,2,-0.75
4,1,1,-0.75
*EL PRINT
S, 
E,  
CS11, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
81.55
*BOUNDARY
2,2,2,-0.05
4,1,1,-0.05
*EL PRINT
S, 
E,  
CS11, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
81.55
*BOUNDARY,LOAD CASE=2
2,2,2,-0.05
4,1,1,-0.05
*EL PRINT
S, 
E,  
CS11, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*END STEP
*STEP,NLGEOM
*STATIC
1.,1.
*BOUNDARY
2,2,2,-0.25
4,1,1,-0.25
*EL PRINT
S, 
E,  
CS11, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
81.55
*BOUNDARY
2,2,2,-0.05
4,1,1,-0.05
*EL PRINT
S, 
E,  
CS11, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
81.55
*BOUNDARY,LOAD CASE=2
2,2,2,-0.05
4,1,1,-0.05
*EL PRINT
S, 
E,  
CS11, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*END STEP