我实现了 Node-To-Surface 接触算法（Wriggers, Peter , Computational contactmechanis. , Berlin: Springer (ISBN 3-540-32608-1/hbk). xii, 518 p. (2006). ZBL1104.74002 .）。代码是通过首先形成接触元件来完成的，接触元件包括球上的从节点和底座的主线。在球穿透软块的示例中，模拟可以运行到接触单元内的投影落在牛顿-拉夫森迭代中相应主单元之外的点。随后，隐式 NR 循环不再收敛。有没有人有处理这个问题的经验？

这是失败前的系统快照

这就是问题所在，注意到黑色圆圈，从节点（蓝色）靠近主节点之一（紫色），导致上述问题。