声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2653|回复: 0

[综合讨论] ANSYS流固耦合分析成功的条件【转】

[复制链接]
发表于 2014-5-25 10:57 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
1.首先在建模和条件设置方面要按照这样的设置顺序:
  (1)选取流体单元,(打开keyopt(4)选项),建立流体模型,注意此处挖去固体所占的空间,然后分区划分流体场网格(好像在ls_dyna里面不要挖去固体所占空间),注意靠近挖去空间的部分网格应该细小些,还有若要采用remesh在计算中重划网格,一定要使用三角形单元(所有流体场)
  (2)流体场模型建立完成后,首先要在流固耦合的边界上施加流体耦合标签FSI,然后在在流体场区域施加必要的边界条件,诸如位移约束,速度、压力等等。然后设置求解流体场的时间步长、求解时间、流体属性,打开ALE选项(瞬态分析)网格重画属性等等
(3)这样的工作完成后,进入/prep7,加入固体单元,设置固体材料属性,在挖去的部分建立固体模型,划分固体网格,在固体网格与流体场接触的固体边界上施加流体耦合标签FSI,注意要和前面的number相同。
(4)为固体实际必要的约束条件(看是固体推动流体还是流体推动固体)
(5)设置固体求解的时间步长和求解结束时间
(6)设置流固耦合属性,(看是固体推动流体还是流体推动固体),求解时间步长和求解时间,收敛准则,迭代次数等等。
(7)保存求解。
-----------------------------
总之,在流固耦合分析中,你最好要按着先流体后固体再耦合的属性设置顺序,流固耦合标签FSI要分别加在流固耦合边界的流体边界上和固体边界上,加在的顺序要按照上面所述。在实际的建模中,流固耦合的边界上由于建模的原因会出现节点的重合现象,注意一定不要使用捏合节点的命令来将重合的节点变成一个,这个很重要。固体单元一定要设置求解的时间步长和求解中止时间,时间步长一般和求解流体场和流固耦合的时间步长相等。
ALE+remesh选项是解决瞬态流固耦合问题的一个很重要的方法,流固耦合一般要伴随着流体的形状改变和位置的移动,因此首先启动ALE选项使流体与固体的耦合边界保持一致并规则化流场内部由于挤压而畸变的网格,其次若网格畸变的过于严重,就要启动remesh选项重新划分网格单元。
看一个流固耦合时间步长设置的是否合理的方法是在求解中若发现软件频频调动remesh进行网格重画,那这个时间步长基本上处于最大的时间步长上,若很久出现一次或不出现,说明你的时间步长太小,可以适当提高,但是正常的情况下,软件只允许在一个求解过程中最多有100次网格重画,若超过,将退出求解。
---------------------------------
先写到这些,匆匆草就,可能有很多缺点和错误,大家指正(注意若固体推动流体且又不考虑固体的形变和应力等等,可以直接按照帮助文档上ALE一章)
还有,若流体推动固体的分析中,固体发生了大位移,那么一定要在/solu中打开大变形分析选项,否则求解不收敛而退出

如果是FSI的话,只需要设置FSI的时间参数就可以了吧?!不用固体和流体分别指定了。另外,如果remesh次数太多,在计算一定时间后,每一个时间步都重新划分,你怎么解决?
先试着回答第一个问题吧。在流固耦合中,只设置FSI的时间参数可能是行不通的,FSI边界是动量能量的传递边界,它的参数设置只是定义动量能量以怎样的方式在流体和固体间传递,而流体和固体的相关时间参数设置是定义了:(1)当流体推动固体时,流体的流动参数计算>FSI耦合边界>固体的变形或移动参数计算(2)当固体推动流体时,固体的变形或移动>FSI耦合边界>流体的流动计算。因此本人认为FSI是流固耦合的一个中间环节,没有它,流体和固体的能量和动量的传递不能进行,同样若只有它,流体和固体的参数设置不知道,能量和动量也没有办法传递。所以我认为必须同时指定流体、固体和FSI的相关参数设置。
------------------------------
第二个问题:好像帮助文档上的ALE一章那个例子就是在每个时间步内都重新进行了remensh,本人认为ALE+remesh的原理是这样的,在任意一个时间步内,流体网格的变形都会通过ALE方法来进行规则化,(1)保持流体耦合的边界
(2)减小内部网格的畸变程度。然而当ALE规则化后的网格不能满足你所定义的remesh的最小限度时,就要启动remesh选项重新划分了。这样我认为若要降低remesh的次数你可以通过降低时间步长来实现,本人曾经计算过一个大位移的例子,是为流体输送固体颗粒的案例,当采用较小的时间步长时,整个计算结束都没有出现一次remesh。
---------------------------------
以上是我个人的一点看法,欢迎大家批评指正。

