*HEADING
:*MODAL DYNAMIC RESPONSE UNDER DYNAMIC FORCE;
SIMULATE 1-D CASE; MATERIAL MODEL FOR THE GASKET BEHAVIOR;
[GK3D6,GK3D6N,GK3D8,GK3D8N,GK3D12M,GK3D12MN,GK3D18,GK3D18N]
*RESTART,WRITE,FREQ=20
*NODE,NSET=N1
  1,0.,0.,0.
  2,1.,0.,0.
 11,0.,1.,0.
 12,1.,1.,0.
101,0.,0.,-1.
102,1.,0.,-1.
111,0.,1.,-1.
112,1.,1.,-1.
*NODE,NSET=N1000
1001,10.,0.,0.
1002,11.,0.,0.
1011,10.,1.,0.
1012,11.,1.,0.
1101,10.,0.,-1.
1102,11.,0.,-1.
1111,10.,1.,-1.
1112,11.,1.,-1.
**
1003,11.,0.,0.
1013,11.,1.,0.
1103,11.,0.,-1.
1113,11.,1.,-1.
*NODE,NSET=N2000
2001,20.,0.,0.
2002,21.,0.,0.
2011,20.,1.,0.
2012,21.,1.,0.
2101,20.,0.,-1.
2102,21.,0.,-1.
2111,20.,1.,-1.
2112,21.,1.,-1.
*NODE,NSET=N3000
3001,30.,0.,0.
3002,31.,0.,0.
3011,30.,1.,0.
3012,31.,1.,0.
3101,30.,0.,-1.
3102,31.,0.,-1.
3111,30.,1.,-1.
3112,31.,1.,-1.
**
3003,31.,0.,0.
3013,31.,1.,0.
3103,31.,0.,-1.
3113,31.,1.,-1.
*NODE,NSET=N4000
4001,40.,0.,0.
4002,41.,0.,0.
4011,40.,1.,0.
4012,41.,1.,0.
4101,40.,0.,-1.
4102,41.,0.,-1.
4111,40.,1.,-1.
4112,41.,1.,-1.
**
4005,40.,0.5, 0.0
4006,40.,0.0,-0.5
4007,40.,0.5,-0.5
4008,40.,0.5,-1.0
4009,40.,1.0,-0.5
4015,41.,0.5, 0.0
4016,41.,0.0,-0.5
4017,41.,0.5,-0.5
4018,41.,0.5,-1.0
4019,41.,1.0,-0.5
*NODE,NSET=N5000
5001,50.,0.,0.
5002,51.,0.,0.
5011,50.,1.,0.
5012,51.,1.,0.
5101,50.,0.,-1.
5102,51.,0.,-1.
5111,50.,1.,-1.
5112,51.,1.,-1.
**
5005,50.,0.5, 0.0
5006,50.,0.0,-0.5
5007,50.,0.5,-0.5
5008,50.,0.5,-1.0
5009,50.,1.0,-0.5
5015,51.,0.5, 0.0
5016,51.,0.0,-0.5
5017,51.,0.5,-0.5
5018,51.,0.5,-1.0
5019,51.,1.0,-0.5
**
5202,51.,0.,0.
5212,51.,1.,0.
5702,51.,0.,-1.
5712,51.,1.,-1.
5215,51.,0.5, 0.0
5216,51.,0.0,-0.5
5217,51.,0.5,-0.5
5218,51.,0.5,-1.0
5219,51.,1.0,-0.5
*NODE,NSET=N6000
6001,60.,0.,0.
6002,61.,0.,0.
6011,60.,1.,0.
6012,61.,1.,0.
6101,60.,0.,-1.
6102,61.,0.,-1.
6111,60.,1.,-1.
6112,61.,1.,-1.
**
6005,60.,0.5, 0.0
6006,60.,0.0,-0.5
6007,60.,0.5,-1.0
6008,60.,1.0,-0.5
6050,61.,0.5, 0.0
6060,61.,0.0,-0.5
6070,61.,0.5,-1.0
6080,61.,1.0,-0.5
**
6600,60.,0.5,-0.5
6601,61.,0.5,-0.5
*NODE,NSET=N7000
7001,70.,0.,0.
7002,71.,0.,0.
7011,70.,1.,0.
7012,71.,1.,0.
7101,70.,0.,-1.
7102,71.,0.,-1.
7111,70.,1.,-1.
7112,71.,1.,-1.
**
7005,70.,0.5, 0.0
7006,70.,0.0,-0.5
7007,70.,0.5,-1.0
7008,70.,1.0,-0.5
7050,71.,0.5, 0.0
7060,71.,0.0,-0.5
7070,71.,0.5,-1.0
7080,71.,1.0,-0.5
**
7600,70.,0.5,-0.5
7601,71.,0.5,-0.5
**
7202,71.,0.,0.
7212,71.,1.,0.
7702,71.,0.,-1.
7712,71.,1.,-1.
7250,71.,0.5, 0.0
7260,71.,0.0,-0.5
7270,71.,0.5,-1.0
7280,71.,1.0,-0.5
7621,71.,0.5,-0.5
*NSET, NSET=NALL
N1, N1000, N2000, N3000, N4000
N5000, N6000, N7000
*ELEMENT,TYPE=GK3D6,ELSET=SP1
 1,11,1,101,12,2,102
