Dear professor,
Thank you very much for your help! But I have another question: according to your answer, the codes of cnts3d.f when the csw =14 are incomplete, is it the reason that the 3D explicit contact cannot be executed correctly? I have made two tests, one is a 2D explicit contact test and it can be run successfully, because the program of cnts2d.f is complete including csw=14, the other is a 3D explicit contact test and it always reports errors, of course you have told me that feap cannot do the explicit contact, but cnts3d.f has the program when csw=6(an indicate of explicit solution), is there something that I have misunderstand?
Could you please give me more details? Because I have read the program in cnts3d.f(csw=6) carefully, but I cannot find a bug. According to my feeling, the 3D explicit contact should be run successfully, however, in fact, it can't.
Thanks for your help in advance.
Best.