Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constitutive models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A Realistic Full-Scale 3D Modeling of Turning Using Coupled Smoothed Particle Hydrodynamics and Finite Element Method for Predicting Cutting Forces

Computational modelling is an effective technique for understanding the complex physics of machining. Large deformations, material separation, and high computational requirements are the key challenges faced while simulating machining. This work introduces a full-scale three-dimensional model of turning operations using a combined approach based on the Smoothed Particle Hydrodynamics (SPH) and Finite Element (FE) methods. By exploiting the advantages of each method, this approach leads to high-fidelity coupled SPH-FE machining models. Cutting forces and chip morphology are the primary results of interest. The machining models are validated with the results of turning experiments. Two-dimensional machining model underpredicts the cutting force and feed force by approximately 49% and 70%, respectively. Moreover, passive force cannot be predicted using the two-dimensional model. On the other hand, with the three-dimensional models developed in this manuscript, the difference between the total simulated force and experimentally measured force is ∼17%. The chip morphologies correlate with experiments in terms of the direction of the chip movement and the “long” continuous chips observed while turning Al 6061. This work expands the realm of machining simulations from two-dimensional orthogonal machining or sectional three-dimensional model to a full-scale realistic simulation. The encouraging simulation results show the potential to study more complex phenomena, such as machining stability and tool path modulation.

42 ENGINEERING↗
DOE OSTI2022

A crystal plasticity finite element method modeling of Zircaloy with hydride phases based on scanning electron micrographs

This work developed a microstructure-based finite element model to predict the stress state of alloys with second phase inclusions. Quantitative microstructural details extracted from scanning electron microscopy (SEM) images were used to generate heterogeneous microstructures including the morphology and spatial distribution of hydrides. Generation of digital microstructure was achieved through two steps of tessellations using software Neper and Matlab. The process is demonstrated using an example of Zircaloy material with secondary phases of hydrides dispersed within and stress-strain response of Zircaloy containing hydrides was predicted. The constitutive material model for Zircaloy in this study was based on crystal plasticity theory which considers the hexagonal close-packed (HCP) atomic structure of Zircaloy material. The hydrides were modeled as isotropic elasto-plastic material. A parametric study had been conducted to understand the effect of volume fraction, and lamellae thickness of the hydride phase on the mechanical properties of the overall material. Results can help designers to alter the manufacturing process to obtain the enhanced mechanical properties for components used in nuclear applications made by Zircaloy material.

Kulkarni, Shank S.↗
NASA NTRS2006

A Thermodynamically Consistent Damage Model for Advanced Composites

A continuum damage model for the prediction of damage onset and structural collapse of structures manufactured in fiber-reinforced plastic laminates is proposed. The principal damage mechanisms occurring in the longitudinal and transverse directions of a ply are represented by a damage tensor that is fixed in space. Crack closure under load reversal effects are taken into account using damage variables established as a function of the sign of the components of the stress tensor. Damage activation functions based on the LaRC04 failure criteria are used to predict the different damage mechanisms occurring at the ply level. The constitutive damage model is implemented in a finite element code. The objectivity of the numerical model is assured by regularizing the dissipated energy at a material point using Bazant's Crack Band Model. To verify the accuracy of the approach, analyses of coupon specimens were performed, and the numerical predictions were compared with experimental data.

Maimi, Pere↗
DOE OSTI2020

An Integrated Multiscale Experimental-Numerical Analysis on Reconsolidation of Salt-Clay Mixture for Disposal of Heat-Generating Waste (Final NEUP Technical Report)

The overall purpose of this research is to improve understanding of THMC coupling effect on the reconsolidation of granular (or crushed) salt-clay mixture used for seal systems of shafts and drifts in salt repositories. This proposed work is partially motivated by the recent work on the Waste Isolation Pilot Plant (WIPP) that shows the promising sealing capability of clay-salt mixture compared to crushed salt. In particular, primary emphasis is to develop a fully integrated multiscale experiment-numerical study to determine and explain what leads to the superior sealing ability of the clay-salt mixture. These research activities are designed to seek further understanding of (1) why clay additives may enhance the fluid trapping and (2) whether this flow barrier effect may prevail under different combinations of temperature, confining pressure, deviatoric stress and other foreseeable environmental factors. If successful, this enhanced flow trapping ability of the seal provides significant improvement to the seal and repository performance and therefore make the repository safer in the long-term. The experiment component includes microstructural investigation and macroscopic tests on a reconsolidated salt-clay mixture. In the microstructural study, the goal is to (1) characterize microscopic distributions of distinct phases (e.g., clay, salt crystal boundaries, trapped brine, and pore) to examine the connectivity of the pore network inside the salt-clay mixture with different amounts of clay additive and moisture content and (2) analyze multiscale imaging data to reconstruct the polycrystalline microstructures for numerical simulations. Meanwhile, macroscopic tests are performed to analyze how clay alters the failure/creep mechanisms in the salt-clay mixture. Microscopic and macroscopic experimental observations will both be used to calibrate and validate a multiscale model that explicitly simulates the capillary and multiphase flow in the connected pores and the deformation due to the presence of intra-crystalline brine at the pore scale via a new polyhedral discrete element–lattice Boltzmann method (DEM-LBM) coupling model. The pore-scale simulations are homogenized via an upscaling procedure that converts pore-scale information (e.g. force exerted on grain boundary, sliding, pressure-solution) to continuum measures (e.g. Cauchy stress, Darcy’s flow) at each integration point in the macroscopic multiphase TMHC model. This multiscale scheme will allow coupling be- tween high-fidelity simulations of brine-salt-clay interaction and the macroscopic TMHC model. The multiscale model helps the understanding of how the trapped brine inclusion affects the pressure-solution mechanism with the presence of clay and moisture. This work brings new insight into the sealing capacity of salt-clay mixture under elevated temperature over a long period of time - a key to evaluating the potential of salt-clay mixture usage for salt repositories.

