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,207 records · Page 67

A direct-adjoint approach for material point model calibration with application to plasticity

Here, this paper proposes a new approach for the calibration of material parameters in local elastoplastic constitutive models. The calibration is posed as a constrained optimization problem, where the constitutive model evolution equations for a single material point serve as constraints. The objective function quantifies the mismatch between the stress predicted by the model and corresponding experimental measurements. To improve calibration efficiency, a novel direct-adjoint approach is presented to compute the Hessian of the objective function, which enables the use of second-order optimization algorithms. Automatic differentiation is used for gradient and Hessian computations. Two numerical examples are employed to validate the Hessian matrices and to demonstrate that the Newton–Raphson algorithm consistently outperforms gradient-based algorithms such as L-BFGS-B.

36 MATERIALS SCIENCE↗

An experimental and numerical investigation of shock-wave induced turbulent boundary-layer separation at hypersonic speeds

An experiment designed to test and guide computations of the interaction of an impinging shock wave with a turbulent boundary layer is described. Detailed mean flow-field and surface data are presented for two shock strengths which resulted in attached and separated flows, respectively. Numerical computations, employing the complete time-averaged Navier-Stokes equations along with algebraic eddy-viscosity and turbulent Prandtl number models to describe shear stress and heat flux, are used to illustrate the dependence of the computations on the particulars of the turbulence models. Models appropriate for zero-pressure-gradient flows predicted the overall features of the flow fields, but were deficient in predicting many of the details of the interaction regions. Improvements to the turbulence model parameters were sought through a combination of detailed data analysis and computer simulations which tested the sensitivity of the solutions to model parameter changes. Computer simulations using these improvements are presented and discussed.

Marvin, J. G.↗

Analysis of the shearout failure mode in composite bolted joints

A semi-empirical shearout strength model has been formulated for the analysis of composite bolted joints with allowance for the effects of joint geometry. The model employs a polynomial stress function in conjunction with a point stress failure criterion to predict strength as a function of fastener size, edge distance, and half spacing. The stress function is obtained by two-dimensional plane-stress finite element analysis using quadrilateral elements with orthotropic material properties. Comparison of experimentally determined shearout strength data with model predicted failures has substantiated the accuracy of the model.

Wilson, D. W.↗

Sensitivity analysis of the non-linear dynamic viscoplastic response of 2-d structures with respect to material parameters

A computational procedure is presented for evaluating the sensitivity coefficients of the viscoplastic response of structures subjected to dynamic loading. A state of plane stress is assumed to exist in the structure, a velocity strain-Cauchy stress formulation is used, and the geometric non-linearities arising from large strains are incorporated. The Jaumann rate is used as a frame indifferent stress rate. The material model is chosen to be isothermal viscoplasticity, and an associated flow rule is used with a von Mises effective stress. The equations of motion emanating from a finite element semi-discretization are integrated using an explicit central difference scheme with an implicit stress update. The sensitivity coefficients are evaluated using a direct differentiation approach. Since the domain of integration is the current configuration, the sensitivity coefficients of the spatial derivatives of the shape functions must be included. Numerical results are presented for a thin plate with a central cutout subjected to an in-plane compressive loading. The sensitivity coefficients are generated by evaluating the derivatives of the response quantities with respect to Young's modulus, and two of the material parameters characterizing the viscoplastic response. Time histories of the response and sensitivity coefficients, and spatial distributions at selected times are presented.

Kulkarni, Makarand↗

The interpretation of crustal dynamics data in terms of plate motions and regional deformation near plate boundaries

