WARP3D - Release 15.9 - April 2008
[ Previous releases of WARP3D can be found here ]

Static and Dynamic Nonlinear Analysis of Fracture in Solids

WARP3D is under continuing development as a research code for the solution of very large-scale, 3-D solid models subjected to static and dynamic loads. Specific features in the code oriented toward the investigation of fracture in metals include a robust finite strain formulation, a general J-integral computation facility (with inertia, thermal, face loading), interaction integrals for computation of linear-elastic fracture parameters, very general element extinction and node release facilities to model crack growth, nonlinear material models including viscoplastic and cyclic, cohesive elements and cohesive constitutive models, and the Gurson-Tvergaard dilatant plasticity model for void growth. We distribute the source code free of charge under the GNU license.


Availability / Downloads

The Source code, makefiles, theory-user guide, and example problems for WARP3D are all freely available below. The theory-user guide is distributed in .pdf format with extensive bookmarking. An interactive "script" in the distribution guides you through the choices to build an executable for a specific platform. The Windows 2000/XP executable, the IA-32 Linux executable, and the AMD64 Linux executable are included in the distribution for users without support software (e.g. compilers) to build the code.

We distribute WARP3D in source form for use under conditions of the Free Software Foundation standard license (a copy of the license agreement is included in the distribution). We provide two sets of files for the distribution, one oriented towards UNIX/Linux users and one for Windows 2000/XP users. The difference is just in the tools used to compress the files.

The following three links give you access to useful information on WARP3D: For fast network connections, the entire WARP3D distribution (including source code) is available for download as one large .tar.gz or .zip file below. For slower network connections, there are separate .tar.gz files and separate .zip files available in the linked directories below. These files can be downloaded one at a time and then unzipped or untarred.
How to Proceed After Downloading

UNIX users
1. Unzip and then untar the distribution:
gzip -dcvf [filename].tar.gz > [filename].tar
tar -xvf [filename].tar
2. Examine the "README" and the "recent_features" files.
3. Compile WARP3D with the ".../src/Makewarp.bash" script.

Linux users
1. Unzip and simultaneously untar the distribution:
tar -xvzf [filename].tar.gz
2. Examine the "README" and the "recent_features" files.
3. Compile WARP3D with the ".../src/Makewarp.bash" script, or use the ready-to-run executable in the ".../run_linux_*" directory corresponding to your CPU architecture.

Windows users
1. Unzip the distribution:
[right-click on [filename].zip and choose "Extract All..."]
2. Examine the "README" and the "recent_features" files. Both can be viewed in WordPad.
3. A ready-to-run executable is supplied in the "run_windows" directory.


Key Algorithms
  • Incremental-iterative, implicit formulation with full Newton iterations
  • Time history integration with Newmark's beta method
  • Adaptive load step logic to assist global Newton iterations
  • Large-displacements, large-strains
  • Polar decomposition based, rotation neutralized finite-strain plasticity
  • Elastic-predictor, radial return constitutive updates
  • Consistent tangent moduli
  • Non-global, absolute coordinate systems at nodes
  • Element block data architectures throughout for vector/parallel
  • Arbitrary rigid body contact (frictionless)
  • Mesh-tieing with automatic generation of multi-point constraints
  • User-specified multi-point constraints
  • Node release and element extinction for crack growth (adaptive damage accumulation)
  • Simplified definition of eigenstrains to model residual stresses
  • Functionally graded materials (liner-nonlinear) supported by definition of all material properties at structure nodes
  • Much of code written in F-90

Equation Solvers
  • Linear-preconditioned conjugate gradient (LPCG) implemented in an element-by-element architecture. Runs in parallel using MPI on all platforms.
  • Sparse-matrix direct solver with various reordering options. The generic sparse solver runs on all platforms. When vendors provide tuned sparse solvers for their hardware, we include an option to invoke that solver. Current vendor specific solvers include:  SGI (serial and parallel), HP (workstations & servers), IBM (RISC 6000 and SP servers), Dec-Alpha (workstations and servers), Windows PCs.

