流固耦合的耦合方式选择（即CTYPE）问题

songhaixian

流固耦合中的耦合方式：常用的是默认情况下的CTYPE=2耦合方式，如果拉格朗日实体中含有刚体时则须采用CTYPE=4或CTYPE=5（考虑侵蚀的solid）的耦合方式。
本人对上述三种耦合方式分别进行了验证，但是结果却让我很不理解，首先介绍一下我的模型：如图
part1（红色）为刚性板（作为地基）；part2（浅蓝）为混凝土框架；part3（绿色）为空气；part4（黄色）为炸药，PART1、PART2之间采用了自动面面接触，并给框架施加了重力！此为1/2模型，对称面为XOZ
由于重点在于分析炸药对框架的影响，因此流固耦合中的拉格朗日实体可以不包含PART1即刚性板（个人理解），因此首先我采用了CTYPE=2进行耦合，计算初期便会有错误提示：
*** Error 5661 (SOL+661)
     ALE FSI coupling type 1,2,3 cannot couple with rigid
     body. Please check *CONSTRAINED_LAGRANGE_IN_SOLID card.

虽有此错误提示，但是程序仍然可以计算，而且最终的计算结果也较为满意，中柱在炸药起爆后便破坏了！
因为有上面的错误提示，随后我又对CTYPE=4和CTYPE=5进行了计算，此时虽然没有了错误提示，然而这两者的计算结果出乎我的意料,特别是CTYPE=5,因为其考虑了材料的侵蚀，在爆炸荷载那么大的作用下没有理由不发生破坏，而CTYPE=4毕竟没有考虑侵蚀，所以不破坏倒是可以理解！
流固耦合关键字：
*ALE_MULTI-MATERIAL_GROUP
3,1
4,1
*SET_PART_LIST
9,0,0,0,0
2
*SET_PART_LIST
10,0,0,0,0
3,4
*CONSTRAINED_LAGRANGE_IN_SOLID
9,10,0,0,3,5,2

$
不知道上述三种耦合方式应该作何选择，哪种选择才是正确的，还请各位指教！

问题补充：
由于存在刚体的影响，我又将刚体去除，只剩下PART1,3,4，此时仍用CTYPE=2，所得结果与存在刚体时的CTYPE=2几乎相同！

tjuabraham

没法回答你的问题。但是对你ctype=2算法有些质疑。流固耦合中为什么没有地面刚体？冲击波和地面之间的相互作用如何体现？

songhaixian

tjuabraham 发表于 2012-2-29 22:25
没法回答你的问题。但是对你ctype=2算法有些质疑。流固耦合中为什么没有地面刚体？冲击波和地面之间的相互 ...

因为我不关心地面刚体的受力，所以我就没有将地面刚体放在拉格朗日实体中...

tjuabraham

songhaixian 发表于 2012-2-29 22:32
因为我不关心地面刚体的受力，所以我就没有将地面刚体放在拉格朗日实体中... ...

流固耦合两种网格能够相互交叉，感觉还是应该考虑上刚性地面的反射作用，把刚性地面加到拉格朗日set里面。刚性地面对冲击波有放大作用。

songhaixian

tjuabraham 发表于 2012-2-29 22:36
流固耦合两种网格能够相互交叉，感觉还是应该考虑上刚性地面的反射作用，把刚性地面加到拉格朗日set里面 ...

嗯，你说的有道理，多谢，但是好像还是不能解决这几种耦合方式之间的差异！

chenyuanjp

在prepost里有有耦合点个数的选项nquad 一般被耦合的网格要细一些 但是增加耦合点会使计算不易安定 一般设置2 or 3 论坛里这方面的讨论很多 你search一下

songhaixian

还请各位大侠指教啊！有人做过相关的对比研究嘛？

chenyuanjp

个人觉得 用2 因为没有和刚体耦合，压力在底部会有泄露。2型在耦合过程里考虑了节点的质量 所以耦合力显得比较大。一般在构造解析里用得比较多。结果正确不正确还要看压力变化过程。个人认为你的算例还是应该用4 或5 （罚函数法） 增加耦合点个数 或细化网格 以加大耦合效果

songhaixian

chenyuanjp 发表于 2012-3-1 10:33
个人觉得 用2 因为没有和刚体耦合，压力在底部会有泄露。2型在耦合过程里考虑了节点的质量 所以耦合力显得 ...

首先非常感谢您的回答，我的网格尺寸采用100mm的，网格的大小应该不会有太大问题。
请问对于增加耦合点个数的有什么相关说明吗？如何控制耦合点的个数呢？如果按您所说的增加耦合点个数可以增大耦合效果，那么如果耦合点个数控制不准确必然会引起结果失真！还请多多指教！

songhaixian

chenyuanjp 发表于 2012-3-2 06:39
在prepost里有有耦合点个数的选项nquad 一般被耦合的网格要细一些 但是增加耦合点会使计算不易安定 一般设 ...

对CTYPE的不同取值，我的耦合点个数都是设置为3，即*CONSTRAINED_LAGRANGE_IN_SOLID中的第5个参数取值

32111258

:victory:

xionganping1985

很值得关注的问题，在过半个月才会重新开始做这个仿真，到时希望能够和您交流下！

吾艾曼

对于  CTYP=5  时要与*MAT_ADD_EROSION  联合起来用  才会出现材料侵蚀

songhaixian

吾艾曼 发表于 2012-3-16 10:38
对于  CTYP=5  时要与*MAT_ADD_EROSION  联合起来用  才会出现材料侵蚀

材料的侵蚀已经定义了！

风浪雨奏

用2的时候爆炸应力特别大，用5特别小时什么原因