*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 = 0
realLoadCase1 = -areaNode
** if *IWI load case=2, then
imagLoadCase2 = areaNode
realLoadCase2 = 0
**  
*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=qa_test
    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=C3D8,elset=_PickedSet3
   1,  583,  607,  190,  189, 1657, 1681,  360,  361
** Section: WATER
*Solid Section, elset=_PickedSet3, material=BARGE_HULL
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
*Nset,nset=qa_test
water-1.qa_test,water-dsload.qa_test
*Nset,nset=BC
water-1.361,water-1.583,water-1.189,water-dsload.361,water-dsload.583,water-dsload.189
*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.
*Boundary
BC,2,3,0.0
** ----------------------------------------------------------------
** ----------------------------------------------------------------
*step, name=Eigen
*frequency,eig=lanc,norm=mass
,,1000.00,
*output,field,freq=0
*output,histo,freq=0
*nodeprint,freq=0
*End Step
** ----------------------------------------------------------------
*Step, name=Subspace
*steadystateDynamic, sub
<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=sphere1,load case=1
acousticLoad2, EXPLOSIVE_PT, STAND_OFF_PT, 1.
*nodeprint
u1,
*Output, field
*Node Output,nset=qa_test
u
*End Step
** ----------------------------------------------------------------
*Step, name=Direct
*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=sphere1,load case=1
acousticLoad2, EXPLOSIVE_PT, STAND_OFF_PT, 1.
**
*nodeprint
u1,
*Output, field
*Node Output,nset=qa_test
u
*End Step