Search NASA⌕ Search

SEARCH · Search NASA

Results for “STRESS CALCULATION”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Transverse Young's modulus of carbon/glass hybrid fiber composites

Hybrid fiber composites offer designers a means of tailoring the stress–strain behavior of lightweight materials used in high-performance structures. While the longitudinal stress–strain behavior of unidirectional hybrid fiber composites has been thoroughly evaluated experimentally and analytically, relatively little information is available on the transverse behavior. The objective of the current investigation is to present data on the transverse modulus of elasticity of unidirectional composites with five different ratios of carbon and glass fiber and to compare the data with predictive and fitted models. The transverse modulus increases monotonically with the proportion of glass fiber in the composite. Finite element analysis was used to evaluate different ways to model voids in the matrix and allowed the unknown transverse properties of the carbon fibers to be backed out using experimental data from the all-carbon composite. The finite element results show that the transverse modulus can be accurately modeled if voids are modeled explicitly in the matrix region and if modulus is calculated based on stress applied along the minimum interfiber distance path between adjacent fibers arranged in a rectangular array. The transverse modulus was under-predicted by the iso-stress model and was well predicted by a modified iso-stress model and a modified Halpin–Tsai model.

Materials Science↗
DOE OSTI2015

Dynamical Theory of Segmental Relaxation and Emergent Elasticity in Supercooled Polymer Melts

We generalize the force-level Elastically Collective Nonlinear Langevin Equation theory of supercooled molecular liquid dynamics to polymer melts based on mapping chains to disconnected and noninterpenetrating Kuhnsized spheres. This allows first-principles, no adjustable parameter calculations to be performed for the temperature-dependent mean segmental relaxation time of chemically diverse van der Waals polymers over a wide range of molecular weights. Despite the simplicity of the mapping, the theory does a good job of a priori predicting the glass transition temperature (T g ), the dynamic fragility, and full temperature dependence of the -relaxation time for some high molecular weight polymers and the chain length dependence of T g as the consequence of the molecular weight dependence of backbone stiffness. The minimalist model does not capture the unusually low and high fragilities of certain long chain polymers which are not typical of van der Waals molecular liquids. This seems likely due to the simple coarse graining adopted which ignores longer range chain connectivity and nonuniversal factors on the sub-Kuhn length scale. Elasticity, not of an entropic single chain origin, emerges in deeply supercooled polymer liquids due to transient segmental localization and is studied at the microscopic stress-tensor level. Calculations of the frequency-dependent dynamic storage modulus close to T g appear to be qualitatively consistent with recent measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Estimation of Errors in Stress Distributions Computed in Finite Element Simulations of Polycrystals

The accuracy of the stresses predicted from crystal plasticity-based finite element formulation depends on estimation and control of the errors associated with the discretization. In the current work, the errors in the stress distribution are estimated in virtual polycrystalline samples of α-phase titanium (hexagonal close-packed phase of Ti–6Al– 4 V). To estimate the error, the stress field, which does not possess inter-element continuity, is smoothed over a grain using an L 2 projection, thereby providing continuous stress distributions with inter-element continuity. The differences between the continuous (smooth) and discontinuous (raw) stress fields are calculated at individual Gauss quadrature points and used to estimate errors for corresponding elements and grains. Error estimations are performed for a Voronoi-tessellated microstructure, an equiaxed microstructure, and two microstructures with varying grain sizes for tensile loading extending into the fully plastic regime (≈ 5% extension). Magnitudes of the errors are found to depend on microstructural characteristics, particularly the shape and size of grains. Samples having variations in grain size or having less spherical grains exhibited higher errors than samples with uniformly sized, equiaxed grains, with the size variations having a more pronounced effect. In conclusion, errors correlate with proximity to grain boundaries at small (elastic) strains and with deformation-induced features (deformation bands) at large strains.

42 ENGINEERING↗
DOE OSTI2018

Radial Inertia Effect on Dynamic Compressive Response of Polymeric Foam Materials

