*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, 51, 0,.5E-2 61, .1,.5E-2 101, 0, 1.E-2 111, .1, 1.E-2 *NGEN,NSET=CHANNEL 1,11 101,111 *NSET,NSET=SAMPLE,GENERATE 1,11 *ELEMENT,TYPE=CPE3 1,1,2,102 2,102,101,1 *ELGEN,ELSET=CHANNEL 1,10,1,2 2,10,1,2 ************************** *SURFACE,NAME=IW_LOAD 2,S2 ************************** ** CROSS-SECTION IS 1.E-4 *SOLID SECTION,MATERIAL=SOLID,ELSET=CHANNEL 1.E-2, ** *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,2 11,1,1 111,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