关于SPH模拟冲击爆炸问题的一些疑问及讨论

2008hmily

目前国内利用SPH来模拟冲击爆炸问题的算例分析不算很多，可借鉴的参考资料实在太少，有些问题拿出来希望各位高手能解答
1、关于SPH模拟在侵彻或爆炸中坑的边界如何精确描述（如图1）
                       
                   图1
如何来精确捕捉两种材料之间的界面？

2、SPH法在边界处的粒子存在缺陷，导致SPH法不能完全适用于整个区域，边界处的粒子只受到边界内粒子的作用而边界外的粒子对这些粒子没有影响。
     对于SPH粒子边界的处理，可参考的资料很少，采用“虚粒子”来表示对称面的对称约束在李裕春的书中有详细描述，还有FEM和SPH耦合处的固联约束也有描述，但对于SPH边界施加固壁边界约束（采用虚粒子）或无反射边界都找不到相关的参考资料，请问如何施加固壁边界约束（采用虚粒子）或无反射边界，用哪些关键字？
3、SPH算法往往需要消耗大量时间，那么粒子的排布是否可以像网格单元那样渐变分布呢，即在变形大的区域分布较多的粒子，在远端分布稀疏？
4、计算中*CONTROL_SPH中的光滑长度h和memory对计算有很大影响，会出现如图2的警告，但是有些情况下还能正常计算并最终得到正常解答，请问如何解决该问题？
参考帖子：http://forum.simwe.com/thread-996633-1-2.html
http://forum.simwe.com/thread-996608-1-2.html

                              图2
5、由于SPH只有粒子，并没有单元，材料物理属性都固定在粒子上，那若要得到应力，应变等单元特有的数据该如何提取？

6、对于介质中爆炸问题，如图3，炸药所处位置的介质是否需要挖去？“炸药爆轰及水下爆炸的SPH数值模拟研究”这篇文章中模拟时貌似并没有将该部分挖去，而是采用炸药粒子直接叠合介质，这样是否不合逻辑呢，通常采用网格单元时都需要将炸药所处位置的介质单元挖去。
出现新的疑问：若需要挖去该部分，即SPH粒子PART在内部称为不完全连续体，在lspp创建SPH粒子时，如何将该部分创建为一个PART？


                            图3



望高手解答，万分感谢!

2008hmily

附上图1中简单的例子的k文件，只是为了描述问题而做的例子，并没有考虑完全具体的情况，比如边界约束等，望见谅！
简单描述：炸药在混凝土板表面爆炸
k文件计算中会出现一行警告，即上述提到的问题4.

希望共同探讨

csuzy

SPH算法比较复杂，文件中可能出现错误。

loaferhit

这个问题真愁人啊，memory调了以后好像解决不了问题，我个人感觉好像不是哪里的问题，不知道咋整，有的时候能继续计算，有的时候就出现：按任意键继续。计算停止回复 2# 2008hmily

2008hmily

回复 4# loaferhit
还有，我单独设置炸药中心起爆，就不会出现上述的警告，但一旦和其他介质同时计算，则会出现问题，接下去我试试加密炸药sph颗粒试试

lovehf

使用AUTODYN

2008hmily

回复 6# lovehf
确实看到几篇文献使用autodyn的，但现在上手的是dyna，其实我也不做sph这方面的，只是最近看到顺便学习学习

lxg0929

关注本贴，学习中

liangsf204

其他问题，无法解答，粒子渐变功能可以使用，只是在网格划分时候，短长边比不要太小即可

fritter123

    刚刚学习SPH，特别想和对这个方法很熟悉的人共同探讨。希望从最基础的理论开始，然后到代码，而不是直接在现有的软件来做。您的几个问题我尝试回答下：
1、不同介质直接在一开始就设置成不同的粒子，就不存在边界识别问题。
2、SPH方法的边缘误差来源于核函数被边缘所截断，现在边缘处理有一些方法，如排斥法、镜像法、周期性边界等，但在软件中的操作未知。
3、粒子的非均匀分布正是我所期望的，但理论上，标准SPH模型非均匀分布必然带来误差，误差来自离散化过程。可能需要修正的SPH模型。
4、光滑长度增加，会导致内部精度提高，而边缘精度降低。
5、应力应力单元都是有的，算法必须，这是运动的基础。软件中如果提取未知。
6、不知道。
希望找到对SPH模型共同探讨，可以通过合作的形式开展工作。我的QQ:273176055，请注明SPH