舱内爆炸流固耦合的问题

xiayewangzhe

计算方舱内部爆炸，舱壁设了五个part,都是用拉格朗日，空气和炸药（图中方框线和圆柱）为多物质ALE，采用流固耦合算法。计算一开始正常炸药起爆舱壁开始变形，到2800us时出现错误有舱壁上的负体积和空气对称面上的node超量程。求大神指导。之前PFAC设置为1，这个问题出现在400us时候，后来改成0.15就推迟到2800us出现。

*SECTION_SOLID
         1         1
*SECTION_SOLID
         2         1
*SECTION_SOLID
         3         1
*SECTION_SOLID
         4         1
*SECTION_SOLID
         5         1
*SECTION_SOLID_ALE
         6        11
         
*SECTION_SOLID_ALE
         7        11
         
*CONTROL_ALE
         2         1         2        -1
         
*ALE_MULTI-MATERIAL_GROUP
         6         1
         7         1
*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   NORMTYP      DAMP
                          0.15                  
$       CQ      HMIN      HMAX     ILEAK     PLEAK   LCIDPOR                              
                                                  
*SET_PART_LIST
         1
         1         2         3         4         5                  
*SET_PART_LIST
         2
         6         7
$
*CONTROL_CPU
  0.00   
*CONTROL_ENERGY
         2         1         2         2
*CONTROL_SHELL
  20.0             1        -1         1         2         2         1
*CONTROL_BULK_VISCOSITY
  1.50     0.600E-01
*CONTROL_TIMESTEP
    0.0000    0.2000         0  0.00      0.00   
*CONTROL_TERMINATION
0.300E+05         0   0.00000   0.00000   0.00000
$

xiayewangzhe

补充：所有模型均为实体164单元，出现负体积的部分全部在上面舱壁上附图如下

xiayewangzhe

疑惑，是不是网格的问题？

chenyuanjp

看看 你的炸药的物性单元统一吗

xiayewangzhe

非常感谢大侠回复，您说的物性单元指的是什么？小弟不太明白，一下是炸药参数
*MAT_HIGH_EXPLOSIVE_BURN
6，1.7,0.693，0.21
*EOS_JWL
6，3.712，3.75E-02，4.15，0.9，0.35，0.04853，1.0     
炸药网格是固体网格的1.5倍左右，空气网格是炸药的三倍左右。空气域太大了，改细了K文件太大

护韶

xiayewangzhe 发表于 2013-5-10 08:25
非常感谢大侠回复，您说的物性单元指的是什么？小弟不太明白，一下是炸药参数
*MAT_HIGH_EXPLOSIVE_BURN
6 ...

你好，我想请教，流固耦合多物质耦合方法的时候，1炸药和空气的网格是不是各自划分没直接关系？？？2还是炸药和空气网格要重合？？？？？3还是炸药和空气网格要共节点？？？？？？？？

xiayewangzhe

护韶 发表于 2013-6-26 10:43
你好，我想请教，流固耦合多物质耦合方法的时候，1炸药和空气的网格是不是各自划分没直接关系？？？2还是 ...

在建模上直接切除炸药就可以，是共节点的。在K文件里设置ALE选项就行了。见上面我的K 文件。上面的K文件有问题，流固耦合的MCOUP参数应该改为0不然就没有耦合空气，看不到冲击波的传播。

护韶

xiayewangzhe 发表于 2013-7-16 12:46
在建模上直接切除炸药就可以，是共节点的。在K文件里设置ALE选项就行了。见上面我的K 文件。上面的K文件 ...

您好，我的K文件里面MCOUP参数是设为0了。然后我想问怎么看冲击波的传播？？？是在后处理器中怎么操作？？？

caoxu666

xiayewangzhe 发表于 2013-5-9 18:51
补充：所有模型均为实体164单元，出现负体积的部分全部在上面舱壁上附图如下 ...

请教一下，我的毕业课题跟你的一样也是研究舱室内爆，爆炸效果一切正常，但是舱室没有发生变形怎么回事，我的单位也是设的g cm 微秒，舱体材料
*MAT_PLASTIC_KINEMATIC
3,7.85,1.98,0.3,0.28,0.136,0.5
40.0,5.00,0.5
还是需要在加什么关键字？

wangzelongsd

caoxu666 发表于 2016-8-5 16:24
请教一下，我的毕业课题跟你的一样也是研究舱室内爆，爆炸效果一切正常，但是舱室没有发生变形怎么回事， ...

硕士？可以一起讨论讨论

caoxu666

wangzelongsd 发表于 2016-8-26 00:43
硕士？可以一起讨论讨论

可以qq305086278