Using the concept of stress intensity factors, we suggest a way to include adhesion into boundary elements simulation of contacts. A local criterion concerning the maximum admissible surface stresses decides whether the adhesive bonds in particular grid points fail or not. By taking into account the grid spacing, a robust methodology is found. Validation is done using the theoretically derived cases of JKR adhesion.