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 91 records · Page 5

A 3D Full Stress Tensor Model for Oklahoma

Abstract The stress tensor is an important property for upper crustal studies such as those that involve pore fluids and earthquake hazards. At tectonic plate scale, plate boundary forces and mantle convection are the primary drivers of the stress field. In many local settings (10–100 s of km and <10 km depth) in tectonic plate interiors, we can simplify by assuming a constant background stress field that is perturbed by local heterogeneity in density and elasticity. Local stress orientation and sometimes magnitude can be estimated from earthquake and borehole‐based observations when available. Modeling of the local stress field often involves interpolating sparse observations. We present a new method to estimate the 3D stress field in the upper crust and demonstrate it for Oklahoma. We created a 3D material model by inverting multiple types of geophysical observations simultaneously. Integrating surface‐wave dispersion, local travel times and gravity observations produces a model of P‐wave velocity, S‐wave velocity, and density. The stress field can then be modeled using finite element simulations. The simulations are performed using our simplified view of the local stress field as the sum of a constant background stress field that is perturbed by local density and elasticity heterogeneity and gravitational body forces. An orientation of N82°E, for the maximum compressive tectonic force, best agrees with previously observed stress orientations and faulting types in Oklahoma. The gravitational contribution of the horizontal stress field has a magnitude comparable to the tectonic contribution for the upper 5 km of the subsurface.

58 GEOSCIENCES↗

Nondestructive Evaluation of Stress Corrosion Cracking in a Welded Steel Plate Using Guided Ultrasonic Waves

Stress corrosion cracking (SCC) had occurred in early-generation high-level nuclear waste tanks constructed by welding carbon steel. This paper describes an ultrasonic inspection system and its fundamental ability to detect and quantify the length of SCC on thick welded steel plates. The finite element method (FEM) was applied to simulate the welding process to estimate the welding residual stress field. Growth of stress corrosion cracks is driven by crack stress intensities exceeding the subcritical cracking threshold intensity. The subject plate was experimentally inspected with ultrasonic nondestructive evaluation (NDE) techniques to characterize the extent of SCC. The NDE system uses a piezoelectric transducer to generate guided waves in the thick steel plate, and a scanning laser Doppler vibrometer (SLDV) to measure multidimensional time–space wavefield data over a user-defined scanning area in the plate surface. The measured wavefield data can show wave interactions in a localized area in the plate due to the presence of the discontinuities of the SCC. To generate an inspection image that can precisely show the crack’s location and/or the dimension, the wavefield data are further processed to generate inspection image that maps the entire sample plate so the crack can be clearly identified in the plate while its length can be readily estimated. In conclusion, the ultrasonic test results for crack length agree well with the visually estimated length and are close to that predicted by the FEM for cracks in the weld residual stress field.

36 MATERIALS SCIENCE↗

Identifying multiple synergistic factors on the susceptibility to stress relaxation cracking in variously heat-treated weldments

The 347H austenitic stainless steel has been widely used for pressure vessels and pipeline (PVP) applications due to its excellent creep and corrosion resistance, which fit ideally to the harsh conditions in petrochemical industries, fossil fuel or nuclear power plants, and modern energy storages. However, a failure mode has been commonly observed with cracks emerging at the heat affected zone (HAZ) of weldments during post-weld heat treatment (PWHT) or under intermediate to high temperature service conditions. This phenomenon is termed as Stress Relaxation Cracking (SRC) since the purpose of PWHT is to relieve the welding-induced residual stress fields, or as Stress Age Cracking (SAC) if failure happens during service. A leading literature explanation of this failure suggests that the residual stress relaxation and the precipitation dissolution and/or re-precipitation occur in the same temperature range, which can lead to locally high strains and thus to crack at the grain boundaries. Since in situ spatial measurements of residual stress fields, microstructural evolution, and failure processes are nearly infeasible, this work recourses to a micromechanical finite element framework that models the high temperature failure as the nucleation and growth of grain boundary cavities, whereas various parameters such as thermomechanical loading history and its evolution, the competition of grain-interior dislocation creep and grain-boundary diffusion in failure lifetime, and microstructural heterogeneities (such as the precipitate free zone near grain boundaries) can be quantitatively incorporated. It can be concluded from these microstructure-explicit simulations that an accurate knowledge of residual stress evolution and a carefully calibrated set of material constitutive parameters are the essential prerequisites for lifetime predictions. The understanding of individual governing factors also leads to a mechanistic interpretation of the observed SRC susceptibility C-curves. In conclusion, these results suggest that the criticality of residual stress evolution, but not the precipitation-induced local strains, be the leading factor for SRC.

