Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Strain 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 739 records · Page 41

Evolution of a moderately short turbulent boundary layer in a severe pressure gradient

The early and intermediate development of a highly accelerated (or decelerated) turbulent boundary layer is analyzed. For sufficiently large accelerations (or pressure gradients) and for total normal strains which are not excessive, the equation for the Reynolds shear stress simplifies to give a stress that remains approximately constant as it is convected along streamlines. The theoretical results for the evolution of the mean velocity in favourable and adverse pressure gradients agree well with experiment for the cases considered. A calculation which includes mass injection at the wall is also given.

Deissler, R. G.↗

A triangular thin shell finite element: Nonlinear analysis

Aspects of the formulation of a triangular thin shell finite element which pertain to geometrically nonlinear (small strain, finite displacement) behavior are described. The procedure for solution of the resulting nonlinear algebraic equations combines a one-step incremental (tangent stiffness) approach with one iteration in the Newton-Raphson mode. A method is presented which permits a rational estimation of step size in this procedure. Limit points are calculated by means of a superposition scheme coupled to the incremental side of the solution procedure while bifurcation points are calculated through a process of interpolation of the determinants of the tangent-stiffness matrix. Numerical results are obtained for a flat plate and two curved shell problems and are compared with alternative solutions.

Thomas, G. R.↗

A single fracture toughness parameter for fibrous composite laminates

A general fracture toughness parameter Qc was previously derived and verified to be a material constant, independent of layup, for centrally cracked boron aluminum composite specimens. The specimens were made with various proportions of 0 and + or - 45 degree plies. A limited amount of data indicated that the ratio Qc/epsilon tuf' where epsilon tuf is the ultimate tensile strain of the fibers, might be a constant for all composite laminates, regardless of material and layup. In that case, a single value of Qc/epsilon tuf could be used to predict the fracture toughness of all fibrous composite laminates from only the elastic constants and epsilon tuf. Values of Qc/epsilon tuf were calculated for centrally cracked specimens made from graphite/polyimide, graphite/epoxy, E glass/epoxy, boron/epoxy, and S glass graphite/epoxy materials with numerous layups. Within ordinary scatter, the data indicate that Qc/epsilon tuf is a constant for all laminates that did not split extensively at the crack tips or have other deviate failure modes.

Poe, C. C., Jr.↗

Numerical and experimental simulations of the postbuckling response of laminated anisotropic panels

Numerical simulations are used in conjuction with experiments to study the buckling and postbuckling responses and failure initiation of flat unstiffened composite panels. The numerical simulations are conducted using two-dimensional shear-flexible finite elements. The effect of the laminate stacking sequence on the buckling and postbuckling responses is studied. Correlation between numerical and experimental results is good through buckling, but the numerical models overestimate the postbuckling stiffness of the panels. To explain the discrepancies in the postbuckling stiffnesses, analytic sensitivity derivatives are calculated and are used to study the sensitivity of the buckling and postbuckling responses to variations in different material and lamination parameters of the panel. Experimental results indicate that failure occurs along a nodal line. Numerical results show that the location of failure initiation corresponds to that of the maximum transverse shear strain energy density in the panel, which occurs at the edge of the panel at a nodal line.

Noor, Ahmed K.↗

Analysis of notched metal matrix composites under tensile loading

This paper presents techniques based on a three-dimensional finite-element analysis for the analysis of continuous fiber reinforced metal matrix composite. Examples are shown for specific metal matrix composites such as boron/aluminum and silicon-carbide/aluminum. Specimen stress-strain behavior and stress at first fiber failure were predicted for boron/aluminum laminates containing circular holes and crack-like slits. The predictions compared very well with test data for laminates containing 0 deg fibers and reasonably well for (+/-45)2s laminates. Mesh configuration was shown to have an effect on the calculation of stresses local to the notch. The presence of thin interface layers of matrix material had a significant influence on the slit tip stress state, causing sharper stress gradients near the notch. Interface layers reduced the slit-tip fibers stress in a (+/-45)2s silicon-carbide/aluminum laminate but increased them in a (0/90)2s laminate.

Bigelow, C. A.↗

Diffraction Studies of the Atomic Vibrations of Bulk and Surface Atoms in the Reciprocal and Real Spaces of Nanocrystalline SiC

