
This information is from here: http://www.oecdnea.org/tools/abstract/detail/nesc9909/
(I ran my PhD models using this software on a Cray XMP at the University of London Computer Centre )
NESC9909 DYNA3D.
DYNA3D, 3D Finite Elements for Dynamic Response of Inelastic Solids
3. DESCRIPTION OF PROGRAM OR FUNCTION
DYNA3D is an explicit, three dimensional, finite element program for analyzing the large deformation dynamic response of inelastic solids and structures. DYNA3D contain 30 material models and 10 equations of state (EOS) to cover a wide range of material behavior. The material models implemented are: elastic, orthotropic elastic, kinematic/isotropic plasticity, thermoelastoplastic, soil and crushable foam, linear viscoelastic, BlatzKo rubber, high explosive burn, hydrodynamic without deviatoric stresses, elastoplastic hydrodynamic, temperature dependent elastoplastic, isotropic elastoplastic, isotropic elastoplastic with failure, soil and crushable foam with failure, Johnson/Cook plasticity model, pseudo TENSOR geological model, elastoplastic with fracture, power law isotropic plasticity, strainÂ Â rate dependent plasticity, rigid, thermal orthotropic, composite damage model, thermal orthotropic with 12 curves, piecewise linear isotropic plasticity, and inviscid two invariant geologic cap, orthotropic crushable model, MoonsyRivlin rubber, resultant plasticity, closed form update shell plasticity, and FrazerNash rubber model. The IBM 3090 version does not contain the last two models mentioned.
The hydrodynamic material models determine only the deviatoric stresses. Pressure is determined by one of ten equations of state including linear polynomial, JWL high explosive, Sack "Tuesday" high explosive, Gruneisen, ratio of polynomials, linear polynomial with energy deposition, ignition and growth of reaction in HE, tabulatedÂ Â compaction, tabulated, and TENSOR pore collapse. DYNA3D generates three binary output databases. One contains information for complete states at infrequent intervals; 50 to 100 states is typical. The second contains information for a subset of nodes and elements at frequent intervals; 1,000 to 10,000 states is typical. The last contains interfaces data for contact surfaces.
This package is distributed by: Energy Science and Technology Software Center P.O. Box 62 1 Science.Gov Way Oak Ridge, TN 37831 (865) 5762606 TEL (865) 5766436 FAX Email: ESTSC@osti.gov
4. METHOD OF SOLUTION
A contactimpact algorithm permits gaps and sliding along material interfaces with friction. All versions except for the IBM3090 include an interface type defining oneway treatment of sliding with voids and friction. By a specialization of this algorithm, such interfaces can be rigidly tied to admit variable zoning with no need for transition regions. Spatial discretization is achieved by implementation of HughesLiu rectangular beams and shells, BelytschkoTsay shells and beams, triangular shell elementsÂ Â based on work by Belytschko and colleagues, and 8node solidshell elements. All element classes can be included as parts of a rigid body. Threedimensional plane stress constitutive subroutines update the stress tensor for the shell elements such that the stress component normal to the shell midsurface is zero. One constitutive evaluation is made for each integration point through the shell thickness. The 8node solid element uses either one point integration or the Flanagan and Belytschko constant stress formulation with exact volume integration. Zero energy modes in theÂ Â shell and solid elements are controlled by either an hourglass viscosity of stiffness. The equations of motion are integrated in time by the central difference method. A Jaumann stress rate formulation is used with the exception of the orthotropic elastic and the rubber material subroutines which use GreenSt.Venant strains to compute second PiolaKirchoff stresses which transform to Gauchy stresses.
5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM
Storage allocation is dynamic. The only limit that exists is the storage capacity of the computer. Typical calculations have 10,000 to 200,000 elements.
6. TYPICAL RUNNING TIME
Execution speeds on the Cray XMP for various DYNA3D element types range from 4 to 48 microseconds per element cycle. These timings are approximate and do not account for the inclusion of sliding interfaces or complex material models. TypicalÂ Â problems require from 4 to 80 minutes on the Cray XMP.
7. UNUSUAL FEATURES OF THE PROGRAM:
8. RELATED AND AUXILIARY PROGRAMS
The soil and crushable foam, linear viscoelastic, and the rubber BlatzKo subroutines were adapted fromÂ Â HONDO and recoded for vectorization, the ignition and growth EOS was adapted from KOVEC. The forms of the first five equations of state are from KOVEC also. The threedimensional contactimpact algorithmÂ Â is an extension of the NIKE2D (NESC9923) twodimensional algorithm.Â Â The tied and sliding only interface options are similar to the twoÂ Â dimensional algorithm used in DYNA2D (NESC9910). DYNA3D relies on standalone mesh generators. INGRID (NESC9649) is recommended for creation of input files, since it provides complete support for all slide surface data, boundary conditions, loads, material properties, and control parameters. The compatible release of TAURUS (NESC9908) processes DYNA3D output plotting contours, fringes, time histories, and deformed shapes. A variety of strain measures, reaction forces along constrained boundaries, and momentaÂ Â can be computed.
9. STATUS
No status
10. REFERENCES
 D.W. Stillman and J.O. Hallquist INGRID: A ThreeDimensional Mesh Generator for Modeling Nonlinear Systems, UCID20506, July 1985.  J.P. Woodruff KOVEC User's Manual, UCID17306, November 23, 1976.  Michael A. Gerhard SLIC: An Interative, Graphic Mesh Generator for FiniteElement and Finite Difference Applications Programs, UCRL52823, September 13, 1979.  DYNA3D, NESC No. 9909.CRA1B, DYNA3D Edition B Cray Tape Description National Energy Software Center Note 9050, March 30, 1990.  DYNA3D, NESC No. 9909.VAX, DYNA3D DEC VAX Version Tape Description National Energy Software Center Note 9051, March 30, 1990.
11. MACHINE REQUIREMENTS:
12. PROGRAMMING LANGUAGE(S) USED
No specified programming language
13. OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED
COS (Cray), VMS (DEC VAX11), VM/MVS (IBM3090), UNIX BSD 4.2 (SUN).
14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS
Presently, only parts of the sliding interface logic are vectorized.
15. NAME AND ESTABLISHMENT OF AUTHORS
R.G. Whirley Lawrence Livermore National Laboratory P.O. Box 808 Livermore, California 94550
T. Kennedy IBM Corporation
16. MATERIAL AVAILABLE
17. CATEGORIES
 I. Deformation and Stress Distributions, Structural Analysis and Engineering Design Studies
Keywords: deformation, dynamic loads, finite element method, hydrodynamics, solids, threedimensional.
