SEARCH · Engineering Papers
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.
The advantages of a newly developed code are demonstrated by comparisons of the analysis with existing theoretical data as well as with other available finite element programs. The new program shows promise to significantly reduce the computer time. It also permits accurate and efficient structural analysis of engine hot section components.
The structural performance of a boron-epoxy-reinforced titanium drag strut, containing a bonded scarf joint and designed to the criteria of a large commercial transport, has been evaluated experimentally and analytically. The strut was exposed to two lifetimes of fatigue loading and was statically loaded to the tensile and compressive design ultimate loads. Throughout the test program no evidence of any damage in the drag strut was detected by strain-gage measurements, ultrasonic inspection, or visual observation. The bonded joint was analyzed using the NASTRAN computer program. A comparison of the strains predicted by the NASTRAN computer program with the experimentally determined values shows excellent agreement. An analytical study indicated that the nonlinear behavior of a structural spacer at each end of the strut could be explained by the inelastic behavior and possible creep of the adhesive.
General simulation equations are derived for the rigid body motion of slung load systems. These systems are viewed as consisting of several rigid bodies connected by straight-line cables or links. The suspension can be assumed to be elastic or inelastic, both cases being of interest in simulation and control studies. Equations for the general system are obtained via D'Alembert's principle and the introduction of generalized velocity coordinates. Three forms are obtained. Two of these generalize previous case-specific results for single helicopter systems with elastic or inelastic suspensions. The third is a new formulation for inelastic suspensions. It is derived from the elastic suspension equations by choosing the generalized coordinates so as to separate motion due to cable stretching from motion with invariant cable lengths. The result is computationally more efficient than the conventional formulation, and is readily integrated with the elastic suspension formulation and readily applied to the complex dual lift and multilift systems. Equations are derived for dual lift systems. Three proposed suspension arrangements can be integrated in a single equation set. The equations are given in terms of the natural vectors and matrices of three-dimensional rigid body mechanics and are tractable for both analysis and programming.
Since most advanced material systems (for example metallic-, polymer-, and ceramic-based systems) being currently researched and evaluated are for high-temperature airframe and propulsion system applications, the required constitutive models must account for both reversible and irreversible time-dependent deformations. Furthermore, since an integral part of continuum-based computational methodologies (be they microscale- or macroscale-based) is an accurate and computationally efficient constitutive model to describe the deformation behavior of the materials of interest, extensive research efforts have been made over the years on the phenomenological representations of constitutive material behavior in the inelastic analysis of structures. From a more recent and comprehensive perspective, the NASA Glenn Research Center in conjunction with the University of Akron has emphasized concurrently addressing three important and related areas: that is, 1) Mathematical formulation; 2) Algorithmic developments for updating (integrating) the external (e.g., stress) and internal state variables; 3) Parameter estimation for characterizing the model. This concurrent perspective to constitutive modeling has enabled the overcoming of the two major obstacles to fully utilizing these sophisticated time-dependent (hereditary) constitutive models in practical engineering analysis. These obstacles are: 1) Lack of efficient and robust integration algorithms; 2) Difficulties associated with characterizing the large number of required material parameters, particularly when many of these parameters lack obvious or direct physical interpretations.
Three dimensional, nonlinear finite element structural analyses were performed for a simulated combustor liner specimen to assess the capability of nonlinear analyses using classical inelastic material models to represent the thermoplastic creep response of the one half scale component. Results indicate continued cyclic hardening and ratcheting while experimental data suggested a stable stress strain response after only a few loading cycles. The computed stress strain history at the critical location was put into two life prediction methods, strainrange partitioning and a Pratt and Whitney combustor life prediction method to evaluate their ability to predict cyclic crack initiation. It is found that the life prediction analyses over predicted the observed cyclic crack initiation life.
Using the concept of an internal time as related to plastic strains, a differential stress-strain relation for elastoplasticity is rederived, such that (1) the concept of a yield-surface is retained; (2) the definitions of elastic and plastic processes are analogous to those in classical plasticity theory; and (3) its computational implementation, via a 'tangent-stiffness' finite element method and a 'generalized-midpoint-radial-return' stress-integration algorithm, is simple and efficient. Also, using the concept of an internal time, as related to both the inelastic strains as well as the Newtonian time, a constitutive model for creep-plasticity interaction, is discussed. The problem of modeling experimental data for plasticity and creep, by the present analytical relations, as accurately as desired, is discussed. Numerical examples which illustrate the validity of the present relations are presented for the cases of cyclic plasticity and creep.
This paper presents an assessment of the solution procedures available for the analysis of inelastic and/or large deflection structural behavior. A literature survey is given which summarized the contribution of other researchers in the analysis of structural problems exhibiting material nonlinearities and combined geometric-material nonlinearities. Attention is focused at evaluating the available computation and solution techniques. Each of the solution techniques is developed from a common equation of equilibrium in terms of pseudo forces. The solution procedures are applied to circular plates and shells of revolution in an attempt to compare and evaluate each with respect to computational accuracy, economy, and efficiency. Based on the numerical studies, observations and comments are made with regard to the accuracy and economy of each solution technique.
The focus of this dissertation is on advanced development of the boundary element method for elastic and inelastic thermal stress analysis. New formulations for the treatment of body forces and nonlinear effects are derived. These formulations, which are based on particular integral theory, eliminate the need for volume integrals or extra surface integrals to account for these effects. The formulations are presented for axisymmetric, two and three dimensional analysis. Also in this dissertation, two dimensional and axisymmetric formulations for elastic and inelastic, inhomogeneous stress analysis are introduced. The derivatives account for inhomogeneities due to spatially dependent material parameters, and thermally induced inhomogeneities. The nonlinear formulation of the present work are based on an incremental initial stress approach. Two inelastic solutions algorithms are implemented: an iterative; and a variable stiffness type approach. The Von Mises yield criterion with variable hardening and the associated flow rules are adopted in these algorithms. All formulations are implemented in a general purpose, multi-region computer code with the capability of local definition of boundary conditions. Quadratic, isoparametric shape functions are used to model the geometry and field variables of the boundary (and domain) of the problem. The multi-region implementation permits a body to be modeled in substructured parts, thus dramatically reducing the cost of analysis. Furthermore, it allows a body consisting of regions of different (homogeneous) material to be studied. To test the program, results obtained for simple test cases are checked against their analytic solutions. Thereafter, a range of problems of practical interest are analyzed. In addition to displacement and traction loads, problems with body forces due to self-weight, centrifugal, and thermal loads are considered.
The collision terms in the Boltzmann equation associated with various processes are derived. For processes having a Fokker-Planck (F-P) limit, the associated F-P operator is derived by means of physical arguments to determine the form of the operator; its multiplying constant is fixed by calculating the total energy exchange rate and comparing with the rate computed by other means. In this manner, the F-P operator is derived for electron-ion scattering, electron-electron scattering in the high-velocity limit, electron-atom elastic scattering, Compton scattering, and the high-velocity limit for inelastic scattering. Other processes considered are bremsstrahlung, radiative-recombination, photoionization, collisional ionization of atoms, and suprathermal-particle ionization of atoms.
Two nonintrusive flowfield diagnostics based on spectrally-resolved elastic (Rayleigh) and inelastic (Raman) laser light scattering were developed for obtaining local flowfield measurements in low-thrust gaseous H2/O2 rocket engines. The objective is to provide an improved understanding of phenomena occurring in small chemical rockets in order to facilitate the development and validation of advanced computational fluid dynamics (CFD) models for analyzing engine performance. The laser Raman scattering diagnostic was developed to measure major polyatomic species number densities and rotational temperatures in the high-density flowfield region extending from the injector through the chamber throat. Initial application of the Raman scattering diagnostic provided O2 number density and rotational temperature measurements in the exit plane of a low area-ratio nozzle and in the combustion chamber of a two-dimensional, optically-accessible rocket engine. In the low-density nozzle exit plane region where the Raman signal is too weak, a Doppler-resolved laser Rayleigh scattering diagnostic was developed to obtain axial and radial mean gas velocities, and in certain cases, H2O translational temperature and number density. The results from these measurements were compared with theoretical predictions from the RPLUS CFD code for analyzing rocket engine performance. Initial conclusions indicate that a detailed and rigorous modeling of the injector is required in order to make direct comparisons between laser diagnostic measurements and CFD predictions at the local level.
Neutrino, gamma-ray, electron and positron production resulting from inelastic proton-proton collisions in a highly relativistic plasma such as may exist in extragalactic radio or gamma-ray burst sources is examined. The source functions of primary (pions, kaons, and neutrons) and secondary (photons, electrons, positrons and neutrinos) products of relativistic nuclear collisions are computed for the cases of power law and Maxwellian relativistic proton distributions. It is shown that in plasma which is optically thin to interactions between the plasma and secondary gamma-rays, electrons and positrons, only a small fraction of the initial energy is emitted in the observable form of neutrinos and gamma rays. In an optically thick plasma on the other hand, most of the energy of the relativistic protons is found to be equally divided between gamma rays and neutrinos, although only the neutrinos may escape freely to be observed.
The results are presented of a viscoplastic stress-strain analysis of a cylindrical thrust chamber used in experiments that were designed to simulate the SSME operating conditions. The inelastic strain was calculated by using a viscoplastic model, and a quasi-three-dimensional structural analysis was performed by using a finite element program MARC. The temperatures and pressures were calculated by using the loading cycles of experiments, and these temperatures and pressures were used in the computations. The deformed shape of the component was predicted after the end of each loading cycle. The predicted shape qualitatively replicated the deformed shape of the component as observed in experiments. The results indicate that the use of viscoplastic models for structural analysis may lead to more realistic life assessments of experimental thrust chambers.
It is often advantageous to account for the microstructure of the material directly using multiscale modeling. For computational tractability, an idealized repeating unit cell (RUC) is used to capture all of the pertinent features of the microstructure. Typically, the RUC is dimensionless and depends only on the relative volume fractions of the different phases in the material. This works well for non-linear and inelastic behavior exhibiting a positive-definite constitutive response. Although, once the material exhibits strain softening, or localization, a mesh objective failure theories, such as smeared fracture theories, nodal and element enrichment theories (XFEM), cohesive elements or virtual crack closure technique (VCCT), can be utilized at the microscale, but the dimensions of the RUC must then be defined. One major challenge in multiscale progressive damage modeling is relating the characteristic lengths across the scales in order to preserve the energy that is dissipated via localization at the microscale. If there is no effort to relate the size of the macroscale element to the microscale RUC, then the energy that is dissipated will remain mesh dependent at the macroscale, even if it is regularized at the microscale. Here, a technique for mapping characteristic lengths across the scales is proposed. The RUC will be modeled using the generalized method of cells (GMC) micromechanics theory, and local failure in the matrix constituent subcells will be modeled using the crack band theory. The subcell characteristic lengths used in the crack band calculations will be mapped to the macroscale finite element in order to regularize the local energy in a manner consistent with the global length scale. Examples will be provided with and without the regularization, and they will be compared to a baseline case where the size and shape of the element and RUC are coincident (ensuring energy is preserved across the scales).
A thermoviscoplastic computational method for hypersonic structures is presented. The method employs unified viscoplastic constitutive model implemented in a finite element approach for quasi-static thermal-structural analysis. Applications of the approach to convectively cooled hypersonic structures illustrate the effectiveness of the approach and provide insight into the transient inelastic structural behavior at elevated temperatures.
An alternative time-dependent wave-packet method for treating three-dimensional gas phase reactive atom-diatom collisions is presented. The method employs a nonreactive body-frame wave packet propagation procedure, made possible by judicious use of absorbing optical potentials, a novel scheme for interpolating the wave function from coordinates in one arrangement to those in another and the fact that the time-dependent Schroedinger equation is an initial-value problem. The last feature makes possible a computationally viable and accurate procedure for changing from one arrangement's coordinates to another. In addition, the method allows the determination of S-matrix elements over a wide range of energies from a single wave-packet propagation. The method is illustrated by carrying out detailed calculations of inelastic and reactive scattering in the H + H2 system using the Liu-Siegbahn-Truhlar-Horowitz potential surface.
General simulation equations are derived for the rigid body motion of slung-load systems. This work is motivated by an interest in trajectory control for slung loads carried by two or more helicopters. An approximation of these systems consists of several rigid bodies connected by straight-line cables or links. The suspension can be assumed elastic or inelastic. Equations for the general system are obtained from the Newton-Euler rigid-body equations with the introduction of generalized velocity coordinates. Three forms are obtained: two generalize previous case-specific results for single-helicopter systems with elastic and inelastic suspensions, respectively; and the third is a new formulation for inelastic suspensions. The latter is derived from the elastic suspension equations by choosing the generalized coordinates so that motion induced by cable stretching is separated from motion with invariant cable lengths, and by then nulling the stretching coordinates to get a relation for the suspension forces. The result is computationally more efficient than the conventional formulation, is readily integrated with the elastic suspension formulation, and is easily applied to the complex dual-lift and multilift systems. Results are given for two-helicopter systems; three configurations are included and these can be integrated in a single simulation. Equations are also given for some single-helicopter systems, for comparison with the previous literature, and for a multilift system. Equations for degenerate-body approximations (point masses, rigid rods) are also formulated and results are given for dual-lift and multilift systems. Finally, linearlized equations of motion are given for general slung-load systems are presented along with results for the two-helicopter system with a spreader bar.
The results presented here are part of an ongoing research program, to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to high strain rate impact loads. A micromechanics approach is employed in this work, in which state variable constitutive equations originally developed for metals have been modified to model the deformation of the polymer matrix, and a strength of materials based micromechanics method is used to predict the effective response of the composite. In the analysis of the inelastic deformation of the polymer matrix, the definitions of the effective stress and effective inelastic strain have been modified in order to account for the effect of hydrostatic stresses, which are significant in polymers. Two representative polymers, a toughened epoxy and a brittle epoxy, are characterized through the use of data from tensile and shear tests across a variety of strain rates. Results computed by using the developed constitutive equations correlate well with data generated via experiments. The procedure used to incorporate the constitutive equations within a micromechanics method is presented, and sample calculations of the deformation response of a composite for various fiber orientations and strain rates are discussed.
Activities reported include: (1) testing of the computer program used to obtain transport properties for the Hulburt-Hirschfelder potential; (2) calculation of transport properties for the C2-C interaction; (3) preliminary calculations for the C2-C2 interaction; (4) calculation of transport properties for the C2H-He interaction; (5) consideration of the effect of inelastic collisions on the transport properties; and (6) the use of the Hulburt-Hirschfelder potential to model ion-atom interactions.