为什么靶板突然消失了？？？？大神指点

yanpenghui

T=16250时间 还好好的
T=15500时间 靶板就突然消失不见了 不知道什么原因。
C30靶板 HJC模型。


大神帮忙分析一下 是什么原因？

lzclsrlxy

yanpenghui 发表于 2018-5-2 09:31
顶啊 我自己顶啊

我帮你测试了一下你的K文件，总共跑了3遍。（没想到你单元这么多，一遍就是4小时）
第一遍：
1.通过查看LS-DYNA说明文档（其实我利用五一两天半时间才把add_erosion翻译完），发现你的add_erosion参数用得有问题。其中，numfip、ncs是满足失效准则的数目，不需要设置成1234，如果是一种准则失效，设置成1就可以了，这样不会有warning。
2.但是到达一定时间后仍会有整块part消失的情况。我猜测是你设置的mnpres(  -4.5E-5)静水拉应力失效准则导致了达到一定压力后全部失效。
（其实这里还有一个关键字mxpres，这是971版本才有的，没有用过，不知道取正值用起来与mnpres能不能达到相同效果，期待大家来讨论）
第二遍：
1.通过查看HJC本构模型，发现其中带有失效关键字FS，但是说明文档非常简单，没有进行详细说明。为了测试单一失效准则对HJC材料的影响，我又将FS先设置为0，Add_erosion保持不变。重新计算，发现材料脆性更明显，而且在更短的时间内就会出现整块part删除。由此证明FS失效并不是影响材料part整体删除的关键。

第三遍：
将FS恢复原来数值（我后来取了0.67），将mnpres最小拉力失效设置为1234，即不添加任何add_eriosion失效准则，仅用材料自带关键字FS进行失效。结果果然没有在一定时间后失效，直至Termination/endtm结束为止。
即mnpres失效准则用得有问题，才会造成整体大面积失效。

然而，这样的设置虽然仿真有了一个完整的结果，但是对于混凝土侵彻的HJC本构模型而言，这样的结果并不能用以模拟真实侵彻下混凝土材料的裂纹生成、断裂及碎石崩裂的效果。
为了解决这一问题，我又用了一天时间搜了百度和论坛各个关于HJC本构模型的帖子，真是各种说法都有，也提供了一些混凝土的参数（主要都是970及之前版本的），但最终让我发现孙其然博士的这篇《 HJC 模型模拟钢筋混凝土侵彻实验的参数研究 》文章写得最为详细和透彻，孙博士借鉴了姜华、王君杰等人《弹体侵彻混凝土数值模拟失效指标 》中的静水拉应力与最大主应变双重失效标准，仿真效果非常不错。不仅如此，孙其然还将钢筋结构纳入混凝土单元中，进行了进一步的侵彻仿真分析。如果你需要进一步的完善你的仿真内容，建议去看看孙其然的论文。另外，你的HJC材料参数也需要再修改一下，看看是不是真的从实验中获取的。否则即使仿真出来，结果也会和试验数据对不上的。


好了，就帮你到这里吧，相信你沿着这个思路做下去，一定能很快解决这个问题的。
（这几天真是苦了我了）

yanpenghui

HJC材料和 材料侵蚀

*MAT_JOHNSON_HOLMQUIST_CONCRETE
         1       2.4     0.134      0.79      1.60     0.007      0.61   3.00E-4
   3.39e-5    1.0e-6      0.01       7.0   1.08e-4   7.18e-4   1.05e-2       0.1
      0.04       1.0     0.174     0.388     0.298      0.45     
*MAT_ADD_EROSION
$#     mid      excl    mxpres     mneps    effeps    voleps    numfip       ncs
         1 1234.0000 1234.0000 1234.0000 1234.0000 1234.0000 1234.0000 1234.0000
$#  mnpres     sigp1     sigvm     mxeps     epssh     sigth   impulse    failtm
   -4.5E-5 1234.0000 1234.0000      0.45 1234.0000 1234.0000 1234.0000 1234.0000

yanpenghui

厘米 克 微秒

单位

lzclsrlxy

能发一下你的K文件吗？我帮你看看，erosion的失效有个跟时间有关的你调节一下，但也有可能在时间步上，需要调试才知道。

yanpenghui

lzclsrlxy 发表于 2018-4-30 17:22
能发一下你的K文件吗？我帮你看看，erosion的失效有个跟时间有关的你调节一下，但也有可能在时间步上，需要 ...

好的，稍等一下啊

yanpenghui

lzclsrlxy 发表于 2018-4-30 17:22
能发一下你的K文件吗？我帮你看看，erosion的失效有个跟时间有关的你调节一下，但也有可能在时间步上，需要 ...

我的qq：742026961
方便加下么

yanpenghui

