There is some information in the Theory Manual. It is based on the Lagrange multipliers (LM) used to enforce the contact constraints. According to Analysis User's manual, Abaqus may do some further processing of LM (perhaps some averaging of the values provided by individual constraints) before converting them to CPRESS which are unique nodal values. However, I could not find much details about this.