load_blast和ale的差别

jiayixia

求解，为什么我定义的Load-blast在message中所提供的压力值和使用ale方法得到的压力值相差那么多呢~求解
炸药距离板1.0m处，炸药约5kg，我使用load-blast得到的最大超压为20.53Mpa，但是我使用ale方法得到炸药投影点为超压峰值只有3.5Mpa左右。但是两者板的响应又差不多一样大，求解答。我用亨里奇公式算出来，比距离等于0.585时，炸药的超压峰值为2.02MPa，再加上反射压力应为入射压力的2~8倍，所以我觉得load-blast所提供的炸药超压是正确的。但是ale为什么错呢，不解

东方之子

楼主的k文件能否上传一下，让大家研究一下

jiayixia

好的，不过ALe的K文件太大了，我上传关键字吧
*initial_volume_fraction_geometry
1,1,1
6,0,2
0.0,20.0,0.0,9
*ALE_MULTI-MATERIAL_GROUP
1,1
4,1
*CONSTRAINED_lagrange_in_SOLID
2,1,1,1,3,4,2,0
0,1e10,0.1,0,0.$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             SECTION DEFINITIONS                              $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*SECTION_SOLID_ale
1,11
*SECTION_SHELL
         2         2    0.8333       5.0       0.0       0.0         0
0.900     0.900     0.900     0.900      0.00   
*SECTION_SHELL
         3         2    0.8333       5.0       0.0       0.0         0
0.650     0.650     0.650     0.650      0.00   
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             MATERIAL DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*MAT_NULL
         1 0.120E-02  0.00      0.00      0.00      0.00      0.00      0.00   
*EOS_LINEAR_POLYNOMIAL
         1  0.00      0.00      0.00      0.00     0.400     0.400      0.00   
  0.00      0.00   
*MAT_PLASTIC_KINEMATIC
         2  7.85      2.06      0.300000 0.345E-02 0.618E-01  0.00   
0.400E-04  5.00     0.500E-01
*MAT_HIGH_EXPLOSIVE_BURN
3,1.63,0.693,0.21,0.0,0,0,0
*EOS_JWL
3,3.7400002,0.0737,4.1500001,0.95,0.3,0.043,1.00
*INITIAL_DETONATION
1,0,20,0,0.0
  
$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              PARTS DEFINITIONS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
*PART
Part          1 for Mat         1 and Elem Type         1
         1         1         1         1         0         0         0
$
*PART
Part          2 for Mat         2 and Elem Type         2
         2         2         2         0         0         0         0
$
*PART
Part          3 for Mat         2 and Elem Type         2
         3         3         2         0         0         0         0
*PART
Part          3 for Mat         3 and Elem Type         1
         4         1         3         3         0         0         03,1
*CONSTRAINED_lagrange_in_SOLID
3,1,1,1,3,4,2,0
0,1e10,0.1,0,0.3,1
[/quote]回复 2# 东方之子

蓝牙

压力取值要取ouler网格上的压力，直接取板上压力是不对的！

jiayixia

我知道的，就是欧拉网格上的压力和load-blast算出来差别好大回复 4# 蓝牙

蓝牙

回复 5# jiayixia
我曾经比对过炸药在水中的压力峰值，采用ALE和经验公式对比差别最大的也不超过20%
你可以尝试
1.把板拿掉，单独比较炸药的压力传播
2.传播介质（如水或者空气）的状态方程你要确保没有问题，这个你可以参考dynalook上的相关文章的参数
3.网格尺寸不要太粗糙，ALE对网格要求比较高，尽可能的把网格划细些
4炸药的形状对计算结果有一定的影响，我的经验在较远的地方炸药的形状对计算结果的影响不是很大，但是在近处，影响还是较大的
5.blast一般要比ALE精确，但是blast的使用受到诸多限制，如果没有应力波的绕行和反射你就采用blast得了，精确并且工作量不大，我一直用它！个人之言

jiayixia

谢谢斑竹的回复，我得到的结论是，欧拉网格的压力是一个单元的平均应力，而load-blast的到的压力是真正作用在板上的压力，这两个还是有所不同的，真正作用到板上的压力要比欧拉单元的压力大很多，这个可以通过经验公式和TM5-1300的比较就知道了回复 6# 蓝牙

蓝牙

回复 7# jiayixia
你可以采用*database_FSI进行实际作用压力 力的检测
也可以在板上布置压力传感器 *database_FSI_sensor进行压力的校对

jiayixia

这个高深了，求教楼主怎么做回复 8# 蓝牙