During our participation in the NASA Crustal Dynamics Project under NASA contract NAS-27339 and grant NAG5-814 for the period 1982-1991, we published or submitted for publication 30 research papers and 52 abstracts of presentations at scientific meetings. In addition, five M.I.T. Ph.D. students (Eric Bergman, Steven Bratt, Dan Davis, Jeanne Sauber, Anne Sheehan) were supported wholly or in part by this project during their thesis research. Highlights of our research progress during this period include the following: application of geodetic data to determine rates of strain in the Mojave block and in central California and to clarify the relation of such strain to the San Andreas fault and Pacific-North American plate motions; application of geodetic data to infer post seismic deformation associated with large earthquakes in the Imperial Valley, Hebgen Lake, Argentina, and Chile; determination of the state of stress in oceanic lithosphere from a systematic study of the centroid depths and source mechanisms of oceanic intraplate earthquakes; development of models for the state of stress in young oceanic regions arising from the differential cooling of the lithosphere; determination of the depth extent and rupture characteristics of oceanic transform earthquakes; improved determination of earthquake slip vectors in the Gulf of California, an important data set for the estimation of Pacific-North American plate motions; development of models for the state of stress and mechanics of fold-and-thrust belts and accretionary wedges; development of procedures to invert geoid height, residual bathymetry, and differential body wave travel time residuals for lateral variations in the characteristic temperature and bulk composition of the oceanic upper mantle; and initial GPS measurements of crustal deformation associated with the Imperial-Cerro Prieto fault system in southern California and northern Mexico. Full descriptions of the research conducted on these topics may be found in the Semi-Annual status Reports submitted regularly to NASA over the course of this project and in the publications listed.

Solomon, Sean C.↗

Transient water stress in a vegetation canopy - Simulations and measurements

Consideration is given to observational and modeling evidence of transient water stress, the effects of the transpiration plateau on the canopy radiometric temperature, and the factors responsible for the onset of the transpiration plateau, such as soil moisture. Attention is also given to the point at which the transient stress can be detected by remote measurement of surface temperature.

Carlson, Toby N.↗

Environmental Barrier Coating (EBC) Durability Modeling; An Overview and Preliminary Analysis

A study outlining a fracture mechanics based model that is being developed to investigate crack growth and spallation of environmental barrier coating (EBC) under thermal cycling conditions is presented. A description of the current plan and a model to estimate thermal residual stresses in the coating and preliminary fracture mechanics concepts for studying crack growth in the coating are also discussed. A road map for modeling life and durability of the EBC and the results of FEA model(s) developed for predicting thermal residual stresses and the cracking behavior of the coating are generated and described. Further initial assessment and preliminary results showed that developing a comprehensive EBC life prediction model incorporating EBC cracking, degradation and spalling mechanism under stress and temperature gradients typically seen in turbine components is difficult. This is basically due to mismatch in thermal expansion difference between sub-layers of EBC as well as between EBC and substrate, diffusion of moisture and oxygen though the coating, and densification of the coating during operating conditions as well as due to foreign object damage, the EBC can also crack and spall from the substrate causing oxidation and recession and reducing the design life of the EBC coated substrate.

Abdul-Aziz, A.↗

A Model of the THUNDER Actuator

A THUNDER actuator is a composite of three thin layers, a metal base, a piezoelectric wafer and a metal top cover, bonded together under pressure and at high temperature with the LaRC SI polyimid adhesive. When a voltage is applied between the metal layers across the PZT the actuator will bend and can generate a force. This document develops and describes an analytical model the transduction properties of THUNDER actuators. The model development is divided into three sections. First, a static model is described that relates internal stresses and strains and external displacements to the thermal pre-stress and applied voltage. Second, a dynamic energy based model is described that allows calculation of the resonance frequencies, developed force and electrical input impedance. Finally, a fully coupled electro-mechanical transducer model is described. The model development proceeds by assuming that both the thermal pre-stress and the piezoelectric actuation cause the actuator to deform in a pure bend in a single plane. It is useful to think of this as a two step process, the actuator is held flat, differential stresses induce a bending moment, the actuator is released and it bends. The thermal pre-stress is caused by the different amounts that the constituent layers shrink due to their different coefficients of thermal expansion. The adhesive between layers sets at a high temperature and as the actuator cools, the metal layers shrink more than the PZT. The PZT layer is put into compression while the metal layers are in tension. The piezoelectric actuation has a similar effect. An applied voltage causes the PZT layer to strain, which in turn strains the two metal layers. If the PZT layer expands it will put the metal layers into tension and PZT layer into compression. In both cases, if shear force effects are neglected, the actuator assembly will experience a uniform in-plane strain. As the materials each have a different elastic modulus, different stresses will develop in each layer and these stresses will induce a bending moment. When the actuator is released from its flat configuration, the differential stresses are relieved as the actuator bends.

