Publications

Abstract not provided.

Sylgard® 184/Glass Microballoon (GMB) potting material is currently used in many NW systems. Analysts need a macroscale constitutive model that can predict material behavior under complex loading and damage evolution. To address this need, ongoing modeling and experimental efforts have focused on study of damage evolution in these materials. Micromechanical finite element simulations that resolve individual GMB and matrix components promote discovery and better understanding of the material behavior. With these simulations, we can study the role of the GMB volume fraction, time-dependent damage, behavior under confined vs. unconfined compression, and the effects of partial damage. These simulations are challenging and push the boundaries of capability even with the high performance computing tools available at Sandia. We summarize the major challenges and the current state of this modeling effort, as an exemplar of micromechanical modeling needs that can motivate advances in future computing efforts.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Polyurethane foams are used widely for encapsulation and structural purposes because they are inexpensive, straightforward to process, amenable to a wide range of density variations (1 lb/ft3 - 50 lb/ft3), and able to fill complex molds quickly and effectively. Computational model of the filling and curing process are needed to reduce defects such as voids, out-of-specification density, density gradients, foam decomposition from high temperatures due to exotherms, and incomplete filling. This paper details the development of a computational fluid dynamics model of a moderate density PMDI structural foam, PMDI-10. PMDI is an isocyanate-based polyurethane foam, which is chemically blown with water. The polyol reacts with isocyanate to produces the polymer. PMDI- 10 is catalyzed giving it a short pot life: it foams and polymerizes to a solid within 5 minutes during normal processing. To achieve a higher density, the foam is over-packed to twice or more of its free rise density of 10 lb/ft3. The goal for modeling is to represent the expansion, filling of molds, and the polymerization of the foam. This will be used to reduce defects, optimize the mold design, troubleshoot the processed, and predict the final foam properties. A homogenized continuum model foaming and curing was developed based on reaction kinetics, documented in a recent paper; it uses a simplified mathematical formalism that decouples these two reactions. The chemo-rheology of PMDI is measured experimentally and fit to a generalized- Newtonian viscosity model that is dependent on the extent of cure, gas fraction, and temperature. The conservation equations, including the equations of motion, an energy balance, and three rate equations are solved via a stabilized finite element method. The equations are combined with a level set method to determine the location of the foam-gas interface as it evolves to fill the mold. Understanding the thermal history and loads on the foam due to exothermicity and oven curing is very important to the results, since the kinetics, viscosity, and other material properties are all sensitive to temperature. Results from the model are compared to experimental flow visualization data and post-test X-ray computed tomography (CT) data for the density. Several geometries are investigated including two configurations of a mock structural part and a bar geometry to specifically test the density model. We have found that the model predicts both average density and filling profiles well. However, it under predicts density gradients, especially in the gravity direction. Further model improvements are also discussed for future work.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Previous numerical studies of Sylgard filled with glass microballoons (GMB) have relied on various microstructure idealizations to achieve a large range of volume fractions with high mesh quality. This study investigates how different microstructure idealizations and constraints affect the apparent homogenized elastic constants in the virgin state of the material, in which all GMBs are intact and perfectly bonded to the Sylgard matrix, and in the fully damaged state of the material in which all GMBs are destroyed. In the latter state, the material behaves as an elastomeric foam. Four microstructure idealizations are considered relating to how GMBs are packed into a representative volume element (RVE): (1) no boundary penetration nor GMB-GMB overlap, (2) GMB-GMB overlap, (3) boundary penetration, and (4) boundary penetration and GMB-GMB overlap. First order computational homogenization with kinematically uniform displacement boundary conditions (KUBCs) was employed to determine the homogenized (apparent) bulk and shear moduli for the four microstructure idealizations in the intact and fully broken GMB material states. It was found that boundary penetration has a significant effect on the shear modulus for microstructures with intact GMBs, but that neither boundary penetration nor GMB overlap have a significant effect on homogenized properties for microstructures with fully broken GMBs. The primary conclusion of the study is that future investigations into Sylgard/GMB micromechanics should either force GMBs to stay within the RVE fully and/or use periodic BCs (PBCs) to eliminate the boundary penetration issues. The implementation of PBCs requires the improvement of existing tools in Sandia’s Sierra/SM code.

