Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stress field”

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

Completed Beltrami-Michell Formulation for Analyzing Radially Symmetrical Bodies

A force method formulation, the completed Beltrami-Michell formulation (CBMF), has been developed for analyzing boundary value problems in elastic continua. The CBMF is obtained by augmenting the classical Beltrami-Michell formulation with novel boundary compatibility conditions. It can analyze general elastic continua with stress, displacement, or mixed boundary conditions. The CBMF alleviates the limitations of the classical formulation, which can solve stress boundary value problems only. In this report, the CBMF is specialized for plates and shells. All equations of the CBMF, including the boundary compatibility conditions, are derived from the variational formulation of the integrated force method (IFM). These equations are defined only in terms of stresses. Their solution for kinematically stable elastic continua provides stress fields without any reference to displacements. In addition, a stress function formulation for plates and shells is developed by augmenting the classical Airy's formulation with boundary compatibility conditions expressed in terms of the stress function. The versatility of the CBMF and the augmented stress function formulation is demonstrated through analytical solutions of several mixed boundary value problems. The example problems include a composite circular plate and a composite circular cylindrical shell under the simultaneous actions of mechanical and thermal loads.

Kaljevic, Igor↗

Crack tip finite element analyses

Computer-oriented finite-element techniques for extracting the required information about the influence of cracks on stress fields are reviewed based on the knowledge that cracks induce stress singularities. Attention is focused on modeling behavior at a crack tip. Also discussed are crack tip analysis, utility of crack tip finite elements, and newer developments for cases where the singularity is not known a priori, i.e. the strength of the singularity is not established in advance of the computation. Available results are briefly summarized.

Swedlow, J. L.↗

Thermocryogenic buckling and stress analyses of a partially filled cryogenic tank subjected to cylindrical strip heating

Thermocryogenic buckling and stress analyses were conducted on a horizontally oriented cryogenic tank using the finite element method. The tank is a finite-length circular cylindrical shell with its two ends capped with hemispherical shells. The tank is subjected to cylindrical strip heating in the region above the liquid-cryogen fill level and to cryogenic cooling below the fill level (i.e., under thermocryogenic loading). The effects of cryogen fill level on the buckling temperature and thermocryogenic stress field were investigated in detail. Both the buckling temperature and stress magnitudes were relatively insensitive to the cryogen fill level. The buckling temperature, however, was quite sensitive to the radius-to-thickness ratio. A mechanical stress analysis of the tank also was conducted when the tank was under: (1) cryogen liquid pressure loading; (2) internal pressure loading; and (3) tank-wall inertia loading. Deformed shapes of the cryogenic tanks under different loading conditions were shown, and high-stress domains were mapped on the tank wall for the strain-gage installations. The accuracies of solutions from different finite element models were compared.

Ko, William L.↗

Coseismic fault slip associated with the 1992 M(sub w) 6.1 Joshua Tree, California, earthquake: Implications for the Joshua Tree-Landers earthquake sequence

Coseismic surface deformation associated with the M(sub w) 6.1, April 23, 1992, Joshua Tree earthquake is well represented by estimates of geodetic monument displacements at 20 locations independently derived from Global Positioning System and trilateration measurements. The rms signal to noise ratio for these inferred displacements is 1.8 with near-fault displacement estimates exceeding 40 mm. In order to determine the long-wavelength distribution of slip over the plane of rupture, a Tikhonov regularization operator is applied to these estimates which minimizes stress variability subject to purely right-lateral slip and zero surface slip constraints. The resulting slip distribution yields a geodetic moment estimate of 1.7 x 10(exp 18) N m with corresponding maximum slip around 0.8 m and compares well with independent and complementary information including seismic moment and source time function estimates and main shock and aftershock locations. From empirical Green's functions analyses, a rupture duration of 5 s is obtained which implies a rupture radius of 6-8 km. Most of the inferred slip lies to the north of the hypocenter, consistent with northward rupture propagation. Stress drop estimates are in the range of 2-4 MPa. In addition, predicted Coulomb stress increases correlate remarkably well with the distribution of aftershock hypocenters; most of the aftershocks occur in areas for which the mainshock rupture produced stress increases larger than about 0.1 MPa. In contrast, predicted stress changes are near zero at the hypocenter of the M(sub w) 7.3, June 28, 1992, Landers earthquake which nucleated about 20 km beyond the northernmost edge of the Joshua Tree rupture. Based on aftershock migrations and the predicted static stress field, we speculate that redistribution of Joshua Tree-induced stress perturbations played a role in the spatio-temporal development of the earth sequence culminating in the Landers event.