Polymeric foams have been extensively used in shock isolation applications because of their superior shock or impact energy absorption capability. However, as a type of soft condensed matter, the highly nonlinear, heterogeneous, and dissipative behavior of polymeric foams may result in an ineffective mitigation or isolation to shock/blast loading. To meet certain desired shock mitigation or isolation requirements, the polymeric foams need to be experimentally characterized to obtain their intrinsic material response. However, radial inertia during dynamic compression has become a severe issue and needs to be fully understood. In this study, we developed an analytical method to calculate the additional stress induced by radial inertia in a polymeric foam specimen. The radial inertia is generally caused by Poisson’s effect and associated with three different mechanisms – axial strain acceleration, large deformation, and Poisson’s ratio change. In conclusion, the effect of Poisson’s ratio change during deformation on radial inertia was specifically investigated for hyperelastic foam materials, and verified with experimental results obtained from Kolsky compression bar tests on a silicone foam.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Long Term Creep Behavior of a CoCrFeNiMn High Entropy Alloy

The potential of high entropy alloys (HEAs) to meet or exceed austenitic stainless-steel performance with the additional benefit of improved hot corrosion/oxidation resistance makes FCC HEAs attractive for use in energy applications. While shorter term creep tests have been reported in the literature on HEAs, not all methodologies utilize repeatable techniques. This manuscript reports on over 23,500 accumulated hours of tensile creep testing with adherence to ASTM standards on a melt solidified ingot of CoCrFeNiMn HEA converted to wrought plate using conventional thermo-mechanical processing techniques. The typical standard creep analyses are reported, i.e., Larson-Miller parameter, Monkman-Grant relationship, activation energy for creep, and creep stress exponents were calculated and compared to previously reported short-term creep tests. Additionally, characteristics of creep fracture and microstructural evolution are reported with cursory dislocation mechanisms investigated.

Rozman, Kyle A.↗
DOE OSTI2012

Crack tip plasticity in single crystal UO2: Atomistic simulations

The fracture behavior of single crystal uranium dioxide is studied using molecular dynamics simulations at room temperature. Initially, an elliptical notch is created on either {111} or {110} planes, and tensile loading is applied normal to the crack planes. For cracks on both planes, shielding of crack tips by plastic deformation is observed, and crack extension occurs for crack on {111} planes only. Two plastic processes, dislocation emission and phase transformation are identified at crack tips. The dislocations have a Burgers vector of ?110?/2, and glide on {100} planes. Two metastable phases, the so-called Rutile and Scrutinyite phases, are identified during the phase transformation, and their relative stability is confirmed by separate density- functional-theory calculations. Examination of stress concentration near crack tips reveals that dislocation emission is not an effective shielding mechanism. The formation of new phases may effectively shield the crack provided all phase interfaces formed near the crack tips are coherent, as in the case of cracks residing on {110} planes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Local residual stresses and microstructure within recrystallizing grains in iron

The intragranular strains and the microstructure within recrystallizing grains in single phase pure iron have been systematically studied in 3D using synchrotron micro-diffraction. In contradiction to common knowledge we observe that residual elastic strains of up to 1–2 × 10 –3 are present in the recrystallizing grains. The local strain variations within individual grains can be similar in magnitude to those between different grains. The effects of processing parameters, including rolling reduction, recrystallized volume fraction, and cooling rate, as well as microstructural parameters, including grain size, orientation and internal misorientation, on the development of local residual elastic strains and microstructure have been quantified. It is suggested that the difference in local defect density between recrystallized grains and the surrounding deformed matrix, as well as the redistribution of medium to long-range residual stresses in the deformed grains, are the main reasons for the development of such strains in the recrystallized grains. The local residual stress levels are calculated and their possible impact on the mechanical properties and recrystallization processes are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2016

The comparison of microstructures and mechanical properties between 14Cr-Al and 14Cr-Ti ferritic ODS alloys

