Hi, I’m also quite interested in the answer to this question. From what I was able to gather from this forum, this should be achievable using one of the following technique:
- Using a DG formulation, although this can get tricky and not very practical
- Using Multiphenicsx such as in this example.
- Using Dolfinx_MPC, although it’s not clear to me how to relate the traction-separation law to one of the implemented constraints.
I still haven’t found or been able to have a simple working example though. If you manage to do it, let me know!