Bennett, Richard A.↗

Implementation of equivalent domain integral method in the two-dimensional analysis of mixed mode problems

An equivalent domain integral (EDI) method for calculating J-intergrals for two-dimensional cracked elastic bodies is presented. The details of the method and its implementation are presented for isoparametric elements. The total and product integrals consist of the sum of an area of domain integral and line integrals on the crack faces. The line integrals vanish only when the crack faces are traction free and the loading is either pure mode 1 or pure mode 2 or a combination of both with only the square-root singular term in the stress field. The EDI method gave accurate values of the J-integrals for two mode I and two mixed mode problems. Numerical studies showed that domains consisting of one layer of elements are sufficient to obtain accurate J-integral values. Two procedures for separating the individual modes from the domain integrals are presented. The procedure that uses the symmetric and antisymmetric components of the stress and displacement fields to calculate the individual modes gave accurate values of the integrals for all problems analyzed. The EDI method when applied to a problem of an interface crack in two different materials showed that the mode 1 and mode 2 components are domain dependent while the total integral is not. This behavior is caused by the presence of the oscillatory part of the singularity in bimaterial crack problems. The EDI method, thus, shows behavior similar to the virtual crack closure method for bimaterial problems.

Raju, I. S.↗

A workflow to assess the efficacy of brine extraction for managing injection-induced seismicity potential using data from a CO 2 injection site near Decatur, Illinois

Injection of CO 2 for storage in a geologic formation increases pore pressure and alters in situ stresses. Depending on the orientation of any existing fault and fracture planes, such as critically stressed planes, this stress alteration will modify normal stresses acting on planes and could result in frictional sliding and release stored energy in the form of seismicity. Brine extraction (BE) is a technique that can be applied prior to, or during, CO 2 injection to reduce pore pressure for increasing storage capacity and, potentially, for reducing the likelihood of frictional sliding. Here a workflow is described to assess the efficacy of BE for mitigating frictional sliding (i.e., seismicity) during injection and entails: site characterization, stress calculations and failure assessment, static and dynamic modeling, and BE operational planning. Site characterization describes the stress field used to calculate the Coulomb Failure Function (CFF) that constrains allowable pore pressure changes and injection rates in the numerical simulation of CO 2 injection scenarios. The inclusion of BE in the workflow allows for determination of the potential need for pressure reduction, and evaluation of the effectiveness of this operation. Example application of the workflow using an injection field dataset near Decatur, IL, provides insight on fracture planes and stresses at the site, formation properties and the impact of variable CO 2 injection-rate targets on whether BE plans are required. The study workflow indicates that BE could enhance CO 2 injection rate by 39% and correspondingly reduce the potential for injection-induced seismicity as indicated by a reduction in CFF.

58 GEOSCIENCES↗

Damage evolution in polymer due to exposure to high-pressure hydrogen gas

The use of hydrogen as a fuel is increasing exponentially, and the most economical way to store and transport hydrogen for fuel use is as a high-pressure gas. Polymers are widely used for hydrogen distribution and storage systems because they are chemically inert towards hydrogen. However, when exposed to high-pressure hydrogen, some hydrogen diffuses through polymers and occupies the preexisting cavities inside the material. Upon depressurization, the hydrogen trapped inside polymer cavities can cause blistering or cracking by expanding these cavities. A continuum mechanics--based deformation model was deployed to predict the stress distribution and damage propagation while the polymer undergoes depressurization after high-pressure hydrogen exposure. The effects of cavity size, cavity location, and pressure inside the cavity on damage initiation and evolution inside the polymer were studied. The stress and damage evolution in the presence of multiple cavities was also studied, because interaction among cavities alters the damage and stress field. It was found that all these factors significantly change the stress state in the polymer, resulting in different paths for damage propagation. The effect of adding carbon black filler particles and plasticizer on the damage was also studied. It was found that damage tolerance of the polymer increases drastically with the addition of carbon black fillers, but decreases with the addition of the plasticizer.

