*HEADING
 CASK DROP PROBLEM
 3D MODEL (refined mesh)
 SECTION CONTROLS USED (KINEMA=ORTHOGONAL, HOURGLASS=COMBINED)
**
** Create liner mesh
**
** Top shell mesh
*PREPRINT, CONTACT=YES, MODEL=YES, HISTORY=YES
*NODE,NSET=N11
1, 0.0,0.0,28.3461
*NODE,NSET=N12
8, 0.0,4.13385,28.3461
*NODE,NSET=N13
57, 4.13385,0.0,28.3461
*NODE,NSET=N14
64, 4.13385,4.13385,28.3461
*NODE,NSET=N15
65, 4.7244,0.0,28.3461
*NODE,NSET=N16
73, 4.7244,4.7244,28.3461
*NODE,NSET=N17
81, 0.0,4.7244,28.3461
*NODE,NSET=N4
184, 11.8109,0.0,28.3461
185, 11.754027,1.1576706,28.3461
186, 11.583957,2.3041923,28.3461
187, 11.302327,3.4285233,28.3461
188,10.911849,4.5198358,28.3461
189,10.416284,5.5676197,28.3461
190,9.8204044,6.5617845,28.3461
191,9.1299492,7.4927556,28.3461
192,8.3515675,8.3515675,28.3461
193,7.4927556,9.1299492,28.3461
194,6.5617845,9.8204044,28.3461
195,5.5676197,10.416284,28.3461
196,4.5198358,10.911849,28.3461
197,3.4185233,11.302327,28.3461
198,2.3041923,11.583957,28.3461
199,1.1576706,11.754027,28.3461
200, 0.0,11.8109,28.3461
*NFILL,NSET=N1
N11,N12,7,1
*NFILL,NSET=N2
N13,N14,7,1
*NFILL
N1,N2,7,8
*NFILL,NSET=N3
N15,N16,8,1
*NFILL,NSET=N3
N16,N17,8,1
*NFILL
N3,N4,7,17
*ELEMENT,TYPE=S4R
1, 1,9,10,2
50, 57,65,66,58
57, 64,72,73,74
58, 56,64,74,75
59, 48,56,75,76
60, 40,48,76,77
61, 32,40,77,78
62, 24,32,78,79
63, 16,24,79,80
64,  8,16,80,81
65, 65,82,83,66
*ELGEN
1,7,1,1,7,8,7
*ELGEN
50,7,1,1
*ELGEN
65,16,1,1,7,17,16
*NSET,NSET=NTOP,GENERATE
1,200,1
*ELSET,ELSET=ETSHELL,GENERATE
1,176,1
**
**
**
*NCOPY,OLD SET=N4, NEW SET=N5, CHANGE NUMBER=357, SHIFT
0.0,0.0,-21.796473
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
*NFILL
N4,N5,21,17
*ELEMENT,TYPE=S4R
177,201,202,185,184
513,741,742,542,541
*ELGEN
177,16,1,1,21,17,16
513,16,1,1
*ELSET,ELSET=W1SHELL,GENERATE
177,528,1
**
**
**
*NCOPY,OLD SET=NTOP, NEW SET=NMID, CHANGE NUMBER=557, SHIFT
0.0,0.0,-22.8344
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
*ELCOPY, ELEMENT SHIFT=528, OLD SET=ETSHELL, SHIFT NODE=557, NEW SET=EMSHELL
**
**
**
*NSET,NSET=N6,GENERATE
741,757,1
*NCOPY,OLD SET=N4, NEW SET=N7, CHANGE NUMBER=625, SHIFT
0.0,0.0,-27.24376
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
*NFILL
N6,N7,4,17
*ELEMENT,TYPE=S4R
705, 758,759,742,741
769, 1009,1010,810,809
*ELGEN
705, 16,1,1,4,17,16
769, 16,1,1
*ELSET,ELSET=W2SHELL,GENERATE
705,784,1
**
**
**
*NCOPY,OLD SET=NTOP, NEW SET=NBOT, CHANGE NUMBER=825, SHIFT
0.0,0.0,-28.3461
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
*ELCOPY, ELEMENT SHIFT=784, OLD SET=ETSHELL, SHIFT NODE=825, NEW SET=EBSHELL
*ELSET,ELSET=SHELL1
ETSHELL,W1SHELL,EMSHELL,W2SHELL,EBSHELL
**
** Shell material definition
**
*SHELL SECTION,MATERIAL=A2,ELSET=ETSHELL, SECTION INTEGRATION=GAUSS
0.3748,5
*SHELL SECTION,MATERIAL=A1,ELSET=W1SHELL, SECTION INTEGRATION=GAUSS
0.1874,5
*SHELL SECTION,MATERIAL=A2,ELSET=EMSHELL, SECTION INTEGRATION=GAUSS
0.5,5
*SHELL SECTION,MATERIAL=A2,ELSET=W2SHELL, SECTION INTEGRATION=GAUSS
0.0591,5
*SHELL SECTION,MATERIAL=A2,ELSET=EBSHELL, 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
**
** Top water nodes
*NCOPY,OLD SET=NTOP, NEW SET=NWTOP, CHANGE NUMBER=2000, SHIFT
0.0,0.0,-1.1811
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
*NCOPY,OLD SET=NTOP, NEW SET=NWBOT, CHANGE NUMBER=6800, SHIFT
0.0,0.0,-22.8344
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
*NFILL,NSET=NWATER
NWTOP,NWBOT,24,200
*ELEMENT, TYPE=C3D8R
2001,2201,2209,2210,2202,2001,2009,2010,2002
2050,2257,2265,2266,2258,2057,2065,2066,2058
2057,2264,2272,2273,2274,2064,2072,2073,2074
2058,2256,2264,2274,2275,2056,2064,2074,2075
2059,2248,2256,2275,2276,2048,2056,2075,2076
2060,2240,2248,2276,2277,2040,2048,2076,2077
2061,2232,2240,2277,2278,2032,2040,2077,2078
2062,2224,2232,2278,2279,2024,2032,2078,2079
2063,2216,2224,2279,2280,2016,2024,2079,2080
2064,2208,2216,2280,2281,2008,2016,2080,2081
2065,2265,2282,2283,2266,2065,2082,2083,2066
*ELGEN
2001,7,1,1,7,8,7
2050,7,1,1
2065,16,1,1,7,17,16
*ELSET,ELSET=E1,GENERATE
2001,2176,1
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=176,SHIFT NODE=200
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=352,SHIFT NODE=400
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=528,SHIFT NODE=600
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=704,SHIFT NODE=800
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=880,SHIFT NODE=1000
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=1056,SHIFT NODE=1200
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=1232,SHIFT NODE=1400
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=1408,SHIFT NODE=1600
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=1584,SHIFT NODE=1800
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=1760,SHIFT NODE=2000
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=1936,SHIFT NODE=2200
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=2112,SHIFT NODE=2400
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=2288,SHIFT NODE=2600
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=2464,SHIFT NODE=2800
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=2640,SHIFT NODE=3000
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=2816,SHIFT NODE=3200
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=2992,SHIFT NODE=3400
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=3168,SHIFT NODE=3600
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=3344,SHIFT NODE=3800
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=3520,SHIFT NODE=4000
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=3696,SHIFT NODE=4200
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=3872,SHIFT NODE=4400
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=4048,SHIFT NODE=4600
*ELSET,ELSET=WATER,GENERATE
2001,6224,1
*NSET,NSET=NWATER,GENERATE
2001,7000,1
*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=C
*SECTION CONTROLS, KINEMA=ORTHOGONAL, HOURGLASS=COMBINED, NAME=C
**
** Create foam mesh
**
*NCOPY,OLD SET=NWBOT,CHANGE NUMBER=1200,NEW SET=NFTOP,SHIFT
0.0,0.0,-0.5
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
*NCOPY,OLD SET=NFTOP,CHANGE NUMBER=1000,NEW SET=NFBOT,SHIFT
0.0,0.0,-5.0117
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
*NFILL
NFTOP,NFBOT,5,200
*ELSET,ELSET=E2,GENERATE
2001,2880,1
*ELCOPY,OLD SET=E2,ELEMENT SHIFT=6000,SHIFT NODE=6000,NEW SET=EFOAM
*ELSET,ELSET=FOAM,GENERATE
 8001,8880,1
