Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Here, a new method for the solution of the non-linear equations forming the core of constitutive model integration is proposed. Specifically, the trust-region method that has been developed in the numerical optimization community is successfully modified for use in implicit integration of elastic-plastic models. Although attention here is restricted to these rate-independent formulations, the proposed approach holds substantial promise for adoption with models incorporating complex physics, multiple inelastic mechanisms, and/or multiphysics. As a first step, the non-quadratic Hosford yield surface is used as a representative case to investigate computationally challenging constitutive models. The theory and implementation are presented, discussed, and compared to other common integration schemes. Multiple boundary value problems are studied and used to verify the proposed algorithm and demonstrate the capabilities of this approach over more common methodologies. Robustness and speed are then investigated and compared to existing algorithms. Through these efforts, it is shown that the utilization of a trust-region approach leads to superior performance versus a traditional closest-point projection Newton-Raphson method and comparable speed and robustness to a line search augmented scheme.

Here, a new method for the solution of the non-linear equations forming the core of constitutive model integration is proposed. Specifically, the trust-region method that has been developed in the numerical optimization community is successfully modified for use in implicit integration of elastic-plastic models. Although attention here is restricted to these rate-independent formulations, the proposed approach holds substantial promise for adoption with models incorporating complex physics, multiple inelastic mechanisms, and/or multiphysics. As a first step, the non-quadratic Hosford yield surface is used as a representative case to investigate computationally challenging constitutive models. The theory and implementation are presented, discussed, and compared to other common integration schemes. Multiple boundary value problems are studied and used to verify the proposed algorithm and demonstrate the capabilities of this approach over more common methodologies. Robustness and speed are then investigated and compared to existing algorithms. Through these efforts, it is shown that the utilization of a trust-region approach leads to superior performance versus a traditional closest-point projection Newton-Raphson method and comparable speed and robustness to a line search augmented scheme.

A new method for the solution of the non-linear equations forming the core of constitutive model integration is proposed. Specifically, the trust-region method that has been developed in the numerical optimization community is successfully modified for use in implicit integration of elastic-plastic models. Although attention here is restricted to these rate-independent formulations, the proposed approach holds substantial promise for adoption with models incorporating complex physics, multiple inelastic mechanisms, and/or multiphysics. As a first step, the non-quadratic Hosford yield surface is used as a representative case to investigate computationally challenging constitutive models. The theory and implementation are presented, discussed, and compared to other common integration schemes. Multiple boundary value problems are studied and used to verify the proposed algorithm and demonstrate the capabilities of this approach over more common methodologies. Robustness and speed are then investigated and compared to existing algorithms. As a result through these efforts, it is shown that the utilization of a trust-region approach leads to superior performance versus a traditional closest-point projection Newton-Raphson method and comparable speed and robustness to a line search augmented scheme.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The numerical integration of constitutive models in computational solid mechanics codes allows for the solution of boundary value problems involving complex material behavior. Metal plasticity models, in particular, have been instrumental in the development of these codes. Here, most plasticity models implemented in computational codes use an isotropic von Mises yield surface. The von Mises, of J₂, yield surface has a simple predictor-corrector algorithm - the radial return algorithm - to integrate the model.

Abstract not provided.

Experiments were performed to characterize the mechanical response of several different rigid polyurethane foams to large deformation. In these experiments, the effects of load path, loading rate, and temperature were investigated. Results from these experiments indicated that rigid polyurethane foams exhibit significant damage, volumetric and deviatoric plasticity when they are compressed. Rigid polyurethane foams were also found to be extremely strain-rate and temperature dependent. These foams are also rather brittle and crack when loaded to small strains in tension or to larger strains in compression. Thus, a phenomenological Unified Creep Plasticity Damage (UCPD) model was developed to describe the mechanical response of these foams to large deformation at a variety of temperatures and strain rates. This paper includes a description of recent experiments and experimental findings. Next, development of a UCPD model for rigid, polyurethane foams is described. Finite element simulations with the new UCPD model are compared with experimental results to show behavior that can be captured with this model.

