*HEADING -- Same as iw_1d_b21_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, spherical and defined through pressure amplitude ** ** C = root (K/rho)= 31.22 m/s (in the solid medium) K=bulk mod, rho=density. ** ** 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 ** *ACOUSTIC WAVE FORMULATION, TYPE=TOTAL *PREPRINT,ECHO=NO,HISTORY=NO,MODEL=YES *NODE,NSET=BEAM 1,0. , 2,0.5, 3,1.0, *ELEMENT,TYPE=B21,ELSET=BEAM 1,1,3 ************************** *SURFACE, NAME=IW_LOAD BEAM,SPOS ************************** *BEAM SECTION,SECTION=RECT,TEMPERATURE=GRADIENTS,MATERIAL=SOLID,ELSET=BEAM 1.E-2, 1.E-2 0., 0., -1. 3, *MATERIAL,NAME=SOLID *DENSITY 1.E3 *ELASTIC 1.E6,0.0 ** ************************** *INCIDENT WAVE PROPERTY,NAME=IWPROP,TYPE=SPHERE 0.5,0.0 0.5,10.0E04 *INCIDENTWAVE FLUID PROPERTY 0.1424E6,1.21,0.0 *AMPLITUDE,DEFINITION=TABULAR,NAME=PRESSUREVTIME 0.0,0.0,0.0005,1.0 ************************** *BOUNDARY BEAM,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=BEAM RF, *OUTPUT,HISTORY,TIME INTERVAL=0.1E-3 *NODE OUTPUT,NSET=BEAM RF, *ENDSTEP ** **