42 ENGINEERING↗
NASA NTRS2006

Database of Properties of Meteors

A database of properties of meteors, and software that provides access to the database, are being developed as a contribution to continuing efforts to model the characteristics of meteors with increasing accuracy. Such modeling is necessary for evaluation of the risk of penetration of spacecraft by meteors. For each meteor in the database, the record will include an identification, date and time, radiant properties, ballistic coefficient, radar cross section, size, density, and orbital elements. The property of primary interest in the present case is density, and one of the primary goals in this case is to derive densities of meteors from their atmospheric decelerations. The database and software are expected to be valid anywhere in the solar system. The database will incorporate new data plus results of meteoroid analyses that, heretofore, have not been readily available to the aerospace community. Taken together, the database and software constitute a model that is expected to provide improved estimates of densities and to result in improved risk analyses for interplanetary spacecraft. It is planned to distribute the database and software on a compact disk.

Suggs, Rob↗
DOE OSTI2025

Modeling of melting and erosion in high velocity microparticle impacts

High velocity collisions can cause melting, a situation involving elasticity, plasticity, and phase transformations in tandem, which is very challenging for mechanical modeling. Here, this paper advances the use of smoothed particle hydrodynamics (SPH) simulations for this complex situation, investigating stainless steel microspheres impacting a tin substrate. After extensive validation against a range of mechanical and geometrical measurements of experimental impacts, the model provides a detailed, quantitative view of melt volume production over a broad range of impact velocities. It also enables direct decomposition of the impact energy partitioning, including analysis of the energy absorbed by plasticity in the impactor, the system kinetic energy, and the adiabatic heating of the substrate in addition to the melting itself. Since these quantities are not available experimentally, SPH modeling substantially enhances the quantitative understanding of such impacts. Finally, the model speaks to the phase characteristics of the ejecta formed during impact, confirming the experimental inference that the ejected material likely contains a significant fraction of solid material.

Constitutive modeling↗
DOE OSTI2021

USAMP Low-Cost Mg Sheet Component Development and Demonstration Project

The overarching objective of this USAMP project was to develop and demonstrate door panels made from magnesium (Mg) sheet with a cost penalty over conventional steel stampings of no more than $\$5.50$/kg saved. The technical approach integrated experiments with advanced computational tools based on Integrated Computational Materials Engineering (ICME) methods to develop new alloy chemistries and their thermomechanical processing that promise improved formability and lower forming temperatures. A penultimate task before finally forming the stampings was to incorporate actual microstructure into models that would enable formability simulations. This approach would, for the first time, account for individual magnesium grains moving in an anisotropic fashion unlike that for aluminum or steel that have isotropic properties upon which the current simulation tools are based. In separate activities, new coatings and lubricants to facilitate forming and improved corrosion protection and joining strategies, were developed to ensure that the door could be produced with stated product requirements. A technical cost model, which included parts production, assembly, and paint for a door specifically designed for Mg sheet, showed the cost penalty to be between $\$4.26$ to $\$6.60$/kg saved, which enveloped the project’s cost targets. The cost of the coated Mg sheet was identified as the key driver for the cost penalty. The mass of the Mg-intensive door was 7.9 kg, which was 54% less than the baseline steel door.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Low-Cost Magnesium Alloy Sheet Component Development and Demonstration Project

