Publications

Abstract not provided.

Wrought Al-Ge-Si alloys were designed and produced to ensure dislocation bypass strengthening (“hard pin” precipitates) without significant precipitate cutting/shearing (“soft pin” precipitates). These unusual alloys were processed from the melt, solution heat treated and aged. Aging curves at temperatures of 120, 160, 200 and 240 °C were established and the corresponding precipitate spacings, sizes, and morphologies were measured using TEM. The role of non-shearable precipitates in determining the magnitude of Bauschinger was revealed using large-strain compression/tension tests. The effect of precipitates on the Bauschinger response was stronger than that of grain boundaries, even for these dilute alloys. The Bauschinger effect increases dramatically from the under-aged to the peak aged condition and remains constant or decreases slowly through over-aging. This is consistent with reported behavior for Al-Cu alloys (maximum effect at peak aging) and for other Al alloys (increasing through over-aging) such as Al-Cu-Li, Al 6111, Al 2524, and Al 6013. The Al-Ge-Si alloy response was simulated with three microstructural models, including a novel SD (SuperDislocation) model, to reveal the origins of the Bauschinger effect in dilute precipitation-hardened / bypass alloys. The dominant mechanism is related to the elastic interaction of polarized dislocation arrays (generalized pile-up or bow-out model) at precipitate obstacles. Such effects are ignored in continuum and crystal plasticity models.

Ductile rupture in metals is generally a multi-step process of void nucleation, growth, and coalescence. Particle decohesion and particle fracture are generally invoked as the primary microstructural mechanisms for room-temperature void nucleation. However, because high-purity materials also fail by void nucleation and coalescence, other microstructural features must also act as sites for void nucleation. Early studies of void initiation in high-purity materials, which included post-mortem fracture surface characterization using scanning electron microscopy (SEM) and high-voltage electron microscopy (HVEM) and in-situ HVEM observations of fracture, established the presence of dislocation cell walls as void initiation sites in high-purity materials. Direct experimental evidence for this contention was obtained during in-situ HVEM tensile tests of Be single crystals. Voids between 0.2 and 1 μm long appeared suddenly along dislocation cell walls during tensile straining. However, subsequent attempts to replicate these results in other materials, particularly α -Fe single crystals, were unsuccessful because of the small size of the dislocation cells, and these remain the only published in-situ HVEM observations of void nucleation at dislocation cell walls in the absence of a growing macrocrack. Despite this challenge, other approaches to studying void nucleation in high-purity metals also indicate that dislocation cell walls are nucleation sites for voids.

A parallel, adaptive overlay grid procedure is proposed for use in generating all-hex meshes for stochastic (SVE) and representative (RVE) volume elements in computational materials modeling. The mesh generation process is outlined including several new advancements such as data filtering to improve mesh quality from voxelated and 3D image sources, improvements to the primal contouring method for constructing material interfaces and pillowing to improve mesh quality at boundaries. We show specific examples in crystal plasticity and syntactic foam modeling that have benefitted from the proposed mesh generation procedure and illustrate results of the procedure with several practical mesh examples.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Traditional singularity-based fracture mechanics theories rely on their ability to infer the crack tip driving force (local field) by surveying macroscopic physical magnitudes far from the crack tip (far field). This key capability allows engineers to employ nominal forces or displacements to estimate the potential for stable or unstable crack growth. In the case of heterogeneous or anisotropic materials, traditional fracture approaches are not fully theoretically sound and applications rely on extrapolating methodologies with ad-hoc corrections. This Express Laboratory Directed Research and Development (ELDRD) program employed mesoscale-sensitive finite element simulations to assess the impact of grain size and texture on the crack tip behavior. A dislocation-based crystal plasticity model conveys grain size effects by computing the constraint on dislocation cell structures. We assessed the effects of microstructural variability on multiple displacement-based measurements of the fracture driving forces for crack opening (Mode I) and sliding (Mode II). We also consider multiple microstructural realizations of single phase metals undergoing ductile failure. The results show that grain size and texture affect the applied fracture driving force and can induce a significant Mode II deformation under force and displacement control, which is completely neglected in homogeneous models. A large variability in driving forces upon identical far field applied conditions is attributed to a buffering effects of the microstructure. Furthermore, crack mouth opening displacement is almost insensitive to microstructure, which suggests that experimental measurements using such a magnitude (e.g., plastic hinge model) may underestimate local crack tip driving force variability.

We have conducted molecular dynamics (MD) simulations of quasi-isentropic ramp-wave compression to very high pressures over a range of strain rates from 1011 down to 108 1/s. Using scaling methods, we collapse wave profiles from various strain rates to a master profile curve, which shows deviations when material response is strain-rate dependent. Thus, we can show with precision where, and how, strain-rate dependence affects the ramp wave. We find that strain rate affects the stress-strain material response most dramatically at strains below 20%, and that above 30% strain the material response is largely independent of strain rate. We show good overall agreement with experimental stress-strain curves up to approximately 30% strain, above which simulated response is somewhat too stiff. We postulate that this could be due to our interatomic potential or to differences in grain structure and/or size between simulation and experiment. Strength is directly measured from per-atom stress tensor and shows significantly enhanced elastic response at the highest strain rates. This enhanced elastic response is less pronounced at higher pressures and at lower strain rates.

