【案例分享】LS-DYNA圆柱体内部爆炸试验仿真模拟

美国劳伦斯利弗莫尔实验室 (LLNL) 测试了一个圆柱体试验,目的是通过测量圆柱体壁在被内部炸药的爆炸产物膨胀时的运动来表征爆炸产物的绝热膨胀。

模拟结果:

以下为.K文件

$ This model is created by Eric Piskula -- eric.piskula@ansys.com
$ received on 04/12/2024
$
$# LS-DYNA Keyword file created by LS-PrePost(R) V4.10.9-08Dec2023
$# Created on Apr-11-2024 (10:49:37)
*KEYWORD
*TITLE
$#                                                                         title
LS-DYNA keyword deck by LS-PrePost
*CONTROL_ACCURACY
$      osu       inn    pidosu                                            unused
         0         1         0         0       0.0
*CONTROL_ALE
        -1         1         2      -1.0       0.0       0.0       0.0       0.0
       0.01.00000E20       1.01.00000E-6         0         01.00000E-4         0
         1        50         0       0.0       0.0         0       0.0       0.0
         0         0         01.00000E-5
*CONTROL_ENERGY
$     hgen      rwen    slnten     rylen                                  unused
         2         2         2         2         2         1         1         1
*CONTROL_MPP_DECOMPOSITION_DISTRIBUTE_ALE_ELEMENTS
*CONTROL_MPP_IO_NODUMP
*CONTROL_OUTPUT
$    npopt    neecho    nrefup    iaccop     opifs    ipnint    ikedit    iflush
         1         0         0         0       0.0         0         0         0
$    iprtf    ierode     tet10    msgmax    ipcurv                        unused
         0         0         2         0         0       0.0         0         0
*CONTROL_PARALLEL
$     ncpu    numrhs    consti      para                                  unused
         0         0         1         0
*CONTROL_TERMINATION
$   endtim    endcyc     dtmin    endeng    endmas                        unused
     0.065  10000000     0.001       0.0  100000.0         0
*CONTROL_TIMESTEP
$   endtim    endcyc     dtmin    endeng    endmas                        unused
       0.0       0.5         0       0.0       0.0         0         1         0
       0.0         0         0                           0.0       0.0         0
*DATABASE_GLSTAT
2.00000E-5         3         0         1
*DATABASE_MATSUM
2.00000E-5         3         0         1
*DATABASE_NODOUT
1.00000E-5         0         0         1       0.0         0
*DATABASE_TRHIST
2.00000E-5         0         0         1
*DATABASE_BINARY_D3PLOT
$       dt      lcdt      beam     npltc    psetid                       unused1
1.00000E-4         0         0         0         0
         0       0.0       0.0       0.0         0         0
*DATABASE_BINARY_D3PROP
$    ifile     imatl     iwall                                           unused2
         2         1         0
*DATABASE_EXTENT_BINARY
$    neiph     neips    maxint    strflg    sigflg    epsflg    rltflg    engflg
         0         0         0        10         1         1         1         1
$   cmpflg    ieverp    beamip     dcomp      shge     stssz    n3thdt   ialemat
         0         0         0         1         1         1         2         0
$  nintsld   pkp_sen      sclp     hydro     msscl     therm    intout    nodout
         0         0       0.0         0         2         0                    
*DATABASE_FORMAT
$    iform   ibinary                                                      unused
         0         1
*DATABASE_HISTORY_NODE
     75945         0         0         0         0         0         0         0
*DATABASE_TRACER
       0.0         0      15.3     244.8       0.0         2         0       0.0
*DEFINE_COORDINATE_SYSTEM
$       ID        xO        yO        zO        xL        yL        zL   unused1
        32       0.0       0.0       0.0       1.0       0.0       0.0         0
$       xP        yP        zP                                           unused2
       0.0       1.0       0.0
*DEFINE_CURVE
         1         0       1.0       1.0       0.0       0.0         0         0
                 0.0                 0.0
               0.026                20.0
*MAT_HIGH_EXPLOSIVE_BURN
$       ID        ro         d       pcj      beta         k         g      sigy
         11.63000E-6    6930.0      21.0       0.0       0.0       0.0       0.0
*MAT_RIGID
         28.94200E-6     126.0     0.345       0.0       0.0       0.0          
       1.0         7         7
       0.0       0.0       0.0       0.0       0.0       0.0
*MAT_JOHNSON_COOK
         38.96000E-6    45.926       0.0      0.35       0.0       0.0       0.0
      0.09     0.292      0.31     0.025      1.09    1656.0     293.0       1.0
     383.0       0.0       2.0       0.0       0.0       0.0       0.0       0.0
       0.0       0.0       0.01.00000E-6       0.0       0.0       0.0
