[ [K_T^i]\Delta u^i = F_ext - F_int^i ]
Newton-Raphson diverged when stiffness became negative. Displacement DOF at crown exceeded 1e7 mm. ansys an internal solution magnitude limit was exceeded
[ [K(u)]u = F(u) ]