lzclsrlxy 发表于 2018-4-30 17:22
能发一下你的K文件吗？我帮你看看，erosion的失效有个跟时间有关的你调节一下，但也有可能在时间步上，需要 ...

*KEYWORD
*TITLE
$# title
LS-DYNA keyword deck by LS-PrePost
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$   Units: gram, cm, microsec, 1e+07 N, 10e11pa,Mbar, 1e+07 N-cm，10e4m/s
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$指定计算结果输出格式为DYNA格式
$数据库格式   二进制文件计算结果输出格式
*DATABASE_FORMAT
$#   iform   ibinary
         0         0
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$     Define Section
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$定义有限元模型单元类型1为体单元，并选择单元算法为常应变单点积分体元
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$     Define Control Ouput
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*CONTROL_TERMINATION
$#  endtim    endcyc     dtmin    endeng    endmas
19800.000         0     0.000     0.000     0.000
$$$$$$ DYNA determines initinal step size
$$$$$$ tssfac Scale factor for computed time step (default=0.9;if high explosive are used,the default is lowered to 0.67)
*CONTROL_TIMESTEP
$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st
     0.000  0.800000         0     0.000     0.000         0         0         0
$#  dt2msf   dt2mslc     imscl    unused    unused     rmscl
     0.000         0         0                         0.000
$$$$$$hourglass energy ,rigidwall energy,sliding interface energy dissipation,rayleigh dissipation
*CONTROL_ENERGY
$#    hgen      rwen    slnten     rylen
         2         2         2         2
*CONTROL_OUTPUT
$#   npopt    neecho    nrefup    iaccop     opifs    ipnint    ikedit    iflush
         1         3         0         0     0.000         0       100      5000
*CONTROL_HOURGLASS
$#     ihq        qh
         1  0.100000
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$     Time History
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*DATABASE_BINARY_D3PLOT
$#      dt      lcdt      beam     npltc    psetid
       250         0         0         0         0
$#   ioopt
         0
*DATABASE_BINARY_D3THDT
$#      dt      lcdt      beam     npltc    psetid
        10         0         0         0         0
*DATABASE_EXTENT_BINARY
$#   neiph     neips    maxint    strflg    sigflg    epsflg    rltflg    engflg
         0         0         3         0         1         1         1         1
$#  cmpflg    ieverp    beamip     dcomp      shge     stssz    n3thdt   ialemat
         0         1         0         1         1         1         2         1
*DATABASE_BINARY_D3DUMP
$#      dt      lcdt      beam     npltc    psetid
2000.0000         0         0         0         0
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$     Define  Materials
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$混凝土（23MPa混凝土）
$取自龚自明 防护工程 BLU-109B侵彻厚混凝土靶体的计算与分析
*MAT_JOHNSON_HOLMQUIST_CONCRETE
         1       2.4     0.134      0.79      1.60     0.007      0.61   3.00E-4
   3.39e-5    1.0e-6      0.01       7.0   1.08e-4   7.18e-4   1.05e-2       0.1
      0.04       1.0     0.174     0.388     0.298      0.55     
*MAT_ADD_EROSION
$#     mid      excl    mxpres     mneps    effeps    voleps    numfip       ncs
         1 1234.0000 1234.0000 1234.0000 1234.0000      0.40 1234.0000 1234.0000
$#  mnpres     sigp1     sigvm     mxeps     epssh     sigth   impulse    failtm
   -4.0E-5 1234.0000 1234.0000 1234.0000 1234.0000 1234.0000 1234.0000 1234.0000
$外壳
*MAT_PLASTIC_KINEMATIC
$#     mid        ro         e        pr      sigy      etan      beta
         2  7.850000  2.060000  0.330000  0.015000  0.020600  1.000000
$#     src       srp        fs        vp
2590.0000  5.558830  3.000000     0.000
$配重
*MAT_PLASTIC_KINEMATIC
$#     mid        ro         e        pr      sigy      etan      beta
         3  2.800000  2.060000  0.330000  0.015000  0.020600  1.000000
$#     src       srp        fs        vp
2590.0000  5.558830  3.000000     0.000
$加速度
*MAT_PLASTIC_KINEMATIC
$#     mid        ro         e        pr      sigy      etan      beta
         4  6.360000  2.060000  0.330000  0.015000  0.020600  1.000000
$#     src       srp        fs        vp
2590.0000  5.558830  3.000000     0.000
$后底盖
*MAT_PLASTIC_KINEMATIC
$#     mid        ro         e        pr      sigy      etan      beta
         5  7.850000  2.060000  0.330000  0.015000  0.020600  1.000000
$#     src       srp        fs        vp
2590.0000  5.558830  3.000000     0.000  
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$     Define Parts
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$     
*PART
concrete
$      pid       sid       mid     eosid      hgid      grav    adpopt      tmid  
         1         1         1
