“…Herein, is the identity, and represents a suitable discrete approximation of the mapping χ between the contact sides, see, eg, Dickopf and Krause and Puso for more details. Discretization of the nonpenetration constraint in yields the discrete gap function at each node j , ie, …”