Software details

Software authors

Javier Segurado, Javier LLorca, Sergio Lucarini, Aitor Cruzado, Sarra Haouala, Vicente Herrera

Intellectual Property Rights

Copyright © 2012-2021 Fundación IMDEA Materiales and Universidad Politécnica de Madrid. All rights reserved

Transfer opportunity

Software license

Software description

CAPSUL is a package of crystal plasticity and polycrystalline homogenization simulation tools. The package includes:

  • A crystal plasticity constitutive model, aimed at predicting the elasto-plastic behaviour at the crystal level. The model incorporates physically-based and phenomenological implementations and both slip and twinning mechanisms are considered. Both monotonic response and cyclic behaviour are considered by the combination of different laws for isotropic hardening, kinematic hardening and cyclic softening. The model has already been successfully applied to simulate the mechanical performance of FCC (Al and Ni-based superalloys) and HCP alloys (Mg and Ti), including stress-strain curves and texture evolution under monotonic loading as well as crack initiation and fatigue life under cyclic loading. The model is programmed as a UMAT subroutine for Abaqus.
  • A tool to generate finite element models representative of the actual alloy microstructure (grain size, shape and orientation distribution) using as input statistical data obtained from microscopy images.
  • An inverse optimization tool to obtain the crystal plasticity model parameters from the result of a set of mechanical tests (both microtests or tests on polycrystals).
  • A set of python scripts to generate cyclic loading conditions and to postprocess the results to obtain fatigue indicator parameters and other measures from microfields.


CAPSUL capabilities

Representative volume elements of polycrystalline Ti. (a) Voxel model with 1000 cubic finite elements in which each one stands for a single crystal. (b) Realistic RVE containing 100 crystals discretized with 64000 cubic finite elements.
Distribution of the third Euler angle in the microstructure of an RVE of polycrystalline Ti. Picture axis corresponds to x = x1,y = x2 and z = x3 and contour plot is drawn on the reference configuration (a) Initial data before deformation showing the grain structure. (b) After drawing along the z = x3 axis.
Contour plot of the accumulated plastic slip Γ in the microstructure of polycrystalline Ti after 100% of tensile strain. (a) RVE containing 1000 crystals, each one represented by one cubic finite element. (b) RVE containing 100 crystals, each one represented with 640 cubic finite elements.


Knowledge & Technology Transfer Department, IMDEA Materials Institute


telephone: +34 91 5493422