根据楼上的提出的一些问题我再作一次详细的解释吧。
----------------------------
1“挖去固体所占的空间”应该这样来解释。我们知道在流固耦合中,固体可能处于流体的空间内部,比如ANSYS帮助文档上那个橡皮垫正是处于流体场的内部。为什么要在流体场的建模时要挖去固体所占的区域呢,这是流体场网格采用ALE算法而决定的,ALE算法首先可以很好的跟踪流体或者是固体的运动边界(这点和固体分析时Lagrange算法是相似的),其次它可以通过自身算法的特点来规则化内部网格的畸变,而Eluer算法则不行,Euler算法在计算流体时必须将网格划分到流体可能流到的所有区域。因此采用ALE方法做流固耦合的流体场建模时是要挖去固体区域的,这样在计算过程中即使固体发生了变形后,采用ALE算法流体网格仍然可以拟合变形后的边界。挖去的固体空间这个过程只是在做流体场的建模时挖去的,我们要在固体区的建模时补上。
2.我所说是一种直接流固耦合的方法。
3.流固耦合的时间步长设置要分作三个部分:流体的、固体的和流固耦合的。流体的就是在Sol-Flotranset up-Execution ctrl里面进行,Sol-Load step opt-Time/Frequenc用来设置固体时间步长,FSI里面设置流固耦合的时间步长。
4.耦合标签是成对出现的,一对流固耦合有两个number,比如1-1,2-2,这两个number是相同的,这样软件才能识别是一个耦合对。
5.不错就是在Solution-Define Load-Field Surface Intr里面设置流固耦合FSI

楼上的看样对ALE和remesh的理解有点和我不怎么相同
打个比方,ALE和remesh对网格的作用如同“改革”和“革命”这两个概念一样,ALE只是在有限的范围内适当调整单元内节点的位置,比如帮助文档上ALE一章就有一个高速物体撞击后的网格变化的比较一部分,它在调整网格的畸变程度时并不会改变网格单元的数目。而remesh则不一样,由于remesh要重新划分网格,因此所占用的计算机资源是很大的,在动态分析中联立这两种网格划分方式可以节省计算机资源。其次在流固耦合里面采用ALE方法有助于捕捉流体的运动边界。这点我想很重要。
一个完整的流固耦合分析包含这个三个过程(流体推动固体)
1.单一的流体场分析
2.通过耦合边界FSI传递动量能量
3.单一的固体(结构)分析
固体分析时只是继承了流场分析时通过FSI传来的动量能量,所有对单一的流体场和固体的计算设置只是普通的流场分析和结构分析设置,并不存在所谓的固体耦合参数。耦合参数设置是对于固体和流体同时而言的。不知楼上的如何理解的。

流固耦合力学一门新型的交叉学科。它研究流体与固体相互作用下的一系列力学问题。流固耦合力学渗透在国民经济的各行各业中。对流固耦合力学问题的探究、认识,为各类工程结构的安全可靠运行提供科学依据,是工程设计中的当务之急。
ADINA可以作到即实双向流固耦合(Real-Time Two-Way Fluid-Structure Interaction (FSI) )功能, 用户可以在流固区间各使用不同的网格, 且在界面用不吻合网格(Dis-similar Mesh), 但程序依然可快速处理流固双向耦合.界面定义, 数据传输及数学模式的耦合, 在ADINA中均有稳定且自动的处理方式.
除此之外, 在流固互动过程, ADINA可以自动重建分布网格, 以避免因过度变型而造成网格扭曲, 导至计算停摆发散.需注意的是, ADINA中的FSI是真实流固耦合, 而非一般CFD软件的单纯moving boundary. 而其网格重建也非一般的ALE方法, 或是mesh morphing, 而是真正的实时重建.

评分

1

查看全部评分

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-29 18:34 , Processed in 0.057920 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表