You are here
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.
(4) support Discontinuous Galerkin methods.
The program is written in such a way that any FE mesh format is supported. Currently only Gmsh format is however implemented. The output can be written to Abaqus compatible *.inp files if needed. Some examples are given below and the algorithm is described in the attached paper which is the preprint of the article
VP Nguyen, An open source program to generate zero-thickness cohesive interface elements, Advances in Engineering Software, 2014(74): 27--39.
The code is available at http://sourceforge.net/projects/cigen/ .