【原创】SPH在入水冲击、机轮滑水和水上迫降中的应用
本帖最后由 cnetuwmj 于 2012-1-18 22:58 编辑本人毕设方向为SPH的流固耦合工程应用
以下就是几个例子,供大家参考
首先是SPH法在飞机水上迫降中的应用
由于项目secrect,飞机方面K文件不予分享,表示抱歉~
但是下面几楼中有楔形体和球入水的K文件,大家可以下载作为参考。
里面有些设置并非完美的,希望同学们继续探索,能比我走的更远一些。
jijingzhisheng 发表于 2011-12-25 20:59 static/image/common/back.gif
楼主很强大。。我也在做sph的,半年了,仍然有不少问题。想请教下楼主关于SPH边界条件的问题。SPH施加边界 ...
过奖了,多个SPH的PART之间不需要接触,我看到过类似的模拟,但是自己没做过。边界问题,我一般对称问题就是对称平面用SPH SYMMETRY,底部刚性墙,四周SPC或者刚性墙均可,但是全部SPC貌似会出问题。如果外围有FE单元耦合,那非反射边界条件自然是上乘之选啦~ 此乃神贴也~~~ 楼主很强大,学习了 本帖最后由 cnetuwmj 于 2011-12-22 22:00 编辑
在球冲击中的应用
前后处理LSPP3.1
这里粒子是加了效果的~见后处理~不过显示的漂亮意味着耗内存
K文件上,很大,压缩了之后还很多,大家慢慢下~
本帖最后由 cnetuwmj 于 2011-12-22 22:01 编辑
K文件来啦~
楔形体入水冲击,ALE和SPH我都做了尝试,其中有一组是不同疏密的SPH粒子间距和ALE网格大小的比较
K文件已传,其中SPH粒子很稀疏,大家当做试练吧~
本帖最后由 cnetuwmj 于 2011-12-22 21:37 编辑
这是机轮溅水的后视图~
这是当年做的一个SPH+FE的唯一例子,很可惜K文件现在都找不到了。
在我看来虽然SPH+FE可以减小计算量,而且FE能加非反射边界,但是计算结果的精度还不如纯SPH模拟单相~在水体内部的SPH粒子和网格的边界上,接触设置再完美也是两种离散方式之间物理量的传递~从目前MIT等大学的博士论文来看,结果不如全SPH的结果~
LS-DYNA中SPH的边界条件只有对称边界是特别针对该法而提的,事实上还应该有更多的边界条件
包括状态方程,MONAGHAN 的EOS为什么不加进去是我很想不明白的。
至此,算例基本分享完毕。
正值毕业季,我也很忙,平时很多朋友的留言,站内信我都没有即时看到~
玩SPH的人很多,我相信多算多想多看才是王道~
妞妞~2011.12.22 @学院路37号NMB
好,不错!顶!这对我们今后使用仿真软件有关键性的决定性的积极的战略性的指导作用! luther6789 发表于 2011-12-22 21:52 static/image/common/back.gif
好,不错!顶!这对我们今后使用仿真软件有关键性的决定性的积极的战略性的指导作用! ...
额,评价太高了,有点受不住了 咦,学院路,难道是北航的童鞋 此贴必火!!!:)
支持楼主,坛子中SPH方面的最高手了!:)
如果时间方便,不知楼主能否以锲形块入水为例做个step by step教程,以及推荐一些SPH方面比较有帮助的参考文献和书籍? 20020824 发表于 2011-12-23 10:53 static/image/common/back.gif
此贴必火!!!
支持楼主,坛子中SPH方面的最高手了!
如果时间方便,不知楼主能否以锲形块入水为例做 ...
下周吧,我下周正准备做教程呢 20020824 发表于 2011-12-23 10:53 static/image/common/back.gif
此贴必火!!!
支持楼主,坛子中SPH方面的最高手了!
如果时间方便,不知楼主能否以锲形块入水为例做 ...
另:高手真算不上啦~很多人都潜伏着呢。。。 happydragon 发表于 2011-12-22 23:15 static/image/common/back.gif
咦,学院路,难道是北航的童鞋
你懂的~ 期待楼主的教程,:D楼主我做的SPH可以运算,就是粒子不见动,只是闪闪而已,不知道是哪里的问题,楼主可以留下联系方式吗?谢谢了:D 远方的小牛 发表于 2011-12-23 14:57 static/image/common/back.gif
期待楼主的教程,楼主我做的SPH可以运算,就是粒子不见动,只是闪闪而已,不知道是哪里的问题,楼主可以 ...
问题描述不清楚,发图或者K文件~有人路过也会你些意见的 高手,我最近要做一个SPH切削模型,希望楼主能指点下 哈哈,还可以,学习了~ 楼主强悍啊! 本帖最后由 jijingzhisheng 于 2011-12-25 21:48 编辑
楼主很强大。。我也在做sph的,半年了,仍然有不少问题。想请教下楼主关于SPH边界条件的问题。SPH施加边界约束时用什么合适呢?*CONSTRAINED_GLOBAL,SPH_SYMMETRY_PLANE还是spc? 后面的两种我都用过,貌似计算结果没多大差别(有人说sph不适合用spc)。对于CONSTRAINED_GLOBAL没弄明白是怎么回事,没用过。。。
还有个问题,我在一篇论文里看到过作者把两个三维的part都用sph划分,让两者接触碰撞,这可能吗?两个三维的SPH好像没办法接触吧?二维的倒是见过,也不需要设置接触。