*PART
waike
$      pid       sid       mid     eosid      hgid      grav    adpopt      tmid  
         2         2         2                                       
*PART
peizong
$      pid       sid       mid     eosid      hgid      grav    adpopt      tmid  
         3         3         3         
*PART
jiasudu
$      pid       sid       mid     eosid      hgid      grav    adpopt      tmid  
         4         4         4        
*PART
houdigai
$      pid       sid       mid     eosid      hgid      grav    adpopt      tmid  
         5         5         5                                   
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$     Initial Condition
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$         
*INITIAL_VELOCITY_GENERATION
$#nsid/pid      styp     omega        vx        vy        vz     ivatn      icid
         1         1     0.000   -250e-4         0         0         0         0
$#      xc        yc        zc        nx        ny        nz     phase    irigid
     0.000     0.000     0.000     0.000     0.000     0.000         0         0
*SET_PART_LIST
$#     sid       da1       da2       da3       da4    solver
         1     0.000     0.000     0.000     0.000MECH
$#    pid1      pid2      pid3      pid4      pid5      pid6      pid7      pid8
         2         3         4         5      
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$     Define Contacts
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE
$#     cid                                                                 title
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         2         3         3         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
       0.5       0.5     0.000     0.000     0.000         0     0.0001.0000E+20
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
  1.000000  1.000000     0.000     0.000  1.000000  1.000000  1.000000  1.000000
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE
$#     cid                                                                 title
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         3         4         3         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
       0.5       0.5     0.000     0.000     0.000         0     0.0001.0000E+20
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
  1.000000  1.000000     0.000     0.000  1.000000  1.000000  1.000000  1.000000
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE
$#     cid                                                                 title
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         4         5         3         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
       0.5       0.5     0.000     0.000     0.000         0     0.0001.0000E+20
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
  1.000000  1.000000     0.000     0.000  1.000000  1.000000  1.000000  1.000000
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE
$#     cid                                                                 title
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         2         4         3         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
       0.5       0.5     0.000     0.000     0.000         0     0.0001.0000E+20
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
  1.000000  1.000000     0.000     0.000  1.000000  1.000000  1.000000  1.000000
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE
$#     cid                                                                 title
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         2         5         3         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
       0.5       0.5     0.000     0.000     0.000         0     0.0001.0000E+20
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
  1.000000  1.000000     0.000     0.000  1.000000  1.000000  1.000000  1.000000  
*CONTACT_ERODING_SURFACE_TO_SURFACE
$#     cid                                                                 title
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         2         1         3         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
       0.5       0.5     0.000     0.000     0.000         0     0.0001.0000E+20
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
       1.2       1.2     0.000     0.000  1.000000  1.000000  1.000000  1.000000
$#    isym    erosop      iadj
         1         1         0
*CONTROL_CONTACT
   4.80000   0.00000         2         0         1         1         1
         0         0        10         0  4.00            
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$     Define HOURGLASS
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$         
*HOURGLASS
$#    hgid       ihq        qm       ibq        q1        q2    qb/vdc        qw
         1         0  0.100000         0  1.500000 6.0000E-2  0.100000  0.100000
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$  included  keyword file   
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*INCLUDE
waike.k
*INCLUDE_TRANSFORM
peizong.k
50000,50000



*INCLUDE_TRANSFORM
jiasudu.k
80000,80000



*INCLUDE_TRANSFORM
houdigai.k
100000,100000



*INCLUDE_TRANSFORM
target1.k
550000,550000



*END

lzclsrlxy

yanpenghui 发表于 2018-4-30 20:49
*KEYWORD
*TITLE
$# title

传我邮箱吧,lzclsrlxy@163.com，我回头帮你试试

yanpenghui

lzclsrlxy 发表于 2018-4-30 21:35
传我邮箱吧,，我回头帮你试试

以及发送  谢谢大神了

yanpenghui

顶啊 我自己顶啊

yanpenghui

lzclsrlxy 发表于 2018-5-4 14:01
我帮你测试了一下你的K文件，总共跑了3遍。（没想到你单元这么多，一遍就是4小时）
第一遍：
1.通过查看L ...

跪谢大神

tianmeng

yanpenghui 发表于 2018-5-4 15:08
跪谢大神

你好，你的k文件怎么跳的，麻烦说一下吗？学习中，跪谢

tianmeng

yanpenghui 发表于 2018-5-2 09:31
顶啊 我自己顶啊

你好，你的k文件怎么跳的，麻烦说一下吗？学习中，跪谢

tianmeng

lzclsrlxy 发表于 2018-4-30 17:22
能发一下你的K文件吗？我帮你看看，erosion的失效有个跟时间有关的你调节一下，但也有可能在时间步上，需要 ...

你好，我在做爆破模拟岩石失效产生裂隙的时候也出现了岩石突然消失的情况，能不能麻烦大神给指点一下