SPH计算总是出现cannot find element of master segment

许一个冬天

做SPH计算时，总是出现这种问题。我见了一个简单模型，就一块刚性板快速压缩一个物体
*** Error 30061 (INI+61)
     cannot find element of master segment 334 of interface 2

*** Error 30061 (INI+61)
     cannot find element of master segment 335 of interface 2

*** Error 30061 (INI+61)
     cannot find element of master segment 336 of interface 2

*** Error 30061 (INI+61)
     cannot find element of master segment 337 of interface 2

*** Error 30061 (INI+61)
     cannot find element of master segment 338 of interface 2

*** Error 30061 (INI+61)
     cannot find element of master segment 339 of interface 2

*** Error 30061 (INI+61)
     cannot find element of master segment 340 of interface 2

*** Error 30061 (INI+61)
     cannot find element of master segment 341 of interface 2

*** Error 30061 (INI+61)
     cannot find element of master segment 342 of interface 2

*** Error 30061 (INI+61)
     cannot find element of master segment 343 of interface 2

*** Error 30061 (INI+61)
     cannot find element of master segment 344 of interface 2

*** Error 30061 (INI+61)
     cannot find element of master segment 345 of interface 2

主要关键字如下：
*CONTROL_ENERGY
$#    hgen      rwen    slnten     rylen     
         2         2         2         2
*CONTROL_SHELL
$#  wrpang     esort     irnxx    istupd    theory       bwc     miter      proj
      20.0         1        -1         1         2         2         1         0
$# rotascl    intgrd    lamsht    cstyp6    tshell      
       1.0         0         0         1         0
$# psstupd   sidt4tu     cntco    itsflg    irquad      
         0         0         0         0         2
$#  nfail1    nfail4   psnfail    keepcs     delfr   drcpsid    drcprm      
         0         0         0         0         0         0       1.0
*CONTROL_SPH
$#    ncbs     boxid        dt      idim    memory      form     start      maxv
         1         01.00000E20         0       150         0       0.01.00000E15
$#    cont     deriv       ini     ishow     ierod     icont     iavis     isymp
         0         0         0         0         0         0         0       100
*CONTROL_TERMINATION
$#  endtim    endcyc     dtmin    endeng    endmas      
     0.001         0       0.0       0.0       0.0
*CONTROL_TIMESTEP
$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st
       0.0       0.9         0       0.0       0.0         0         0         0
$#  dt2msf   dt2mslc     imscl    unused    unused     rmscl     
       0.0         0         0                           0.0
*DATABASE_BNDOUT
$#      dt    binary      lcur     ioopt     
1.00000E-5         0         0         1
*DATABASE_DEFGEO
$#      dt    binary      lcur     ioopt     
1.00000E-5         0         0         1
*DATABASE_DEFORC
$#      dt    binary      lcur     ioopt     
1.00000E-5         0         0         1
*DATABASE_GCEOUT
$#      dt    binary      lcur     ioopt     
1.00000E-5         0         0         1
*DATABASE_GLSTAT
$#      dt    binary      lcur     ioopt     
1.00000E-5         0         0         1
*DATABASE_JNTFORC
$#      dt    binary      lcur     ioopt     
1.00000E-5         0         0         1
*DATABASE_MATSUM
$#      dt    binary      lcur     ioopt     
1.00000E-5         0         0         1
*DATABASE_NCFORC
$#      dt    binary      lcur     ioopt     
1.00000E-5         0         0         1
*DATABASE_RBDOUT
$#      dt    binary      lcur     ioopt     
1.00000E-5         0         0         1
*DATABASE_RCFORC
$#      dt    binary      lcur     ioopt     
1.00000E-5         0         0         1
$




*CONTACT_AUTOMATIC_NODES_TO_SURFACE_ID
$#     cid                                                                 title
         2                                                                     
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         2         3         4         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
       0.0       0.0       0.0       0.0       0.0         0       0.01.00000E20
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
       1.0       1.0       0.0       0.0       1.0       1.0       1.0       1.0
*PART
$#                                                                         title
Part          2 for Mat         2 and Elem Type         1
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         1         1         2         0         0         0         0         0
*SECTION_SOLID
$#   secid    elform       aet   
         1         1         0
*MAT_ELASTIC
$#     mid        ro         e        pr        da        db  not used        
         27.86000E-92.000000E7      0.27       0.0       0.0         0
*PART
$#                                                                         title
Part          3 for Mat         3 and Elem Type         1
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         2         1         3         0         0         0         0         0
*MAT_ELASTIC
$#     mid        ro         e        pr        da        db  not used        
         37.86000E-92.000000E7      0.27       0.0       0.0         0
*PART
$#                                                                         title
SphNode
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         3         2         1         0         0         0         0         0
*SECTION_SPH
$#   secid      cslh      hmin      hmax    sphini     death     start     
         2       1.2       0.2       2.0       0.01.00000E20       0.0

niuhon

这种一般是接触contact的问题，调整接触关系试试看，不知道楼主有没有模型图片