Kulkarni, Shank S.↗

Estimation of the orientation of stress in the Earth’s crust without earthquake or borehole data

Abstract Mechanical stress acting in the Earth’s crust is a fundamental property that is important for a wide range of scientific and engineering applications. The orientation of maximum horizontal compressive stress can be estimated by inverting earthquake source mechanisms and measured directly from borehole-based measurements, but large regions of the continents have few or no observations. Here we present an approach to determine the orientation of maximum horizontal compressive stress by measuring stress-induced anisotropy of nonlinear susceptibility, which is the derivative of elastic modulus with respect to strain. Laboratory and Earth experiments show that nonlinear susceptibility is azimuthally dependent in an anisotropic stress field and is maximum in the orientation of maximum horizontal compressive stress. We observe this behavior in the Earth—in Oklahoma and New Mexico, U.S.A, where maximum nonlinear susceptibility coincides with the orientation of maximum horizontal compressive stress measured using traditional methods. Our measurements use empirical Green’s functions and solid-earth tides and can be applied at different temporal and spatial scales.

58 GEOSCIENCES↗

Near-zero-field magnetoresistance measurements: A simple method to track atomic-scale defects involved in metal-oxide-semiconductor device reliability

We demonstrate the ability of a relatively new analytical technique, near-zero-field magnetoresistance (NZFMR), to track atomic-scale phenomena involved in the high-field stressing damage of fully processed Si metal-oxide-semiconductor field-effect transistors. We show that the technique is sensitive to both the Pb0 and Pb1 dangling bond centers and that the presence of both centers can be inferred through NZFMR via hyperfine interactions with the central 29Si atoms of the dangling bonds. The NZFMR results also provide evidence for the redistribution of mobile hydrogen atoms at the Si/SiO2 interface and also a potential change in the average dipolar coupling constant between electrons in neighboring defects. This work shows that NZFMR offers significant analytical power for studying technologically relevant semiconductor device reliability problems and has advantages in experimental simplicity over comparable techniques.

Moxim, Stephen J. (ORCID:0000000309659904)↗

Measuring the residual stress and stress-corrosion cracking susceptibility of additively manufactured 316L by ASTM G36-94

Residual stress is a contributor to stress-corrosion cracking (SCC) and a common byproduct of additive manufacturing (AM). In this work, the relationship between residual stress and SCC susceptibility in laser powder bed fusion AM 316L stainless steel was studied through immersion in saturated boiling magnesium chloride per ASTM G36-94. The residual stress was varied by changing the sample height for the as-built condition and additionally by heat treatments at 600, 800, and 1200 °C to control, and in some cases reduce, residual stress. In general, all samples in the as-built condition showed susceptibility to SCC with the thinner, lower residual stress samples showing shallower cracks and crack propagation occurring perpendicular to melt tracks due to local residual stress fields. The heat-treated samples showed a reduction in residual stress for the 800 and 1200 °C samples. Both were free of cracks after >300 hours of immersion in MgCl 2 , while the 600 °C sample showed similar cracking to their as-built counterpart. Geometrically necessary dislocation (GND) density analysis indicates that the dislocation density may play a major role in the SCC susceptibility.

36 MATERIALS SCIENCE↗

Influence of transverse shear on an axial crack in a cylindrical shell

An axial crack in a cylindrical shell is investigated by use of a 10th order shell theory, which accounts for transverse shear deformations as well as a special kind of orthotropy. The symmetric problem is formulated in terms of two coupled singular integral equations, which are solved numerically. The asymptotic membrane and bending stress fields ahead of the crack are found to be self similar. Stress intensity factors are given as a function of the shell parameter for various values of the ratio crack length to shell thickness. considerable differences from 8th order shell theory results are found for the bending stresses, while the membrane stresses of the 8th order theory seem to be a lower limit reached for very thin shells.

Krenk, S.↗

Influence of transverse shear on an axial crack in a cylindrical shell

