注入分析

LS-DYNA实例显示用离散元法向漏斗注入粒子

*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.海龙腾博士18luck快乐彩Nils Karajan(DYNAMERE)元日期:Tue Mar2619:42:53CET2013$E-Mail:info@dyamore.de$We$ $------------------------------------------------------------------------------- $# 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