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,009 records · Page 56

Analysis of multifastener composite joints

A new numerical procedure for determining load proportioning in multiple-hole mechanical joints in composite plates is presented. The joints are loaded in double lap fashion in tension with pins through the holes. The commercial finite-element program ABAQUS is used to predict the load proportioning among holes using two independent plane stress finite-element models, one representing the composite inner lap and one representing the two steel outer laps, interacting through rigid circular surfaces. The circular surfaces effectively represent rigid pins. Load proportioning is predicted for a number of geometries. Excellent correlation with experimental data is obtained. Experimental and computed surface strains are also found to compare well. The assumption of a radial cosine distribution of contact stress between pin and hole is shown to be substantially in error for some holes.

Griffin, O. H., Jr.↗

Factors Influencing Dwell Fatigue Cracking in Notches of Powder Metallurgy Superalloys

The influences of heat treatment and cyclic dwells on the notch fatigue resistance of powder metallurgy disk superalloys were investigated for low solvus high refractory (LSHR) and ME3 disk alloys. Disks were processed to produce material conditions with varied microstructures and associated mechanical properties. Notched specimens were first subjected to baseline dwell fatigue cycles having a dwell at maximum load, as well as tensile, stress relaxation, creep rupture, and dwell fatigue crack growth tests at 704 C. Several material heat treatments displayed a bimodal distribution of fatigue life with the lives varying by two orders-of-magnitude, while others had more consistent fatigue lives. This response was compared to other mechanical properties, in search of correlations. The wide scatter in baseline dwell fatigue life was observed only for material conditions resistant to stress relaxation. For selected materials and conditions, additional tests were then performed with the dwells shifted in part or in total to minimum tensile load. The tests performed with dwells at minimum load exhibited lower fatigue lives than max dwell tests, and also exhibited early crack initiation and a substantial increase in the number of initiation sites. These results could be explained in part by modeling evolution of peak stresses in the notch with continued dwell fatigue cycling. Fatigue-environment interactions were determined to limit life for the fatigue cycles with dwells.

Gabb, T. P.↗

Scanning electron acoustic microscopy of residual stresses in ceramics - Theory and experiment

The paper presents a three-dimensional mathematical model of signal generation and contrast in brittle materials and uses the model to simulate the effect of residual stress fields on the scanning electron acoustic microscopy (SEAM)-generated electron acoustic signal. According to the model, a positive (tensile) strain produces an increase in the output signal, whereas a negative (compressive) strain produces a decrease in the ouput signal. Dark field contrast conditions occur at a chopping frequency at which V2 - V1 is greater than 0 (where V2 = V is the SEAM output in a region of residual stresses, and V1 is the output in a stress-free region of the sample). Under ideal conditions (maximum contrast) V1 approaches zero. It was found that tensile strains of the order 0.2-0.3 percent, possible in brittle materials, would produce a variation of the acoustic output signal of the order 10 nV (about 1 percent), well within the image contrast and signal processing capability of the SEAM electronics.

Cantrell, John H.↗

Computational Analyses of Propulsion Aeroacoustics for Mixed Flow Nozzle Pylon Installation at Takeoff

A CFD analyses is presented for a set of baseline and noise suppression mixed flow nozzles with and without a pylon installation. The five model configurations are as follows; a baselinecore/fan dual-stream nozzle with an external plug, a chevron mixer nozzle with a peak on the symmetry plane with external plug, both of the above nozzles with an installed bifurcatingpylon and lastly a clocked chevron mixer nozzle such that a trough is aligned with the center of the pylon. The fluid flow is simulated by solving the asymptotically steady, compressible, Reynolds-averaged Navier-Stokes equations using an implicit, up-wind, flux-difference splitting finite volume scheme and standard two equation k-epsilon turbulence model with a linear stress representation. All computations are performed using the multiblock, parallel, structuredcode PAB3D. Results indicate that the clocked chevron with pylon case achieves the most optimal levels of average and peak turbulence kinetic energy and vorticity and therefore is expected to be the quietest of the five configurations tested. Further study is required to refine expressions which are indicative of noise and mate these with rigorous noise prediction models.

Massey, Steven J.↗

Physical Activity and Bone Density in Women