Abstract not provided.

Abstract not provided.

In this work, we develop a tantalum strength model that incorporates effects of temperature, strain rate and pressure. Dislocation kink-pair theory is used to incorporate temperature and strain rate effects while the pressure dependent yield is obtained through the pressure dependent shear modulus. Material constants used in the model are parameterized from tantalum single crystal tests and polycrystalline ramp compression experiments. It is shown that the proposed strength model agrees well with the temperature and strain rate dependent yield obtained from polycrystalline tantalum experiments. Furthermore, the model accurately reproduces the pressure dependent yield stresses up to 250 GPa. The proposed strength model is then used to conduct simulations of a Taylor cylinder impact test and validated with experiments. This approach provides a physically-based multi-scale strength model that is able to predict the plastic deformation of polycrystalline tantalum through a wide range of temperature, strain and pressure regimes.

Abstract not provided.

Abstract not provided.

In this work, a crystal plasticity-finite element (CP-FE) model is used to investigate the effects of microstructural variability at a notch tip in tantalum single crystals and polycrystals. It is shown that at the macroscopic scale, the mechanical response of single crystals is sensitive to the crystallographic orientation while the response of polycrystals shows relatively small susceptibility to it. However, at the microscopic scale, the local stress and strain fields in the vicinity of the crack tip are completely determined by the local crystallographic orientation at the crack tip for both single and polycrystalline specimens with similar mechanical field distributions. Variability in the local metrics used (maximum von Mises stress and equivalent plastic strain at 3% deformation) for 100 different realizations of polycrystals fluctuates by up to a factor of 2-7 depending on the local crystallographic texture. Comparison with experimental data shows that the CP model captures variability in stress-strain response of polycrystals that can be attributed to the grain-scale microstructural variability. This work provides a convenient approach to investigate fluctuations in the mechanical behavior of polycrystalline materials induced by grain morphology and crystallographic orientations.

Abstract not provided.

The mechanical properties of materials systems are highly influenced by various features at the microstructural level. The ability to capture these heterogeneities and incorporate them into continuum-scale frameworks of the deformation behavior is considered a key step in the development of complex non-local models of failure. In this study, we present a modeling framework that incorporates physically-based realizations of polycrystalline aggregates from a phase field (PF) model into a crystal plasticity finite element (CP-FE) framework. Simulated annealing via the PF model yields ensembles of materials microstructures with various grain sizes and shapes. With the aid of a novel FE meshing technique, FE discretizations of these microstructures are generated, where several key features, such as conformity to interfaces, and triple junction angles, are preserved. The discretizations are then used in the CP-FE framework to simulate the mechanical response of polycrystalline α-iron. It is shown that the conformal discretization across interfaces reduces artificial stress localization commonly observed in non-conformal FE discretizations. The work presented herein is a first step towards incorporating physically-based microstructures in lieu of the overly simplified representations that are commonly used. In broader terms, the proposed framework provides future avenues to explore bridging models of materials processes, e.g. additive manufacturing and microstructure evolution of multi-phase multi-component systems, into continuum-scale frameworks of the mechanical properties.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

In this report, we present a multi-scale computational model to simulate plastic deformation of tantalum and validating experiments. In atomistic/ dislocation level, dislocation kink- pair theory is used to formulate temperature and strain rate dependent constitutive equations. The kink-pair theory is calibrated to available data from single crystal experiments to produce accurate and convenient constitutive laws. The model is then implemented into a BCC crystal plasticity finite element method (CP-FEM) model to predict temperature and strain rate dependent yield stresses of single and polycrystalline tantalum and compared with existing experimental data from the literature. Furthermore, classical continuum constitutive models describing temperature and strain rate dependent flow behaviors are fit to the yield stresses obtained from the CP-FEM polycrystal predictions. The model is then used to conduct hydro- dynamic simulations of Taylor cylinder impact test and compared with experiments. In order to validate the proposed tantalum CP-FEM model with experiments, we introduce a method for quantitative comparison of CP-FEM models with various experimental techniques. To mitigate the effects of unknown subsurface microstructure, tantalum tensile specimens with a pseudo-two-dimensional grain structure and grain sizes on the order of millimeters are used. A technique combining an electron back scatter diffraction (EBSD) and high resolution digital image correlation (HR-DIC) is used to measure the texture and sub-grain strain fields upon uniaxial tensile loading at various applied strains. Deformed specimens are also analyzed with optical profilometry measurements to obtain out-of- plane strain fields. These high resolution measurements are directly compared with large-scale CP-FEM predictions. This computational method directly links fundamental dislocation physics to plastic deformations in the grain-scale and to the engineering-scale applications. Furthermore, direct and quantitative comparisons between experimental measurements and simulation show that the proposed model accurately captures plasticity in deformation of polycrystalline tantalum.