*HEADING
 CASK DROP PROBLEM
 3D MODEL (coarse mesh)
 SECTION CONTROLS USED (KINEMA=ORTHOGONAL,
 HOURGLASS=COMBINED)
**
** Create liner mesh
**
** Top shell mesh
*PREPRINT, CONTACT=YES, MODEL=YES, HISTORY=YES
*NODE,NSET=N1
1, 0.0,0.0,28.3461
2, 2.3622,0.0,28.3461
3, 0.0,2.3622,28.3461
4, 2.3622,2.3622,28.3461
*NODE,NSET=N2
5, 4.7244,0.0,28.3461
6, 4.7244,2.3622,28.3461
7, 4.7244,4.7244,28.3461
8, 2.3622,4.7244,28.3461
9, 0.0,4.7244,28.3461
*NODE,NSET=N3
20, 11.8109,0.0,28.3461
21,10.911849,4.5198358,28.3461
22,8.3515675,8.3515675,28.3461
23,4.5198358,10.911849,28.3461
24, 0.0,11.8109,28.3461
*NFILL,NSET=N31
N2,N3,3,5
*ELEMENT,TYPE=S4R,ELSET=ETSHELL
1, 1,2,4,3
2, 2,5,6,4
3, 3,4,8,9
4, 4,6,7,8
5, 5,10,11,6
*NSET,NSET=NTSHELL
N1,N2,N3,N31
*ELGEN, ELSET=ETSHELL
5,4,1,1,3,5,4
** First shell wall and middle shell mesh
*NCOPY,OLD SET=N3, NEW SET=N4, CHANGE NUMBER=5,
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=N4, NEW SET=N5, CHANGE NUMBER=45,
SHIFT
0.0,0.0,-17.7163
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
*NCOPY,OLD SET=NTSHELL, NEW SET=NMSHELL,
CHANGE NUMBER=89, 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
N4,N5,9,5
*NCOPY,OLD SET=N4, NEW SET=N6, CHANGE NUMBER=60,
SHIFT
0.0,0.0,-20.669
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
*NFILL
N5,N6,3,5
** Second shell wall and bottom shell mesh
*NCOPY,OLD SET=NTSHELL, NEW SET=NBSHELL,
CHANGE NUMBER=133, 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
*NCOPY,OLD SET=N4, NEW SET=N7, CHANGE NUMBER=104,
SHIFT
0.0,0.0,-26.0627
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
*NSET,NSET=N61
109,110,111,112,113
*NFILL
N61,N7,4,5
*ELEMENT,TYPE=S4R
17, 25,26,21,20
69, 109,110,86,85
89, 114,115,110,109
105,153,154,130,129
*ELGEN, ELSET=W1SHELL
17, 4,1,1,13,5,4
*ELGEN, ELSET=W1SHELL
69, 4,1,1
*ELCOPY, ELEMENT SHIFT=72, OLD SET=ETSHELL,
SHIFT NODE=89, NEW SET=EMSHELL
*ELCOPY, ELEMENT SHIFT=36, OLD SET=EMSHELL,
SHIFT NODE=44, NEW SET=EBSHELL
*ELGEN,ELSET=W2SHELL
89,4,1,1,4,5,4
*ELGEN,ELSET=W2SHELL
105, 4,1,1
*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
**
** Stainless 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
**
*NCOPY,OLD SET=NTSHELL,CHANGE NUMBER=200,
NEW SET=WTOP,SHIFT
0.0,0.0,-2.2834
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
*NCOPY,OLD SET=WTOP, NEW SET=WMID,
CHANGE NUMBER=216, SHIFT
0.0,0.0,-16.8145
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
*NCOPY,OLD SET=WTOP, NEW SET=WBOT,
CHANGE NUMBER=312, SHIFT
0.0,0.0,-20.551
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
*NFILL
WTOP,WMID,9,24
*NFILL
WMID,WBOT,4,24
*ELEMENT,TYPE=C3D8R,ELSET=E1
201, 225,226,228,227,201,202,204,203
202, 226,229,230,228,202,205,206,204
203, 227,228,232,233,203,204,208,209
204, 228,230,231,232,204,206,207,208
205, 229,234,235,230,205,210,211,206
206, 230,235,236,231,206,211,212,207
207, 231,236,237,232,207,212,213,208
208, 232,237,238,233,208,213,214,209
209, 234,239,240,235,210,215,216,211
210, 235,240,241,236,211,216,217,212
211, 236,241,242,237,212,217,218,213
212, 237,242,243,238,213,218,219,214
213, 239,244,245,240,215,220,221,216
214, 240,245,246,241,216,221,222,217
215, 241,246,247,242,217,222,223,218
216, 242,247,248,243,218,223,224,219
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=16,
SHIFT NODE=24
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=32,
SHIFT NODE=48
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=48,
SHIFT NODE=72
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=64,
SHIFT NODE=96
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=80,
SHIFT NODE=120
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=96,
SHIFT NODE=144
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=112,
SHIFT NODE=168
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=128,
SHIFT NODE=192
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=144,
SHIFT NODE=216
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=160,
SHIFT NODE=240
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=176,
SHIFT NODE=264
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=192,
SHIFT NODE=288
*ELSET,ELSET=WATER,GENERATE
201,408,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=NMSHELL,CHANGE NUMBER=511,
NEW SET=FTOP,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=NBSHELL,CHANGE NUMBER=587,
NEW SET=FBOT,SHIFT
0.0,0.0,0.0
0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
*NFILL
FTOP,FBOT,5,24
*ELCOPY,OLD SET=E1,ELEMENT SHIFT=300,
SHIFT NODE=400,NEW SET=E2
*ELCOPY,OLD SET=E2,ELEMENT SHIFT=16,
SHIFT NODE=24
*ELCOPY,OLD SET=E2,ELEMENT SHIFT=32,
SHIFT NODE=48
*ELCOPY,OLD SET=E2,ELEMENT SHIFT=48,
SHIFT NODE=72
*ELCOPY,OLD SET=E2,ELEMENT SHIFT=64,
SHIFT NODE=96
*ELSET,ELSET=FOAM,GENERATE
501,580,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
201,224,1
220,532,24
221,533,24
222,534,24
223,535,24
224,536,24
513,536,1
*NSET, NSET=FCONT, GENERATE
601,624,1
620,740,24
621,741,24
622,742,24
623,743,24
624,744,24
721,744,1
**
** Define boundary condition
**
*NSET,NSET=FEND, GENERATE
134,157,1
*NSET,NSET=YANODES
1,2,5,10,15,20,25,30,35,40,45,50,55,60,65,70,
75,80,85,90,91,94,99,104,109,114,119,124,129,
134,135,138,143,148,153,513,514,517,522,527,
532,721,722,725,730,735,740,
*NSET,NSET=YANODES,GENERATE
201,513,24
202,514,24
205,517,24
210,522,24
215,527,24
220,532,24
601,721,24
602,722,24
605,725,24
610,730,24
615,735,24
620,740,24
*NSET,NSET=XANODES
1,3,9,14,19,24,29,34,39,44,49,54,59,64,69,74,
79,84,89,90,92,98,103,108,113,118,123,128,133,
134,136,142,147,152,157,513,515,521,526,531,
536,721,723,729,734,739,744,
*NSET,NSET=XANODES,GENERATE
201,513,24
203,515,24
209,521,24
214,526,24
219,531,24
224,536,24
601,721,24
603,723,24
609,729,24
614,734,24
619,739,24
624,744,24
*BOUNDARY
FEND,ZSYMM
XANODES,XSYMM
YANODES,YSYMM
**
** Define initial condition
**
*NSET,NSET=NALL,GENERATE
1,133,1
201,536,1
601,720,1
*INITIAL CONDITIONS, TYPE=VELOCITY
NALL, 3, -523.3
**
** Define history output element set
**
*ELSET,ELSET=FLUID
393, 
*ELSET,ELSET=SOLID
501, 
*ELSET,ELSET=MIDEL
393,501
*ELSET,ELSET=SHELL
97, 
*ELSET,ELSET=ELOUT
MIDEL, 
SHELL, 
*NSET, NSET=NOUT
119,601
*NSET, NSET=NOUT, ELSET=ELOUT
*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