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

Aftershocks of the 2012 Off-Coast of Sumatra Earthquake Sequence

Aftershocks of the 2012 Off-Coast of Sumatra Earthquake Sequence exhibit a complex and diffuse spatial distribution. The first-order complexity in aftershock distribution is clear and well beyond the influence of typical earthquake location uncertainty. The sequence included rupture of multiple fault segments, spatially separated. We use surface-wave based relative centroid locations to examine whether, at the small scale, the distribution of the aftershocks was influenced by location errors. Surface-wave based relative location has delineated precise oceanic transform fault earthquake locations in multiple regions. However, the relocated aftershocks off the coast of Sumatra seldom align along simple linear trends that are compatible with the corresponding fault strikes as estimated for the GCMT catalog. The relocation of roughly 60 moderate-earthquake epicentroids suggests that the faulting involved in the 2012 earthquake aftershock sequence included strain release along many short fault segments. Statistical analysis and temporal variations of aftershocks show a typical decay of the aftershocks but a relatively low number of aftershocks, as is common for intraplate oceanic earthquakes. Coulomb stress calculations indicate that most of the moderate-magnitude aftershocks are compatible with stress changes predicted by the large-event slip models. Furthermore the patterns in the aftershocks suggest that the formation of the boundary and eventual localization of deformation between the Indian and Australian plate is a complicated process.

2012 Sumatra Earthquake↗
DOE OSTI2020

Characterization of elastic mechanical properties of Tuscaloosa Marine Shale from well logs using the vertical transversely isotropic model

To avoid steep declines in the Tuscaloosa Marine Shale (TMS) production, wells are fracture-stimulated to release the hydrocarbons trapped in the matrix of the formation. An accurate estimation of Young’s modulus and Poisson’s ratio is essential for hydraulic fracture propagation. In addition, ignoring the highly heterogeneous and anisotropic character of TMS can lead to erroneous stress values, which subsequently affect hydraulic fracture width estimates and the overall hydraulic fracturing process. We have developed an empirical 1D geomechanical model that takes into account VTI anisotropy, and it is used to characterize the elastic mechanical properties of TMS in two wells. In the analyzed formation, the vertical Poisson’s ratio is less than the horizontal Poisson’s ratio, which suggests the necessity of an alternative to the ANNIE equations. The stiffness coefficients [Formula: see text] and [Formula: see text] were estimated using the relationships developed from the ultrasonic core data available for the two TMS. Further, correlations between the static and dynamic properties from laboratory tests were used to improve the minimum horizontal stress calculation. We compare VTI Young’s moduli, Poisson’s ratios, and minimum horizontal stress with the isotropic solution. VTI modeling improves the estimation of the elastic mechanical properties. The isotropic solution underestimates the minimum horizontal stress in the formation. Moreover, it was shown that the 20 ft shale interval below the TMS base is characterized by a low Young’s modulus (the vertical Young’s modulus is equal to 20 GPa, whereas the horizontal Young’s modulus is equal to 40 GPa) and may be a frac barrier.

Geochemistry & Geophysics↗
DOE OSTI2025

Inelastic deformation of diamond single crystals shock compressed to multimegabar stresses: Wave profile calculations

As the archetypal strong solid, the response of diamond shock compressed to multimegabar stresses is important for fundamental science and for numerical simulations of wave profiles for applications in high energy density physics experiments. Previous experiments and analysis have shown that the commonly used hydrodynamic assumption is invalid for diamond shock compressed to stresses below melt and an elastic–inelastic description is needed. Here, we present a phenomenological material model for calculating wave profiles in shock compressed diamond single crystals that incorporates this description. Also, to support the modeling effort, we carried out wave profile measurements on shock compressed diamond single crystals at the Sandia Z facility to augment previous measurements. Wave profiles for [100] and [111] diamond calculated using the material model provide a good match to the elastic–inelastic response (observed two-wave structure) measured at ∼325 and ∼360 GPa. Furthermore, the calculated peak stresses for single (overdriven) waves provide a good match to the measured Hugoniot states for stresses reaching ∼700 GPa, which is near melting conditions. The present results show that the diamond single crystal response at multimegabar shock stresses is characteristic of a brittle solid—pressure-dependent strength and strength loss due to inelastic deformation.

Deformation↗
DOE OSTI2022

Interlaboratory Reproducibility of Contour Method Data in a High Strength Aluminum Alloy

