*HEADING
 Fluid sphere splashing in a tank
*preprint,model=yes
*parameter
** mesh
scale = 1.0
xLen = 15.0*scale
yLen = 5.0*scale
zLen = 15.0*scale
** tank rigid walls
xRight = -xLen/2. 
yBot   = -yLen/2. 
zFront = -zLen/2. 
xLeft  = -xRight
yTop   = -yBot
zBack  = -zFront
thickWalls = 0.00001
particleRadius = 0.2
particleRadiusVol = 0.25
gravMag      = -.00025*scale
**------------------------------------------------
** sphere
*system
0,14,0,   1,14,0
0,9,0
*include, input=ver_prc_watersplashinpan_sphere.inp
*surface, type=node, name=sphereSPHNodes
sphereNodes,
**------------------------------------------------
*system
0,0,0,   1,0,0
0,1,0
*include, input=ver_prc_watersplashinpan_tank.inp
*surface, type=node, name=tankSPHNodes
tankNodes,
**------------------------------------------------
*elset, elset=allParticles
tankPC3D, spherePC3D
*solid section, elset=allParticles, material=water, controls=SPH_fac
<particleRadiusVol>,
*sectioncontrols, name=SPH_fac


-8, -3.5 ,-8 ,8 ,19 , 8, 
**------------------------------------------------
*material,name=water
*Density
 1.0
*Eos, type=USUP
0.111803398875,0.,0.
**------------------------------------------------
*node, nset=wallNodes
500001, <xLeft>,  <yBot>, <zFront>
500002, <xRight>, <yBot>, <zFront>
500003, <xRight>, <yBot>, <zBack>
500004, <xLeft>,  <yBot>, <zBack>
600001, <xLeft>,  <yTop>, <zFront>
600002, <xRight>, <yTop>, <zFront>
600003, <xRight>, <yTop>, <zBack>
600004, <xLeft>,  <yTop>, <zBack>
700001, 0,        <yBot>, 0
*element, type =s4R, elset=tankWalls
500001, 500001, 600001, 600004, 500004
500002, 500001, 600001, 600002, 500002
500003, 500002, 600002, 600003, 500003
500004, 500004, 600004, 600003, 500003
500005, 500001, 500002, 500003, 500004
*shell section, elset=tankWalls, material=steel
<thickWalls>,
*material, name=steel
*density
7.8e-9
*elastic
2.1e5,0.0
*Rigid body, elset=tankWalls, ref=700001
*Surface, type=element, name=rigidWallsInteriorSurf
500001, SNEG
500002, SPOS
500003, SPOS
500004, SNEG
500005, SPOS
**------------------------------------------------
*Surface Interaction, name=nodeToFace
**------------------------------------------------
*elset,elset=watersplashinpan_QATEST
6519,
*nset,nset=watersplashinpan_QATEST
2975,
*nset,nset=tank_QA_TEST
700001,
**
*Step,nlgeom=yes
*dynamic,explicit
**,element
,5000
*bulk viscosity
0.06,1.2
**
** fix the rigid tank
*boundary
700001, 1,6
**
*cload
sphereNodes, 2,  -0.000005
tankNodes  , 2,  -0.000005
*contact
*contact inclusions
sphereSPHNodes, rigidWallsInteriorSurf
tankSPHNodes  , rigidWallsInteriorSurf
***Contact Property Assignment
**sphereSPHNodes, rigidWallsInteriorSurf, nodeToFace
**tankSPHNodes, rigidWallsInteriorSurf, nodeToFace
***Surface Property Assignment, property=THICKNESS
**sphereSPHNodes, <particleRadius>
**tankSPHNodes, <particleRadius>
**rigidWallsInteriorSurf, 0.0
**
*OUTPUT,FIELD, num=20
*NODE OUTPUT
 U,V,A
*element output
 S,
*CONTACT OUTPUT 
 CFORCE,
*Output, history, freq=10
*Energy Output
ETOTAL,
*****************************************************
*OUTPUT, HISTORY, FREQ=10
*ELEMENT OUTPUT, ELSET =watersplashinpan_QATEST
S,
*NODE OUTPUT, NSET =watersplashinpan_QATEST
U,
*NODE OUTPUT, NSET =tank_QA_TEST
U,
**
*end step