Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stress Modeling”

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 1,027 records · Page 57

An Integrated Computational Materials Engineering (ICME) Approach to Design Nonlinear Transition Zones Between Dissimilar Metals

Current approaches to designing graded transition joints (GTJs) between dissimilar metals often rely on linear changes in both composition profiles and thickness of each sublayer. This increases fabrication cost and may not be optimal with respect to residual stress or the formation of undesirable phases. Here, in this study, GTJs between P91 ferritic/martensitic steel and 347H austenitic stainless steel were designed using Integrated Computational Materials Engineering (ICME) principles with nonlinear composition and length profiles. Guided by inputs from classical mechanics and CALPHAD predictions of carbon chemical potential, a novel transition zone consisting of five discrete compositions was proposed, with the thickness of each sublayer varying according to a brachistochrone-inspired distribution. In addition to carbon potential gradients, CALPHAD was used to predict coefficients of thermal expansion, which were incorporated into finite element models to evaluate stress evolution. The proposed nonlinear design resulted in a smoother carbon potential gradient, lower carbon depletion at the P91 interface, and a comparable residual stress under long-term thermal exposure, compared to a conventional linear design using ten sublayers with equal thickness. This work introduces a brachistochrone-inspired distribution for GTJ design, offering a general framework for optimizing graded interfaces between dissimilar metals.

Directed Energy Deposition↗

Peridynamics and surrogate modeling of pressure-driven well stimulation

In this work we use the peridynamics theory of solid mechanics to simulate fracture in an annular rock domain subject to an in-situ stress and create surrogate models that predict the area of the resulting cracks. Peridynamics is a non-local formulation of continuum mechanics that naturally accommodates material discontinuities. Furthermore, unlike other fracture modeling techniques there is no need to provide information about the crack path. We utilize the peridynamics code Peridigm and take a two-stage approach to fracture modeling. First an implicit solve is performed to compute the in-situ stress state. We then execute an explicit solve where a pressure loading designed to emulate fluid-driven hydraulic fracture is applied at the borehole and transmitted to the pre-stressed rock. We present results from polynomial and single and multi-level Gaussian process surrogate models constructed from a sampling study of the peridynamics model. The surrogates predict crack area given a measure of the in-situ stress anisotropy and rise time and amplitude of the pressure loading. These surrogates take a minuscule fraction of peridynamics model's running time to evaluate and are a step towards enabling advanced optimization and uncertainty quantification workflows that require many model evaluations.

42 ENGINEERING↗

Evaluation of subgrid-scale models using an accurately simulated turbulent flow

A calculation of periodic homogeneous isotropic turbulence is used to simulate the experimental decay of grid turbulence. The calculation is found to match the experiment in a number of important aspects and the computed flow field is then treated as a realization of a physical turbulent flow. From this flow, a calculation is conducted of the large eddy field and the various averages of the subgrid-scale turbulence that occur in the large eddy simulation equations. These quantities are compared with the predictions of the models that are usually applied in large eddy simulation. The results show that the terms which involve the large-scale field are accurately modeled but the subgrid-scale Reynolds stresses are only moderately well modeled. It is also possible to use the method to predict the constants of the models without reference to experiment. Attempts to find improved models have not met with success.

Clark, R. A.↗

Micro-mechanical analysis of damage growth and fracture in discontinuous fiber reinforced metal matrix composites

An experimental verification is presented for a new two parameter fracture model based on the equivalent remote biaxial stresses (ERBS). A detailed comparison is made between the new theory and the constant K(sub IC) approach of linear elastic fracture mechanics (LEFM). Fracture is predicted through a failure curve representing the change in a variable fracture toughness K(sub IC) with the ERBS ratio B(sub E). The nonsingular term (T) in the series expansion of the near crack-tip transverse stress is included in the model. Experimental results for polymethyl methacrylate (PPMA) show that the theory can account for the effects of geometry on fracture toughness as well as indicate the initiation of crack branching. It is shown that the new criterion predicts failure for PMMA with a 95 percent confidence zone which is nearly three times smaller than that of the LEFM K(sub IC) approach.

Goree, James G.↗

Low-angle normal faults-low differential stress at mid crustal levels

