Search NASA⌕ Search

SEARCH · Search NASA

Results for “anisotropic elasticity”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Crystal plasticity modeling of strain-induced martensitic transformations to predict strain rate and temperature sensitive behavior of 304 L steels: Applications to tension, compression, torsion, and impact

This paper advances crystallographically-based Olson-Cohen (direct γ → α’) and deformation mechanism (indirect γ→ε→α’) phase transformation models for predicting strain-induced austenite to martensite transformation. Here, the advanced transformation models enable predictions of not only strain-path sensitive, but also of strain-rate and temperature sensitive deformation of polycrystalline stainless steels (SSs). The deformation of constituent grains in SSs is modeled as a combination of anisotropic elasticity, crystallographic slip, and phase transformation, while the hardening is based on the evolution of dislocation density and explicit shifts in phase fractions. Such grain-scale deformation is implemented within the meso-scale elasto-plastic self-consistent (EPSC) homogenization model, which is coupled with the implicit finite element (FE) method to provide a constitutive response at each FE integration point for solving boundary value problems at the macro-scale. Parameters pertaining to the hardening and transformation models within FEEPSC are calibrated and validated on a suite of data including flow curves and phase fractions for monotonic compression, tension, and torsion as a function of strain-rate and temperature for wrought and additively manufactured (AM) SS304L. To illustrate the potential and accuracy of the integrated multi-level FE-EPSC simulation framework, geometry, mechanical response, phase fractions, and texture evolution are simulated during gas-gun impact deformation of a cylinder and quasi-static tension of a notched specimen made of AM SS304L. Details of the simulation framework, comparison between experimental and simulation results, and insights from the results are presented and discussed.

304L steels↗
DOE OSTI2026

Cross-slip and easy-glide CRSS in titanium: Theoretical predictions and in-situ TEM measurements

This study investigates the mechanics of prismatic and first-order pyramidal $\langle$a$\rangle$ slip in titanium (Ti), elucidating the physics of easy-glide and cross-slip through a combination of theory and experiments. Screw-character prismatic (Pr) dislocations in Ti are of particular interest because their complex cores can be stable or unstable, leading to activation by either cross-slip or planar glide. To investigate these mechanisms, site-specific micro tensile samples were prepared using focused ion beam (FIB) milling and mounted on a push-to-pull (PTP) device for in-situ transmission electron microscopy (TEM) tensile testing. The in-situ experiments provide direct observations of the onset of dislocation motion and the precise determination of the critical resolved shear stress (CRSS) for the activated mechanisms, and its evolution with load cycling. A comprehensive theory has been developed to predict the CRSS values for easy glide, cross-slip, and multiplication of dislocations. Predicted critical stresses for pyramidal (π)-to-Pr and reverse cross-slip agree closely with the experimental measurements. The latter cross-slip stress is a factor of two higher than that of unobstructed planar slip. The model accounts for overlapping dislocation cores and employs a Wigner-Seitz based cell to evaluate misfit energies. By combining ab initio density functional theory (DFT) with anisotropic elasticity, the framework identifies minimum energy pathways for dislocation glide, which can be intermittent and zig-zag. A simplified expression utilizing (π) and Pr Schmid factor ratios is proposed for critical stress corresponding to (π)-to-Pr cross-slip transition. The results are strongly dependent on crystal orientation, underscoring non-Schmid behavior. Overall, this study explores key critical stress parameters essential for informing higher-scale simulations of plasticity in Ti.

Critical stress↗
DOE OSTI2021

Investigation of δ zirconium hydride morphology in a single crystal using quantitative phase field simulations supported by experiments