A mathematical model of bone density regulation as a function of the daily tissue "effective" stress has been derived. Using the model, the influence of daily activity in the form of a daily loading history has been related to bone density of the calcaneus. The theory incorporates a stress exponent m to account for differences in the importance of magnitude and number of load cycles experienced during daily activity. We have derived a parameter from the model, the "Bone Density Index" (BDI). We have developed a method of collecting daily habitual loading histories using an insole force sensor interfaced to a portable digital data logger carried in a fanny pack. Our goal for this study was to determine a stress exponent, m, relating GRFz history to Calcaneal Bone Mineral Density (CBMD).

Bowley, Susan M.↗

Modeling the Nonlinear, Strain Rate Dependent Deformation of Shuttle Leading Edge Materials with Hydrostatic Stress Effects Included

An analysis method based on a deformation (as opposed to damage) approach has been developed to model the strain rate dependent, nonlinear deformation of woven ceramic matrix composites, such as the Reinforced Carbon Carbon (RCC) material used on the leading edges of the Space Shuttle. In the developed model, the differences in the tension and compression deformation behaviors have also been accounted for. State variable viscoplastic equations originally developed for metals have been modified to analyze the ceramic matrix composites. To account for the tension/compression asymmetry in the material, the effective stress and effective inelastic strain definitions have been modified. The equations have also been modified to account for the fact that in an orthotropic composite the in-plane shear response is independent of the stiffness in the normal directions. The developed equations have been implemented into LS-DYNA through the use of user defined subroutines (UMATs). Several sample qualitative calculations have been conducted, which demonstrate the ability of the model to qualitatively capture the features of the deformation response present in woven ceramic matrix composites.

Goldberg, Robert K.↗

Navier-Stokes Computations With One-Equation Turbulence Model for Flows Along Concave Wall Surfaces

This report presents the use of a time-marching three-dimensional compressible Navier-Stokes equation numerical solver with a one-equation turbulence model to simulate the flow fields developed along concave wall surfaces without and with a downstream extension flat wall surface. The 3-D Navier- Stokes numerical solver came from the NASA Glenn-HT code. The one-equation turbulence model was derived from the Spalart and Allmaras model. The computational approach was first calibrated with the computations of the velocity and Reynolds shear stress profiles of a steady flat plate boundary layer flow. The computational approach was then used to simulate developing boundary layer flows along concave wall surfaces without and with a downstream extension wall. The author investigated the computational results of surface friction factors, near surface velocity components, near wall temperatures, and a turbulent shear stress component in terms of turbulence modeling, computational mesh configurations, inlet turbulence level, and time iteration step. The computational results were compared with existing measurements of skin friction factors, velocity components, and shear stresses of the developing boundary layer flows. With a fine computational mesh and a one-equation model, the computational approach could predict accurately the skin friction factors, near surface velocity and temperature, and shear stress within the flows. The computed velocity components and shear stresses also showed the vortices effect on the velocity variations over a concave wall. The computed eddy viscosities at the near wall locations were also compared with the results from a two equation turbulence modeling technique. The inlet turbulence length scale was found to have little effect on the eddy viscosities at locations near the concave wall surface. The eddy viscosities, from the one-equation and two-equation modeling, were comparable at most stream-wise stations. The present one-equation turbulence model is an effective approach for turbulence modeling in the near solid wall surface region of flow over a concave wall.

Wang, Chi R.↗

Modeling the Nonlinear, Strain Rate Dependent Deformation of Woven Ceramic Matrix Composites With Hydrostatic Stress Effects Included

An analysis method based on a deformation (as opposed to damage) approach has been developed to model the strain rate dependent, nonlinear deformation of woven ceramic matrix composites with a plain weave fiber architecture. In the developed model, the differences in the tension and compression response have also been considered. State variable based viscoplastic equations originally developed for metals have been modified to analyze the ceramic matrix composites. To account for the tension/compression asymmetry in the material, the effective stress and effective inelastic strain definitions have been modified. The equations have also been modified to account for the fact that in an orthotropic composite the in-plane shear stiffness is independent of the stiffness in the normal directions. The developed equations have been implemented into a commercially available transient dynamic finite element code, LS-DYNA, through the use of user defined subroutines (UMATs). The tensile, compressive, and shear deformation of a representative plain weave woven ceramic matrix composite are computed and compared to experimental results. The computed values correlate well to the experimental data, demonstrating the ability of the model to accurately compute the deformation response of woven ceramic matrix composites.

