*HEADING
 1D transient wave propagation, solid-acoustic coupling (between ACAX4 and SAX1 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
 1, 0, 0
 2, 1, 0
 3, 1, 1
 4, 0, 1
*NODE, NSET=SOLID
 1001, 0, 1
 1002, 1, 1
*ELEMENT,TYPE=ACAX4,ELSET=ACOUSTIC
 1, 1,2,3,4
*ELEMENT,TYPE=SAX1,ELSET=SOLID
 1001, 1002,1001
*NSET,NSET=SYMM
 1001
**************************
*SURFACE,NAME=IW_LOAD_ACOUSTIC
1,S3
*SURFACE,NAME=IW_LOAD_SOLID
1001,SPOS
*SURFACE,NAME=TIE_ACOUSTIC
1,S3
*SURFACE,NAME=TIE_SOLID
1001,SPOS
**************************
** CROSS-SECTION IS Pi*1.E-4  
*SOLID SECTION,MATERIAL=AIR,ELSET=ACOUSTIC
 1.E-4, 
**
*MATERIAL,NAME=AIR
*DENSITY
 1.21,
*ACOUSTIC MEDIUM,BULK MODULUS
 0.1424E6,
**
*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
0.0,0.0,
0.0,-10.0,
*INCIDENT WAVE FLUID PROPERTY
0.1424E6,1.21,0.0
*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
SYMM,XSYMM
1001,2,2
1002,1,2
**************************
**
*TIE,NAME=COUPLING
TIE_ACOUSTIC,TIE_SOLID
**************************
*RESTART,WRITE
**
*STEP,INC=1000
*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