In light water nuclear reactors, waterside corrosion of the cladding material leads to the production of hydrogen, a fraction of which is picked up by the zirconium cladding and precipitates into brittle hydride particles. These nanoscale hydride particles aggregate into mesoscale hydride clusters. The principal stacking direction of the nanoscale hydrides precipitated in the cladding tube changes from circumferential in the absence of applied stress to radial under circumferential applied stress. A quantitative phase field model has been developed to predict the hydride morphology observed experimentally and identify the mechanisms responsible for nanoscale hydride stacking. The model focuses on nanoscale hydride precipitation in a single zirconium grain with a detailed description of the anisotropic elastic contribution. The model predictions concerning the shape, orientation, and stacking behavior of nanoscale hydride are analyzed and compared with experimental observations. The model accurately accounts for the experimentally observed elongated nanoscale hydride shape and the stacking of hydrides along the basal plane of the hexagonal zirconium matrix. When investigating the role of applied stress in hydride morphology, the model challenges some of the mechanisms previously proposed to explain hydride reorientation. Finally, although hydride reorientation has been hypothesized to be caused by a change in nanoscale hydride shape, the current study shows that these mechanisms are unlikely to occur.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Sea Ice Rheology Experiment (SIREx): 1. Scaling and Statistical Properties of Sea-Ice Deformation Fields

As the sea-ice modeling community is shifting to advanced numerical frameworks, developing new sea-ice rheologies, and increasing model spatial resolution, ubiquitous deformation features in the Arctic sea ice are now being resolved by sea-ice models. Initiated at the Forum for Arctic Modeling and Observational Synthesis, the Sea Ice Rheology Experiment (SIREx) aims at evaluating state-of-the-art sea-ice models using existing and new metrics to understand how the simulated deformation fields are affected by different representations of sea-ice physics (rheology) and by model configuration. Part 1 of the SIREx analysis is concerned with evaluation of the statistical distribution and scaling properties of sea-ice deformation fields from 35 different simulations against those from the RADARSAT Geophysical Processor System (RGPS). For the first time, the viscous-plastic (and the elastic-viscous-plastic variant), elastic-anisotropic-plastic, and Maxwell-elasto-brittle rheologies are compared in a single study. We find that both plastic and brittle sea-ice rheologies have the potential to reproduce the observed RGPS deformation statistics, including multi-fractality. Model configuration (e.g., numerical convergence, atmospheric representation, spatial resolution) and physical parameterizations (e.g., ice strength parameters and ice thickness distribution) both have effects as important as the choice of sea-ice rheology on the deformation statistics. It is therefore not straightforward to attribute model performance to a specific rheological framework using current deformation metrics. In light of these results, we further evaluate the statistical properties of simulated Linear Kinematic Features in a SIREx Part 2 companion paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Machine learning unifies flexibility and efficiency of spinodal structure generation for stochastic biomaterial design

Abstract Porous biomaterials design for bone repair is still largely limited to regular structures (e.g. rod-based lattices), due to their easy parameterization and high controllability. The capability of designing stochastic structure can redefine the boundary of our explorable structure–property space for synthesizing next-generation biomaterials. We hereby propose a convolutional neural network (CNN) approach for efficient generation and design of spinodal structure—an intriguing structure with stochastic yet interconnected, smooth, and constant pore channel conducive to bio-transport. Our CNN-based approach simultaneously possesses the tremendous flexibility of physics-based model in generating various spinodal structures (e.g. periodic, anisotropic, gradient, and arbitrarily large ones) and comparable computational efficiency to mathematical approximation model. We thus successfully design spinodal bone structures with target anisotropic elasticity via high-throughput screening, and directly generate large spinodal orthopedic implants with desired gradient porosity. This work significantly advances stochastic biomaterials development by offering an optimal solution to spinodal structure generation and design.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

A new woven composite constitutive model validated by shock wave experiments

