为什么用H-J-C模型进行混凝土侵彻仿真没有崩裂，而是像海绵似的被挤下去啦

youli

K文件如下所示，请各位帮忙看看错在哪儿？
*KEYWORD
*TITLE
                                                                                
$
*DATABASE_FORMAT
         0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                               NODE DEFINITIONS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*NODE
       1 0.000000000E+00 0.000000000E+00 0.000000000E+00       0       0
       2 7.500000000E-01 0.000000000E+00 0.000000000E+00       0       0
       3 1.250000000E-01 0.000000000E+00 0.000000000E+00       0       0
       4 2.500000000E-01 0.000000000E+00 0.000000000E+00       0       0
       5 3.750000000E-01 0.000000000E+00 0.000000000E+00       0       0
       6 5.000000000E-01 0.000000000E+00 0.000000000E+00       0       0
       7 6.250000000E-01 0.000000000E+00 0.000000000E+00       0       0
       8 0.000000000E+00 7.500000000E-01 0.000000000E+00       0       0
       9 0.000000000E+00 6.250000000E-01 0.000000000E+00       0       0
      10 0.000000000E+00 5.000000000E-01 0.000000000E+00       0       0
  
  .......................................................................
  116867 1.920000000E+01 4.400000000E+00-7.800000000E+00       0       0
  116868 1.960000000E+01 4.400000000E+00-7.800000000E+00       0       0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             SECTION DEFINITIONS                              $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*SECTION_SOLID
         1         1
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             MATERIAL DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*MAT_JOHNSON_COOK
         1  7.83      0.77     
12.61e-03 5.107e-03     0.018 0.30      0.94     0.1793e+04       294 0.100e-05               
0.477e-05 -9.00e+00      3.00       0.0      4.00      0.00      0.00      0.00         
      0.00   
*EOS_GRUNEISEN
         1    0.4569     1.490      0.00      0.00      2.17      0.46      0.00         
  1.00  
*MAT_JOHNSON_HOLMQUIST_CONCRETE
         2       2.4    0.1486      0.79      1.60     0.007      0.61   0.00048
      4E-5      1E-6      0.01       7.0   0.00016     0.001    0.0080      0.10
      0.04       1.0      0.85     -1.71      2.08
$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              PARTS DEFINITIONS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
*PART
Part          1 for Mat         1 and Elem Type         1
         1         1         1         1         0         0         0
$
*PART
Part          2 for Mat         2 and Elem Type         1
         2         1         2         0         0         0         0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             ELEMENT DEFINITIONS                              $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*ELEMENT_SOLID
       1       1      16      24      29      15     231     816    1050     192
       2       1      24      25      28      29     816     855    1128    1050
       3       1      25      19      21      28     855     777    1011    1128
       4       1      15      29      27      14     192    1050    1089     153
       5       1      29      28      26      27    1050    1128    1167    1089
       6       1      28      21      20      26    1128    1011     972    1167
       7       1      14      27       7       2     153    1089     387     114
       8       1      27      26       6       7    1089    1167     388     387
       9       1      26      20       5       6    1167     972     389     388
      
   ............................................................................
  109079       2  115737  116868   70068   70067   71658   71629   68547   68546
  109080       2  116868   41366   41308   70068   71629   41309   41269   68547
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              COORDINATE SYSTEMS                              $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                               LOAD DEFINITIONS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                                RIGID BOUNDRIES                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             BOUNDARY DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*SET_NODE_LIST
         1     0.000     0.000     0.000     0.000
         1        45       348       349       350       351       352       353
       354       355       356       357       358       359       360       361
       362       363       364       365       366       367       368       369
       370       371       372       373       374       375       376       377
       378       379       380       381       382       383       384       385
       386      2169      2189      2190      2191      2192      2193      2194
      2195      2196      2197      2198      2199      2200      2201      2202
      2203      2204      2205      2206      2207      2208      2209      2210
      2211      2212      2213      2214      2215      2216      2217      2218
*BOUNDARY_SPC_SET
         1         0         1         1         0         1         1         1
*SET_NODE_LIST
         2     0.000     0.000     0.000     0.000
         2         3         4         5         6         7        38        51
        52        53        54        55       114       115       116       117
       118       119       120       121       122       123       124       125
       126       127       128       129       130       131       132       133
       134       135       136       137       138       139       140       141
       142       143       144       145       146       147       148       149
       150       151       152       387       388       389       390       391
       392       393       394       395       396       397       398       399
    ............................................................................
     45737     45738     45739     45740     45741     45742     45743     45744
     45745     45746     45747     45748     45749     45750     45751     45752
     45753     45754     45755     45756     45757     45758     45759     45760
     45761     45762     45763     45764     45765     45766     45767     45768
     45769     45770     45771     45772     45773     45774     45775     45776
     45777     45778     45779
*BOUNDARY_SPC_SET
         2         0         0         1         0         1         0         1
*SET_NODE_LIST
         3     0.000     0.000     0.000     0.000
         8         9        10        11        12        13        39        46
        47        48        49        50        75        76        77        78
        79        80        81        82        83        84        85        86
        87        88        89        90        91        92        93        94
      ..........................................................................
     19713     19714     19715     19716     19717     19718     19719     19720
     19721     19722     19723     19724     19725     19726     19727     19728
     19729     19730     19731     19732     19733     19734     19735     19736
     19737     19738     19739
*BOUNDARY_SPC_SET
         3         0         1         0         0         0         1         1