Curtis, Alan R. D.↗

Langley facility for tests at Mach 7 of subscale, hydrogen-burning, airframe-integratable, scramjet models

Modifications to a 20-megawatt arc-heated facility for testing a hydrogen-burning, airframe-integratable, subscale, scramjet model are described. Arc-heated flow is mixed with unheated air to furnish a test flow duplicating Mach 7 flight. (Stagnation temperature is 2220 K.) Modifications to the commercially available heater to improve survivability and smoothness are described. Pitot profiles show uniform flow and a slightly thinner nozzle boundary layer than predicted. Comparison of the tunnel boundary layer, which will be ingested by the engine model, with the boundary layer that a flight engine might ingest from its vehicle forebody shows a difference in the density distribution through the boundary layer. Calculations of wall heating and transient wall temperatures of the engine model show that for a 30-sec burn, the heat sink model requires cooling at selected locations to avoid thermal-stress, cycle-life problems. Model performance predictions show that fuel equivalence ratio and nozzle exit area both have large effects on thrust. Average inlet entrance Mach number (as affected by boundary-layer ingestion) has little effect on thrust.

Boatright, W. B.↗

Variational mechanics analysis of the stresses in microdrop debond specimens

A recently derived variational mechanics analysis of stresses in single-fiber model composites has been applied to the analysis of the stresses in the microdrop debond specimen. The new analysis is more accurate than the commonly applied shear-lag or elastic-plastic analyses. The results from a sample stress state calculation suggest that interfacial failure between the fiber and the microdrop is by mode I or opening mode failure at the beginning of the microdrop. The opening mode failure is caused by a large tensile radial stress at the fiber/matrix interface. Previous analyses of microdrop debond data have been in terms of a shear strength. It is suggested that these analyses misrepresent microdrop debond results and recommend instead a failure analysis based on energy release rate and interfacial fracture toughness. A procedure for calculating the energy release rate for the growth of an interfacial crack is described.

Scheer, Robert J.↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling↗

Mission Cycle Life Prediction for C/SiC Composites

The residual strength and life of C/SiC is dominated by carbon interface and fiber oxidation if thermal and mechanical stresses combine with thermal expansion mismatch in a manner such that seal coat and matrix cracks are open to allow oxygen ingress. Under these circumstances life can be predicted by simple oxidation based models based on reaction controlled kinetics at low temperature and gas phase diffusion controlled kinetics at high temperature. Key life governing variables in these models include temperature, applied stress, initial strength, oxygen partial pressure and total pressure. In this paper these models are described and extended to mission cycle life prediction.

Levine, Stanley R.↗

Wind-Stress Simulations and Equatorial Dynamics in an AGCM: Basic Results from a 1979-1999 Forced SST Experiment - Part 1

This is the first of a two part study examining the connection of the equatorial momentum budget in an AGCM (Atmospheric General Circulation Model), with simulated equatorial surface wind stresses over the Pacific. The AGCM used in this study forms part of a newly developed coupled forecasting system used at NASA's Seasonal- to-Interannual Prediction Project. Here we describe the model and present results from a 20-year (1979-1999) AMIP-type experiment forced with observed SSTs (Sea Surface Temperatures). Model results are compared them with available observational data sets. The climatological pattern of extra-tropical planetary waves as well as their ENSO-related variability is found to agree quite well with re-analysis estimates. The model's surface wind stress is examined in detail, and reveals a reasonable overall simulation of seasonal interannual variability, as well as seasonal mean distributions. However, an excessive annual oscillation in wind stress over the equatorial central Pacific is found. We examine the model's divergent circulation over the tropical Pacific and compare it with estimates based on re-analysis data. These comparisons are generally good, but reveal excessive upper-level convergence in the central Pacific. In Part II of this study a direct examination of individual terms in the AGCM's momentum budget is presented. We relate the results of this analysis to the model's simulation of surface wind stress.