347H stainless steel weldments↗

Statistics of internal stress fluctuations in dislocated crystals and relevance to density-based dislocation dynamics models

A statistical analysis of internal stress fluctuations, defined as the difference between the local mean stress and stress on dislocations, is presented for deforming crystals with 3D discrete dislocation systems. Dislocation realizations are generated using dislocation dynamics simulations and the associated stress field is computed as a superposition of a regularized stress field of dislocation lines within the domain of the solution and a complementary stress field computed via a finite-element boundary value problem. The internal stress fluctuations of interest are defined by an ensemble of the difference between the stress on dislocation lines and the local mean field stress in the crystal. The latter is established in a piecewise fashion over small voxels in the crystal thus allowing the difference between the local average stress and stress on segments to be easily estimated. The results show that the Schmid stress (resolved shear stress) and Escaig stress fluctuations on various slip systems sampled over a random set of points follow a Cauchy (Lorentz) distribution at all strain levels, with the amplitude and width of the distribution being dependent on the strain. Finally, the implications of the Schmid and Escaig internal stress fluctuations are discussed from the points of view of dislocation cross-slip and the dislocation motion in continuum dislocation dynamics.

36 MATERIALS SCIENCE↗

Stress Recovery and Error Estimation for 3-D Shell Structures

The C1-continuous stress fields obtained from finite element analyses are in general lower- order accurate than are the corresponding displacement fields. Much effort has focussed on increasing their accuracy and/or their continuity, both for improved stress prediction and especially error estimation. A previous project developed a penalized, discrete least squares variational procedure that increases the accuracy and continuity of the stress field. The variational problem is solved by a post-processing, 'finite-element-type' analysis to recover a smooth, more accurate, C1-continuous stress field given the 'raw' finite element stresses. This analysis has been named the SEA/PDLS. The recovered stress field can be used in a posteriori error estimators, such as the Zienkiewicz-Zhu error estimator or equilibrium error estimators. The procedure was well-developed for the two-dimensional (plane) case involving low-order finite elements. It has been demonstrated that, if optimal finite element stresses are used for the post-processing, the recovered stress field is globally superconvergent. Extension of this work to three dimensional solids is straightforward. Attachment: Stress recovery and error estimation for shell structure (abstract only). A 4-node, shear-deformable flat shell element developed via explicit Kirchhoff constraints (abstract only). A novel four-node quadrilateral smoothing element for stress enhancement and error estimation (abstract only).

Riggs, H. R.↗

Thermoelastic stress in oceanic lithosphere due to hotspot reheating

The effect of hotspot reheating on the intraplate stress field is investigated by modeling the three-dimensional thermal stress field produced by nonuniform temperature changes in an elastic plate. Temperature perturbations are calculated assuming that the lithosphere is heated by a source in the lower part of the thermal lithosphere. A thermal stress model for the elastic lithosphere is calculated by superposing the stress fields resulting from temperature changes in small individual elements. The stress in an elastic plate resulting from a temperature change in each small element is expressed as an infinite series, wherein each term is a source or an image modified from a closed-from half-space solution. The thermal stress solution is applied to midplate swells in oceanic lithosphere with various thermal structures and plate velocities. The results predict a stress field with a maximum deviatoric stress on the order of 100 MPa covering a broad area around the hotspot plume. The predicted principal stress orientations show a complicated geographical pattern, with horizontal extension perpendicular to the hotspot track at shallow depths and compression along the track near the bottom of the elastic lithosphere.

Zhu, Anning↗

Implementation of J-A Methodology Elastic-Plastic Crack Instability Analysis Capability into the WARP-3D Code

