三维sixDoFRigidBodyMotion算例发散



  • 尝试使用sixDoFRigidBodyMotion计算一个弹簧连接的三维圆柱在来流下的振动情况,但是结果发散。

    在算到发散的前几步之前,库朗数开始急剧增加。

    Courant Number mean: 0.097781 max: 12712.2
    

    然后后处理发现本来应该随着刚体一起运动的近壁网格发生了变形,图如下(红色是库朗数大的地方,同样是网格变形的地方)
    77e82604-f266-4283-81ff-19f3c442bf43-image.png

    二维圆柱动的算例之前已经跑通过了,所以我先以为是constraints约束不够,后来加了几个约束后还是发散。 然后考虑到是不是icem没有关联好,cylinder边界重新关联了一下还是没有解决。

    现在不知道问题出在哪里:有可能是感觉是innerDistance 0.1没有起到作用,但是使用slice切片观察Co发现圆柱中部的库朗数没有边界处这么大。所以不知道是不是约束没有设置合理。:136:
    5cbf4059-88ee-4c7d-af7d-93afc6d81b3f-image.png
    [0_1607610864136_pointDisplacement](正在上传 100%)
    附件中是相关的文件,不知道有没有人能帮忙解答一下。。网格太大了传不上来 ,,,
    dynamicDictandpointdisplacement.zip


  • 版主 副教授

    网格变形成这样,肯定算不对啊。
    所以找问题得找:为什么网格变形成这样。
    如何不让网格变形成这样。

    至于发散、库朗数之类的,得在网格正常的基础上讨论。

    看起来是圆柱的移动逼近了极限。所以要么限制圆柱的运动(离深蓝色和浅蓝色交界处原点),要么扩大运动的范围(就是可以变形的深蓝色区域)。:chigua2:



  • @bestucan

    最后找到了解决办法,是边界条件的原因。本来上下圆柱连接的上下两个边界用的symmetry边界条件,改成cylic边界之后就能算了。。。具体原因还没研究出来:136:


Log in to reply
 


CFD中文网 | 东岳流体学术 | 东岳流体商业 | 吉ICP备20003622号-1