*HEADING 1D transient wave propagation, solid medium, incident wave load, planar and defined through pressure amplitude ** ** Analytical solution to prescribed constant acceleration of a0 at z=0 ** ** p(z,t) = (K a0 / c^2) (ct - z) for z < ct ** = 0 for z > ct ** *PREPRINT,ECHO=NO,HISTORY=NO,MODEL=YES *NODE, NSET=PLATE 1, 0. , 0. , 2, 1. , 0. , *ELEMENT,TYPE=SAX1,ELSET=PLATE 1,2,1 ************************** *SURFACE,NAME=IW_LOAD PLATE,SPOS ************************** **CROSS-SECTION IS Pi *SHELL SECTION,MATERIAL=SOLID,ELSET=PLATE 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 *INCIDENTWAVE FLUID PROPERTY 0.1424E6,1.21,0.0 *AMPLITUDE,DEFINITION=TABULAR,NAME=PRESSUREVTIME 0.0,0.0,0.0005,1.0 ************************** *BOUNDARY PLATE,ENCASTRE ** ************************** *RESTART,WRITE ** *STEP *DYNAMIC,EXPLICIT, DIRECT USER 0.2E-4,.0005, ** *INCIDENTWAVE,PRESSUREAMPLITUDE=PRESSUREVTIME,PROPERTY=IWPROP IW_LOAD,1.E3 ** *OUTPUT,FIELD,NUMBER INTERVAL=2 *NODE OUTPUT,NSET=PLATE RF, *OUTPUT,HISTORY,TIME INTERVAL=0.1E-3 *NODE OUTPUT,NSET=PLATE RF, *ENDSTEP ** **