Publications

Results 1–50 of 69
Skip to search filters

Additive manufacturing: Toward holistic design

Scripta Materialia

Jared, Bradley H.; Aguilo, Miguel A.; Beghini, Lauren L.; Boyce, Brad B.; Clark, Brett W.; Cook, Adam W.; Kaehr, Bryan J.; Robbins, Joshua R.

Additive manufacturing offers unprecedented opportunities to design complex structures optimized for performance envelopes inaccessible under conventional manufacturing constraints. Additive processes also promote realization of engineered materials with microstructures and properties that are impossible via traditional synthesis techniques. Enthused by these capabilities, optimization design tools have experienced a recent revival. The current capabilities of additive processes and optimization tools are summarized briefly, while an emerging opportunity is discussed to achieve a holistic design paradigm whereby computational tools are integrated with stochastic process and material awareness to enable the concurrent optimization of design topologies, material constructs and fabrication processes.

More Details

ALEGRA : version 4.6

Wong, Michael K.; Brunner, Thomas A.; Garasi, Christopher J.; Haill, Thomas A.; Mehlhorn, Thomas A.; Drake, Richard R.; Hensinger, David M.; Robbins, Joshua R.; Robinson, Allen C.; Summers, Randall M.; Voth, Thomas E.

ALEGRA is an arbitrary Lagrangian-Eulerian multi-material finite element code used for modeling solid dynamics problems involving large distortion and shock propagation. This document describes the basic user input language and instructions for using the software.

More Details

ALEGRA Update: Modernization and Resilience Progress

Robinson, Allen C.; Petney, Sharon P.; Drake, Richard R.; Weirs, Vincent G.; Adams, Brian M.; Vigil, Dena V.; Carpenter, John H.; Garasi, Christopher J.; Wong, Michael K.; Robbins, Joshua R.; Siefert, Christopher S.; Strack, Otto E.; Wills, Ann E.; Trucano, Timothy G.; Bochev, Pavel B.; Summers, Randall M.; Stewart, James R.; Ober, Curtis C.; Rider, William J.; Haill, Thomas A.; Lemke, Raymond W.; Cochrane, Kyle C.; Desjarlais, Michael P.; Love, Edward L.; Voth, Thomas E.; Mosso, Stewart J.; Niederhaus, John H.

Abstract not provided.

ALEGRA: User Input and Physics Descriptions Version 4.2

Boucheron, Edward A.; Haill, Thomas A.; Peery, James S.; Petney, Sharon P.; Robbins, Joshua R.; Robinson, Allen C.; Summers, Randall M.; Voth, Thomas E.; Wong, Michael K.; Brown, Kevin H.; Budge, Kent G.; Burns, Shawn P.; Carroll, Daniel E.; Carroll, Susan K.; Christon, Mark A.; Drake, Richard R.; Garasi, Christopher J.

ALEGRA is an arbitrary Lagrangian-Eulerian finite element code that emphasizes large distortion and shock propagation. This document describes the user input language for the code.

More Details

An extended finite element method formulation for modeling the response of polycrystalline materials to dynamic loading

AIP Conference Proceedings

Robbins, Joshua R.; Voth, Thomas E.

The extended Finite Element Method (X-FEM) is a finite-element based discretization technique developed originally to model dynamic crack propagation [1]. Since that time the method has been used for modeling physics ranging from static meso-scale material failure to dendrite growth. Here we adapt the recent advances of Vitali and Benson [2] and Song et. al. [3] to model dynamic loading of a polycry stalline material. We use demonstration problems to examine the method's efficacy for modeling the dynamic response of polycrystalline materials at the meso-scale. Specifically, we use the X-FEM to model grain boundaries. This approach allows us to i) eliminate ad-hoc mixture rules for multi-material elements and ii) avoid explicitly meshing grain boundaries. © 2007 American Institute of Physics.

More Details

Concurrent Shape and Topology Optimization

Robbins, Joshua R.; Alberdi, Ryan A.; Clark, Brett W.

The typical topology optimization workflow uses a design domain that does not change during the optimization process. Consequently, features of the design domain, such as the location of loads and constraints, must be determined in advance and are not optimizable. A method is proposed herein that allows the design domain to be optimized along with the topology. This approach uses topology and shape derivatives to guide nested optimizers to the optimal topology and design domain. The details of the method are discussed, and examples are provided that demonstrate the utility of this approach.

More Details

Coupling volume-of-fluid based interface reconstructions with the extended finite element method

Computer Methods in Applied Mechanics and Engineering

Voth, Thomas E.; Mosso, Stewart J.; Robbins, Joshua R.