*SET_NODE_LIST
         4     0.000     0.000     0.000     0.000
     15189     15229     15230     15231     15232     15233     15234     15235
     ..........................................................................
     72729     72730     72731     72732     72733     72734     72735     72736
     72737     72738     72739     72740     72741     72742     72743     72744
     72745     72746     72747     72748     72749     72750     72751     72752
     72753     72754     72755     72756     72757     72758     72759
*BOUNDARY_SPC_SET
         4         0         1         1         1         0         0         0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              CONTACT DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*CONTACT_ERODING_SURFACE_TO_SURFACE
         1         2         3         3         0         0         0         0
0.000     0.000     0.000     0.000     0.000             1  
1.000     1.000     0.000     0.000     1.000     1.000     1.000     1.000   
         1         1         1
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              CONTROL OPTIONS                                 $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*CONTROL_CONTACT
         2      
            
*CONTROL_ENERGY
         1         2         2         1
*CONTROL_BULK_VISCOSITY
  1.00     0.600E-01
*CONTROL_TIMESTEP
    0.0000    0.9000         0  0.00      0.00   
*CONTROL_TERMINATION
  50.0             0   0.00000   0.00000   0.00000
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                                 TIME HISTORY                                 $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_BINARY_D3PLOT
1.000   
*DATABASE_BINARY_D3THDT
1.000   
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              DATABASE OPTIONS                                $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_EXTENT_BINARY
         0         0         3         1         0         0         0         0
         0         0         4         0         0         0
*INITIAL_VELOCITY_GENERATION
         1         2                                   -0.15
                  
*END

[ 本帖最后由 youli 于 2006-11-12 22:19 编辑 ]

redboy202

我把你的k文件重新运算了一遍，弹体被穿透了。这就不对了！我觉得应该加*MAT_ADD_EROSION。但我对这个也不熟。大家一起想想

yulyu

能贴出最后侵彻完后的图吗.?
网格密度?失效的定义?

weizq

速度太大了吧
1500m/s
这个也挺正常的

redboy202

我是有看到加*MAT_ADD_EROSION的文章，不过我也不清楚这样作对不对。怎么样建新模型？

youli

*MAT_JOHNSON_COOK
         1  7.85      0.7275     2.01       0.3
9.5e-03  6.000e-03     0.260 0.140e-01      1.03 0.1793e+04       294 0.100e-05               
0.477e-05 -9.00e+00      3.00       0.0      0.80      0.00      0.00      0.00         
      0.00   
*EOS_GRUNEISEN
         1    0.4569     1.490      0.00      0.00      2.17      0.46      0.00         
  1.00   
*MAT_JOHNSON_HOLMQUIST_CONCRETE
         2       2.4    0.1486      0.79      1.60     0.007      0.61   0.00048
      4E-5      1E-6      0.01       7.0   0.00016     0.001    0.0080      0.10
      0.04       1.0      0.85     -1.71      2.08
$*MAT_ADD_EROSION                                                      
$        2     12345                                                  
$    -4e-5     12345     12345     12345     12345     12345     12345
*MAT_ADD_EROSION
         2     12345
     -4e-5     12345     12345     12345      0.18     12345     12345
*HOURGLASS
         1         3 0.100             0  1.00     0.600E-01  0.00      0.00   
*HOURGLASS
         2         3 0.100             0  1.00     0.600E-01  0.00      0.00  
材料参数设成这样后，侵彻时混凝土就会产生崩落现象。
参看970手册，就有各项参数的意义。
只是不明白*MAT_ADD_EROSION对材料2要设两次？

redboy202

这个好像我也看到过。但我个人觉得它并不能很好的模拟。那个$后面的*MAT_ADD_EROSION是用来解释的。没有实际的意义。

yulyu

只是不明白*MAT_ADD_EROSION对材料2要设两次？ ,,


加了add_erosin.的材料号.,

就是说,针对此材料号的材料如果满足所定义的任意一种失效方式,单元就删除.

redboy202

什么是单点全积分。在选择单元中是不是用Full  Int S/R

youli

为什么弹体经过的地方混凝土还在，是定义接触的问题还是材料定义的问题？

redboy202

问一下：你们有没有A computational constitutive model for concrete subjective to large
strain，high strain rates，and high pressure这篇原文。有的话能否再贴出来。万分感谢！我也下了一些外文。大家一起互换吧

elmund

正在做这个例子，希望能帮上忙

静水流深11

不是定义两次。实际上只定义了一次。有一次定义前面有$,这个符号的意思是注释，电脑是不读入这个信息的。

fuyuyl

弹体换个材料模型试试，如3号材料

tianjitingyu

youli 发表于 2006-11-19 22:47
为什么弹体经过的地方混凝土还在，是定义接触的问题还是材料定义的问题？ ...

这个是穿透了，你修改一下材料参数或者是罚函数因子增大一下

yingongzi

youli 发表于 2006-11-19 22:47
为什么弹体经过的地方混凝土还在，是定义接触的问题还是材料定义的问题？ ...

单位不统一，你检查单位制~

zhaidx101

yingongzi 发表于 2017-10-26 10:53
单位不统一，你检查单位制~

你也不瞅瞅发帖时间的？~~

mry

zhaidx101 发表于 2017-12-5 21:02
你也不瞅瞅发帖时间的？~~

哈哈哈，忍不住了

Mr.YANG

zhaidx101 发表于 2017-12-5 21:02
你也不瞅瞅发帖时间的？~~

可以给后来人一点提示，也挺好的