*HEADING
 1D transient wave propagation in solid medium with incident wave load
 defined through bubble-load amplitude
**
*PREPRINT,ECHO=NO,HISTORY=NO,MODEL=YES
*NODE, NSET=ODB_NSET
   1, 0,
  11, .1,
  51, 0, .5E-2
  61, .1, .5E-2
 101, 0, 1.E-2
 111, .1, 1.E-2
*NGEN,NSET=CHANNEL
   1,11
  51,61
 101,111
*ELEMENT,TYPE=CPE4R
 1,1,2,102,101
*ELGEN,ELSET=CHANNEL1
 1,10
*NCOPY,OLDSET=CHANNEL,NEWSET=CHANNEL,SHIFT,CHANGENUMBER=10000
0.0,0.03,0.0
1,0,0, 0,1,0
*NSET,NSET=QA_TEST,GENERATE
 1,11
10001,10010
*ELCOPY,ELEMENTSHIFT=1000,OLDSET=CHANNEL1,NEWSET=CHANNEL2,SHIFTNODES=10000
**************************
*SURFACE,NAME=IW_LOAD1
1,S4
*SURFACE,NAME=IW_LOAD2
1001,S4
**************************
**
*SOLID SECTION,MATERIAL=SOLID,ELSET=CHANNEL1
 1.E-2, 
*SOLID SECTION,MATERIAL=SOLID,ELSET=CHANNEL2
 1.E-2, 
**
*MATERIAL,NAME=SOLID
*DENSITY
 1.E3,
*ELASTIC
 1.E6,0.0
**************************
*INCIDENT WAVE PROPERTY,NAME=IWPROP,TYPE=SPHERE
0.0,0.5E-2,
-10.0E04,0.0,
*INCIDENT WAVE FLUID PROPERTY
2.25E9,1.0E3,0.0
*AMPLITUDE,DEFINITION=BUBBLE,NAME=WAVES_ON
52.1E6, 0.09E-3, 0.18, 0.185, 839.6E6, 1.27, 1.5E3, 226.8, 137.16
1.0E3, 1.5E3, 1.0, 0.0, 0.0
9.8E0, 9.8E4,1.0
2.5,1400,1e-12,1e-12,0.2
*AMPLITUDE,DEFINITION=BUBBLE,NAME=WAVES_OFF
52.1E6, 0.09E-3, 0.18, 0.185, 839.6E6, 1.27, 1.5E3, 226.8, 137.16
1.0E3, 1.5E3, 1.0, 0.0, 0.0
9.8E0, 9.8E4,0.0
2.5,1500,1e-12,1e-12,0.2
**************************
*RESTART,WRITE,FREQUENCY=1000
*BOUNDARY
CHANNEL,2,2
 11,1,1
111,1,1
10011,1,1
10111,1,1
**************************
**
*STEP,INC=10000
*DYNAMIC,alpha=0.0,haftol=10.0
.0001,0.5,,0.005
**
*INCIDENT WAVE,PRESSUREAMPLITUDE=WAVES_ON, PROPERTY=IWPROP
IW_LOAD1,1.0
*INCIDENT WAVE,PRESSUREAMPLITUDE=WAVES_OFF,PROPERTY=IWPROP
IW_LOAD2,1.0
**
*EL PRINT,FREQ=1000
 S,
 E,
*NODE PRINT,NSET=QA_TEST,FREQ=1000
 U,
*OUTPUT,VAR=PRESELECT,FIELD,FREQ=1
*OUTPUT,VAR=PRESELECT,HISTORY,FREQ=1
*OUTPUT,FIELD, FREQUENCY=1000
*NODE OUTPUT
 U
*ELEMENT OUTPUT
 S,E
*OUTPUT,HISTORY,FREQUENCY=10
*NODE OUTPUT,NSET=QA_TEST
 U,
*ENDSTEP