“…Following Puso and Laursen [7,Section 2], the potential function of the normal contact constraints pertaining to the mortar method can be derived from an integral form of the contact complementarity condition. Accordingly, (1) , t)·(u (1),h (X (1) , t)−u (2),h (X (2) , t)) d =:…”