A free program to generate interface elements in an existing FE mesh
Last year when I started implementing interface elements to model material failure, I realized that the formulation is easy except how to generate a mesh with interface elements. I did a googling to search for such a free program. Amazingly, I did not find any although there are many researchers working on the fracture mechanics field.
So, I wrote a small object-oriented C++ program which reads a FE mesh, duplicates nodes and insert interface elements where asked. The program is able to
(1) insert 1D/2D interface elements everywhere in a FE mesh
(2) insert i1D/2D nterface elements along material interfaces, useful for delamination analysis or interface debonding
(3) insert 1D/2D interface elements along grain boundaries for a pollycrytalline solid, useful for inter-granular crack modelling in pollycrystals.
The program is written in such a way that any FE mesh format is supported. Currently only Gmsh format is however implemented. Some examples are given below.
If you are interested in this program, feel free to email me, email@example.com. I will find a place to put it on the web.