*HEADING
 CASK DROP PROBLEM
 3D MODEL (refined mesh)
**
** Create liner mesh
**
** Top shell mesh
*NODE,NSET=N11
*NODE,NSET=N12
*NODE,NSET=N13
*NODE,NSET=N14
*NODE,NSET=N15
*NODE,NSET=N16
*NODE,NSET=N17
*NODE,NSET=N4
*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
*ELGEN
*ELGEN
*ELGEN
*NSET,NSET=NTOP,GENERATE
*ELSET,ELSET=ETSHELL,GENERATE
**
**
**
*NCOPY,OLD SET=N4, NEW SET=N5, CHANGE NUMBER=357, SHIFT
*NFILL
N4,N5,21,17
*ELEMENT,TYPE=S4R
*ELGEN
*ELSET,ELSET=W1SHELL,GENERATE
**
**
**
*NCOPY,OLD SET=NTOP, NEW SET=NMID, CHANGE NUMBER=557, SHIFT
*ELCOPY, ELEMENT SHIFT=528, OLD SET=ETSHELL, SHIFT NODE=557, NEW SET=EMSHELL
**
**
**
*NSET,NSET=N6,GENERATE
*NCOPY,OLD SET=N4, NEW SET=N7, CHANGE NUMBER=625, SHIFT
*NFILL
N6,N7,4,17
*ELEMENT,TYPE=S4R
*ELGEN
*ELSET,ELSET=W2SHELL,GENERATE
**
**
**
*NCOPY,OLD SET=NTOP, NEW SET=NBOT, CHANGE NUMBER=825, SHIFT
*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
*NCOPY,OLD SET=NTOP, NEW SET=NWBOT, CHANGE NUMBER=6800, SHIFT
*NFILL,NSET=NWATER
NWTOP,NWBOT,24,200
*ELEMENT, TYPE=C3D8R
*ELGEN
*ELSET,ELSET=E1,GENERATE
*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
*NSET,NSET=NWATER,GENERATE
*MATERIAL, NAME=MWATER
*DENSITY
9.2E-5,
*EOS, TYPE=USUP
*TENSILE FAILURE,ELEMENT DELETION=NO,PRESSURE=BRITTLE,SHEAR=BRITTLE
*SOLID SECTION, ELSET=WATER, MATERIAL=MWATER
**
** Create foam mesh
**
*NCOPY,OLD SET=NWBOT,CHANGE NUMBER=1200,NEW SET=NFTOP,SHIFT
*NCOPY,OLD SET=NFTOP,CHANGE NUMBER=1000,NEW SET=NFBOT,SHIFT
*NFILL
NFTOP,NFBOT,5,200
*ELSET,ELSET=E2,GENERATE
*ELCOPY,OLD SET=E2,ELEMENT SHIFT=6000,SHIFT NODE=6000,NEW SET=EFOAM
*ELSET,ELSET=FOAM,GENERATE
**
*SOLID SECTION,ELSET=FOAM,MATERIAL=MFOAM
*MATERIAL,NAME=MFOAM
*ELASTIC
 18738., 0.0
*CRUSHABLE FOAM,HARDENING=VOLUMETRIC
** 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
** yield stress is in psi unit
*DENSITY
 2.85E-5,
**
** Define fluid and foam contact node set
**
*NSET, NSET=WCONT, GENERATE
*NSET, NSET=FCONT, GENERATE
**
** Define boundary condition
**
*NSET,NSET=YANODES
*NSET,NSET=YANODES,GENERATE
*NSET,NSET=XANODES
*NSET,NSET=XANODES,GENERATE
*NSET,NSET=FEND,GENERATE
*BOUNDARY
FEND,ZSYMM
XANODES,XSYMM
YANODES,YSYMM
*NSET,NSET=NALL,GENERATE
*INITIAL CONDITIONS, TYPE=VELOCITY
NALL, 3, -523.3
*ELSET,ELSET=FLUID
*ELSET,ELSET=SOLID
*ELSET,ELSET=MIDEL
*ELSET,ELSET=SHELL
*ELSET,ELSET=ELOUT
MIDEL,
SHELL,
*NSET, NSET=NOUT
*RESTART, WRITE, NUM=1, TIMEMARKS=NO
*SURFACE,TYPE=ELEMENT, NAME=SSIDE1, NO THICK
ETSHELL,SNEG
W1SHELL,SNEG
EMSHELL,SPOS
*SURFACE,TYPE=NODE,NAME=FCONT
 FCONT,
*SURFACE,TYPE=ELEMENT, NAME=SSIDE2, NO THICK
EMSHELL,SNEG
W2SHELL,SNEG
EBSHELL,SPOS
*SURFACE,TYPE=NODE,NAME=WCONT
 WCONT,
*STEP
*DYNAMIC, EXPLICIT
,0.005
*CONTACT PAIR
 SSIDE1,WCONT
*CONTACT PAIR
 SSIDE2,FCONT
*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
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
**S33,PEEQ,MISES,ERV
***ELEMENT OUTPUT,ELSET=ELOUT
**PRESS,
***ELEMENT OUTPUT,ELSET=SHELL
**SF,STH
***ELEMENT OUTPUT,ELSET=SHELL,SECTIONPOINT=1
**S33,PEEQ,MISES
***ELEMENT OUTPUT,ELSET=SHELL,SECTIONPOINT=2
**S33,PEEQ,MISES
***ELEMENT OUTPUT,ELSET=SHELL,SECTIONPOINT=3
**S33,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,NSET=NOUT
 U,
*OUTPUT,HISTORY,OP=NEW,TIME INTERVAL=1.E-5
*ENERGY OUTPUT
 ALLAE,ALLIE,ALLKE,ALLPD,ALLSE,ALLVD,ALLWK,ETOTAL
*ELEMENT OUTPUT,ELSET=MIDEL
 S33,PEEQ,MISES,ERV
*END STEP