关于RPIM形状功能的问题

我正在用无网格方法求解一个简单的固体力学问题。使用RPIM计算形状函数。

当我计算形状函数的导数时，我发现当点在问题域的边界上时，在计算点的导数不接近于零。如果计算点在问题域中，则可以得到较好的求导结果。

有没有人有使用RPIM形状函数的经验来帮助我解决这个问题?