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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1984

On constitutive relations at finite strain - Hypo-elasticity and elasto-plasticity with isotropic or kinematic hardening

Nagtegaal and de Jong (1982) have studied stresses generated by simple finite shear in the case of elastic-plastic and rigid-plastic materials which exhibit anisotropic hardening. They reported that the shear stress is oscillatory in time. It was found that the occurrence of such an 'anomaly' is not restricted to anisotropic plasticity. Similar behavior in finite shear may result even in the case of hypoelasticity and classical isotropic hardening plasticity theory. The present investigation is concerned with the central problem of 'generalizing' with respect to the finite strain case, taking into account the constitutive relations of infinitesimal strain theories of classical plasticity with isotropic or kinematic hardening. The problem of hypoelasticity is also considered. It is shown that current controversies surrounding the choice of stress rate in the finite-strain generalizations of the constitutive relations and the anomalies surrounding kinematic hardening plasticity theory are easily resolvable.

Atluri, S. N.↗
NASA NTRS1992

Failure analysis of thick composite cylinders under external pressure

Failure of thick section composites due to local compression strength and overall structural instability is treated. Effects of material nonlinearity, imperfect fiber architecture, and structural imperfections upon anticipated failure stresses are determined. Comparisons with experimental data for a series of test cylinders are described. Predicting the failure strength of composite structures requires consideration of stability and material strength modes of failure using linear and nonlinear analysis techniques. Material strength prediction requires the accurate definition of the local multiaxial stress state in the material. An elasticity solution for the linear static analysis of thick anisotropic cylinders and rings is used herein to predict the axisymmetric stress state in the cylinders. Asymmetric nonlinear behavior due to initial cylinder out of roundness and the effects of end closure structure are treated using finite element methods. It is assumed that local fiber or ply waviness is an important factor in the initiation of material failure. An analytical model for the prediction of compression failure of fiber composites, which includes the effects of fiber misalignments, matrix inelasticity, and multiaxial applied stresses is used for material strength calculations. Analytical results are compared to experimental data for a series of glass and carbon fiber reinforced epoxy cylinders subjected to external pressure. Recommendations for pretest characterization and other experimental issues are presented. Implications for material and structural design are discussed.

Caiazzo, A.↗
DOE OSTI2023

Effects of anisotropic shale creep on the stress and permeability evolution of a geological nuclear waste repository

To ensure long-term safety and performance, geological nuclear waste repositories require low-permeability barriers such as bentonite buffers and/or shale host rock. Shale is not only known for its low permeability but also for its trend to undergo time-dependent deformation (i.e., creep), which could heal damage, but the effects of shale creep on the long-term performance of nuclear waste repositories have not been clearly understood. In particular, the anisotropic nature of shale (i.e., bedding) could have a significant effect on its creep behavior, and consequently, on the long-term performance of nuclear waste repositories. In this research, numerical simulations were carried out with the objective of showing the effects of anisotropic shale creep on the stress and permeability evolution of a generic geological nuclear waste repository in shale. The TOUGH-FLAC simulator, a thermo-hydromechanically (THM) coupled numerical code, was used for the simulations. To achieve the objective, comparisons were performed between the results of anisotropic shale creep simulations and those of different simulation cases, namely, no creep (i.e., elastic), isotropic creep, and long-term creep shale cases. Results of the comparisons show that the elastic and isotropic creep shale cases respectively led to the overestimation and underestimation of stress and permeability in the repository, whereas the long-term creep shale case, which accumulated greater creep in later periods than in earlier periods, helped to keep large shear and tensile stresses from developing while maintaining compressive spherical stress, resulting in consistently low permeability levels. These results indicate that performance assessments with elastic and isotropic creep formation models will provide the upper and lower bound estimates of stress and permeability, while more reasonable estimates will be given by an anisotropic creep formation model, and that shale with long-term creep characteristics will be beneficial in many aspects of the safety and performance of nuclear waste repositories.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Sliding dynamics for bubble phases on periodic modulated substrates

