内爆炸数值模拟结果与实际不符

ycy_0517

我对大空间钢结构内爆炸情况进行数值模拟，ALE空间为62米直径，21米高的圆柱体，边界条件设定为无反射边界条件，技算结果与实际超不符，提取了炸药上方一点的压力时程曲线（附件里面，不知道怎么直接粘图片 ），不知道什么原因，是边界条件不对还是有什么东西没考虑，希望有经验的你们给我点意见，谢谢啊

wyc412721

是有点问题，刚才始超压就很大，楼主能否介绍下建模思路

ycy_0517

#[/url] wyc412721
你好，我的结构形式是双层球壳结构，壳体杆件选用LINK160，地面和墙壁用SHELL163,炸药和空气都是SOLID164单元，用ALE法，边界条件采用无反射边界条件，我昨天晚上回去想了一下，觉得自己犯了个错误，计算终止时间设的太短了，今天设为1S之后整个趋势还是一样的，部分K文件如下：$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*SECTION_BEAM
1    ,     3
0.04   
*SECTION_SHELL
2     ,    2  ,  1.0,  2.0 ,   0.0 ,   0.0  ,   0
0.001  ,0.001 , 0.001 ,  0.001 ,     0.00   
*SECTION_SOLID_ALE
         3   , 11
         0,0,0,0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             MATERIAL DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*MAT_NULL
1,1.29,0,0,0,0,0,0
*EOS_LINEAR_POLYNOMIAL
$    EOSID        C0        C1        C2        C3        C4        C5        C6
        1,  0.000E+00, 0.000E+00 ,0.000E+00, 0.000E+00 ,    0.400,     0.400 ,0.000E+00
$       E0        V0
2.5e5,1   
*MAT_HIGH_EXPLOSIVE_BURN
$      MID        RO         D       PCJ      BETA         K         G      SIGY
2,1.6310E+3,0.6718E+4,1e10,0.0000000,0.0000000,0.0000000,0.0000000
*EOS_JWL
$    EOSID         A         B        R1        R2      OMEG        E0        V0
2,5.409E+11,0.094E+11,4.5000000,1.1000000,0.3500000,0.080E+11,1.0000000
*MAT_rigid
$      MID        RO         E        PR      SIGY      ETAN      BETA
3,7830,2.07e30,0.300,0,0.000E+00,0.000E+00
1,7,7
0,0,0,0,0,0
*ALE_MULTI-MATERIAL_GROUP
3,1
4,1
*INITIAL_DETONATION
$      PID         X         Y         Z        LT
4,0.0,0.0,19.0,0
*CONSTRAINED_LAGRANGE_IN_SOLID
$    SLAVE    MASTER     SSTYP     MSTYP     NQUAD     CTYPE     DIREC     MCOUP
       1,       2,        0  ,      0,         3 ,        4   ,      2,         1
$    START       END      PFAC      FRIC    FRCMIN      NORM
       0.0 ,     0.0  ,    0.0   ,    0.0  ,    0.3  ,    0   
$       CQ      HMIN      HMAX     ILEAK
       0.0                             1
*SET_PART_LIST
1
1,2
*SET_PART_LIST
2
3,4   
$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              PARTS DEFINITIONS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
*PART
Part          1 for Mat         3 and Elem Type         2
         1         1         3         0         0         0         0
$
*PART
Part          2 for Mat         3 and Elem Type         3
         2         2         3         0         0         0         0
$
*PART
Part          3 for Mat         1 and Elem Type         1
         3         3         1         1         0         0         0
$
*PART
Part          4 for Mat         2 and Elem Type         1
         4         3         2         2         0         0         0
$
*CONSTRAINED_RIGID_BODIES
1,2
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             ELEMENT DEFINITIONS                              $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*ELEMENT_BEAM
       1       1       1       2       0
       2       1       1       3       0
............................................................................

     30637     30638     30639     30640     30641     30642     30643     30644
     30645     30646     30647     30648     30649
*BOUNDARY_SPC_SET
         1         0         0         1         0         1         0         1
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              CONTACT DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              CONTROL OPTIONS                                 $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*CONTROL_ENERGY
         2         2         2         2
*CONTROL_SHELL
  20.0             1        -1         1         2         2         1
*CONTROL_ALE
         3         1         2 -1.00      0.00      0.00      0.00      0.00   
  0.00     0.100E+21  1.00      0.00      0.00             0
*CONTROL_TIMESTEP
    0.0000    0.4000         0  0.00      0.00   
*CONTROL_TERMINATION
0.800E-01         0   0.00000   0.00000   0.00000
*DATABASE_HISTORY_NODE
         1
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                                 TIME HISTORY                                 $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_GLSTAT  
0.4E-04
*DATABASE_BINARY_D3PLOT
0.4E-04
*DATABASE_BINARY_D3THDT
0.4E-04
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              DATABASE OPTIONS                                $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_EXTENT_BINARY
         0         0         3         1         0         0         0         0
         0         0         4         0         0         0
*END

ycy_0517

还有就是我屋面没有设阻挡结构，以壳顶点为零点坐标往上建立结构，就整个结构倒过来了，这应该不会影响模拟结果吧

jyqqcm

个人觉得只要结构的相对位置对了，以何处为零坐标应该不会影响结果吧，
可能是其他原因，你再找找