关于水下爆炸不爆炸的原因？求解答!

chenzecong

有关水下爆炸炸药没爆炸的原因？觉得问题应该是和关键字*INITIAL_VOLUME_FRACTION_GEOMETRY有关？一般可能是什么原因，另外出现炸药的表面节点上out-of-range velocities是什么原因，应如何解决？

chenzecong

附件是删掉节点和单元信息后的不完整k文件。
*CONTROL_ALE
$#     dct      nadv      meth      afac      bfac      cfac      dfac      efac
         2         1         2 -1.000000     0.000     0.000     0.000     0.000
$#   start       end     aafac     vfact      prit       ebc      pref   nsidebc
     0.000     0.000     0.000     0.000         0         0     0.000         0
$#    ncpl      nbkt    imascl    checkr
         1        50         0     0.000
*CONTROL_CPU
$#  cputim
     0.000
*CONTROL_TERMINATION
$#  endtim    endcyc     dtmin    endeng    endmas
2500.0000         0     0.000     0.000     0.000
*CONTROL_TIMESTEP
$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st
     0.000  0.010000         0     0.000     0.000         0         1         0
$#  dt2msf   dt2mslc     imscl
     0.000         0         0
*DATABASE_GLSTAT
$#      dt    binary      lcur     ioopt
50.000000         0         0         1
*DATABASE_MATSUM
$#      dt    binary      lcur     ioopt
50.000000         0         0         1
*DATABASE_BINARY_D3PLOT
$#      dt      lcdt      beam     npltc    psetid
50.000000         0         0         0         0
$#   ioopt
         0
*DATABASE_BINARY_D3THDT
$#      dt      lcdt      beam     npltc    psetid
50.000000         0         0         0         0
*DATABASE_FORMAT
$#   iform   ibinary
         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         0
$# nintsld   pkp_sen      sclp    unused     msscl     therm    intout    nodout
         0         0  1.000000         0         0         0STRESS    STRESS
*BOUNDARY_NON_REFLECTING
$#    ssid        ad        as
         1     0.000     0.000
*BOUNDARY_SPC_SET
$#    nsid       cid      dofx      dofy      dofz     dofrx     dofry     dofrz
         1         0         1         1         1         1         1         1
*SET_NODE_LIST_TITLE
FIX
$#     sid       da1       da2       da3       da4    solver
         1     0.000     0.000     0.000     0.000MECH
$#    nid1      nid2      nid3      nid4      nid5      nid6      nid7      nid8
       200         1       199       198       197       196       195       194
       193       192       191       190       189       188       187       186
*PART
$# title
Steel
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         1         1         1         0         1         0         0         0
*SECTION_SHELL
$#   secid    elform      shrf       nip     propt   qr/irid     icomp     setyp
         1         2  1.000000         3         1         0         0         1
$#      t1        t2        t3        t4      nloc     marea      idof    edgset
  0.600000  0.600000  0.600000  0.600000     0.000     0.000     0.000         0
*MAT_PLASTIC_KINEMATIC_TITLE
Steel
$#     mid        ro         e        pr      sigy      etan      beta
         1  7.830000  2.700000  0.300000  0.003000  0.002500     0.000
$#     src       srp        fs        vp
40.000000  5.000000  0.200000     0.000
*HOURGLASS
$#    hgid       ihq        qm       ibq        q1        q2    qb/vdc        qw
         1         0     0.000         0     0.000     0.000     0.000     0.000
*PART
$# title
Water
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         2         2         2         2         1         0         0         0
*SECTION_SOLID_ALE
$#   secid    elform       aet
         2        11         1
$#    afac      bfac      cfac      dfac     start       end     aafac
     0.000     0.000     0.000     0.000     0.000     0.000     0.000
*MAT_NULL_TITLE
Water
$#     mid        ro        pc        mu     terod     cerod        ym        pr
         2  1.000000     0.000     0.000     0.000     0.000     0.000     0.000
*EOS_GRUNEISEN_TITLE
Water
$#   eosid         c        s1        s2        s3     gamao         a        e0
         2  0.148000  1.750000     0.000     0.000  0.280000     0.000     0.000
$#      v0
  1.000000
*PART
$# title
Air
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         3         3         3         3         1         0         0         0
*SECTION_SOLID_ALE
$#   secid    elform       aet
         3        11         1
$#    afac      bfac      cfac      dfac     start       end     aafac
     0.000     0.000     0.000     0.000     0.000     0.000     0.000
*MAT_NULL_TITLE
Air
$#     mid        ro        pc        mu     terod     cerod        ym        pr
         3  0.001250     0.000     0.000     0.000     0.000     0.000     0.000
*EOS_LINEAR_POLYNOMIAL_TITLE
Air
$#   eosid        c0        c1        c2        c3        c4        c5        c6
         3     0.000     0.000     0.000     0.000     0.000     0.000     0.000
$#      e0        v0
     0.000     0.000
*PART
$# title
TNT
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         4         4         4         4         1         0         0         0
*SECTION_SOLID_ALE
$#   secid    elform       aet
         4        11         1
$#    afac      bfac      cfac      dfac     start       end     aafac
     0.000     0.000     0.000     0.000     0.000     0.000     0.000
*MAT_HIGH_EXPLOSIVE_BURN_TITLE
TNT
$#     mid        ro         d       pcj      beta         k         g      sigy
         4  1.630000  0.693000  0.270000     0.000     0.000     0.000     0.000
*EOS_JWL_TITLE
TNT
$#   eosid         a         b        r1        r2      omeg        e0        vo
         4     0.000     0.000     0.000     0.000     0.000     0.000     0.000
*INITIAL_VOLUME_FRACTION_GEOMETRY
$#   fmsid   fmidtyp     bammg    ntrace
         2         1         2         3
$# conttyp   fillopt     fammg        vx        xy        xz    radvel    unused
         5         0         1     0.000     0.000     0.000         0         0
$#    xmin      ymin      zmin      xmax      ymax      zmax    unused    unused
-0.200000 -0.200000 -0.800000  0.200000  0.200000 -0.400000     0.000     0.000
*INITIAL_DETONATION
$#     pid         x         y         z        lt
         4     0.000     0.000 -0.600000     0.000
*SET_SEGMENT_TITLE
NON_RFE
$#     sid       da1       da2       da3       da4    solver
         1     0.000     0.000     0.000     0.000MECH
$#      n1        n2        n3        n4        a1        a2        a3        a4
     22719     22720     22696     22695     0.000     0.000     0.000     0.000
     22455     22441     22440     22454     0.000     0.000     0.000     0.000

chenzecong

已解决，原因是水域和空气域的网格之间没有共节点

chenzecong

已解决，原因是水域和空气域的网格之间没有共节点

xiao0015

chenzecong 发表于 2016-6-23 13:24
已解决，原因是水域和空气域的网格之间没有共节点

水下爆炸没有空气域吧，是和炸药网格没有公共节点吧？