Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.
Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.
Constitutive model parameterizations for the General Plastics EF4003 low density 3 pound per cubic foot are needed for design and qualification purposes in normal and abnormal mechanical simulations. The material is expected to be deformed in two ways: first during preloading, and second under impact conditions of the system (transient dynamic). All analyses are to be performed at room temperature. The goal is to provide the analysis community a robust constitutive model parameterization to represent the compression behavior of the EF4003 foam from small deformations up to massive compressive deformations when the foam is densifying. It is worth noting the EF4003 exhibits anisotropy in its stress-strain behavior between the rise and transverse directions (See figure 2.8c-d) as well as plateau behavior that is very likely to cause material stability issues, due to the buckling transition, (and has historically done so) when using Sandia’s current workhorse models for flexible foams, Hyperfoam and Flex Foam. A Stability-informed Hyperfoam parameterization procedure is developed and executed to calibrate a hyperfoam model for the EF4003 room temperature, transversely loaded data. A rise orientation parameterization was not attempted due to localization in the experiments.
Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.
Luo, Chaoqian; Chung, Christopher; Yakacki, Christopher M.; Long, Kevin N.; Yu, Kai
Liquid crystal elastomers (LCEs) exhibit soft elasticity due to the alignment and reorientation of mesogens upon mechanical loading, which provides additional mechanisms to absorb and dissipate energy. This enhanced response makes LCEs potentially transformative materials for biomedical devices, tissue replacements, and protective equipment. However, there is a critical knowledge gap in understanding the highly rate-dependent dissipative behaviors of LCEs due to the lack of real-time characterization techniques that probe the microscale network structure and link it to the mechanical deformation of LCEs. In this work, we employ in situ optical measurements to evaluate the alignment and reorientation degree of mesogens in LCEs. The data are correlated to the quantitative physical analysis using polarized Fourier-transform infrared spectroscopy. The time scale of mesogen alignment is determined at different strain levels and loading rates. The mesogen reorientation kinetics is characterized to establish its relationship with the macroscale tensile strain, and compared to theoretical predictions. Overall, this work provides the first detailed study on the time-dependent evolution of mesogen alignment and reorientation in deformed LCEs. It also provides an effective and more accessible approach for other researchers to investigate the structural-property relationships of different types of polymers.
Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.
Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.
This SAND report fulfills the completion requirements for the ASC Physics and Engineering Modeling Level 2 Milestone 7836 during Fiscal Year 2021. The Sandia Simplified potential energy clock (SPEC) non-linear viscoelastic constitutive model was developed to predict a whole host of polymer glass physical behaviors in order to provide a tool to assess the effects of stress on these materials over their lifecycle. Polymer glasses are used extensively in applications such as electronics packaging, where encapsulants and adhesives can be critical to device performance. In this work, the focus is on assessing the performance of the model in predicting material evolution associated with long-term physical aging, an area that the model has not been fully vetted in. These predictions are key to utilizing models to help demonstrate electronics packaging component reliability over decades long service lives, a task that is very costly and time consuming to execute experimentally. The initiating hypothesis for the work was that a model calibration process can be defined that enables confidence in physical aging predictions under ND relevant environments and timescales without sacrificing other predictive capabilities. To test the hypothesis, an extensive suite of calibration and aging data was assembled from a combination of prior work and collaborating projects (Aging and Lifetimes as well as the DoD Joint Munitions Program) for two mission relevant epoxy encapsulants, 828DGEBA/DEA and 828DGEBA/T403. Multiple model calibration processes were developed and evaluated against the entire set of data for each material. A qualitative assessment of each calibration's ability to predict the wide range of aging responses was key to ranking the calibrations against each other. During this evaluation, predictions that were identified as non-physical, i.e., demonstrated something that was qualitatively different than known material behavior, were heavily weighted against the calibration performance. Thus, unphysical predictions for one aspect of aging response could generate a lower overall rating for a calibration process even if that process generated better quantitative predictions for another aspect of aging response. This insurance that all predictions are qualitatively correct is important to the overall aim of utilizing the model to predict residual stress evolution, which will depend on the interplay amongst the different material aging responses. The DSC-focused calibration procedure generated the best all-around aging predictions for both materials, demonstrating material models that can qualitatively predict the whole host of different physical aging responses that have been measured. This step forward in predictive capability comes from an unanticipated source, utilization of calorimetry measurements to specify model parameters. The DSC-focused calibration technique performed better than compression-focused techniques that more heavily weigh measurements more closely related to the structural responses to be predicted. Indeed, the DSC-focused calibration procedure was only possible due to recent incorporation of the enthalpy and heat capacity features into SPEC that was newly verified during this L2 milestone. Fundamentally similar aspects of the two material model calibrations as well as parametric studies to assess sensitives of the aging predictions are discussed within the report. A perspective on the next steps to the overall goal of residual stress evolution predictions under stockpile conditions closes the report.
We have characterized the three-dimensional evolution of microstructural anisotropy of a family of elastomeric foams during uniaxial compression via in-situ X-ray computed tomography. Flexible polyurethane foam specimens with densities of 136, 160 and 240 kg/m3 were compressed in uniaxial stress tests both parallel and perpendicular to the foam rise direction, to engineering strains exceeding 70%. The uncompressed microstructures show slightly elongated ellipsoidal pores, with elongation aligned parallel to the foam rise direction. The evolution of this microstructural anisotropy during deformation is quantified based on the autocorrelation of the image intensity, and verified via the mean intercept length as well as the shape of individual pores. Trends are consistent across all three methods. In the rise direction, the material remains transversely anisotropic throughout compression. Anisotropy initially decreases with compression, reaches a minimum, then increases up to large strains, followed by a small decrease in anisotropy at the largest strains as pores collapse. Compression perpendicular to the foam rise direction induces secondary anisotropy with respect to the compression axis, in addition to primary anisotropy associated with the foam rise direction. In contrast to compression in the rise direction, primary anisotropy initially increases with compression, and shows a slight decrease at large strains. These surprising non-monotonic trends and qualitative differences in rise and transverse loading are explained based on the compression of initially ellipsoidal pores. Microstructural anisotropy trends reflect macroscopic stress-strain and lateral strain response. These findings provide novel quantitative connections between three-dimensional microstructure and anisotropy in moderate density polymer foams up to large deformation, with important implications for understanding complex three-dimensional states of deformation.
Mechanical impact protection is an important consideration in many applications, ranging from product transportation to sports. Cellular materials are typically used due to their desirable energy absorption properties and light weight. However, their large deformation and rate dependent responses (especially of polymer foams) are challenging to consider in design. Additionally, the use of foams with uniform properties, such as uniform density and uniform stiffness, often restricts the designed foams to only be suitable for a narrow range of mechanical impact conditions whereas real applications commonly face unpredictable situations. 3D printing offers fabrication flexibility and thus opens the door to create foams with tailored properties. In this work, we investigate the feasibility of using 3D printing for functionally graded foams (FGFs) that are optimal over a broad range of mechanical environments. The foams are fabricated by the recently developed grayscale digital light processing (g-DLP) method which can print parts with locally designed properties. These foams are tested under drop test conditions and with slower displacement control. We also model the large deformation behavior of FGFs using finite element analysis in which we account for the different viscoelastic behaviors of the distinct grayscale regions. We then use the model to examine the impact mitigation capabilities of FGFs in different loading scenarios. Finally, we show how FGFs can be used to satisfy real-world design goals using the case study of a motorcycle helmet. In contrast to prior work, we investigate continuous, functionally graded foams of a single density that differ in their viscoelastic responses. This work provides further insight into the benefits of viscoelastic properties and modulus graded foams and presents a manufacturing approach that can be used to produce the next generation of flexible lattice foams as mechanical absorbers.
Luo, Chaoqian; Chung, Christopher; Traugutt, Nicholas A.; Yakacki, Christopher M.; Long, Kevin N.; Yu, Kai
Polymer foams are an essential class of lightweight materials used to protect assets against mechanical insults, such as shock and vibration. Two features are important to enhance their energy absorption characteristics: the foam structure and the matrix phase mechanical behavior. This study investigates novel approaches to control both of these features to enhance the energy absorption capability of flexible lattice foams. First, we consider 3D printing via digital light processing (DLP) as a method to control the foam mesostructure across a suite of periodic unit cells. Second, we introduce an additional energy dissipation mechanism in the solid matrix phase material by 3D printing the lattice foams with polydomain liquid crystal elastomer (LCE), which undergo a mechanically induced phase transition under large strains. This phase transition is associated with LC mesogen rotation and alignment and provides a second mechanism for mechanical energy dissipation in addition to the viscoelastic relaxation of the polymer network. We contrast the 3D printed LCE lattices with conventional, thermomechanically near-equivalent elastomer lattice foams to quantify the energy-absorbing enhancement the LCE matrix phase provides. Under cyclic quasi-static uniaxial compression conditions, the LCE lattices show dramatically enhanced energy dissipation in uniaxial compression compared to the non-LCE equivalent foams printed with a commercially available photocurable elastomer resin. The lattice geometry also plays a prominent role in determining the energy dissipation ratio between the LCE and non-LCE foams. We show that when increasing the lattice connectivity, the foam deformation transitions from bending-dominated to stretching-dominated deformations, which generates higher axial strains in the struts and higher energy dissipation in the lattice foam, as stretching allows greater mesogen rotation than bending. The LCE foams demonstrate superior energy absorption during the repeated dynamic loading during drop testing compared with the non-LCE equivalent foams, demonstrating the potential of LCEs to enhance physical protection systems against mechanical impact.
Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented
The purpose of this document is to provide evidence for assessing the adequacy of parameterized material models for a collection of materials used in a finite element analyses setting. “Adequacy” is relative to the intended use of the material in particular analyses. The intended application of the material models covered within this document is for system level abnormal mechanical solid mechanics analyses. Generally, material model parameterizations should be valid from temperatures of approximately -50 to 70° C, across a range of strain rates, and (depending on details of the parts involved) large deformations. Each material covered in this document is presented in its own chapter with a common format across materials. Model assumptions, limitations, existing validation results, readiness for use with uncertainty quantification, general usage guidance, and failure considerations are all provided along with specific parameterization inputs suitable for the finite element analysis code Sierra/Solid Mechanics.
In this paper we introduce a method to compare sets of full-field data using Alpert tree-wavelet transforms. The Alpert tree-wavelet methods transform the data into a spectral space allowing the comparison of all points in the fields by comparing spectral amplitudes. The methods are insensitive to translation, scale and discretization and can be applied to arbitrary geometries. This makes them especially well suited for comparison of field data sets coming from two different sources such as when comparing simulation field data to experimental field data. We have developed both global and local error metrics to quantify the error between two fields. We verify the methods on two-dimensional and three-dimensional discretizations of analytical functions. We then deploy the methods to compare full-field strain data from a simulation of elastomeric syntactic foam.
Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.
Time-dependent, viscoelastic responses of materials like polymers and glasses have long been studied. As such, a variety of models have been put forth to describe the behavior including simple rheological models (e.g. Maxwell, Kelvin), linear "fading memory" theories, and hereditary integral based linear thermal viscoelastic approaches as well as more recent nonlinear theories that are either integral, fictive temperature, or differential internal state variable based. The current work details a new LINEAR_THERMOVISCOELASTIC model that has been added to LAME. This formulation represents a viscoelastic theory that neglects some of the phenomenological details of the PEC/SPEC models in favor of efficiency and simplicity. Furthermore, this new model is a first step towards developing "modular" viscoelastic capabilities akin to those available with hardening descriptions for plasticity models in LAME. Specifically, multiple different (including user-defined) shift-factor forms are implemented with each being easily selected via parameter specification rather than requiring distinct material models.
In this work, we investigated microstructural features of elastomeric foam with the goal of identifying descriptors other than porosity that have a significant effect on the macroscale mechanical response. X-ray computed tomography (XCT) provided three-dimensional images of several flexible polyurethane foam samples prior to mechanical testing. The samples were then compressed to approximately 80% engineering strain. Stereo digital image correlation was used to measure the three-dimensional surface displacement data, from which strain was determined. The strain data, which were calculated with respect to the undeformed coordinates, were then overlaid on the corresponding surface generated from XCT. Heterogeneities in the strain-field were cross-correlated with topological quantities such as pore size distribution. A statistically significant correlation was identified between the distance transform of the pore phase and strain fluctuations.
In this work, we applied the in-situ X-ray Computed Tomography (XCT) mechanical testing method that coupled the in-situ mechanical loading with the XCT imaging to study the damage mechanism of GMBs inside the Sylgard as the material was subject to mechanical loading. We studied Sylgard specimens with different volume fraction of GMBs to understand how they behave differently under compression loading and how the volume fraction of GMBs affect the Sylgard failure. Both high resolution (1.5 μm/voxel) and low resolution (10 μm/voxel) XCT imaging were performed at different loading levels to visualize the GMB collapse during the compression of Sylgard with different volume fraction of GMBs. Feret shape of GMBs were calculated from the high resolution XCT images to determine whether the GMBs were intact or fractured, as well as the relationship between the size distribution of GMBs and their Feret shapes. Through these quantitative analysis of the high resolution XCT data, we were able to understand how the size and volume fraction of GMBs affect their failure behavior. The Digital volume correlation (DVC) technique was applied to the low resolution XCT images to calculate the local deformation of Sylgard specimen, which enabled us to understand the different failure propagation and failure mechanisms of Sylgard with different volume fraction of GMBs.
To elucidate the damage mechanisms in syntactic foams with hollow glass microballoon (GMB) reinforcement and elastomer matrices, in situ X-ray computed tomography mechanical testing was performed on syntactic foams with increasing GMB volume fraction. Image processing and digital volume correlation techniques identified very different damage mechanisms compared to syntactic foams with brittle matrices. In particular, the prevailing mechanism transitioned from dispersed GMB collapse at low volume fraction to clustered GMB collapse at high volume fraction. Moreover, damage initiated and propagated earlier in closely-packed GMBs for all specimens. Both of these trends were attributed to increased interaction between closely-packed GMBs. This was confirmed by statistical analysis of GMB damage, which identified a consistent, inverse relationship between the probability of survival and the local coordination number (Nneighbor) across all specimens.
Silicone elastomer filled with glass micro balloons (GMB) is an elastomeric syntactic foam used in electronics and component packaging for encapsulation, potting, stress-relief layer, and electrical insulation purposes. Under mechanical loading, the reinforcing phase, namely the GMBs embedded in the elastomer matrix, may break or delaminate, leading to internal damage and macroscale stiffness degradation, which can alter the material's protective capacity against mechanical shock and vibration. The degree of damage is controlled by the loading history, delamination, and failure behavior of the GMBs. We investigate the GMB failure behavior in this work wherein we present an indentation experiment to measure the force required to fail individual GMBs that are either embedded in the elastomer matrix or adhered to the surface of an elastomer layer. The indentation apparatus is augmented with an inverted optical microscope to enable in situ imaging of the GMB. Failure modes for the embedded or non-embedded GMBs are discussed based on the morphology of the broken GMBs and the measured failure forces. We also measure the adhesion energy between the glass balloon and the elastomer, based on which the possibility of delamination between the GMB and the surrounding elastomer matrix during the failure process is evaluated. Our results can facilitate the development of a failure criterion of GMBs which is necessary for establishing a physics-based constitutive model to describe the macroscopic damage mechanics of elastomeric syntactic foams.
In this report, we investigate how manufacturing conditions result in the warpage of moderate density PMDI polyurethane foam (12-50 lb/ft 3 ) when they are released from a mold. We have developed a multiphysics modeling framework to simulate the manufacturing process including resin injection, foaming and mold filling, gelation of the matrix, elevated cure, vitrification, cool down, and demolding. We have implemented this framework within the Sierra Mechanics Finite Element Code Suite. We couple Aria for flow, energy conservation, and foaming/curing kinetics with Adagio for the nonlinear viscoelastic solid response in a multi-staged simulation process flow. We calibrate a model for the PMDI-10S (10 lb/ft 3 free rise foam) through a suite of characterization data presented here to calibrate the solid cure behavior of the foam. The model is then used and compared to a benchmark experiment, the manufacturing and warpage over 1 year of a 10 cm by 10 cm by 2.5 cm foam "staple'. This component features both slender and thick regions that warp considerably differently over time. Qualitative agreement between the model and the experiment is achieved but quantitative accuracy is not. 2
The ability to relax a macroscopically applied stress is often associated with molecular mobility, or the possibility for a molecule to move outside the confines of its current position, within the material of which the stress is applied. Here, a viscoelastic constitutive analysis is used to investigate the counter-intuitive experimental observation of “mobility decrease with increased deformation through yield” [1] for a glass forming polymer during stress relaxation while under compressive and tensile loading conditions. The behavior of an epoxy thermoset is examined using an extensively validated, thermorheologically simple, material “clock” model, the Simplified Potential Energy Clock (SPEC) model.[2] This methodology allows for a comparison between the linear viscoelastic (LVE) limit and the true non-linear viscoelastic (NLVE) representation and enables exploration of a wide range of conditions that are not practical to investigate experimentally. The model predicts the behavior previously described as “mobility decrease with increased deformation” in the LVE limit and at low strain rates for NLVE. Only when loading rates are sufficient to decrease the material shift factor by multiple orders of magnitude is the anticipated deformation induced mobility or “mobility increase with increased deformation” observed. While the model has not been “trained” for these behaviors, it also predicts that the normalized stress relaxation response is indistinguishable amongst strain levels in the “post-yield” region, as has been experimentally reported. At long time, which has not been examined experimentally, the model predicts that even the normalized relaxation curves that exhibit “mobility increase with increased deformation” “cross back over” and return to the LVE ordering. These findings demonstrate the ability of rheologically simple models to represent the counter-intuitive experimentally measured material response and present predictions at long time scales that could be tested experimentally.
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.
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.
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.
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.
We are studying PMDI polyurethane with a fast catalyst, such that filling and polymerization occur simultaneously. The foam is over-packed to tw ice or more of its free rise density to reach the density of interest. Our approach is to co mbine model development closely with experiments to discover new physics, to parameterize models and to validate the models once they have been developed. The model must be able to repres ent the expansion, filling, curing, and final foam properties. PMDI is chemically blown foam, wh ere carbon dioxide is pr oduced via the reaction of water and isocyanate. The isocyanate also re acts with polyol in a competing reaction, which produces the polymer. A new kinetic model is developed and implemented, which follows a simplified mathematical formalism that decouple s these two reactions. The model predicts the polymerization reaction via condensation chemis try, where vitrification and glass transition temperature evolution must be included to correctly predict this quantity. The foam gas generation kinetics are determined by tracking the molar concentration of both water and carbon dioxide. Understanding the therma l history and loads on the foam due to exothermicity and oven heating is very important to the results, since the kinetics and ma terial properties are all very sensitive to temperature. The conservation eq uations, including the e quations of motion, an energy balance, and thr ee rate equations are solved via a stabilized finite element method. We assume generalized-Newtonian rheology that is dependent on the cure, gas fraction, and temperature. The conservation equations are comb ined with a level set method to determine the location of the free surface over time. Results from the model are compared to experimental flow visualization data and post-te st CT data for the density. Seve ral geometries are investigated including a mock encapsulation part, two configur ations of a mock stru ctural 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. Thoughts on m odel improvements are also discussed.
Full-field axial deformation within molten-salt batteries was measured using x-ray imaging with a sampling moiré technique. This method worked for in situ testing of the batteries because of the inherent grid pattern of the battery layers when imaged with x-rays. High-speed x-ray imaging acquired movies of the layer deformation during battery activation. Numerical validation of the technique, as implemented in this paper, was done using synthetic and numerically shifted images. Typical results of a battery are shown for one test. Ongoing work on validation and more test results are in progress.
We are developing computational models to help understand manufacturing processes, final properties and aging of structural foam, polyurethane PMDI. Th e resulting model predictions of density and cure gradients from the manufacturing process will be used as input to foam heat transfer and mechanical models. BKC 44306 PMDI-10 and BKC 44307 PMDI-18 are the most prevalent foams used in structural parts. Experiments needed to parameterize models of the reaction kinetics and the equations of motion during the foam blowing stages were described for BKC 44306 PMDI-10 in the first of this report series (Mondy et al. 2014). BKC 44307 PMDI-18 is a new foam that will be used to make relatively dense structural supports via over packing. It uses a different catalyst than those in the BKC 44306 family of foams; hence, we expect that the reaction kineti cs models must be modified. Here we detail the experiments needed to characteriz e the reaction kinetics of BKC 44307 PMDI-18 and suggest parameters for the model based on these experiments. In additi on, the second part of this report describes data taken to provide input to the preliminary nonlinear visco elastic structural response model developed for BKC 44306 PMDI-10 foam. We show that the standard cu re schedule used by KCP does not fully cure the material, and, upon temperature elevation above 150°C, oxidation or decomposition reactions occur that alter the composition of the foam. These findings suggest that achieving a fully cured foam part with this formulation may be not be possible through therma l curing. As such, visco elastic characterization procedures developed for curing thermosets can provide only approximate material properties, since the state of the material continuously evolves during tests.
Lithium-ion battery electrodes rely on a percolated network of solid particles and binder that must maintain a high electronic conductivity in order to function. Coupled mechanical and electrochemical simulations may be able to elucidate the mechanisms for capacity fade. We present a framework for coupled simulations of electrode mechanics that includes swelling, deformation, and stress generation driven by lithium intercalation. These simulations are performed at the mesoscale, which requires 3D reconstruction of the electrode microstructure from experimental imaging or particle size distributions. We present a novel approach for utilizing these complex reconstructions within a finite element code. A mechanical model that involves anisotropic swelling in response to lithium intercalation drives the deformation. Stresses arise from small-scale particle features and lithium concentration gradients. However, we demonstrate, for the first time, that the largest stresses arise from particle-to-particle contacts, making it important to accurately represent the electrode microstructure on the multi-particle scale. Including anisotropy in the swelling mechanics adds considerably more complexity to the stresses and can significantly enhance peak particle stresses. Shear forces arise at contacts due to the misorientation of the lattice structure. These simulations will be used to study mechanical degradation of the electrode structure through charge/discharge cycles.
Light activated polymers, which are capable of mechanically responding to light, promise to offer exciting, innovative, and unique material capabilities. Such materials include: photo-radical mediated cleavage and reformation of the polymer backbone in cross-linked elastomers that results in local stress relaxation; photo-switching cross-links in shape memory polymers; and photo-isomerization of azobenzene groups contained in liquid crystal elastomers. In this paper, using our recent material model that couples multiphysical processes involved in light-activated polymers, we demonstrate that a variety of patterns can be created on light activated polymer thin films when coupling mechanical deformation with light irradiation. Here, the polymer thin film is first stretched uniaxially or biaxially. Light is then irradiated on the surface of the thin film. After light irradiation, removal external load partially recovers the initial stretching of the polymer thin film and induces patterns. The variation of the geometry of the patterns can be controlled by a variety of parameters such as initial stretching, light intensity, etc. Photo-patterning with light activated polymer therefore offers a novel way to create surface patterns.