Axisymmetric cannon
At t=0, an air chamber with high temperature and high pressure gaz is released propelling a projectile at high velocity. This case uses the axisymmetric CESE solver is order to reduce the calculation times.
https://www.dynaexamples.com/cese/oricese/cese_axisymmetric
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Axisymmetric cannon
At t=0, an air chamber with high temperature and high pressure gaz is released propelling a projectile at high velocity. This case uses the axisymmetric CESE solver is order to reduce the calculation times.
At t=0, an air chamber with high temperature and high pressure gaz is released propelling a projectile at high velocity. This case uses the axisymmetric CESE solver is order to reduce the calculation times.
CESE Velocity fringes
*CESE_BOUNDARY_AXISYMMETRIC_SET *CESE_BOUNDARY_REFLECTIVE_SET *CESE_CONTROL_LIMITER *CESE_CONTROL_SOLVER *CESE_CONTROL_TIMESTEP *CESE_EOS_IDEAL_GAS *CESE_INITIAL *CESE_INITIAL_SET *CESE_PART *CONTROL_TERMINATION *DATABASE_BINARY_D3PLOT *DEFINE_CURVE *ELEMENT_SHELL *EOS_GRUNEISEN_TITLE *KEYWORD *MAT_CRUSHABLE_FOAM_TITLE *MAT_ELASTIC_PLASTIC_HYDRO_TITLE *MAT_LINEAR_ELASTIC_DISCRETE_BEAM *MAT_RIGID *NODE *PARAMETER *PART *TITLE *SECTION_SHELL *SET_NODE_LIST_TITLE *SET_PART_LIST_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 R8.0.0 (or higher) with double precision
$X
$X------------------------------------------------------------------------------
$# UNITS: (g/cm/mus)
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
CESE Axisymmetric example
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ PARAMETERS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
R T_end 5000
R dt_plot 50
$
$--- Fluid
$
R dt_fluid 0.0001
Rcfl_fluid 0.5
R T_1 293
R P_1 1e-6
R T_2 353
R P_2 3e-4
R cv 717.5e-8
R cp 1004.5e-8
R W 130
*INCLUDE
mesh.k
*INCLUDE
struc.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
200 1 101
*CESE_CONTROL_TIMESTEP
$ iddt cfl dtint
2&cfl_fluid &dt_fluid
*CESE_CONTROL_LIMITER
$ idlmt alfa beta epsr
0 0.0 1.0 0.5
$---+----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
1 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
&P_1 &T_1
*CESE_INITIAL_SET
$ setID Duic vic wic rhoic pic tic
111 &P_2 &T_2
*CESE_BOUNDARY_NON_REFLECTIVE_SET
$ ssid
5
*CESE_BOUNDARY_REFLECTIVE_SET
1
2
3
4
6
*CESE_BOUNDARY_AXISYMMETRIC_SET
7
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ DATABASE (OUTPUT) $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DATABASE_BINARY_D3PLOT
&dt_plot
*END