*Heading
basic test of incident wave interaction in ssdd
*Preprint, echo=NO, model=yes, history=y, contact=NO
**
** analytical spherical wave cloads (approximate)
*Parameter
areaNode=(25.2*25.4363)/4
hertz = 100
omega =2*3.1415*hertz
rhof  = 9.3427e-05
cccc  = 59059.2
** if *IWI load case=1, then
imagLoadCase1 = -areaNode*omega/rhof/cccc
realLoadCase1 = -areaNode/rhof/400
** if *IWI load case=2, then
imagLoadCase2 = -areaNode/rhof/400
realLoadCase2 = -areaNode*omega/rhof/cccc
**
*Part, name=EXPLOSIVE_PT
*Node
*Nset, nset=EXPLOSIVE_PT-RefPt_, internal
*End Part
**
*Part, name=STAND_OFF_PT
*Node
*Nset, nset=STAND_OFF_PT-RefPt_, internal
*End Part
**
*Part, name=Water
*Node,nset=nall
*Element, type=AC3D8,elset=_PickedSet3
** Section: WATER
*Solid Section, elset=_PickedSet3, material=WATER
1.,
*End Part
**
**
** ASSEMBLY
**
*Assembly, name=Assembly
**
*Instance, name=Water-1, part=Water
*End Instance
*Instance, name=Water-DSLOAD, part=Water
*End Instance
**
*Instance, name=EXPLOSIVE_PT, part=EXPLOSIVE_PT
*End Instance
**
*Instance, name=STAND_OFF_PT-1, part=STAND_OFF_PT
*End Instance
**
*Nset, nset=EXPLOSIVE_PT, instance=EXPLOSIVE_PT
*Nset, nset=STAND_OFF_PT, instance=STAND_OFF_PT-1
*Nset, nset=cload, instance=water-dsload
*Elset, elset=__PickedSurf15_S6, internal, instance=Water-1
*Surface, type=ELEMENT, name=acousticLoad
__PickedSurf15_S6, S6
*Elset, elset=qickedSurf15_S6, internal, instance=Water-DSLOAD
*Surface, type=ELEMENT, name=acousticCLoad
qickedSurf15_S6, S6
*End Assembly
***
***
** MATERIALS
**
*Material, name=BARGE_HULL
*Density
 0.000253623,
*Elastic
 2e+11, 0.34
*Material, name=WATER
*Acoustic Medium
*Density
 9.3427e-05,
**
** INTERACTION PROPERTIES
**
*Incident Wave Interaction Property, name=Sphere1, type=sphere
*Incident Wave Interaction Property, name=plane, type=plane
*AMPLITUDE, DEFINITION=TABULAR, NAME=PRESSUREVTIME
0.0, 1.0,0.5,1.,1.0,1.
** ----------------------------------------------------------------
**
** STEP: Step-1
**
*Step, name=Step-1
*steadystateDynamic, direct
<hertz>,
**
** INTERACTIONS
**
** Interaction: INCIDENT_WAVE
*incident wave interaction, PRESSURE AMPLITUDE=PRESSUREVTIME,property=Sphere1,load case=2
acousticLoad, EXPLOSIVE_PT, STAND_OFF_PT, 1.
*incident wave interaction, PRESSURE AMPLITUDE=PRESSUREVTIME,property=Sphere1,load case=1
acousticLoad, EXPLOSIVE_PT, STAND_OFF_PT, 1.
**
** IWI causes real (loadcase=1) and imag (loadcase=2) loads for either specified loadcase
*cload,loadcase=1
cload,8,<realLoadCase1>
*cload,loadcase=2
cload,8,<imagLoadCase1>
**
*cload,loadcase=1
cload,8,<realLoadCase2>
*cload,loadcase=2
cload,8,<imagLoadCase2>
**
*nodeprint
por,
*Output, field
*Node Output
por
*End Step
** ----------------------------------------------------------------
*Step, name=Step-2
*steadystateDynamic, direct
<hertz>,
**
** INTERACTIONS
**
** Interaction: INCIDENT_WAVE
*incident wave interaction, PRESSURE AMPLITUDE=PRESSUREVTIME,property=plane,load case=1
acousticLoad, EXPLOSIVE_PT, STAND_OFF_PT, 1.
*incident wave interaction, PRESSURE AMPLITUDE=PRESSUREVTIME,property=plane,load case=2
acousticLoad, EXPLOSIVE_PT, STAND_OFF_PT, 1.
**
** same comment pertains as in previous step
** also: here for plane wave case, take lim rj->inf to find correct cloads.
*cload,loadcase=2
cload,8,<imagLoadCase1>
*cload,loadcase=1
cload,8,<realLoadCase2>
**
*nodeprint
por
*Output, field
*Node Output
por
*End Step