*HEADING
:*MODAL DYNAMIC RESPONSE UNDER DYNAMIC FORCE;
SIMULATE 1-D CASE; GASKET BEHAVIOR MODEL FOR THE GASKET BEHAVIOR;
[GK3D4L,GK3D4LN,GK3D6L,GK3D6LN]
*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.
**
1021,11.,1.
1022,11.,0.
*NODE,NSET=N2000
2001,20.,0.
2002,20.,0.5
2003,20.,1.
2011,21.,0.
2012,21.,0.5
2013,21.,1.
*NODE,NSET=N3000
3001,30.,0.
3002,30.,0.5
3003,30.,1.
3011,31.,0.
3012,31.,0.5
3013,31.,1.
**
3110,31.,0.
3120,31.,0.5
3130,31.,1.
*NSET, NSET=NALL
N1, N1000, N2000, N3000
*ELEMENT,TYPE=GK3D4L,ELSET=SP1
1,11,1,12,2
*ELEMENT,TYPE=GK3D4LN,ELSET=SP1000
1001,1011,1001,1012,1002
*ELEMENT,TYPE=GK3D6L,ELSET=SP2000
2001,2003,2002,2001,2013,2012,2011
*ELEMENT,TYPE=GK3D6LN,ELSET=SP3000
3001,3003,3002,3001,3013,3012,3011
*ELSET,ELSET=SP
SP1,SP1000,SP2000,SP3000
*GASKET SECTION,ELSET=SP,BEHAVIOR=BE
1.2,
*GASKET BEHAVIOR,NAME=BE
*GASKET THICKNESS BEHAVIOR,TYPE=DAMAGE,DIRECTION=LOADING,VARIABLE=FORCE
0.,0.
200.,1.
*GASKET CONTACT AREA
1.5,0.
1.8,0.05
*ELEMENT,TYPE=MASS,ELSET=MA1
2,2
3,12
*ELEMENT,TYPE=MASS,ELSET=MA1000
1002,1021
1003,1022
*ELSET,ELSET=MA
MA1,MA1000
*MASS,ELSET=MA
0.25, 
*ELEMENT,TYPE=MASS,ELSET=MA2000
2002,2011
2003,2012
2004,2013
*ELEMENT,TYPE=MASS,ELSET=MA3000
3002,3110
3003,3120
3004,3130
*ELSET,ELSET=MASS
MA2000,MA3000
*MASS,ELSET=MASS
0.16666, 
*AMPLITUDE,NAME=COS,DEFINITION=PERIODIC
1,10.,0.,0.
10.,0.
** N,w,t_0,A_0 | A_1,B_1  ==>  P = P_0 cos(wt) = 10 cos(10t)
****************************************
*MPC
TIE,1022,1002
TIE,1021,1012
**
TIE,3110,3011
TIE,3120,3012
TIE,3130,3013
****************************************
*EQUATION
2, 
12,1,1.,2,1,-1.
**
2, 
1012,1,1.,1002,1,-1
**
2, 
2012,1,1.,2011,1,-1.
2, 
2013,1,1.,2011,1,-1.
**
2, 
3012,1,1.,3011,1,-1.
2, 
3013,1,1.,3011,1,-1.
*BOUNDARY
1,1,3,0
2,2,3,0
11,1,3,0
12,2,3,0
**
1001,1,1,0
1011,1,1,0
**
2001,1,3,0
2002,1,3,0
2003,1,3,0
2011,2,3,0
2012,2,3,0
2013,2,3,0
**
3001,1,1,0
3002,1,1,0
3003,1,1,0
****************************************
*STEP
*FREQUENCY, EIGENSOLVER=LANCZOS, SIM=NO
4,1.0
*OUTPUT,FIELD,FREQ=0
*MODAL FILE
*END STEP
****************************************
*STEP
*MODAL DYNAMIC
1e-3,1e0,
*MODAL DAMPING
1,4
*CLOAD, AMP=COS
2,1,-5.
12,1,-5.
**
1002,1,-5.
1012,1,-5.
**
2011,1,-3.3333
2012,1,-3.3333
2013,1,-3.3333
**
3011,1,-3.3333
3012,1,-3.3333
3013,1,-3.3333
*EL PRINT,FREQ=500
S,
E,
*EL FILE,FREQ=50
S,
E,
*NODE PRINT,FREQ=500
U,
RF, 
*NODE FILE,FREQ=50
U,
RF, 
*ENERGY PRINT,FREQ=500
*OUTPUT,FIELD,OP=NEW,FREQ=50
*ELEMENT OUTPUT,ELSET=SP
S,E
*NODE OUTPUT, NSET=NALL
U,RF
*OUTPUT,HISTORY,FREQ=50
*ELEMENT OUTPUT,ELSET=SP
S,E
*NODE OUTPUT, NSET=NALL
U,RF
*END STEP