11,101,111,11,102,112,12
*ELEMENT,TYPE=GK3D6N,ELSET=SP1000
1001,1011,1001,1101,1012,1002,1102
1011,1101,1111,1011,1102,1112,1012
*ELEMENT,TYPE=GK3D8,ELSET=SP2000
2001,2011,2001,2101,2111,2012,2002,2102,2112
*ELEMENT,TYPE=GK3D8N,ELSET=SP3000
3001,3011,3001,3101,3111,3012,3002,3102,3112
*ELEMENT,TYPE=GK3D12M,ELSET=SP4000
4001,4011,4001,4101,4005,4006,4007,4012,4002,4102,4015,4016,4017
4011,4101,4111,4011,4008,4009,4007,4102,4112,4012,4018,4019,4017
*ELEMENT,TYPE=GK3D12MN,ELSET=SP5000
5001,5011,5001,5101,5005,5006,5007,5012,5002,5102,5015,5016,5017
5011,5101,5111,5011,5008,5009,5007,5102,5112,5012,5018,5019,5017
*ELEMENT,TYPE=GK3D18,ELSET=SP6000
6001,6011,6001,6101,6111,6005,6006,6007,6008,6012,6002,6102,6112,
6050,6060,6070,6080,6600,6601
*ELEMENT,TYPE=GK3D18N,ELSET=SP7000
7001,7011,7001,7101,7111,7005,7006,7007,7008,7012,7002,7102,7112,
7050,7060,7070,7080,7600,7601
*ELSET,ELSET=SP
SP1,SP1000,SP2000,SP3000,SP4000,SP5000,SP6000,SP7000
*GASKET SECTION,ELSET=SP,MATERIAL=BE
1.2,
*MATERIAL,NAME=BE
*ELASTIC
240.,0
*ELEMENT,TYPE=MASS,ELSET=MA1
2,2
3,12
4,102
5,112
*ELEMENT,TYPE=MASS,ELSET=MA1000
1002,1003
1003,1013
1004,1103
1005,1113
*ELEMENT,TYPE=MASS,ELSET=MA2000
2002,2002
2003,2012
2004,2102
2005,2112
*ELEMENT,TYPE=MASS,ELSET=MA3000
3002,3003
3003,3013
3004,3103
3005,3113
*ELSET,ELSET=MA
MA1,MA1000,MA2000,MA3000
*MASS,ELSET=MA
0.125, 
*ELEMENT,TYPE=MASS,ELSET=MA4000
4002,4002
4003,4012
4004,4102
4005,4112
4006,4015
4007,4016
4008,4017
4009,4018
4010,4019
*ELEMENT,TYPE=MASS,ELSET=MA5000
5002,5202
5003,5212
5004,5702
5005,5712
5006,5215
5007,5216
5008,5217
5009,5218
5010,5219
*ELEMENT,TYPE=MASS,ELSET=MA6000
6002,6002
6003,6012
6004,6102
6005,6112
6006,6050
6007,6060
6008,6070
6009,6080
6010,6601
*ELEMENT,TYPE=MASS,ELSET=MA7000
7002,7202
7003,7212
7004,7702
7005,7712
7006,7250
7007,7260
7008,7270
7009,7280
7010,7621
*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
2, 
 12,1,1., 2,1,-1.