In this study, two kinds of 14Cr ODS alloys (14Cr-Al and 14Cr-Ti) were investigated to reveal the different effects between Al and Ti on the microstructures and mechanical properties of 14Cr ferritic ODS alloys. The microstructure information such as grains, minor phases of these two alloys has been investigated by high-energy X-ray diffraction and transmission electron microscopy (TEM). The in situ synchrotron X-ray diffraction tensile test was applied to investigate the mechanical properties of these two alloys. The lattice strains of different phases through the entire tensile deformation process in these two alloys were analyzed to calculate their elastic stresses. From the comparison of elastic stress, the strengthening capability of Y 2 Ti 2 O 7 is better than TiN in 14Cr-Ti, and the strengthening capability of YAH is much better than YAM and AlN in 14Cr-Al ODS. The dislocation densities of 14Cr-Ti and 14Cr-Al ODS alloys during tensile deformation were also examined by modified Williamson-Hall analyses of peak broadening, respectively. In conclusion, the different increasing speed of dislocation density with plastic deformation reveals the better strengthening effect of Y-Ti-O particles in 14Cr-Ti ODS than that of Y-Al-O particles in 14Cr-Al ODS alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Understanding the origin of softness in structurally tailored and engineered macromolecular (STEM) gels: A DPD study

Structurally tailored and engineered macromolecular (STEM) gels are polymer networks containing a primary network (the “STEM-0” gel) and side chains that are anchored to this network (forming the “STEM-1” gel). The ability to vary the features of the side chains provides a robust means of tailoring the macroscopic properties of the material. We use dissipative particle dynamics (DPD) simulations to determine the mechanical response of STEM-1 gels to uniaxial compression for various values of the side chain length, n sc . While the Young’s modulus of the material is significantly decreased with increasing n sc (at a fixed density of grafted side chains), above a certain saturation value, increasing n sc does not lead to any further softening of the sample. Using the simulations, we calculate the relevant stresses, which are related to the Young’s modulus and the number of entanglements in the network. We show that the backbone chains become more spread out upon addition of sidechains, leading to a decrease in the physical entanglements between backbone segments. This observation accounts for the decrease of stress in the backbone network and the softening of the networks through the addition of sidechains. For long chains, however, the number density of physical crosslinks between side chains is independent of side chain length; this observation explains the observed stress saturation in the STEM-1 gels. Our approach allows us to correlate the molecular architecture of the gels to the resultant macroscopic mechanical behavior and provide guidelines for fabricating STEM gels with well-defined mechanical properties, which allow the materials to be used for a variety of applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Three-dimensional paganica fault morphology obtained from hypocenter clustering (L'Aquila 2009 seismic sequence, Central Italy)

In seismic modelling, fault planes are normally assumed to be flat due to the lack of data which can constrain fault morphology. However, incorporating 3D fault morphology is important for modelling several phenomena, for example calculating mainshock induced stress changes. Here we utilize a data-analytical method to unveil the 3D rupture morphology of faults using unsupervised clustering techniques applied to earthquake hypocenters in seismic sequences. We apply this method to the 2009 L'Aquila seismic sequence which involved a M W 6.1 mainshock on April 6th. We use a dataset of about 50,000 relocated events, mostly microearthquakes, reaching magnitude of completeness equal to 0.7. Clustering distinguishes the earthquakes as occurring in three main clusters along with other minor fault segments. We then represent the morphology of the main Paganica fault system (responsible for the largest mainshock) using splines. This method shows promise as a step toward robustly and quickly obtaining 3D rupture morphologies where earthquake sequences have been monitored. The 3D model is presented interactively online, and the processing is presented in an interactive Jupyter Notebook (https://bit.ly/2MnCFdj).

58 GEOSCIENCES↗
DOE OSTI2014

Numerical exploration of spontaneous broken symmetries in multiorbital Hubbard models

Here, we study three proposals for broken symmetry in the cuprate pseudogap-oxygen antiferromagnetism, Theta(II) orbital loop currents, and circulating currents involving apex oxygens-through numerical exploration of multiorbital Hubbard models. Our numerically exact results show no evidence for the existence of oxygen antiferromagnetic order or the Theta(II) phase in the three-orbital Hubbard model. The model also fails to sustain an ordered current pattern even with the presence of additional apex oxygen orbitals. Thus, we conclude that it is difficult to stabilize the aforementioned phases in the multiorbital Hubbard models for parameters relevant to cuprate superconductors. However, the Theta(II) phase might be stabilized through explicit flux terms. We also found an enhanced propensity for circulating currents with such terms in calculations simulating applied stress or strain, which skew the copper-oxygen plane to resemble a kagome lattice. We propose an experimental viewpoint to shed additional light on this problem.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Prediction of the stability and electronic properties of carbon nanotori synthesized by a high-voltage pulsed discharge in ethanol vapor

An experimental technique for increasing the yield of carbon-nanotube nanotori using the modified arc synthesis method is proposed. New physical knowledge on the systematic features of the interrelation between the properties of nanotori and atomic-network topology are theoretically established for the first time. The experiments are performed based on new technology for synthesizing nanotori on nickel-catalyst particles by a high-voltage pulsed discharge in ethanol vapor and using atomic-force microscopy. Stability is predicted using an original procedure for calculating local atomic stresses. Simulation shows that the zigzag chirality corresponds to the most stable topology of nanotori. Using the tight binding method, it is shown that, depending on the chirality type, nanotori are divided into two classes, i.e., those with metal and semiconductor conductivity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Multi-Objective Optimization of a Wrought Magnesium Alloy for High Strength and Ductility

An optimization technique is coupled with crystal plasticity based finite element (CPFE) computations to aid the microstructural design of a wrought magnesium alloy for improved strength and ductility. The initial microstructure consists of a collection of sub-micron sized grains containing deformation twins. The variables used in the simulations are crystallographic texture, and twin spacing within the grains. It is assumed that plastic deformation occurs mainly by dislocation slip on two sets of slip systems classified as hard and soft modes. The hard modes are those slip systems that are inclined to the twin planes and the soft mode consists of dislocation glide along the twin plane. The CPFE code calculates the stress-strain response of the microstructure as a function of the microstructural parameters and the length-scale of the features. A failure criterion based on a critical shear strain and a critical hydrostatic stress is used to define ductility. The optimization is based on the sequential generation of an initial population defined by the texture and twin spacing variables. The CPFE code and the optimizer are coupled in parallel so that new generations are created and analyzed dynamically. In each successive generation, microstructures that satisfy at least 90% of the mean strength and mean ductility in the current generation are retained. Multiple generation runs based on the above procedure are carried out in order to obtain maximum strength-ductility combinations. Finally, the implications of the computations for the design of a wrought magnesium alloy are discussed. Research sponsored by the Laboratory Directed Research and Development Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the U. S. Department of Energy.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Multi-Objective Optimization of a Wrought Magnesium Alloy for High Strength and Ductility

An optimization technique is coupled with crystal plasticity based finite element (CPFE) computations to aid the microstructural design of a wrought magnesium alloy for improved strength and ductility. The initial microstructure consists of a collection of sub-micron sized grains containing deformation twins. The variables used in the simulations are crystallographic texture, and twin spacing within the grains. It is assumed that plastic deformation occurs mainly by dislocation slip on two sets of slip systems classified as hard and soft modes. The hard modes are those slip systems that are inclined to the twin planes and the soft mode consists of dislocation glide along the twin plane. The CPFE code calculates the stress-strain response of the microstructure as a function of the microstructural parameters and the length-scale of the features. A failure criterion based on a critical shear strain and a critical hydrostatic stress is used to define ductility. The optimization is based on the sequential generation of an initial population defined by the texture and twin spacing variables. The CPFE code and the optimizer are coupled in parallel so that new generations are created and analyzed dynamically. In each successive generation, microstructures that satisfy at least 90% of the mean strength and mean ductility in the current generation are retained. Multiple generation runs based on the above procedure are carried out in order to obtain maximum strength-ductility combinations. The implications of the computations for the design of a wrought magnesium alloy are discussed. Research sponsored by the Laboratory Directed Research and Development Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the U. S. Department of Energy.

Radhakrishnan, Balasubramaniam↗
DOE OSTI2015

Residual stress within nanoscale metallic multilayer systems during thermal cycling

Projected applications for nanoscale metallic multilayers will include wide temperature ranges. Since film residual stress has been known to alter system reliability, stress development within new film structures with high interfacial densities should be characterized to identify potential long-term performance barriers. To understand factors contributing to thermal stress evolution within nanoscale metallic multilayers, stress in Cu/Nb systems adhered to Si substrates was calculated from curvature measurements collected during cycling between 25 °C and 400 °C. Additionally, stress within each type of component layers was calculated from shifts in the primary peak position from in-situ heated X-ray diffraction. The effects of both film architecture (layer thickness) and layer order in metallic multilayers were tracked and compared with monolithic Cu and Nb films. Analysis indicated that the thermoelastic slope of nanoscale metallic multilayer films depends on thermal expansion mismatch, elastic modulus of the components, and also interfacial density. The layer thickness (i.e. interfacial density) affected thermoelastic slope magnitude while layer order had minimal impact on stress responses after the initial thermal cycle. When comparing stress responses of monolithic Cu and Nb films to those of the Cu/Nb systems, the nanoscale metallic multilayers show a similar increase in stress above 200 °C to the Nb monolithic films, indicating that Nb components play a larger role in stress development than Cu. Local stress calculations from X-ray diffraction peak shifts collected during heating reveal that the component layers within a multilayer film respond similarly to their monolithic counterparts.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Hybrid immersed interface-immersed boundary methods for AC dielectrophoresis

Dielectrophoresis, a nonlinear electrokinetic transport mechanism, has become popular in many engineering applications including manipulation, characterization and actuation of biomaterials, particles and biological cells. In this paper, we present a hybrid immersed interface–immersed boundary method to study AC dielectrophoresis where an algorithm is developed to solve the complex Poisson equation using a real variable formulation. An immersed interface method is employed to obtain the AC electric field in a fluid media with suspended particles and an immersed boundary method is used for the fluid equations and particle transport. The convergence of the proposed algorithm as well as validation of the hybrid scheme with experimental results is presented. In this paper, the Maxwell stress tensor is used to calculate the dielectrophoretic force acting on particles by considering the physical effect of particles in the computational domain. Thus, this study eliminates the approximations used in point dipole methods for calculating dielectrophoretic force. A comparative study between Maxwell stress tensor and point dipole methods for computing dielectrophoretic forces are presented. The hybrid method is used to investigate the physics of dielectrophoresis in microfluidic devices using an AC electric field. The numerical results show that with proper design and appropriate selection of applied potential and frequency, global electric field minima can be obtained to facilitate multiple particle trapping by exploiting the mechanism of negative dielectrophoresis. Our numerical results also show that electrically neutral particles form a chain parallel to the applied electric field irrespective of their initial orientation when an AC electric field is applied. This proposed hybrid numerical scheme will help to better understand dielectrophoresis and to design and optimize microfluidic devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Insights into the physical and chemical properties of a cement-polymer composite developed for geothermal wellbore applications

To isolate injection and production zones from overlying formations and aquifers during geothermal operations, cement is placed in the annulus between well casing and formations. Wellbore cement can fracture due to chemical and physical stress. A polymer-cement (composite) with self-healing properties, stable at high temperatures, was recently developed by our group. Here we present data on the elemental distribution, mineral composition, internal microstructure, and tensile elasticity of this composite and how it compares to cement without polymer (control cement). Entrapped air porosity was similar between control cement (0.44% by volume) and composite (0.22% by volume). Total voids increased with polymer content in the cement (10% and 25% composite), and polymer was found to fill most of these voids. Polymer appears to act as a (cement hydration) retarder increasing the amount of unhydrated silicate minerals in the matrix. Tensile elasticity was improved for composites as compared to control cement, suggesting that the composite materials perform better under mechanical stress. Density-functional calculations reproduce the drop in Young’s modulus by about 60% observed experimentally when polymer is added to the cement. In addition, calculations indicate that the bonding interactions between the cement and polymer remain stable over the range of strain studied.

Rod, Kenton A. (ORCID:0000000318189094)↗