I have compared contact analysis models using Abaqus, Lusas (using spring gap elements) , MSC/Nastran (also using CGAP elements) and FEMdesigner , on a pin within a spherical bearing which itself was within a spherical lug, i.e. two regions of contact. All solvers gave very close results, however FEmdesigner was more than 50% quicker than any of the other three !