*HEADING -- Same as iw_1d_cax3_xpl_p_pp.inp but now the total wave formulation option is chosen. Note however that this option is meaningful only when the wave propagation in fluids is modeled as well. Hence, the presence of the total wave formulation option should not affect the current analysis. This input deck checks to make sure that load application on solid surfaces is not affected. -- 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 ** *ACOUSTIC WAVE FORMULATION, TYPE=TOTAL *PREPRINT,ECHO=NO,HISTORY=NO,MODEL=YES *NODE, NSET=ODB_NSET 1, 0, 11, 0, .1, 51, .5E-2,0. 61, .5E-2, .1 101, 1.E-2,0. 111, 1.E-2, .1 *NGEN,NSET=CHANNEL 1,11 101,111 *NSET,NSET=SAMPLE,GENERATE 1,11 *ELEMENT,TYPE=CAX3 1,1,102,2 2,1,101,102 *ELGEN,ELSET=CHANNEL 1,10,1,2 2,10,1,2 ************************** *SURFACE,NAME=IW_LOAD 2,S1 ************************** ** CROSS-SECTION IS Pi*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 0.0,-10.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 SAMPLE,XSYMM 11,2,2 111,2,2 ************************** *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