lsdyna爆破模拟求助~炸药单元最外一层网格总是会过度拉伸,然后岩体突然消失,请问怎么回事?
建模比较简单,就是一个长方体,用工作面且出来炸药,空气间隔和堵塞的部分,用的ms-mm-g-MPa-N的单位制
材料,堵塞是*MAT_SOIL_AND_FOAM,岩体用的是*MAT_PLASTIC_KINEMATIC
再贴出来主要的k文件
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ SECTION DEFINITIONS $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*SECTION_SOLID
1 2
*SECTION_SOLID
2 11
*SECTION_SOLID
3 11
*SECTION_SOLID
4 11
*ALE_MULTI-MATERIAL_GROUP
2 1
3 1
4 1
*CONTROL_ALE
2 1 2 -1.00 0.00 0.00 0.00 0.00
0.00 100 0
*CONSTRAINED_LAGRANGE_IN_SOLID
1 2 0 0 0 5 3 0
0 0 0.1 0.00 0.2 0 0 0.0
0.0 0 0 0 0.0 0 0 0
*SET_PART_LIST
1
1
*SET_PART_LIST
2
2 3 4
*HOURGLASS
1 1 0.1 0 2.0 0
*HOURGLASS
2 1 1.0E-6 0 2.0 0
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ MATERIAL DEFINITIONS $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*MAT_PLASTIC_KINEMATIC
$ MID RO E PR SIGY(屈服) ETAN BETA
1 2.27E-3 45000 0.25 60 45 0
1 3 0
*MAT_NULL
2 1.293E-6 -1.0 1.81E-8 0.00 0.00 0.00 0.00
*EOS_LINEAR_POLYNOMIAL
2 0.00 0.00 0.00 0.00 0.40 0.40 0.00
254374E-6 0.75
*MAT_HIGH_EXPLOSIVE_BURN
3 1.1E-3 3.2E3 5.733E3 0.00 0.00 0.00 0.00
*EOS_JWL
3 3.71E5 3230.6 4.150 0.95 0.30 7000.0 1.00
*INITIAL_DETONATION
0 0 0 7500.0 0.0
*MAT_SOIL_AND_FOAM
4 2E-3 16 10E7 0 0 133.3 0.0
0.0 0.0
0.0 0.050 0.090 0.11 0.15 0.19 0.21 0.22
0.25 0.30
0.0 3.420E-2 4.530E-2 6.760E-2 1.270E-1 2.080E-1 2.710E-1 3.920E-1
5.660E-1 1.230
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ PARTS DEFINITIONS $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
*PART
Part 1 for Mat 1 and Elem Type 1 这是岩体
1 1 1 0 2 0 0
$
*PART
Part 2 for Mat 2 and Elem Type 1 这是空气间隔
2 2 2 2 1 0 0
$
*PART
Part 3 for Mat 3 and Elem Type 1 这是炸药
3 3 3 3 1 0 0
$
*PART
Part 4 for Mat 4 and Elem Type 1 这是堵塞
4 4 4 0 1 0 0
如图所示中间明显被拉长的两列网格是炸药
再下一步,就只剩下中间的炸药,空气,填塞了岩体就消失了
不知道是材料 的问题,建模的问题,还是流固耦合的问题,求前辈们指导~
另外,还有一些其他问题,想请教资深大神们,可以有偿,但是还是学生没有收入,没有太多钱,或者可以和其他同学们一起交流,我的wx号是:MRC41312032~~