DETONATION

Material properties

Beta command

This command is in the beta stage and the format may change over time.

*DETONATION
"Optional title"
coid
entype, enid, $x_d$, $y_d$, $z_d$, $t_d$, $R$, $\phi$
gid

Parameter definition

Variable
Description
coid
Command ID (required)
entype
Entity type
options: P, PS
enid
Entity identification number
$x_d$, $y_d$, $z_d$
Detonation point or function (fcn)
quantity: Length
$t_d$
Detonation time or function (fcn)
quantity: Time
$R$
Optional detonation radius
quantity: Length
$\phi$
Optional curve taking delay parameter
gid
ID of detonation geometry
default: not used

Description

Detonation point definition, needed by materials such as MAT_EXPLOSIVE_JWL. The detonation radius $R$ limits the distance the detonation front is allowed to propagate through programmed burn. $\phi$ is a delay parameter (unit time). It slows the detonation front when navigating corners (such as wave shapers). The $\phi$ value quantifies the time lost in a full 90-degree turn.

A detonation function (fcn) may be used in place of the fixed time ($t_d$). The event is triggered when returning a positive value, as illustrated in the example below.

If specified, gid causes all elements enclosed by the referenced geometry to detonate simulataneously, replacing initiation at the detonation point.

Example

Detonation functions

Example where the detonation point and time are not defined a priori, but with functions (fcn). In this specific case, detonation is assumed to take place when the material pressure exceeds $1GPa$

*UNIT_SYSTEM SI # # --- TIME --- # *TIME 5.0e-5 # # --- MESH --- # *COMPONENT_BOX "C4" 1, 1, 10, 10, 10 -0.05, -0.05, -0.1, 0.05, 0.05, 0.0 *COMPONENT_SPHERE "Impactor" 2, 2, 6 -0.03, -0.02, 0.04, 0.03 # # --- MATERIAL --- # *MAT_OBJECT "C4" "{34b95a1a-c7fa-4c63-a8c7-4c4f729063ef}", 1 *MAT_RIGID "Steel" 2, 7800.0 # # --- PART --- # *PART "C4" 1, 1 "Impactor" 2, 2 # # --- DETONATION --- # *DETONATION 1 P, 1, fcn(1), fcn(2), fcn(3), fcn(4) *FUNCTION "xd" 1 pres_part(1,1) *FUNCTION "yd" 2 pres_part(1,2) *FUNCTION "zd" 3 pres_part(1,3) *FUNCTION "td" 4 pres_part(1,0) - 1.0e9 # # --- INITIAL VELOCITY --- # *INITIAL_VELOCITY "Impactor" P, 2, 0, 0, -500.0 # # --- CONTACT --- # *CONTACT 1 ALL, 0, ALL, 0, 0.0, 1.0e14 *END