Search NASA⌕ Search

SEARCH · Search NASA

Results for “Anisotropic elastic materials”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modeling brittle fracture due to anisotropic thermal expansion in polycrystalline materials

Here, this work investigated brittle fracture of polycrystalline materials due to thermal stresses arising from anisotropic thermal expansion. We used phase-field fracture simulations with the properties of alpha-uranium (α-U) and assumed a linear elastic mechanical response. Three-dimensional simulations were used to predict fracture for various conditions and crystallographic textures. We found that fracture was more pronounced during cooling than during heating because the anisotropy increased with temperature. We also found that the total crack surface area increased with increasing average misorientation, while the net shape change of the material decreased with increasing misorientation. Two-dimensional simulations in which one crystallographic coefficient of thermal expansion (CTE) was set to zero indicated that the expansion behavior in the crystallographic direction with the smallest CTE was the primary cause of fracture.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1983

The mechanics of delamination in fiber-reinforced composite materials. Part 1: Stress singularities and solution structure

The fundamental mechanics of delamination in fiber composite laminates is studied. Mathematical formulation of the problem is based on laminate anisotropic elasticity theory and interlaminar fracture mechanics concepts. Stress singularities and complete solution structures associated with general composite delaminations are determined. For a fully open delamination with traction-free surfaces, oscillatory stress singularities always appear, leading to physically inadmissible field solutions. A refined model is introduced by considering a partially closed delamination with crack surfaces in finite-length contact. Stress singularities associated with a partially closed delamination having frictional crack-surface contact are determined, and are found to be diferent from the inverse square-root one of the frictionless-contact case. In the case of a delamination with very small area of crack closure, a simplified model having a square-root stress singularity is employed by taking the limit of the partially closed delamination. The possible presence of logarithmic-type stress singularity is examined; no logarithmic singularity of any kind is found in the composite delamination problem. Numerical examples of dominant stress singularities are shown for delaminations having crack-tip closure with different frictional coefficients between general (1) and (2) graphite-epoxy composites.

Wang, S. S.↗
DOE OSTI2024

Experimental Observation of a New Attenuation Mechanism in hcp ‐Metals That May Operate in the Earth's Inner Core

Abstract Seismic observations show the Earth's inner core has significant and unexplained variation in seismic attenuation with position, depth and direction. Interpreting these observations is difficult without knowledge of the visco‐ or anelastic dissipation processes active in iron under inner core conditions. Here, a previously unconsidered attenuation mechanism is observed in zinc, a low pressure analog ofhcp‐iron, during small strain sinusoidal deformation experiments. The experiments were performed in a deformation‐DIA combined with X‐radiography, at seismic frequencies (∼0.003–0.1 Hz), high pressure and temperatures up to ∼80% of melting temperature. Significant dissipation (0.077 ≤ Q −1 (ω) ≤ 0.488) is observed along with frequency dependent softening of zinc's Young's modulus and an extremely small activation energy for creep (⩽7 kJ mol −1 ). In addition, during sinusoidal deformation the original microstructure is replaced by one with a reduced dislocation density and small, uniform, grain size. This combination of behavior collectively reflects a mode of deformation called “internal stress superplasticity”; this deformation mechanism is unique to anisotropic materials and activated by cyclic loading generating large internal stresses. Here we observe a new form of internal stress superplasticity, which we name as “elastic strain mismatch superplasticity.” In it the large stresses are caused by the compressional anisotropy. If this mechanism is also active inhcp‐iron and the Earth's inner‐core it will be a contributor to inner‐core observed seismic attenuation and constrain the maximum inner‐core grain‐size to ≲10 km.

Geochemistry & Geophysics↗
NASA NTRS1973

A study of fracture phenomena in fiber composite laminates

