Blast wave
It is possible to combine the keyword LOAD_BLAST_ENHANCED with the CESE solver in order to achieve an accurate tracking of the pressure loads and pressure wave propagation due to the detonation of a conventional explosive.
https://www.dynaexamples.com/cese/oricese/cese_blast
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Blast wave
It is possible to combine the keyword LOAD_BLAST_ENHANCED with the CESE solver in order to achieve an accurate tracking of the pressure loads and pressure wave propagation due to the detonation of a conventional explosive.
It is possible to combine the keyword LOAD_BLAST_ENHANCED with the CESE solver in order to achieve an accurate tracking of the pressure loads and pressure wave propagation due to the detonation of a conventional explosive.
CESE pressure fringes
*CESE_BOUNDARY_BLAST_LOAD_SET *CESE_BOUNDARY_REFLECTIVE_SET *CESE_DATABASE_ELOUT *CESE_CONTROL_LIMITER *CESE_CONTROL_SOLVER *CESE_CONTROL_TIMESTEP *CESE_EOS_IDEAL_GAS *CESE_INITIAL *CESE_PART *CONTROL_TERMINATION *DATABASE_BINARY_D3PLOT *ELEMENT_SOLID *KEYWORD *LOAD_BLAST_ENHANCED *NODE *PARAMETER *PARAMETER_EXPRESSION *TITLE *SET_SEGMENT
$-----------------------------------------------------------------------------
$
$ Example provided by Iñaki (LSTC)
$
$ E-Mail: info@dynamore.de
$ Web: http://www.dynamore.de
$
$ Copyright, 2015 DYNAmore GmbH
$ Copying for non-commercial usage allowed if
$ copy bears this notice completely.
$
$X------------------------------------------------------------------------------
$X
$X 1. Run file as is.
$X Requires LS-DYNA MPP R9.0.0 (or higher) with double precision
$X
$X------------------------------------------------------------------------------
$# UNITS: (kg/m/s)
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
CESE Load Blast example
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ PARAMETERS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
R T_end 5e-5
R dt_plot 1e-6
$
$--- Fluid
$
R dt_fluid 0.0001
Rcfl_fluid 0.5
R ro_1 1.00
R P_1 1.000e+5
R cv 717.5
R cp 1004.5
R W 130
*INCLUDE
mesh.k
*INCLUDE
load_blast.k
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ CESE CONTROL CARDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CONTROL_TERMINATION
$ endtim endcyc dtmin endeng endmas
&T_end
*CESE_CONTROL_SOLVER
$ iframe iflow igeom
0 1 3
*CESE_CONTROL_TIMESTEP
$ iddt cfl dtint
2&cfl_fluid &dt_fluid
*CESE_CONTROL_LIMITER
$# idlmt alfa beta epsr
2 10.0 1.0 .1
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ CESE PARTS/ EOS/ MATERIAL $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CESE_PART
$ pid mid eosid
5 3
*CESE_EOS_IDEAL_GAS
$ eosid cv cp
3 &Cv &Cp
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ CESE BOUNDARY/INITIAL CONDITIONS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CESE_INITIAL
$# u v w rho p t
0.0 0.0 0.0 1.29 100000.0 0.0
*CESE_BOUNDARY_REFLECTIVE_SET
$# ssid
3
*CESE_BOUNDARY_BLAST_LOAD_SET
$# bid ssid
1 1
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ DATABASE (OUTPUT) $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DATABASE_BINARY_D3PLOT
&dt_plot
*CESE_DATABASE_ELOUT
$# outlv dtout
1 0.0
$# elsid
111
*END