Command manual

Command list

VULNERABLE_AREA_ASSESSMENT

Loads
Attention: This command is in the beta stage and the format may change over time.
*VULNERABLE_AREA_ASSESSMENT
"Optional title"
coid
entype${}_a$, enid${}_a$, entype${}_v$, enid${}_v$, tid, $RHA_c$, cid${}_{RHAc}$, $v_{imp}$
$\psi_0$, $\psi_1$, $N_{\psi}$, $\theta_0$, $\theta_1$, $N_{\theta}$
cid${}_{\phi def}$, $\Delta_{map}$, $\Delta_{hit}$

Parameter definition

VariableDescription
coid Command ID
entype${}_a$ Armour entity type
options: P, PS, ALL
enid${}_a$ Armour entity identification number
entype${}_v$ Protected volume entity type
options: P, PS, ALL
enid${}_v$ Protected volume entity identification number
tid ID of table with material properties
$RHA_c$ Critical RHA equivalent armour thickness at normal impact
cid${}_{RHAc}$ ID of a CURVE defining the critical RHA equivalent armour thickness versus impact velocity
default: not used
$v_{imp}$ Impact velocity (only needed if cid${}_{RHAc}$ is defined)
default: not used
$\psi_0$ Horizontal start angle
$\psi_1$ Horizontal end angle
$N_{\psi}$ Number of horizontal frames
$\theta_0$ Vertical start angle
$\theta_1$ Vertical end angle
$N_{\theta}$ Number of vertical frames
cid${}_{\phi def}$ ID of a CURVE defining default reduction scale factor for the threat penetration capacity versus impact angle (degrees). Can be overruled by material specific curves in then material property TABLE.
default: not used
$\Delta_{map}$ Grid size for cartesian mapping of armour mesh
default: $\Delta_{map} = 0.1 \cdot RHA_c$
$\Delta_{hit}$ Hit spot grid size
default: $\Delta_{hit} = 2 \cdot \Delta_{map}$

Description

This command is primarily a tool that can be used to evaluate the ballistic protection level of armoured vehicles. It is not based on any Finite Element calculations, but uses simplified material and geometric assumptions to quickly estimate the ballistic protection at many thousands of hit spots. It is assumed that the threat follows a straight path through the armour and that that the protective performance of the different materials can be expressed in terms of specific relative RHA equivalence.

The RHA equivalent penetration capacity of the threat is either defined with the parameter $RHA_c$, or with a penetration capacity CURVE cid${}_{RHAc}$ + an impact velocity $v_{imp}$. (The threat specific curve specifies the RHA equivalent penetraton capacity versus impact velocity.)

By default $RHA_c$ defines the distance the threat can penetrate the target material before being stopped. The threat is assumed to follow a straight line and, hence, the effective armour thickness versus impact angle $\phi$ is scaled by a factor $1/\mathrm{cos}(\phi)$.


The larger impact angle, the less armour is needed to stop the threat. The default critical armour thickness versus impact angle $h(\phi)$ becomes:

$\displaystyle{h(\phi) = RHA_c \cdot \mathrm{cos}(\phi)}$

In reality the impact angle effect can be quite different, especially if the the threat is disturbed enough to ricochet or significantly change its path. This tool always assumes a straight penetration path. However, the real impact angle influence on the penetration capacity can still be taken into account. It is done by providing a CURVE cid${}_{\phi def}$ specifying the relative reduction in penetration capacity versus impact angle. This CURVE will then replace the default $\mathrm{cos}(\phi)$ relationship between impact angle and penetration capacity.

User defined impact angle dependency
User defined impact angle dependency

The command references to a TABLE listing the specific RHA equivalence (RHA equivialence per unit thickness) of the armour materials. The specific RHA equivalence of each material can be complemented with optional parameters defining edge effects, local softening in heat affected zones (see INITIAL_STATE_HAZ) and an optional material specific impact angle dependency.

The local specific RHA equivalence at a point on (or inside) the armour is obtained by averaging the material properties in a circular region perpendicular to the shooting direction. The radius of this region is determined by the impacted material's sensitivity to edge effects. Four simple examples are given below. The target material ID's are 100 (grey) and 200 (orange).

