OUTPUT_DAMAGE_MAP

Output
Attention: This command is in the beta stage and the format may change over time.
*OUTPUT_DAMAGE_MAP
"Optional title"
coid
pid, N${}_x$, N${}_y$, N${}_z$, csysid
$x_0$, $y_0$, $z_0$, $x_1$, $y_1$, $z_1$
Parameter definition
VariableDescription
coid Command ID
pid ID of part from which data will be mapped at termination time
N${}_x$ Number of data points in local x-direction
N${}_y$ Number of data points in local y-direction
N${}_z$ Number of data points in local z-direction
csysid Optional local coordinate system ID
default: global coordinates are used
$x_0$ Mapping domain minimum x-coordinate
$y_0$ Mapping domain minimum y-coordinate
$z_0$ Mapping domain minimum z-coordinate
$x_1$ Mapping domain maximum x-coordinate
$y_1$ Mapping domain maximum y-coordinate
$z_1$ Mapping domain maximum z-coordinate
Description

This command outputs the damage distribution inside a user defined domain.

The command instructs the solver to generate the file damage_map1.k at termination time which stores the damage values in a cartesian map of $n_1 \times n_2 \times n_3$ data points, see MAP. The stored data can later be used as input data in subsequent simulations by including the file and map the damage (and erosion) to the new part with the command INITIAL_DAMAGE_MAP.

Example
Local model of fragment impacting concrete slab (step 1)

A simple model of a fragment impacting a concrete slab. The damage (and erosion) is mapped from the local model (step 1) to a global model (step 2). The global model can be found under examples in INITIAL_DAMAGE_MAP.

Domain where damage is automatically sampled (and stored with *MAP command) at termination time. Negative damage values indicates empty space.

*UNIT_SYSTEM
SI
#
# --- TIME and OUTPUT ---
#
*TIME
0.0003
*OUTPUT
0.0001
#
# --- MESH ---
#
*COMPONENT_BOX
1, 1, 80, 80, 20
-0.4, -0.4, -0.2, 0.4, 0.4, 0.0
*COMPONENT_SPHERE
2, 2, 4
0, 0, 0.016, 0.015
*CHANGE_P-ORDER
ALL, 0, 3, 123
*GEOMETRY_PIPE
123
0, 0, -0.2, 0, 0, 0.2, 0.1
*SMOOTH_MESH
ALL, 0, 45.0
*REDISTRIBUTE_MESH_CARTESIAN
1
P, 1, XYZ, 0, 1
0, 0, 0
#
# --- MATERIAL ---
#
*MAT_HJC_CONCRETE
"concrete"
1, 2310.0, 14.86e9
0.79, 1.60, 0.61, 0.007, 50.5e6, 4.0e6, 1.0, 0.01
7.0, 16.0e6, 0.001, 800.0e6, 0.10, 0.04, 1.0, 85.0e9
-171.0e9, 208.0e9, 0
*MAT_RIGID
2, 7800.0
#
# --- PART ---
#
*PART
"concrete block"
1, 1, 0, 0, 0, 5.0e-8, 2.0, 2.0
"impactor"
2, 2
#
# --- INITIAL VELOCITY ---
#
*INITIAL_VELOCITY
P, 2, 0, 0, -800.0
#
# --- CONTACT ---
#
*CONTACT
1
ALL, 0, ALL, 0
#
# --- DAMAGE MAP ---
#
*OUTPUT_DAMAGE_MAP
"test 1"
1
1, 40, 40, 40
-0.1, -0.1, -0.2, 0.1, 0.1, 0.0
*END