Stochastic local FEM for computational homogenization of heterogeneous materials exhibiting large plastic deformations

Computational homogenization is a powerful tool allowing to obtain homogenized properties of materials on the macroscale from simulations of the underlying microstructure. The response of the microstructure is, however, strongly affected by variations in the microstructure geometry. In particular, we consider heterogeneous materials with randomly distributed non-overlapping inclusions, which radii are also random. In this work we extend the earlier proposed non-deterministic computational homogenization framework to plastic materials, thereby increasing the model versatility and overall realism. We apply novel soft periodic boundary conditions and estimate their effect in case of non-periodic material microstructures. We study macroscopic plasticity signatures like the macroscopic von-Mises stress and make useful conclusions for further constitutive modeling. Simulations demonstrate the effect of the novel boundary conditions, which significantly differ from the standard periodic boundary conditions, and the large influence of parameter variations and hence the importance of the stochastic modeling.

Fully-coupled micro–macro finite element simulations of the Nakajima test using parallel computational homogenization

The Nakajima test is a well-known material test from the steel and metal industry to determine the forming limit of sheet metal. It is demonstrated how FE2TI, our highly parallel scalable implementation of the computational homogenization method FE$$^2$$ 2 , can be used for the simulation of the Nakajima test. In this test, a sample sheet geometry is clamped between a blank holder and a die. Then, a hemispherical punch is driven into the specimen until material failure occurs. For the simulation of the Nakajima test, our software package FE2TI has been enhanced with a frictionless contact formulation on the macroscopic level using the penalty method. The appropriate choice of suitable boundary conditions as well as the influence of symmetry assumptions regarding the symmetric test setup are discussed. In order to be able to solve larger macroscopic problems more efficiently, the balancing domain decomposition by constraints (BDDC) approach has been implemented on the macroscopic level as an alternative to a sparse direct solver. To improve the computational efficiency of FE2TI even further, additionally, an adaptive load step approach has been implemented and different extrapolation strategies are compared. Both strategies yield a significant reduction of the overall computing time. Furthermore, a strategy to dynamically increase the penalty parameter is presented which allows to resolve the contact conditions more accurately without increasing the overall computing time too much. Numerically computed forming limit diagrams based on virtual Nakajima tests are presented.