We analyze a bubble-forming system composed of particles with competing long-range repulsive and short-range attractive interactions driven over a quasi-one-dimensional periodic substrate. We find various pinned and sliding phases as a function of substrate strength and drive amplitude. When the substrate is weak, a pinned bubble phase appears that depins elastically into a sliding bubble lattice. For stronger substrates, we find anisotropic bubbles, disordered bubbles, and stripe phases. Plastic depinning occurs via the hopping of individual particles from one bubble to the next in a pinned bubble lattice, and as the drive increases, there is a transition to a state where all of the bubbles are moving but are continuously shedding and absorbing individual particles. This is followed at high drives by a moving bubble lattice in which the particles can no longer escape their individual bubbles. The transition between the plastic and elastic sliding phases can be detected via signatures in the velocity-force curves, differential conductivity, and noise. When the bubbles shrink due to an increase in the attractive interaction term, they fit better inside the pinning troughs and become more strongly pinned, leading to a reentrant pinning phase. For weaker attractive terms, the size of the bubbles becomes greater than the width of the pinning troughs and the depinning becomes elastic with a reduced depinning threshold. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Sliding Dynamics for Bubble Phases on Periodic Modulated Substrates

We analyze a bubble-forming system composed of particles with competing long-range repulsive and short range attractive interactions driven over a quasi-one-dimensional periodic substrate. We find various pinned and sliding phases as a function of substrate strength and drive amplitude. When the substrate is weak, a pinned bubble phase appears that depins elastically into a sliding bubble lattice. For stronger substrates, we find anisotropic bubbles, disordered bubbles, and stripe phases. Plastic depinning occurs via the hopping of individual particles from one bubble to the next in a pinned bubble lattice, and as the drive increases, there is a transition to a state where all of the bubbles are moving but are continuously shedding and absorbing individual particles. This is followed at high drives by a moving bubble lattice in which the particles can no longer escape their individual bubbles. The transition between the plastic and elastic sliding phases can be detected via signatures in the velocity-force curves, differential conductivity, and noise. When the bubbles shrink due to an increase in the attractive interaction term, they fit better inside the pinning troughs and become more strongly pinned, leading to a reentrant pinning phase. For weaker attractive terms, the size of the bubbles becomes greater than the width of the pinning troughs and the depinning becomes elastic with a reduced depinning threshold.

36 MATERIALS SCIENCE↗
NASA NTRS1998

Finite Element and Plate Theory Modeling of Acoustic Emission Waveforms

A comparison was made between two approaches to predict acoustic emission waveforms in thin plates. A normal mode solution method for Mindlin plate theory was used to predict the response of the flexural plate mode to a point source, step-function load, applied on the plate surface. The second approach used a dynamic finite element method to model the problem using equations of motion based on exact linear elasticity. Calculations were made using properties for both isotropic (aluminum) and anisotropic (unidirectional graphite/epoxy composite) materials. For simulations of anisotropic plates, propagation along multiple directions was evaluated. In general, agreement between the two theoretical approaches was good. Discrepancies in the waveforms at longer times were caused by differences in reflections from the lateral plate boundaries. These differences resulted from the fact that the two methods used different boundary conditions. At shorter times in the signals, before reflections, the slight discrepancies in the waveforms were attributed to limitations of Mindlin plate theory, which is an approximate plate theory. The advantages of the finite element method are that it used the exact linear elasticity solutions, and that it can be used to model real source conditions and complicated, finite specimen geometries as well as thick plates. These advantages come at a cost of increased computational difficulty, requiring lengthy calculations on workstations or supercomputers. The Mindlin plate theory solutions, meanwhile, can be quickly generated on personal computers. Specimens with finite geometry can also be modeled. However, only limited simple geometries such as circular or rectangular plates can easily be accommodated with the normal mode solution technique. Likewise, very limited source configurations can be modeled and plate theory is applicable only to thin plates.

Prosser, W. H.↗
DOE OSTI2024

Design of biphenylene-derived tunable dirac materials

