*HEADING
:*STEADY STATE DYNAMIC,DIRECT RESPONSE WITH VISCOELASTICITY;
SIMULATE 1-D CASE; GASKET BEHAVIOR MODEL FOR THE GASKET BEHAVIOR;
[GKPE4,GKPS4,GKPS4N,GKPE6,GKPS6,GKPS6N]
*PREPRINT,MODEL=YES
*RESTART,WRITE,FREQ=20
*NODE,NSET=N1
1,0.,0.
2,1.,0.
11,0.,1.
12,1.,1.
*NODE,NSET=N1000
1001,10.,0.
1002,11.,0.
1011,10.,1.
1012,11.,1.
*NODE,NSET=N2000
2001,20.,0.
2002,21.,0.
2011,20.,1.
2012,21.,1.
**
2020,21.,0.
2120,21.,1.
*NODE,NSET=N3000
3001,30.,0.
3002,31.,0.
3011,30.,0.5
3012,31.,0.5
3021,30.,1.
3022,31.,1.
*NODE,NSET=N4000
4001,40.,0.
4002,41.,0.
4011,40.,0.5
4012,41.,0.5
4021,40.,1.
4022,41.,1.
*NODE,NSET=N5000
5001,50.,0.
5002,51.,0.
5011,50.,0.5
5012,51.,0.5
5021,50.,1.
5022,51.,1.
**
5020,51.,0.
5120,51.,0.5
5220,51.,1.
*ELEMENT,TYPE=GKPE4,ELSET=SP1
1,11,1,12,2
*ELEMENT,TYPE=GKPS4,ELSET=SP1000
1001,1011,1001,1012,1002
*ELEMENT,TYPE=GKPS4N,ELSET=SP2000
2001,2011,2001,2012,2002
*ELEMENT,TYPE=GKPE6,ELSET=SP3000
3001,3021,3011,3001,3022,3012,3002
*ELEMENT,TYPE=GKPS6,ELSET=SP4000
4001,4021,4011,4001,4022,4012,4002
*ELEMENT,TYPE=GKPS6N,ELSET=SP5000
5001,5021,5011,5001,5022,5012,5002
*ELSET,ELSET=SP
SP1,SP1000,SP2000,SP3000,SP4000,SP5000
*GASKET SECTION,ELSET=SP,BEHAVIOR=BE
1.2,
*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
3,12
*ELEMENT,TYPE=MASS,ELSET=MA1000
1002,1002
1003,1012
*ELEMENT,TYPE=MASS,ELSET=MA2000
2002,2020
2003,2120
*ELSET,ELSET=MA
MA1,MA1000,MA2000
*MASS,ELSET=MA
0.25, 
*ELEMENT,TYPE=MASS,ELSET=MA3000
3002,3002
3003,3012
3004,3022
*ELEMENT,TYPE=MASS,ELSET=MA4000
4002,4002
4003,4012
4004,4022
*ELEMENT,TYPE=MASS,ELSET=MA5000
5002,5020
5003,5120
5004,5220
*ELSET,ELSET=MASS
MA3000,MA4000,MA5000
*MASS,ELSET=MASS
0.16666, 
****************************************
*MPC
TIE,2020,2002
TIE,2120,2012
**
TIE,5020,5022
TIE,5120,5022
TIE,5220,5022
****************************************
*EQUATION
2, 
12,1,1.,2,1,-1.
**
2, 
1012,1,1.,1002,1,-1.
**
2, 
2012,1,1.,2002,1,-1.
**
2, 
3002,1,1.,3022,1,-1.
2, 
3012,1,1.,3022,1,-1.
**
2, 
4002,1,1.,4022,1,-1.
2, 
4012,1,1.,4022,1,-1.
**
2, 
5002,1,1.,5022,1,-1.
2, 
5012,1,1.,5022,1,-1.
*BOUNDARY
1,1,2,0
2,2,2,0
11,1,2,0
12,2,2,0
**
1001,1,2,0
1002,2,2,0
1011,1,2,0
1012,2,2,0
**
2001,1,1,0
2011,1,1,0
**
3001,1,2,0
3002,2,2,0
3011,1,2,0
3012,2,2,0
3021,1,2,0
3022,2,2,0
**
4001,1,2,0
4002,2,2,0
4011,1,2,0
4012,2,2,0
4021,1,2,0
4022,2,2,0
**
5001,1,1,0
5011,1,1,0
5021,1,1,0
****************************************
*STEP,NLGEOM
*STATIC
1.,1.
*BOUNDARY
2,1,1,-0.25
1002,1,1,-0.25
2002,1,1,-0.25
3022,1,1,-0.25
4022,1,1,-0.25
5022,1,1,-0.25
*EL PRINT
E,
S,
*OUTPUT,FIELD
*ELEMENT OUTPUT
E,
S,
*END STEP
****************************************
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
8.4
*BOUNDARY
2,1,1,-0.05
1002,1,1,-0.05
2002,1,1,-0.05
3022,1,1,-0.05
4022,1,1,-0.05
5022,1,1,-0.05
*EL PRINT
E,
S,
*OUTPUT,FIELD
*ELEMENT OUTPUT
E,
S,
*END STEP
****************************************
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
8.4
*BOUNDARY,LOAD CASE=2
2,1,1,-0.05
1002,1,1,-0.05
2002,1,1,-0.05
3022,1,1,-0.05
4022,1,1,-0.05
5022,1,1,-0.05
*EL PRINT
E,
S,
*OUTPUT,FIELD
*ELEMENT OUTPUT
E,
S,
*END STEP
****************************************
*STEP,NLGEOM
*STATIC
1.,1.
*BOUNDARY
2,1,1,-0.75
1002,1,1,-0.75
2002,1,1,-0.75
3022,1,1,-0.75
4022,1,1,-0.75
5022,1,1,-0.75
*EL PRINT
E,
S,
*OUTPUT,FIELD
*ELEMENT OUTPUT
E,
S,
*END STEP
****************************************
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
81.55
*BOUNDARY
2,1,1,-0.05
1002,1,1,-0.05
2002,1,1,-0.05
3022,1,1,-0.05
4022,1,1,-0.05
5022,1,1,-0.05
*EL PRINT
E,
S,
*OUTPUT,FIELD
*ELEMENT OUTPUT
E,
S,
*END STEP
****************************************
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
81.55
*BOUNDARY,LOAD CASE=2
2,1,1,-0.05
1002,1,1,-0.05
2002,1,1,-0.05
3022,1,1,-0.05
4022,1,1,-0.05
5022,1,1,-0.05
*EL PRINT
E,
S,
*OUTPUT,FIELD
*ELEMENT OUTPUT
E,
S,
*END STEP
****************************************
*STEP,NLGEOM
*STATIC
1.,1.
*BOUNDARY
2,1,1,-0.25
1002,1,1,-0.25
2002,1,1,-0.25
3022,1,1,-0.25
4022,1,1,-0.25
5022,1,1,-0.25
*EL PRINT
E,
S,
*OUTPUT,FIELD
*ELEMENT OUTPUT
E,
S,
*END STEP
****************************************
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
81.55
*BOUNDARY
2,1,1,-0.05
1002,1,1,-0.05
2002,1,1,-0.05
3022,1,1,-0.05
4022,1,1,-0.05
5022,1,1,-0.05
*EL PRINT
E,
S,
*OUTPUT,FIELD
*ELEMENT OUTPUT
E,
S,
*END STEP
****************************************
*STEP,PERTURBATION
*STEADY STATE DYNAMICS,DIRECT
81.55
*BOUNDARY,LOAD CASE=2
2,1,1,-0.05
1002,1,1,-0.05
2002,1,1,-0.05
3022,1,1,-0.05
4022,1,1,-0.05
5022,1,1,-0.05
*EL PRINT
E,
S,
*OUTPUT,FIELD
*ELEMENT OUTPUT
E,
S,
*END STEP