OUTPUT_DAMAGE_MAP

Output

Beta command

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

Variable
Description
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