**
*SOLID SECTION,ELSET=FOAM,MATERIAL=MFOAM,CONTROL=C
*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=WCONT, GENERATE
2001,2183,1
2184,6984,200
2185,6985,200
2186,6986,200
2187,6987,200
2188,6988,200
2189,6989,200
2190,6990,200
2191,6991,200
2192,6992,200
2193,6993,200
2194,6994,200
2195,6995,200
2196,6996,200
2197,6997,200
2198,6998,200
2199,6999,200
2200,7000,200
6801,6983,1
*NSET, NSET=FCONT, GENERATE
8001,8183,1
8184,9184,200
8185,9185,200
8186,9186,200
8187,9187,200
8188,9188,200
8189,9189,200
8190,9190,200
8191,9191,200
8192,9192,200
8193,9193,200
8194,9194,200
8195,9195,200
8196,9196,200
8197,9197,200
8198,9198,200
8199,9199,200
8200,9200,200
9001,9183,1
**
** Define boundary condition
**
*NSET,NSET=YANODES
1,9,17,25,33,41,49,57,65,82,99,116,133,150,167,184,
558,566,574,582,590,598,606,614,622,639,656,673,690,707,724,741,
826,834,842,850,858,866,874,882,890,907,924,941,958,975,992,1009,
*NSET,NSET=YANODES,GENERATE
184,541,17
741,809,17
2001,6801,200
2009,6809,200
2017,6817,200
2025,6825,200
2033,6833,200
2041,6841,200
2049,6849,200
2057,6857,200
2065,6865,200
2082,6882,200
2099,6899,200
2116,6916,200
2133,6933,200
2150,6950,200
2167,6967,200
2184,6984,200
8001,9001,200
8009,9009,200
8017,9017,200
8025,9025,200
8033,9033,200
8041,9041,200
8049,9049,200
8057,9057,200
8065,9065,200
8082,9082,200
8099,9099,200
8116,9116,200
8133,9133,200
8150,9150,200
8167,9167,200
8184,9184,200
*NSET,NSET=XANODES
1,2,3,4,5,6,7,8,81,98,115,132,149,166,183,200,
558,559,560,561,562,563,564,565,638,655,672,689,706,723,740,757,
826,827,828,829,830,831,832,833,906,923,940,957,974,991,1008,1025,
*NSET,NSET=XANODES,GENERATE
200,557,17
757,825,17
2001,6801,200
2002,6802,200
2003,6803,200
2004,6804,200
2005,6805,200
2006,6806,200
2007,6807,200
2008,6808,200
2081,6881,200
2098,6898,200
2115,6915,200
2132,6932,200
2149,6949,200
2166,6966,200
2183,6983,200
2200,7000,200
8001,9001,200
8002,9002,200
8003,9003,200
8004,9004,200
8005,9005,200
8006,9006,200
8007,9007,200
8008,9008,200
8081,9081,200
8098,9098,200
8115,9115,200
8132,9132,200
8149,9149,200
8166,9166,200
8183,9183,200
8200,9200,200
*NSET,NSET=FEND,GENERATE
826,1025,1
*BOUNDARY
FEND,ZSYMM
XANODES,XSYMM
YANODES,YSYMM
*NSET,NSET=NALL,GENERATE
1,825,1
2001,7000,1
8001,9200,1
*INITIAL CONDITIONS, TYPE=VELOCITY
NALL, 3, -523.3
*ELSET,ELSET=FLUID
6049, 
*ELSET,ELSET=SOLID
8001, 
*ELSET,ELSET=MIDEL
6049,8001
*ELSET,ELSET=SHELL
737, 
*ELSET,ELSET=ELOUT
MIDEL, 
SHELL, 
*NSET, NSET=NOUT
791,8001
*RESTART, WRITE, NUM=1, TIMEMARKS=NO
*SURFACE,TYPE=ELEMENT, NAME=SSIDE1
ETSHELL,SNEG
W1SHELL,SNEG
EMSHELL,SPOS
*SURFACE,TYPE=ELEMENT, NAME=SSIDE2
EMSHELL,SNEG
W2SHELL,SNEG
EBSHELL,SPOS
*SURFACE,TYPE=ELEMENT, NAME=SWATER
 WATER,
