*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=ODB_NSET 1, 0., 11, .1, 101, 0., 1E-2 111, .1, 1E-2 151, 0. , 1E-2, -1E-2 161, .1, 1E-2, -1E-2 251, 0., 0, -1E-2 261, .1, 0, -1E-2 *NGEN,NSET=CHANNEL 1,11 101,111 151,161 251,261 *NSET,NSET=SAMPLE,GENERATE 1,11 *ELEMENT,TYPE=C3D6 1,1,102,2,251,152,252 2,1,101,102,251,151,152 *ELGEN,ELSET=CHANNEL 1,10,1,2 2,10,1,2 ************************** *SURFACE,NAME=IW_LOAD 2,S3 ************************** *SOLID SECTION,MATERIAL=SOLID,ELSET=CHANNEL ** *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.0,0.0 *INCIDENT WAVE FLUID PROPERTY 0.1424E6,1.21,0.0 *AMPLITUDE,DEFINITION=TABULAR,NAME=PRESSUREVTIME 0.0,0.0,0.0015,-1.0 ************************** ** *BOUNDARY CHANNEL,2,3 11,1,1 111,1,1 161,1,1 261,1,1 ************************** *RESTART,WRITE ** *STEP *DYNAMIC,EXPLICIT ,.0015, ** *INCIDENT WAVE,PRESSUREAMPLITUDE=PRESSUREVTIME,PROPERTY=IWPROP IW_LOAD,415.09517 ** *NSET,NSET=QA_TEST SAMPLE, *OUTPUT,FIELD,NUMBER INTERVAL=2 *NODE OUTPUT U, *NODE OUTPUT,NSET=QA_TEST U, *ELEMENT OUTPUT S, *OUTPUT,HISTORY,TIME INTERVAL=0.1E-3 *NODE OUTPUT,NSET=SAMPLE U, *ENDSTEP