The overarching objective of this USAMP project was to develop and demonstrate door panels made from magnesium (Mg) sheet with a cost penalty over conventional steel stampings of no more than $5.50/kg saved. The technical approach integrated experiments with advanced computational tools based on Integrated Computational Materials Engineering (ICME) methods to develop new alloy chemistries and their thermomechanical processing that promise improved formability and lower forming temperatures. A penultimate task before finally forming the stampings was to incorporate actual microstructure into models that would enable formability simulations. This approach would, for the first time, account for individual magnesium grains moving in an anisotropic fashion unlike that for aluminum or steel that have isotropic properties upon which the current simulation tools are based. In separate activities, new coatings and lubricants to facilitate forming and improved corrosion protection and joining strategies, were developed to ensure that the door could be produced with stated product requirements. A technical cost model, which included parts production, assembly, and paint for a door specifically designed for Mg sheet, showed the cost penalty to be between 4.26 USD to 6.60 USD/kg saved, which enveloped the project’s cost targets. The cost of the coated Mg sheet was identified as the key driver for the cost penalty. The mass of the Mg-intensive door was 7.9 kg, which was 54% less than the baseline steel door.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Model for the Softening Factor within Stages of Work Hardening

A formulation for a dimensionless coefficient c b is derived that represents a scale of microstructural softening for alloys which follow Kocks-Mecking (K-M) work hardening behavior. The variation of the true plastic strain ε p between the proportional limit σ y and the strength σ u at the instability is determined using the Considère criterion. Parameterization of the model is limited to variables expressly measured within tensile experiments. Further development is now made for the softening factor c bi through the individual and sequential stages 3 and 4 of plastic deformation during K-M work hardening Q behavior. Application is shown for tensile test results of Ti-6Al-4V made by different additively manufactured (AM) processes. It is found that the variation in plastic strain as a function of c bi produces a continuous curve representative of the alloy system. The results of data analysis indicate that Θ o3 and c b3 increase, while Θ o4 and c b4 decrease, as the total plastic strain ε p increases. Furthermore, formulations derived for K-M stages 3 and 4 enable the evaluation of other material parameters such as the activation volume ν* for the onset of plastic deformation. This activation volume is found to be near constant, at a ν*-value of 0.353±0.036 nm 3 as computed using a strain-rate sensitivity of strength exponent m of 0.014, irrespective of the AM method used to produce the Ti-6Al-4V alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Non-linear elastic-plastic behavior of the invert glass lithium phosphorous oxynitride (LiPON)

Micro and nano-scale mechanical behavior of network/ionic glasses is dictated by their composition and the number of constraints per structural unit. From this perspective, the glass LiPON presents the opportunity for enhancement of the microscale ductility due to its rather unconstrained orthophosphate structure. Here, we use instrumented nanoindentation with different tip geometries to investigate the mechanical response of LiPON glass. The results reveal that the elastic modulus of LiPON is not constant and depends on pressure. With the method utilizing spherical nanoindentation and continuous stiffness measurement (CSM) we determine the yield point of LiPON. We propose the Drucker-Prager type of yield criterion for LiPON and estimate its yield stress in compression as 2.4 GPa. There exists another stress threshold, however, around 490 MPa, at which the elastic deformation becomes non-linear, and this can be mistaken for the yield stress.

Constitutive model↗
DOE OSTI2020

Discrepancy in the Critical State Void Ratio of Poorly Graded Sand due to Shear Strain Localization

The critical state (CS) concept is a theoretical framework that models the constitutive behavior of soils, including sand and other granular materials. It supports the notion of a unique postfailure state, where the soil ultimately experiences continuous shearing with no change in the plastic volumetric strain. However, the published literature has frequently noted the nonconvergence of sand specimens with different initial densities to a unique CS in the compression plane due to many factors such as specimen fabric, particle morphology, breakage, and grain size distribution. This paper examines the CS for poorly graded (uniform) glass beads and 3 different types of silica sands using 50 conventional triaxial compression (CTC) experiments, 12 oedometer tests, and in situ synchrotron microcomputed tomography (SMT) scans for 10 CTC experiments. The results of the 50 CTC experiments revealed a diffused CS zone in the compression plane, which was further examined using the in situ SMT scans. A thorough three-dimensional image analysis of the SMT scans accurately quantified the evolution of the local void ratio (elocal) versus axial compression within zones of intensive shearing toward the center of the specimen. The evolution of the void ratio was also measured using the entire volume of the specimen (eglobal). At the CS, the elocal/eglobal ratio was assessed to be ~1.25 when a single shear band developed within the scanned specimens and ~1.1-1.15 for specimens that failed via external bulging that was internally manifested by the development of multiple shear bands. This finding suggests that the CS zone in the compression plane can be attributed to the common wrong consideration of eglobal evolution in lieu of elocal within the developing shear bands. Furthermore, the lack of shear band development in uniaxial compression has made the results of the oedometer test reliable in quantifying the CS parameters in the compression plane.

58 GEOSCIENCES↗
DOE OSTI2021

Genomic Signatures of a Major Adaptive Event in the Pathogenic Fungus Melampsora larici-populina