In this paper, we present results of plate impact simulations of shock compressed woven glass fiber-reinforced plastic (GRP) performed using the Arbitrary Lagrangian–Eulerian three-dimensional finite element code. A hyperelastic large-strain-based empirical Continuum Damage Mechanics (CDM) formulation is employed to describe damage initiation and growth in the shock-compressed GRP. The model parameters calibration scheme utilizes the Velocity Interferometer System for Any Reflector normal particle velocity measurements at the free surface of the GRP target plates. The impact velocity in the experiments ranged from 8.5 to 418 m/s. The finite element model considered planar 0°/90° bidirectional plies with an individual ply thickness of 0.68 mm, stacked to reach a total laminate thickness of 6.8 mm. The anisotropic elastic strains were estimated from the experimentally determined tetragonal symmetry stiffness matrix for the GRP. The strain-based damage model captures several salient features observed in the measured free surface particle wave profiles, including the shock rise time, onset of Elastic—Elastic Cracking, and the shape of the nonlinear portion of the experimental particle velocity profiles. Furthermore, the CDM model predicts the dominant damage mode to be matrix microcracking due to shear and the associated bulk expansion (bulking) under the global compressive loading in the plate impact configuration.

42 ENGINEERING↗
DOE OSTI2025

Shear strength of a refractory high entropy alloy MoNbTaVW under high pressure

Radial X-ray diffraction (R-XRD) was performed in situ using a Panoramic Diamond Anvil Cell on the refractory high entropy alloy MoNbTaVW. The lattice parameters were determined through a Le Bail fit using a body centered cubic (BCC) lattice of symmetry $Im$$\overline{3}$$m$ (international space group number 229). Upper and lower bounds to the shear strength were determined up to 80 GPa nonhydrostatically using copper as a pressure standard. The equation of state was derived at the magic angle ψ = 54.7° and yielded a bulk modulus of K 0 = 220.8 ± 1.65 GPa. The experimental lattice parameters and bulk modulus match closely with corresponding density functional theory (DFT) calculations. The BCC phase remains stable up to the highest pressure of 80 GPa studied and is shown to be elastically anisotropic. The shear strength was found to saturate around 70 GPa with a value of τ = 1.75 GPa, and the shear moduli are presented in different limits of iso-strain and iso-stress.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Kinematics and dynamics of disclination lines in three-dimensional nematics

An exact kinematic law for the motion of disclination lines in nematic liquid crystals as a function of the tensor order parameter Q is derived. Unlike other order parameter fields that become singular at their respective defect cores, the tensor order parameter remains regular. Following earlier experimental and theoretical work, the disclination core is defined to be the line where the uniaxial and biaxial order parameters are equal, or equivalently, where the two largest eigenvalues of Q cross. This allows an exact expression relating the velocity of the line to spatial and temporal derivatives of Q on the line, to be specified by a dynamical model for the evolution of the nematic. By introducing a linear core approximation for Q, analytical results are given for several prototypical configurations, including line interactions and motion, loop annihilation, and the response to external fields and shear flows. Behaviour that follows from topological constraints or defect geometry is highlighted. Finally, the analytic results are shown to be in agreement with three-dimensional numerical calculations based on a singular Maier–Saupe free energy that allows for anisotropic elasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Highly Extensible X-ray Diffraction Toolkit

This code has been previous released as open source under LLNL-CODE-529294. The HEXRD program comprises a library of X-ray diffraction analysis tools and a stand-alone GUI application. While there are many software packages that provide similar capabilities, HEXRD was designed to provide an extremely flexible and extensible basis for describing a broad array of X-ray diffraction instruments; this includes both mono- and poly-chromatic modalities as well as poly- and single-crystal samples. The high-level objects include Material: an abstraction for defining crystalline materials, including unit cell parameters and anisotropic elastic moduli Imageseries: an interface for series of diffraction images containing one or more frames, including provisions for metadata. Instrument: an abstraction of an X-ray diffraction instrument including support for the definition of multiple detector elements as well as powder, Laue, and mono-chromatic single/multi-grain diffraction analysis modalities (including attendant calibration routines). The hexrdgui package is a standalone PySide2/Qt5 GUI built on top of hexrd. The 0.8.1 release includes various bug fixes as well as some new functionality in the GUI, including: o Enhancements to the import tools for LLNL diffraction instruemtns (PXRDIP and TARDIS) o The ability to save 1-d LeBail fits of integrated diffraction spectra.

