求助大佬给解决下计算不收敛问题！

最近在做基于数字图像细观建模相关方面的课题，利用1比1像素还原花岗岩真实细观结构，但计算一直出现不收敛问题，采用的显式算法求解也不行，应该是接触出现问题，但自己是个小白，不是特别懂，希望平台大佬给予一些参考意见，谢谢！图一是放大后的压头附近的应力云图，图二和图三是垫块附近的受力云图。为什么左边的垫块出现了严重的接触穿透现象？

                                                                    图一

                                                                       图二

                                                                    图三