A simple model for frictional slip on pre-existing faults that considers the local stress state near the fault and the effect of non-hydrostatic fluid pressures predicts that low-angle normal faulting is restricted to areas of the crust characterized by low differential stress and nearly lithostatic fluid pressures. The model considers frictional slip on a cohesionless low-angle normal fault governed by the failure criterion tau = mu sub f (sigma (*) sub n) =mu sub f (sigma sub n - P sub f) where tau and sigma sub n are the shear and normal stresses across the fault plane, mu sub f is the static coefficient of friction, and P sub f is the pore fluid pressure. As a first approximation, the model considers a vertical greatest principal compressive stress, sigma sub 1. It is apparent that if slip on low-angle normal faults is governed by the avove frictional failure criterion, slip on the low-angle normal fault occurs only if the least effective principal stress, sigma (*) sub 3 = sigma sub 3 - P sub f, is tensile, whenever tan superscrip -1(mu sub f d, where d is the dip of the fault. If detachment faulting occurs at any significant depth in the crust, P sub f sigma sub 3 is required. In light of this conclusion I allow P sub f to vary as necessary to allow slip on the low-angle normal fault.

Power, W. L.↗

Structure-Property Relationships of Bismaleimides

The purpose of this research was to control and systematically vary the network topology of bismaleimides through cure temperature and chemistry (addition of various coreactants) and subsequently attempt to determine structure-mechanical property relationships. Characterization of the bismaleimide structures by dielectric, rheological, and thermal analyses, and density measurements was subsequently correlated with mechanical properties such as modulus, yield strength, fracture energy, and stress relaxation. The model material used in this investigation was 4,4'-BismaleiMidodIphenyl methane (BMI). BMI was coreacted with either 4,4'-Methylene Dianiline (MDA), o,o'-diallyl bisphenol A (DABA) from Ciba Geigy, or Diamino Diphenyl Sulfone (DDS). Three cure paths were employed: a low- temperature cure of 140 C where chain extension should predominate, a high-temperature cure of 220 C where both chain extension and crosslinking should occur simultaneously, and a low-temperature (140 C) cure followed immediately by a high-temperature (220 C) cure where the chain extension reaction or amine addition precedes BMI homopolymerization or crosslinking. Samples of cured and postcured PMR-15 were also tested to determine the effects of postcuring on the mechanical properties. The low-temperature cure condition of BMI/MDA exhibited the highest modulus values for a given mole fraction of BMI with the modulus decreasing with decreasing concentration of BMI. The higher elastic modulus is the result of steric hindrance by unreacted BMI molecules in the glassy state. The moduli values for the high- and low/high-temperature cure conditions of BMI/MDA decreased as the amount of diamine increased. All the moduli values mimic the yield strength and density trends. For the high-temperature cure condition, the room- temperature modulus remained constant with decreasing mole fraction of BMT for the BMI/DABA and BMI/DDS systems. Postcuring PMR-15 increases the modulus over that of the cured material even though density values of cured and postcured PMR were essentially the same. Preliminary results of a continuous and intermittent stress relaxation experiment for BMI:MDA in a 2:1 molar ratio indicate that crosslinking is occurring when the sample is in the undeformed state. Computer simulation of properties such as density, glass transition temperature, and modulus for the low- temperature cure conditions of BMI/MDA and BMI/DABA were completed. The computer modeling was used to help further understand and confirm the structure characterization results. The simulations correctly predicted the trends of these properties versus mole fraction BMI and were extended to other BMI/diamine systems.

Tenteris-Noebe, Anita D.↗

Free-edge effects in laminates under extension, bending and twisting. II - Sublaminate/layer modeling and analysis

The stress-function-based variational method of Yin (1991) is extended and modified into a combined layer/sublaminate approach applicable to a laminated strip composed of a large number of differently orientated, anisotropic elastic plies. Lekhnitskii's (1963) stress functions are introduced into two interior layers adjacent to a particular interface. The remaining layers are grouped into an upper sublaminate and a lower sublaminate. The stress functions are expanded in truncated power series of the thickness coordinate, and the differential equations governing the coefficient functions are derived by using the complementary virtual work principle. The layer/sublaminate approach limits the dimension of the eigenvalue problem to a fixed number irrespective of the number of layers in the sublaminate, so that reasonably accurate solutions of the interlaminar stresses can be computed with extreme ease. For symmetric, four-layer, angle-ply and cross-ply laminates, a comparison of the previous analysis results based on the pure layer model and new results based on two different layer/sublaminate models indicates reasonable over-all agreement in the interlaminar stresses and superior agreement in the total peeling and shearing force.

Yin, Wan-Lee↗

Thermal expansion mismatch and plasticity in thermal barrier coating

The basic objective of this investigation is the quantitative determination of stress states in a model thermal barrier coating (TBC) as it cools in the air to 600 C from an assumed stress-free state at 700 C. This model is intended to represent a thin plasma-sprayed zirconia-yttria ceramic layer with a nickel chromium-aluminum-yttrium bond coat on a cylindrical substrate made of nickel-based superalloys typically found in gas turbines.

Chang, George C.↗

Two Surface Temperature Retrieval Methods Compared Over Agricultural Lands

Accurate, spatially distributed surface temperatures are required for modeling evapotranspiration (ET) over agricultural fields under wide ranging conditions, including stressed and unstressed vegetation. Modeling approaches that use surface temperature observations, however, have the burden of estimating surface emissivities. Emissivity estimation, the subject of much recent research, is facilitated by observations in multiple thermal infrared bands. But it is nevertheless a difficult task. Using observations from a multiband thermal sensor, the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), estimated surface emissivities and temperatures are retrieved in two different ways: the temperature emissivity separation approach (TES) and the normalized emissivity approach (NEM). Both rely upon empirical relationships, but the assumed relationships are different. TES relies upon a relationship between the minimum spectral emissivity and the range of observed emissivities. NEM relies upon an assumption that at least one thermal band has a pre-determined emissivity (close to 1.0). The benefits and consequences of each approach will be demonstrated for two different landscapes: one in central Oklahoma, USA and another in southern New Mexico.

