Please find attached the C++ code needed to add contact limit forces to the model.