Abstract not provided.

Abstract not provided.

Numerous experiments were performed to characterize the mechanical response of several different rigid polyurethane foams (FR3712, PMDI10, PMDI20, and TufFoam35) to large deformation. In these experiments, the effects of load path, loading rate, and temperature were investigated. Results from these experiments indicated that rigid polyurethane foams exhibit significant volumetric and deviatoric plasticity when they are compressed. Rigid polyurethane foams were also found to be very strain-rate and temperature dependent. These foams are also rather brittle and crack when loaded to small strains in tension or to larger strains in compression. Thus, a new Unified Creep Plasticity Damage (UCPD) model was developed and implemented into SIERRA with the name Foam Damage to describe the mechanical response of these foams to large deformation at a variety of temperatures and strain rates. This report includes a description of recent experiments and experimental findings. Next, development of a UCPD model for rigid, polyurethane foams is described. Selection of material parameters for a variety of rigid polyurethane foams is then discussed and finite element simulations with the new UCPD model are compared with experimental results to show behavior that can be captured with this model.

Abstract not provided.

This memo derives a particular use of a rate-dependent, perfectly plastic Von Mises Plasticity constitutive model that represents a Newtonian (or Bingham) fluid. Only two parameters are needed: the viscosity and bulk modulus to complete the model parameterization. Under certain deformations, the model exactly represents a Newtonian fluid. Under others, the model shear thins. A Newtonian fluid representation may be used in confined flow applications in which a fluid-solid interaction is desired but code coupling is not.

Abstract not provided.

The Virtual Fields Method (VFM) is an inverse method for constitutive model parameter identication that relies on full-eld experimental measurements of displacements. VFM is an alternative to standard approaches that require several experiments of simple geometries to calibrate a constitutive model. VFM is one of several techniques that use full-eld exper- imental data, including Finite Element Method Updating (FEMU) techniques, but VFM is computationally fast, not requiring iterative FEM analyses. This report describes the im- plementation and evaluation of VFM primarily for nite-deformation plasticity constitutive models. VFM was successfully implemented in MATLAB and evaluated using simulated FEM data that included representative experimental noise found in the Digital Image Cor- relation (DIC) optical technique that provides full-eld displacement measurements. VFM was able to identify constitutive model parameters for the BCJ plasticity model even in the presence of simulated DIC noise, demonstrating VFM as a viable alternative inverse method. Further research is required before VFM can be adopted as a standard method for constitu- tive model parameter identication, but this study is a foundation for ongoing research at Sandia for improving constitutive model calibration.

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.

Abstract not provided.

Instrumented, fully coupled thermal-mechanical experiments were conducted to provide validation data for finite element simulations of failure in pressurized, high temperature systems. The design and implementation of the experimental methodology is described in another paper of this conference. Experimental coupling was accomplished on tubular 304L stainless steel specimens by mechanical loading imparted by internal pressurization and thermal loading by side radiant heating. Experimental parameters, including temperature and pressurization ramp rates, maximum temperature and pressure, phasing of the thermal and mechanical loading and specimen geometry details were studied. Experiments were conducted to increasing degrees of deformation, up to and including failure. Mechanical characterization experiments of the 304L stainless steel tube material was also completed for development of a thermal elastic-plastic material constitutive model used in the finite element simulations of the validation experiments. The material was characterized in tension at a strain rate of 0.001/s from room temperature to 800 C. The tensile behavior of the tube material was found to differ substantially from 304L bar stock material, with the plasticity characteristics and strain to failure differing at every test temperature.

