HFEM

HFEM is a simple teaching code to illustrate the basic features of the finite element method. HFEM is written in Fortran 77. It is also used for an easy tool in student's programming exercises. The present version solves only stationary problems.

hfem.tar.gz, version MAY-2002 gzipped tar-archive. You can unpack the archive by executing gtar -xvzf hfem.tar.gz .

The archive contains:

• Source files (suffix .f)
• Makefile hf-ownblas.mk (or hfem.mk)
• four input files (suffix .inp) and the corresponding output files (suffices .out , .mes , .res, .neu ).

Program characteristics

• Heat transfer, also with convective and reactive terms (1,2 and 3 D)
• Stress analysis of 2 (plane stress and strain) and 3 D solids
• 2D Euler-Bernoulli beam (frame)
• Truss element (1, 2 and 3 D)
• 2D isoparametric arch element (2 to 4 nodes)
• 3 and 4 node stabilized MITC plate elements
• Standard isoparametric (bilinear, biquadratic) RM-plate elements

• 1D: linear, quadratic and cubic Lagrange interpolation and the first order Hermitian polynomials.
  Hierarchical Lagrange (C0) and Hermite (C1) interpolation of arbitrary degree.
• 2D: linear, quadratic for triangular elements and bilinear and biquadratic (also 8-node Serendipity) for quadrilateral elements.
• 3D: trilinear and reduced triquadratic (20-node Serendipity) for hexahedral elements.

• Direct band Cholesky for symmetric positive definite systems and band LU solver for unsymmetric problems. Block band Cholesky for SPD problems.
• Conjugate gradient and bi-conjugate gradient algorithms for iterative solution of the linear system. Incomplete Cholesky (IC) and incomplete LU (ILU) preconditioners included. Storage modes: Compressed Sparse Row (CSR) and sparse DIAgonal (DIA) (only IC coded with DIA mode).
• Gibbs, Poole, Stockmeyer (GPS) bandwidth reduction algorithm for band solvers or the Gibbs, King profile reduction algorithm for skyline or envelope solvers.