Engineering Papers⌕ Search

Engineering topics

Zecevic, Miroslav

Publications and source records attributed to Zecevic, Miroslav.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Green's function fast multipole method for continuum mechanics (SM-FMM)

Solid Mechanics Fast Multipole Method based on elastic Green's function and accelerated using FFTs The code calculates the mechanical fields (stress and strain) for a heterogeneous elasto-plastic material unit cell under quasi-static conditions (no dynamic effects). The fields are calculated using a Green's function method, where the strain field is given by a discrete convolution of Green's operator with an auxiliary stress field. The convolution is calculated using the fast multipole method.

Zecevic, Miroslav↗
DOE OSTI2024

Generalized grain boundary constitutive description implemented in a strain-gradient large-strain FFT-based formulation: Application to nano-metallic laminates

Herein this paper presents a general treatment of grain boundary constitutive behavior in the context of strain-gradient (SG) plasticity, and its numerical implementation in a large-strain (LS) elasto-viscoplastic (EVP) fast Fourier transform (FFT)-based micromechanical model. Two novel grain boundary constitutive equations are proposed, allowing for more accurate description of the Burgers vector flow at the grain boundary. The capabilities of the generalized SG-LS-EVPFFT formulation are illustrated for the case of kink-band formation during layer-parallel compression of nano-metallic laminates (NMLs), requiring consideration of the interaction between dislocations and interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A parallel and performance portable implementation of a full-field crystal plasticity model

We have developed a parallel implementation of an Elasto-Viscoplastic Fast Fourier Transform-based (EVPFFT) micromechanical solver to enable computationally efficient crystal plasticity modeling for polycrystalline materials. Our primary focus lies in achieving performance portability, allowing a single EVPFFT implementation to run optimally on various homogeneous architectures, including multi-core Central Processing Units (CPUs), as well as on heterogeneous computer architectures comprising multi-core CPUs and Graphics Processing Units (GPUs) from different vendors. To accomplish this goal, we have leveraged MATAR, a C++ software library that simplifies the creation and utilization of multidimensional dense or sparse matrix and array data structures. These data structures are designed to be portable across diverse architectures through the use of Kokkos, a performance-portable library. Additionally, we have employed the Message Passing Interface (MPI) to efficiently distribute the computational workload among processors. The heFFTe (Highly Efficient FFT for Exascale) library is used to facilitate the performance portability of the fast Fourier transforms (FFTs) computation. The computational performance of EVPFFT is evaluated and presented in terms of parallel scalability and simulation runtime on different high-performance computing (HPC) architectures. As a result, the utility of the developed framework to efficiently simulate the micro-mechanical fields in polycrystalline microstructures in engineering applications is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Role of grain boundary crystallography on void growth in FCC metals

The growth of an isolated spherical void at and near a grain boundary or a triple junction is studied by means of full-field crystal plasticity simulations with explicit representation of the void. To examine the fundamental aspects of void/grain boundary interactions in these cases, a three-dimensional large-strain elasto-viscoplastic fast-Fourier transform (LS-EVP-FFT) model with axisymmetric tensile periodic boundary conditions was developed. We employ this model to reveal the role of crystallography (grain orientations, grain boundary misorientation and inclination angle with respect to the largest component of the applied stress) and type of loading on void growth. It was found that grain orientation was the primary crystallographic feature governing void growth in strain-rate controlled simulations. For strains small enough to not significantly alter the local geometry, the overall growth of the void in bicrystal boundaries was not strongly affected by grain boundary inclination angle. The relative amount of void growth experienced within each grain in a bicrystal or tri-crystal was strongly dependent on the grain boundary inclination angle. We found a subset of orientations that exhibit faster void growth in bicrystal boundaries than inside the constituent single crystals. Last, we show that when placed primarily in one grain but near the grain boundary, void growth was similar to void growth when placed at the grain boundary. Furthermore, the growing void remained in the original grain and did not grow into the other grain, regardless of crystallographic character.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Unraveling kinking: A plasticity enhancing failure mode in high strength nano metallic laminates

