采用cel计算空爆，炸药穿透平板

41439937 · 发表于 2015-5-20 10:20:04

本帖最后由 41439937 于 2019-1-7 15:17 编辑

zho123855 · 发表于 2015-5-20 10:20:05

本帖最后由 zho123855 于 2015-6-3 22:28 编辑

41439937 发表于 2015-5-20 10:23
空气材料
*Material, name=GAS
*Density

强烈关注你的这个问题，我的经验是：尝试把Euler网格画得细一些。流体穿透流固耦合边界是CEL的顽疾，或者更准确的说，是Euler相关算法的通病，它源自Euler算法对于边界的处理（Euler算法没法模拟出明确的流固耦合边界，只能是Fraction）。
另，请教：你选择了理想气体作为其本构模型，其合理性你是否考察过？根据Abaqus的帮助文档，在CEL分析中要谨慎地使用理想气体模型（只要有一小点真空（void）出现，真空就会不断地扩展、造成计算结果完全不物理）。欢迎交流CEL空气相关的问题。再补充几句，Abaqus在结构动力学、计算固体力学方面的能力是毋庸置疑的，但对流体的关注也就是最近几年才开始的（我记得CFD模块加进来没几年），CEL算法可以解决一部分流固耦合问题，但我一直觉得好像Abaqus并没有打算把CEL作为一个强有力的流固耦合工具，最明显的就是CEL里面没有流体力学里面经典的本构关系（CEL里面都是用EOS来表达流体的本构关系的），也没有湍流模型。不知道大家对于这个问题是怎么看的。欢迎讨论啊。

41439937 · 发表于 2015-5-20 10:23:49

本帖最后由 41439937 于 2015-5-27 14:52 编辑

空气材料
*Material, name=GAS
*Density
1.225,
*Eos, type=IDEALGAS
280,100000.
*Specific Heat
717.6,288.
炸药材料
*Material, name=TNT
*Density
1630.,
*Eos, type=JWL
6930., 3.737e+11, 3.7471e+09, 0.35, 4.15, 0.9, 5.99999e+09, 0.
*Detonation Point
0., 0., -0.15, 0.
接触属性
*Surface Interaction, name=INTPROP-1
*Friction
0.1,
*Surface Behavior, pressure-overclosure=HARD
接触设置
*Contact, op=NEW
*Contact Inclusions, ALL EXTERIOR
*Contact Property Assignment
, , INTPROP-1

41439937 · 发表于 2015-5-20 14:53:24

求斑竹帮忙看看啊，

41439937 · 发表于 2015-5-27 14:52:41

没人看看吗？

leo0072008 · 发表于 2015-6-1 20:08:15

本帖最后由 leo0072008 于 2015-6-1 20:16 编辑

5.99999e+09是不是量级错了？？？
个人觉得应该是2~10之间这个量级，inp能看看吗？leo2002leo2002@sina.com

41439937 · 发表于 2015-6-2 21:59:25

leo0072008 发表于 2015-6-1 20:08
5.99999e+09是不是量级错了？？？
个人觉得应该是2~10之间这个量级，inp能看看吗？leo2002leo2002@sina.com ...

之前发的这个是有问题，e0是3.6e6左右，后来改过。穿透问题的话，主要是网格匹配这块。

半数花开 · 发表于 2015-6-2 23:30:55

楼主上传inp 我们看看怎么回事吧这样光看也不容易看出是哪里的问题啊

41439937 · 发表于 2015-6-3 19:25:12

本帖最后由 41439937 于 2015-6-3 19:36 编辑

半数花开发表于 2015-6-2 23:30
楼主上传inp 我们看看怎么回事吧这样光看也不容易看出是哪里的问题啊

后来问题解决得差不多了，穿透是由于欧拉网格太大，规定是在拉格朗日欧拉耦合处差不多板厚方向要3到5个欧拉网格，当占比小于0.5的时候。会不计入耦合，所以产生了穿透。不过总的计算结过还是比dyna与autodyn算得差。耦合参考这个资料http://wenku.baidu.com/view/66b6 ... e?qq-pf-to=pcqq.c2c

