batmanabcdefg / c1aba

Tetrahedral 4-node C1-continuous finite element for Abaqus

Incompatible mode C₁-continuous tetrahedral 4-node 3D user element for Abaqus

This repository contains all the necessary tools to use a tetrahedral incompatible mode C₁-continuous user element in Abaqus that has the following features:

1. Fully compatible with C3D4-meshes
2. Large strain formulation
3. Simple user material interface
4. 4 degrees of freedom (DOF, field value and gradient) at each corner node, i.e.~16 total DOF per scalar field
5. Full 5^th order polynomial shape function
6. Numerical integration exact up to a full 6^th order polynomial

Walkthrough

1. Element formulation
2. Mathematica: Shape function derivation and code generation
3. Fortran UEL
4. Abaqus example
5. Visualisation with Paraview