Coupled thermal-mechanical experiments with well-defined, controlled boundary conditions were designed through an iterative process involving a team of experimentalists, material modelers and computational analysts. First the basic experimental premise was selected: an axisymmetric tubular specimen mechanically loaded by internal pressurization and thermally loaded asymmetrically by side radiant heating. Then several integrated experimental-analytical steps were taken to determine the experimental details. The boundary conditions were mostly thermally driven and were chosen so they could be modeled accurately; the experimental fixtures were designed to ensure that the boundary conditions were met. Preliminary, uncoupled analyses were used to size the specimen diameter, height and thickness with experimental consideration of maximum pressure loads and fixture design. Iterations of analyses and experiments were used to efficiently determine heating parameters including lamp and heating shroud design, set off distance between the lamps and shroud and between the shroud and specimen, obtainable ramp rates, and the number and spatial placement of thermocouples. The design process and the experimental implementation of the final coupled thermomechanical failure experiment design will be presented.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Constitutive models for computational solid mechanics codes are in LAME--the Library of Advanced Materials for Engineering. These models describe complex material behavior and are used in our finite deformation solid mechanics codes. To ensure the correct implementation of these models, regression tests have been created for constitutive models in LAME. A selection of these tests is documented here. Constitutive models are an important part of any solid mechanics code. If an analysis code is meant to provide accurate results, the constitutive models that describe the material behavior need to be implemented correctly. Ensuring the correct implementation of constitutive models is the goal of a testing procedure that is used with the Library of Advanced Materials for Engineering (LAME) (see [1] and [2]). A test suite for constitutive models can serve three purposes. First, the test problems provide the constitutive model developer a means to test the model implementation. This is an activity that is always done by any responsible constitutive model developer. Retaining the test problem in a repository where the problem can be run periodically is an excellent means of ensuring that the model continues to behave correctly. A second purpose of a test suite for constitutive models is that it gives application code developers confidence that the constitutive models work correctly. This is extremely important since any analyst that uses an application code for an engineering analysis will associate a constitutive model in LAME with the application code, not LAME. Therefore, ensuring the correct implementation of constitutive models is essential for application code teams. A third purpose of a constitutive model test suite is that it provides analysts with example problems that they can look at to understand the behavior of a specific model. Since the choice of a constitutive model, and the properties that are used in that model, have an enormous effect on the results of an analysis, providing problems that highlight the behavior of various constitutive models to the engineer can be of great benefit. LAME is currently implemented in the Sierra based solid mechanics codes Adagio [3] and Presto [4]. The constitutive models in LAME are available in both codes. Due to the nature of a transient dynamics code--e.g. Presto--it is difficult to test a constitutive model due to inertia effects that show up in the solution. Therefore the testing of constitutive models is primarily done in Adagio. All of the test problems detailed in this report are run in Adagio. It is the goal of the constitutive model test suite to provide a useful service for the constitutive model developer, application code developer and engineer that uses the application code. Due to the conflicting needs and tight time constraints on solid mechanics code development, no requirements exist for implementing test problems for constitutive models. Model developers are strongly encouraged to provide test problems and document those problems, but given the choice of having a model without a test problem or no model at all, certain requirements must be kept loose. A flexible code development environment, especially with regards to research and development in constitutive modeling, is essential to the success of such an environment. This report provides documentation of a number of tests for the constitutive models in LAME. Each section documents a separate test with a brief description of the model, the test problem and the results. This report is meant to be updated periodically as more test problems are created and put into the test suite.

