WARP3D Changelog

What's new in WARP3D 17.5.6

May 22, 2015
  • The plane option of the constraints command now offers more control over the distance proximity checks and provides additional output on the distance checking. Similar improvements are included in the list command (Section 2.16).
  • A new output capability that simplifies significantly the specification of analyses with many load(time) steps. A file of output commands may be defined that will be executed automatically after each step in a specified list of steps. Such automatic output can also be mixed with additional output commands given after specific load steps as in the past.
  • Asymmetric equation solvers are now available for the systems of equations that derive from the crystal plasticity model. Asymmetric systems have symmetry of the non-zero locations in the as- sembled equations but non-symmetric numerical values. The asymmetric solver requires twice the amount of memory and twice the equation solving time for each global Newton iteration { the savings in solution time can arise from the reduced number of global Newton iterations. The asymmetric option is available in threads only execution via the Intel Pardiso solver direct and iterative options Section 2.10.2).
  • Output strain values for solid element in finite strain/rotation analyses are now given by the summed increments of unrotated strain, delta d rotated to the fixed-global axes using R from the polar de-composition, F = RU. These values match closely those output by Abaqus which uses summed, incrementally rotated values over the loading history. See Section 1.8.6 for additional details.
  • A new output capability to write at files of nodal and element results. These files have the simplest possible structure of a 2D array with no other readable information. The at files are provided in text or stream binary) forms. The text files may be read directly using Python ( loadtxt imported into Excel, Matlab, Mathematica or read easily with C, C++, Fortran programs. The stream files may be read with fromfile in Python; Matlab and Mathematica using binary stream I/O, and using standard stream operators in C, C++ and Fortran (2003 and later).
  • Similar to the at results files, a new output capability is available to write a simple, flat fille for the model description. The .text file includes the number of nodes/elements, coordinates of nodes and element data (type, material assigned to element and nodes to which the element connects). The file structure is designed to be read very efficiently by Python, imported into spreadsheets, read by C, C++, ... programs. Appendix K includes sample Python programs to read this model description file and the various at results files.
  • Manual Chapter 1 has been re-written, updated and expanded with a wider range of example analyses/discussion and an expanded discussion of theoretical developments implemented in the code. A new capability is available to release absolute constraints on nodes during the solution. The corresponding reactions are relaxed uniformly to zero over the specified number of load (time) steps.
  • Node releases are continued across analysis restarts if required and are integrated into the adaptive step subdivision of the nonlinear solution.
  • The capability to write the equilibrium equations to a file at any time during a nonlinear analysis is described in Section 2.10.16. These files may become quite large (10s GB) for large models. Moreover, users may want to read the equations for processing with C and C++ programs. The solver file now uses the access=stream capability now available in Fortran.
  • The manual section on multipoint constraints is updated to clarify the more general command syntax, limitations that the equations must be homogeneous, and assumptions the WARP3D equation solver makes about dependent vs independent displacement components.

New in WARP3D 17.5.0 (Jan 6, 2014)

  • A new capability is available during the global Newton iterations to detect a likely diverging solution. This option is β€œon” by default. When a divergent sequence of iterations is detected, the adaptive solution procedure is triggered immediately (if enabled).
  • A user solution parameters routine is now available. See Section 2.10 and Appendix I for full details. This routine is called before the analysis of each load (time) step. The routine may change (1) values of variables that control the nonlinear solution process and (2) the incremental loading definition or the upcoming step. May be valuable, for example, in complex nonlinear simulations to modify load (time) step definitions on-the-fly based on the convergence history of the last few load steps.
  • All user-routines are now collected into the user routines.f file replacing the older umats.f file. The distributed version of this .f file has example implementations of all user routines.
  • The crystal plasticity (CP) model now has an exact consistent tangent to accelerate convergence of global Newton iterations. Computational efficiency in the CP routines is also improved.
  • The Pardiso equation solver now (automatically) uses threaded reordering to increase performance.
  • The routines to write a restart file and to read a restart file now use a segmental structuring of data
  • which improves performance and eliminates occasional runtime issues with file buffering.
  • Increased to 5,000 the maximum number of loading conditions (patterns).

New in WARP3D 17.4.0 (Jun 13, 2013)

  • New crystal plasticity material model for rate and temperature dependent simulation of microscale plastic flow in metals. Includes options for simple gradient-based geometric hardening, which incorporates the effect of necessary dislocations on the hardening properties of the material.
  • Expanded solution capabilities for MPI execution: (1) the stiffness assembly process is now performed using a distributed approach across the MPI ranks, (2) the hypre solver now includes an option for the BoomerAMG pre-conditioner (parallel implementation of algebraic multigrid).
  • The interface-cohesive elements now support hardcopy and binary packet file results for tractions and displacement jumps.
  • Nodal temperatures may now be requested in printed output, Patran compatible result files and in binary packet files.
  • Older solvers removed. User input requirements for solution parameters are much simpler.
  • Startup process (Bash shell script) for MPI-based (+ threads) executions has improved reliability and speed in starting MPI processes.

New in WARP3D 17.3.2 (Oct 17, 2012)

  • The blocking command has been extended to provide automatic assignment of elements to blocks for the most common analyses performed with WARP3D. The automatic blocking feature is applicable for: (1) threads-only parallel execution on Windows, Linux, Mac OS X, (2) use of the Pardiso sparse equation solver (either direct or iterative). For other, much less common situations (EBE solver and/or MPI+threads execution), the blocking information must be provided in the input as in prior versions of the code.
  • The recently added UMAT feature provides a very convenient means to incorporate new material constitutive behaviors. WARP3D expects the UMAT routine and all supporting routines to be writ- ten using thread-safe principles – parallel processing of element blocks will invariably invoke multiple instances of the UMAT concurrently. Many older UMATs use COMMON blocks (not thread-safe), are often exceedingly complex and cannot be easily re-written to make them thread-safe. An option in the nonlinear analysis parameters command is now available to request serial (1 thread) execution of element blocks using the UMAT during only the stress update for non thread-safe UMATs. See Appendix J.

New in WARP3D 17.3.1 (Oct 2, 2012)

  • Includes UMAT interface compatible with Abaqus, a user-nodal loads routine, mean-dilatation (F-bar) for 8-node brick element (needed by UMAT and coming crystal plasticity models), improved displacement extrapolation algorithms to reduce number of global Newton iterations, parallel-distributed stiffness assembly with MPI version, manual updates, new verification problems, bug fixes. Version 17.2 added the hypre linear solver, cyclic plasticity model with temperature dependent material properties, new parallel execution method

New in WARP3D 17.3 (Aug 29, 2012)

  • WARP3D can now integrate existing Abaqus compatible umat routines to support user-defined material behavior. The umat routines are invoked from within (thread) parallel processing of element blocks. The umat thus has full benefit of parallel processing already designed into WARP3D. See Section 3.11 and Appendix J.
  • Major changes have been made to internal data structures and algorithms that impose non-zero dis- placement and temperature increments over a load step. Testing reveals a generally reduced number of equilibrium iterations for convergence. No changes in the manual.
  • Updated descriptions of element blocking requirements. See Section 2.6 specifically if you make use of computational material models cyclic, mises hydrogen, and umat.

New in WARP3D 17.2.0 (Mar 14, 2012)

  • Includes hypre linear solver, cyclic plasticity model with temperature dependent material properties, new parallel execution method