French, Andrew N.↗

A finite element overlay technique for modeling pinned composite joints

A finite element technique using overlaid plane stress elements is presented for modeling a pinned composite joint. The technique allows for discrete modeling of the pin and web regions of the joint in a two-dimensional (2-D) finite element analysis. Thus, predictive capability is substantially increased without the added complexity of a three-dimensional (3-D) analysis. The overlay technique requires a thorough evaluation of the model constraints between the joint components. Bearing panels are used to account for a reduced compressive modulus under bearing loads. Application of the plane stress overlay technique to the Filament Wound Case (FWC) Program is discussed. The experimental approach to determine bearing panel moduli using a double lap shear test is detailed, and finite element results are correlated with full-scale FWC joint test results.

Colvin, G. E.↗

Generalized Reliability Methodology Applied to Brittle Anisotropic Single Crystals

A generalized reliability model was developed for use in the design of structural components made from brittle, homogeneous anisotropic materials such as single crystals. The model is based on the Weibull distribution and incorporates a variable strength distribution and any equivalent stress failure criteria. In addition to the reliability model, an energy based failure criterion for elastically anisotropic materials was formulated. The model is different from typical Weibull-based models in that it accounts for strength anisotropy arising from fracture toughness anisotropy and thereby allows for strength and reliability predictions of brittle, anisotropic single crystals subjected to multiaxial stresses. The model is also applicable to elastically isotropic materials exhibiting strength anisotropy due to an anisotropic distribution of flaws. In order to develop and experimentally verify the model, the uniaxial and biaxial strengths of a single crystal nickel aluminide were measured. The uniaxial strengths of the <100> and <110> crystal directions were measured in three and four-point flexure. The biaxial strength was measured by subjecting <100> plates to a uniform pressure in a test apparatus that was developed and experimentally verified. The biaxial strengths of the single crystal plates were estimated by extending and verifying the displacement solution for a circular, anisotropic plate to the case of a variable radius and thickness. The best correlation between the experimental strength data and the model predictions occurred when an anisotropic stress analysis was combined with the normal stress criterion and the strength parameters associated with the <110> crystal direction.

Salem, Jonathan A.↗

Thermal shock resistance of lightweight cements developed for geothermal conditions

Cements are a critical component in well construction, as they act to prevent well fluid and gas escape, prevent corrosion of the casing, and strengthen the wellbore to prevent deformation. Under the high temperature/pressure conditions common in geothermal systems, the injection of cold water for energy production is expected to induce cyclic damage to the borehole cement through the rapid temperature fluctuations. These “thermal shocks” are expected to cause casing shrinkage, annulus formation, and cement tensile stresses. To understand the effect of cold water injection on the wellbore environment, a set of rock-cement-steel samples were created to simulate the structure of a geothermal well. Lightweight thermally-insulating cement blends were tested under thermal shock conditions in this study. In each test, the samples were pressurized to an effective pressure of ~3.5 MPa and placed at high temperatures. Thermal shocks were performed by injecting cold water (~10-15 °C) through the samples at a constant rate while keeping the samples at high temperatures until the sample temperature stopped decreasing and deformation ceased. Eight thermal shock tests were conducted with each sample – two at 100 °C and six at 200 °C. Post-tests analysis was then conducted by cutting open each sample to examine the damage in each component of the simulated wellbore. Experimental results suggest that all samples experienced similar degrees of axial and lateral contraction during cold water injection, but for the most part this contraction is recoverable when injection halts. Post-test analysis revealed that fly ash cenosphere pre-treatment had the best effect on improving thermal shock resistance in the cement blends. Thermomechanical modeling of likely stress paths experienced by the cements during heating/cooling cycles shows that elasto-plastic cement constitutive behavior results in most plastic strain occurring during the initial heating steps, with mostly elastic strain occurring during the thermal shock cycles. In conclusion, this agrees with experimental evidence, suggesting that cement damage from shocking occurs via other mechanisms such as chemical alteration, corrosion, and fatigue.

