*HEADING
:*STEADY STATE DYNAMIC,DIRECT RESPONSE WITH VISCOELASTICITY;
SIMULATE 1-D CASE; GASKET BEHAVIOR MODEL FOR THE GASKET BEHAVIOR;
[GKAX4,GKAX4N,GKAX6,GKAX6N]
*RESTART,WRITE,FREQ=20
*NODE,NSET=N2000
2001,0.     ,20.
2002,0.56419,20.
2011,0.     ,21.
2012,0.56419,21.
*NODE,NSET=N3000
3001,0.     ,30.
3002,0.56419,30.
3011,0.     ,31.
3012,0.56419,31.
**
3010,0.     ,30.
3020,0.56419,30.
3110,0.     ,31.
3120,0.56419,31.
*NODE,NSET=N4000
4001,0.      ,40.
4002,0.282095,40.
4003,0.56419 ,40.
4011,0.      ,41.
4012,0.282095,41.
4013,0.56419 ,41.
*NODE,NSET=N5000
5001,0.      ,50.
5002,0.282095,50.
5003,0.56419 ,50.
5011,0.      ,51.
5012,0.282095,51.
5013,0.56419 ,51.
**
5010,0.      ,50.
5020,0.282095,50.
5030,0.56419 ,50.
5110,0.      ,51.
5120,0.282095,51.
5130,0.56419 ,51.
*ELEMENT,TYPE=GKAX4,ELSET=SP2000
2001,2001,2002,2011,2012
*ELEMENT,TYPE=GKAX4N,ELSET=SP3000
3001,3001,3002,3011,3012
*ELEMENT,TYPE=GKAX6,ELSET=SP4000
4001,4001,4002,4003,4011,4012,4013
*ELEMENT,TYPE=GKAX6N,ELSET=SP5000
5001,5001,5002,5003,5011,5012,5013
*ELSET,ELSET=SPR
SP2000,SP3000,SP4000,SP5000
*GASKET SECTION,ELSET=SP,BEHAVIOR=BE
1.2,,,0.56419,0,1,0
*GASKET SECTION,ELSET=SPR,BEHAVIOR=BE
1.2,,,,0,1,0
*GASKET BEHAVIOR,NAME=BE
*GASKET THICKNESS BEHAVIOR,TYPE=DAMAGE,DIRECTION=LOADING,VARIABLE=STRESS
0.,0.
200.,0.5
300.,1.
*GASKET THICKNESS BEHAVIOR,TYPE=DAMAGE,DIRECTION=UNLOADING,VARIABLE=STRESS
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
*ELEMENT,TYPE=MASS,ELSET=MA1
2,2
*ELEMENT,TYPE=MASS,ELSET=MA1000
1002,1012
*ELSET,ELSET=MA
MA1,MA1000
*MASS,ELSET=MA
0.5, 
*ELEMENT,TYPE=MASS,ELSET=MA2000
2002,2011
2003,2012
*ELEMENT,TYPE=MASS,ELSET=MA3000
3002,3110
3003,3120
*ELSET,ELSET=MAS
MA2000,MA3000
*MASS,ELSET=MAS
0.25, 
*ELEMENT,TYPE=MASS,ELSET=MA4000
4002,4011
4003,4012
4004,4013
*ELEMENT,TYPE=MASS,ELSET=MA5000
5002,5110
5003,5120
5004,5130
*ELSET,ELSET=MASS
MA4000,MA5000
*MASS,ELSET=MASS
0.16666, 
*MPC
TIE,1011,1001
TIE,1012,1002
**
TIE,3010,3001
TIE,3020,3002
TIE,3110,3011
TIE,3120,3012
**
TIE,5010,5001
TIE,5020,5002
TIE,5030,5003
TIE,5110,5011
TIE,5120,5012
TIE,5130,5013
****************************************
*EQUATION
2, 
2012,2,1.,2011,2,-1.
**
2, 
3012,1,1.,3011,1,-1.
**
2, 
4012,2,1.,4011,2,-1.
2, 
4013,2,1.,4011,2,-1.
**
2, 
5012,1,1.,5011,1,-1.
2, 
5013,1,1.,5011,1,-1.
****************************************
*BOUNDARY
2001,1,2,0
2002,1,2,0
2011,1,1,0
2012,1,1,0
**
3001,1,1,0
3002,1,1,0
**
4001,1,2,0
4002,1,2,0
4003,1,2,0
4011,1,1,0
4012,1,1,0
4013,1,1,0
**
5001,1,1,0
5002,1,1,0
5003,1,1,0
****************************************
*STEP
*STATIC
1.,1.
*BOUNDARY
2011,2,2,-0.25
3011,1,1,-0.25
4011,2,2,-0.25
5011,1,1,-0.25
*EL PRINT,ELSET=SPR
S,
E,
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMIC,DIRECT
8.4
*BOUNDARY
2011,2,2,-0.05
3011,1,1,-0.05
4011,2,2,-0.05
5011,1,1,-0.05
*EL PRINT,ELSET=SPR
S,
E,
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMIC,DIRECT
8.4
*BOUNDARY,LOAD CASE=2
2011,2,2,-0.05
3011,1,1,-0.05
4011,2,2,-0.05
5011,1,1,-0.05
*EL PRINT,ELSET=SPR
S,
E,
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*END STEP
*STEP
*STATIC
1.,1.
*BOUNDARY
2011,2,2,-0.75
3011,1,1,-0.75
4011,2,2,-0.75
5011,1,1,-0.75
*EL PRINT,ELSET=SPR
S,
E,
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMIC,DIRECT
81.55
*BOUNDARY
2011,2,2,-0.05
3011,1,1,-0.05
4011,2,2,-0.05
5011,1,1,-0.05
*EL PRINT,ELSET=SPR
S,
E,
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMIC,DIRECT
81.55
*BOUNDARY,LOAD CASE=2
2011,2,2,-0.05
3011,1,1,-0.05
4011,2,2,-0.05
5011,1,1,-0.05
*EL PRINT,ELSET=SPR
S,
E,
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*END STEP
*STEP
*STATIC
1.,1.
*BOUNDARY
2011,2,2,-0.25
3011,1,1,-0.25
4011,2,2,-0.25
5011,1,1,-0.25
*EL PRINT,ELSET=SPR
S,
E,
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMIC,DIRECT
81.55
*BOUNDARY
2011,2,2,-0.05
3011,1,1,-0.05
4011,2,2,-0.05
5011,1,1,-0.05
*EL PRINT,ELSET=SPR
S,
E,
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*END STEP
*STEP,PERTURBATION
*STEADY STATE DYNAMIC,DIRECT
81.55
*BOUNDARY,LOAD CASE=2
2011,2,2,-0.05
3011,1,1,-0.05
4011,2,2,-0.05
5011,1,1,-0.05
*EL PRINT,ELSET=SPR
S,
E,
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*END STEP