The extension of linear elastic fracture mechanics from ostensibly homogeneous isotropic metallic alloys to heterogeneous anisotropic advanced fiber composites is considered. It is analytically demonstrated that the effects of material anisotropy do not alter the principal characteristics exhibited by a crack in an isotropic material. The heterogeneity of fiber composites is experimentally shown to have a negligible effect on the behavior of a sufficiently long crack. A method is proposed for predicting the fracture strengths of a large class of composite laminates; the values predicted by this method show good agreement with limited experimental data. The limits imposed by material heterogeneity are briefly discussed, and areas for further study are recommended.

Konish, H. J., Jr.↗
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 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.↗
DOE OSTI2021

3D orthorhombic earth model effects on seismic source characterization

Most earth materials are anisotropic with regard to seismic wave-speeds, especially materials such as shales, or where oriented fractures are present. However, the base assumption for many numerical simulations is to treat earth materials as isotropic media. This is done for simplicity, the apparent weakness of anisotropy in the far field, and the lack of well-characterized anisotropic material properties for input into numerical simulations. One approach for addressing the higher complexity of actual geologic regions is to model the material as an orthorhombic medium. We have developed an explicit time-domain, finite-difference (FD) algorithm for simulating three-dimensional (3D) elastic wave propagation in a heterogeneous orthorhombic medium. The objective of this research is to investigate the errors and biases that result from modeling a non-isotropic medium as an isotropic medium. This is done by computing “observed data” by using synthetic, anisotropic simulations with the assumption of an orthorhombic, anisotropic earth model. Green’s functions for an assumed isotropic earth model are computed and then used an inversion designed to estimate moment tensors with the “observed” data. One specific area of interest is how shear waves, which are introduced in an anisotropic model even for an isotropic explosion, affect the characterization of seismic sources when isotropic earth assumptions are made. This work is done in support of the modeling component of the Source Physics Experiment (SPE), a series of underground chemical explosions at the Nevada National Security Site (NNSS).

58 GEOSCIENCES↗
NASA NTRSDate not supplied

Porous Microstructure Analysis (PuMA) software

The Porous Microstructure Analysis (PuMA) software was developed to provide a robust and efficient framework for computing material properties based on their microstructures. The development was motivated by advancements in X-ray microtomography, an imaging technology that can resolve the structure of a material at a sub-micron scale, in 3D and even in 4D (over time). PuMA provides the capability of computing a comprehensive spectrum of properties, from the most fundamental geometric features of a microstructure, to advanced anisotropic thermo-elastic properties. In addition, the software can generate artificial microstructures, ranging from simple analytical shapes to complex fibrous woven and non- woven geometries, which can be used in performance optimization studies. This presentation will highlight many of the capabilities of the recent open-source release.

microtomography↗
NASA NTRS1998

Determination of Stress Coefficient Terms in Cracked Solids for Monoclinic Materials with Plane Symmetry at x3 = 0

Determination of all the coefficients in the crack tip field expansion for monoclinic materials under two-dimensional deformation is presented in this report. For monoclinic materials with a plane of material symmetry at x(sub 3) = 0, the in-plane deformation is decoupled from the anti-plane deformation. In the case of in-plane deformation, utilizing conservation laws of elasticity and Betti's reciprocal theorem, together with selected auxiliary fields, T-stress and third-order stress coefficients near the crack tip are evaluated first from path-independent line integrals. To determine the T-stress terms using the J-integral and Betti's reciprocal work theorem, auxiliary fields under a concentrated force and moment acting at the crack tip are used respectively. Through the use of Stroh formalism in anisotropic elasticity, analytical expressions for all the coefficients including the stress intensity factors are derived in a compact form that has surprisingly simple structure in terms of the Barnett-Lothe tensors, L. The solution forms for degenerated materials, orthotropic, and isotropic materials are presented.

Yuan, F. G.↗
DOE OSTI2021

Effect of interfacial structures on phonon transport across atomically precise Si/Al heterojunctions