The contour method for residual stress measurement has seen significant development, but an experimental reproducibility study utilizing physical samples has not been published. A double-blind reproducibly study is reported, having scope beginning with EDM cutting and ending with residual stress calculation. A reinforced I-beam sample geometry is identified for its unique residual stress profile when extracted from residual stress bearing quenched aluminum bar (7050-T74). Contour measurements are prescribed on a midplane of symmetry with dimensions 24.0 mm by 50.0 mm. Fourteen identically prepared samples are fabricated from a single long bar with well characterized and uniform residual stress. Five samples throughout the bar are identified for planning measurements to validate sample uniformity and overall suitability of the residual stress field. The planning measurements employ a range of techniques: contour method, neutron diffraction, and hole-drilling. Eight samples are distributed to an international group of participants to execute their standard measurement practice. A double-blind process is followed to provide anonymity. Results are provided by eight participants: six being self-similar and two being quite different, the latter set aside as outliers. An average residual stress field is established from non-outlying results and the spatial distribution of reproducibility standard deviation is determined. The average stress field ranges from -60 to 70 MPa and the reproducibility standard deviation averages 8.1 MPa on the measurement plane. The average reproducibility standard deviation is about 3 × larger for points within 1.0 mm of plane boundaries (17.6 MPa) than for the remaining points (6.1 MPa). Reproducibility standard deviation (among different labs) for contour method residual stress measurement is found to be very similar to repeatability standard deviation (in a single lab) reported in prior work. The reproducibility observed here, for the entire measurement process, is also similar to that found in a prior reproducibility study limited to contour method data analysis.

36 MATERIALS SCIENCE↗
NASA NTRS1975

Fracture analysis of various cracked configurations in sheet and plate materials

A two-parameter fracture criterion was derived which relates the linear-elastic stress-intensity factor at failure, the elastic nominal failure stress, and two material parameters. The fracture criterion was used previously to analyze fracture data for surface- and through-cracked sheet and plate specimens under tensile loading. The fracture criterion was applied to center-crack tension, compact, and notch-band fracture specimens made of steel, titanium, or aluminum alloy materials tested at room temperature. The fracture data included a wide range of crack lengths, specimen widths, and thicknesses. The materials analyzed had a wide range of tensile properties. Failure stresses calculated using the criterion agreed well (plus or minus 10 percent) with experimental failure stresses. The criterion was also found to correlate fracture data from different specimen types (such as center-crack tension and compact specimens), within plus or minus 10 percent for the same material, thickness, and test temperature.

Newman, J. C., Jr.↗
NASA NTRS1976

Fracture analysis of various cracked configurations in sheet and plate materials

A two-parameter fracture criterion has been derived which relates the linear-elastic stress-intensity factor at failure, the elastic nominal failure stress, and two material parameters. The fracture criterion was applied to center-crack tension, compact, and notch-bend fracture specimens made of steel, titanium, or aluminum alloy materials tested at room temperature. The fracture data included a wide range of crack lengths, specimen widths, and thicknesses. The materials analyzed had a wide range of tensile properties. Failure stresses calculated using the criterion agreed well (plus or minus 10 percent) with experimental failure stresses. The criterion was also found to correlate fracture data from different specimen types (such as center-crack tension and compact specimens), within plus or minus 10 percent for the same material, thickness, and test temperature.

Newman, J. C., Jr.↗
NASA NTRS1994

Interlaminar tension strength of graphite/epoxy composite laminates

An L-shaped curved beam specimen and a tension loading fixture were proposed to measure the interlaminar tension strength of laminated and textile composites. The specimen size was 2 X 2 in. (51 X 51 mm). The use of a standard tension test machine and the introduction of load nearly at the specimen midthickness were the advantages of the proposed specimen. Modified Lekhnitskii and beam theory equations for calculating interlaminar stresses of an L-beam were verified by finite element analysis. The beam theory equation is simple and accurate for mean radius to thickness ratios greater than 1.5. The modified Lekhnitskii equations can be used for detailed stress field calculation. AS4/3501-6 graphite/epoxy unidirectional specimens with thicknesses of 16, 24, and 32 piles were fabricated and tested. The delamination initiation site agreed with the calculated maximum interlaminar tension stress location for all three thicknesses. Average interlaminar tension strengths of 16-, 24-, and 32-ply laminates were 47.6, 40.9, and 23.4 MPa, respectively.

Shivakumar, Kunigal N.↗
NASA NTRS1993

State of stress, faulting, and eruption characteristics of large volcanoes on Mars

