*Heading
Underwater Explosion Simulation Benchmark: Spherical Shell and Step Wave
**
** PARTS
**
*Part, name=fluid
*End Part
*Part, name=sphere
*End Part
**
*PART,NAME=DUMMY
**************************
** THIS IS A DUMMY PART USED TO INCLUDE
** THE ANALYTIC SOLUTION TO
** THE OUTPUT DATABASE (.ODB) FILE. THE ANALYTIC
** SOLUTION IS FOUND IN:
** H. HUANG, JASA VOL 45, NO 3, PP661-670, 1969.
** THESE ANALYTIC RESULTS ARE
** FOR VELOCITY (V2) AT THE FINITE ELEMENT NODES 4 AND 5:
** THE DUMMY V1 CORRESPONDS TO NODE 4, V2 TO NODE 5.
**************************
*NODE,NSET=ANALYTIC
*ELEMENT,TYPE=MASS,ELSET=ANALYTIC
*MASS,ELSET=ANALYTIC
1.0
*END PART
**
** ASSEMBLY
**
*Assembly, name=Assembly
**
*Instance, name=Dummy, part=dummy
*End Instance
**
*Instance, name=sphere-1, part=sphere
*Node
*Element, type=S4R
** Region: (shell:Picked)
*Elset, elset=_I1, generate
** Section: shell
*Shell Section, elset=_I1, material=shell
0.02, 5
*End Instance
**
*Instance, name=fluid-1, part=fluid
*Node
*Element, type=AC3D8R
** Region: (fluid:Picked)
*Elset, elset=_I1, generate
** Section: fluid
*Solid Section, elset=_I1, material=fluid
1.,
*End Instance
*Nset, nset=sphere, instance=sphere-1, generate
*Elset, elset=sphere, instance=sphere-1, generate
*Nset, nset=fluid, instance=fluid-1, generate
*Elset, elset=fluid, instance=fluid-1, generate
*Nset, nset=zsymm, instance=sphere-1
*Elset, elset=zsymm, instance=sphere-1, generate
*Nset, nset=xsymm, instance=sphere-1
*Elset, elset=xsymm, instance=sphere-1
*Nset, nset=enda, instance=fluid-1
*Nset, nset=ends, instance=sphere-1
*Elset, elset=_shell-fluid_SPOS, instance=sphere-1, generate
*Surface, type=ELEMENT, name=shell-fluid
_shell-fluid_SPOS, SPOS
*Elset, elset=_fluid-shell_S2, instance=fluid-1
*Surface, type=ELEMENT, name=fluid-shell
_fluid-shell_S2, S2
*Elset, elset=_radiation_S3, instance=fluid-1
*Elset, elset=_radiation_S4, instance=fluid-1
*Elset, elset=_radiation_S2, instance=fluid-1
*Surface, type=ELEMENT, name=radiation
_radiation_S3, S3
_radiation_S4, S4
_radiation_S2, S2
** Constraint: Constraint-1
*Tie, name=Constraint-1, adjust=yes
fluid-shell, shell-fluid
*End Assembly
**
** MATERIALS
**
*Material, name=fluid
*Acoustic Medium
*Density
1.,
*Material, name=shell
*Density
 7.79,
*Elastic
97., 0.3
**
** ANALYTICAL SOLUTION AMPLITUDES
**
*Amplitude,name=Analytic_V2_node4,definition=tabular
   0.000000E+00,   1.031800E-13,   1.000000E-03,   3.979600E-07,   2.000000E-03,   1.285600E-06,   3.000000E-03,   2.367600E-06,
   4.000000E-03,   3.497200E-06,   5.000000E-03,   4.613700E-06,   6.000000E-03,   5.701800E-06,   7.000000E-03,   6.767200E-06,
   8.000000E-03,   7.822600E-06,   9.000000E-03,   8.879700E-06,   1.000000E-02,   9.946900E-06,   1.100000E-02,   1.102700E-05,
   1.200000E-02,   1.212100E-05,   1.300000E-02,   1.322500E-05,   1.400000E-02,   1.433300E-05,   1.500000E-02,   1.544200E-05,
   1.600000E-02,   1.654500E-05,   1.700000E-02,   1.763900E-05,   1.800000E-02,   1.872100E-05,   1.900000E-02,   1.979000E-05,
   2.000000E-02,   2.084500E-05,   3.000000E-02,   3.093100E-05,   4.000000E-02,   4.059300E-05,   5.000000E-02,   4.938300E-05,
   6.000000E-02,   5.755800E-05,   7.000000E-02,   6.519700E-05,   8.000000E-02,   7.207300E-05,   9.000000E-02,   7.833900E-05,
   0.100000    ,   8.417200E-05,   0.110000    ,   8.946900E-05,   0.120000    ,   9.415000E-05,   0.130000    ,   9.835100E-05,
