*HEADING 1D transient wave propagation, solid medium, incident wave load, planar and defined through pressure amplitude Use *BEAM FLUID INERTIA with B32, *BEAM GENERAL SECTION, ROTARY INERTIA=ISOTROPIC ** *PREPRINT,ECHO=NO,HISTORY=NO,MODEL=YES *NODE,NSET=BEAM 1,0.0 , 2,0.5, 3,1.0,0.0 *ELEMENT,TYPE=B32,ELSET=BEAM 1,1,2,3 ************************** *SURFACE, NAME=IW_LOAD BEAM,S1 ************************** *BEAM GENERAL SECTION,SECTION=RECT,DENSITY=1.E3,ELSET=BEAM, ROTARY INERTIA=ISOTROPIC 2.E-2, 1.E-2 0.0, 0.0, -1. 1.E6 *BEAM FLUID INERTIA, FULL 1.21,0.0,0.0,0.01,1.0,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,.0010,0.0 ************************** *BOUNDARY BEAM,ENCASTRE ** ************************** *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=BEAM RF, *ENDSTEP ** **