For impacts using a penalty approach, contact laws can be implemented in user functions as described in the tutorial in “External forces” part.
For impacts using a non-smooth approach, a coupling with lmgc90 has been developed. More information here (lien web)