36 MATERIALS SCIENCE↗

Analysis of Bonded Joints Between the Facesheet and Flange of Corrugated Composite Panels

This paper outlines a method for the stress analysis of bonded composite corrugated panel facesheet to flange joints. The method relies on the existing HyperSizer Joints software, which analyzes the bonded joint, along with a beam analogy model that provides the necessary boundary loading conditions to the joint analysis. The method is capable of predicting the full multiaxial stress and strain fields within the flange to facesheet joint and thus can determine ply-level margins and evaluate delamination. Results comparing the method to NASTRAN finite element model stress fields are provided illustrating the accuracy of the method.

Yarrington, Phillip W.↗

Proliferated Resilient Economical Half-Meter Aperture Space Telescopes (PREEMPT)

The PREEMPT LDRD was motivated by a major national need for better space-based imaging systems that are both high performing and affordable for the US Government. Current optical payloads for intelligence, surveillance and reconnaissance (ISR) and space domain awareness (SDA) missions can cost hundreds of millions of dollars and take years to develop, which makes it difficult to build the large constellations required for persistent, world-wide coverage. To address this, the team aimed to advance a different kind of large-aperture (>25 cm) telescope, called a monolithic telescope, in which key optical surfaces are built into a single piece of fused silica. This design greatly reduces payload and spacecraft complexity, the need for precision focus actuators, improves mechanical and thermal robustness, and lowers cost when compared with traditional Cassegrain telescopes that rely on many precisely aligned components. The project focused on five primary thrusts. The first thrust was to advance the concept of a V10, 25 cm, monolithic telescope forward from optical design to flight-ready stage. This was accomplished in partnership with Optimax, who delivered the first test unit in the early stages of the LDRD. The team developed several technologies necessary for this optic to be integrated into a flight demonstration. These include carbon fiber housings, highly detailed structural and thermal models and stress-reducing elastic averaging Hirth groove designs. These technologies resulted in a successful maturation of the optic, which is now slated to fly in late 2026/early 2027 for a demonstration mission. The second thrust was to advance the manufacturability of these optics. In collaboration with NIF’s optical manufacturing shop, we reduced polishing time from 480 hours to 65 hours through the implementation of optimized processes and new tools. The NIF team utilized a conceptual V8 (18 cm) optic to demonstrate this optimization, though it can be applied to the rest of the monolithic optic portfolio. Third, the team developed the first conceptual V20 (50 cm) payload, which is slated to be the next generation of LLNL optical payload systems. A set of structural, dynamic and thermal simulations were performed to identify potential challenges in the future development of this payload. Early-stage simulations suggest the payload is feasible, though thermal management will be the key focus area to maintain optimal performance. Fourth, a non-linear model of Viton was developed, to further enhance the reliability of our structural and dynamic models for future payloads. Viton acts as the primary interface material between the optic and its housing. Lastly, the team focused on successfully displaying the feasibility of using additively manufactured metal composites for optical space payloads. The team successfully demonstrated layer by layer deposition of Al-SiC composites, which have highly tunable structural and coefficient of thermal expansion (CTE) properties. These are crucial for optical payloads because CTE mismatch is one of the causes for degraded optical performance for telescopes in orbit. Overall, the work showed that monolithic telescopes could become a practical, lower-cost path to high-resolution space imaging for both national security and scientific missions.

42 ENGINEERING↗

Postseismic viscoelastic surface deformation and stress. Part 1: Theoretical considerations, displacement and strain calculations

A model of viscoelastic deformations associated with earthquakes is presented. A strike-slip fault is represented by a rectangular dislocation in a viscoelastic layer (lithosphere) lying over a viscoelastic half-space (asthenosphere). Deformations occur on three time scales. The initial response is governed by the instantaneous elastic properties of the earth. A slower response is associated with viscoelastic relaxation of the lithosphere and a yet slower response is due to viscoelastic relaxation of the asthenosphere. The major conceptual contribution is the inclusion of lithospheric viscoelastic properties into a dislocation model of earthquake related deformations and stresses. Numerical calculations using typical fault parameters reveal that the postseismic displacements and strains are small compared to the coseismic ones near the fault, but become significant further away. Moreover, the directional sense of the deformations attributable to the elastic response, the lithospheric viscoelastic softening, and the asthenospheric viscoelastic flow may differ and depend on location and model details. The results and theoretical arguments suggest that the stress changes accompanying lithospheric relaxation may also be in a different sense than and be larger than the strain changes.

Cohen, S. C.↗

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.↗