Goldberg, Robert K.↗

Large area sheet task: Advanced dendritic web growth development

The growth of silicon dendritic web for photovoltaic applications was investigated. The application of a thermal model for calculating buckling stresses as a function of temperature profile in the web is discussed. Lid and shield concepts were evaluated to provide the data base for enhancing growth velocity. An experimental web growth machine which embodies in one unit the mechanical and electronic features developed in previous work was developed. In addition, evaluation of a melt level control system was begun, along with preliminary tests of an elongated crucible design. The economic analysis was also updated to incorporate some minor cost changes. The initial applications of the thermal model to a specific configuration gave results consistent with experimental observation in terms of the initiation of buckling vs. width for a given crystal thickness.

Duncan, C. S.↗

Investigating Some Technical Issues on Cohesive Zone Modeling of Fracture

This study investigates some technical issues related to the use of cohesive zone models (CZMs) in modeling fracture processes. These issues include: why cohesive laws of different shapes can produce similar fracture predictions; under what conditions CZM predictions have a high degree of agreement with linear elastic fracture mechanics (LEFM) analysis results; when the shape of cohesive laws becomes important in the fracture predictions; and why the opening profile along the cohesive zone length needs to be accurately predicted. Two cohesive models were used in this study to address these technical issues. They are the linear softening cohesive model and the Dugdale perfectly plastic cohesive model. Each cohesive model constitutes five cohesive laws of different maximum tractions. All cohesive laws have the same cohesive work rate (CWR) which is defined by the area under the traction-separation curve. The effects of the maximum traction on the cohesive zone length and the critical remote applied stress are investigated for both models. For a CZM to predict a fracture load similar to that obtained by an LEFM analysis, the cohesive zone length needs to be much smaller than the crack length, which reflects the small scale yielding condition requirement for LEFM analysis to be valid. For large-scale cohesive zone cases, the predicted critical remote applied stresses depend on the shape of cohesive models used and can significantly deviate from LEFM results. Furthermore, this study also reveals the importance of accurately predicting the cohesive zone profile in determining the critical remote applied load.

Wang, John T.↗

Improvements in Simulating a Mach 0.80 Transonic Truss-Braced Wing Configuration using the Spalart-Allmaras and k-ω SST Turbulence

Teams from the NASA Ames Research Centers (ARC) and Langley Research Center (LaRC) have been working on validating their computational fluid dynamics (CFD) results for the Boeing Mach 0.80 Transonic Truss-Braced Wing (TTBW) configuration. Experimental data used for the validation were gathered from a test conducted of a 4.5% scale Mach 0.80 TTBW model in the NASA Ames Research Center 11-by 11-Foot Transonic Wind Tunnel. The CFD simulations were initially run with both LAVA and USM3D Mixed Element solvers utilizing the Spalart-Allmaras (SA) turbulence model. A discrepancy was observed between CFD and experimental loads and moments ranging on average from 0.047 to 0.063 for lift coefficient, 16.5 to 27 in drag counts, and -0.02 to -0.032 for pitching moment, varying with the solver used. With introducing the refactored version of LAVA these values dropped to 0.034 to 0.047 for lift coefficient, 7.74 to 16.5 in drag counts, and lastly -0.012 to -0.02 for pitching moment. Based on findings from Boeing, who observed an improved comparison to experimental data when using the k-ω Shear Stress Transport (SST) turbulence model, the NASA teams conducted simulations with SST to investigate turbulence modeling sensitivities. OVERFLOW and USM3D V6 solvers were used for the comparison of SST and SA simulations. CFD results using the SST turbulence model demonstrate an improvement in matching with experimental CL values, reducing the discrepancies seen by 0.021-0.055 (62-93% reduction in discrepancy respective to the OVERFLOW and USM3D V6). The SST model has varying effects on CD, based on the solver / grid paradigm; for USM3D this prediction is only improved at higher angles of attack (above the cruise design point) resulting in a reduced discrepancy of 8.8–25 drag counts However, at the mid-to-low angles of attack, SST increases the CD discrepancy by 4-55 counts, increasing at lower angles of attack. For OVERFLOW this discrepancy only exists at the low angles of attack and maxes out at 10 counts of drag. Above α = 1.5° OVERFLOW SST reduces discrepancy by 9.3 drag counts on average.

AATT↗