The Library of Advanced Materials for Engineering (LAME) provides a common repository for constitutive models that can be used in computational solid mechanics codes. A number of models including both hypoelastic (rate) and hyperelastic (total strain) constitutive forms have been implemented in LAME. The structure and testing of LAME is described in Scherzinger and Hammerand ([3] and [4]). The purpose of the present report is to describe the material models which have already been implemented into LAME. The descriptions are designed to give useful information to both analysts and code developers. Thus far, 33 non-ITAR/non-CRADA protected material models have been incorporated. These include everything from the simple isotropic linear elastic models to a number of elastic-plastic models for metals to models for honeycomb, foams, potting epoxies and rubber. A complete description of each model is outside the scope of the current report. Rather, the aim here is to delineate the properties, state variables, functions, and methods for each model. However, a brief description of some of the constitutive details is provided for a number of the material models. Where appropriate, the SAND reports available for each model have been cited. Many models have state variable aliases for some or all of their state variables. These alias names can be used for outputting desired quantities. The state variable aliases available for results output have been listed in this report. However, not all models use these aliases. For those models, no state variable names are listed. Nevertheless, the number of state variables employed by each model is always given. Currently, there are four possible functions for a material model. This report lists which of these four methods are employed in each material model. As far as analysts are concerned, this information is included only for the awareness purposes. The analyst can take confidence in the fact that model has been properly implemented and the methods necessary for achieving accurate and efficient solutions have been incorporated. The most important method is the getStress function where the actual material model evaluation takes place. Obviously, all material models incorporate this function. The initialize function is included in most material models. The initialize function is called once at the beginning of an analysis and its primary purpose is to initialize the material state variables associated with the model. Many times, there is some information which can be set once per load step. For instance, we may have temperature dependent material properties in an analysis where temperature is prescribed. Instead of setting those parameters at each iteration in a time step, it is much more efficient to set them once per time step at the beginning of the step. These types of load step initializations are performed in the loadStepInit method. The final function used by many models is the pcElasticModuli method which changes the moduli that are to be used by the elastic preconditioner in Adagio. The moduli for the elastic preconditioner are set during the initialization of Adagio. Sometimes, better convergence can be achieved by changing these moduli for the elastic preconditioner. For instance, it typically helps to modify the preconditioner when the material model has temperature dependent moduli. For many material models, it is not necessary to change the values of the moduli that are set initially in the code. Hence, those models do not have pcElasticModuli functions. All four of these methods receive information from the matParams structure as described by Scherzinger and Hammerand.

Abstract not provided.

Constitutive modeling is an important aspect of computational solid mechanics. Sandia National Laboratories has always had a considerable effort in the development of constitutive models for complex material behavior. However, for this development to be of use the models need to be implemented in our solid mechanics application codes. In support of this important role, the Library of Advanced Materials for Engineering (LAME) has been developed in Engineering Sciences. The library allows for simple implementation of constitutive models by model developers and access to these models by application codes. The library is written in C++ and has a very simple object oriented programming structure. This report summarizes the current status of LAME.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

It has been shown that thermal energy imparted to a metallic substrate by laser heating induces a transient temperature gradient through the thickness of the sample. In favorable conditions of laser fluence and absorptivity, the resulting inhomogeneous thermal strain leads to a measurable permanent deflection. This project established parameters for laser micro forming of thin materials that are relevant to MESA generation weapon system components and confirmed methods for producing micrometer displacements with repeatable bend direction and magnitude. Precise micro forming vectors were realized through computational finite element analysis (FEA) of laser-induced transient heating that indicated the optimal combination of laser heat input relative to the material being heated and its thermal mass. Precise laser micro forming was demonstrated in two practical manufacturing operations of importance to the DOE complex: micrometer gap adjustments of precious metal alloy contacts and forming of meso scale cones.

Abstract not provided.

Honeycomb is a structure that consists of two-dimensional regular arrays of open cells. High-density aluminum honeycomb has been used in weapon assemblies to mitigate shock and protect payload because of its excellent crush properties. In order to use honeycomb efficiently and to certify the payload is protected by the honeycomb under various loading conditions, a validated honeycomb crush model is required and the mechanical properties of the honeycombs need to be fully characterized. Volume I of this report documents an experimental study of the crush behavior of high-density honeycombs. Two sets of honeycombs were included in this investigation: commercial grade for initial exploratory experiments, and weapon grade, which satisfied B61 specifications. This investigation also includes developing proper experimental methods for crush characterization, conducting discovery experiments to explore crush behaviors for model improvement, and identifying experimental and material uncertainties.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.