The formation of a large volcano loads the underlying lithospheric plate and can lead to lithospheric flexure and faulting. In turn, lithospheric deformation affects the stress field beneath and within the volcanic edifice and can influence magma transport. Modeling the interaction of these processes is crucial to an understanding of the history of eruption characteristics and tectonic deformation of large volcanoes. We develop models of time-dependent stress and deformation for the Tharsis volcanoes on Mars. By means of a finite element code, we calculate stresses and displacements due to a volcano-shaped load emplaced on an elastic plate overlying a viscoelastic mantle. Models variously incorporate growth of the volcanic load with time and a detachment between volcano and lithosphere. The models illustrate the manner in which time-dependent stresses induced by lithospheric plate flexure beneath the volcanic load may affect eruption histories, and the derived stress fields can be related to tectonic features on and surrounding Martian volcanoes. As a result of flexure there are three regions where stresses become sufficiently large to cause failure by faulting, according to the Mohr-Coulomb criterion: at the surface of the plate just outward of the volcano, near the base of the elastic lithosphere beneath the center of the volcano, and on the upper flanks of the volcano early in its growth history.

Mcgovern, Patrick J.↗
NASA NTRS1993

Wrinkle ridges, reverse faulting, and the depth penetration of lithospheric stress in lunae planum, Mars

Tectonic features on a planetary surface are commonly used as constraints on models to determine the state of stress at the time the features formed. Quantitative global stress models applied to understand the formation of the Tharsis province on Mars constrained by observed tectonics have calculated stresses at the surface of a thin elastic shell and have neglected the role of vertical structure in influencing the predicted pattern of surface deformation. Wrinkle ridges in the Lunae Planum region of Mars form a conentric pattern of regularly spaced features in the eastern and southeastern part of Tharsis; they are formed due to compressional stresses related to the response of the Martian lithosphere to the Tharsis bulge. As observed in the exposures of valley walls in areas such as the Kasei Valles, the surface plains unit is underlain by an unconsolidated impact-generated megaregolith that grades with depth into structurally competent lithospheric basement. The ridges have alternatively been hypothesized to reflect deformation restricted to the surface plains unit ('thin skinned deformation') and deformation that includes the surface unit, megaregolith and basement lithosphere ('thick skinned deformation'). We have adopted a finite element approach to quantify the nature of deformation associated with the development of wrinkle ridges in a vertically stratified elastic lithosphere. We used the program TECTON, which contains a slippery node capability that allowed us to explicitly take into account the presence of reverse faults believed to be associated with the ridges. In this study we focused on the strain field in the vicinity of a single ridge when slip occurs along the fault. We considered two initial model geometries. In the first, the reverse fault was assumed to be in the surface plains unit, and in the second the initial fault was located in lithospheric basement, immediately beneath the weak megaregolith. We are interested in the conditions underwhich strain in the surface layer and basement either penetrates or fails to penetrate through the megaregolith. We thus address the conditions required for an initial basement fault to propagate through the megaregolith to the surface, as well as the effect of the megareolith on the strain tensor in the vicinity of a fault that nucleates in the surface plains unit.

Zuber, M. T.↗
DOE OSTI2017

Inclusion of Structural Flexibility in Design Load Analysis for Wave Energy Converters: Preprint

Hydroelastic interactions, caused by ocean wave loading on wave energy devices with deformable structures, are studied in the time domain. A midfidelity, hybrid modeling approach of rigid-body and flexible-body dynamics is developed and implemented in an open-source simulation tool for wave energy converters (WEC-Sim) to simulate the dynamic responses of wave energy converter component structural deformations under wave loading. A generalized coordinate system, including degrees of freedom associated with rigid bodies, structural modes, and constraints connecting multiple bodies, is utilized. A simplified method of calculating stress loads and sectional bending moments is implemented, with the purpose of sizing and designing wave energy converters. Results calculated using the method presented are verified with those of high-fidelity fluid-structure interaction simulations, as well as low-fidelity, frequency-domain, boundary element method analysis.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

WHOLESCALE: Coordinates of wells at San Emidio, Nevada

This dataset includes position coordinates and elevation information for wells at the WHOLESCALE San Emidio project location. Well positions in the attached file are characterized by UTM coordinates (Easting, Northing) in meters, and WHOLESCALE coordinates (Easting, Northing) relative in meters from a chosen reference point. The elevation of the top and bottom of open intervals within each well (representing locations of perforated or open-hole sections) are measured in meters (positive) above WGS84 geoid (mean sea level). The WHOLESCALE acronym stands for Water & Hole Observations Leverage Effective Stress Calculations and Lessen Expenses. The goal of the WHOLESCALE project is to simulate the spatial distribution and temporal evolution of stress in the geothermal system at San Emidio in Nevada, United States.