Orion Handling Qualities During ISS Rendezvous and Docking

The Orion spacecraft was designed to rendezvous with multiple vehicles in low earth orbit (LEO) and beyond. To perform the required rendezvous and docking task, Orion must provide enough control authority to perform coarse translational maneuvers while maintaining precision to perform the delicate docking corrections. While Orion has autonomous docking capabilities, it is expected that final approach and docking operations with the International Space Station (ISS) will initially be performed in a manual mode. A series of evaluations was conducted by NASA and Lockheed Martin at the Johnson Space Center to determine the handling qualities (HQ) of the Orion spacecraft during different docking and rendezvous conditions using the Cooper-Harper scale. This paper will address the specifics of the handling qualities methodology, vehicle configuration, scenarios flown, data collection tools, and subject ratings and comments. The initial Orion HQ assessment examined Orion docking to the ISS. This scenario demonstrates the Translational Hand Controller (THC) handling qualities of Orion. During this initial assessment, two different scenarios were evaluated. The first was a nominal docking approach to a stable ISS, with Orion initializing with relative position dispersions and a closing rate of approximately 0.1 ft/sec. The second docking scenario was identical to the first, except the attitude motion of the ISS was modeled to simulate a stress case ( 1 degree deadband per axis and 0.01 deg/sec rate deadband per axis). For both scenarios, subjects started each run on final approach at a docking port-to-port range of 20 ft. Subjects used the THC in pulse mode with cues from the docking camera image, window views, and range and range rate data displayed on the Orion display units. As in the actual design, the attitude of the Orion vehicle was held by the automated flight control system at 0.5 degree deadband per axis. Several error sources were modeled including Reaction Control System (RCS) jet angular and position misalignment, RCS thrust magnitude uncertainty, RCS jet force direction uncertainty due to self plume impingement, and Orion center of mass uncertainty.

Hart, Jeremy J.↗

Failure Modes in Capacitors When Tested Under a Time-Varying Stress

Power-on failure has been the prevalent failure mechanism for solid tantalum capacitors in decoupling applications. A surge step stress test (SSST) has been previously applied to identify the critical stress level of a capacitor batch to give some predictability to the power-on failure mechanism [1]. But SSST can also be viewed as an electrically destructive test under a time-varying stress (voltage). It consists of rapidly charging the capacitor with incremental voltage increases, through a low resistance in series, until the capacitor under test is electrically shorted. When the reliability of capacitors is evaluated, a highly accelerated life test (HALT) is usually adopted since it is a time-efficient method of determining the failure mechanism; however, a destructive test under a time-varying stress such as SSST is even more time efficient. It usually takes days or weeks to complete a HALT test, but it only takes minutes for a time-varying stress test to produce failures. The advantage of incorporating a specific time-varying stress profile into a statistical model is significant in providing an alternative life test method for quickly revealing the failure mechanism in capacitors. In this paper, a time-varying stress that mimics a typical SSST has been incorporated into the Weibull model to characterize the failure mechanism in different types of capacitors. The SSST circuit and transient conditions for correctly surge testing capacitors are discussed. Finally, the SSST was applied for testing Ta capacitors, polymer aluminum capacitors (PA capacitors), and multi-layer ceramic (MLC) capacitors with both precious metal electrodes (PME) and base metal electrodes (BME). The test results are found to be directly associated with the dielectric layer breakdown in Ta and PA capacitors and are independent of the capacitor values, the way the capacitors were built, and the capacitors manufacturers. The test results also show that MLC capacitors exhibit surge breakdown voltages much higher than the rated voltage and that the breakdown field is inversely proportional to the dielectric layer thickness. The SSST data can also be used to comparatively evaluate the voltage robustness of capacitors for decoupling applications.

Liu, David (Donhang)↗

The prediction of two- and three-dimensional asymmetrical turbulent wakes - A comparision of the performance of three turbulence models for the effects of streamline curvature and rotation

Two- and three-dimensional turbulent wakes which develop under the influence of streamline curvature and rotation were calculated with three turbulence closure models. The first model comprised of transport equations for the turbulent kinetic energy and the rate of energy dissipation. The second model comprised of equations for the rate of turbulent kinetic energy dissipation and Reynolds stresses, but the effects of the convection and diffusion in the Reynolds stress transport equation were handled collectively. The third model utilizes equations of turbulent kinetic energy dissipation and Reynolds stresses in nearly exact form. All of the three models were modified for the effects of streamline curvature. The second and third models include the effects of rotation through the rotation-originated redistribution term in the transport equation of Reynolds stresses. The numerical results for the wakes of airfoils, cascades and turbomachinery rotor blades demonstrate that the second and third models provide accurate predictions, but computer time and storage can be considerably saved with the second model.

