找回密码
 注册
Simdroid-非首页
查看: 1075|回复: 17

[09.流固耦合] 气泡在液体中动态运动的流固耦合结果第一周期振荡

[复制链接]
发表于 2011-10-8 11:06:00 | 显示全部楼层 |阅读模式 来自 广东深圳
悬赏10仿真币已解决
本帖最后由 yuanzitou 于 2011-10-12 14:00 编辑

物理问题介绍: 小气泡在弹性软管中振荡(超声激励下),想看软管壁的受力问题。

模块:NS,Solid,ALE

用到了弱约束解决表面曲率问题,在这里对soliton同学和kinggun表示深深的感谢!

动画:


原以为问题圆满解决了。查看内壁的受力,出问题了:在开始的时间段,结果振荡得非常厉害,如下图:

再查看血管壁面的pressure,也出现振荡,如下图


有什么能解决的办法吗??

谢谢大家! 感激不尽!!卡在这里很久了。

附件: 您需要 登录 才可以下载或查看,没有账号?注册

最佳答案

查看完整内容

很好奇这个现象出现的原因到底是什么。 目前楼主采用的这个办法从结果的角度当然没有,但从数值计算的角度,到底是问题本身就是这样,还是数值计算产生的?
发表于 2011-10-8 11:06:01 | 显示全部楼层 来自 上海
Simdroid开发平台
很好奇这个现象出现的原因到底是什么。

目前楼主采用的这个办法从结果的角度当然没有,但从数值计算的角度,到底是问题本身就是这样,还是数值计算产生的?
回复

使用道具 举报

发表于 2011-10-9 10:56:11 | 显示全部楼层 来自 上海
没有太明白这个,问题到底在哪里?

直觉是对的啊!
回复

使用道具 举报

 楼主| 发表于 2011-10-9 14:25:54 | 显示全部楼层 来自 广东深圳
soliton 发表于 2011-10-9 10:56
没有太明白这个,问题到底在哪里?

直觉是对的啊!

就是 那个力的曲线 在一开始的时间段,不平滑,不像后面几个周期的结果。而我要的恰恰就是最开始的。。。
回复

使用道具 举报

发表于 2011-10-9 15:50:25 | 显示全部楼层 来自 上海
yuanzitou 发表于 2011-10-9 14:25
就是 那个力的曲线 在一开始的时间段,不平滑,不像后面几个周期的结果。而我要的恰恰就是最开始的。。。 ...

至少基本的趋势是对的了。看起来是数值计算过程中的一个振荡。

记得你提过,控制时间步长,计算的容差等等都试过了。不知道这个结果和你取结果的“点”有没有关系,看看周围的点是什么趋势?(或者细化你感兴趣区域的网格?)
回复

使用道具 举报

 楼主| 发表于 2011-10-9 17:03:44 | 显示全部楼层 来自 广东深圳
soliton 发表于 2011-10-9 15:50
至少基本的趋势是对的了。看起来是数值计算过程中的一个振荡。

记得你提过,控制时间步长,计算的容差等 ...

液体中的点,取pressure值,也是振荡的。看我第三张图片。

所以不知道什么原因
回复

使用道具 举报

 楼主| 发表于 2011-10-12 13:58:48 | 显示全部楼层 来自 广东深圳
问题解决了!!

在前面加个时间段,然后不施加激励的超声条件,就解决了,多谢kinggun!

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×

评分

1

查看全部评分

回复

使用道具 举报

 楼主| 发表于 2011-10-12 13:59:56 | 显示全部楼层 来自 广东深圳
悬赏结束
回复

使用道具 举报

 楼主| 发表于 2011-10-17 15:47:20 | 显示全部楼层 来自 广东深圳
我现在正在调参数。我觉得不是问题的本身引起的。因为后续几个周期就很平滑。这个应该是数值振荡。我发现这样只是权宜之计,还是得仔细调程序
回复

使用道具 举报

 楼主| 发表于 2011-10-20 17:36:53 | 显示全部楼层 来自 广东深圳
解决不了。。。请高手回答!~~

我在结构场用 准静态分析,发现不会有这个振荡。但是我这个问题貌似不适合准静态分析!!

回复

使用道具 举报

 楼主| 发表于 2011-10-20 18:33:13 | 显示全部楼层 来自 广东深圳
我发现分析类型为transient时,在equation system里查看subdomain setting,dweak这个选项卡第一栏是
rho_smaxi*(uor2_test*uor2tt*R^3+w2_test*w2tt*R)

这是代表什么意思?



        我把这一项去掉后,发现得出来的结果和准静态分析是一模一样的,而且之前的数值振荡现象也没有了。

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
回复

使用道具 举报

 楼主| 发表于 2011-10-20 18:44:08 | 显示全部楼层 来自 广东深圳
yuanzitou 发表于 2011-10-20 18:33
我发现分析类型为transient时,在equation system里查看subdomain setting,dweak这个选项卡第一栏是
rho ...

哦。我明白了。在transient的分析类型中,考虑了加速度的惯性力这个因素,而quasi-static的控制方程里面根本没有这一项。

还是对解决问题没有思路啊
回复

使用道具 举报

 楼主| 发表于 2011-10-27 09:31:05 | 显示全部楼层 来自 广东深圳
问题已解决。给结构场加上合适的阻尼即可!

评分

1

查看全部评分

回复

使用道具 举报

发表于 2011-12-25 20:30:47 | 显示全部楼层 来自 甘肃兰州
不太明白

评分

1

查看全部评分

回复

使用道具 举报

发表于 2012-5-16 10:07:22 | 显示全部楼层 来自 四川成都
学习!厉害....
回复

使用道具 举报

发表于 2012-5-23 10:23:41 | 显示全部楼层 来自 湖北武汉
应该是最开始阶段的数值震荡
回复

使用道具 举报

发表于 2012-5-23 16:01:09 | 显示全部楼层 来自 重庆沙坪坝区
问题:1 ,LZ用的是什么软件?
          2 ,结构里面的阻尼在哪里加?结构的阻尼好像是由材料产生的。
回复

使用道具 举报

发表于 2012-6-13 14:28:24 | 显示全部楼层 来自 日本
学习了~~~

评分

1

查看全部评分

回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-28 21:04 , Processed in 0.054147 second(s), 17 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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