*SURFACE,TYPE=ELEMENT, NAME=SFOAM
 FOAM
*************************************************************
**General contact definition
*CONTACT
*SURFACEPROPERTYASSIGNMENT, PROPERTY=THICKNESS
 SSide1, 0.0, 
 SSide2, 0.0, 
*CONTACT INCLUSIONS, ALLELEMENTBASED
*CONTACT FORMULATION, TYPE=PURE MASTER-SLAVE
 SSIDE1, SWATER, MASTER
 SSIDE2, SFOAM,  MASTER
*************************************************************
*STEP
*DYNAMIC, EXPLICIT
,0.005
*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
STH
*NODE OUTPUT, NSET=QA_TEST
U,
*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,NSET=NOUT
 U,V
*OUTPUT,HISTORY,OP=NEW,TIME INTERVAL=0.005
*ENERGY OUTPUT
 ALLAE,ALLIE,ALLKE,ALLPD,ALLSE,ALLVD,ALLWK,ETOTAL
*OUTPUT,HISTORY,OP=NEW,TIME INTERVAL=0.0
*NODE OUTPUT,NSET=NOUT
 U,V
*ELEMENT OUTPUT,ELSET=MIDEL
 S33,PEEQ,MISES,ERV
*ELEMENT OUTPUT,ELSET=ELOUT
 PRESS,
*ELEMENT OUTPUT,ELSET=SHELL
 SF,STH
*ELEMENT OUTPUT,ELSET=SHELL
 S33,PEEQ,MISES
*ELEMENT OUTPUT,ELSET=SHELL
 S33,PEEQ,MISES
*ELEMENT OUTPUT,ELSET=SHELL
 S33,PEEQ,MISES
*END STEP