An axial crack in a cylindrical shell is investigated by use of a 10th order shell theory, which accounts for transverse shear deformations as well as a special kind of orthotropy. The symmetric problem is formulated in terms of two coupled singular integral equations, which are solved numerically. The asymptotic membrane and bending stress fields ahead of the crack are found to be self similar. Stress intensity factors are given as a function of the shell parameter for various values of the ratio crack length to shell thickness. Considerable differences from 8th order shell theory results are found for the bending stresses, while the membrane stresses of the 8th order theory seems to be a lower limit reached for very thin shells.

Krenk, S.↗

Transverse shear effects on the stress-intensity factor for a circumferentially cracked, specially orthotropic cylindrical shell

The problem of a cylindrical shell containing a circumferential through crack is considered by taking into account the effect of transverse shear deformations. The formulation is given for a specially orthotropic material within the confines of a linearized shallow shell theory. The particular theory used permits the consideration of all five boundary conditions regarding moment and stress resultants on the crack surface. Consequently, aside from multiplicative constants representing the stress intensity factors, the membrane and bending components of the asymptotic stress fields near the crack tip are found to be identical. The stress intensity factors are calculated separately for a cylinder under a uniform membrane load, and that under a uniform bending moment. Sample results showing the nature of the out-of-plane crack surface displacement and the effect of the Poisson's ratio are presented.

Delale, F.↗

Influence of transverse shear on an axial crack in a cylindrical shell

An axial crack in a cylindrical shell is investigated by use of a 10th order shell theory, which accounts for transverse shear deformations as well as a special kind of orthotropy. The symmetric problem is formulated in terms of two coupled singular integral equations, which are solved numerically. The asymptotic membrane and bending stress fields ahead of the crack are found to be self similar. Stress intensity factors are given as a function of the shell parameter for various values of the ratio shell radius to shell thickness. Considerable differences from 8th order shell theory results are found for the bending stresses, while the membrane stresses of the 8th order theory seems to be a lower limit reached for very thin shells.

Krenk, S.↗

On the Application of Deformation Kinetics to Nonlinear Constitutive Relations at Higher Temperatures

A single phenomenological constitutive equation is derived theoretically from first principles and applied to aluminum, tin and lead. The theory is based on deformation kinetics of steady creep in which the fundamental mechanism is atomic transport over potential barriers whose conformation is distorted by the application of a stress field. The form of the functional dependence of barrier distortion and stress over the entire temperature range is found to be a sigmoidal curve which tends to straight lines of a unit slope in the small and high stress regions. With this form of barrier distortion, the constitutive equation prediction the steady creep behavior of aluminum, tin and lead over a wide range of temperature and stress.

Valanis, K. C.↗

On the numerical performance of three-dimensional thick shell elements using a hybrid/mixed formulation

Three-dimensional thick shell elements with 8, 16, and 18 nodes are formulated by using the hybrid/mixed method. In bending applications, these elements are free from locking effect and give improved stress predictions. Finite element equations are derived from the Hellinger-Reissner variational principle in which both the displacement and stress fields are approximated by independent interpolation functions. For the assumption of stress parameters, three guidelines are followed: (1) suppression of kinematic deformation modes, (2) invariant element property, and (3) the constraint index exhibited by the element, when applied to constrained-media problems, must be greater than or equal to one. Numerical results are presented to show the element's behavior characteristics regarding sensitivity to locking, distortion effect (patch tests), mesh convergence and the accuracy of stress evaluation.

Graf, W.↗

Resolved shear stress intensity coefficient and fatigue crack growth in large crystals

Fatigue crack growth tests were carried out on large-grain Al 7029 aluminum alloy and the finite element method was used to calculate the stress field near the tip of a zigzag crack. The resolved shear stresses on all 12 slip systems were computed, and the resolved shear stress intensity coefficient (RSSIC) was defined. The RSSIC was used to analyze the irregular crack path and was correlated with the rate of single-slip-plane shear crack growth. Fatigue crack growth was found to be caused primarily by shear decohesion at a crack tip. When the RSSIC on a single-slip system was much larger than all the others, the crack followed a single-slip plane. When the RSSICs on two conjugate slip systems were comparable, a crack grew in a zigzag manner on these planes and the macrocrack-plane bisected the two active slip planes. The maximum RSSIC on the most active slip system is proposed as a parameter to correlate with the shear fatigue crack growth rate in large crystals.

Chen, Q.↗