Pair Nodes in contact automatically. (Python script or in Abaqus CAE-ODB GUI)

Hi all,

I am predicting wear depth in Abaqus using Python scripts.

Initially, I need to find all the nodes which are in contact , and extract the label of those nodes. Afterwards, for a node on a part, find the closest node on the other part (which both are in contact). For instance if we have node No. 100 on part 1, in contact, I want to pair that with the closest node on part 2 which is in contact as well, say node No.200. The out put requested should be (100, 200). (I just need the node labels as an output)

I have done this procedure manually for an axisymmetric model. It takes ages to pair all the nodes manually for a 3D model with a refined mesh. I just wounder, whether there is any way to find the nodes in contact and pair them automatically to show all pair nodes label separately, either by using Abaqus CAE/ODB or Python script?

Any suggestions would be greatly appreciated.

Ariyan