*HEADING
 CASK DROP PROBLEM
 AXISYMMETRIC MODEL (coarse mesh)
**
** Create liner mesh
**
*NODE,NSET=NSHE1
1, 0.0,28.3461
2, 2.3622,28.3461
3, 4.7244,28.3461
4, 7.0865,28.3461
5, 9.4487,28.3461
6, 11.8109,28.3461
7, 11.8109,27.1650
8, 11.8109,25.1965
9, 11.8109,23.2281
10, 11.8109,21.2596
11, 11.8109,19.2911
12, 11.8109,17.3226
13, 11.8109,15.3541
14, 11.8109,13.3857
15, 11.8109,11.4172
16,11.8109,9.4487
17,11.8109,8.4645
18,11.8109,7.4802
19,11.8109,6.4960
20,11.8109,5.5117
21,9.4487,5.5117
22,7.0865,5.5117
23,4.7244,5.5117
24,2.3622,5.5117
25, 0.0,5.5117
26,11.8109,4.4094
27,11.8109,3.3070
28,11.8109,2.2047
29,11.8109,1.1023
30,11.8109,0.0
31,9.4487,0.0
32,7.0865,0.0
33,4.7244,0.0
34,2.3622,0.0
35,0.0,0.0
*ELEMENT, TYPE=SAX1
1, 1,2
6, 6,7
20, 20,26
21, 26,27
30, 20,21
*ELGEN, ELSET=S1
1,5,1,1
*ELGEN, ELSET=S2
6,14,1,1
*ELGEN, ELSET=S3
30, 5,1,1
*ELSET,ELSET=S4
20,
*ELGEN, ELSET=S4
21,9,1,1
*ELSET,ELSET=SSTEEL
S1,S2,S3,
*ELSET,ELSET=MSTEEL
S4,
**
** Shell material definition
**
*SHELL SECTION,MATERIAL=A2,ELSET=S1, SECTION INTEGRATION=GAUSS
0.3748,5
*SHELL SECTION,MATERIAL=A1,ELSET=S2, SECTION INTEGRATION=GAUSS
0.1874,5
*SHELL SECTION,MATERIAL=A2,ELSET=S3, SECTION INTEGRATION=GAUSS
0.5,5
*SHELL SECTION,MATERIAL=A2,ELSET=S4, SECTION INTEGRATION=GAUSS
0.0591,5
**
** Stainess steel (304L)
**
*MATERIAL,NAME=A1
*DENSITY
0.00075163, 
*ELASTIC
28E6,0.28
*PLASTIC
44300.0,0.
66300.0,0.1
**
** Mild steel (A36)
**
*MATERIAL,NAME=A2
*DENSITY
0.00075163, 
*ELASTIC
28E6,0.28
*PLASTIC
30000.0,0.
**
** Create fluid mesh
**
*NODE, NSET=WTOP
101, 0.0,26.0627
102, 2.36218,26.0627
103, 4.72436,26.0627
104, 7.08654,26.0627
105, 9.44872,26.0627
106, 11.8109,26.0627
*NODE, NSET=WMID
155, 0.0,9.2482
156, 2.36218,9.2482
157, 4.72436,9.2482
158, 7.08654,9.2482
159, 9.44872,9.2482
160, 11.8109,9.2482
*NODE, NSET=WBOT
179, 0.0,5.5117
180, 2.36218,5.5117
181, 4.72436,5.5117
182, 7.08654,5.5117
183, 9.44872,5.5117
184, 11.8109,5.5117
*NFILL,NSET=NWATER
WTOP,WMID,9,6
*NFILL,NSET=NWATER
WMID,WBOT,4,6
*ELEMENT, TYPE=CAX4R
101, 102,101,107,108
*ELGEN, ELSET=WATER
101, 5,1,1,13,6,5
*MATERIAL, NAME=MWATER
*DENSITY
9.2E-5, 
*EOS, TYPE=USUP
57100.874,0.0,0.0
*TENSILE FAILURE,ELEMENT DELETION=NO,PRESSURE=BRITTLE,SHEAR=BRITTLE
 0.,
*SOLID SECTION, ELSET=WATER, MATERIAL=MWATER
**
** Create foam mesh
**
*NODE, NSET=FTOP
201, 0.0,5.0117
202, 2.36218,5.0117
203, 4.72436,5.0117
204, 7.08654,5.0117
205, 9.44872,5.0117
206, 11.8109,5.0117
*NODE,NSET=FBOT
231, 0.0,0.0
232, 2.36218,0.0
233, 4.72436,0.0
234, 7.08654,0.0
235, 9.44872,0.0
236, 11.8109,0.0
*NFILL, NSET=NFOAM
FTOP,FBOT,5,6
*ELEMENT, TYPE=CAX4R
201, 202,201,207,208
*ELGEN, ELSET=FOAM
201, 5,1,1,5,6,5
**
*SOLID SECTION,ELSET=FOAM,MATERIAL=MFOAM
*MATERIAL,NAME=MFOAM
*ELASTIC
 18738., 0.0
*CRUSHABLE FOAM,HARDENING=VOLUMETRIC
 0.3916, 0.2798
** k, kt 
** k = sigma_0 / p_c0 = 313.3 / 800.0
** kt = p_t / p_c0 = 223.8 / 800.0 = 0.2798
**
*CRUSHABLE FOAM HARDENING
    313.3,  0.000
    325.3,  0.010
    337.4,  0.020
    468.6,  0.345
    712.8,  0.440
   1188.9,  0.490
   2128.3,  0.510
 110067.7,  2.000