Characterization of the near crack-tip stress/strain fields is the foundation of fracture mechanics. The description of the near tip stress field and the prediction of when fracture occurs is well established for brittle materials that exhibit linear elastic behavior. However, in ductile materials or conditions that violate linear elastic assumptions (Aluminum alloys, Al 2024-T3, Al 2024- T351 etc.), the elastic-plastic crack-tip stress fields are characterized by the Hutchison-Rice-Rosengren (HRR) field. The J-integral is commonly used to characterize amplitude of the HRR field under elastic-plastic conditions. The J-integral has been demonstrated for crack-tip fields that are under high constraint conditions (i.e., small-scale plasticity where the J-dominance is maintained). However, as the external load increases, yielding changes from small- to largescale plasticity and usually a loss of constraint (i.e., reduction in the triaxial stress field along the crack front). The loss of constraint leads to the deviation of the crack-tip stress fields from that given by the HRR field. Hence, the J-dominance will be gradually lost and additional parameter(s) are required to quantify the crack-tip stress fields and predict fracture behavior. The assessment objectives were to: 1) implement a two-parameter (i.e., J-A) fracture criterion into an elastic-plastic three-dimensional (3D) finite element analysis (FEA), 2) validate the implementation by comparison with the A parameter from literature data, 3) conduct material characterization tests to quantify the material behavior and provide fracture data for validation of the J-A fracture criteria, and (4) perform evaluations to establish if the J-A criteria can be used to predict fracture in a ductile metallic material (e.g., aluminum alloys). The A parameter in these criteria is the second parameter in a three-term elastic-plastic asymptotic expansion of the neartip stress behavior. A series of extensive FEAs were performed using WARP3D software package to obtain solutions for the A parameter for different specimen configurations. The methodology needed for the estimation of the A parameter in the asymptotic expansion was developed and implemented using Matlab®. A user material (UMAT) routine was used to model the material stress-strain response using a Ramberg-Osgood power law with a hardening exponent (n) and a material coefficient (alpha). This UMAT routine was successfully implemented in WARP3D software and validated through comparison with the experimental data. Three configurations were extracted from published results: 1) center cracked plate (CCP), 2) single edge-cracked plate (SECP), and 3) double edge-cracked plate (DECP). These configurations and four other configurations (three-hole tension (THT)), three-point bend (3PTB), three-hole compact tension (3PCT), and compact tension (CT)) were analyzed to verify the methodology that was developed and implemented into WARP3D. Solutions of the A parameter were obtained for remote tension loading conditions that started with small-scale yielding and continued into the large-scale plasticity regime. The results indicate that the methodology developed can be used to calculate the elastic-plastic J-A parameters for test specimens with a range of crack geometries, material strain hardening behaviors, and loading conditions. The J-A parameters were implemented as fracture criteria and used to predict the test results. For comparison, other fracture criteria were used to predict the same test results. Major findings include: The A constraint parameter A varies with specimen type and applied load thus accurate determination is crucial in predicting the failure load, and the A parameter is asymptotic as the failure load is approached, making an accurate determination difficult (i.e., small differences in the A parameter can cause large variations in failure load) for materials exhibiting elastic-plastic behavior. The failure predictions from J-A methodology were more accurate than the traditionally used KC and J methods, and have comparable scatter to that observed when using the crack-tip opening angle (CTOA) method. However, the J-A methodology requires considerable effort (expertise level and labor) to implement and to evaluate the A parameter for different specimen types and materials, or to apply this methodology to part-through crack (e.g., 3D problems) structural applications.

Hamm, Kenneth R., Jr.↗

Variations of Stress Parameters in the Southern California Plate Boundary Around the South Central Transverse Ranges