Kinking is an important and plasticity-enhancing deformation/failure mode in numerous mechanically anisotropic materials including high-strength nano metallic laminates (NMLs). However, our current limited understanding of the mechanics of kinking and its dependence on microstructural attributes is insufficient for thoroughly comprehending and eventually being able to control failure behaviors of materials. In this study, we investigate kinking dependencies on microstructural attributes in NMLs via in situ micropillar compression, multiscale microstructure characterization, dislocation dynamic simulations, and crystal plasticity modeling. Additionally, by examining several NML systems (Cu/Fe, Ag/Fe, Al-4Mg/Fe), we demonstrate that the development of internal stresses during loading activates local layer-parallel glide triggering kinking in NMLs. Furthermore, this work reveals the effect of key microstructural features including layer thickness, layer waviness, interface barrier strength, and work hardening capacity on kink band formation in NMLs. More broadly, our efforts represent a generically applicable approach for probing large-strain deformation behavior of complex materials via synergetic modeling and experimental efforts.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Implementation and experimental validation of nonlocal damage in a large-strain elasto-viscoplastic FFT-based framework for predicting ductile fracture in 3D polycrystalline materials

Ductile materials, such as metal alloys, can undergo substantial deformation before failure. Additionally, these materials are usually of polycrystalline composition and exhibit strongly anisotropic behavior at small length scales. Previously developed fast Fourier transform (FFT)-based models can model ductile fracture of isotropic materials or the elastic–plastic behavior of anisotropic polycrystalline materials; however, there remains a need to couple both capabilities. This work extends a large-strain FFT-based crystal plasticity model to simulate ductile fracture of polycrystalline materials. In this work, a triaxiality-based continuum damage mechanics (CDM) formulation is incorporated into a large-strain elasto-viscoplastic FFT (LS-EVPFFT) framework. The CDM formulation is augmented with an integral-based nonlocal regularization approach that correctly handles gas-phase material necessary to model unconstrained surfaces. To validate the damage-enabled LS-EVPFFT framework, mesoscale copper tensile coupons were machined using microwire electrical discharge machining and experimentally characterized using electron backscatter diffraction. In-situ optical digital image correlation was performed during uniaxial testing to provide a side-by-side comparison of the experimental and computational strain fields and stress–strain responses. The damage-enabled LS-EVPFFT framework can simulate the complete macroscopic stress–strain response of ductile polycrystals to failure. The model reproduces necking behavior that qualitatively agrees with experimental observations. By leveraging the relatively low computational cost of the damage-enabled LS-EVPFFT framework, the framework presented here allows the ductile fracture response of 3D polycrystalline materials to be tractably predicted.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Non-local large-strain FFT-based formulation and its application to interface-dominated plasticity of nano-metallic laminates

Here, this paper presents a novel formulation and its robust numerical implementation of strain-gradient (SG) crystal plasticity within a large-strain (LS) elasto-viscoplastic (EVP) fast Fourier transform (FFT)-based micromechanical model. The resulting non-local SG-LS-EVPFFT formulation is used to model and understand the process of kink band formation during layer-parallel compression of nano-metallic laminates (NMLs). NMLs are layered composites with nanoscale thicknesses, thus requiring consideration of the interaction between dislocations and interfaces within the micromechanical model. The length-scale parameter of the SG model is calibrated by simulating a double pile-up and comparing predictions to analytical solution. This required new expressions for the defect energy, resulting in more accurate double pile-up predictions. The calibrated SG-LS-EVPFFT model is then used to simulate layer-parallel compression of copper-niobium NML. Formation of kink bands is predicted, and the model is used to rationalize the microscopic mechanisms enabling the formation process. It is found that accumulation of dislocations at interfaces leads to activation of layer-parallel slip, which in turn leads to kink band formation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Predicting extreme anisotropy and shape variations in impact testing of tantalum single crystals