** yield stress is in psi unit
*DENSITY
 2.85E-5,
**
** Define fluid and foam contact node set
**
*NSET, NSET=WCON1
101,102,103,104,105,106,112,118,124,130,136,142,148,154,160,166,
172,178,179,180,181,182,183,184
*NSET,NSET=FCON1
201,202,203,204,205,206,212,218,224,230,231,232,233,234,235,236
**
** Define boundary condition
**
*NSET, NSET=FEND
30,31,32,33,34,35
*NSET, NSET=WAXI, GENERATE
101,179,6
*NSET, NSET=FAXI, GENERATE
201,231,6
*BOUNDARY
FEND,2,2
35,XSYMM
25,XSYMM
1,XSYMM
WAXI,XSYMM
FAXI,XSYMM
**
** Define initial condition
**
*NSET,NSET=NALL, GENERATE
1,29,1
101,184,1
201,236,1
*INITIAL CONDITIONS, TYPE=VELOCITY
NALL, 2, -523.3
**
** Define history output element set
**
*ELSET,ELSET=FLUID
161, 
*ELSET,ELSET=SOLID
201, 
*ELSET,ELSET=MIDEL
161,201
*ELSET,ELSET=SHELL
22, 
*ELSET,ELSET=ELOUT
MIDEL, 
SHELL, 
*NSET, NSET=NOUT
27,201
*RESTART, WRITE, NUM=1, TIMEMARKS=NO
*SURFACE,TYPE=NODE,NAME=FCON1
 FCON1,
*SURFACE,TYPE=NODE,NAME=WCON1
 WCON1,
*SURFACE,TYPE=ELEMENT, NAME=SSIDE1, NO THICK
S1,SNEG
S2,SNEG
S3,SNEG
*SURFACE,TYPE=ELEMENT, NAME=SSIDE2, NO THICK
S3,SPOS
S4,SNEG
*STEP
*DYNAMIC, EXPLICIT
,0.005
*CONTACT PAIR
SSIDE1,WCON1
*CONTACT PAIR
 SSIDE2,FCON1
*FILE OUTPUT, NUM=4, TIMEMARKS=NO
*EL FILE, ELSET=FLUID
 PRESS,
*EL FILE, ELSET=SOLID
 PEEQ,
*EL FILE, ELSET=SHELL
 PEEQ,
 STH,
*NODE FILE, NSET=NOUT
 U,
**********
*ELSET, ELSET=QA_TEST_FLUID
FLUID,
*ELSET, ELSET=QA_TEST
SOLID, SHELL
*NSET, NSET=QA_TEST
FEND,
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=4
*ELEMENT OUTPUT, ELSET=QA_TEST_FLUID

PRESS,
*ELEMENT OUTPUT, ELSET=QA_TEST

PEEQ, STH
*NODE OUTPUT, NSET=QA_TEST
U,
*OUTPUT, HIST, FREQ=999
*ENERGY OUTPUT, VAR=PRESELECT 
***********
***OUTPUT,HISTORY,TIME INTERVAL=0.0
***ENERGY OUTPUT
**ALLAE,ALLIE,ALLKE,ALLPD,ALLSE,ALLVD,ALLWK,ETOTAL
***NODE OUTPUT,NSET=NOUT
**U,V
***ELEMENT OUTPUT,ELSET=MIDEL
**S22,PEEQ,MISES,ERV
***ELEMENT OUTPUT,ELSET=ELOUT
**PRESS,
***ELEMENT OUTPUT,ELSET=SHELL
**SF,STH
***ELEMENT OUTPUT,ELSET=SHELL,SECTIONPOINT=1
**S22,PEEQ,MISES
***ELEMENT OUTPUT,ELSET=SHELL,SECTIONPOINT=2
**S22,PEEQ,MISES
***ELEMENT OUTPUT,ELSET=SHELL,SECTIONPOINT=3
**S22,PEEQ,MISES
**
***OUTPUT,FIELD,VARIABLE=PRESELECT,OP=NEW,NUM=1,TIMEMARKS=NO
*OUTPUT,FIELD,OP=NEW,NUM=4,TIMEMARKS=NO
*ELEMENT OUTPUT,ELSET=FLUID
 PRESS,
*ELEMENT OUTPUT,ELSET=SOLID
 PEEQ,PE,LE
 MISES,PRESS
*ELEMENT OUTPUT,ELSET=SHELL
 PEEQ,PE,LE
 MISES,PRESS
 SF,STH
*NODE OUTPUT
 U,V
*OUTPUT,HISTORY,OP=NEW,TIME INTERVAL=0.0
*ENERGY OUTPUT
 ALLAE,ALLIE,ALLKE,ALLPD,ALLSE,ALLVD,ALLWK,ETOTAL
*NODE OUTPUT,NSET=NOUT
 U,V
*ELEMENT OUTPUT,ELSET=MIDEL
 S22,PEEQ,MISES,ERV
*ELEMENT OUTPUT,ELSET=ELOUT
 PRESS,
*ELEMENT OUTPUT,ELSET=SHELL
 SF,STH
*ELEMENT OUTPUT,ELSET=SHELL
 S22,PEEQ,MISES
*ELEMENT OUTPUT,ELSET=SHELL
 S22,PEEQ,MISES
*ELEMENT OUTPUT,ELSET=SHELL
 S22,PEEQ,MISES
*END STEP