最近在进行算法的验证工作，发现案例中的mls-mpm边界处的粒子会出现应力分布不均匀的现象，如下图所示，mpm128案例中果冻弹性体的竖向应力分布云图以及应力路径图：
这应该是边界网格内的质量点在计算仿射矩阵（或速度梯度）时只能用到6个网格的节点速度，导致最靠近边界处的粒子放射矩阵过大，而第二层粒子的放射矩阵过小。请问各位大佬有没有什么解决办法呢？