Abstract not provided.

Kinetic models have been developed to understand the manufacturing of polymeric foams, which evolve from low viscosity Newtonian liquids, to bubbly liquids, finally producing solid foam. Closed-form kinetics are formulated and parameterized for PMDI-10, a fast curing polyurethane, including polymerization and foaming. PMDI-10 is chemically blown, where water and isocyanate react to form carbon dioxide. The isocyanate reacts with polyol in a competing reaction, producing polymer. Our approach is unique, although it builds on our previous work and the polymerization literature. This kinetic model follows a simplified mathematical formalism that decouples foaming and curing, including an evolving glass transition temperature to represent vitrification. This approach is based on IR, DSC, and volume evolution data, where we observed that the isocyanate is always in excess and does not affect the kinetics. The kinetics are suitable for implementation into a computational fluid dynamics framework, which will be explored in subsequent articles. © 2017 American Institute of Chemical Engineers AIChE J, 63: 2945–2957, 2017.

This work was done to support customer questions about whether a Sylgard/Glass Microballoon (GMB) potting material in current use could be replaced with pure Sylgard and if this would significantly change stresses imparted to internal components under thermal cycling conditions. To address these questions, we provide micromechanics analysis of Sylgard/GMB materials using both analytic composite theory and finite element simulations to better understand the role of the GMB volume fraction in determining thermal expansion coefficient, elastic constants, and behavior in both confined and unconfined compression boundary value problems. A key finding is that damage accumulation in the material from breakage of GMBs significantly limits the global stress magnitude and results in a plateau stress behavior over large ranges of compressive strain. The magnitude of this plateau stress is reduced with higher volume fractions of GMBs. This effect is particularly pronounced in confined compression, which we estimate bears the most similarity to the application of interest. This stress-limiting damage mechanism is not present in pure Sylgard, however, and the result is much higher stresses under confined compression. Thus, we recommend that some volume fraction greater than 10% GMBs be used for confined deformation applications.

We present selected results from a series of Open Stack thermal battery tests performed in FY14 and FY15 and discuss our findings. These tests were meant to provide validation data for the comprehensive thermal battery simulation tools currently under development in Sierra/Aria under known conditions compared with as-manufactured batteries. We are able to satisfy this original objective in the present study for some test conditions. Measurements from each test include: nominal stack pressure (axial stress) vs. time in the cold state and during battery ignition, battery voltage vs. time against a prescribed current draw with periodic pulses, and images transverse to the battery axis from which cell displacements are computed. Six battery configurations were evaluated: 3, 5, and 10 cell stacks sandwiched between 4 layers of the materials used for axial thermal insulation, either Fiberfrax Board or MinK. In addition to the results from 3, 5, and 10 cell stacks with either in-line Fiberfrax Board or MinK insulation, a series of cell-free “control” tests were performed that show the inherent settling and stress relaxation based on the interaction between the insulation and heat pellets alone.

Abstract not provided.

Abstract not provided.

We calibrate a linear thermoviscoelastic model for solid Sylgard 184 (90-10 formulation), a lightly cross-linked, highly flexible isotropic elastomer for use both in Sierra / Solid Mechanics via the Universal Polymer Model as well as in Sierra / Structural Dynamics (Salinas) for use as an isotropic viscoelastic material. Material inputs for the calibration in both codes are provided. The frequency domain master curve of oscillatory shear was obtained from a report from Los Alamos National Laboratory (LANL). However, because the form of that data is different from the constitutive models in Sierra, we also present the mapping of the LANL data onto Sandia’s constitutive models. Finally, blind predictions of cyclic tension and compression out to moderate strains of 40 and 20% respectively are compared with Sandia’s legacy cure schedule material. Although the strain rate of the data is unknown, the linear thermoviscoelastic model accurately predicts the experiments out to moderate strains for the slower strain rates, which is consistent with the expectation that quasistatic test procedures were likely followed. This good agreement comes despite the different cure schedules between the Sandia and LANL data.

Abstract not provided.