Hex Elements (8, 9, 12, 15, 20-node isoparametric solids)
  • Finite strains
  • B-bar formulation with hourglass stabilization
  • Lumped mass
  • Body and face loadings, temperature loadings
  • Extensive output options

Tet Elements (4 and 10-node isoparametric solids)
  • Finite strains
  • Lumped mass
  • Body and face loadings, temperature loadings
  • Extensive output options

Cohesive Elements
  • 8 and 16 node, quadrilateral elements to connect 8 and 20 node hex elements
  • 6 and 12 node, triangular elements to connect 4 and 10-node tet elements
  • Large displacements

Material Models
  • Mises Plasticity
    • Bilinear, power-law, segmental uniaxial response
    • Power-law viscoplasticity
    • Finite-strain, small-strain options
    • Mixed isotropic-kinematic hardening
    • Cyclic plasticity with mixed nonlinear kinematic-isotropic hardening
    • Hydrogen effects on local plastic flow
  • Deformation Plasticity
    • Nonlinear (J2) elasticity
    • Linear, then power-law uniaxial response
    • Small-strains, inviscid response
  • Gurson-Tvergaard Dilatant Plasticity
    • Bilinear, power-law, segmental uniaxial matrix response
    • Power-law viscoplasticity
    • Finite-strain, small-strain options
    • Nucleation models (stress, strain controlled)
    • Adaptive stress-update algorithms
  • Cohesive (coupled shear & normal, linear, nonlinear 2 parameter exponential)

Domain Integral Computation
  • Domain Interaction Integral Formulations for J
  • Interaction integrals to compute all 3 stress-intensity factors for linear-elastic models
  • Interaction integrals to compute T-stresses for linear-elastic models
  • Automatic and Manual q-function Definition
  • Inertia Loading of Near-Tip Material Included
  • Crack-Face Tractions Included
  • Temperature Effects Included (supports anisotropic expansion coefficients)
  • Collapsed or blunt crack tip shapes
  • All domain and interaction integrals support models with FGMs

Grack Growth Procedures
  • Element Extinction
    • Uses critical porosity to trigger extinction, or
    • Critical plastic strain defined by Stress Modified Critical Strain criterion
    • Force releases by traction-separation model
    • Force releases by fixed number of steps
    • Adaptive load stepping to control increments of plastic strain or porosity
  • Node Release
    • Use CTOA to trigger node release
    • Any number of initial crack fronts supported
    • Crack fronts may coalesce during solution
    • Initial and growth CTOAs user specified
    • Option to enforce uniform growth over front based on critical node
    • Automatic location/tracking of crack front in 3-D
    • Force releases by traction-separation model
    • Force releases by fixed number of steps
    • Adaptive load stepping to control increments of CTOA

Program Input
  • English-like, fully free-form commands
  • Interactive and batch execution modes
  • Translator from MSC-Patran neutral file format to WARP3D provided (performs element blocking operations, domain decomposition and automatic definition of transition elements between 8 and 20-node brick elements)

Program Output
  • Printed by node/element user defined lists
  • Patran compatbile node/element ASCII and binary files
  • Packet files of requested values in binary format

Large Analysis Features
  • CPU Time Limit - code stops analysis and writes restart file if next step will exceed time limit
  • Fully Automated Restart From A Single File
  • Status Message File For Jobs Running In Batch/Background
  • Stacks of Input Files

Supported Computer Platforms
  • HP B,C, J-Class Systems Running HPUX 11.xx
  • HP Itanium-2 and Higher Systems Running HPUX 11.xx
  • SGI Origin/Octane Systems Running IRIX
  • SGI Altix Systems Running SGI Linux
  • DEC-Alpha Systems Running Digital UNIX
  • IBM RISC-6000 & SP-x Systems Running AIX
  • Intel IA-32 and AMD Systems Running Linux
  • AMD64 and Intel EM64T Systems Running Linux
  • Intel IA-32 Systems Running Win2K/WinXP