The exploration of carbon allotropes has unveiled a series of two-dimensional (2D) materials with unique electronic and mechanical properties, yet the need for stable structures with tailored electronic properties persists. Here, in this study, we introduce a new class of 2D carbon allotropes derived from the biphenylene network (BPN), incorporating acetylenic linkages to tune their structural and electronic characteristics. Through density functional theory calculations, we identified ten novel BPN-derived structures that exhibit both energetic and dynamic stability, confirmed by cohesive energy and phonon spectrum analyses. Among them, BPN-02 and BPN04 are metallic, featuring critically-tilted type-III Dirac cones under ~ 5 % biaxial strain, while BPN-22 is a semiconductor with a band gap of 0.95 eV and exhibits highly anisotropic carrier mobility. Additionally, these structures demonstrate significant anisotropy in their elastic properties, further distinguishing them from other 2D carbon materials like graphene. Our findings suggest that these novel BPN-based structures have strong potential for next-generation electronic and optoelectronic applications, providing new avenues for the design and synthesis of advanced carbon materials.

2D carbon allotropes↗
DOE OSTI2021

Effect of neighboring grain orientation on strain localization in slip bands in HCP materials

Particularly in plastically anisotropic crystals, such as hexagonal close packed (HCP) materials, plastic deformation is realized by slip acting in the small volumes within individual crystals. Here we extend a full field fast Fourier transform (FFT)-based elasto-viscoplastic formulation to simulate the development of a single slip band on either prismatic or basal planes spanning a crystal. Calculations of the strain and stress fields induced locally within the band and parent crystal, and ahead of the band/grain boundary junction in the neighboring crystal are analyzed as the slip band intensifies under increasing applied strain. We report a substantial influence of the crystallographic orientation of the nearest neighboring grain on the rate of slip band localization. Performing the analysis on two materials, CP-Ti and Mg, indicates that the strength of the material affects the rate of localization, with stronger materials tending to localize more easily. A slip band tip stress-based criterion is proposed for identifying the nearest neighbor orientations in which slip band transmission is possible and the likely slip system for which it occurs. This indicator is validated against experimental studies on commercially pure Ti, an Mg–Y alloy, and Ti–6Al–4V. Finally, we show that for low GB misorientations, the slip band is likely to transmit into another slip band of the same type in the neighbor grain, while for high GB misorientations, it is likely to transmit into one of a different type or to not transmit at all.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Electronic, magnetic, and optical properties of bulk and (1 1 1)-surfaces of CoMnZnSi quaternary Heusler alloy

The structural stability, electronic structure, and magnetic and optical properties of CoMnZnSi quaternary Heusler alloy (QHA) in the bulk and (1 1 1)-slab forms have been investigated by performing density functional theory (DFT) calculations. Formation and cohesive energies, and elastic constants confirm that the bulk CoMnZnSi is chemically and mechanically stable at an equilibrium lattice parameter of 5.81 Å. The calculated elastic constants also indicate that this QHA has ductile and anisotropic features. We investigate and discuss the bonding behavior from charge density distribution and density of states. Furthermore, these calculations in the bulk phase show a perfect half-metallic behavior with an integer value of magnetic moment (4 µB) and a large spin-flip gap of 0.39 eV. On the other hand, the Co (1 1 1) and Si (1 1 1)-slabs exhibit semi-metallic nature at the Fermi level. These results indicate that CoMnZnSi qualifies for spintronic applications. In the bulk and Co (1 1 1), Zn (1 1 1), and Si (1 1 1)-surfaces, the magnetic moment of Co and Mn atoms are ferromagnetically aligned, while the ferrimagnetic alignment of Co with Mn has been found in Mn (1 1 1)-surface. We also investigate the optimal conductivity, dielectric functions, reflectivity, absorption coefficient, refractive index, and loss function to understand the underlying optical properties of this alloy.

36 MATERIALS SCIENCE↗
NASA NTRS2006

Magneto-thermo-elastokinetics of Geometrically Nonlinear Laminated Composite Plates. Part 1: Foundation of the Theory

