@FEA way, thank you, that's a great point.
I am trying replicating another example, with no contact, and prescribred displacement. In doing that, I have a side question: when exporting the stiffness matrix through the input lines in the original post, I always get some nodes with negative...