*HEADING Use reaction forces to compute CONWEP incident and reflected pressures ** ** DEFORMABLE ELEMENT ** *NODE 1001, -25., -25., 250.0 1002, -25., 25., 250.0 1003, 25., -25., 250.0 1004, 25., 25., 250.0 2001, -25., -25., -250.0 2002, -25., 25., -250.0 2003, 25., -25., -250.0 2004, 25., 25., -250.0 ** *NODE, NSET=NSOURCE 1, 0., 0., 0. ** *NSET, NSET=NBC 1001, 1002, 1003, 1004, 2001, 2002, 2003, 2004 ** *ELEMENT, TYPE=S4RS, ELSET=PLATE 401, 1001, 1002, 1004, 1003 402, 2001, 2002, 2004, 2003 *Shell section, Elset=PLATE, Material=METAL 5.0, 3 ** *MATERIAL, NAME=METAL *DENSITY 2.78e-9, *ELASTIC 6.9e+4, 0.3 ** *BOUNDARY NBC, 1, 6 NSOURCE, 1, 6 ** *SURFACE,Type=element,NAME=Front 401, SPOS ** *SURFACE,Type=element,NAME=Back 402, SPOS ** *INCIDENT WAVE INTERACTION PROPERTY, NAME=Small_Charge, TYPE=AIR BLAST *CONWEP Charge Property ** ** Conwep units: TNT mass=ton, length=mm, time=sec, pressure=MPa ** 5.e-4, 1000. 1.e-3, 1., 1.e6 ** *STEP *DYNAMIC,EXPLICIT, DIRECT USER CONTROL 2.e-08, 2.e-04 ** *INCIDENT WAVE INTERACTION, PROPERTY=Small_Charge, CONWEP Front, NSOURCE, , 0.0, 1.0 ** *INCIDENT WAVE INTERACTION, PROPERTY=Small_Charge, CONWEP Back, NSOURCE ** *OUTPUT, FIELD, NUMBER=2 *NODE OUTPUT,NSET=NBC RT, *OUTPUT,HISTORY,FREQ=1,OP=NEW *NODE OUTPUT,NSET=NBC RT *END STEP