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

1.JPG

2.JPG


炸药和薄板均是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
(3条)
默认 最新
这个你板子材料是有失效的,你可以采用一号单元试下,或者24号单元
评论 3 点赞 1
回复
多谢,在下初涉SPH,不明白1号单元,24号单元是指什么??
评论 1 点赞
回复
我说错了,是板料的材料可以采用不失效的单元,jc是带有失效的
评论 1 点赞
查看其它1条回复 >
如问题解决了,请“采纳答案”,老师能获得额外的技术分
评论 点赞

查看更多评论 >

点赞 评论 6 收藏 2
关注