The Düsseldorf Advanced Material Simulation Kit: DAMASK
F. Roters*, C. Reuber*,1, P. Shanthraj*,2, M. Diehl*, P. Eisenlohr*,3, D. Raabe**Department of Microstructure Physics and Alloy Design 1now at: FEV GmbH, Alsdorf, Germany 2Aachen Institute for Advanced Study in Computational Engineering Science, RWTH Aachen University, Germany 3now at: Chemical Engineering & Materials Science, Michigan State University, USA
Crystal plasticity modelling has matured into a cornerstone of computational materials science . Developing this field from its original mean-field homogenization approach in conjunction with viscoplastic constitutive hardening rules into an advanced multiphysics continuum field solution strategy requires a long term initiative. The “Theory and Simulation” group is working in this field since 2000. Starting originally from Kalidindi’s numerical crystal plasticity integration scheme  user subroutines for the commercial FEM packages Abaqus and MSC.Marc were developed. The scientific focus of the group was, however, the development of advanced constitutive models based on dislocation densities as internal state variables, which are capable of providing microstructure based predictions . As simulations of different material classes and on different length scales necessitate the use of constitutive descriptions of varying degree of sophistication, it soon turned out that the numerical implementation was not flexible enough to incorporate different constitutive models. Therefore, within the project “Computational Modelling of Polycrystals” (CMCn, the first joined MPG-FhG project ever, established in 2006) the development of a flexible framework, for crystal plasticity simulations on all scales from the single crystal up to the engineering component, was initiated. The new code is strictly modularized to allow easy incorporation of additional models on all length scales. In this way, two different homogenization schemes and seven different constitutive models for plasticity, including a non-local one treating dislocation fluxes [4, 5], have been incorporated into the code package. The ability to use different constitutive models within a single simulation is a unique feature of the code. In 2010, in collaboration with Prof. R. Lebensohn (LANL, Humboldt awardee at MPIE at that time), a spectral method based boundary value problem solver was added to complement the commercial FEM solvers. In 2011, it was decided to release the code, by now called DAMASK, as free software to the public domain. For that purpose, a website (damask.mpie.de) was launched in September 2011. The science community very well received the idea of a free crystal plasticity code. At least 15 groups, including universities like UCLA, research facilities like LANL as well as companies like Tata Steel (IJmuiden), today use DAMASK worldwide (as it can be freely downloaded we do not know the real number). These groups contribute to the further code development as well as adding new features to DAMASK including new/modified constitutive models . Currently DAMASK is extended towards a multi-field solver. This will allow treating coupled problems, e.g. thermo-mechanical-chemical, in a fully consistent way.
 F. Roters, P. Eisenlohr, L. Hantcherli, D. D. Tjahjanto, T. R. Bieler, D. Raabe: Acta Materialia 58 (2010) 1152 – 1211
 S. R. Kalidindi, C. A. Bronkhorst, L. Anand: J. Mech. Phys. Solids, 40 (1992) 537 – 569.
 A. Ma, F. Roters: Acta Materialia 52 (2004) 3603 – 3612
 C. Kords: PhD Thesis RWTH Aachen 2014
 C. Reuber, P. Eisenlohr, F. Roters, D. Raabe: Acta Materialia 71 (2014) 333 – 348
 D. Cereceda, M. Diehl, F. Roters, P. Shanthraj, D. Raabe, J. M. Perlado, J. Marian: GAMM-Mitt. 38 (2015) 213 – 227