*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 x=0 ** ** p(x,t) = (K a0 / c^2) (ct - x) for x < ct ** = 0 for x > ct ** *PREPRINT,ECHO=NO,HISTORY=NO,MODEL=YES *NODE, NSET=PLATE 1, 0. , 0. , 2, 1. , 0. , 3, 1. , 1. , 4, 0. , 1. , 5, 0.5, 0. , 6, 1.0, 0.5, 7, 0.5, 1.0, 8, 0. , 0.5, 9, 0.5, 0.5, *ELEMENT,TYPE=S4R,ELSET=PLATE 1,1,2,3,4 ************************** *SURFACE,NAME=IW_LOAD PLATE,SPOS ************************** *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.5,0.5,0.0 0.5,0.5,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 ** **