Phonons are important carriers of energy and information in many cryogenic devices used for quantum information science and in fundamental physics experiments such as dark matter detectors. In these systems phonon behavior can be dominated by interfaces, hence there is increasing demand for a more detailed understanding of interfacial phonon transport in relevant material systems. Previous studies have focused on understanding thermal transport over the entire phonon spectrum at and above room temperature. At ultra-low temperatures, however, knowledge is missing regarding athermal phonon behavior due to the challenge in modeling the extreme conditions in micro-scale, heterogeneous cryogenic systems, as well as extracting single-phonon information from a large ensemble. In this work, we delineate interfacial phonon transport using a combination of classical molecular dynamics (MD) and phonon wave packet simulations, to illustrate the consistency and differences between the ensemble- and single-phonon dynamics. We consider three single-crystal Si surface reconstructions --- (1×1), (v3×v3) and (7×7) --- and model both experimentally observed Si(1×1)/Al interfaces and hypothesized Si(v3×v3)/Al and Si(7×7)/Al interfaces. The overall interfacial thermal conductance calculated from non-equilibrium MD shows that for the Si(1×1)/Al system, the presence of Al twin grain boundaries can hinder phonon transport and reduce thermal conductance by 2--12% relative to single-crystal Al; whereas the Si (v3×v3) and (7×7) reconstructions can enhance it by 6--19%. Normal mode decomposition reveals that both the increase and decrease in conductance are related to inelastic phonon scattering. Single-phonon wave packet simulations show results consistent with NEMD, while further suggesting that phonon polarization conversion is significant even when elastic transmission dominates, and that the interfacial structures have anisotropic impacts on atomic vibrations along different lattice directions. Our findings suggest alternative perspectives in achieving selective phonon transport through controlling interfacial structures of materials with atomically-precise fabrication techniques, and that the phonon wave packet formalism is a potentially powerful method for developing a detailed understanding of non-equilibrium phenomena in the low-temperature limit.

Lu, Zexi↗
DOE OSTI2022

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Computationally efficient method for determining limiting velocities of edge dislocations in anisotropic crystals

The continuum-limit theory of dislocations in crystals predicts divergences in the elastic energy at crystal-geometry dependent limiting velocities vL, which separate subsonic, transsonic, and supersonic dislocation glide regimes and are therefore import for material strength models at high strain rates. Although it is known how to calculate those limiting velocities, there is one special case - edge dislocations with reflection symmetry, but non-vanishing elastic constants c16 or c26 - where previous methods have been notoriously slow. In this letter, we address this deficiency by deriving a computationally efficient method for determining the limiting velocities of edge dislocations with reflection symmetry which is two orders of magnitude faster than the previous method.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Structure and stability of hydrous minerals at high pressure

The presence of even small amounts of hydrogen in the Earth's deep interior may have profound effects on mantle melting, rheology, and electrical conductivity. The recent discovery of a large class of high-pressure H-bearing silicates further underscores the potentially important role for hydrous minerals in the Earth's mantle. Hydrogen may also be a significant component of the Earth's core, as has been recently documented by studies of iron hydride at high pressure. In this study, we explore the role of H in crystal structures at high pressure through detailed Raman spectroscopic and x ray diffraction studies of hydrous minerals compressed in diamond anvil cells. Brucite, Mg(OH)2, has a simple structure and serves as an analogue for the more complex hydrous silicates. Over the past five years, this material has been studied at high pressure using shock-compression, powder x ray diffraction, infrared spectroscopy, Raman spectroscopy, and neutron diffraction. In addition, we have recently carried out single-crystal synchrotron x-ray diffraction on Mg(OH)2 and Raman spectroscopy on Mg(OD)2 at elevated pressure. From all these studies, an interesting picture of the crystal chemical behavior of this material at high pressure is beginning to emerge. Some of the primary conclusions are as follows: First, hydrogen bonding is enhanced by the application of pressure. Second, layered minerals which are elastically anisotropic at low pressure may not be so at high pressure. Furthermore, the brucite data place constraints on the effect of hydrogen on seismic velocities and density at very high pressure. Third, the stability of hydrous minerals may be enhanced at high P by subtle structural rearrangements that are difficult to detect using traditional probes and require detailed spectroscopic analyses. Finally, brucite appears to be unique in that it undergoes pressure-induced disordering that is confined solely to the H-containing layers of the structure.

