Learning Forum

[kavya_p230097ce]

I am simulating the impact of discrete element (DEM) particles on a wire mesh structure in LS-DYNA. The particles are generated using *DEFINE_DE_INJECTION. The structure consists of a frame modeled with solid elements (*MAT_ELASTIC) and wires modeled with beam elements (*MAT_PLASTIC_KINEMATIC). The keyword file was initially generated using ANSYS LS-DYNA (where meshing was also done) and finalized in LS-PrePost.

For structural contacts, I defined *CONTACT_AUTOMATIC_GENERAL and *CONTACT_AUTOMATIC_NODES_TO_SURFACE. For the DEM-to-structure interaction, I implemented *DEFINE_DE_TO_BEAM_COUPLING and *DEFINE_DE_TO_SURFACE_COUPLING.

However, during the simulation, the DEM particles completely penetrate both the wires and the frame as if the structure does not exist. Are there any missing keywords or specific parameter adjustments required to properly enforce the physical interaction between the DEM particles and the structure?

[erko]

In principal yes, this can be done in theory.

Try using only the beam_coupling keyword only.

Thank you

A sample below -

see that 24 is used here for sid in both cards (23 are nr of beams). PART set id 55 contains all beams.

*DEFINE_DE_TO_BEAM_COUPLING
$#   slave    master     stype     mtype     
        24        55         0         0
$#   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
        24        24       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