15 GEOTHERMAL ENERGY↗
DOE OSTI2016

WHOLESCALE: Mass Flux Rates for Wells at San Emidio in December 2016

This dataset provides mass flux rates in kg/s from six (production and injection) wells at San Emidio at minute intervals from December 1, 2016 - December 15, 2016. Files for injection wells are named with "IW", for instance "WellIW42-21SI.csv", and include negative flux rates. Files for production wells are named with "PW", for instance "WellPW61-21SI.csv", and include positive flux rates. The time interval studied includes the shutdown of pumping at the site, so flux rates are at or near zero around December 9. The WHOLESCALE acronym stands for Water & Hole Observations Leverage Effective Stress Calculations and Lessen Expenses. The goal of the WHOLESCALE project is to simulate the spatial distribution and temporal evolution of stress in the geothermal system at San Emidio in Nevada, United States.

15 GEOTHERMAL ENERGY↗
NASA NTRS1976

On the thermoelastic analysis of solar cell arrays and related material properties

Accurate prediction of failure of solar cell arrays requires accuracy in the computation of thermally induced stresses. This was accomplished by using the finite element technique. Improved procedures for stress calculation were introduced together with failure criteria capable of describing a wide range of ductile and brittle material behavior. The stress distribution and associated failure mechanisms in the N-interconnect junction of two solar cell designs were then studied. In such stress and failure analysis, it is essential to know the thermomechanical properties of the materials involved. Measurements were made of properties of materials suitable for the design of lightweight arrays: microsheet-0211 glass material for the solar cell filter, and Kapton-H, Kapton F, Teflon, Tedlar, and Mica Ply PG-402 for lightweight substrates. The temperature-dependence of the thermal coefficient of expansion for these materials was determined together with other properties such as the elastic moduli, Poisson's ratio, and the stress-strain behavior up to failure.

Salama, M. A.↗
DOE OSTI2024

Uniaxial compressive creep tests by spark plasma sintering of 70% theoretical density α -uranium and U-10Zr

Metallic fuels hold numerous advantages over conventional uranium dioxide fuels and are a key component of several liquid metal-cooled advanced reactor concepts including sodium fast reactors. These fuels undergo rapid swelling during early burnup; consequently, they spend most of their reactor lifetime in a porous state. The presence of this porosity alters many of the mechanical properties of the fuel including creep impacting fuel deformation during axial swelling. This work investigates the creep behavior of the porous fuel using a spark plasma sintering technique. Creep tests were performed for the first time on porous α-phase uranium and uranium with 10 wt. % zirconium (U-10Zr) samples. The samples of α-phase uranium and U-10Zr were fabricated from depleted uranium by spark plasma sintering and subjected to uniaxial compressive creep testing. Calculated stress exponents were found to be 2.6±1.6 and 5.7±1.4 for α-U and U-10Zr, respectively, and calculated activation energies were found to be 61.6±1.1kJ/mol for α-U. The creep data were also used to evaluate existing porosity inclusive in creep models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Coupled thermal–hydrological–mechanical modeling of CO 2 -enhanced coalbed methane recovery

CO 2 -enhanced coalbed methane recovery, also known as CO 2 -ECBM, is a potential win-win approach for enhanced methane production while simultaneously sequestering injected anthropogenic CO 2 to decrease CO 2 emissions into the atmosphere. Here, CO 2 -ECBM is simulated using a coupled thermal–hydrological–mechanical (THM) numerical model that considers multiphase (gas and water) flow and solubility, multicomponent (CO 2 and CH 4 ) diffusion and adsorption, heat transfer and coal deformation. The coupled model is based on the TOUGH-FLAC simulator, which is applied here for the first time to model CO 2 -ECBM. The capacity of the simulator for modeling methane production is verified by a code-to-code comparison with the general-purpose finite-element solver COMSOL. Then, the TOUGH-FLAC simulator is applied in an isothermal simulation to study the variations in permeability evolution during a CO 2 -ECBM operation while considering four different stress-dependent permeability models that have been implemented into the simulator. Finally, the TOUGH-FLAC simulator is applied in non-isothermal simulations to model THM responses during a CO 2 -ECBM operation.Our simulations show that the permeability evolution, mechanical stress, and deformation are all affected by changes in pressure, temperature and adsorption swelling, with adsorption swelling having the largest effect. The calculated stress changes do not induce any mechanical failure in the coal seam, except near the injection well in one case of a very unfavorable stress field.

01 COAL, LIGNITE, AND PEAT↗