Duffy, T. S.↗
DOE OSTI2020

Cell alignment by smectic liquid crystal elastomer coatings with nanogrooves

Control of cells behavior through topography of substrates is a vital theme in biomedical applications. Among many materials used as substrates, polymers show advantages since they can be tailored by chemical functionalization. Fabrication of polymer substrates with nano- and microscale topography requires processing by lithography, microprinting, etching, etc. In this work, we present a different approach based on anisotropic elastic properties of polymerized smectic A (SmA) liquid crystal elastomer (LCE). When the SmA liquid crystal coating is deposited onto a substrate with planar alignment of the molecules, it develops nanogrooves at its free surface. After photopolymerization, these nanogrooves show an excellent ability to align human dermal fibroblasts (hDFs) over large areas. The alignment quality is good for both bare SmA LCE substrates and for substrates coated with fibronectin. The SmA LCE nano-topographies show a high potential for tissue engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Fractonic gauge theory of smectics

Motivated by striped correlated quantum matter, and the recently developed duality between elasticity of a two-dimensional (2D) crystal and a gauge theory, we derive a dual coupled U(1) vector gauge theory for a two-dimensional (2D) quantum smectic, where the disclination is mapped onto the fractonic charge, that we demonstrate can only move transversely to smectic layers. This smectic gauge theory dual also emerges from a gauge dual of a quantum crystal through a Higgs transition corresponding to a single flavor of its dipole condensation, an anisotropic quantum melting via dislocation proliferation. A condensation of the second flavor of dislocations a corresponds to another Higgs transition describing the smectic-to-nematic melting. We also utilize the electrostatic limit of this duality to formulate a melting of a 2D classical smectic in terms of a higher derivative sine–Gordon model, demonstrating its instability to a nematic at any nonzero temperature. Generalizing this classical duality to a 3D smectic, gives formulation of a 3D nematic-to-smectic transition in terms of an anisotropic Abelian-Higgs model.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A machine learning approach to determine the elastic properties of printed fiber-reinforced polymers

This work focuses on the simultaneous determination of the elastic constants and the fiber orientation state for a short fiber-reinforced polymer composite by performing a minimum of experimental tests. Here we introduce a methodology that enables the inverse determination of fiber orientation state and the in-situ polymer properties by performing tensile tests at the composite coupon level. We demonstrate the approach for the extrusion deposition additive manufacturing (EDAM) process to illustrate one application of the methodology, but the development is such that it can be applied to short fiber-reinforced polymer (SFRP) systems processed via other methods. Currently, developing composites additive manufacturing digital twins require extensive material characterization. In particular, the mechanical characterization of the orthotropic elastic properties of a composite involves extensive sample preparation and testing, therefore the elasticity tensor is generally populated using a micromechanics model. This, however, requires measuring the fiber orientation state in addition to knowing the constituent material properties. Experimentally measuring the fiber orientation state can be tedious and time consuming. Further, optical methods are limited to resolving the orientation of cylindrical fibers or cluster of non-cylindrical fibers, and computed tomography (CT) methods scan regions of volume that are much smaller than a full printed bead. Therefore, we propose a methodology, accelerated by machine learning, to identify the anisotropic mechanical properties and fiber orientation state at the same time. Early results show that inference of the fiber orientation and composite properties is possible with as few as three tensile tests. Our results show that a combination of the choice of the micromechanics model and reliable set of experiments can yield the nine elastic constants, as well as, the fiber orientation state.

36 MATERIALS SCIENCE↗