关于SPH粒子飞散的问题

ericmuse

各位大神好，我最近刚开始学习ls dyna的SPH模拟，参考了很多论坛上的算例和教程。
我是想做一个刚体在JHC模型上的划痕/犁沟的模拟，有点类似切削的过程，单位是g-mm-ms，底面和两侧SPC约束(试过rigidwal)，材料参数是按照文献数据换算单位以后设定的。

但是试了好多遍，每次都是刚体运动2-3步之后SPH模型就飞散掉，完全没法算。搜了搜论坛里之前有人遇到过类似的情况，但是参考他的方法，不管是加大SPH粒子的MAXV还是检查sph的mass，都不能解决。


K文件发上来，请各位前辈大神指点。如果还有其他错误特别欢迎一并指出。
==========更新问题=====
我后来尝试吧sph的mass增大了e7倍，飞散问题貌似解决了，但是刚体运动会受到影响，接触到sph的时候会突然被弹飞，然后下落穿透模型。而且粒子质量增加这么多倍是否会影响计算结果？


补张飞散的图




另外还想问一个关于后处理的问题，怎么让sph粒子显示体积增大，结果动画的渲染效果更好一点呢？

蝶恋花

看一看单位是否统一

ericmuse

蝶恋花 发表于 2016-3-21 16:35
看一看单位是否统一

恩，我应该很仔细检查过了的，主要检查了材料参数那一块，我再重新都查一遍吧。除了单位，还会有什么原因可能造成粒子飞散？

除了粒子飞散，我给刚体小球加载，不管是用load rigid还是load segment，只要加载了力，小球运动第二步就会立马起飞。。力越大飞的越快（力方向肯定没错）
这种情况又是因为啥呢？

十分感谢帮助

蝶恋花

粒子的间距

ericmuse

蝶恋花 发表于 2016-3-22 15:02
粒子的间距

谢谢解答，我后来把单位换成 g-cm-μs 之后解决了飞散的问题。虽然还是不明白换个单位制为什么有这个效果。。

我生成SPH的方法是ansys划好网格导出K文件之后，在prepost里直接用solid node生成SPH粒子，所以这个粒子的间距应该怎么调整呢？

zqpcgm

ericmuse 发表于 2016-3-22 18:54
谢谢解答，我后来把单位换成 g-cm-μs 之后解决了飞散的问题。虽然还是不明白换个单位制为什么有这个效果 ...

粒子间距就按网格尺寸来的，网格密，粒子密。

破ui

ericmuse 发表于 2016-3-22 18:54
谢谢解答，我后来把单位换成 g-cm-μs 之后解决了飞散的问题。虽然还是不明白换个单位制为什么有这个效果 ...

楼主  你现在知道  为什么改掉单位后  就能解决粒子飞散问题了吗