*MAT_NULL
         41.00000E-9       0.0       0.0       0.0       0.0       0.0       0.0
*MAT_VACUUM
         51.00000E-9
*ELEMENT_SHELL*PART
$                                                             partname   unused1
SYS-1\Surface
         1         1         2         0         0         0         0         0
*SECTION_SHELL
$       ID    elform      shrf       nip     propt        qr     icomp     setyp
         1        14  0.833333         3       1.0         0         0         1
    1000.0    1000.0    1000.0    1000.0       0.0       0.0       0.0         0
$Material Id = {2,{73725CD1-DA9D-465D-9E21-6FA3A00F7328}, Copper, pure, C10100, hard}
*PART
$                                                             partname   unused1
SYS-1\Surface
         2         2         3         2         0         0         0         0
*SECTION_SHELL
$       ID    elform      shrf       nip     propt        qr     icomp     setyp
         2        14  0.833333         3       1.0         0         0         1
    1000.0    1000.0    1000.0    1000.0       0.0       0.0       0.0         0
$   Keyword Snippet
$   Do not leave any unintentional empty lines in this editor
*EOS_GRUNEISEN
         2    3940.0     1.489       0.0       0.0      1.97      0.47       0.0
       1.0                   0
*ALE_STRUCTURED_MESH
         1      1000         0         0                              1.00000E16
      1001      1002         0         0         0
*ALE_STRUCTURED_MESH_CONTROL_POINTS
      1001                   1       1.0                 0.0
                   1                 0.0                0.26
                 230                60.0                0.26
                 300               120.0                 0.0
*ALE_STRUCTURED_MESH_CONTROL_POINTS
      1002                   0       1.0               -13.0
                   1                 0.0                 0.0
                1365               355.0                 0.0
*ALE_STRUCTURED_MULTI-MATERIAL_GROUP_AXISYM
Air                4         3                                               0.0
TNT                1         1                                               0.0
*ALE_STRUCTURED_MESH_VOLUME_FILLING
         1          1                            3                             0
ALL                0       0.0       0.0       0.0       0.0       0.0
*ALE_STRUCTURED_MESH_VOLUME_FILLING
         1          3                            3                             0
PART               0         1       0.0         0         0         0
*ALE_STRUCTURED_MESH_REFINE
         1         1         1         1
*ALE_STRUCTURED_FSI
         2      1000         1         1                                      -1
       0.01.00000E10       1.0       0.0                   0
*BOUNDARY_SALE_MESH_FACE
NOFLOW             1         0         0         1         0         0         0
*EOS_JWL
$       ID         A         B        R1        R2      omeg        E0        V0
         1    371.21     3.723      4.15      0.95       0.3       7.0       0.0
*EOS_LINEAR_POLYNOMIAL
         3       0.0       0.0       0.0       0.0       0.4       0.4       0.0
2.50000E-4       0.0
*INITIAL_DETONATION
         0       0.0       0.0       0.0       0.0                   0
         0       0.5       0.0       0.0       0.0                   0
         0       1.0       0.0       0.0       0.0                   0
         0       1.5       0.0       0.0       0.0                   0
         0       2.0       0.0       0.0       0.0                   0
         0       2.5       0.0       0.0       0.0                   0
         0       3.0       0.0       0.0       0.0                   0
         0       3.5       0.0       0.0       0.0                   0
         0       4.0       0.0       0.0       0.0                   0
         0       4.5       0.0       0.0       0.0                   0
         0       5.0       0.0       0.0       0.0                   0
         0       5.5       0.0       0.0       0.0                   0
         0       6.0       0.0       0.0       0.0                   0
         0       6.5       0.0       0.0       0.0                   0
         0       7.0       0.0       0.0       0.0                   0
         0       8.0       0.0       0.0       0.0                   0
         0       8.5       0.0       0.0       0.0                   0
         0       9.0       0.0       0.0       0.0                   0
         0       9.5       0.0       0.0       0.0                   0
         0      10.0       0.0       0.0       0.0                   0
         0      10.5       0.0       0.0       0.0                   0
         0      11.0       0.0       0.0       0.0                   0
         0      11.5       0.0       0.0       0.0                   0
         0      12.0       0.0       0.0       0.0                   0
         0      12.5       0.0       0.0       0.0                   0
         0       7.5       0.0       0.0       0.0                   0
*SET_NODE_LIST_TITLE
NODESET(SPC) 1
         1       0.0       0.0       0.0       0.0MECH      1                   
     63494     63496     74468     74469     74470     74471     74472     74473
     74474     74475     74476         0         0         0         0         0
*SET_MULTI-MATERIAL_GROUP_LIST
         1
3         0         0         0         0         0         0         0
*END

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值