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