I have a question about newton solver convergence problem. I am simulation a domain with two different material properties. The first part is soft material, and Young’s modulus E_s is 3MPa, the left is hard area, where E_h is 800MPa, when I run it in dolfinx, the algorithm doesn’t converge.
But when I set the hard material Young’s modulus E_h = 80MPa(1 order low than the original parameter), and soft area E_s = 3MPa (doesn’t change). Then the model converges.
Could someone tell me why and how I can fix the problem?