We report recent Taylor cylinder impact tests carried out for Ta single crystals showed strong variations in dimensional changes for different crystallographic directions aligned with the cylindrical axis. In order to capture the effect of crystallography on the deformation characteristics and final shapes of the impacted cylinders, a single crystal material subroutine is adapted and embedded in the solid mechanics/dynamics Finite Element solver Abaqus to simulate the aforementioned single crystal Ta Taylor impact experiments. Details of the coupled model implementation, and insights on the role played by single crystal anisotropic flow on the deformation behavior across a broad range of strain rates and temperatures for different single crystal orientations are presented and discussed. We demonstrate the predictive capability of the adopted crystal plasticity model to capture the significant role played by crystal orientation-induced anisotropy, as well as strain hardening and adiabatic heating, on the dynamic deformation response of crystalline materials. This re-emphasizes the need of microstructure-aware models to improve the accuracy of simulations for high-consequence engineering design.

36 MATERIALS SCIENCE↗
DOE OSTI2022

New large-strain FFT-based formulation and its application to model strain localization in nano-metallic laminates and other strongly anisotropic crystalline materials

This paper presents a new robust large-strain (LS) elasto-viscoplastic (EVP) formulation based on Fast Fourier Transforms (FFTs) for the prediction of the micro-mechanical response and microstructure evolution of polycrystalline and multiphase materials, with emphasis on the effect of strong crystallographic and/or morphologic anisotropy on localization of plastic deformation. In this work, the novel LS-EVPFFT formulation allows treatment of complex initial geometries and large deformations considering three grids of material points: a regular grid in the reference configuration, where FFTs can be performed; an irregular grid in the initial configuration, created by applying a stress-free displacement field to the reference regular grid; and an irregular grid in the current configuration, undergoing large strains and rotations as the material is loaded. Further numerical stability of the new formulation also required the use of a novel expression for the discrete modified Green’s operator, which reduces spurious field oscillations. After presenting and validating the new formulation by comparison with preexisting implementations and analytical solutions, LS-EVPFFT is applied to the prediction of slip and kink bands formation in polycrystalline columnar ice, and kink bands in single crystal zinc wires, showing good agreement with classic experiments. Finally, the model is used to study kink band formation during compression of Cu–Nb nano-metallic laminates (NMLs), in which accurate treatment of the complex geometry associated with the tortuosity of interfaces and large deformations become critical, showing consistency with corresponding micropillar experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Approximation of periodic Green's operator in real space using numerical integration and its use in fast Fourier transform-based micromechanical models

In this work, we propose an expression for the periodic first derivative of Green's function in real space. The proposed expression allows an alternative way of computing the periodic Green's operator based on periodically summing the free-space Green's operator in terms of an appropriate quadrature rule. We provide computational examples, which show the accuracy of the proposed approach, together with reduced spurious oscillations in the solution fields.

42 ENGINEERING↗
DOE OSTI2020

New self-consistent homogenization for thermo-elastic polycrystals with imperfect interfaces

Here we present a new thermo-elastic self-consistent formulation to calculate the thermo-mechanical response of polycrystalline aggregates with imperfect interfaces. The new model is based on an improved treatment of the effect of imperfect interfaces on the micromechanical response of the grains compared to pre-existing models. New expressions/algorithms are derived for the effective elastic and thermal properties, localization relations, and calculation of intragranular field fluctuations. In addition, derivation of a novel integral equation for the strain field in a heterogeneous elastic medium subjected to arbitrary thermal strain field and displacement jumps is presented. The new model is applied to study the effective thermo-elastic properties and average field fluctuations as a function of the interface properties for imperfectly-bonded triaminotrinitrobenzene (TATB) polycrystalline granular aggregates. The proposed formulation is general in the sense that it can be extended to other material systems and deformation regimes that require a proper consideration of the effect of imperfect interface properties on the material's thermo-mechanical response.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Defects in epitaxial Ru(0001) on Al2O3(0001): Dislocations, stacking faults, and deformation twins

