关于爆炸的问题

18761687231

我做了一个关于爆炸的数值模拟，用manager求解完之后，只是炸药所在的地方有爆炸的现象，炸药周围的土壤空气以及其他介质都没有反应，这是怎么回事啊，求大神指导啊！*KEYWORD
*TITLE
$# title
LS-DYNA keyword deck by LS-PrePost
*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.0001.0000E+20  1.000000 1.0000E-6         0         0     0.000         0
$#    ncpl      nbkt    imascl    checkr
         1        50         0     0.000
*CONTROL_BULK_VISCOSITY
$#      q1        q2      type     btype
  1.500000  0.060000         1         0
*CONTROL_CPU
$#  cputim
     0.000
*CONTROL_ENERGY
$#    hgen      rwen    slnten     rylen
         2         1         2         2
*CONTROL_SHELL
$#  wrpang     esort     irnxx    istupd    theory       bwc     miter      proj
20.000000         1        -1         1         2         2         1         0
$# rotascl    intgrd    lamsht    cstyp6    tshell
  1.000000         0         0         1         0
$# psstupd   sidt4tu     cntco    itsflg    irquad
         0         0         0         0         2
$#  nfail1    nfail4   psnfail    keepcs     delfr   drcpsid    drcprm
         0         0         0         0         0         0  1.000000
*CONTROL_TERMINATION
$#  endtim    endcyc     dtmin    endeng    endmas
1500.0000         0     0.000     0.000     0.000
*CONTROL_TIMESTEP
$#  dtinit    tssfac      isdo    tslimt     dt2ms      lctm     erode     ms1st
     0.000  0.670000         0     0.000     0.000         0         0         0
$#  dt2msf   dt2mslc     imscl    unused    unused     rmscl
     0.000         0         0                         0.000
*DATABASE_BINARY_D3PLOT
$#      dt      lcdt      beam     npltc    psetid
15.000000         0         0         0         0
$#   ioopt
         1
*DATABASE_BINARY_D3THDT
$#      dt      lcdt      beam     npltc    psetid
  1.500000         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         0         0         1         1         1         2         1
$# nintsld   pkp_sen      sclp    unused     msscl     therm    intout    nodout
         0         0  1.000000                   0         0STRESS    STRESS
$#    dtdt    resplt
         0         0
*BOUNDARY_NON_REFLECTING
$#    ssid        ad        as
         1     0.000     0.000
         2     0.000     0.000
*BOUNDARY_SPC_SET
$#    nsid       cid      dofx      dofy      dofz     dofrx     dofry     dofrz
         1         0         0         1         0         0         0         0
*SET_NODE_LIST_TITLE
NODESET(SPC) 1
$#     sid       da1       da2       da3       da4    solver
         1     0.000     0.000     0.000     0.000MECH
$#    nid1      nid2      nid3      nid4      nid5      nid6      nid7      nid8
         1         2         3         4         5         6         7        22
        23        24        25        26        27        28        43        44
        45        46        47        48        49        64        65        66
        67        68        69        70        85        86        87        88

.........
     47047     47048     47049     47050     47261     47262     47263     47264
     47265     47266     47267     47478     47479     47480     47481     47482
     47483     47484     47695     47696     47697     47698     47699     47700
     47701     47912     47913     47914     47915     47916     47917     47918
*BOUNDARY_SPC_SET
$#    nsid       cid      dofx      dofy      dofz     dofrx     dofry     dofrz
         2         0         1         0         0         0         0         0
*SET_NODE_LIST_TITLE
NODESET(SPC) 2
$#     sid       da1       da2       da3       da4    solver
         2     0.000     0.000     0.000     0.000MECH
$#    nid1      nid2      nid3      nid4      nid5      nid6      nid7      nid8
         1         2         3         4         5         6         7         8
         9        10        11        12        13        14        15        16
        17        18        19        20        21        64        65        66
        67        68        69        70        71        72        73        74
        75        76        77        78        79        80        81        82
        83        84       127       128       129       130       131       132
       133       134       135       136       137       138       13