Hah, C.↗

Exploring the Langtry-Menter Transition Model for High Speed Applications Using FUN3D

A series of Reynolds-averaged Navier-Stokes (RANS) simulations were performed using the FUN3D flow solver to explore the capabilities of the Langtry-Menter Shear-Stress Transport (LM-SST) transition model for predicting transition for aircraft inlet applications. Two geometries were simulated: a zeropressure-gradient flat plate and an axisymmetric cone exposed to hypersonic flow. In addition to the transition-sensitized LM-SST model investigations, simulations were run with the one-equation SpalartAllmaras (SA) and the two-equation Menter Shear-Stress Transport (SST-V) RANS models in fully turbulent mode to identify the natural RANS model transition behavior as a function of Mach number when executed in fully turbulent mode. The flat plate simulations showed that (1) the transition model was able to predict rapid transition at a freestream Mach number of 0.2, which is expected but (2) the predicted transition location moved downstream as the freestream Mach number was increased for the simulations that used the SST-V turbulence model. The latter is significant as it is usually assumed that one- and two-equation turbulence models will produce fully turbulent flow very near the boundary layer origin. The flat plate simulation freestream Mach number trend was confirmed with simulations using the Wind-US code, which also saw a similar trend when employing the SA turbulence model. For the axisymmetric cone simulations, the transition location was highly sensitive to the inflow turbulence levels. This is significant as the prediction of the transition location is crucial when trying to predict inlet performance, especially for hypersonic vehicle applications. It was also noted that the predicted transition location for the cone when using the SST-V turbulence model agreed well with the predicted transition location from the equivalent zero-pressure-gradient cold wall flat plate case.

Transition Model↗

Thermal Stresses in Silicon-Web Growth

Mathematical model represents thermomechanical behavior of solid web of rectangular cross section as pulled from melt of same material. Intended for application to growth of silicon webs or ribbons, model general enough to be applicable to other materials as well. Shows suitable temperature distribution results in stress-free ribbon. Conceptual ribbon of rectangular cross section pulled from melt. Simplified configuration used to calculate thermal conditions for stress-free ribbon.

Wade, Ben K.↗

Second-order closure models for supersonic turbulent flows

Recent work on the development of a second-order closure model for high-speed compressible flows is reviewed. This turbulent closure is based on the solution of modeled transport equations for the Favre-averaged Reynolds stress tensor and the solenoidal part of the turbulent dissipation rate. A new model for the compressible dissipation is used along with traditional gradient transport models for the Reynolds heat flux and mass flux terms. Consistent with simple asymptotic analyses, the deviatoric part of the remaining higher-order correlations in the Reynolds stress transport equations are modeled by a variable density extension of the newest incompressible models. The resulting second-order closure model is tested in a variety of compressible turbulent flows which include the decay of isotropic turbulence, homogeneous shear flow, the supersonic mixing layer, and the supersonic flat-plate turbulent boundary layer. Comparisons between the model predictions and the results of physical and numerical experiments are quite encouraging.

Speziale, Charles G.↗

Second-order closure models for supersonic turbulent flows

Recent work by the authors on the development of a second-order closure model for high-speed compressible flows is reviewed. This turbulence closure is based on the solution of modeled transport equations for the Favre-averaged Reynolds stress tensor and the solenoidal part of the turbulent dissipation rate. A new model for the compressible dissipation is used along with traditional gradient transport models for the Reynolds heat flux and mass flux terms. Consistent with simple asymptotic analyses, the deviatoric part of the remaining higher-order correlations in the Reynolds stress transport equation are modeled by a variable density extension of the newest incompressible models. The resulting second-order closure model is tested in a variety of compressible turbulent flows which include the decay of isotropic turbulence, homogeneous shear flow, the supersonic mixing layer, and the supersonic flat-plate turbulent boundary layer. Comparisons between the model predictions and the results of physical and numerical experiments are quite encouraging.

Speziale, Charles G.↗