Search NASA⌕ Search

SEARCH · Search NASA

Results for “STRESS CALCULATION”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

A nonlinear fracture mechanics approach to the growth of small cracks

An analytical model of crack closure is used to study the crack growth and closure behavior of small cracks in plates and at notches. The calculated crack opening stresses for small and large cracks, together with elastic and elastic plastic fracture mechanics analyses, are used to correlate crack growth rate data. At equivalent elastic stress intensity factor levels, calculations predict that small cracks in plates and at notches should grow faster than large cracks because the applied stress needed to open a small crack is less than that needed to open a large crack. These predictions agree with observed trends in test data. The calculations from the model also imply that many of the stress intensity factor thresholds that are developed in tests with large cracks and with load reduction schemes do not apply to the growth of small cracks. The current calculations are based upon continuum mechanics principles and, thus, some crack size and grain structure exist where the underlying fracture mechanics assumptions become invalid because of material inhomogeneity (grains, inclusions, etc.). Admittedly, much more effort is needed to develop the mechanics of a noncontinuum. Nevertheless, these results indicate the importance of crack closure in predicting the growth of small cracks from large crack data.

Newman, J. C., Jr.↗

Calculation of turbulence-driven secondary motion in ducts with arbitrary cross section

Calculation methods for turbulent duct flows are generalized for ducts with arbitrary cross-sections. The irregular physical geometry is transformed into a regular one in computational space, and the flow equations are solved with a finite-volume numerical procedure. The turbulent stresses are calculated with an algebraic stress model derived by simplifying model transport equations for the individual Reynolds stresses. Two variants of such a model are considered. These procedures enable the prediction of both the turbulence-driven secondary flow and the anisotropy of the Reynolds stresses, in contrast to some of the earlier calculation methods. Model predictions are compared to experimental data for developed flow in triangular duct, trapezoidal duct and a rod-bundle geometry. The correct trends are predicted, and the quantitative agreement is mostly fair. The simpler variant of the algebraic stress model procured better agreement with the measured data.

Demuren, A. O.↗

Calculation of turbulence-driven secondary motion in ducts with arbitrary cross-section

Calculation methods for turbulent duct flows are generalized for ducts with arbitrary cross-sections. The irregular physical geometry is transformed into a regular one in computational space, and the flow equations are solved with a finite-volume numerical procedure. The turbulent stresses are calculated with an algebraic stress model derived by simplifying model transport equations for the individual Reynolds stresses. Two variants of such a model are considered. These procedures enable the prediction of both the turbulence-driven secondary flow and the anisotropy of the Reynolds stresses, in contrast to some of the earlier calculation methods. Model predictions are compared to experimental data for developed flow in triangular duct, trapezoidal duct and a rod-bundle geometry. The correct trends are predicted, and the quantitative agreement is mostly fair. The simpler variant of the algebraic stress model procured better agreement with the measured data.

Demuren, A. O.↗

The axisymmetric elasticity problem for a laminated plate containing a circular hole

The elasticity problem for a laminated thick plate which consists of two bonded dissimilar layers and which contains a circular hole is considered. The problem is formulated for arbitrary axisymmetric tractions on the hole surface by using the Love strain function. Through the expansion of the boundary conditions into Fourier series the problem is reduced to an infinite system of algebraic equations which is solved by the method of reduction. Of particular interest in the problem are the stresses along the interface as they relate to the question of delamination failure of the composite plate. These stresses are calculated and are observed to become unbounded at the hole boundary. An approximate treatment of the singular behavior of the stress state is presented and the stress intensity factors are calculated.

Delale, F.↗

Elevated temperature compressive steady state deformation and failure in the oxide dispersion strengthened alloy MA 6000E

