CPFEM-VUMAT

Crystal plasticity finite element code, VUMAT files for Abaqus

Introduction

The repository contains Abaqus crystal plasticity VUMAT code for the following papers,

1. AM-316steel-NC/ Chen, Wen, et al. "Microscale residual stresses in additively manufactured stainless steel." Nature Communications, 10, 1-13, (2019)

2. AM-HEA/ Ren, Jie, et al. "Strong yet ductile nanolamellar high-entropy alloys by additive manufacturing." Nature 608, 62-68, (2022)

3. AM-dislocation-cell/Zhang, Yin, et al. "Modeling of microscale internal stresses in additively manufactured stainless steel." Modelling and Simulation in Materials Science and Engineering, 30, 074001 (2022).

4. Crack/ Baolin, Wang, et al. in preparation.

Code Structure

The VUMAT codes are written in Fortran, other postprocessing scripts are using Python3 and MATLAB. The main crystal plasticity code is vumat_*.f, this VUMAT code requires external input files such as slip system definition and crystal orientation. Please remember to change the path of these files (variables of FILE1, FILE2 and etc.) when using this code. The sample running script can be found in the 'run.sh'

After the calculation, the postprocessing script could be used to get the lattice strain response, such as the commands in 'fit.sh'. The outputs will be written into 'fitting.dat', 'lattice_strain_LD.dat' and 'lattice_strain_TD.dat'.

Contact

Please contact Dr. Yin Zhang [email protected] for more technical details.