*HEADING 1D transient wave propagation, solid medium, incident wave load, planar and defined through pressure amplitude Use *BEAM FLUID INERTIA,FULL and HALF with B31, *BEAM GENERAL SECTION ** *PREPRINT,ECHO=NO,HISTORY=NO,MODEL=NO *NODE,NSET=BEAM 1,0.0 , 2,0.5, 3,1.0,0.0 *NCOPY,OLDSET=BEAM,NEWSET=BEAMHALF,CHANGE=100,SHIFT 0,0,1 1,0,0, 0,1,0 *NCOPY,OLDSET=BEAM,NEWSET=BEAMDRY,CHANGE=200,SHIFT 0,0,2 1,0,0, 0,1,0 *NSET,NSET=ALL,GENERATE 1,99999 *ELEMENT,TYPE=B31,ELSET=BEAM 1,1,3 *ELEMENT,TYPE=B31,ELSET=BEAMHALF 2,101,103 *ELEMENT,TYPE=B31,ELSET=BEAMDRY 3,201,203 ************************** *SURFACE, NAME=IW_LOAD BEAM,S1 BEAMHALF,S1 BEAMDRY,S1 ************************** *BEAM GENERAL SECTION,SECTION=RECT,DENSITY=1.E3,ELSET=BEAM 2.E-2, 1.E-2 0.0, 0.0, -1.0 1.E6 *BEAM FLUID INERTIA, FULL 1000.0,0.0,0.0,0.01,1.0,0.0 *BEAM GENERAL SECTION,SECTION=RECT,DENSITY=1.E3,ELSET=BEAMHALF 2.E-2, 1.E-2 0.0, 0.0, -1. 1.E6 *BEAM FLUID INERTIA,HALF 1000.0,0.0,0.0,0.01,1.0,0.0 *BEAM GENERAL SECTION,SECTION=RECT,DENSITY=1.E3,ELSET=BEAMDRY 2.E-2, 1.E-2 0.0, 0.0, -1. 1.E6 *BEAM FLUID INERTIA,HALF 1.0,0.0,0.0,0.01,1.0e-3,0.0 ************************** *INCIDENT WAVE PROPERTY,NAME=IWPROP,TYPE=PLANE 0.5,0.0 0.5,10.0 *INCIDENTWAVE FLUID PROPERTY 0.1424E6,1.21,0.0 *AMPLITUDE,DEFINITION=SMOOTH STEP,NAME=PRESSUREVTIME 0.0,0.0,0.0005,1.0,0.0010,0.0 ************************** *BOUNDARY ALL,1,1,0.0 ALL,3,3,0.0 ** ************************** *RESTART,WRITE ** *STEP *DYNAMIC, EXPLICIT, DIRECT .0001,.0010 *INCIDENTWAVE,PRESSUREAMPLITUDE=PRESSUREVTIME,PROPERTY=IWPROP IW_LOAD,1.E3 ** *OUTPUT,FIELD, OP=ADD,NUMBER=1 *NODE OUTPUT RF, *OUTPUT,HISTORY,OP=ADD,FREQUENCY=1 *NODE OUTPUT,NSET=ALL RF, *ENDSTEP ** **