The compressive flow strength-strain rate behavior of the oxide-dispersion-strengthened alloy MA 6000E has been studied in the temperature range 1144-1366 K, with strain rates ranging from 2.1 x 10 to the -5th to 2.1 x 10 to the -7th per s. It is found that the inherent strength of the alloy is essentially the same in all test directions and that the low strength observed in tensile tests results from the inability of grain boundaries to support high tensile stresses. The failure of MA 6000E under high-temperature, slow plastic flow conditions is shown to be the result of concentrated slip. Slow plastic deformation in MA 6000E can be described by a threshold stress model of creep where threshold stresses are calculated from relatively fast testing procedures and the effective stress exponent for creep is assumed to be 3.5.

Whittenberger, J. D.↗

Analysis of silicon stress/strain relationships

In the study of stress-strain relationships in silicon ribbon, numerous solutions were calculated for stresses, strain rates, and dislocation densities through the use of the Sumino model. It was concluded that many cases of failure of computer solutions to converge are analytical manifestations of shear bands (Luder's band) observed in experiments.

Dillon, O.↗

Analysis of an Interface Crack for a Functionally Graded Strip Sandwiched between Two Homogeneous Layers of Finite Thickness

The interface crack problem for a composite layer that consists of a homogeneous substrate, coating and a non-homogeneous interface was formulated for singular integral equations with Cauchy kernels and integrated using the Lobatto-Chebyshev collocation technique. Mixed-mode Stress Intensity Factors and Strain Energy Release Rates were calculated. The Stress Intensity Factors were compared for accuracy with relevant results previously published. The parametric studies were conducted for the various thickness of each layer and for various non-homogeneity ratios. Particular application to the Zirconia thermal barrier on steel substrate is demonstrated.

Shbeeh, N. I.↗

Nonlinear temperature dependent failure analysis of finite width composite laminates

A quasi-three dimensional, nonlinear elastic finite element stress analysis of finite width composite laminates including curing stresses is presented. Cross-ply, angle-ply, and two quasi-isotropic graphite/epoxy laminates are studied. Curing stresses are calculated using temperature dependent elastic properties that are input as percent retention curves, and stresses due to mechanical loading in the form of an axial strain are calculated using tangent modulii obtained by Ramberg-Osgood parameters. It is shown that curing stresses and stresses due to tensile loading are significant as edge effects in all types of laminate studies. The tensor polynomial failure criterion is used to predict the initiation of failure. The mode of failure is predicted by examining individual stress contributions to the tensor polynomial.

Nagarkar, A. P.↗

Thermal residual stresses in silicon-carbide/titanium (0/90) laminate

The current work formulated a micromechanical analysis of a cross-ply laminate and calculated the thermal residual stress in a very thick (0/90)(sub 2n) silicon-carbide/titanium laminate. Results were also shown for a unidirectional laminate of the same material. Discrete fiber-matrix models assuming a rectangular array of fibers with a fiber volume fraction of 32.5 percent and a three-dimensional, finite-element analysis were used. Significant differences in the trends and magnitudes for the fiber, matrix, and interface stresses were calculated for unidirectional and (0/90) models. Larger hoop stresses calculated for the (0/90) model indicate that it may be more susceptible to radial cracking when subjected to mechanical loading than the unidirectional model. The axial stresses in the matrix were calculated to be slightly larger for the (0/90) model. The compressive axial stresses in the fiber were significantly larger in the (0/90) model. The presence of the cross-ply in the (0/90) model reduced the constraint on the fiber, producing radial interface stresses that were less compressive, which could lead to earlier failure of the fiber-matrix interface.

Bigelow, C. A.↗

Sensitivity of Lumped Constraints Using the Adjoint Method

Adjoint sensitivity calculation of stress, buckling and displacement constraints may be much less expensive than direct sensitivity calculation when the number of load cases is large. Adjoint stress and displacement sensitivities are available in the literature. Expressions for local buckling sensitivity of isotropic plate elements are derived in this study. Computational efficiency of the adjoint method is sensitive to the number of constraints and, therefore, the method benefits from constraint lumping. A continuum version of the Kreisselmeier-Steinhauser (KS) function is chosen to lump constraints. The adjoint and direct methods are compared for three examples: a truss structure, a simple HSCT wing model, and a large HSCT model. These sensitivity derivatives are then used in optimization.

Akgun, Mehmet A.↗

Stress-concentration factors for finite orthotropic laminates with a pin-loaded hole

Stresses were calculated for finite size orthotropic laminates loaded by a frictionless steel pin in a circular hole of the same diameter. The calculations were based on finite element analyses for six laminates. Stress concentration factors, based on nominal bearing stress, were determined for wide ranges of the ratios of width to diameter, w/d and edge distance to diameter, e/d. An infinite laminate case was analyzed for each laminate. Orthotropy had a significant influence on the tensile stress concentration at the hole. For example, the stress concentration factors for the infinite laminate cases ranged from 0.82 to 2.16, compared with 0.87 for the quasi-isotropic laminate. The finite widths and edge distances strongly influenced the tensile stress concentration. For the practical range w/d or = 3, the peak tensile stresses were as much as 80% larger than the infinite laminate reference value. For e/d or = 3, these stresses were amplified by as much as 50%. In contrast, the finite width and edge distance had little effect on shear-out and bearing stress concentrations.

Crews, J. H., Jr.↗

Prediction of Fretting Crack Location and Orientation in a Single Crystal Nickel Alloy

Fretting is a structural damage mechanism arising between two nominally clamped surfaces subjected to an oscillatory loading. A critical location for fretting induced damage has been identified at the blade/disk and blade/damper interfaces of gas turbine engine turbomachinery and space propulsion components. The high- temperature, high-frequency loading environment seen by these components lead to severe stress gradients at the edge-of-contact that could potentially foster crack growth leading to component failure. These contact stresses drive crack nucleation in fretting and are very sensitive to the geometry of the contacting bodies, the contact loads, materials, temperature, and contact surface tribology (friction). Recently, a high-frequency, high-temperature load frame has been designed for experimentally investigating fretting damage of single crystal nickel materials employed in aircraft and spacecraft turbomachinery. A modeling method for characterizing the fretting stresses of the spherical fretting contact stress behavior in this experiment is developed and described. The calculated fretting stresses for a series of experiments are then correlated to the observed fretting damage. Results show that knowledge of the normal stresses and resolved shear stresses on each crystal plane can aid in predicting crack locations and orientations.

Matlik, J. F.↗

TURBOMAT-FR: Turbomachinery Aeroelastic Forced Response Analysis Automation Using MATLAB®: PART 1

The primary aim of this effort is to develop a software tool named “Turbomachinery Aeroelastic Analysis Tool for Forced Response (TURBOMAT-FR)” using MATLAB® (The MathWorks, Inc.) scripts. The tool will aid in performing routine forced response analysis of turbomachinery fan blades with the aid of the modal summation method (MSUM). The MSUM approach requires a modal solution consisting of modal displacements, modal frequencies, and modal stresses. The MATLAB® scripts facilitate easy extraction of the required quantities from the computational fluid dynamics (CFD) and MSC Nastran (Hexagon AB) solution files. These quantities are subsequently utilized by the MSUM approach to perform a complete forced response analysis.The MATLAB® scripts in TURBOMAT-FR are utilized in postprocessing of the MSC Nastran finite element analysis outputs in three types of solutions or steps, namely (1) stresses and displacements from static and steady response solution, (2) unsteady stresses and displacements from direct forced response solution, and (3) modal solution containing modal displacements and modal stresses. In step (4), separate MATLAB® scripts are also developed to read the processed outputs in step (3) and perform forced response analysis using MSUM. The TURBOMAT-FR MATLAB® scripts mentioned previously are presented in two reports designated “Part 1” and “Part 2”. The current report, Part 1, presents the scripts needed to accomplish steps (1) and (2), noted previously. This will be exclusively useful for performing only the static and steady response analysis and the forced response analysis using the direct forced response solution approach. Part 2 contains the TURBOMAT-FR MATLAB® scripts needed to accomplish steps (3) and (4). The MATLAB® scripts given in Part 2 facilitate the calculation of unsteady stresses using MSUM and provide an automated procedure associated with aeroelastic forced response analysis with different pressure files for the same structural model. The unsteady stresses calculated using MSUM in Part 2 can be compared with those presented in Part 1. A cantilevered blade-alone model subjected to a realistic pressure loading to a boundary was used for static and unsteady direct forced response analysis to demonstrate the effectiveness of the TURBOMATFR MATLAB® scripts.

MATLAB-Scripts↗

TURBOMAT-FR: Turbomachinery Aeroelastic Forced Response Analysis Automation Using MATLAB® Part 2

The primary aim of this effort is to develop a software tool named “Turbomachinery Aeroelastic Analysis Tool for Forced Response (TURBOMAT-FR)” using MATLAB® (The MathWorks, Inc.) scripts. The tool will aid in performing routine forced response analysis of turbomachinery fan blades with the aid of the modal summation method (MSUM). The MSUM approach requires a modal solution consisting of modal displacements, frequencies, and stresses. The MATLAB® scripts facilitate easy extraction of the required quantities from the computational fluid dynamics (CFD) and MSC Nastran (Hexagon AB) solution files. These quantities are subsequently utilized by the MSUM approach to perform a complete forced response analysis using the MATLAB® scripts. The MATLAB® scripts in TURBOMAT-FR are utilized in postprocessing of the MSC Nastran finite element analysis outputs in three types of solutions or steps, namely (1) stresses and displacements from static or steady response solution, (2) unsteady stresses and displacements from direct forced response solution, and (3) modal solution containing modal displacements and modal stresses. In step (4), separate MATLAB® scripts are also developed to read the processed outputs in step (3) and perform forced response analysis using MSUM. It should be noted that the direct forced response solution from step (2) will provide for a comparison with those obtained from the MSUM approach in steps (3) and (4). The TURBOMAT-FR MATLAB® scripts mentioned previously are presented in two reports designated Parts 1 and 2. Part 1 presents the scripts needed to accomplish steps (1) and (2). The current report, Part 2, presents the scripts needed to accomplish steps (3) and (4), noted previously. The TURBOMAT-FR MATLAB® scripts provided here in Part 2 facilitate the calculation of unsteady stresses using MSUM and provide an automated procedure associated with aeroelastic forced response analysis with different pressure files. The unsteady stresses calculated using MSUM can be compared with those presented in Part 1. A cantilevered blade-alone model subjected to a realistic pressure loading to a boundary was used for static and unsteady direct forced response analysis to demonstrate the effectiveness of the TURBOMATFR MATLAB® scripts.

T S R Reddy↗

Analysis of a Generally Oriented Crack in a Functionally Graded Strip Sandwiched Between Two Homogeneous Half Planes

The driving forces for a generally oriented crack embedded in a Functionally Graded strip sandwiched between two half planes are analyzed using singular integral equations with Cauchy kernels, and integrated using Lobatto-Chebyshev collocation. Mixed-mode Stress Intensity Factors (SIF) and Strain Energy Release Rates (SERR) are calculated. The Stress Intensity Factors are compared for accuracy with previously published results. Parametric studies are conducted for various nonhomogeneity ratios, crack lengths. crack orientation and thickness of the strip. It is shown that the SERR is more complete and should be used for crack propagation analysis.

Shbeeb, N.↗

Contact and Bending Durability Calculation for Spiral-Bevel Gears

The objective of this project is to extend the capabilities of the gear contact analysis solver Calyx, and associated packages Transmission3D, HypoidFaceMilled, HypoidFaceHobbed. A calculation process for the surface durability was implemented using the Dowson-Higginson correlation for fluid film thickness. Comparisons to failure data from NASA's Spiral Bevel Gear Fatigue rig were carried out. A bending fatigue calculation has been implemented that allows the use of the stress-life calculation at each individual fillet point. The gears in the NASA test rig did not exhibit any bending fatigue failure, so the bending fatigue calculations are presented in this report by using significantly lowered strength numbers.

life↗