To describe and evaluate the vibrational properties of nanoparticles it is necessary to distinguish between the surface and the core of the particles. Theoretical calculations show that vibrational density of states of the inner atoms of nanograins is similar to bulk material but shifted to higher energies which can be explained by the fact that the gain core is stressed (hardened) due to the presence of internal pressure. Theoretical calculations also show that there is a difference between vibrational properties of a crystal lattice of the grain interior in isolated particles and in a dense (sintered) nanocrystalline material. This is probably due to a coupling of the modes inside the grains via the grain boundaries in dense nanocrystalline bodies. We examined strains present in the surface shell based on examination of diamond and Sic nanocrystals in reciprocal (Bragg-type scattering) and real (PDF analysis) space analysis of neutron diffraction data. Recently we examined the atomic thermal motions in nanocrystalline Sic based on the assumption of a simple Einstein model for uncorrelated atomic notions. According to this model, the Bragg intensity is attenuated as a function of scattering angle by the Debye-Waller factor. Based on this assumption overall temperature factors were determined from the Wilson plots.

Stelmakh, S.↗

Calculation of Dynamic Loads Due to Random Vibration Environments in Rocket Engine Systems

An important part of rocket engine design is the calculation of random dynamic loads resulting from internal engine "self-induced" sources. These loads are random in nature and can greatly influence the weight of many engine components. Several methodologies for calculating random loads are discussed and then compared to test results using a dynamic testbed consisting of a 60K thrust engine. The engine was tested in a free-free condition with known random force inputs from shakers attached to three locations near the main noise sources on the engine. Accelerations and strains were measured at several critical locations on the engines and then compared to the analytical results using two different random response methodologies.

Christensen, Eric R.↗

Stress wave calculations in composite plates using the fast Fourier transform.

The protection of composite turbine fan blades against impact forces has prompted the study of dynamic stresses in composites due to transient loads. The mathematical model treats the laminated plate as an equivalent anisotropic material. The use of Mindlin's approximate theory of crystal plates results in five two-dimensional stress waves. Three of the waves are flexural and two involve in-plane extensional strains. The initial value problem due to a transient distributed transverse force on the plate is solved using Laplace and Fourier transforms. A fast computer program for inverting the two-dimensional Fourier transform is used. Stress contours for various stresses and times after application of load are obtained for a graphite fiber-epoxy matrix composite plate. Results indicate that the points of maximum stress travel along the fiber directions.

Moon, F. C.↗

ANISAP: A three-dimensional finite element program for laminated composites subjected to mechanical loading

ANISAP is a 3-D finite element FORTRAN 77 computer code for linear elastic, small strain, analysis of laminated composites with arbitrary geometry including free edges and holes. Individual layers may be isotropic or transversely isotropic in material principal coordinates; individual layers may be rotated off-axis about a global z-axis. The laminate may be a hybrid. Three different isoparametric elements, variable order of gaussian integration, calculation of stresses at element boundaries, and loading by either nodal displacement of forces are included in the program capability. Post processing capability includes failure analysis using the tensor polynominal failure criterion.

Burns, S. W.↗

Thermal/Mechanical Durability of Polymer-Matrix Composites in Cryogenic Environments

In order to increase the reliability of the next generation of space transportation systems, the mechanical behavior of polymeric-matrix composite (PMC) materials at cryogenic temperatures must be investigated. This paper presents experimental data on the residual mechanical properties of a carbon fiber polymeric composite, IM7/PETI-5 as a function of temperature and aging. Tension modulus and strength were measured at room temperature, -196 C, and -269 C on five different specimens ply lay-ups. Specimens were preconditioned with one set of coupons being isothermally aged for 576 hours at -184 C, in an unloaded state. Another set of corresponding coupons were mounted in constant strain fixtures such that a constant uniaxial strain was applied to the specimens for 576 hours at -184 C. A third set was mechanically cycled in tension at -184 C. The measured properties indicated that temperature, aging, and loading mode can all have significant influence on performance. Moreover, this influence is a strong function of laminate stacking sequence. Thermal-stress calculations based on lamination theory predicted that the transverse tensile ply stresses could be quite high for cryogenic test temperatures. Microscopic examination of the surface morphology showed evidence of degradation along the exposed edges of the material because of aging at cryogenic temperatures. ________________

Gates, Thomas S.↗

Strain Gage Load Calibration of the Wing Interface Fittings for the Adaptive Compliant Trailing Edge Flap Flight Test