We examine the stress field in Southern California with a focus on the region near the South Central Transverse Ranges (SCTR), using a refined stress inversion methodology. We independently invert declustered and aftershocks focal mechanisms from the years 1981 to 2017. Comparisons between the stress field orientations and the stress ratios from these inversions provide information on local sources of the stress field. Regionally, the SHmax tends to be toward NNE and the stress ratios vary from transtensional stress regime near the Eastern California Shear Zone (ECSZ), to strike-slip faulting near the SCTR, and toward transpression near the Western Transverse Ranges. Detailed analysis of stress parameters near the SCTR indicates deviations from the regional strike-slip faulting including the following: (1) The San Bernardino Mountain area shows S Hmax direction toward NNW and transpressional stress components likely associated with the motions of the San Andreas Fault and ECSZ; (2) the Cajon Pass and San Gorgonio Pass show transpressional stress regime near the bottom of the seismogenic zone likely associated with the elevated topography; and (3) in Crafton Hills, rotation of the principal stress plunges and SHmax direction and transtensional stress regime below ~10 km, along with lower estimated apparent friction coefficient, suggest a weak fault possibly associated with deep creep. The performed multiscale analysis resolves effects of regional and local loadings. The average stress parameters in the analyzed 37 years do not show significant temporal variations of stress field near the SCTR region.

58 GEOSCIENCES↗

Magnetic stress anisotropy field in plated cylindrical Permalloy films.

An analysis is made of the magnetic stress anisotropy field (Hks) arising from internal and external stress sources in plated-wire memory elements. The analysis takes into consideration circumferential composition variation and cylindrical geometry of the Permalloy film. Expressions are derived relating Hks to uniaxial film stress, average composition, and amplitude of composition variation. A result of particular importance is that even for average zeromagnetostrictive composition (ZMC) films, Hks may still make an appreciable contribution to the total anisotropy field if the composition is not uniform. Calculated Hks characteristics are shown to correlate with anisotropy field changes observed in annealing experiments. Examples are given to show the importance of composition uniformity in determining the stability of the anisotropy field. The utility of the analysis is extended by the inclusion of data expressing the inverse relation between anisotropy field and easy-axis dispersion in the film.

Lutes, O. S.↗

A variational method for finite element stress recovery: Applications in one-dimension

It is well-known that stresses (and strains) calculated by a displacement-based finite element analysis are generally not as accurate as the displacements. In addition, the calculated stress field is typically discontinuous at element interfaces. Because the stresses are typically of more interest than the displacements, several procedures have been proposed to obtain a smooth stress field, given the finite element stresses, and to improve the accuracy. Hinton and Irons introduced global least squares smoothing of discrete data defined on a plane using a finite element formulation. Tessler and co-workers recently developed a conceptually similar formulation for smoothing of two-dimensional data based on a discrete least square approximation with a penalty constraint. The penalty constraint results in a stress field which is C(exp 1)-continuous, a result not previously obtained. The approach requires additional, 'smoothing' finite element analysis and for their two-dimensional application, they used a conforming C(exp 0)-continuous triangular finite element based on a conforming plate element. This paper presents the results of a detailed investigation into the application of Tessler's smoothing procedure to the smoothing of finite element stresses from one-dimensional problems. Although the one-dimensional formulation has some practical applicability, such as in truss, beam, axisymmetric mechanics, and one-dimensional heat conduction, the primary motivation for developing the one-dimensional smoothing case is to explore the characteristics of the general smoothing strategy. In particular, it is used to describe the behavior of the method and to explore the suitability of criteria proposed for the smoothing analysis. Prior to presenting numerical results, the variational formulation of the smoothing strategy is presented and a criterion for the smoothing analysis is described.

Riggs, H. Ronald↗

Wind stress over the Arabian Sea from ship reports and Seasat scatterometer data

Seasat scatterometer data over the Arabian Sea are used to build wind-stress fields during July and August 1978. They are first compared with 3-day wind analyses from ship data along the Somali coast. Seasat scatterometer specifications of 2-m/s and 20-deg accuracy are fulfilled in almost all cases. The exceptions are for winds stronger than 14 m/s, which are underestimated by the scatterometer by 15 percent. Wind stress is derived from these wind data using a bulk formula with a drag coefficient depending on the wind intensity. A successive-correction objective analysis is used to build the wind-stress field over the Arabian Sea with 2 x 2-deg and 6-day resolution. The final wind-stress fields are not significantly dependent on the objective analysis because of the dense coverage of the scatterometer. The combination of scatterometer and coastal ship data gives the best coverage to resolve monsoon wind structures even close to the coast. The final wind stress fields show wind features consistent with other monthly mean wind stress field. However, a high variability is observed on the 6-day time scale.

Perigaud, C.↗

$S$ Hmax orientation in the Alpine region from observations of stress-induced anisotropy of nonlinear elasticity

