有偿求助，关于SPHl粒子以及Solid单元穿透的问题？ 135

不同part之间设置相同本构，相互穿透的问题。

工况描述：不同part之间采用完全相同的本构、材料参数以及失效准则，但是相互穿透。具体情况如下：

图1 Solid之间接触穿透

图2 失效转化后SPH与Solid接触穿透

出现的问题：

1. SPH以及solid之间的穿透问题。

2. SPH粒子穿越对称约束面。

补充：

1. 这套材料参数应该没有问题，用相同参数跑通了很多模型。

2. 尝试过更换网格尺寸、接触参数以及罚函数因子等方法，均无法解决。

3. 有偿求助。

部分关键字：

*CONTACT_ERODING_SURFACE_TO_SURFACE_ID

$#     cid                                                                 title

        24ceraJIAO-PU

$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr

         5         6         3         3         1         1         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

$#    isym    erosop      iadj    

         1         1         1

$#    soft    sofscl    lcidab    maxpar     sbopt     depth     bsort    frcfrq

         1       0.1         0     1.025       2.0         1         0         1

*CONTACT_ERODING_NODES_TO_SURFACE_ID

$#     cid                                                                 title

        25ceraZHONGXINsph-ceraBIAN1

$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr

         9         3         3         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      10.0       0.0       0.0       1.0       1.0       1.0       1.0

$#    isym    erosop      iadj    

         1         1         1

$#    soft    sofscl    lcidab    maxpar     sbopt     depth     bsort    frcfrq

         1       0.1         0     1.025       2.0         2         0         1

*CONTROL_CONTACT

$#  slsfac    rwpnal    islchk    shlthk    penopt    thkchg     orien    enmass

       0.1       0.0         1         0         0         0         1         0

$#  usrstr    usrfrc     nsbcs    interm     xpene     ssthk      ecdt   tiedprj

         0         0         0         0       4.0         0         0         0

$#   sfric     dfric       edc       vfc        th     th_sf    pen_sf      

       0.0       0.0       0.0       0.0       0.0       0.0       0.0

$#  ignore    frceng   skiprwg    outseg   spotstp   spotdel   spothin       

         0         0         0         0         0         0       0.0

$#    isym    nserod    rwgaps    rwgdth     rwksf      icov    swradf    ithoff

         1         0         1       0.0       1.0         0       0.0         0

$#  shledg    pstiff    ithcnt    tdcnof     ftall    unused    shltrw      

         1         0         0         0         0                 0.0

*CONTROL_TIMESTEP

$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st

       0.0       0.6         0       0.0       0.0         0         0         0

$#  dt2msf   dt2mslc     imscl    unused    unused     rmscl     

       0.0         0         0                           0.0