SPH接触爆炸，粒子飞散讨论

zqpcgm

炸药和薄板均是SPH建模，在爆炸过程中，薄板应该向下弯曲，但是在模拟时出现如图所示的破裂，粒子快速飞散，炸药粒子甚至直接穿透薄板飞出。炸药和薄板参数均选自AUTODYN模型库。后来我发现，建模后，将薄板向下移动一个粒子尺寸的距离，这种现象就消失了。但这不是真实的原因，因为在autodyn里接触建模是没有任何问题的。请教造成这种现象的真实原因是什么？应该如何改进？附上我的K文件，单位cm-g-us
模型是二维的，Shell填充粒子生成的。
文件太大，上传不了，我就把关键字段附上，删去了NODE 和ELEMENT SPH两部分的信息。
*CONTROL_BULK_VISCOSITY
$#      q1        q2      type     btype     
       1.5      0.06         1         0
*CONTROL_ENERGY
$#    hgen      rwen    slnten     rylen     
         2         2         2         2
*CONTROL_SPH
$#    ncbs     boxid        dt      idim    memory      form     start      maxv
         0         01.00000E20         2       150         5       0.00.69999999
$#    cont     deriv       ini     ishow     ierod     icont     iavis     
         0         0         0         0         0         0         0
*CONTROL_TERMINATION
$#  endtim    endcyc     dtmin    endeng    endmas      
      20.0         0       0.0       0.01.000000E8
*CONTROL_TIMESTEP
$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st
       0.00.67000002         0       0.0       0.0         0         0         0
$#  dt2msf   dt2mslc     imscl    unused    unused     rmscl     
       0.0         0         0                           0.0
*DATABASE_GLSTAT
$#      dt    binary      lcur     ioopt     
       0.2         0         0         1
*DATABASE_MATSUM
$#      dt    binary      lcur     ioopt     
       0.2         0         0         1
*DATABASE_SPHOUT
$#      dt    binary      lcur     ioopt     
       0.2         0         0         1
*DATABASE_BINARY_D3PLOT
$#      dt      lcdt      beam     npltc    psetid      
      0.02         0         0         0         0
$#   ioopt     
         0
*DATABASE_EXTENT_BINARY
$#   neiph     neips    maxint    strflg    sigflg    epsflg    rltflg    engflg
         0         0         3         1         1         1         1         1
$#  cmpflg    ieverp    beamip     dcomp      shge     stssz    n3thdt   ialemat
         0         0         4         1         1         1         2         1
$# nintsld   pkp_sen      sclp     hydro     msscl     therm    intout    nodout
         0         0       1.0         0         0         0STRESS              
$#    dtdt    resplt     neipb     
         0         0         0
*PART
$#                                                                         title
SphNode
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         3         2         3         3         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
*MAT_HIGH_EXPLOSIVE_BURN
$#     mid        ro         d       pcj      beta         k         g      sigy
         30.930999990.41600001    0.0515       0.0       0.0       0.0       0.0
*EOS_JWL
$#   eosid         a         b        r1        r2      omeg        e0        vo
         3    0.4946   0.01891 3.9070001     1.1180.33333334   0.02484       1.0
*PART
$#                                                                         title
SphNode
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         4         2         4         4         0         0         0         0
*MAT_JOHNSON_COOK
$#     mid        ro         g         e        pr       dtf        vp    rateop
         4 8.93999960.46000001       0.0       0.0       0.0       0.0       0.0
$#       a         b         n         c         m        tm        tr      epso
9.00000E-4   0.00292      0.31       0.0      1.09    1356.0     293.0       0.0
$#      cp        pc     spall        it        d1        d2        d3        d4
4.60000E-6      -9.0       3.0       0.00.54000002 4.8899999     -3.03     0.014
$#      d5      c2/p      erod     efmin     
      1.12       0.0         01.00000E-6
*EOS_GRUNEISEN
$#   eosid         c        s1        s2        s3     gamao         a        e0
         4     0.391      1.51       0.0       0.0      1.97       0.0       0.0
$#      v0  
   0.11206
*INITIAL_DETONATION
$#     pid         x         y         z        lt  
         3       0.0       1.1       0.0       0.0

*ELEMENT_SPH
$#   nid     pid            mass   
   18313       3   3.723993e-004
   18314       3   3.723999e-004

   35806       4   8.939996e-004
   35807       4   8.939996e-004
   35808       4   8.939999e-004
   35809       4   8.939999e-004
   35810       4   8.939999e-004
   35811       4   8.939999e-004
   35812       4   8.939999e-004

guojunhang

这个原因也很简单，楼主不是已经用自己的实践摸索出来了吗？一个sph粒子在空间占据的空间由它的半径决定，那么它的实际外表面离粒子的中心是有间距的，等于粒子半径。所以把实体离散为粒子时应该适当将其外边界缩小一个粒子半径那么大再划分网格。只不过你这里把板材往外移动了一些。

liangsf204

炸药爆炸的威力足以破坏底部挡板，建议对地板做约束

xiaocui_1234

请问炸药sph粒子与薄板之间是否设置接触或者耦合以达到荷载的传递？☺

xiaocui_1234

是不是设置了可以穿透？

guojunhang

就是这个示意图啦。这些粒子之间要保留合适的距离，

hanxiaoxiao

请问楼主，不需要接触么

P-wenson

hanxiaoxiao 发表于 2016-8-2 14:36
请问楼主，不需要接触么

sph粒子与sph粒子间不需要接触

hawaha1113

厉害了，往外移这都能想出来

工大傻小子

楼主您好，小弟非常想学习sph在dyna中实现爆破的模拟，可是小弟愚钝，对sph技术一无所知，恳请楼主能推荐我一些关于这方面的资料，不胜感激