Injection Analysis
This LS-DYNA example shows the injection of a granular material into a funnel using the discrete element method (DEM).
https://www.dynaexamples.com/dem/injection-analysis
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Injection Analysis
This LS-DYNA example shows the injection of a granular material into a funnel using the discrete element method (DEM).
*CONTROL_DISCRETE_ELEMENT *CONTROL_ENERGY *CONTROL_OUTPUT *CONTROL_TERMINATION *CONTROL_TIMESTEP *DATABASE_BINARY_D3PLOT *DATABASE_BINARY_DEMFOR *DEFINE_BOX *DEFINE_CURVE *DEFINE_DE_ACTIVE_REGION *DEFINE_DE_INJECTION *DEFINE_DE_TO_SURFACE_COUPLING *ELEMENT_SHELL *END *KEYWORD *LOAD_BODY_Z *MAT_ELASTIC_TITLE *MAT_RIGID *NODE *PART *SECTION_SOLID_TITLE *SECTION_SHELL *TITLE
$----------------------------------------------------------------------------- $ $ Example provided by Dr. Hailong Teng (LSTC) and Dr. Nils Karajan (DYNAmore) $ $ Date: Tue Mar 26 19:42:53 CET 2013 $ $ 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. $ $------------------------------------------------------------------------------- $# UNITS: ton / mm / s / N / MPa / Nmm $------------------------------------------------------------------------------- *KEYWORD *TITLE Funnel Injection $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ $ $ CONTROL CARDS $ $ $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *CONTROL_TERMINATION $# endtim endcyc dtmin endeng endmas 3.0 *CONTROL_TIMESTEP $# dtinit tssfac isdo tslimt dt2ms lctm erode ms1st 2 $# dt2msf dt2mslc imscl unused unused rmscl *DEFINE_CURVE $# lcid sidr sfa sfo offa offo dattyp lcint 2 $# a1 o1 1.000000e-004 1000.0 1.000000e-004 *CONTROL_ENERGY $# hgen rwen slnten rylen 1 2 2 1 *CONTROL_OUTPUT $# npopt neecho nrefup iaccop opifs ipnint ikedit iflush 1 $# iprtf ierode tet10 msgmax ipcurv gmdt ip1dblt eocs 2 50 *DATABASE_BINARY_D3PLOT $# dt lcdt beam npltc psetid 0.02 *DATABASE_BINARY_DEMFOR $# dt lcdt beam npltc psetid 0.02 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ $ $ CONTACT DEFINITION $ $ $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *CONTROL_DISCRETE_ELEMENT $# ndamp tdamp fric fricr normk sheark cap mxnsc 0.5 0.5 0.57 0.1 *DEFINE_DE_ACTIVE_REGION $# id type xm ym zm 1 1 *DEFINE_BOX $# boxid xmn xmx ymn ymx zmn zmx 1 -450.0 450.0 -350.0 350.0 -300.0 700.0 $ $--- $ *DEFINE_DE_TO_SURFACE_COUPLING $# slave master stype mtype 1 1 1 $# frics fricd damp bsort lcvx lcvy lcvz wearc 0.3 0.3 0.8 10 1.5 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ $ $ BOUNDARY/INITIAL CONDITIONS $ $ $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *DEFINE_DE_INJECTION $# pid sid xc yc zc xl yl cid 2 1 0.0 0.0 600.0 300.0 400.0 $# rmass rmin rmax vx vy vz tbeg tend 0.0015 2.0 3.0 2.0 *LOAD_BODY_Z $# lcid sf lciddr xc yc zc cid 1 9810.0 *DEFINE_CURVE $# lcid sidr sfa sfo offa offo dattyp lcint 1 $# a1 o1 1.0 1000.0 1.0 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ $ $ PARTS/ SECTION/ MATERIAL $ $ $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *PART Steel funnel and plate $# pid secid mid eosid hgid grav adpopt tmid 1 1 2 *SECTION_SHELL $# secid elform shrf nip propt qr/irid icomp setyp 1 $# t1 t2 t3 t4 nloc marea idof edgset 1.0 *MAT_RIGID $# mid ro e pr n couple m alias 2 7.85E-9 2.10e+5 0.3 $# cmo con1 con2 1.0 7 7 $#lco or a1 a2 a3 v1 v2 v3 $ $--- $ *PART DES $# pid secid mid eosid hgid grav adpopt tmid 2 2 1 *SECTION_SOLID_TITLE DES $# secid elform aet 2 *MAT_ELASTIC_TITLE DES $# mid ro e pr da db not used 1 1.390E-9 0.70e+5 0.3
