sph模拟液体晃动计算过程中sph粒子穿透刚性单元

mutong

本人现在做的sph方法模拟的液体晃荡——计算过程中出现sph粒子穿透sph单元：
腔体采用的是shell单元，腔体材料选择的是rigid。计算时间7s。
计算过程中1.65s以前未出现sph粒子穿透腔体单，但是在1.65s之后，1.75s时有少量sph粒子穿透腔体单元到腔体外面，1.85s时有大量粒子穿透腔体单元跑到腔体外面。图片中是计算时输出的结果，和腔体施加的运动曲线。
发生sph粒子穿透后计算时输出d3plot文时间间隔边长。未出现sph粒子穿透d3plot输出时间间隔是20分钟，出现sph粒子穿透后d3plot文件的输出时间间隔是120分钟。
k文件中的设置
1.材料：
1.1腔体
*MAT_RIGID
$#     mid        ro         e        pr         n    couple         m     alias
         2 1600.00001.2000E+12  0.300000     0.000     0.000     0.000         
$#     cmo      con1      con2   
  1.000000        6.        0.
$#lco or a1        a2        a3        v1        v2        v3  
     0.000     0.000     0.000     0.000     0.000     0.000

  1.2 水
*MAT_NULL
$#     mid        ro        pc        mu     terod     cerod        ym        pr
         1 775.00000-10.000000 8.6000E-4     0.000     0.000     0.000     0.000
*EOS_GRUNEISEN
$#   eosid         c        s1        s2        s3     gamao         a        e0
         1 1647.0000  1.921000     0.000     0.000  0.350000     0.000     0.000

    1.3 水的状态方程
*EOS_GRUNEISEN
$#   eosid         c        s1        s2        s3     gamao         a        e0
         1 1647.0000  1.921000     0.000     0.000  0.350000     0.000     0.000
$#      v0  
     0.000

2.接触设置
*CONTACT_AUTOMATIC_NODES_TO_SURFACE_ID
$#     cid                                                                 title
         1                                                                     
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         1         1         4         2         0         0         1         1
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
  0.200000  0.200000     0.000     0.000 10.000000         1     0.0001.0000E+20
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
  1.200000  1.200000     0.000     0.000  1.000000  1.000000  1.000000  1.000000
$#    soft    sofscl    lcidab    maxpar     sbopt     depth     bsort    frcfrq
         1  0.100000         0  1.025000  2.000000         2         3         1

3. 控制设置
*CONTROL_BULK_VISCOSITY
$#      q1        q2      type     btype     
  1.500000 6.0000E-2         1         0
*CONTROL_CONTACT
$#  slsfac    rwpnal    islchk    shlthk    penopt    thkchg     orien    enmass
  0.100000     0.000         1         0         1         0         1         0
$#  usrstr    usrfrc     nsbcs    interm     xpene     ssthk      ecdt   tiedprj
         0         0         0         0  4.000000         0         0         0
$#   sfric     dfric       edc       vfc        th     th_sf    pen_sf      
     0.000     0.000     0.000     0.000     0.000     0.000     0.000
$#  ignore    frceng   skiprwg    outseg   spotstp   spotdel   spothin      
         0         0         0         0         0         0     0.000
$#    isym    nserod    rwgaps    rwgdth     rwksf      icov    swradf    ithoff
         0         0         1     0.000  1.000000         0     0.000         0
$#  shledg    pstiff    ithcnt    tdcnof     ftall    unused    shltrw      
         0         0         0         0         0               0.000
*CONTROL_ENERGY
$#    hgen      rwen    slnten     rylen     
         2         2         2         2
*CONTROL_HOURGLASS
$#     ihq        qh  
         0     0.000
*CONTROL_SHELL
$#  wrpang     esort     irnxx    istupd    theory       bwc     miter      proj
20.000000         0        -1         0         2         2         1         0
$# rotascl    intgrd    lamsht    cstyp6    tshell      
  1.000000         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.000000
*CONTROL_SPH
$#    ncbs     boxid        dt      idim    memory      form     start      maxv
         1         01.0000E+20         3       800         0     0.0001.0000E+15
$#    cont     deriv       ini     ishow     ierod     icont     iavis     
         0         0         0         0         0         0         0
*CONTROL_TERMINATION
$#  endtim    endcyc     dtmin    endeng    endmas      
19.000000         0     0.000     0.000 1.0000E+8
*CONTROL_TIMESTEP
$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st
     0.000  0.900000         0     0.000     0.000         0         0         0
$#  dt2msf   dt2mslc     imscl    unused    unused     rmscl     
     0.000         0         0   

mutong

重新上传图片。

mutong

怎么没人回答，请大神回答。

mutong

sph粒子模拟液体晃动，当计算到一定时间时，突然出现粒子“爆炸”现象，出现大量sph粒子穿透腔体现象。下面两幅图分别为t=4.2s 和t=4.25s的输出结果。

zdsimwe

接触问题，shell的厚度？

欧阳飞

1、增大接触厚度
2、增大罚刚度因子
3、SPH方法本身存在拉伸不稳定性的固有缺陷，如果此缺陷不幸在此CASE里被激发，唯一能做尝试就是加密粒子（但无法根除）。

Try 1 & 2， keep watch on 3.

mutong

欧阳飞 发表于 2015-12-9 11:56
1、增大接触厚度
2、增大罚刚度因子
3、SPH方法本身存在拉伸不稳定性的固有缺陷，如果此缺陷不幸在此CASE里 ...

谢谢，是有这方面的原因

mutong

zdsimwe 发表于 2015-8-20 12:36
接触问题，shell的厚度？

谢谢啊，

小威学子

谢谢楼主，学习一下

小威学子

学习一下:):):):)

anak2006

最后修改效果怎么样？？、

小威学子

学习一下

枫叶1

你好，我也在做这个类型，在ansys里是个什么流程啊，也是油箱晃动问题，一开始建模，要画网格吗，在dyna没找到sph这个模块啊，多谢了！

枫叶1

楼主，这个求解要花费多长时间啊！