Search NASA⌕ Search

SEARCH · Search NASA

Results for “Deformation Gradient”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Simultaneous concentration and purification through gradient deformation chromatography

Mobile-phase additives, commonly used to modulate absorbate retention in gradient elution chromatography, are usually assumed to be either linearly retained or unretained. Previous theoretical work from our laboratory has shown that these modulators, such as salts in ion-exchange and hydrophobic interaction chromatography and organic modifiers in reversed-phase chromatography, can absorb nonlinearly, giving rise to gradient deformation. Consequently, adsorbate peaks that elute in the vicinity of the head of the deformed gradient may exhibit unusual shapes, form shoulders, and/or be concentrated. These effects for a reversed-phase sorbent with aqueous acetonitrile (ACN) as the modulator are verified experimentally. Gradient deformation is demonstrated experimentally and agrees with simulations based on ACN isotherm parameters that are independently determined from batch equilibrium studies using the layer model. Unusual absorbate peak shapes were found experimentally for single-component injections of phenylalanine, similar to those calculated by the simulations. A binary mixture of tryptophan and phenylalanine is used to demonstrate simultaneous concentration and separation, again in agreement with simulations. The possibility of gradient deformation in ion-exchange and hydrophobic interaction chromatography is discussed.

NASA Discipline Number 93-10↗
NASA NTRS2015

Representing Matrix Cracks Through Decomposition of the Deformation Gradient Tensor in Continuum Damage Mechanics Methods

A method is presented to represent the large-deformation kinematics of intraply matrix cracks and delaminations in continuum damage mechanics (CDM) constitutive material models. The method involves the additive decomposition of the deformation gradient tensor into 'crack' and 'bulk material' components. The response of the intact bulk material is represented by a reduced deformation gradient tensor, and the opening of an embedded cohesive interface is represented by a normalized cohesive displacement-jump vector. The rotation of the embedded interface is tracked as the material deforms and as the crack opens. The distribution of the total local deformation between the bulk material and the cohesive interface components is determined by minimizing the difference between the cohesive stress and the bulk material stress projected onto the cohesive interface. The improvements to the accuracy of CDM models that incorporate the presented method over existing approaches are demonstrated for a single element subjected to simple shear deformation and for a finite element model of a unidirectional open-hole tension specimen. The material model is implemented as a VUMAT user subroutine for the Abaqus/Explicit finite element software. The presented deformation gradient decomposition method reduces the artificial load transfer across matrix cracks subjected to large shearing deformations, and avoids the spurious secondary failure modes that often occur in analyses based on conventional progressive damage models.

Leone, Frank A., Jr.↗
NASA NTRS2016

A Continuum Damage Mechanics Model to Predict Kink-Band Propagation Using Deformation Gradient Tensor Decomposition

A new model is proposed that represents the kinematics of kink-band formation and propagation within the framework of a mesoscale continuum damage mechanics (CDM) model. The model uses the recently proposed deformation gradient decomposition approach to represent a kink band as a displacement jump via a cohesive interface that is embedded in an elastic bulk material. The model is capable of representing the combination of matrix failure in the frame of a misaligned fiber and instability due to shear nonlinearity. In contrast to conventional linear or bilinear strain softening laws used in most mesoscale CDM models for longitudinal compression, the constitutive response of the proposed model includes features predicted by detailed micromechanical models. These features include: 1) the rotational kinematics of the kink band, 2) an instability when the peak load is reached, and 3) a nonzero plateau stress under large strains.

Bergan, Andrew C.↗
DOE OSTI2021

A finite element formulation for deformation twinning induced strain localization in polycrystal magnesium alloys

Deformation twinning induces shear strain localization in hexagonal close-packed crystals and is critical for the material’s ductility and failure. Cracks often occur at twin-twin or twin-grain boundary intersections and propagate along twin bands. However, most crystal plasticity models for deformation twinning are based on a “pseudo-slip” approach and do not capture the localized deformation associated with the formation of each discrete twin band. The few exceptions are discrete twin models that involve very complex numerical algorithms and are often compromised in accuracy due to the numerical convergence. These factors make the discrete twin models hard to adopt. This paper proposes a modification to the conventional finite element weak form, to fully incorporate a twin-induced heterogeneous deformation that does not depend on the “pseudo-slip” assumption. The model starts by splitting the deformation gradient into elastic-slip-twinning components. The twin-induced deformation gradient component is computed separately by solving a microstructural evolution problem and then implemented into finite element weak form by constructing a global “twin-force” vector. The constitutive update (e.g., in the user-defined material subroutine, or UMAT, for ABAQUS) therefore avoids dealing with the twinning and recovers to the form of a regular slip-based crystal plasticity model. The results presented here indicate that the twin-induced strain localization and the associated stress-reversal phenomena near the twin band were naturally captured in the model, which was validated against an in-situ synchrotron X-ray micro-diffraction experiment.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Higher-Order Approximations for Stabilizing Zero-Energy Modes in Peridynamics Crystal Plasticity Models with Large Horizon Interactions