The safety-of-flight parameters for the Adaptive Compliant Trailing Edge (ACTE) flap experiment require that flap-to-wing interface loads be sensed and monitored in real time to ensure that the structural load limits of the wing are not exceeded. This paper discusses the strain gage load calibration testing and load equation derivation methodology for the ACTE interface fittings. Both the left and right wing flap interfaces were monitored; each contained four uniquely designed and instrumented flap interface fittings. The interface hardware design and instrumentation layout are discussed. Twenty-one applied test load cases were developed using the predicted in-flight loads. Pre-test predictions of strain gage responses were produced using finite element method models of the interface fittings. Predicted and measured test strains are presented. A load testing rig and three hydraulic jacks were used to apply combinations of shear, bending, and axial loads to the interface fittings. Hardware deflections under load were measured using photogrammetry and transducers. Due to deflections in the interface fitting hardware and test rig, finite element model techniques were used to calculate the reaction loads throughout the applied load range, taking into account the elastically-deformed geometry. The primary load equations were selected based on multiple calibration metrics. An independent set of validation cases was used to validate each derived equation. The 2-sigma residual errors for the shear loads were less than eight percent of the full-scale calibration load; the 2-sigma residual errors for the bending moment loads were less than three percent of the full-scale calibration load. The derived load equations for shear, bending, and axial loads are presented, with the calculated errors for both the calibration cases and the independent validation load cases.

Miller, Eric J.↗

Method for calculating alloy energetics

A semiempirical method for the computation of alloy energies is introduced. It is based on the equivalent-crystal theory of defect-formation energies in elemental solids. The method is both simple and accurate. Heats of formation as a function of composition are computed for some binary alloys of Cu, Ni, Al, Ag, Pd, Pt, and Au using the heats of solution in the dilute limit as experimental input. The separation of heats into strain and chemical components helps in understanding the energetics. In addition, lattice-parameter contractions seen in solid solutions of Ag and Au are accurately predicted. Good agreement with experiment is obtained in all cases.

Bozzolo, Guillermo↗

Fracture modes of high modulus graphite/epoxy angleplied laminates subjected to off-axis tensile loads

Angleplied laminates of high modulus graphite fiber/epoxy were examined in several ply configurations at various tensile loading angles to the zero ply direction to determine the effects of ply orientations on tensile properties, fracture modes, and fracture surface characteristics of the various plies. Experimental results consist of stress-strain data, selected plots, fracture stresses and strains, and scanning electron microscope (SEM) photographs of fracture surfaces. It was found that the stress-strain curves were linear to fracture, and that although fracture surface characteristics for a given fracture mode are similar to those for the same fracture mode in uniaxial specimens, no simple load angle range can be associated with a given fracture mode. It was also concluded that SEM results must be supplemented with ply stress calculations in order to identify ranges of fracture modes occurring as a function of ply orientation with respect to the load direction.

Sinclair, J. H.↗

Setting the Bar for the Replacement of the Probability of Collision Metric

To date, satellite conjunction assessment (CA) risk analysis has largely embraced the probability of collision (Pc) as the omnibus metric to evaluate collision likelihood, and its use in such assessments has mostly been straightforward: at the point at which a conjunction mitigation decision is required, the calculated Pc is compared to a threshold; and if the calculated Pc exceeds that threshold, then a mitigation action is warranted. With only minor variation, this approach is employed by major CA risk assessment centers (e.g., NASA, EUSST, CNES, JAXA) and is advanced as the preferred method in the published CA best practices handbooks. Despite this near unanimity of operational practice, there is a major strain of secondary literature critical of the Pc and willing to propose alternatives. Alfano (2005) pointed out the ability of the Pc to underrepresent the risk in certain situations and counselled a maximum Pc construct. Carpenter (2017, 2019) reiterated this criticism and proposed using instead a confidence interval on the miss distance. Balch et al. (2019) identified what they argued was a defect in the entire Bayesian Pc construct and believed that the use of a more conservative methodology based on covariance ellipsoid overlap was necessary. Delande (2022) introduced the framework of collision “plausibility” to the risk assessment process and sketched out how this might be used operationally. Elkantassi (2022) published a full development of the miss distance confidence interval approach and applied it to several worked examples. While these different approaches to collision risk assessment do differ in their details, they all converge on two central points: first, the Pc’s failure to give an adequate expression of the risk in dilution region situations is a fatal flaw; and second, a conjunction should be presumed risky and in need of mitigation until the evidence of the situation can establish otherwise. These criticisms, if correct, would counsel a number of modifications to current CA operational practice; as such, they force a re-examination of fundamental aspects of the CA problem, including the following: 1. Is the CA risk assessment a probability problem, a statistics problem, or something else? 2. If it is a statistics problem, does it lend itself naturally to a hypothesis test construction? 3. If it can be construed as a hypothesis test, what form should the null hypothesis take, to wit: what constraints exist on the choice of the null hypothesis, what selections are in best alignment with all of the attendant parameters of the problem, and what is implied philosophically by different choices? 4. What are the implications of using the different proposed risk assessment parameters for CA? This question should be answered both in determining how frequently the dilution region situation cited by the critics of the Pc actually appears in an operationally significant manner and the missed detection and false alarm rates of all of the proposed risk assessment metrics, compared both to the Pc and to each other. This paper explores and offers preliminary answers to the above questions, presenting a researched treatment of the philosophical nature of the CA problem and the null hypothesis choice that achieves the greatest consistency with all of the different aspects of operational CA conduct. It then profiles all of the different proposed risk assessment metrics enumerated in the earlier paragraph against an extremely large database of conjunction events at both the 550km and 700km altitudes. The combination of the philosophical exploration of the CA problem and the results of the profiling activity articulates what a risk assessment metric will need to demonstrate, in terms of both innate construction and performance, in order to be a true competitor to the Pc.

