[吐血奉献]流固耦合分析成功的条件

zzw420

在命令流中，设置fluid 141（142）单元之后，加一句keyopt,1（注意：是单元编号，不是单元号）,4,1.就可以了。只有打开这个，才能让141（142）单元实现ALE的功能。
本人一般都用命令流，所以在划分网格时候都是直接指定划分的单元类型和材料类型。至于材料的属性，只要是在/solu之前就可以。所以我说不需要分顺序。

eversoon

那固体的设定单元以后，还用keyopt命令吗？因为我的固体单元编号是2，
所以是不是也得设置keyopt，2，4，1？
   还有一个问题是你在划分流体网格的时候指定的材料类型是什么？
   因为一般是不定义流体的材料类型的，只有固体才有定义材料类型，而流体的只是在Fluid Property里面设定。

zzw420

固体单元似乎没有ALE这个说法。一般来说，划分流体网格的时候，都默认流体为1号材料。不信你试试看，改为2号就会出错。

eversoon

谢谢楼主的细心解答^_^！
我现在可以算通了，但是我是在材料属性里面另外添加了一个材料编号，设置了流体的材料特性（密度和粘度）。然后在流体网格划分的时候材料项就选择了流体的材料编号2，不知道这样会不会和在Fluid Property中设置的密度和粘度重合。

zzw420

可以去试试看。从帮助中的例子（活塞）看，似乎它在流场网格划分中就根本没有指定材料号。固体材料编号为1。看来也是可以的！我认为，流体网格划分最好不要指定材料号，ansys应该会按照你的流体单元去找流体属性的。

mbrslydddd

看了那么多网友热情的讨论，我也从中受益匪浅。个人对于使用ansys做流固耦合也只是摸索着前进，感觉很麻烦，因为我不是流体专业的，我是固体力学方面的，主要做弹塑性方面和热、结构方面的分析，对流体力学了解的不是太多，且ansys帮助文档上的说明也不怎么合意。近来看到ansys10宣传说她们将在新版本中改善流固耦合的分析的特点，将在7月份上市，强烈期待中...。
－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－
上面好多网友讨论材料属性的问题，这个问题是这样解决的，若使用命令流的方式就没有什么好说的了，若采用GUI方式，若你同时定义了fluid单元和固体单元，且定义了固体单元的材料属性，这时若划分流体网格时，固体材料属性也显示在单元属性里面，我告诉你这时此固体单元属性对流体单元是无任何意义的，即使说你无论定义与否，都不会影响流体单元的属性，流体单元的属性只有在fluid Property里面定义才有效，且流体单元本身也不具备弹性模量、泊松比和屈服强度等等固体的材料属性，所以这个问题大家不要担心。
－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－
流固耦合分析中流体单元的编号必为一，具体情况可以参考帮助文档的相关章节。
－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－
我说的流固耦合分析顺序基本上是按照帮助文档上的顺序，这样做比较不容易出错，若顺序改变，应该也可以的。
－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－
ALE方法既可以用于流体分析又可以用于固体分析，这点在ansys里面体现的不是太好，（不知道我说的对不对），若有人用过ls_dyna就会知道了，再做固体的大变形分析时
（例如高速撞击、压力加工等）ALE方法是一个非常有用的方法。大家可以参考相应的理论书籍看看Lagrange、Euler和ALE方法的区别。我介绍一本书，
《连续体和结构的非线性有限元》 Ted Belytschko, Wing Kam Liu, Brian Moran 庄茁 译。 清华大学出版社。

mbrslydddd

流固耦合分析中固体推动流体还是流体推动固体的差别应该是很大的。不说设置方面，单说建模方面，若是流体推动固体的分析，那么你的模型中必须有固体模型，即使固体属于刚体，否则你怎么定义固体的材料属性。而固体推动流体则相对较为好处理些，若固体处于流体的内部，且固体属于刚体的话，那么固体模型就可以省去，此时流体耦合就可以变成单独的流体分析。因此谁先谁后差别还是很大的。

mbrslydddd

哪位有流体推动固体，而固体不受约束，发生大范围移动的例子，小弟做了几次结果不是很理想，若哪位网友有能否传上来，

zzw420

mbrslydddd wrote:
哪位有流体推动固体，而固体不受约束，发生大范围移动的例子，小弟做了几次结果不是很理想，若哪位网友有能否传上来，

我给楼主一点建议，做这种例子，一定要注意载荷的施加方式。否则很容易发生单元负体积，或者结果不收敛的情况！

eversoon

ANSYS的例子中的那个活塞的例子就是属于固体推动流体的流固耦合例子，不过全是命令流的形式。没有GUI的。我将没个命令用gui形式的做了，基本上还能做出来。

eversoon

楼上的朋友，一般的当发生单元负体积的时候，你是怎么解决这个问题的？
丛哪开始着手去寻找问题的根源？

mbrslydddd

zzw420 wrote:
我给楼主一点建议，做这种例子，一定要注意载荷的施加方式。否则很容易发生单元负体积，或者结果不收敛的情况！

我这样的问题到没有怎么出现，一些很奇怪的问题到是出现了，看样zzw420老兄有这方面的经验了，有没有现成的算例可供参考，小弟急需。

zzw420

对于负体积的问题，可以从很多方面找原因：边界条件，载荷施加方式，单元划分，时间步长，remesh参数，…
楼主不妨留一个联系方式啊！我也是自己摸索罢了，可以共同探讨！

mbrslydddd

QQ:379296331
e-mail:mbrslyd@163.com

eversoon

zzw420 wrote:
对于负体积的问题，可以从很多方面找原因：边界条件，载荷施加方式，单元划分，时间步长，remesh参数，…

你说的remesh的参数设置，不知道具体怎么说。该如何设置呢？

ztao03

先试试

gus1977

先顶一下再看

xujiazhu

不知仁兄mbrslydddd能不能上传一下实例呢？小弟不胜感谢