The non-ordinary state-based peridynamics theory combines non-local dynamic techniques with a desirable correspondence material principle, allowing for the use of continuum mechanics constitutive models. Such an approach presents a unique capability for solving problems involving discontinuities (e.g., strain localization, fracture, and fragmentation). However, the correspondence-based peridynamics models often suffer from zero-energy mode instabilities in numerical implementation, primarily due to the approximations of the non-local deformation gradient tensor. This paper focuses on a computational scheme for eliminating the zero-energy mode oscillations using a choice of influence functions that improve the truncation error in a higher-order Taylor series expansion of the deformation gradient. The novelty here is a tensor-based derivation of the linear constraint equations, which can be used to systematically identify the particle interaction weight functions for various user-specified horizon radii. In this paper, the proposed higher-order stabilization scheme is demonstrated for multi-dimensional examples involving polycrystalline and composite microstructures, along with comparisons against conventional finite element methods. The proposed stabilization scheme is shown to be highly effective in suppressing the spurious zero-energy mode oscillations in all numerical examples while enabling efficient simulations of strain localizations across material interfaces.

Non-Ordinary State-Based Peridynamics↗
DOE OSTI2023

Eulerian finite element implementations of a dislocation density-based continuum model

In Eulerian finite element simulations, the mesh moves relative to the material. After every change of position between the mesh and the material, the state variables are interpolated to the new mesh position, which is referred to as advection. Large strain crystal plasticity models are based on the multiplicative decomposition of the total deformation gradient. Further, the stress is evaluated as a function of the thermoelastic strain, temperature, and other state variables. Advection of tensor quantities, such as the strain, is coupled with possibly significant advection errors. In an effort to reduce the advection errors, we develop two rate forms of an established dislocation density-based continuum model. To that end, we replace the multiplicative decomposition of the deformation gradient with the additive decomposition of the velocity gradient, and define the stress rate instead of the total stress. The Eulerian implementation is compared with Lagrangian calculations, and two numerical examples with severe deformation levels are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Comparative analysis of plasticity-based GND density estimation methods in crystal plasticity finite element models

In crystal plasticity finite element (CPFE) simulations, accurately quantifying geometrically necessary dislocations (GNDs) is critical for capturing strain gradients in polycrystals. We compare different methods for quantifying GNDs, all of which originate from the Nye tensor, which is computed as the curl of the plastic deformation gradient. The projection technique directly decomposes the Nye tensor onto individual screw and edge dislocation components to compute GNDs. This approach requires converting a nine-component Nye tensor into densities for a larger number of dislocation systems, a fundamentally underdetermined (non-unique) process, which is resolved using L2 minimization. In contrast, when employing CPFE analysis, one could directly compute dislocation densities on each slip system using shear gradients. Projection and slip gradient methods are compared with respect to their prediction of GNDs with changing grain size, strain, and grain neighborhoods, including multigrain junctions. Although these techniques match analytical GND densities for single slip, single crystal deformation, and are consistent with anticipated overall GND trends, we find that the GND densities from projection techniques are significantly lower than those predicted from CPFE-based slip gradients in polycrystals. A suggested improvement of only using the active dislocation systems in the projection technique almost entirely resolved this mismatch.

Crystal plasticity↗
NASA NTRS1976

Dynamic response of membranes with finite elements

A general set of finite-element equations is developed for treating the nonlinear dynamic response of membranes consisting of complex nonlinear materials. The approach proposed by Key (1974) for a most complete treatment of the finite-element method to the nonlinear dynamic analysis of solids is used in developing the nonlinear equations of motion for membranes with cubic quadrilateral finite elements. The computations are kept general with respect to material modeling by computing the deformation gradient and requiring a constitutive subroutine to supply the Cauchy stress. An example calculation of the vibration of a square membrane demonstrates the accuracy of the formulation for small deformation theory and shows the nonlinear effect of large deformations.

Benzley, S. E.↗
NASA NTRS1982

The effect of small temperature gradients on flow in a continuous flow electrophoresis chamber

