Bar II — Welcome to LS-DYNA Examples

Home » Knowledge Base » Mesh-Free Methods » Bar II

This example shows the impact of two equivalent bars made of plastic material against a planar rigid wall. The upper bar is calculated using Lagrangian discretization, the lower one is calculated using SPH method.

Keywords

*CONTROL_CONTACT
*CONTROL_CPU
*CONTROL_DAMPING
*CONTROL_ENERGY
*CONTROL_OUTPUT
*CONTROL_SHELL
*CONTROL_SPH
*CONTROL_TERMINATION
*CONTROL_TIMESTEP
*DATABASE_BINARY_D3PLOT
*DATABASE_EXTENT_BINARY
*DATABASE_GLSTAT
*DATABASE_MATSUM
*DATABASE_RCFORC
*DATABASE_RWFORC
*DATABASE_SLEOUT
*ELEMENT_SHELL
*ELEMENT_SPH
*END
*INITIAL_VELOCITY
*KEYWORD
*MAT_PLASTIC_KINEMATIC
*NODE
*PART
*RIGIDWALL_PLANAR
*SECTION_SHELL
*SECTION_SPH
*SET_NODE_LIST_GENERATE
*TITLE

Reduced Input

  *KEYWORD
  $* LS-DYNA Keyword file: LSTC LS-PRE/POST Ver 1.0(Beta)- 31OCT2002(1553)
  $* Total Number of Nodes  =     6081,Largest node ID  = 1004800
  $* Total Number of Shells =     1200,Largest Shell ID = 1200
  $* Total Number of Parts  =        2,Largest Part ID  = 101
  $* Total Number of Mats   =        2 Largest Mat ID   = 2
  *TITLE
  SPH / TP1 bis Impact barre 2D contre SW                                         
  *CONTROL_CONTACT
  $*  slsfac    rwpnal    islchk    shlthk    penopt    thkchg     orien    enmass
         0.0       0.0
  $*  usrstr    usrfrc     nsbcs    interm     xpene     ssthk      ecdt   tiedprj
           0         0         0         0       0.0
  *CONTROL_CPU
  $*  cputim    
         0.0
  *CONTROL_DAMPING
  $*  nrcyck     drtol    drfctr    drterm    tssfdr    irelal     edttl    idrflg
           0       0.0       0.0       0.0       0.0         0       0.0
  *CONTROL_ENERGY
  $*    hgen      rwen    slnten     rylen     
           2         2         2         1
  *CONTROL_OUTPUT
  $*   npopt    neecho    nrefup    iaccop     opifs    ipnint    ikedit    iflush
           0         0         0         0       0.0
  *CONTROL_SHELL
  $*  wrpang     esort     irnxx    istupd    theory       bwc     miter      proj
         0.0
  *CONTROL_SPH
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$   *CONTROL_SPH  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $*    ncbs     boxid        dt       dim    memory      form     start      maxv
           0         0       0.0         2         0         0       0.0       0.0
  *CONTROL_TERMINATION
  $*  endtim    endcyc     dtmin    endeng    endmas    
    0.035000         0       0.0       0.0       0.0
  *CONTROL_TIMESTEP
  $*  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st
         0.0  0.900000         0       0.0       0.0
  *DATABASE_GLSTAT
  4.0000e-04
  *DATABASE_MATSUM
  4.0000e-04
  *DATABASE_RCFORC
  4.0000e-04
  *DATABASE_RWFORC
  4.0000e-04
  *DATABASE_SLEOUT
  4.0000e-04
  *DATABASE_BINARY_D3PLOT
    0.004000
  *DATABASE_EXTENT_BINARY
           0         0         0         0         0         0         0         0
           0         0         0         0         0         0         0
  *ELEMENT_SHELL
  $
  *ELEMENT_SPH
  $
  *NODE
  $
  *INITIAL_VELOCITY
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $$$$     INITIAL CONDITION        
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
           0         0         0
   -50.00000       0.0       0.0       0.0       0.0       0.0
  *PART
  $* title                                                                        
                                                                                  
  $*     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
           1         1         1
  *SECTION_SHELL
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ *LAGRANGIAN PART $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $*   secid    elform      shrf       nip     propt   qr/irid     icomp     setyp
           1        13       0.0         0         0       0.0
  $*      t1        t2        t3        t4      nloc     marea     
    1.000000  1.000000  1.000000  1.000000         0       0.0
  *MAT_PLASTIC_KINEMATIC
  $*     mid        ro         e        pr      sigy      etan      beta      
           1  0.0027106.8900e+04  0.330000  117.0000  6.800000  1.000000
  $*     src       srp        fs        vp        
         0.0       0.0       0.0       0.0
  *PART
  $* title                                                                        
                                                                                  
  $*     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         101       101         1
  *SECTION_SPH
  $
  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$  *SPH PART $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
  $
  $*   secid      cslh      hmin      hmax    sphini     death     start     
         101  1.300000       0.0       0.0       0.0       0.0       0.0
  *MAT_PLASTIC_KINEMATIC
  $*     mid        ro         e        pr      sigy      etan      beta      
           2  0.0027106.8900e+04  0.330000  117.0000  6.800000  1.000000
  $*     src       srp        fs        vp        
         0.0       0.0       0.0       0.0
  *RIGIDWALL_PLANAR
           0         0         0       0.0
         0.0       0.0       0.0  1.000000       0.0       0.0       0.0       0.0
  *SET_NODE_LIST_GENERATE
         101       0.0       0.0       0.0       0.0
     1000001   1004800
  *END

Figures

image01.gif

image02.gif

Animated Result

Result

Download