Bacmeister, Julio T.↗

Thermal stress effects in intermetallic matrix composites

Intermetallic matrix composites develop residual stresses from the large thermal expansion mismatch (delta-alpha) between the fibers and matrix. This work was undertaken to: establish improved techniques to measure these thermal stresses in IMC's; determine residual stresses in a variety of IMC systems by experiments and modeling; and, determine the effect of residual stresses on selected mechanical properties of an IMC. X ray diffraction (XRD), neutron diffraction (ND), synchrotron XRD (SXRD), and ultrasonics (US) techniques for measuring thermal stresses in IMC were examined and ND was selected as the most promising technique. ND was demonstrated on a variety of IMC systems encompassing Ti- and Ni-base matrices, SiC, W, and Al2O3 fibers, and different fiber fractions (Vf). Experimental results on these systems agreed with predictions of a concentric cylinder model. In SiC/Ti-base systems, little yielding was found and stresses were controlled primarily by delta-alpha and Vf. In Ni-base matrix systems, yield strength of the matrix and Vf controlled stress levels. The longitudinal residual stresses in SCS-6/Ti-24Al-llNb composite were modified by thermomechanical processing. Increasing residual stress decreased ultimate tensile strength in agreement with model predictions. Fiber pushout strength showed an unexpected inverse correlation with residual stress. In-plane shear yield strength showed no dependence on residual stress. Higher levels of residual tension led to higher fatigue crack growth rates, as suggested by matrix mean stress effects.

Wright, P. K.↗

Predicting transcriptional responses to cold stress across plant species

Although genome-sequence assemblies are available for a growing number of plant species, gene-expression responses to stimuli have been cataloged for only a subset of these species. Many genes show altered transcription patterns in response to abiotic stresses. However, orthologous genes in related species often exhibit different responses to a given stress. Accordingly, data on the regulation of gene expression in one species are not reliable predictors of orthologous gene responses in a related species. Here, we trained a supervised classification model to identify genes that transcriptionally respond to cold stress. A model trained with only features calculated directly from genome assemblies exhibited only modest decreases in performance relative to models trained by using genomic, chromatin, and evolution/diversity features. Models trained with data from one species successfully predicted which genes would respond to cold stress in other related species. Cross-species predictions remained accurate when training was performed in cold-sensitive species and predictions were performed in cold-tolerant species and vice versa. Models trained with data on gene expression in multiple species provided at least equivalent performance to models trained and tested in a single species and outperformed single-species models in cross-species prediction. These results suggest that classifiers trained on stress data from well-studied species may suffice for predicting gene-expression patterns in related, less-studied species with sequenced genomes.

54 ENVIRONMENTAL SCIENCES↗

Second-Order Modeling of Low-Reynolds-Number Turbulence Near Walls

This paper presents a set of second-order closure models for low-Reynolds-number turbulence near the wall. Existing closure models for the Reynolds-stress equations were modified to show proper near-wall behavior. A dissipation-rate equation for the turbulent kinetic energy is also reformulated. The proposed models satisfy realizability and will not produce unphysical behavior. Fully developed channel flows are used for model testing. The equations are solved for the mean velocity, the Reynolds stresses, and the dissipation rate of the turbulent kinetic energy. The calculations are compared with both direct numerical simulations and with measurements. It is shown that the present models perform well in predicting the behavior of the turbulence near a wall. Significant improvements over previous models in predicting the components of the Reynolds stress tensor are obtained in the present models.

Shih, T.-H.↗

Considerations on the moving contact-line singularity, with application to frictional drag on a slender drop

It has previously been shown that the no-slip boundary conditions leads to a singularity at a moving contact line and that this presumes some form of slip. Present considerations on the energetics of slip due to shear stress lead to a yield stress boundary condition. A model for the distortion of the liquid state near solid boundaries gives a physical basis for this boundary condition. The yield stress condition is illustrated by an analysis of a slender drop rolling down an incline. That analysis provides a formula for the frictional drag resisting the drop movement. With the present boundary condition, the length of the slip region becomes a property of the fluid flow.

Durbin, P. A.↗