Continuous flow electrophoresis employs an electric field to separate biological cells suspended in a flowing liquid buffer solution. Good separations based on differences in electrophoretic mobility are obtained only when a unidirectional flow is maintained. The desired flow has a parabolic structure in the narrow dimension of the chamber and is uniform acros the width, except near the edges where the no-slip condition prevails. However, because of buoyancy, very small laterall or axial temperature gradients deform the flow significantly. The results of experiments conducted with a specially instrumented chamber show the origin and structure of the buoyancy-driven perturbations. It is found that very small temperature gradients can disturb the flow significantly, as was predicted by earlier theoretical work.

Rhodes, P. H.↗
NASA NTRS1982

On finite deformation elasto-plasticity

Lee (1969) has proposed a theory based on the decomposition of the total deformation gradient to an elastic and plastic part, and from it has concluded that the additive decomposition of the strain rates holds only approximately. Lubarda and Lee (1981) have declared that Lee's 'exact finite-deformation kinematics shows the almost universal assumption that the total velocity strain or rate of deformation is the sum of elastic and plastic rates to be in error'. Hence questions are raised regarding the validity of essentially all finite deformation elasto-plasticity theories. The present investigation is concerned with these questions. It is shown that the additive decomposition of the strain rates follows from all common finite elasto-plasticity concepts. Lee's theory is examined, and it is shown that this theory also leads to an additive strain rate decomposition, and therefore, his conclusion stems from misinterpretation. It is found that the elastic and the plastic strain rates considered by Lee do not correspond to the same configuration. They are, therefore, not compatible measures.

Nemat-Nasser, S.↗
NASA NTRS1995

Effect of modulator sorption on gradient shape in ion-exchange chromatography

Mobile phase additives, or modulators, are used in gradient elution chromatography to facilitate separation and reduce separation time. The modulators are usually assumed to be linearly adsorbed or unadsorbed. Here, the consequences of nonlinear modulator adsorption are examined for ion-exchange gradient elution through a series of simulations. Even when the buffer salt is identical to the modulator salt, gradient deformation is observed; the extent of deformation increases as the volume of the feed is increased. When the modulator salt is different from the buffer salt, unusual effects are observed, and the chromatograms are quite different from those predicted by classical gradient elution theory. In particular, local increases in the buffer concentration are found between feed bands, and serve to improve the separation. These effects become more pronounced as the feed volume increases, and could therefore prove valuable in preparative applications.

NASA Discipline Life Support Systems↗
NASA NTRSDate not supplied

Nonlinear Homogenization of Finitely Deformed Viscoelastic-Viscoplastic Composites Using Mechanics of Structure Genome

The objective of this paper is to develop a micromechanics approach to homogenizing finitely deformed viscoelastic-viscoplastic composites using the mechanics of structure genome. The incremental constitutive relation for glassy polymers, formulated in the spatial configuration, is implemented in the present approach.This involves (1) pulling-back the constitutive model to the material configuration and (2)choosing the deformation gradient tensor and the first Piola–Kirchhoff stress tensor as the strain and the stress measures during homogenization, respectively. An Euler–Newton predictor–corrector method is developed for homogenization. Each step involves formulating a variational statement using the mechanics of structure genome, discretizing the statement in a finite-dimensional space, and solving the problem using an Euler/multilevel Newton method. The present approach is demonstrated by homogenizing fiber- and particle-reinforced composites undergoing uniaxial, biaxial, or shear deformation, at different stain rates.

Multi-scale modeling, High Strain Composites, Visc↗
DOE OSTI2024

Sequential formulation of all‐way coupled finite strain thermoporomechanics for largely deformable gas hydrate deposits

We develop a numerically stable sequential formulation of thermoporomechanics for largely deformable gas hydrate deposits, extended from the fixed stress split of infinitesimal transformation. Constitutive equations are based on the total Lagrangian approach for both flow and geomechanics, including dynamic full tensor permeability and thermal conductivity updated from the deformation gradient. For space discretization, we take the cell-centered finite volume and node-based finite element method for flow and geomechanics, respectively. Then, we propose a sequential implicit method for all-way coupled thermoporomechanics, where the nonisothermal multiphase flow problem of gas hydrates is solved implicitly first and then the geomechanics problem is solved implicitly at the next step. During solution of the flow problem, we fix the rate of first Pioal total stress for numerical stability as well as apply porosity correction and entropy correction to account for geomechanical effects. We test numerical examples where flow and geomechanics parameters are based on deep oceanic gas hydrate deposits. When applying depressurization, even though the results between the infinitesimal transformation and finite strain geomechanics are similar in the early stages due to small deformation, we find differences between them in the late times as deformation becomes large. Accordingly, permeability and thermal conductivity tensors become nonisotropic full tensors although they are initially isotropic. Furthermore, we identify numerical stability of the developed sequential method from the test cases that exhibit the highly complex coupled gas hydrate systems with large deformation. Thus, the proposed sequential formulation can be applied in largely deformable gas hydrate systems.