........
    41587     41588     41589     41590     41591     41592     41593     41594
     41595     41596     41597     41598     41599     41600     41601     41602
     41603     41604     41605     41606     41607     41608     41609     41610
     41611     41612     41613     41614     41615     41616     41617     41618
*PART
$# title
air
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         1         1         1         1         0         0         0         0
*SECTION_SOLID_TITLE
air
$#   secid    elform       aet
         1        11         0
*MAT_NULL
$#     mid        ro        pc        mu     terod     cerod        ym        pr
         1  0.001250     0.000     0.000     0.000     0.000     0.000     0.000
*EOS_GRUNEISEN
$#   eosid         c        s1        s2        s3     gamao         a        e0
         1  0.344000     0.000     0.000     0.000  1.400000     0.000     0.000
$#      v0
     0.000
*PART
$# title
burn
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         2         2         2         2         2         0         0         0
*SECTION_SOLID_TITLE
burn
$#   secid    elform       aet
         2        11         0
*MAT_HIGH_EXPLOSIVE_BURN
$#     mid        ro         d       pcj      beta         k         g      sigy
         2  1.630000  0.693000  0.272000     0.000     0.000     0.000     0.000
*EOS_JWL
$#   eosid         a         b        r1        r2      omeg        e0        vo
         2  3.710000  0.074300  4.150000  0.950000  0.300000  0.070000  1.000000
*HOURGLASS
$#    hgid       ihq        qm       ibq        q1        q2    qb/vdc        qw
         2         0  0.100000         0  1.500000  0.060000  0.100000  0.100000
*PART
$# title
soil
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         3         3         3         0         0         0         0         0
*SECTION_SOLID_TITLE
soil
$#   secid    elform       aet
         3        11         0
*MAT_SOIL_AND_FOAM
$#     mid        ro         g      bulk        a0        a1        a2        pc
         3  1.800000 1.6015E-4 132.80000  0.003300 1.3100E-7  0.123200     0.000
$#     vcr       ref
     0.000     0.000
$#    eps1      eps2      eps3      eps4      eps5      eps6      eps7      eps8
     0.000  0.050000  0.090000  0.110000  0.150000  0.190000  0.210000  0.220000
$#    eps9     eps10
  0.250000  0.300000
$#      p1        p2        p3        p4        p5        p6        p7        p8
     0.000  0.034200  0.045300  0.067600  0.127000  0.208000  0.271000  0.392000
$#      p9       p10
  0.566000  1.230000
*PART
$# title
JHC
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         4         4         4         0         4         0         0         0
*SECTION_SOLID_TITLE
JHC
$#   secid    elform       aet
         4        11         0
*MAT_JOHNSON_HOLMQUIST_CONCRETE_TITLE
JHC
$#     mid        ro         g         a         b         c         n        fc
         4  2.400000  0.148600  0.790000  1.600000  0.007000  0.610000 4.8000E-4
$#       t      eps0     efmin     sfmax        pc        uc        pl        ul
4.0000E-5 1.0000E-6  0.010000  7.000000 1.6000E-4  0.001000  0.008000  0.100000
$#      d1        d2        k1        k2        k3        fs
  0.040000  1.000000  0.850000 -1.710000  2.080000  0.010000
*HOURGLASS
$#    hgid       ihq        qm       ibq        q1        q2    qb/vdc        qw
         4         3  0.100000         0  1.500000  0.060000     0.000     0.000
*PART
$# title
foam
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         5         5         5         0         0         0         0         0
*SECTION_SOLID_TITLE
foam
$#   secid    elform       aet
         5         1         0
*MAT_FU_CHANG_FOAM_TITLE
foam
$#     mid        ro         e        ed        tc      fail      damp      tbid
         5  0.055000 1.1500E-4     0.000 -4.482E-6     0.000  0.100000         1
$#  bvflag     sflag     rflag     tflag      pvid      sraf       ref        hu
  1.000000  1.000000  1.000000  1.000000         0     0.000     0.000  0.050000
