# Have someone met a strange behavior of Load-Displacement curve in a simple 2D indentation problem

I am modelling a very simple 2D contact problem between a rigid wedge indenter and a deformable squared-shape specimen (general steel material) in frictionless contact case. I used an implicit function f to describe the rigid indenter.

I implemented Lagrange multiplier method for this contact problem and the contact condition is: Inside the contact zone of the deformable specimen, the node set is active if f(node)<=0 and the Lagrange multiplier lambda<=0.

I also generated an exactly the same model on ABAQUS in order to make a code2code comparison of the ouputs. And these are the things I see:

a. Displacement: I must say it is 99% similar considering displacement distribution and value between my implementation and ABAQUS results. (There is only a little bit tolerance at the displacement value).

b. Load-Displacement curve: the ABAQUS curve is very reasonable with a linearly increasing. However, it is so strange in my implementation due to the fact that the load is equal to zero in 3 first steps of an increasing displacement u.

I attached my model, and the 2 load-displacement curve outputs: one is my implementation and the other is ABAQUS result, so that you can see how strange it is.

Thank you very much,
Best regards,

Minh

### This problem i see

This problem i see before in my work.

This Zero force at start may be due to check the contact at the start.

In almost contact problem you have to define a small space between the indenter and the specimen to check contact and in a second step you can apply your load.

So, you can define very small distance between the indenter and the specimen say 0.005 mm and apply displacemnt equal to this value. in aseoned step of loading you can define your load.

Regards

Ahmed

### The problem is still there. :-((

Thank you very much for your feedback. I did try to follow your idea. However, the problem is still there.

Best, Minh 