Can someone explain the strange physical behavior in contact problem using Lagrange method in modeling a simple 2D indentation?
I am modelling a very simple 2D contact problem between a rigid indenter and a deformable squared-shape specimen. I used a implicit function f to describe the rigid indenter.
The contact condition is: Inside the contact zone of the deformable specimen, a node n is outside of rigid indenter for f(n)>0, and inside for f(n)<0. In case f(n)=0, the contacting node lies on the surface of the rigid body.