*Amplitude,name=Analytic_V2_node5,definition=tabular
   0.000000E+00,  -2.586500E-13,   1.000000E-03,  -6.638000E-11,   2.000000E-03,  -6.771100E-11,   3.000000E-03,   4.940000E-11,
   4.000000E-03,   2.058000E-10,   5.000000E-03,   3.178900E-10,   6.000000E-03,   3.430600E-10,   7.000000E-03,   2.844500E-10,
   8.000000E-03,   1.765000E-10,   9.000000E-03,   6.701000E-11,   1.000000E-02,   1.213000E-12,   1.100000E-02,   1.162800E-11,
   1.200000E-02,   1.129000E-10,   1.300000E-02,   3.019600E-10,   1.400000E-02,   5.610900E-10,   1.500000E-02,   8.628200E-10,
   1.600000E-02,   1.175500E-09,   1.700000E-02,   1.468500E-09,   1.800000E-02,   1.716300E-09,   1.900000E-02,   1.901300E-09,
   2.000000E-02,   2.015200E-09,   3.000000E-02,   2.116600E-09,   4.000000E-02,   6.444000E-09,   5.000000E-02,   7.350400E-09,
   6.000000E-02,   1.011900E-08,   7.000000E-02,   1.690100E-08,   8.000000E-02,   1.900400E-08,   9.000000E-02,   2.032700E-08,
   0.100000    ,   2.811300E-08,   0.110000    ,   3.640200E-08,   0.120000    ,   3.801800E-08,   0.130000    ,   3.857100E-08,
**
** BOUNDARY CONDITIONS FOR FINITE ELEMENT MODEL
**
** Name: BC-1 Type: Symmetry/Antisymmetry/Encastre
*Boundary
xsymm, XSYMM
** Name: BC-2 Type: Symmetry/Antisymmetry/Encastre
*Boundary
zsymm, ZSYMM
**
**
*IMPEDANCE PROPERTY, TYPE=SPHERE, NAME=SPHERE
*INCIDENTWAVEPROPERTY,NAME=PLANE
*INCIDENTWAVEFLUIDPROPERTY
*amplitude,definition=smoothstep,name=smoothstep
-3.0,0.0, 0.0,0.0,  0.1,1.0,  15.0,1.0
*amplitude,definition=tabular,name=step,time=totaltime
-3.0,0.0, 0.0,0.0,  0.1,1.0,  15.0,1.0
**
** ----------------------------------------------------------------
**
** STEP: xpl
**
*Step, name=xpl
*Dynamic, Explicit, direct user control
0.01, 10.
*Bulk Viscosity
*incidentwave,pressureamplitude=step,property=plane
assembly_shell-fluid,1.0e-4
assembly_fluid-shell,1.0e-4
*SIMPEDANCE
assembly_radiation,SPHERE
**
** BOUNDARY CONDITIONS FOR DUMMY (ANALYTIC) MODEL
**
*Boundary,Amplitude=Analytic_v2_node4,type=velocity
Assembly.DUMMY.ANALYTIC,1,1,1.0
*Boundary,Amplitude=Analytic_v2_node5,type=velocity
Assembly.DUMMY.ANALYTIC,2,2,1.0
*Boundary
Assembly.DUMMY.ANALYTIC,3,3,0.0
**
*Output, field,numberinterval=25
*Node output
u,v,a
POR,
*Output, history
*Node Output, nset=Assembly.DUMMY.ANALYTIC
v,
*Node Output, nset=ends
V,
*Node Output, nset=enda
POR,
*End Step