Noah
Subscriber

Hello,

I am trying to simulate a simple crash test where simple geometry is deformed using a drop weight.
I am using a shell plate as a drop weight.
See Keywords.

LS-DYNA keyword deck by LS-PrePost
*CONTROL_ACCURACY
$#     osu       inn    pidosu      iacc     exacc     
         1         4         0         0          
*CONTROL_BULK_VISCOSITY
$#      q1        q2      type     btype    tstype      
       1.5      0.06        -2         0         0
*CONTROL_ENERGY
$#    hgen      rwen    slnten     rylen     irgen     maten     drlen     disen
         2         2         2         2         2                              
*CONTROL_SHELL
$#  wrpang     esort     irnxx    istupd    theory       bwc     miter      proj
      20.0         1        -1         4        16         2         1         0
$# rotascl    intgrd    lamsht    cstyp6    thshel      
       1.0         0         0         1         0
$# psstupd   sidt4tu     cntco    itsflg    irquad    w-mode   stretch      icrq
                   0         0         0         3                             0
$#  nfail1    nfail4   psnfail    keepcs     delfr   drcpsid    drcprm   intperr
                             0         0         0         0       1.0         0
*CONTROL_TERMINATION
$#  endtim    endcyc     dtmin    endeng    endmas     nosol     
      0.03         0       0.0       0.01.000000E8         0
*CONTROL_TIMESTEP
$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st
       0.0       0.0         0       0.0       0.0         0         0         0
$#  dt2msf   dt2mslc     imscl    unused    unused     rmscl     emscl      ihdo
       0.0         0         0                           0.0                   0
*DATABASE_NODOUT
$#      dt    binary      lcur     ioopt   option1   option2       
     0.001         0         0         1       0.0         0
*DATABASE_RWFORC
$#      dt    binary      lcur     ioopt     
     0.001         0         0         1
*DATABASE_BINARY_D3PLOT
$#      dt      lcdt      beam     npltc    psetid      
2.00000E-4         0         0         0         0
$#   ioopt      rate    cutoff    window      type      pset    
         0       0.0       0.0       0.0         0         0
*CONTACT_FORMING_SURFACE_TO_SURFACE_ID
$#     cid                                                                 title
         2mass
$#   surfa     surfb  surfatyp  surfbtyp   saboxid   sbboxid      sapr      sbpr
        11         8         3         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
       0.2      0.18       0.0                20.0         0       0.01.00000E20
$#    sfsa      sfsb      sast      sbst     sfsat     sfsbt       fsf       vsf
       1.0       1.0       0.0       0.0       1.0       1.0       1.0       1.0
*CONTACT_AUTOMATIC_SINGLE_SURFACE_ID
$#     cid                                                                 title
         3Selfcontact
$#   surfa     surfb  surfatyp  surfbtyp   saboxid   sbboxid      sapr      sbpr
        11                   3         0                             0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
       0.2       0.2       0.0       0.0       0.0                 0.01.00000E20
$#    sfsa      sfsb      sast      sbst     sfsat     sfsbt       fsf       vsf
       1.5       1.0                           1.0       1.0       1.0       1.0
$#    soft    sofscl    lcidab    maxpar     sbopt     depth     bsort    frcfrq
         1       0.1         0     1.025       1.0         2         0         1
$#  penmax    thkopt    shlthk     snlog      isym     i2d3d    sldthk    sldstf
       0.0         1         2         0         0         0       0.0       0.0
$#    igap    ignore    dprfac    dtstif    unused    unused    flangl   cid_rcf
         1         0       0.0       0.0                           0.0         0
*MAT_MODIFIED_JOHNSON_COOK_TITLE
MAT_107_V2A
$#     mid        ro         e        pr      beta       xsi        cp     alpha
         17.90000E-9  200000.0       0.3       1.0       0.95.000000E81.60000E-5
$#   e0dot        tr        tm        t0     flag1     flag2     
       1.0     293.0    1673.0     293.0       0.0       1.0
$#  a/siga       b/b   n/beta0   c/beta1      m/na    
     310.0    1000.0      0.65     0.007       1.0
$#    q1/a      c1/n q2/alpha0 c2/alpha1         
       0.0       0.0       0.0       0.0
$#   dc/dc     pd/wc     d1/na     d2/na     d3/na     d4/na     d5/na     
       1.0     150.0       0.0       0.0       0.0       0.0       0.0
$#      tc      tauc    
       0.0     850.0
*MAT_RIGID_TITLE
Fallgewicht
$#     mid        ro         e        pr         n    couple         m     alias
         22.00000E-5  200000.0       0.3       0.0       0.0       0.0          
$#     cmo      con1      con2    
       1.0         4         7
$#lco or a1        a2        a3        v1        v2        v3  
       0.0       0.0       0.0       0.0       0.0       0.0
*PART
$#                                                                         title
mass
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         8         2         2         0         0         0         0         0
*SECTION_SHELL_TITLE
mass
$#   secid    elform      shrf       nip     propt   qr/irid     icomp     setyp
         2         2    0.8333         2       1.0         0         0         1
$#      t1        t2        t3        t4      nloc     marea      idof    edgset
       1.0       1.0       1.0       1.0       0.0       0.0       0.0         0
*PART
$#                                                                         title
cylindershell
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
        11         1         1         0         1         0         0         0
*SECTION_SHELL_TITLE
Shell_2
$#   secid    elform      shrf       nip     propt   qr/irid     icomp     setyp
         1        16    0.8333         5       1.0         0         0         1
$#      t1        t2        t3        t4      nloc     marea      idof    edgset
       2.0       2.0       2.0       2.0       0.0       0.0       0.0         0
*HOURGLASS
$#    hgid       ihq        qm       ibq        q1        q2    qb/vdc        qw
         1         8      0.03         0       1.5      0.06      0.03      0.03

 

Unfortunately, this error always occurs and I do not understand what the problem is.