The orientation of $S$ Hmax is commonly estimated from in situ borehole breakouts and earthquake focal mechanisms. Borehole measurements are expensive, and therefore sparse, and earthquake measurements can only be made in regions with many well-characterized earthquakes. Here, we derive the stress-field orientation using stress-induced anisotropy in nonlinear elasticity. In this method, we measure the strain derivative of velocity as a function of azimuth. We use a natural pump-probe (NPP) approach which consists of measuring elastic wave speed using empirical Green’s functions (probe) at different points of the earth tidal strain cycle (pump). The approach is validated using a larger data set in the Northern Alpine Foreland region where the orientation of maximum horizontal compressive stress is known from borehole breakouts and drilling-induced fractures. The technique resolves NNW-SSW to N-S directed $S$ Hmax which is in good agreement with conventional methods and the recent crustal stress model. We confirm that the NPP method can be applied to dense large-scale seismic arrays. The technique is then applied to the Southern Alps to understand the contemporary stress pattern associated with the ongoing deformation due to counterclockwise rotation of the Adriatic plate with respect to the European plate. Our results explain why the two major faults in Northeastern Italy, the Giudicarie Fault and the Periadriatic Line (Pustertal–Gailtal Fault) are currently inactive, while the currently acting stress field allows faults in Slovenia to deform actively. We have demonstrated that the pump-probe method has the potential to fill in the measurement gap left by conventional approaches, both in terms of regional coverage and in depth.

58 GEOSCIENCES↗

A Comparison of Radiation-Induced and High-Field Electrically Stress-Induced Interface Defects in Si/SiO 2 MOSFETs via Electrically Detected Magnetic Resonance

Here, we utilize electrically detected magnetic resonance (EDMR) measurements to compare high-field stressed, and gamma irradiated Si/SiO 2 metal–oxide–silicon (MOS) structures. We utilize spin-dependent recombination (SDR) EDMR detected using the Fitzgerald and Grove dc I-V approach to compare the effects of high-field electrical stressing and gamma irradiation on defect formation at and near the Si/SiO 2 interface. As anticipated, both greatly increase the concentration of P b centers (silicon dangling bonds at the interface) densities. The irradiation also generated a significant increase in the dc I-V EDMR response of E' centers (oxygen vacancies in the SiO 2 films), whereas the generation of an E' EDMR response in high-field stressing is much weaker than in the gamma irradiation case. These results likely suggest a difference in their physical distribution resulting from radiation damage and high electric field stressing.

42 ENGINEERING↗

Weld stresses beyond elastic limit: Materials discontinuity

When welded structures depend on properties beyond the elastic limit to qualify their ultimate safety factor, and weld-parent materials abruptly change at the interface, then stress discontinuity is inevitable. The stress concentration is mildly sensitive to material relative strain hardening and acutely sensitive to applied stress fields. Peak stresses occur on the weld surface, at the interface, and dissipate within a 0.01-inch band. When the stress is intense, the weld will always fracture at the interface. The analysis incorporates a classical mechanics model to more sharply define stress spikes within the bandwidth, and suggests a relative material index and Poisson's ratio related to strain hardening. Implications are discussed which are applicable to industries of high performance structures.

Verderaime, V.↗

Effect of contact stresses in four-point bend testing of graphite/epoxy and graphite/PMR-15 composite beams

The results of in-plane four-point bend experiments on unidirectionally reinforced composite beams are presented for graphite/epoxy (T300/934) and graphite/polyimide (G30-500/PMR-15) composites. The maximum load and the location of cracks formed during failure were measured for testpieces with fibers oriented at various angles to the beam axis. Since most of the beams failed near one or more of the load points, the strength of the beams was evaluated in terms of a proposed model for the local stress distribution. In this model, an exact solution to the problem of a localized contact force acting on a unidirectionally reinforced half plane is used to describe the local stress field. The stress singularity at the load points is treated in a manner similar to the stress singularity at a crack tip in fracture mechanisms problems. Using this approach, the effect of fiber angle and elastic material properties on the strength of the beam is described in terms of a load intensity factor. For fiber angles less than 45 deg from the beam axis, a single crack is initiated near one of the load points at a critical value of the load intensity factor. The critical load intensity factor decreases with increasing fiber angle. For larger fiber angles, multiple cracks occur at locations both near and away from the load points, and the load intensity factor at failure increases sharply with increasing fiber angle.

