*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=BEAM
1,0.0 ,
2,0.5,
3,1.0,0.0
*ELEMENT,TYPE=B21,ELSET=BEAM1
1,1,3
*NCOPY,OLDSET=BEAM,NEWSET=BEAM,SHIFT,CHANGENUMBER=10000
1.1, 0.0, 0.0
1,0,0, 0,1,0
*ELCOPY,ELEMENTSHIFT=1000,OLDSET=BEAM1,NEWSET=BEAM2,SHIFTNODES=10000
**************************
*SURFACE, NAME=IW_LOAD1
BEAM1,SNEG
*SURFACE,NAME=IW_LOAD2
BEAM2,SNEG
**************************
*BEAM SECTION,SECTION=RECT,TEMPERATURE=GRADIENTS,MATERIAL=SOLID,ELSET=BEAM1
1.E-2, 1.E-2
0.0, 0.0,  -1.
3,
*BEAM SECTION,SECTION=RECT,TEMPERATURE=GRADIENTS,MATERIAL=SOLID,ELSET=BEAM2
1.E-2, 1.E-2
0.0, 0.0,  -1.
3,
*MATERIAL,NAME=SOLID
*DENSITY
 1.E3,
*ELASTIC
 1.E6,0.0
**
**************************
*INCIDENT WAVE PROPERTY,NAME=IWPROP,TYPE=SPHERE
0.5,0.0
1.05,-10.0E04
*INCIDENTWAVE 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,
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, 0.0, 1.0, 0.0
9.8E0, 9.8E4, 0.0
2.5,1500,1e-12,1e-12,0.2
**************************
*BOUNDARY
BEAM,ENCASTRE
**
**************************
*RESTART,WRITE
**
*STEP
*DYNAMIC, EXPLICIT, DIRECT USER
1.E-3,0.5,
*INCIDENTWAVE,PRESSUREAMPLITUDE=WAVES_ON,PROPERTY=IWPROP
IW_LOAD1,1.0
*INCIDENTWAVE,PRESSUREAMPLITUDE=WAVES_OFF,PROPERTY=IWPROP
IW_LOAD2,1.0
**
*OUTPUT,FIELD, NUMBER INTERVAL=2
*NODE OUTPUT,NSET=BEAM
RF,
*OUTPUT,HISTORY,TIME INTERVAL=0.1E-3
*NODE OUTPUT,NSET=BEAM
RF,
*ENDSTEP