*HEADING
:*MODAL DYNAMIC RESPONSE UNDER DYNAMIC FORCE;
SIMULATE 1-D CASE; GASKET BEHAVIOR MODEL FOR THE GASKET BEHAVIOR;
[GK3D6,GK3D6N,GK3D8,GK3D8N,GK3D12M,GK3D12MN,GK3D18,GK3D18N]
*RESTART,WRITE,FREQ=20
*NODE,NSET=N1
*NODE,NSET=N1000
**
*NODE,NSET=N2000
*NODE,NSET=N3000
**
*NODE,NSET=N4000
**
*NODE,NSET=N5000
**
**
*NODE,NSET=N6000
**
**
*NODE,NSET=N7000
**
**
**
*NSET, NSET=NALL
N1, N1000, N2000, N3000, N4000
N5000, N6000, N7000
*ELEMENT,TYPE=GK3D6,ELSET=SP1
*ELEMENT,TYPE=GK3D6N,ELSET=SP1000
*ELEMENT,TYPE=GK3D8,ELSET=SP2000
*ELEMENT,TYPE=GK3D8N,ELSET=SP3000
*ELEMENT,TYPE=GK3D12M,ELSET=SP4000
*ELEMENT,TYPE=GK3D12MN,ELSET=SP5000
*ELEMENT,TYPE=GK3D18,ELSET=SP6000
*ELEMENT,TYPE=GK3D18N,ELSET=SP7000
*ELSET,ELSET=SP
SP1,SP1000,SP2000,SP3000,SP4000,SP5000,SP6000,SP7000
*GASKET SECTION,ELSET=SP,BEHAVIOR=BE
*GASKET BEHAVIOR,NAME=BE
*GASKET THICKNESS BEHAVIOR,TYPE=DAMAGE,DIRECTION=LOADING,VARIABLE=STRESS
*ELEMENT,TYPE=MASS,ELSET=MA1
*ELEMENT,TYPE=MASS,ELSET=MA1000
*ELEMENT,TYPE=MASS,ELSET=MA2000
*ELEMENT,TYPE=MASS,ELSET=MA3000
*ELSET,ELSET=MA
MA1,MA1000,MA2000,MA3000
*MASS,ELSET=MA
0.125,
*ELEMENT,TYPE=MASS,ELSET=MA4000
*ELEMENT,TYPE=MASS,ELSET=MA5000
*ELEMENT,TYPE=MASS,ELSET=MA6000
*ELEMENT,TYPE=MASS,ELSET=MA7000
*ELSET,ELSET=MASS
MA4000,MA5000,MA6000,MA7000
*MASS,ELSET=MASS
0.0555,
*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,1003,1002
TIE,1013,1012
TIE,1103,1102
TIE,1113,1112
**
TIE,3003,3002
TIE,3013,3012
TIE,3103,3102
TIE,3113,3112
**
TIE,5202,5002
TIE,5212,5012
TIE,5702,5102
TIE,5712,5112
TIE,5215,5015
TIE,5216,5016
TIE,5217,5017
TIE,5218,5018
TIE,5219,5019
**
TIE,7202,7002
TIE,7212,7012
TIE,7702,7102
TIE,7712,7112
TIE,7250,7050
TIE,7260,7060
TIE,7270,7070
TIE,7280,7080
TIE,7621,7601
**********************************
*EQUATION
**
**
**
**********************************
*BOUNDARY
  1,1,3
  2,2,3
 11,1,3
 12,2,3
101,1,3
102,2,3
111,1,3
112,2,3
**
**
**
**********************************
*STEP
*FREQUENCY, eigensolver=lanczos, SIM=NO
8, , ,10.
*OUTPUT,FIELD,FREQ=0
*MODAL FILE
*END STEP
****************************************
*STEP
*MODAL DYNAMIC
1e-3,1e0,
*MODAL DAMPING
*SELECT EIGENMODES, GENERATE
*CLOAD, AMP=COS
  2,1,-2.5
 12,1,-2.5
102,1,-2.5
112,1,-2.5
**
1002,1,-2.5
1012,1,-2.5
1102,1,-2.5
1112,1,-2.5
**
**
**
*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