Here, we examine the coupling of the patterned-interface-reconstruction (PIR) algorithm with the extended finite element method (X-FEM) for general multi-material problems over structured and unstructured meshes. The coupled method offers the advantages of allowing for local, element-based reconstructions of the interface, and facilitates the imposition of discrete conservation laws. Of particular note is the use of an interface representation that is volume-of-fluid based, giving rise to a segmented interface representation that is not continuous across element boundaries. In conjunction with such a representation, we employ enrichment with the ridge function for treating material interfaces and an analog to Heaviside enrichment for treating free surfaces. We examine a series of benchmark problems that quantify the convergence aspects of the coupled method and examine the sensitivity to noise in the interface reconstruction. Finally, the fidelity of a remapping strategy is also examined for a moving interface problem.

More Details

Exploring wave propagation in heterogeneous metastructures using the relaxed micromorphic model

Journal of the Mechanics and Physics of Solids

Alberdi, Ryan A.; Robbins, Joshua R.; Walsh, Timothy W.; Dingreville, Remi P.

Metamaterials are artificial structures that can manipulate and control sound waves in ways not possible with conventional materials. While much effort has been undertaken to widen the bandgaps produced by these materials through design of heterogeneities within unit cells, comparatively little work has considered the effect of engineering heterogeneities at the structural scale by combining different types of unit cells. In this paper, we use the relaxed micromorphic model to study wave propagation in heterogeneous metastructures composed of different unit cells. We first establish the efficacy of the relaxed micromorphic model for capturing the salient characteristics of dispersive wave propagation through comparisons with direct numerical simulations for two classes of metamaterial unit cells: namely phononic crystals and locally resonant metamaterials. We then use this model to demonstrate how spatially arranging multiple unit cells into metastructures can lead to tailored and unique properties such as spatially-dependent broadband wave attenuation, rainbow trapping, and pulse shaping. In the case of the broadband wave attenuation application, we show that by building layered metastructures from different metamaterial unit cells, we can slow down or stop wave packets in an enlarged frequency range, while letting other frequencies through. In the case of the rainbow-trapping application, we show that spatial arrangements of different unit cells can be designed to progressively slow down and eventually stop waves with different frequencies at different spatial locations. Finally, in the case of the pulse-shaping application, our results show that heterogeneous metastructures can be designed to tailor the spatial profile of a propagating wave packet. Collectively, these results show the versatility of the relaxed micromorphic model for effectively and accurately simulating wave propagation in heterogeneous metastructures, and how this model can be used to design heterogeneous metastructures with tailored wave propagation functionalities.

More Details

LDRD final report : mesoscale modeling of dynamic loading of heterogeneous materials

Robbins, Joshua R.; Dingreville, Remi P.; Voth, Thomas E.; Furnish, Michael D.

Material response to dynamic loading is often dominated by microstructure (grain structure, porosity, inclusions, defects). An example critically important to Sandia's mission is dynamic strength of polycrystalline metals where heterogeneities lead to localization of deformation and loss of shear strength. Microstructural effects are of broad importance to the scientific community and several institutions within DoD and DOE; however, current models rely on inaccurate assumptions about mechanisms at the sub-continuum or mesoscale. Consequently, there is a critical need for accurate and robust methods for modeling heterogeneous material response at this lower length scale. This report summarizes work performed as part of an LDRD effort (FY11 to FY13; project number 151364) to meet these needs.

More Details

LDRD Final Report: Topology Optimization for Nonlinear Transient Applications Using a Minimally Invasive Approach

Robbins, Joshua R.

The purpose of this project was to devise, implement, and demonstrate a method that can use Sandias existing analysis codes (e.g., Sierra, Alegra, the CTH hydro code) with minimal modification to generate objective function gradients for optimization-based design in tran- sient, non-linear, coupled-physics applications. The approach uses a Moving Least Squares representation of the geometry to substantially reduce the number of geometric degrees of freedom. A Multiple-Program Multiple-Data computing model is then used to compute ob- jective gradients via finite differencing. Details of the formulation and implementation are provided, and example applications are presented that show effectiveness and scalability of the approach.

More Details

Microstructural modeling of ferroic switching and phase transitions in PZT

Proceedings of SPIE - The International Society for Optical Engineering

Robbins, Joshua R.; Khraishi, Tariq A.; Chaplya, Pavel

