爆破漏斗SPH模拟计算就停在某一时刻，计算也不中止，也不提示错误是怎么回事？ 40

浏览：1356 回答：5

我做的是炸药在土壤中的爆炸，炸药和近区岩石采用的是SPH粒子，但是计算没有多久，计算就停在了83us，也不中止，也不提示错误。我把k文件中的除了单元节点外的信息写了，麻烦各位老师帮我看看是怎么回事。

*KEYWORD 200000000
*TITLE
$#                                                                         title

*CONTROL_ENERGY
         2         2         2         2
*CONTROL_SPH
         0         0 1.000E+20         3
*CONTROL_TERMINATION
400.0000
*CONTROL_TIMESTEP
     0.000 0.900000
     0.000         0
*DATABASE_BINARY_D3PLOT
2.000000
*DATABASE_FORMAT
         0         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.000000         0         0         0STRESS    STRESS
$#    dtdt    resplt
         0         0
$
$
$外层岩石
*SECTION_SOLID
         2         1
*PART
$#                                                                         title
Part          1 for Mat         2 and Elem Type         2
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         1         2         2         0         0         0         0         0
*MAT_JOHNSON_HOLMQUIST_CONCRETE
         2       2.3       0.1      0.55      1.23     0.097      0.896.0000E-04
4.0000E-051.0000E-06      0.01        202.0000E-041.2500E-03      0.02     0.174
      0.04       1.0      0.39     -2.23       5.5       0.0
$
$内层岩石SPH粒子
*SECTION_SPH
         1 1.2000000
*PART
$#                                                                         title
SphNode
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         2         1         1         0         0         0         0         0
*MAT_JOHNSON_HOLMQUIST_CONCRETE
         1       2.3       0.1      0.55      1.23     0.097      0.896.0000E-04
4.0000E-051.0000E-06      0.01        202.0000E-041.2500E-03      0.02     0.174
      0.04       1.0      0.39     -2.23       5.5       0.0
$
$
$炸药SPH粒子
*SECTION_SPH
         3 1.2000000
*PART
$#                                                                         title
SphNode
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         3         3         3         3         0         0         0         0
*MAT_HIGH_EXPLOSIVE_BURN
         3       1.3       0.4     0.099       0.0
*EOS_JWL
         3     2.144   0.00182       4.2       0.9      0.15    0.0419
*INITIAL_DETONATION
         3      0.00     40.00      0.00      0.00
$
$
$
$固连接触
*CONTACT_TIED_NODES_TO_SURFACE_OFFSET
         1         1         4         0
     0.000     0.000     0.000     0.000     0.000         1
1.000000 1.000000     0.000     0.000 1.000000 1.000000 1.000000 1.000000
$
$对称设置
*CONSTRAINED_GLOBAL
         2         6         2
         3         4         3
*BOUNDARY_SPH_SYMMETRY_PLANE
     0.000     0.000     0.000 1.000000
*BOUNDARY_SPH_SYMMETRY_PLANE
     0.000     0.000     0.000     0.000     0.000 1.000000