CE 521: SOLID MECHANICS
FORTRAN PROGRAMMING
(Last Modified 3-31-99)
Instructor: Chris L. Mullen
Office: 202 Carrier Hall / Phone: 232-5370
Email: cvchris@olemiss.edu
Select Fortran program by "file.f" name given below:
ve-umat.f
You may access Fortran from from sweetgum.mcsr.olemiss.edu through a telnet
session or from the class accounts setup on zinger.ce.olemiss.edu and temblor.ce.olemiss.edu.
The file, ve-umat.f, contains the ABAQUS UMAT subroutine written by Prof. Hackett which is
based on the Maxwell-type 3D viscoelasticity model he presented in class.
Download and try to write a simple flow chart in order to help you relate
it to the class development. Also, read the appropriate chapters on User-Material
subroutines in the ABAQUS User Manual found in Room 206 before attempting
to implement it. You will need to see Dr. Mullen if you do not yet have
access to ABAQUS through temblor or magnolia.
An ABAQUS input file for testing the UMAT is given in umat.inp with excerpted output in umat.dat. This file computes the purely elastic response corresponding to the idealized
rod case with properties given in class. Verify that it gives the correct
top displacement(PL/EA) at full load (P=2.5 k) that are given by ABAQUS
elastic material models using 1D truss elements ( trus.inp,trus.dat ) and 3D hexahedral elements ( hex.inp,hex.dat).
After you have run the above files and verified accuracy, change the properties
in umat.inp to Dr. Hackett's values and run your own solution. Submit evidence
of all results and prepare a little write-up for credit on Hwk. #2