42 ENGINEERING↗
DOE OSTI2022

A smeared crack modeling framework accommodating multi-directional fracture at finite strains

A generic smeared crack modeling framework predicated on the deformation gradient decomposition (DGD) approach is proposed for use in dynamic fracture problems at finite strains, accommodating failure along multiple mutually orthogonal fracture planes embedded within an independently defined bulk material model. Within this constitutive framework, the traction equilibrium conditions imposed at each failure surface are used to determine the associated crack displacements stored as internal state variables. In general, the enforcement of interfacial equilibrium entails the implicit solution of a non-linear system of equations within the constitutive update procedure. However, if inertial effects arising due to the relative motion of the fractured material are incorporated within the model, the traction equilibrium conditions are shown to give rise to corresponding dynamic equations of motion governing the time-evolution of the crack opening displacements. For dynamic problems, an explicit time-integration procedure is devised to efficiently update the material state, subject to a set of internal frictionless contact constraints to prevent material inversion. Finally, the efficacy of the proposed modeling framework is investigated through several benchmark dynamic fracture problems run within the explicit finite element code DYNA3D.

42 ENGINEERING↗
DOE OSTI2023

Automated extraction of interfacial dislocations and disconnections from atomistic data

Here we introduce Interfacial Line Defect Analysis (ILDA), a method for identifying and extracting interfacial dislocations and disconnections with little-to-no input from the user on the nature of the interface. By simply providing the orientations and coherency strains for the crystals in a coherent reference state of the interface, ILDA provides exact Burgers vectors. Alternatively, these orientations and strains can be estimated using local atomic deformation gradients, making ILDA fully automated and providing estimated Burgers vectors. ILDA also determines the step height associated with each defect line segment, in case the associated defect is a disconnection. The heart of the method is the identification of atoms residing at co-incidence sites between the two crystals and the construction of a surface mesh connecting these sites that is used to compose Burgers circuits and insert defect line segments. We demonstrate the performance of ILDA in two test cases: a twist grain boundary and a phase boundary.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A multiscale phase field fracture approach based on the non-affine microsphere model for rubber-like materials

Rubber-like materials have a broad scope of applications due to their unique properties like high stretchability and increased toughness. Hence, computational models for simulating their fracture behavior are paramount for designing them against failures. In this study, the phase field fracture approach is integrated with a multiscale polymer model for predicting the fracture behavior in elastomers. At the microscale, damaged polymer chains are modeled to be made up of a number of elastic chain segments pinned together. Using the phase field approach, the damage in the chains is represented using a continuous variable. Both the bond stretch internal energy and the entropic free energy of the chain are assumed to drive the damage, and the advantages of this assumption are expounded. A framework for utilizing the non-affine microsphere model for damaged systems is proposed here by considering the minimization of a hypothetical undamaged free energy, ultimately connecting the chain stretch to the macroscale deformation gradient. At the macroscale, a thermodynamically consistent formulation is derived in which the total dissipation is assumed to be mainly due to the rupture of molecular bonds. Using a monolithic scheme, the proposed model is numerically implemented and the resulting three-dimensional simulation predictions are compared with existing experimental data. The capability of the model to qualitatively predict the propagation of complex crack paths and quantitatively estimate the overall fracture behavior is verified. Additionally, the effect of the length scale parameter on the predicted fracture behavior is studied for an inhomogeneous system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Multiscale simulation of spatially correlated microstructure via a latent space representation

When deformation gradients act on the scale of the microstructure of a part due to geometry and loading, spatial correlations and finite-size effects in simulation cells cannot be neglected. We propose a multiscale method that accounts for these effects using a variational autoencoder to encode the structure–property map of the stochastic volume elements making up the statistical description of the part. In this paradigm the autoencoder can be used to directly encode the microstructure or, alternatively, its latent space can be sampled to provide likely realizations. Furthermore, we demonstrate the method on three examples using the common additively manufactured material AlSi10Mg in: (a) a comparison with direct numerical simulation of the part microstructure, (b) a push forward of microstructural uncertainty to performance quantities of interest, and (c) a simulation of functional gradation of a part with stochastic microstructure.

Elastoplasticity↗