*Heading
basic test of incident wave interaction in ssds
*Preprint, echo=NO, model=yes, history=y, contact=NO
**
** analytical spherical wave cloads (approximate)
*Parameter
areaNode=(25.2*25.4363)/4
**hertz = 0.15915494
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
      1,   22.8551331,  -601.008118,   274.799988
*Nset, nset=EXPLOSIVE_PT-RefPt_, internal
1, 
*End Part
**  
*Part, name=STAND_OFF_PT
*Node
      1,   22.8551331,  -601.008118,   274.799988
*Nset, nset=STAND_OFF_PT-RefPt_, internal
1, 
*End Part
**  
*Part, name=Water
*Node,nset=nall
    189,   22.8551331,  -601.008118,   274.799988
    190,  -2.29828906,  -601.008118,   274.799988
    360,  -2.29828906,  -601.008118,         300.
    361,   22.8551331,  -601.008118,         300.
    583,   22.8099365,  -575.571899,   274.799988
    607,  -2.30958748,  -575.609863,   274.799988
   1657,   22.8099365,  -575.571899,         300.
   1681,  -2.30958748,  -575.609863,         300.
*Element, type=AC3D8,elset=_PickedSet3
   1,  583,  607,  190,  189, 1657, 1681,  360,  361
** Section: WATER
*Solid Section, elset=_PickedSet3, material=WATER
1.,
*End Part
**  
**
** ASSEMBLY
**
*Assembly, name=Assembly
**  
*Instance, name=Water-1, part=Water
          0.,           0.,           0.
*End Instance
*Instance, name=Water-DSLOAD, part=Water
          0.,           0.,           40.
*End Instance
**  
*Instance, name=EXPLOSIVE_PT, part=EXPLOSIVE_PT
         400,0,0
*End Instance
**  
*Instance, name=STAND_OFF_PT-1, part=STAND_OFF_PT
         0,0,0
*End Instance
**  
*Nset, nset=EXPLOSIVE_PT, instance=EXPLOSIVE_PT
 1,
*Nset, nset=STAND_OFF_PT, instance=STAND_OFF_PT-1
 1,
*Nset, nset=cload, instance=water-dsload
 189, 361, 583, 1657
*Elset, elset=__PickedSurf15_S6, internal, instance=Water-1
 1,
*Surface, type=ELEMENT, name=acousticLoad
__PickedSurf15_S6, S6
*Elset, elset=qickedSurf15_S6, internal, instance=Water-DSLOAD
 1,
*Surface, type=ELEMENT, name=acousticLoad2
qickedSurf15_S6, S6
*End Assembly
***
***
** MATERIALS
** 
*Material, name=BARGE_HULL
*Density
 0.000253623,
*Elastic
 2e+11, 0.34
*Material, name=WATER
*Acoustic Medium
317633.,
*Density
 9.3427e-05,
** 
** INTERACTION PROPERTIES
** 
*Incident Wave Interaction Property, name=Sphere1, type=sphere
59059.2, 9.3427e-05
*Incident Wave Interaction Property, name=plane, type=plane
59059.2, 9.3427e-05
*AMPLITUDE, DEFINITION=TABULAR, NAME=PRESSUREVTIME
0.0, 1.0,0.5,1.,1.0,1.
** ----------------------------------------------------------------
** ----------------------------------------------------------------
*step, name=Eigen
*frequency,eig=lanc,norm=mass, SIM=NO
,,1000.00,
*output,field,freq=0
*output,histo,freq=0
*nodeprint,freq=0
*End Step
** ----------------------------------------------------------------
*Step, name=SphereLoad
*steadystateDynamic, sub
<hertz>,
** 
** INTERACTIONS
** 
** Interaction: INCIDENT_WAVE
*incident wave interaction, PRESSURE AMPLITUDE=PRESSUREVTIME,property=Sphere1,load case=2
acousticLoad, EXPLOSIVE_PT, STAND_OFF_PT, 1.
**
*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=PlaneLoad
*steadystateDynamic, sub
<hertz>,
** 
** INTERACTIONS
** 
** Interaction: INCIDENT_WAVE
*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