conjunction assessment↗

Development of procedures for calculating stiffness and damping properties of elastomers in engineering applications. Part 2: Elastomer characteristics at constant temperature

Dynamic properties of a commerical polybutadiene compound were determined at a constant temperature of 32 C by a forced-vibration resonant mass type of apparatus. The constant thermal state of the elastomer was ensured by keeping the ambient temperature constant and by limiting the power dissipation in the specimen. Experiments were performed with both compression and shear specimens at several preloads (nominal strain varying from 0 to 5 percent), and the results are reported in terms of a complex stiffness as a function of frequency. Very weak frequency dependence is observed and a simple power law type of correlation is shown to represent the data well. Variations in the complex stiffness as a function of preload are also found to be small for both compression and shear specimens.

Gupta, P. K.↗

An X-ray photoelectron spectroscopy study of the thermal nitridation of SiO2/Si

The dependence of the nitrogen distribution in thermally nitrided SiO2 films on the nitridation time and temperature has been studied by means of X-ray photoelectron spectroscopy (XPS). The photoelectron peak intensities were measured by fitting Voigt profiles to the XPS spectra and were used to calculate the film composition as a function of film depth, applying an analytical method described in detail. The times of appearance of the maxima in interfacial nitrogen concentration are shown for 800, 1000, and 1150 C, and the data are related to a kinetic model of Vasquez and Madhukar (1985), which considers the effect of interfacial strain on the nitridation kinetics. In addition, the intensity of a fluorine marker (from the HF used in the etching step) was found to correlate with the nitrogen concentration. It is postulated that the F bonds preferentially to defects. This hypothesis and the measured F intensities are consistent with the proposed strain-dependent energy of defect formation.

Vasquez, R. P.↗

LES/RANS Simulation of a Supersonic Reacting Wall Jet

This work presents results from large-eddy / Reynolds-averaged Navier-Stokes (LES/RANS) simulations of the well-known Burrows-Kurkov supersonic reacting wall-jet experiment. Generally good agreement with experimental mole fraction, stagnation temperature, and Pitot pressure profiles is obtained for non-reactive mixing of the hydrogen jet with a non-vitiated air stream. A lifted flame, stabilized between 10 and 22 cm downstream of the hydrogen jet, is formed for hydrogen injected into a vitiated air stream. Flame stabilization occurs closer to the hydrogen injection location when a three-dimensional combustor geometry (with boundary layer development resolved on all walls) is considered. Volumetric expansion of the reactive shear layer is accompanied by the formation of large eddies which interact strongly with the reaction zone. Time averaged predictions of the reaction zone structure show an under-prediction of the peak water concentration and stagnation temperature, relative to experimental data and to results from a Reynolds-averaged Navier-Stokes calculation. If the experimental data can be considered as being accurate, this result indicates that the present LES/RANS method does not correctly capture the cascade of turbulence scales that should be resolvable on the present mesh. Instead, energy is concentrated in the very largest scales, which provide an over-mixing effect that excessively cools and strains the flame. Predictions improve with the use of a low-dissipation version of the baseline piecewise parabolic advection scheme, which captures the formation of smaller-scale structures superimposed on larger structures of the order of the shear-layer width.

Edwards, Jack R.↗

An ILLIAC program for the numerical simulation of homogeneous incompressible turbulence

An algorithm and ILLIAC computer program, developed for the simulation of homogeneous incompressible turbulence in the presence of an applied mean strain, are described. The turbulence field is represented spatially by a truncated triple Fourier series (spectral method) and followed in time using a fourth-order Runge-Kutta algorithm. These include: (1) transformation of variables suggested by Taylor's sudden-distortion theory; (2) implicit viscous diffusion by use of an integrating factor; (3) implicit pressure calculation suggested by Taylor's sudden-distortion theory, and (4) inexpensive control of aliasing by random and phased coordinate shifts.

Rogallo, R. S.↗