$#      d0        n0        n1        n2        n3        c0        c1        c2
     0.000     0.000     0.000     0.000     0.000     0.000     0.000     0.000
$#      c3        c4        c5       aij       sij      minr      maxr     shape
     0.000     0.000     0.000     0.000     0.000     0.000     0.000     0.000
*PART
$# title
concrete
$#     pid     secid       mid     eosid      hgid      grav    adpopt      tmid
         6         6         6         0         0         0         0         0
*SECTION_SOLID_TITLE
concrete
$#   secid    elform       aet
         6         1         0
*MAT_SOIL_CONCRETE
$#     mid        ro         g         k      lcpv      lcyp      lcfp      lcrp
         6  2.500000  0.125000  0.166700         2         3         4         5
$#      pc       out         b      fail
-2.000E-5     0.000  0.300000  1.000000
*INITIAL_DETONATION
$#     pid         x         y         z        lt
         2     0.000     0.000-120.00000     0.000
*DEFINE_CURVE
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp
         1         0  1.000000  1.000000     0.000     0.000         0
$#                a1                  o1
               0.000               0.000
       5.000000e-008            0.010000
       1.000000e-007            0.020000
       2.500000e-007            0.028000
       5.000000e-007            0.032000
       7.500000e-007            0.037000
       1.000000e-006            0.040000
       1.250000e-006            0.045000
       1.500000e-006            0.047000
       1.750000e-006            0.050000
       2.000000e-006            0.054000
       2.250000e-006            0.060000
       2.500000e-006            0.100000
       2.500000e-006            0.180000
*DEFINE_CURVE
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp
         2         0     0.000     0.000     0.000     0.000         0
$#                a1                  o1
               0.000               0.000
            0.095300            0.023000
            0.182300            0.065000
            0.262400            0.130000
            0.336500            0.175000
            0.405500            0.225000
            0.470000            0.350000
*DEFINE_CURVE
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp
         3         0     0.000     0.000     0.000     0.000         0
$#                a1                  o1
               0.000               0.000
            0.004000            0.002400
            0.008000            0.005500
            0.012000            0.009500
            0.014000            0.012600
            0.018000            0.015500
            0.100000            0.015500
*DEFINE_CURVE
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp
         4         0     0.000     0.000     0.000     0.000         0
$#                a1                  o1
               0.000               0.000
       1.500000e-004               0.000
            0.001000            0.004000
            0.010000            0.490000
            0.100000            4.900000
*DEFINE_CURVE
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp
         5         0     0.000     0.000     0.000     0.000         0
$#                a1                  o1
               0.000            0.010000
       1.500000e-004            0.010000
            0.001000            0.050000
            0.010000            0.500000
            0.100000            5.000000
*SET_SEGMENT
$#     sid       da1       da2       da3       da4    solver
         1     0.000     0.000     0.000     0.000MECH
$#      n1        n2        n3        n4        a1        a2        a3        a4
     47941     47948     47947     47940     0.000     0.000     0.000     0.000
     41311     41317     41316     41310     0.000     0.000     0.000     0.000
     35446     35470     35469     35445     0.000     0.000     0.000     0.000
     35191     35215     35214     35190     0.000     0.000     0.000     0.000

.......
    35190     35214     35213     35189     0.000     0.000     0.000     0.000
*SET_SEGMENT
$#     sid       da1       da2       da3       da4    solver
         2     0.000     0.000     0.000     0.000MECH
$#      n1        n2        n3        n4        a1        a2        a3        a4
     43351     43350     43567     43568     0.000     0.000     0.000     0.000
     29581     29580     30252     30.......

    43350     43349     43566     43567     0.000     0.000     0.000     0.000
     29580     29579     30251     30252     0.000     0.000     0.000     0.000
     22185     22184     22856     22857     0.000     0.000     0.000     0.000
*ALE_MULTI-MATERIAL_GROUP
$#     sid    idtype    gpname
         1         1         
         2         1         
         3         1         
         4         1         
