To download these files, rightclick the
link and Save As... or down load a
zipped file of all the codes

FEM_1D_Static.m Simple
static analysis of a 1D bar subjected to axial body force

FEM_1D_newmark.m Simple
dynamic analysis of a 1D bar subjected to axial body force, using
Newmark time integration.

FEM_1D_modal.m Simple dynamic analysis of a 1D bar subjected to
axial body force, using modal time integration.

FEM_conststrain.m Simple 2D FEA code with constant strain
triangles. Should be run with the input fileFEM_conststrain.txt or
(for a larger problem) FEM_conststrain_holeplate.txt

The
following files all solve 2D or 3D static linear elastic problems, but
illustrate various refinements of the finite element method:

FEM_2Dor3D_linelast_standard.m 2D(plane
strain/stress) or 3D static linear elasticity code with fully
integrated elements. The code can be run with the following input
files.

Linear_elastic_triangles.txt: 2D plane strain problem
with two triangular elements;

Linear_elastic_quad4.txt: 2D
plane strain problem with 4 noded quadrilateral elements;

Linear_elastic_quad8.txt: 2D plane strain problem with 8
noded quadrilateral elements;

Linear_elastic_brick8.txt: 3D problem with 8 noded brick
elements;

Linear_elastic_brick20.txt: 3D problem with 20 noded
brick elements;

Linear_elastic_pressurized_cylinder.txt:
2D simulation of a pressurized
cylinder

FEM_shear_locking_demo.m 
Solves the beam bending problem discussed in Section 8.6.2, and
compares the FEM solution with the exact solution to illustrate
shear locking. This version of the code must be run withshear_locking_demo_linear.txt (solution
with 4 noded quad elements).

FEM_incompatible_modes.m 
Solves the beam bending problem discussed in Section 8.6.2 using
incompatible mode elements, and compares the FEM solution with the
exact solution to demonstrate that the elements avoid shear locking.
This version of the code must be run with shear_locking_demo_linear.txt (solution
with 4 noded quad elements).

FEM_volumetric_locking_demo.m 
Solves the pressurized cylindrical cavity problem discussed in
Section 8.6.2, and compares the FEM solution with the exact
solution. This version of the code must be run withvolumetric_locking_demo_linear.txt (solution
with 4 noded quad elements) or volumetric_locking_demo_quadratic.txt(solution
with 8 noded quadrilateral elements).

FEM_hourglassing_demo.m 
Solves the pressurized cylindrical cavity problem discussed in
Section 8.6.2 with reduced integration elements, demonstrating
hourglassing. This version of the code must be run withvolumetric_locking_demo_linear.txt (solution
with 4 noded quad elements).

FEM_selective_reduced_integration.m 
Solves the pressurized cylindrical cavity problem discussed in
Section 8.6.2 using selectively reduced integration, and compares
the FEM solution with the exact solution. This version of the code
must be run with volumetric_locking_demo_quadratic.txt (solution
with 8 noded quadrilateral elements).

FEM_hourglasscontrol.m 
illustrates use of hourglass control to eliminate hourglassing in 4
noded quadrilateral elements. This version of the code must be run
with volumetric_locking_demo_linear.txt (solution
with 4 noded quad elements)

FEM_Bbar.m –
Solves the pressurized cylinder problem discussed in Section 8.6.2
using the Bbar method, and compares the solution with the exact
solution. This version of the code must be run withvolumetric_locking_demo_linear.txt or volumetric_locking_demo_quadratic.txt.

FEM_hybrid.m – Solves
the pressurized cylinder problem discussed in Section 8.6.2 using
hybrid elements, and compares the FEM solution with the exact
solution. This version of the code must be run withvolumetric_locking_demo_linear.txt or volumetric_locking_demo_quadratic.txt.

FEM_2Dor3D_linelast_dynamic.m: Solves 2D or 3D dynamic linear
elasticity problems, using Newmark time integration. The code can be run
with the input file Linear_elastic_dynamic_beam.txt.

FEM_2Dor3D_modeshapes.m: Calculates
mode shapes and natural frequencies for a linear elastic solid. The code
can be run with the input file Linear_elastic_dynamic_beam.txt.

FEM_2Dor3D_hypoelastic_static.m: Solves 2D (plane strain only) or 3D
static problems for a hypoelastic material, as discussed in Section
8.3.9. The input file is Hypoelastic_quad4.txt.

FEM_2Dor3D_hyperelastic_static.m: Solves 2D (plane strain only) or
3D static problems for a hyperelastic (NeoHookean) material. An input
file is provided in Hyperelastic_quad4.txt.

FEM_2Dor3D_viscoplastic_static.m: Solves 2D (plane strain only) or
3D static problems for a small strain viscoplastic material. An input
file is provided in Viscoplastic_quad4.txt
This web site was originally
developed by
Charles Camp for
CIVL
7117.
This site is
Maintained by the
Department of Civil Engineering
at the University of Memphis.
Your comments and questions are welcomed.