半数花开 · 发表于 2015-6-3 19:37:15

41439937 发表于 2015-6-3 19:25
后来问题解决得差不多了，穿透是由于欧拉网格太大，规定是在拉格朗日欧拉耦合处差不多板厚方向要3到5个欧 ...

那也就是说楼主觉得算爆炸这块还是dyna要好用一点吗现在我们基本上所有爆炸的都是用的dyna，没有用abaqus的。。。

41439937 · 发表于 2015-6-3 19:50:21

半数花开发表于 2015-6-3 19:37
那也就是说楼主觉得算爆炸这块还是dyna要好用一点吗现在我们基本上所有爆炸的都是用的dyna，没有用abaq ...

ale法算的话还是dyna和autodyn好吧，不过abaqus可以用声固耦合和conwep计算水爆和空爆。这个看情况吧。而且也可能是我自己对abaqus的cel算法了解不深，所以没算好。

半数花开 · 发表于 2015-6-3 20:04:06

41439937 发表于 2015-6-3 19:50
ale法算的话还是dyna和autodyn好吧，不过abaqus可以用声固耦合和conwep计算水爆和空爆。这个看情况吧。而 ...

谢谢楼主的回答感谢了！

41439937 · 发表于 2015-6-4 15:03:53

本帖最后由 41439937 于 2015-6-4 15:08 编辑

zho123855 发表于 2015-6-3 22:20
强烈关注你的这个问题，我的经验是：尝试把Euler网格画得细一些。流体穿透流固耦合边界是CEL的顽疾，或者 ...

因为我不是主要搞流体的，这块采用理想气体，是参照别人的论文来给的。材料穿透的话，考虑加密耦合处的网格密度会有很大的改善，推荐厚度方向流体比结构位3-5.多了也出现了一些穿透。。。这个我自己也弄不明白，按理来说应该是网格越密的话，约不会出现Fraction小于0.5的情况的。

41439937 · 发表于 2015-6-4 15:05:58

41439937 发表于 2015-6-4 15:03
因为我不是主要搞流体的，这块采用理想气体，是参照别人的论文来给的。材料穿透的话，考虑加密耦合处的网 ...

另外新版本有流域自适应网格的划分。可以自动加密耦合处的网格，但是我试了一下，并没有什么卵用、、、、该穿透还是穿透，手动加密过后对穿透问题的改善特别大、

leo0072008 · 发表于 2015-6-5 11:29:24

楼主方便发一下input文件吗？大家一起交流学习一下，leo2002leo2002@sina.com

41439937 · 发表于 2015-6-5 14:13:33

leo0072008 发表于 2015-6-5 11:29
楼主方便发一下input文件吗？大家一起交流学习一下，

发你邮箱了

johnnymanson · 发表于 2015-6-19 03:31:19

楼主方便发一下input文件吗？sanyan6868@hotmail.com

41439937 · 发表于 2015-6-19 10:44:54

本帖最后由 41439937 于 2015-6-19 10:47 编辑

johnnymanson 发表于 2015-6-19 03:31
楼主方便发一下input文件吗？

给你发了，但是又被退回来了。。。。好像是你的邮箱不支持。。。不好意思。

johnnymanson · 发表于 2015-6-19 14:30:27

41439937 发表于 2015-6-19 10:44
给你发了，但是又被退回来了。。。。好像是你的邮箱不支持。。。不好意思。 ...

请再发一次感谢
johnnymanson68@qq.com

lmz_714 · 发表于 2015-10-24 15:25:58

楼主，你好！我也在一直在用CEL，最近老是碰到穿透问题！根据你的解答是否可以这样理解：在耦合处欧拉网格尺寸比拉格朗日网格尺寸小3-5倍就可有效改善穿透现象。还是说欧拉网格尺寸必须小到一定值？

账号		自动登录	找回密码
密码			注册

账号		自动登录	找回密码
密码			立即注册

[流固耦合] 采用cel计算空爆，炸药穿透平板

最佳答案