Barton, NathanR.↗
DOE OSTI2021

Efficient Failure Probability Calculations and Modeling Interface Debonding in TRISO Particles with BISON

This document constitutes completion of the NEAMS milestone, which is titled Demonstrate efficiency improvements in TRISO failure probability calculations and the effect of debonding on fission product transport in TRISO fuel particles in Bison . In this report we present the development of the following: (1) the anisotropic elasticity model of pyrolytic carbon and extended material models using the local coordinate system for aspherical particle geometry; (2) a capability of modeling interface debonding in TRISO particles; and (3) an efficient high-fidelity approach to calculate the failure probability of TRISO particles. Benchmark problems and AGR-2 tests are used to verify the models and demonstrate these new capabilities.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Signatures of Velocity-dependent Dark Matter Self-interactions in Milky Way-mass Halos

In this work, we explore the impact of elastic, anisotropic, velocity-dependent dark matter (DM) self-interactions on the host halo and subhalos of Milky Way (MW)--mass systems. We consider a generic self-interacting dark matter (SIDM) model parameterized by the masses of a light mediator and the DM particle. The ratio of these masses, w, sets the velocity scale above which momentum transfer due to DM self-interactions becomes inefficient. We perform high-resolution zoom-in simulations of an MW-mass halo for values of w that span scenarios in which self-interactions either between the host and its subhalos or only within subhalos efficiently transfer momentum, and we study the effects of self-interactions on the host halo and on the abundance, radial distribution, orbital dynamics, and density profiles of subhalos in each case. The abundance and properties of surviving subhalos are consistent with being determined primarily by subhalo--host halo interactions. In particular, subhalos on radial orbits in models with larger values of the cross section at the host halo velocity scale are more susceptible to tidal disruption owing to mass loss from ram pressure stripping caused by self-interactions with the host. This mechanism suppresses the abundance of surviving subhalos relative to collisionless DM simulations, with stronger suppression for larger values of w. Thus, probes of subhalo abundance around MW-mass hosts can be used to place upper limits on the self-interaction cross section at velocity scales of ~200 km s-1, and combining these measurements with the orbital properties and internal dynamics of subhalos may break degeneracies among velocity-dependent SIDM models.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1974

Elastic moduli of precompressed pyrophyllite used in ultrahigh pressure research

The propagation of ultrasonic pulses in pyrophyllite specimens was studied to determine the effect of specimen precompression on the measured elastic moduli. Measurements were made at room pressure and, for the precompressed specimens, to pressures of 3 kbar. Pyrophyllite was found to be elastically anisotropic, apparently the result of the fabric present in our material. The room pressure adiabatic bulk modulus as measured on specimens made of isostatically compacted powered pyrophyllite was determined to be 96.1 kbar. The wave speeds of ultrasonic pulses in pyrophyllite were found to decrease with increasing specimen precompression. A limiting value of precompression was found, above which no further decrease in wave speed was observed. For the shear wave speeds this occurs at 10 kbar while for the longitudinal wave at 25 kbar. In the limit, the shear waves propagate 20% slower than in the unprecompressed samples; for the longitudinal wave the difference is 30%.

Sachse, W.↗
NASA NTRS1975

Structural design and analysis. Part 1

Topics in anisotropic elasticity are discussed along with the analysis of truss, beam, frame, and membrane components. Attention is given to failure criteria and failure analysis of composite structural components. Other subjects considered are related to the analysis of plates, the analysis of shells, and questions of wave propagation and impact in composite materials. Individual items are announced in this issue.

Chamis, C. C.↗
NASA NTRS1976

Fundamental analysis of the failure of polymer-based fiber reinforced composites