Niobium doped Lead Zirconate Titanate (PZT) with a Zr/Ti ratio of 95/5 (i.e., PZT 95/5-2Nb) is a ferroelectric with a rhombohedral structure at room temperature. A crystal (or a subdomain within a crystal) exhibits a spontaneous polarization in any one of eight crystallographically equivalent directions. Such a material becomes polarized when subjected to a large electric field. When the electric field is removed, a remanent polarization remains and a bound charge is stored. A displacive phase transition from a rhombohedral ferroelectric phase to an orthorhombic anti-ferroelectric phase can be induced with the application of a mechanical load. When this occurs, the material becomes depoled and the bound charge is released. The polycrystalline character of PZT 95/5-2Nb leads to highly non-uniform fields at the grain scale. These local fields lead to very complex material behavior during mechanical depoling that has important implications to device design and performance. This paper presents a microstructurally based numerical model that describes the 3D non-linear behavior of ferroelectric ceramics. The model resolves the structure of polycrystals directly in the topology of the problem domain and uses the extended finite element method (X-FEM) to solve the governing equations of electromechanics. The material response is computed from anisotropic single crystal constants and the volume fractions of the various polarization variants (i.e., three variants for rhombohedral anti-ferroelectric and eight for rhomobohedral ferroelectric ceramic). Evolution of the variant volume fractions is governed by the minimization of internally stored energy and accounts for ferroelectric and ferroelastic domain switching and phase transitions in response to the applied loads. The developed model is used to examine hydrostatic depoling in PZT 95/5-2Nb.

More Details

Modeling dislocations in a polycrystal using the generalized finite element method

Robbins, Joshua R.

Modeling the interaction of dislocations with internal boundaries and free surfaces is essential to understanding the effect of material microstructure on dislocation motion. However, discrete dislocation dynamics methods rely on infinite domain solutions of dislocation fields which makes modeling of heterogeneous materials difficult. A finite domain dislocation dynamics capability is under development that resolves both the dislocation array and polycrystalline structure in a compatible manner so that free surfaces and material interfaces are easily treated. In this approach the polycrystalline structure is accommodated using the GFEM, and the displacement due to the dislocation array is added to the displacement approximation. Shown in figure 1 are representative results from simulations of randomly placed and oriented dislocation sources in a cubic nickel polycrystal. Each grain has a randomly assigned (unique) material basis, and available glide planes are chosen accordingly. The change in basis between neighboring grains has an important effect on the motion of dislocations since the resolved shear on available glide planes can change dramatically. Dislocation transmission through high angle grain boundaries is assumed to occur by absorption into the boundary and subsequent nucleation in the neighboring grain. Such behavior is illustrated in figure 1d. Nucleation from the vertically oriented source in the bottom right grain is due to local stresses from dislocation pile-up in the neighboring grain. In this talk, the method and implementation is presented as well as some representative results from large scale (i.e., massively parallel) simulations of dislocation motion in cubic nano-domain nickel alloy. Particular attention will be paid to the effect of grain size on polycrystalline strength.

More Details

Multimode Metastructures: Novel Hybrid 3D Lattice Topologies

Boyce, Brad B.; Garland, Anthony G.; White, Benjamin C.; Jared, Bradley H.; Conway, Kaitlynn C.; Adstedt, Katerina A.; Dingreville, Remi P.; Robbins, Joshua R.; Walsh, Timothy W.; Alvis, Timothy A.; Branch, Brittany A.; Kaehr, Bryan J.; Kunka, Cody; Leathe, Nicholas L.

With the rapid proliferation of additive manufacturing and 3D printing technologies, architected cellular solids including truss-like 3D lattice topologies offer the opportunity to program the effective material response through topological design at the mesoscale. The present report summarizes several of the key findings from a 3-year Laboratory Directed Research and Development Program. The program set out to explore novel lattice topologies that can be designed to control, redirect, or dissipate energy from one or multiple insult environments relevant to Sandia missions, including crush, shock/impact, vibration, thermal, etc. In the first 4 sections, we document four novel lattice topologies stemming from this study: coulombic lattices, multi-morphology lattices, interpenetrating lattices, and pore-modified gyroid cellular solids, each with unique properties that had not been achieved by existing cellular/lattice metamaterials. The fifth section explores how unintentional lattice imperfections stemming from the manufacturing process, primarily sur face roughness in the case of laser powder bed fusion, serve to cause stochastic response but that in some cases such as elastic response the stochastic behavior is homogenized through the adoption of lattices. In the sixth section we explore a novel neural network screening process that allows such stocastic variability to be predicted. In the last three sections, we explore considerations of computational design of lattices. Specifically, in section 7 using a novel generative optimization scheme to design novel pareto-optimal lattices for multi-objective environments. In section 8, we use computational design to optimize a metallic lattice structure to absorb impact energy for a 1000 ft/s impact. And in section 9, we develop a modified micromorphic continuum model to solve wave propagation problems in lattices efficiently.

More Details
Results 1–50 of 69
Results 1–50 of 69