*HEADING
 CASK DROP PROBLEM
 AXISYMMETRIC MODEL (refined mesh)
 SECTION CONTROLS USED (HOURGLASS=COMBINED)
**
** Create liner mesh
**
*NODE,NSET=NSHE1
1, 0.0,28.3461
2, 0.59055,28.3461
3, 1.1811,28.3461
4, 1.7716,28.3461
5, 2.3622,28.3461
6, 2.9527,28.3461
7, 3.5433,28.3461
8, 4.1338,28.3461
9, 4.7244,28.3461
10, 5.7368,28.3461
11, 6.7491,28.3461
12, 7.7615,28.3461
13, 8.7738,28.3461
14, 9.7862,28.3461
15, 10.799,28.3461
16, 11.8109,28.3461
17, 11.8109,27.308
18, 11.8109,26.270
19, 11.8109,25.232
20, 11.8109,24.194
21, 11.8109,23.156
22, 11.8109,22.119
23, 11.8109,21.081
24, 11.8109,20.043
25, 11.8109,19.005
26, 11.8109,17.967
27, 11.8109,16.929
28, 11.8109,15.891
29, 11.8109,14.853
30, 11.8109,13.815
31, 11.8109,12.777
32, 11.8109,11.739
33, 11.8109,10.701
34, 11.8109,9.6634
35, 11.8109,8.6255
36, 11.8109,7.5876
37, 11.8109,6.5496
38, 11.8109,5.5117
39, 10.799,5.5117
40, 9.7862,5.5117
41, 8.7738,5.5117
42, 7.7615,5.5117
43, 6.7491,5.5117
44, 5.7368,5.5117
45, 4.7244,5.5117
46, 4.1338,5.5117
47, 3.5433,5.5117
48, 2.9527,5.5117
49, 2.3622,5.5117
50, 1.7716,5.5117
51, 1.1811,5.5117
52, 0.59055,5.5117
53, 0.0,5.5117
54,11.8109,4.4094
55,11.8109,3.3070
56,11.8109,2.2047
57,11.8109,1.1023
58,11.8109,0.0
59, 10.799,0.0
60, 9.7862,0.0
61, 8.7738,0.0
62, 7.7615,0.0
63, 6.7491,0.0
64, 5.7368,0.0
65, 4.7244,0.0
66, 4.1338,0.0
67, 3.5433,0.0
68, 2.9527,0.0
69, 2.3622,0.0
70, 1.7716,0.0
71, 1.1811,0.0
72, 0.59055,0.0
73, 0.0,0.0
*ELEMENT, TYPE=SAX1
1, 1,2
16, 16,17
38, 38,39
53, 38,54
54, 54,55
*ELGEN, ELSET=S1
1,15,1,1
*ELGEN, ELSET=S2
16,22,1,1
*ELGEN, ELSET=S3
38,15,1,1
*ELSET, ELSET=S4
53,
*ELGEN, ELSET=S4
54,19,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=N1
201, 0.0,26.0627
202, 0.59055,26.0627
203, 1.1811,26.0627
204, 1.7716,26.0627
205, 2.3622,26.0627
206, 2.9527,26.0627
207, 3.5433,26.0627
208, 4.1338,26.0627
209, 4.7244,26.0627
210, 5.7368,26.0627
211, 6.7491,26.0627
212, 7.7615,26.0627
213, 8.7738,26.0627
214, 9.7862,26.0627
215, 10.799,26.0627
216, 11.8109,26.0627
*NODE,NSET=N2
585, 0.0,5.5117
586, 0.59055,5.5117
587, 1.1811,5.5117
588, 1.7716,5.5117
589, 2.3622,5.5117
590, 2.9527,5.5117
591, 3.5433,5.5117
592, 4.1338,5.5117
593, 4.7244,5.5117
594, 5.7368,5.5117
595, 6.7491,5.5117
596, 7.7615,5.5117
597, 8.7738,5.5117
598, 9.7862,5.5117
599, 10.799,5.5117
600, 11.8109,5.5117
*NFILL,NSET=NWATER
N1,N2,24,16
*ELEMENT, TYPE=CAX4R
201, 202,201,217,218
*ELGEN, ELSET=WATER
201, 15,1,1,24,16,15
*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,CONTROL=B
*SECTION CONTROLS, HOURGLASS=COMBINED, NAME=B
**
** Create foam mesh
**
*NODE,NSET=N3
601, 0.0,5.0117
602, 0.59055,5.0117
603, 1.1811,5.0117
604, 1.7716,5.0117
605, 2.3622,5.0117
606, 2.9527,5.0117
607, 3.5433,5.0117
608, 4.1338,5.0117
609, 4.7244,5.0117
610, 5.7368,5.0117
611, 6.7491,5.0117
612, 7.7615,5.0117
613, 8.7738,5.0117
614, 9.7862,5.0117
615, 10.799,5.0117
616, 11.8109,5.0117
*NODE,NSET=N4
681, 0.0,0.0
682, 0.59055,0.0
683, 1.1811,0.0
684, 1.7716,0.0
685, 2.3622,0.0
686, 2.9527,0.0
687, 3.5433,0.0
688, 4.1338,0.0
689, 4.7244,0.0
690, 5.7368,0.0
691, 6.7491,0.0
692, 7.7615,0.0
693, 8.7738,0.0
694, 9.7862,0.0
695, 10.799,0.0
696, 11.8109,0.0
*NFILL, NSET=NFOAM
N3,N4,5,16
*ELEMENT, TYPE=CAX4R
601, 602,601,617,618
*ELGEN
601, 15,1,1,5,16,15
*ELSET, ELSET=FOAM, GENERATE
601,675,1
**
*SOLID SECTION,ELSET=FOAM,MATERIAL=MFOAM,CONTROL=B
*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
201,202,203,204,205,206,207,208,209,210,211,212,213,214,215,216,
585,586,587,588,589,590,591,592,593,594,595,596,597,598,599,600
*NSET, NSET=WCON1, GENERATE
216,600,16
*NSET,NSET=FCON1, GENERATE
601,616,1
616,696,16
681,696,1
**
** Define boundary condition
**
*NSET, NSET=FEND, GENERATE
58,73,1
*NSET, NSET=WAXI, GENERATE
201,585,16
*NSET, NSET=FAXI, GENERATE
601,681,16
*BOUNDARY
FEND,2,2
1,XSYMM
53,XSYMM
73,XSYMM
WAXI,XSYMM
FAXI,XSYMM
**
** Define initial condition
**
*NSET,NSET=NALL, GENERATE
1,57,1
201,696,1
*INITIAL CONDITIONS, TYPE=VELOCITY
NALL, 2, -523.3
**
** Define history output element set
**
*ELSET,ELSET=FLUID
546, 
*ELSET,ELSET=SOLID
601, 
*ELSET,ELSET=MIDEL
546,601
*ELSET,ELSET=SHELL
55, 
*ELSET,ELSET=ELOUT
MIDEL, 
SHELL, 
*NSET, NSET=NOUT
55,601
*NSET, NSET=NOUT, ELSET=ELOUT
*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,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