A fully coupled magneto-thermo-elastokinetic model of laminated composite, finitely electroconductive plates incorporating geometrical nonlinearities and subjected to a combination of magnetic and thermal fields, as well as carrying an electrical current is developed, In this context. the first-order transversely shearable plate theory in conjunction with von-Karman geometrically nonlinear strain concept is adopted. Related to the distribution of electric and magnetic field disturbances within the plate, the assumptions proposed by Ambartsumyan and his collaborators are adopted. Based on the electromagnetic equations (i.e. the ones by Faraday, Ampere, Ohm, Maxwell and Lorentz), the modified Fourier's law of heat conduction and on the elastokinetic field equations, the 3-D coupled problem is reduced to an equivalent 2- D one. The theory developed herein provides a foundation for the investigation, both analytical and numerical, of the interacting effects among the magnetic, thermal and elastic fields in multi-layered thin plates made of anisotropic materials.

Hasanyan, Davresh↗
NASA NTRS2003

Buckling Of Long Compression-Loaded Anisotropic Plates Restrained Against Inplane Lateral and Shear Deformations

An approach for synthesizing buckling results and behavior for thin balanced and unbalanced symmetric laminates that are subjected to uniform axial compression loads and elastically restrained against inplane expansion, contraction, and shear deformation is presented. This approach uses a nondimensional analysis for infinitely long, flexurally anisotropic plates (coupling between bending and twisting) that are subjected to combined mechanical loads and is based on nondimensional parameters. In addition, nondimensional loading parameters are derived that account for the effects of the elastic inplane deformation restraints, membrane orthotropy, and membrane anisotropy on the induced prebuckling stress state. The loading parameters are used to determine buckling coefficients that include the effects of flexural orthotropy and flexural anisotropy. Many results are presented, for some selected laminates, that are intended to facilitate a structural designer's transition to the use of the generic buckling design curves that are presented and discussed in the paper. Several buckling response curves are presented that provide physical insight into the behavior for combined loads, in addition to providing useful design data. An example is presented that demonstrates the use of the generic design curves, which are applicable to a wide range of laminate constructions. The analysis approach and generic results indicate the effects and characteristics of laminate orthotropy and anisotropy in a very general and unifying manner.

Nemeth, Michael P.↗
DOE OSTI2023

A new approach to model geomaterials with heterogeneous properties in thermo-hydro-mechanical coupled problems

The main objective of this article is to present a new approach to model coupled thermo-hydro-mechanical problems considering geomaterials with heterogeneous properties. This approach has been implemented in the software CODE_BRIGHT and it provides the possibility of considering geomaterials with a spatially correlated heterogeneous field of porosity, following a normal distribution. This spatial correlation can be isotropic or anisotropic. An important feature of this approach is that material properties such as intrinsic permeability, thermal conductivity, diffusivity, retention curve, elastic modulus or cohesion are defined as a function of porosity and, thus, they become heterogeneous with spatial correlation and, eventually, anisotropic. A validation exercise and other basic numerical examples have been carried out to illustrate the possibilities of the proposed approach. The results, which have been compared with a homogeneous case, show that considering heterogeneous fields can be relevant in different modelling problems, especially coupled thermo-hydro-mechanical problems.

42 ENGINEERING↗
DOE OSTI2026

Short and medium range structure in elastic deformation of metallic and covalent glasses

Here, we present a concise methodology to analyze structural response to the applied stress in amorphous solids, including metallic glasses (MG), glassy selenium, silica and polycarbonate, using high energy x-ray diffraction and atomic pair distribution function (PDF) analysis. To assess the structural anisotropy induced by applied axial stress, diffraction data were expanded into spherical harmonics. Using Bessel transformation, components of the structure function were converted into isotropic and anisotropic PDFs. The PDFs were compared to the expected model behavior for ideal elastic deformation to separate homogeneous affine strain from local non-affine strains. In metallic glass the range of non-affine deformation is limited to the nearest neighbor shell, suggesting local strain relaxation under stress that occurs even in the elastic regime. Beyond the second atomic shell strain is uniform. However, in glassy silica, polycarbonate and selenium strong local bonding inhibits local displacements and strain in short range order is accommodated by rotation of local units. Interestingly, beyond a molecular unit, deformation in covalent systems is similar to MG, and response of the medium range order scales with the macroscopic stress.

