All,
I am trying to leverage the work from @bleyerj from Cohesive zone modelling of debonding and bulk fracture to approach traction separation laws using dolfinx. I have worked to replicate the original code in the updated framework but am facing numerous challenges in its implementation, primarily in selecting appropriate function spaces and solver settings.
Has anyone done similar work in dolfinx?
Hi, sorry I did not have time yet to port this demo to dolfinx. I will try to do it in the near future though…
Can you be more specific on what is your exact issue ?
@bleyerj , Thanks for the reply. My intent was just to ping the community, happy to provide specific details if you are interested.
Selection of Function Spaces
I have done little work with quadrature elements and spaces other than CG/DG. As per your comment in the above example
The ideal solution would be here to use Quadrature elements on the mesh facets, which are unfortunately not available in FEniCS. We will instead use a
Discontinuous Lagrange Trace(also known asHDiv Trace) space which represents polynomials living on the facet of the meshes and which are discontinuous at the vertices.
The issue hence is deciding the element family for V_f. I see from hdiv-trace-elements that DGT is not supported yet in Dolfinx. I am not familiar with using Quadrature elements in Dolfinx though I am trying to learn. Any advice is appreciated.
Thanks,
-Sven