Modeling grain boundaries in polycrystals using cohesive elements: Qualitative and quantitative analysis

We propose and demonstrate several tests to estimate the performance of the cohesive elements in ABAQUS for modeling grain boundaries in complex spatial structures such as polycrystalline aggregates. The performance of the cohesive elements is checked by comparing the computed stresses with the theoretically predicted values for a homogeneous material under uniaxial tensile loading. Statistical analyses are performed under different loading conditions for two elasto-plastic models of the grains: isotropic elasticity with isotropic hardening plasticity and anisotropic elasticity with crystal plasticity. Tests are conducted on an analytical finite element model generated from Voronoi tessellation as well as on a realistic finite element model of a stainless steel wire. The results of the analyses highlight several issues related to the computation of normal and shear stresses. The most severe issue is found within the plastic grain response where the computed normal stresses on a particularly oriented cohesive elements are significantly underestimated. Other issues are found to be related to topological constraints in the modeling space and result in the increased scatter of the computed stresses.