*HEADING
:*MODAL DYNAMIC RESPONSE UNDER DYNAMIC FORCE;
SIMULATE 1-D CASE; GASKET BEHAVIOR MODEL FOR THE GASKET BEHAVIOR;
[GK2D2,GK2D2N,GK3D2,GK3D2N]
*RESTART,WRITE,FREQ=20
*NODE,NSET=N1
1,0.,0.
2,1.,0.
*NODE,NSET=N100
101,10.,0.
102,11.,0.
**
120,11.,0.
*NODE,NSET=N200
201,20.,0.
202,21.,0.
*NODE,NSET=N300
301,30.,0.
302,31.,0.
**
320,31.,0.
*NSET, NSET=NALL
N1, N100, N200, N300
*ELEMENT,TYPE=GK2D2,ELSET=SP1
1,1,2
*ELEMENT,TYPE=GK2D2N,ELSET=SP100
101,101,102
*ELEMENT,TYPE=GK3D2,ELSET=SP200
201,201,202
*ELEMENT,TYPE=GK3D2N,ELSET=SP300
301,301,302
*ELSET,ELSET=SP
SP1,SP100,SP200,SP300
*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=MA
2,2
102,120
202,202
302,320
*MASS,ELSET=MA
0.5, 
*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,120,102
TIE,320,302
*BOUNDARY
1,1,2
2,2,2
101,1,1
201,1,3
202,2,3
301,1,1
****************************************
*STEP
*FREQUENCY, EIGENSOLVER=LANCZOS
4,1.
*OUTPUT,FIELD,FREQ=0
*MODAL FILE
*END STEP
****************************************
*STEP
*MODAL DYNAMIC
1e-3,1e0,
*CLOAD, AMP=COS
2,1,-10
102,1,-10
202,1,-10
302,1,-10
*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