Home

Welcome to the FEM-Explicit2D wiki!

I am trying to separate the subroutines and combine them into modules to achieve data encapsulation. This also avoids the trouble of using common blocks. Right now I have connectivity module, element module, explicit module, assembly module and BC module.

The Connectivity module reads in the node, element, and BC data. The subroutines inside this module can determine the number of nodes, elements, and BCs. I used LS-Prepost to generate the mesh data so my code only supports that format, in which data are separated by space and there are 8 columns per line for BC data. It can also provide node coordinates of an element providing the element id.

The Element module calculates the Jacobian, shape function, the derivative of shape function, B matrix and mass of the element. I use 2 by 2 Gauss integration point and first-order quadrilateral element.

The explicit module calculates the critical time step, initializes the internal, external force vector, displacement, velocity, and acceleration datatype.