glassy structure↗
DOE OSTI2021

Evolution of anisotropy with saturation and its implications for the elastoplastic responses of clay rocks

Many clay rocks have distinct bedding planes. Experimental studies have shown that their mechanical properties evolve with the degree of saturation (DOS), often with higher stiffness and strength after drying. For transversely isotropic rocks, the effects of saturation can differ between the bed-normal (BN) and bed-parallel (BP) directions, which gives rise to saturation-dependent stiffness and strength anisotropy. Accurate prediction of the mechanical behavior of clay rocks under partially saturated conditions requires numerical models that can capture the evolving elastic and plastic anisotropy with DOS. In this study, we present an anisotropy framework for coupled solid deformation-fluid flow in unsaturated elastoplastic media. Here we incorporate saturation-dependent strength anisotropy into an anisotropic modified Cam-Clay (MCC) model and consider the evolving anisotropy in both the elastic and plastic responses. The model was calibrated using experimental data from triaxial tests to demonstrate its capability in capturing strength anisotropy at various levels of saturation. Through numerical simulations, we demonstrate the role of evolving stiffness and strength anisotropy in the mechanical behavior of clay rocks. Plane strain simulations of triaxial compression tests were also conducted to demonstrate the impacts of material anisotropy and DOS on the mechanical and fluid flow responses.

36 MATERIALS SCIENCE↗
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↗
DOE OSTI2026

Effect of magneto-mechanical synergism in the process-structure correlation in Fe–C alloys: A phase-field modeling approach

Applied magnetic fields can alter phase equilibria and kinetics in steels; however, quantitatively resolving how magnetic, chemical, and elastic driving forces jointly influence the microstructure remains challenging. We develop a quantitative magneto-mechanically coupled phase-field model for the Fe–C system that couples a CALPHAD-based chemical free energy with demagnetization-field magnetostatics and microelasticity. Here, the model reproduces single- and multi-particle evolution during the α → γ inverse transformation at 1023 K under external fields up to 20 T, including ellipsoidal morphologies observed experimentally at 8 T. Chemically driven growth is isotropic; a magnetic interaction introduces an anisotropic driving force that elongates γ precipitates along the field into ellipsoids, while elastic coherency promotes faceting, yielding elongated cuboidal or “brick-like” particles under combined magneto-elastic coupling. Growth kinetics increase with C content, and decrease with field strength and misfit strain. Multi-particle simulations reveal dipolar interaction-mediated coalescence for field-parallel neighbors and ripening for field-perpendicular neighbors. Incorporating field-dependent diffusivity from experiment slows kinetics as expected; a first-principles-motivated anisotropic diffusivity correction is estimated to be small (<2%). These results establish a process-structure link for magnetically assisted heat treatments of Fe–C alloys and provide guidance for microstructure control via chemo-magneto-mechanical synergism.

Magnetic field↗
DOE OSTI2021

How to determine limiting velocities of dislocations in anisotropic crystals

In the continuum limit, the theory of dislocations in crystals predicts a divergence in the elastic energy of the host material at a crystal geometry dependent limiting (or critical) velocity vc. Explicit expressions for vc are scattered throughout the literature and are available in analytic form only for special cases with a high degree of symmetry. The fact that in some cases (like pure edge dislocations in fcc) v c happens to coincide with the lowest shear wave speed of a sound wave traveling parallel to the dislocation's gliding direction has led to further confusion in the more recent literature. The aim of this short review therefore is to provide a concise overview of the limiting velocities for dislocations of arbitrary character in general anisotropic crystals, and how to efficiently compute them, either analytically or numerically.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