ten1917 发表于 2009-3-2 13:36:30

关于流固耦合的三个问题(已附.in)

本帖最后由 ten1917 于 2009-5-16 15:47 编辑

最近建了一个模型:一个圆柱形液池中有一个梁(用ADINA-M建的)在底面的基底上振动。尺寸最小4微米(修改了Geometry下的Tolerance,才允许用按米单位建),用了自适应网格,这过程中发现以下三个问题:

1. 给梁一个位移载荷,然而打开Load Plot,却说没有载荷,后来是在生成网格时,可以显示,显示为0;然后就改magnitude与Time Function中的Value,发现反复改大这两个值,Load Plot显示出的值在e-23量级。

2. 出现上面的情况后,求解,在刚计算Time Number= 1 时,求解自动关闭,没有任何提示。之后,把Time Function的Value改为1E17,然后求解FSI,结果是求解就停在Time Step 1,算了一夜没动静!单元10000个多一点,显示的求解器才用了20多MB。
不过后来在另一个模型上,又在梁体上加了位移,却可以求解,区别仅在于这个模型载荷加的是集中力在一条Edge上(另外这个模型的单位不协调,我是新手)。

3. 我想让这个梁的一端振动,当幅值较大时,会与下面的基底接触,即振动接触,请问这种情况如何设接触。试过一次,好像接触对要求二者原先就已接触。

恳请高手能指点一二,非常感谢!

后来,把模型放大了100倍,可以运行到第三步,然后按照提示改小了Displacement reflexation Factor,可以运行到第5步,再改小,可以运行完全的20步。在可以运行的基础上,我把Displacement又改大(只是1E-7m量级),再次出现,在第一步不收敛的问题。

到底是我的位移载荷太大,还是模型太小,还是二者都有的问题?

附件是今天下午新建的模型,问题是在求解Time Step 1时停住,晚上才导出的.in文件,不知有没有什么影响。

mgcomr 发表于 2009-3-5 00:27:57

请问搂主,

如何改小 Displacement reflexation Factor? 请教教我可以吗? 我在想我自己的问题就是出现在这个上面。我想改变Displacement reflexation Factor. 但是不会。

谢谢你哦!

ten1917 发表于 2009-3-5 11:04:01

2# mgcomr

这个就在CFD下Program module条,当选中FSI时,右边会有一个蓝色的“FSI”(流体类型,如incompressible的左边),点击打开,最下面的两个就是了。看了手册,说这个是在运算时将什么矩阵除以这个因子,该值小于等于1。

对于上面的的三个问题,最近一直在试,怎么都不行。然后就在那个能够求解的model的基础上一点点改,发现:上面说得第一个问题没再出现;如果一个周期设置的时间步过多会不收敛;载荷过大也会不收敛。现在还在试错中……

期待有人详细介绍一下关于收敛的问题,手册上那些建议如果能够给出物理含义,以后这样的问题就会注意了

mgcomr 发表于 2009-3-5 12:47:52

3# ten1917 谢谢你阿!

jacobi 发表于 2009-3-5 22:10:08

尺寸最小4微米--怎么用m做建模单位??

太小的模型怎么不用小的单位系统做呢??
比如:mm,N,Mpa, t/mm^3
等更小的单位制?

感觉模型太小的话,数值计算精度损失殆尽,
比如32系统下,数值的有效位数才15位吧(记得不太清)。

其他的就不懂了,fsi是一窍不通.....

ten1917 发表于 2009-3-6 10:06:57

本帖最后由 ten1917 于 2009-3-6 10:14 编辑

1# ten1917

第一个问题我想已经找到原因了:是由于每次点击“Load Plot”时,显示的载荷加载在最后一个时间步,于是会显示一个很小的值;如果移动时间步再显示载荷,就没问题了。
很奇怪自己竟没注意到,抱歉。

ten1917 发表于 2009-3-6 10:12:59

本帖最后由 ten1917 于 2009-3-6 10:14 编辑

5# jacobi

因为开始时对协调单位制一窍不通(以前没用过有限元软件),后来我用过SI(mm),还是出错,不过也许那时载荷也有错误,而且我试的不多,正准备再试验。
而我前两天在放大的模型上试,用1Hz频率振动,时间步、步长设为20×0.05,是可以求解的。问题是一旦高频,比如角频率设为3600,就不行,或者时间步长小到0.001(1Hz),也提示不收敛。

很感谢。

ten1917 发表于 2009-3-8 09:06:36

4# mgcomr

你说得对,是单位的问题。周五周六我分别用mm、kg、s和um、g、s建了模型,在时间步长不太小的情况下可以求解。问题是一旦时间步长小于0.00001就不行,错误类型主要是ADF3045,而且在time Step 1,看手册,应该还是我哪里有错误,继续摸索
再次感谢!

前面的问题误导了,抱歉!

shegaoxiang123 发表于 2009-4-24 14:54:24

高频的用有限元没法搞,一般超过500HZ的就不用有限元来解决了!

ten1917 发表于 2009-4-24 20:56:31

9# shegaoxiang123

还是能做的,用ms的单位,我也用s的单位做过100kHz的。

另外上面的问题基本都解决了,主要注意的是在CFD中,FSI的计算选用Direct,不要用Iterative,这是前几天看到ADINA中国的一个网页“ADINA FSI Capabilities”中的一句话,然后试了,不再有如ADF3034、ADF3045的错误。

waynekey 发表于 2009-5-14 19:16:09

建议给楼主加分,解决了问题还把解决方法弄出来,大家才有提高的机会。

Crespo97 发表于 2009-9-9 14:10:19

非常感谢楼主ten1917 ,我的模拟跟你的有些相似,也跟你遇到了相似的困难。但是看了你的建议后,问题也解决了,非常感谢!
页: [1]
查看完整版本: 关于流固耦合的三个问题(已附.in)