$ ELEMENT CARDS FOR SOLID ELEMENTS
*ELEMENT_SOLID
$#   eid     pid      n1      n2      n3      n4      n5      n6      n7      n8
       1       1       1      22      29       8       2      23      30       9
       2       1      22      43      50      29      23      44      51      30
       3       1       8      29 。。。。。。

41404       6   47904   48121   48128   47911   47903   48120   48127   47910
$
$ NODES
$
*NODE
$#   nid               x               y               z      tc      rc
       1           0.000           0.000      -40.000000       0       0
       2           0.000           0.000           0.000       0       0
       3           0.000     .........

  48126       1200.0000       1200.0000      -1560.0000       0       0
   48127       1200.0000       1200.0000      -1600.0000       0       0
   48128       1200.0000       1200.0000      -1640.0000       0       0
*END

babymice

真乱！没看明白，但是根据个人经验，总结如下：
（1）设置炸药的爆炸坐标，不要忘记设置了，即需要起爆控制
（2）如果设置了起爆控制的话，炸药一定会传播的，不管其他材料是啥，你再看看。

920706938@qq.co

我也遇到相同的问题，请问你解决了吗

18761687231

920706938@qq.co 发表于 2015-10-25 09:50
我也遇到相同的问题，请问你解决了吗

解决了，我最后用的流固耦合， 炸药和空气共节点，空气把其他的材料都包住

fys

楼主，你也在做炸药爆炸吗？？能不能交流一下，我做炸药接触钢板爆炸，刚开始，能不能请教一下您？？

18761687231

fys 发表于 2016-1-12 10:01
楼主，你也在做炸药爆炸吗？？能不能交流一下，我做炸药接触钢板爆炸，刚开始，能不能请教一下您？？ ...

你有什么问题直接说吧  我也有一段时间没有做了  不知道能不能给你解决了

fys

对您来说应该是比较简单的吧，我现在按照石少卿老师的那本书做实例，但是发现他用的老版本，现在修改材料属性需要用Change Material Number 现在么有了

fys

18761687231 发表于 2016-1-12 10:05
你有什么问题直接说吧  我也有一段时间没有做了  不知道能不能给你解决了 ...

哥，能不能加一下您的qq什么的。。。感谢

18761687231

fys 发表于 2016-1-12 10:53
哥，能不能加一下您的qq什么的。。。感谢

我现在用的是lsdyna  没用ANSYS 我的qq是 522195945

fys

18761687231 发表于 2016-1-12 16:16
我现在用的是lsdyna  没用ANSYS 我的qq是 522195945

哥，就是用lsdyna

zhoust0707

fys 发表于 2016-1-12 10:43
对您来说应该是比较简单的吧，我现在按照石少卿老师的那本书做实例，但是发现他用的老版本，现在修改材料属 ...

我用的也是14.5的，你可以输入命令流mpchg,2,all更改材料属性。2为所选单元

zhoust0707

18761687231 发表于 2016-1-11 19:38
解决了，我最后用的流固耦合， 炸药和空气共节点，空气把其他的材料都包住 ...

哈哈，刚想回答你，可能是结点的问题，你已经解决了。我问一下，DEFINE CURVE怎么用，因为我现在要在模型中加重力，需要定义重力加速度

18761687231

zhoust0707 发表于 2016-1-13 11:00
哈哈，刚想回答你，可能是结点的问题，你已经解决了。我问一下，DEFINE CURVE怎么用，因为我现在要在模型 ...

定义荷载的话  就是直接定义你想要加的荷载曲线 就可以了     如果需要定义重力加速度的话好像是有定力加速度的关键字 。define curve 不能定义加速度的，你可以看一下英文关键字手册

鱼丫丫

18761687231 发表于 2016-1-11 19:38
解决了，我最后用的流固耦合， 炸药和空气共节点，空气把其他的材料都包住 ...

请问楼主做的是三维的吗  我三维的流固耦合后 岩体依然没有反应 请问楼主如何解决，谢谢

huangyh

成功的k文件能否上传学习，谢谢