/BEM/DAA
Block Format Keyword Doubly Asymptotic Approximation for Underwater Explosion, where the fluid mass matrix is computed by boundary element method.
Format
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
/BEM/DAA/daa_ID/unit_ID | |||||||||
daa_title | |||||||||
surf_ID | grav_ID | ||||||||
C | |||||||||
Xs | Ys | Zs | |||||||
Iform | Ipri | Ipres | Kform | Freesurf | Afterflow | Integr |
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
Pm |
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
fct_IDP | FscaleP |
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
Xc | Yc | Zc |
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
XA | YA | ZA | |||||||
Dir-X | Dir-Y | Dir-Z |
Definitions
Field | Contents | SI Unit Example |
---|---|---|
daa_ID | DAA block
identifier (Integer, maximum 10 digits) |
|
unit_ID | Unit Identifier (Integer, maximum 10 digits) |
|
daa_title | DAA block
title (Character, maximum 100 characters) |
|
surf_ID | Wet surface identifier
2
3 (Integer) |
|
grav_ID | /GRAV
option identifier (Integer) |
|
Fluid
density (Real) |
||
C | Fluid sound
speed (Real) |
|
Xs | X-coordinate of standoff
point 3 (Real) |
|
Ys | Y-coordinate of standoff
point 3 (Real) |
|
Zs | Z-coordinate of standoff
point 3 (Real) |
|
Iform | BEM solution flag
(Integer) |
|
Ipri | Printout flag level
|
|
Ipres | Pressure loading flag
6
(Integer) |
|
Kform | Analysis flag
(Integer) |
|
Freesurf | Free surface flag 6
(Integer) |
|
Afterflow | Afterflow computation
7
|
|
Integr | Time integer flag
|
|
Pm | Maximum pressure 5 (Real) |
|
Decay
time (Real) |
||
fct_IDP | Incident pressure function
identifier (Integer) |
|
FscaleP | Ordinate (pressure) scale
factor for
fct_IDP (Real) |
|
XC | X-coordinate of explosive
charge (Real) |
|
YC | Y-coordinate of explosive
charge (Real) |
|
ZC | Z-coordinate of explosive
charge (Real) |
|
XA | X-coordinate of point A on
the free surface (Real) |
|
YA | Y-coordinate of point A on
the free surface (Real) |
|
ZA | Z-coordinate of point A on
the free surface (Real) |
|
Dir-X | X-component of the normal
to the free surface plane (Real) |
|
Dir-Y | Y-component of the normal
to the free surface plane (Real) |
|
Dir-Z | Z-component of the normal
to the free surface plane (Real) |
Comments
- The entire structure must be modeled. Symmetric analysis is not supported.
- The surface normal should be pointed into the fluid.
- Standoff point defined with (Xs, Ys, Zs) is the location where the incident pressure wave is given at time t=0:
- A plane wave can be simulated using a spherical wave and putting the explosive charge far enough away.
- Pressure at the standoff point as
a function of time is:
(1) Where,- Maximum pressure
- Time
- Decay time
The maximum pressure and decay time can be calculated using:(2) (3) - Explosive mass
- Distance to the explosion
- , α, and
- Constants depending on the explosive
If in kg, in meter, in MPa and in ms.α TNT 52.12 1.180 0.0895 -0.185 PETN 56.21 1.194 0.0860 -0.257 HBX 53.51 1.144 0.0920 -0.247 - A free surface is a plane defined by a point and its normal vector.
- The afterflow normal velocity is
calculated as:
(4) - P
- Fluid point
- C
- Explosive charge point
- S
- Standoff point