*HEADING
:*MODAL DYNAMIC RESPONSE UNDER DYNAMIC FORCE;
SIMULATE 1-D CASE; MATERIAL MODEL FOR THE GASKET BEHAVIOR;
[GKAX2,GKAX2N,GKAX4,GKAX4N,GKAX6,GKAX6N]
*RESTART,WRITE,FREQ=20
*NODE,NSET=N1
*NODE,NSET=N1000
**
*NODE,NSET=N2000
*NODE,NSET=N3000
**
*NODE,NSET=N4000
*NODE,NSET=N5000
**
*NSET, NSET=NALL
N1, N1000, N2000, N3000, N4000, N5000
*ELEMENT,TYPE=GKAX2,ELSET=SP1
*ELEMENT,TYPE=GKAX2N,ELSET=SP1000
*ELSET,ELSET=SP
SP1,SP1000
*ELEMENT,TYPE=GKAX4,ELSET=SP2000
*ELEMENT,TYPE=GKAX4N,ELSET=SP3000
*ELEMENT,TYPE=GKAX6,ELSET=SP4000
*ELEMENT,TYPE=GKAX6N,ELSET=SP5000
*ELSET,ELSET=SPR
SP2000,SP3000,SP4000,SP5000
*GASKET SECTION,ELSET=SP,MATERIAL=BE
*GASKET SECTION,ELSET=SPR,MATERIAL=BE
*MATERIAL,NAME=BE
*ELASTIC
240.,0.
*ELEMENT,TYPE=MASS,ELSET=MA1
*ELEMENT,TYPE=MASS,ELSET=MA1000
*ELSET,ELSET=MA
MA1,MA1000
*MASS,ELSET=MA
0.5,
*ELEMENT,TYPE=MASS,ELSET=MA2000
*ELEMENT,TYPE=MASS,ELSET=MA3000
*ELSET,ELSET=MAS
MA2000,MA3000
*MASS,ELSET=MAS
0.25,
*ELEMENT,TYPE=MASS,ELSET=MA4000
*ELEMENT,TYPE=MASS,ELSET=MA5000
*ELSET,ELSET=MASS
MA4000,MA5000
*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,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
**
**
**
****************************************
*BOUNDARY
1,1,2,0
2,1,1,0
**
1001,1,1,0
**
2001,1,2,0
2002,1,2,0
2011,1,1,0
2012,1,1,0
**
3001,1,1,0
**
**
****************************************
*STEP
*FREQUENCY, EIGENSOLVER=LANCZOS
6,1.
*OUTPUT,FIELD,FREQ=0
*MODAL FILE
*END STEP
****************************************
*STEP
*MODAL DYNAMIC
1e-3,1e0,
*CLOAD, AMP=COS
2,2,-10.
**
1002,1,-10
**
2011,2,-5.
2012,2,-5.
**
3011,1,-5.
3012,1,-5.
**
4011,2,-3.3333
4012,2,-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