*HEADING
 1D transient wave propagation, solid-acoustic coupling (between AC3D8 and S4R elements),
 incident wave load, planar and defined through pressure amplitude
**
** C = root (K/rho)= 343.05 m/s   (in acoustic medium)
**
** C = root (K/rho)=  31.62 m/s   (in solid medium) K=bulk mod, rho=density.
**
**
*PREPRINT,ECHO=NO,HISTORY=NO,MODEL=YES
*NODE,NSET=ACOUSTIC
*NODE, NSET=SOLID
*ELEMENT,TYPE=AC3D8,ELSET=ACOUSTIC
*ELEMENT,TYPE=S4R,ELSET=SOLID
**************************
*SURFACE,NAME=IW_LOAD_ACOUSTIC
1,S1
*SURFACE,NAME=IW_LOAD_SOLID
SOLID,SPOS
*SURFACE,NAME=TIE_ACOUSTIC
1,S1
*SURFACE,NAME=TIE_SOLID
SOLID,SPOS
**************************
*SOLID SECTION,MATERIAL=AIR,ELSET=ACOUSTIC
**
*MATERIAL,NAME=AIR
*DENSITY
 1.21,
*ACOUSTIC MEDIUM,BULK MODULUS
**
*SHELL SECTION,MATERIAL=SOLID,ELSET=SOLID
 1.E-4,
**
*MATERIAL,NAME=SOLID
*DENSITY
 1.E3,
*ELASTIC
 1.E6,0.0
**
**************************
*INCIDENT WAVE PROPERTY,NAME=IWPROP,TYPE=PLANE
*INCIDENT WAVE FLUID PROPERTY
*AMPLITUDE,DEFINITION=TABULAR,NAME=PRESSUREVTIME_ACOUSTIC
0.0,0.0,1.0,-1.0
*AMPLITUDE,DEFINITION=TABULAR,NAME=PRESSUREVTIME_SOLID
0.0,0.0,1.0,1.0
**************************
**
*BOUNDARY
SOLID,ENCASTRE
**************************
**
*TIE,NAME=COUPLING
TIE_ACOUSTIC,TIE_SOLID
**************************
*RESTART,WRITE
**
*STEP,INC=1000,NLGEOM
*DYNAMIC,ALPHA=0.0
.0001,.005
**
*INCIDENT WAVE,PRESSUREAMPLITUDE=PRESSUREVTIME_ACOUSTIC,PROPERTY=IWPROP
IW_LOAD_ACOUSTIC,415.09517
*INCIDENT WAVE,PRESSUREAMPLITUDE=PRESSUREVTIME_SOLID,PROPERTY=IWPROP
IW_LOAD_SOLID,1.E3
**
*NODE PRINT,NSET=ACOUSTIC,FREQ=20
 POR,
*NODE PRINT,NSET=SOLID,FREQ=20
 RF,
*OUTPUT,FIELD,VAR=PRESELECT,FREQ=1
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT,NSET=ACOUSTIC
 POR,
*NODE OUTPUT,NSET=SOLID
 RF,
*OUTPUT,HISTORY,FREQ=1
*NODE OUTPUT,NSET=ACOUSTIC
 POR,
*NODE OUTPUT,NSET=SOLID
 RF,
*ENDSTEP