求助：模拟土中炸药爆炸对地下结构影响

要治本

材料有ale的土壤和炸药，结构采用*MAT_PLASTIC_KINEMATIC，用的拉格朗日算法，两者用的流固耦合。
计算40us时就出现畸变。如图中畸变的那一面设置的无反射边界条件。模型厚度为100cm。

现在不知道问题是什么：1 不能用无反射边界条件？
2 炸药参数设置不合理？
3 被爆结构设置不合理？
4 计算的时候我的slave 和master换了说无法计算，这里是指slave必须是拉格朗日单元吗？

要治本

bichengcheng 发表于 2017-6-10 10:44
你这模型K文件存在的问题比较多，例如：
1.单元属性如果是拉格朗日型的，即通过*SECTION_SOLID来定义，通常 ...

谢谢您的回复，这里还有几点疑问：
1 对于您说的第三点，我看过一篇论文里面没有填充空气，有没有什么合理的方法让结构上部，或让整个结构类似于刚体不产生变形，我的参数都是和那篇论文一样的，这里应该没问题。
2 我这里爆速取得0.71cm/us,然后整个模型的宽度为1200cm。然后我在misc-pressure里查看爆炸波时，只用了80us爆炸波就到了模型边缘，这里应该也是有问题的。
再次谢谢您上面的回复，收益很多，谢谢了。

bichengcheng

你这模型K文件存在的问题比较多，例如：
1.单元属性如果是拉格朗日型的，即通过*SECTION_SOLID来定义，通常只有输入一行数据，而你却又多空了一行，这是不合理的，通常情况下计算不了；
2.你采用*MAT_SOIL_AND_FOAM来定义土体，体积应变与压力之间的关系是成对出现的，通常为10组关键参数，而你只给出了压力参数。所以建议要么都不给，要么给全。
3.你的模型中结构为框架结构，中间没有任何填充物，由于炸药距离结构物上表面较近，结构物受到冲击波的作用发生变形，ALE网格与Lagrange网格出现没有重叠的部分，因此该部位流固耦合失效，产生不可控制的计算错误，为此建议该部分用空气填充
4.流固耦合关键字*CONSTRAINED_LAGRANGE_IN_SOLID第二行第2个关键字表示耦合失效时间，很多例子给的都是0，个人觉得还是采用默认值1E10比较合理。
以上四点均是本人做了大量的数值模拟得出的浅薄的经验之谈，供你参考。

要治本

部分K文件，第一次发，不太会弄。求大神不吝赐教
*SECTION_SOLID
         1         1
         
*SECTION_SOLID_ALE
         2        11
         
*SECTION_SOLID_ALE
         3        11
         
$
$
*ALE_MULTI-MATERIAL_GROUP
         2         1
         3         1  
*CONSTRAINED_LAGRANGE_IN_SOLID
         2         1         0         0         3         5         3         0
         0         0         0         0         
*SET_PART_LIST
         1
         2         
*SET_PART_LIST
         2
         1
*CONTROL_ALE
         2         1         2        -1         0         0         0         0   
  0.00       0.1E+21         1         0         0         0   
$      
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             MATERIAL DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*MAT_PLASTIC_KINEMATIC
         1       2.5       0.3      0.21   2.01E-4      3E-3       0.1

*MAT_SOIL_AND_FOAM
         2       1.8   6.38e-4   1.328e2  3.33E-03   7.03e-7       0.3         0
         0         0
         
         
       0.0 3.420E-02 4.530E-02 6.760E-02 1.270E-01 2.080E-01 2.710E-01 3.920E-01
5.660E-01 1.230E+00              
*MAT_HIGH_EXPLOSIVE_BURN
         3  0.163E+1      0.71      0.21       0.0       0.0       0.0
*EOS_JWL
         1      3.74     0.074      4.15      0.95       0.3     0.043       1.0
*INITIAL_DETONATION
1,0,1320,0,0