The recent availability of genome-wide sequencing techniques has allowed systematic screening for molecular signatures of adaptation, including in nonmodel organisms. Host–pathogen interactions constitute good models due to the strong selective pressures that they entail. We focused on an adaptive event which affected the poplar rust fungus Melampsora larici-populina when it overcame a resistance gene borne by its host, cultivated poplar. Based on 76 virulent and avirulent isolates framing narrowly the estimated date of the adaptive event, we examined the molecular signatures of selection. Using an array of genome scan methods based on different features of nucleotide diversity, we detected a single locus exhibiting a consistent pattern suggestive of a selective sweep in virulent individuals (excess of differentiation between virulent and avirulent samples, linkage disequilibrium, genotype–phenotype statistical association, and long-range haplotypes). Our study pinpoints a single gene and further a single amino acid replacement which may have allowed the adaptive event. Although our samples are nearly contemporary to the selective sweep, it does not seem to have affected genome diversity further than the immediate vicinity of the causal locus, which can be explained by a soft selective sweep (where selection acts on standing variation) and by the impact of recombination in mitigating the impact of selection. Therefore, it seems that properties of the life cycle of M. larici-populina, which entails both high genetic diversity and outbreeding, has facilitated its adaptation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

pnnl/MCRASTA

McRasta (Markov Chain Rate and State Analysis) was developed to estimate parameter uncertainty in constitutive friction models via Bayesian inverse and Markov Chain Monte Carlo (MCMC) methods.

Fichera, Marissa [Pacific Northwest National Labor↗
DOE OSTI2021

A Constitutive Structural Parameter for the Work Hardening Behavior of Additively Manufactured Ti-6Al-4V

The mechanical behavior of Ti-6Al-4V produced by additive manufacturing processes is assessed based on a formulation developed from the Kocks-Mecking relationship. A constitutive parameter derived for the microstructure is characteristic of the work hardening behavior determined by the plastic strain between the proportional limit and the strength at the instability point. The varied plasticity behavior associated with surface and build direction effects can be evaluated with this approach as presented for the case of Ti-6Al-4V under uniaxial tension.

36 MATERIALS SCIENCE↗
NASA NTRS1976

Nonlinear coupling between pulsation and convection in late type stars

A simple idealized nonlinear model applicable to long-period variable stars has been formulated which assumes that the convective envelope of M giants is composed of giant convection cells comparable in size to the stellar radius. The simplicity of this model essentially constitutes a physical analog to the strong dynamic coupling that occurs if the convective envelope of the star undergoes both modes of motion. As shown implicitly in the time scales associated with these motions, the coupling produces asymmetrical fluctuations of the entire star, the mean velocity of which is comparable to the escape velocity of the star at particular values of the ratio of the pulsation and convection time scales. It is suggested that this can account for the mass loss from late-type stars and the circumstellar dust shells that are associated extensively with long-period variables. For critical values of the pulsation and convection time scales, the solutions correspond to the rapid expansion of the entire convective envelope and are the basis of a mechanism that simulates the manner in which pulsating stars ballistically accelerate their convective shells to form planetary nebulae.

Anand, S. P. S.↗
NASA NTRS1983

Engine cyclic durability by analysis and material testing

The problem of calculating turbine engine component durability is addressed. Nonlinear, finite-element structural analyses, cyclic constitutive behavior models, and an advanced creep-fatigue life prediction method called strainrange partitioning were assessed for their applicability to the solution of durability problems in hot-section components of gas turbine engines. Three different component or subcomponent geometries are examined: a stress concentration in a turbine disk; a louver lip of a half-scale combustor liner; and a squealer tip of a first-stage high-pressure turbine blade. Cyclic structural analyses were performed for all three problems. The computed strain-temperature histories at the critical locations of the combustor linear and turbine blade components were imposed on smooth specimens in uniaxial, strain-controlled, thermomechanical fatigue tests of evaluate the structural and life analysis methods.

Kaufman, A.↗
NASA NTRS1982

Predicting viscoelastic response and delayed failures in general laminated composites

Although graphite fibers behave in an essentially elastic manner, the polymeric matrix of graphite/epoxy composites is a viscoelastic material which exhibits creep and delayed failures. The creep process is quite slow at room temperature, but may be accelerated by higher temperatures, moisture absorption, and other factors. Techniques are being studied to predict long-term behavior of general laminates based on short-term observations of the unidirectional material at elevated temperatures. A preliminary numerical procedure based on lamination theory is developed for predicting creep and delayed failures in laminated composites. A modification of the Findley nonlinear power law is used to model the constitutive behavior of a lamina. An adaptation of the Tsai-Hill failure criterion is used to predict the time-dependent strength of a lamina. Predicted creep and delayed failure results are compared with typical experimental data.

Dillard, D. A.↗