Binienda, W. K.↗

The effect of contact stresses in four-point bend testing of graphite/epoxy and graphite/PMR-15 composite beams

The results of in-plane four-point bend experiments on unidirectionally reinforced composite beams are presented for graphite/epoxy (T300/934) and graphite/polyimide (G30-500/PMR-15) composites. The maximum load and the location of cracks formed during failure were measured for testpieces with fibers oriented at various angles to the beam axis. Since most of the beams failed near one or more of the load points, the strength of the beams was evaluated in terms of a proposed model, for the local stress distribution. In this model, an exact solution to the problem of a localized contact force acting on a unidirectionally reinforced half plane is used to describe the local stress field. The stress singularity at the load points is treated in a manner similar to the stress singularity at a crack tip in fracture mechanisms problems. Using this approach, the effect of fiber angle and elastic material properties on the strength of the beam is described in terms of a load intensity factor. For fiber angles less than 45 deg from the beam axis, a single crack is initiated near one of the load points at a critical value of the load intensity factor. The critical load intensity factor decreases with the increasing fiber angle. For larger fiber angles, multiple cracks occur at locations both near and away from the load points, and the load intensity factor at failure increases sharply with increasing fiber angle.

Binienda, Wieslaw K.↗

Interlaboratory Reproducibility of Contour Method Data Analysis and Residual Stress Calculation

While the contour method for residual stress assessment has developed rapidly, no published study documents its interlaboratory reproducibility. Objective: Here we report an initial reproducibility experiment focused on contour method data analysis and residual stress calculation. The experiment uses surface topography data from a physical process simulation of elastic-plastic beam bending. The simulation provides surface topography, for input to the contour method data analysis, as well as a known residual stress field with 130 MPa peak magnitude. Additionally, to increase realism, noise and specific artifacts are added to the topography data. A group of participants received the topography data (without the known residual stress), independently analyzed the data, and submitted results as a two-dimensional residual stress field. Analysis of submissions provides a group average residual stress field and the spatial distribution of reproducibility standard deviation. The group average residual stress agrees with the known stress in magnitude and spatial trend. The reproducibility standard deviation ranges from 2 to 54 MPa over the measurement plane, with an average of 5.4 MPa. Reproducibility standard deviation is smaller in the cross-section interior (≤ 5 MPa), modest near local extrema in the stress field (5 to 10 MPa), and larger near the cross-section boundaries (10 to 30 MPa). Overall, the largest values of reproducibility standard deviation (up to 54 MPa) occur in limited areas where artifacts had been added to the topography data; while some participants identified and removed these artifacts, some did not, leading to systematic differences that elevated the standard deviation.

42 ENGINEERING↗

Three Dimensional Numerical Simulation and Characterization of Crack Growth in the Weld Region of a Friction Stir Welded Structure

Friction stir welding (FSW) fabrication technology is being adopted in aerospace applications. The use of this technology can reduce production cost, lead-times, reduce structural weight and need for fasteners and lap joints, which are typically the primary locations of crack initiation and multi-site fatigue damage in aerospace structures. FSW is a solid state welding process that is well-suited for joining aluminum alloy components; however, the process introduces residual stresses (both tensile and compressive) in joined components. The propagation of fatigue cracks in a residual stress field and the resulting redistribution of the residual stress field and its effect on crack closure have to be estimated. To insure the safe insertion of complex integral structures, an accurate understanding of the fatigue crack growth behavior and the complex crack path process must be understood. A life prediction methodology for fatigue crack growth through the weld under the influence of residual stresses in aluminum alloy structures fabricated using FSW will be detailed. The effects and significance of the magnitude of residual stress at a crack tip on the estimated crack tip driving force are highlighted. The location of the crack tip relative to the FSW and the effect of microstructure on fatigue crack growth are considered. A damage tolerant life prediction methodology accounting for microstructural variation in the weld zone and residual stress field will lead to the design of lighter and more reliable aerospace structures

Seshadri, Banavara R.↗