2, 
102,1,1., 2,1,-1.
2, 
112,1,1., 2,1,-1.
**
2, 
1012,1,1.,1002,1,-1.
2, 
1102,1,1.,1002,1,-1.
2, 
1112,1,1.,1002,1,-1.
**
2, 
2012,1,1.,2002,1,-1.
2, 
2102,1,1.,2002,1,-1.
2, 
2112,1,1.,2002,1,-1.
**
2, 
3012,1,1.,3002,1,-1.
2, 
3102,1,1.,3002,1,-1.
2, 
3112,1,1.,3002,1,-1.
**
2, 
4012,1,1.,4002,1,-1.
2, 
4102,1,1.,4002,1,-1.
2, 
4112,1,1.,4002,1,-1.
2, 
4015,1,1.,4002,1,-1.
2, 
4016,1,1.,4002,1,-1.
2, 
4017,1,1.,4002,1,-1.
2, 
4018,1,1.,4002,1,-1.
2, 
4019,1,1.,4002,1,-1.
**
2, 
5012,1,1.,5002,1,-1.
2, 
5102,1,1.,5002,1,-1.
2, 
5112,1,1.,5002,1,-1.
2, 
5015,1,1.,5002,1,-1.
2, 
5016,1,1.,5002,1,-1.
2, 
5017,1,1.,5002,1,-1.
2, 
5018,1,1.,5002,1,-1.
2, 
5019,1,1.,5002,1,-1.
**
2, 
6012,1,1.,6002,1,-1.
2, 
6102,1,1.,6002,1,-1.
2, 
6112,1,1.,6002,1,-1.
2, 
6050,1,1.,6002,1,-1.
2, 
6060,1,1.,6002,1,-1.
2, 
6070,1,1.,6002,1,-1.
2, 
6080,1,1.,6002,1,-1.
2, 
6601,1,1.,6002,1,-1
**
2, 
7012,1,1.,7002,1,-1.
2, 
7102,1,1.,7002,1,-1.
2, 
7112,1,1.,7002,1,-1.
2, 
7050,1,1.,7002,1,-1.
2, 
7060,1,1.,7002,1,-1.
2, 
7070,1,1.,7002,1,-1.
2, 
7080,1,1.,7002,1,-1.
2, 
7601,1,1.,7002,1,-1
**********************************
*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
**
1001,1,1
1011,1,1
1101,1,1
1111,1,1
**
2001,1,3
2002,2,3
2011,1,3
2012,2,3
2101,1,3
2102,2,3
2111,1,3
2112,2,3
**
3001,1,1
3011,1,1
3101,1,1
3111,1,1
**
4001,1,3
4002,2,3
4011,1,3
4012,2,3
4101,1,3
4102,2,3
4111,1,3
4112,2,3
**
4005,1,3
4006,1,3
4007,1,3
4008,1,3
4009,1,3
4015,2,3
4016,2,3
4017,2,3
4018,2,3
4019,2,3
**
5001,1,1
5011,1,1
5101,1,1
5111,1,1
5005,1,1
5006,1,1
5007,1,1
5008,1,1
5009,1,1
**
6001,1,3
6002,2,3
6011,1,3
6012,2,3
6101,1,3
6102,2,3
6111,1,3
6112,2,3
**
6005,1,3
6006,1,3
6007,1,3
6008,1,3
6600,1,3
6050,2,3
6060,2,3
6070,2,3
6080,2,3
6601,2,3
**
7001,1,1
7011,1,1
7101,1,1
7111,1,1
7005,1,1
7006,1,1
7007,1,1
7008,1,1
7600,1,1
**********************************
*STEP
*FREQUENCY, eigensolver=lanczos, SIM=NO
8, , ,10. 
*OUTPUT,FIELD,FREQ=0
*MODAL FILE
*END STEP
****************************************
*STEP
*MODAL DYNAMIC
1e-3,1e0,
*MODAL DAMPING
** skip low-frequency mode
1,8
*SELECT EIGENMODES, GENERATE
1,8,1
*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
**
2002,1,-2.5
2012,1,-2.5
2102,1,-2.5
2112,1,-2.5
**
3002,1,-2.5
3012,1,-2.5
3102,1,-2.5
3112,1,-2.5
**
4002,1,-1.111111111111
4012,1,-1.111111111111
4102,1,-1.111111111111
4112,1,-1.111111111111
4015,1,-1.111111111111
4016,1,-1.111111111111
4017,1,-1.111111111111
4018,1,-1.111111111111
4019,1,-1.111111111111
**
5002,1,-1.111111111111
5012,1,-1.111111111111
5102,1,-1.111111111111
5112,1,-1.111111111111
5015,1,-1.111111111111
5016,1,-1.111111111111
5017,1,-1.111111111111
5018,1,-1.111111111111
5019,1,-1.111111111111
**
6002,1,-1.111111111111
6012,1,-1.111111111111
6102,1,-1.111111111111
6112,1,-1.111111111111
6050,1,-1.111111111111
6060,1,-1.111111111111
6070,1,-1.111111111111
6080,1,-1.111111111111
6601,1,-1.111111111111
**
7002,1,-1.111111111111
7012,1,-1.111111111111
7102,1,-1.111111111111
7112,1,-1.111111111111
7050,1,-1.111111111111
7060,1,-1.111111111111
7070,1,-1.111111111111
7080,1,-1.111111111111
7601,1,-1.111111111111
*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