In this study, defects in epitaxial Ru(0001) films on c-plane sapphire, with nominal thicknesses of 10–80 nm, deposited at 350 °C and step-annealed to 950 °C, were characterized using transmission electron microscopy. The variation of Ru and sapphire lattice parameters with temperature is such that the misfit strain for the observed 30° rotated-honeycomb epitaxial relationship is essentially constant with temperature at 1.5%, resulting in a biaxial stress of 10.0 GPa and an energy density of 150 MJ m –3 in unrelaxed films. Stress relaxation occurs by the formation of defects. For the 20–80 nm thick films, the defects are a- and c-type dislocations and stacking faults, argued to be of I 2 type. In addition, the films show the surprising presence of $\{11\bar{2}1\}1/3\langle11\bar{2}\bar{6}\rangle$ deformation twins. The 10 nm-thick films were found to be defect free. The critical thickness for misfit strain relaxation via the formation of threading and misfit dislocations is computed as 7±2 nm, depending on the choice of the dislocation core radius. Energetic analysis of twin formation, using both the infinite-matrix and the finite-matrix (Mori–Tanaka) approaches, provides values of the twin aspect ratios, assumed to be ellipsoidal, and shows that the latter but not the former approach can qualitatively explain the formation of the observed twins. In addition to providing the maximum strain relief compared to other potential twin types, $\{11\bar{2}1\}1/3\langle11\bar{2}\bar{6}\rangle$ twins do not require lattice shuffles and have a boundary that is a special boundary, namely, a 35° tilt boundary with a-type dislocations every other {0002} plane, that may also favor their formation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental verification of a crystal plasticity-based simulation framework for predicting microstructure and geometric shape changes: Application to bending and Taylor impact testing of Zr

This paper is concerned with experimental verification of a recently developed multi-scale simulation framework for plastic deformation of metallic materials from quasi-static to impact deformation conditions. The framework is a visco-plastic self-consistent (VPSC) polycrystalline model embedded in an implicit finite element method (FE-VPSC) to provide a microstructure-sensitive constitutive response at each material point. Each material point of the FEM model is a polycrystalline aggregate with crystallographic deformation mechanisms operating at the single crystal scale with their evolving activity based on a dislocation density-based hardening law and texture. Four beams and three cylinders machined in different orientations from a textured plate of high-purity zirconium are tested quasi-statically in 4-point bending and at speeds of 100 m/s, 170 m/s and 243 m/s during Taylor impact tests, respectively. The variation in dimensional changes resulting from different sample orientations in the plate with respect to loading directions is measured for each sample. Moreover, texture and twinning characterization is performed using electron backscattered diffraction (EBSD). The deformation processes and underlying evolution of microstructure are successfully simulated using the FE-VPSC framework. In doing so, the model parameters are optimized and validated across a broad range of strain rates and temperatures. Simulation results in terms of geometrical changes and microstructural evolution are compared with the experimental measurements. The model predicts anisotropic material flow resulting from the hard-to-deform crystallographic directions, the development of gradients in texture and twinning through the geometries, tension–compression asymmetry, as well as the extent of plasticity under impact.

42 ENGINEERING↗
DOE OSTI2020

Modeling of the thermo-mechanical response and texture evolution of WE43 Mg alloy in the dynamic recrystallization regime using a viscoplastic self-consistent formulation

This paper presents a microstructure sensitive model for predicting mechanical response and texture evolution of metals in the dynamic recrystallization regime. A recently proposed viscoplastic self-consistent (VPSC) formulation for the prediction of recrystallization driven by strain energy and intragranular misorientation is extended to hexagonal close-packed (hcp) metals. The model is applied to the dynamic recrystallization of magnesium alloy WE43 at different temperatures and strain rates. Model predictions in terms of stress-strain response and texture evolution are compared to the experimental measurements and acceptable agreement is achieved. According to the model predictions, superplastic behavior of nuclei was found to be the dominant softening mechanism at high temperatures and low strain rates. High concentration of precipitates at the grain boundaries and presence of alloying elements are the likely causes of low boundary mobility, resulting in nucleation dominated dynamic recrystallization. Relatively strong basal compression textures indicate dominant activity of basal slip, which can be achieved only through large difference in slip resistance between soft basal and hard prismatic and pyramidal modes.

36 MATERIALS SCIENCE↗