Specific RHA equivalence at a point on the armour (single material)
Specific RHA equivalence at a point on the armour (single material)
Specific RHA equivalence at a point on the armour (two materials with a gap in between)
Specific RHA equivalence at a point on the armour (two materials with a gap in between)
Specific RHA equivalence at a point on the armour (single material, partially damaged (HAZ)
Specific RHA equivalence at a point on the armour (single material, partially damaged (HAZ)

The specific RHA equivalence integrated along the threat path through the armour renders a local RHA equivalent armour thickness $RHA_l$.

Local RHA equivalent armour thickness, a simple example
Local RHA equivalent armour thickness, a simple example

This value is compared to the critical armour thickness $RHA_c$ and a damage parameter is defined as:

$\displaystyle{D = \frac{RHA_l}{RHA_c}}$

$D \geq 1$ means that the armour will be perforated at the specific location for the given firing direction.

To speed up the analysis the original Finite Element mesh is mapped to cartesian background grid. All penetration calculations are done on this background grid. The grid size is $\Delta_{map}$. The finer grid the more accurate results, but also longer calculation time.

Example

Ballistically protected teapot

A simple example of a protected teapot. The critical armour thickness is 20mm (threat dependent parameter). That is, the teapot will be damaged if it is shot at from a direction and location where the RHA equivalent armour thickness is less than 20mm.

A protected teapot
A protected teapot
*PARAMETER
RHA_c = 0.02, "Critical armour RHA equivalence 20mm"
tid = 100, "ID of table with material properties"
psi_0 = 0, "Horizontal start angle"
psi_1 = 360, "Horizontal end angle"
N_psi = 10, "Number of horizontal frames"
theta_0 = 0, "Vertical start angle"
theta_1 = 20, "Vertical end angle"
N_theta = 3, "Number of vertical frames"
#
RHAeq_10 = 1.00, "RHA equivalence (material 10)"
Redge_10 = 0.04, "Edge effect radius (material 10)"
sfedge_10 = 0.50, "Edge scale factor (material 10)"
sfgap_10 = 0.50, "Gap scale factor (material 10)"
sfdmg_10 = 1.00, "Damage scale factor (material 10)"
#
RHAeq_11 = 0.50, "RHA equivalence (material 11)"
Redge_11 = 0.04, "Edge effect radius (material 11)"
sfedge_11 = 0.50, "Edge scale factor (material 11)"
sfgap_11 = 0.50, "Gap scale factor (material 11)"
sfdmg_11 = 1.00, "Damage scale factor (material 11)"
*VULNERABLE_AREA_ASSESSMENT
"protected teapot example"
1
PS, 1014, PS, 12, [%tid], [%RHA_c]
[%psi_0], [%psi_1], [%N_psi], [%theta_0], [%theta_1], [%N_theta]
*SET_PART
"armour"
1014
10..14
*SET_PART
"teapot (protected volume)"
12
1, 2
*TABLE
"material properties"
[%tid], 7
10, [%RHAeq_10], [%Redge_10], [%sfedge_10], [%sfgap_10], [%sfdmg_10], 0
11, [%RHAeq_11], [%Redge_11], [%sfedge_11], [%sfgap_11], [%sfdmg_11], 0
#
# --- PART ---
#
*PART
"lid"
1, 1
"pot"
2, 1
"plate 1"
10, 10
"plate 2"
11, 11
"plate 3"
12, 11
"plate 4"
13, 11
"plate 5"
14, 10
*END

We choose to test the protection in 10 horizontal and 3 in vertical angles. This brings us to a total of 30 different shooting directions. In each of these directions a dense grid of hit spots will be tested. A damage value will be calculated for each one of these hit spots. $D=0$ means that the armour is infinitiely thick, $D=0.5$ corresponds to an RHA equivalent thickness that is twice as large as needed to protect the teapot and $D=1$ is the limit where the protection is perforated.

The grids of damage values are used to generate damage maps. These maps visualize the armour performance in the different firing directions.

Teapot, armour and 30 damage maps
Teapot, armour and 30 damage maps
One specific damage map projected onto the armour
One specific damage map projected onto the armour