A mathematical model is described which will permit predictions of the strength of fiber reinforced composites containing known flaws to be made from the basic properties of their constituents. The approach was to embed a local heterogeneous region (LHR) surrounding the crack tip into an anisotropic elastic continuum. The model should (1) permit an explicit analysis of the micromechanical processes involved in the fracture process, and (2) remain simple enough to be useful in practical computations. Computations for arbitrary flaw size and orientation under arbitrary applied load combinations were performed from unidirectional composites with linear elastic-brittle constituent behavior. The mechanical properties were nominally those of graphite epoxy. With the rupture properties arbitrarily varied to test the capability of the model to reflect real fracture modes in fiber composites, it was shown that fiber breakage, matrix crazing, crack bridging, matrix-fiber debonding, and axial splitting can all occur during a period of (gradually) increasing load prior to catastrophic fracture. The computations reveal qualitatively the sequential nature of the stable crack process that precedes fracture.

Kanninen, M. F.↗
NASA NTRS1977

Preliminary development of a fundamental analysis model for crack growth in a fiber reinforced composite material

A mathematical model for the strength of fiber reinforced composites containing specific flaws is described. The approach is to embed a local heterogeneous region surrounding the crack tip in an anisotropic elastic continuum. By consideration of the individual failure events activated near the flaw tip, a strength prediction can be made from basic properties of the composite constituents. Computations for arbitrary flaw size and orientation have been performed for unidirectional composites with linear elastic-brittle constituent behavior. The mechanical properties were those of graphite epoxy. With the rupture properties arbitrarily varied to test the capability of the model to reflect real fracture modes in fiber composites, it is shown that fiber breakage, matrix crazing, crack bridging, matrix-fiber debonding, and axial splitting all can occur during gradually increasing load prior to catastrophic fracture. Qualitative comparisons with experimental results on edge-notched unidirectional graphite epoxy specimens have also been made.

Kanninen, M. F.↗
NASA NTRS1979

Boundary layer thermal stresses in angle-ply composite laminates

Boundary-layer thermal stress singularities and distributions of angle-ply composite laminates under uniform thermal loading are investigated through a system of sixth-order governing partial differential equations developed with the aid of the anisotropic elasticity field equations and Lekhnitskii's complex stress functions. Results are presented for cases of various angle-ply graphite/epoxy laminates, and it is shown that the boundary-layer thickness depends on the degree of anisotropy of each individual lamina, thermomechanical properties of each ply, and the relative thickness of adjacent layers.

Wang, S. S.↗
NASA NTRS1981

Edge delamination in angle-ply composite laminates

Edge delamination has caused severe concern in the design and analysis of advanced composite materials and structures. Due to its complex nature, very limited knowledge for the problem is currently available. It involves not only geometric and material discontinuities but also inherently coupled mode I, II and III fracture in the layered anisotropic system. Based on complex-variable stress potentials in the anisotropic elasticity theory and eigenfunction expansion, exact orders of the crack-tip stress singularity and complete field solutions are obtained. Results are given for edge delaminated, angle-ply composites subjected to uniform axial extension for illustrative purposes. Effects of geometric, lamination, and crack variables are determined.

Wang, S. S.↗
NASA NTRS1981

Edge delamination in angle-ply composite laminates, part 5

A theoretical method was developed for describing the edge delamination stress intensity characteristics in angle-ply composite laminates. The method is based on the theory of anisotropic elasticity. The edge delamination problem is formulated using Lekhnitskii's complex-variable stress potentials and an especially developed eigenfunction expansion method. The method predicts exact orders of the three-dimensional stress singularity in a delamination crack tip region. With the aid of boundary collocation, the method predicts the complete stress and displacement fields in a finite-dimensional, delaminated composite. Fracture mechanics parameters such as the mixed-mode stress intensity factors and associated energy release rates for edge delamination can be calculated explicity. Solutions are obtained for edge delaminated (theta/-theta theta/-theta) angle-ply composites under uniform axial extension. Effects of delamination lengths, fiber orientations, lamination and geometric variables are studied.

Wang, S. S.↗