Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational inelasticity”

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 109 records · Page 6

3-D inelastic analysis methods for hot section components. Volume 2: Advanced special functions models

This Annual Status Report presents the results of work performed during the third year of the 3-D Inelastic Analysis Methods for Hot Sections Components program (NASA Contract NAS3-23697). The objective of the program is to produce a series of computer codes that permit more accurate and efficient three-dimensional analyses of selected hot section components, i.e., combustor liners, turbine blades, and turbine vanes. The computer codes embody a progression of mathematical models and are streamlined to take advantage of geometrical features, loading conditions, and forms of material response that distinguish each group of selected components.

Wilson, R. B.↗

Probabilistic Evaluation of Advanced Ceramic Matrix Composite Structures

The objective of this report is to summarize the deterministic and probabilistic structural evaluation results of two structures made with advanced ceramic composites (CMC): internally pressurized tube and uniformly loaded flange. The deterministic structural evaluation includes stress, displacement, and buckling analyses. It is carried out using the finite element code MHOST, developed for the 3-D inelastic analysis of structures that are made with advanced materials. The probabilistic evaluation is performed using the integrated probabilistic assessment of composite structures computer code IPACS. The affects of uncertainties in primitive variables related to the material, fabrication process, and loadings on the material property and structural response behavior are quantified. The primitive variables considered are: thermo-mechanical properties of fiber and matrix, fiber and void volume ratios, use temperature, and pressure. The probabilistic structural analysis and probabilistic strength results are used by IPACS to perform reliability and risk evaluation of the two structures. The results will show that the sensitivity information obtained for the two composite structures from the computational simulation can be used to alter the design process to meet desired service requirements. In addition to detailed probabilistic analysis of the two structures, the following were performed specifically on the CMC tube: (1) predicted the failure load and the buckling load, (2) performed coupled non-deterministic multi-disciplinary structural analysis, and (3) demonstrated that probabilistic sensitivities can be used to select a reduced set of design variables for optimization.

Abumeri, Galib H.↗

Three-dimensional inelastic approximate analysis code (MOMM)

The Mechanics of Materials Model (MOMM) is one of a series of new stand-alone three dimensional nonlinear structural analysis codes. Incorporation of a general purpose finite element computer code into the hot section design process was severely limited by the high costs involved. MOMM is a stiffness method finite element code that uses an internally generated network of beams to characterize hot section component behavior. The method was proposed as a fast, easy to use, computationally efficient tool for approximate analyses. MOMM incorporates a wide variety of analysis capabilities, material models, and load type specifiers instrumental for the analysis of hot section components.

Meister, Jeffrey P.↗

On 3-D inelastic analysis methods for hot section components. Volume 1: Special finite element models

This Annual Status Report presents the results of work performed during the third year of the 3-D Inelastic Analysis Methods for Hot Section Components program (NASA Contract NAS3-23697). The objective of the program is to produce a series of new computer codes that permit more accurate and efficient three-dimensional analysis of selected hot section components, i.e., combustor liners, turbine blades, and turbine vanes. The computer codes embody a progression of mathematical models and are streamlined to take advantage of geometrical features, loading conditions, and forms of material response that distinguish each group of selected components. This report is presented in two volumes. Volume 1 describes effort performed under Task 4B, Special Finite Element Special Function Models, while Volume 2 concentrates on Task 4C, Advanced Special Functions Models.

Nakazawa, S.↗

On 3-D inelastic analysis methods for hot section components. Volume 1: Special finite element models

This annual status report presents the results of work performed during the fourth year of the 3-D Inelastic Analysis Methods for Hot Section Components program (NASA Contract NAS3-23697). The objective of the program is to produce a series of new computer codes permitting more accurate and efficient 3-D analysis of selected hot section components, i.e., combustor liners, turbine blades and turbine vanes. The computer codes embody a progression of math models and are streamlined to take advantage of geometrical features, loading conditions, and forms of material response that distinguish each group of selected components. Volume 1 of this report discusses the special finite element models developed during the fourth year of the contract.

Nakazawa, S.↗

Body frame close coupling wave packet approach to gas phase atom-rigid rotor inelastic collisions

The close coupling wave packet (CCWP) method is formulated in a body-fixed representation for atom-rigid rotor inelastic scattering. For J greater than j-max (where J is the total angular momentum and j is the rotational quantum number), the computational cost of propagating the coupled channel wave packets in the body frame is shown to scale approximately as N exp 3/2, where N is the total number of channels. For large numbers of channels, this will be much more efficient than the space frame CCWP method previously developed which scales approximately as N-squared under the same conditions.

Sun, Y.↗

Perspectives on hypersonic viscous and nonequilibrium flow research

An attempt is made to reflect on current focuses in certain areas of hypersonic flow research by examining recent works and their issues. Aspects of viscous interaction, flow instability, and nonequilibrium aerothermodynamics pertaining to theoretical interest are focused upon. The field is a diverse one, and many exciting works may have either escaped the writer's notice or been abandoned for the sake of space. Students of hypersonic viscous flow must face the transition problems towards the two opposite ends of the Reynolds or Knudsen number range, which represents two regimes where unresolved fluid/gas dynamic problems abound. Central to the hypersonic flow studies is high-temperature physical gas dynamics; here, a number of issues on modelling the intermolecular potentials and inelastic collisions remain the obstacles to quantitative predictions. Research in combustion and scramjet propulsion will certainly be benefitted by advances in turbulent mixing and new computational fluid dynamics (CFD) strategies on multi-scaled complex reactions. Even for the sake of theoretical development, the lack of pertinent experimental data in the right energy and density ranges is believed to be among the major obstacles to progress in aerothermodynamic research for hypersonic flight. To enable laboratory simulation of nonequilibrium effects anticipated for transatmospheric flight, facilities capable of generating high enthalpy flow at density levels higher than in existing laboratories are needed (Hornung 1988). A new free-piston shock tunnel capable of realizing a test-section stagnation temperature of 10(exp 5) at Reynolds number 50 x 10(exp 6)/cm is being completed and preliminary tests has begun (H. Hornung et al. 1992). Another laboratory study worthy of note as well as theoretical support is the nonequilibrium flow experiment of iodine vapor which has low activation energies for vibrational excitation and dissociation, and can be studied in a laboratory with modest resources (Pham-Van-Diep et al. 1992).

Cheng, H. K.↗

Implementation of Fiber Substructuring Into Strain Rate Dependent Micromechanics Analysis of Polymer Matrix Composites

A research program is in progress to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to impact loads. Previously, strain rate dependent inelastic constitutive equations developed to model the polymer matrix were incorporated into a mechanics of materials based micromechanics method. In the current work, the micromechanics method is revised such that the composite unit cell is divided into a number of slices. Micromechanics equations are then developed for each slice, with laminate theory applied to determine the elastic properties, effective stresses and effective inelastic strains for the unit cell. Verification studies are conducted using two representative polymer matrix composites with a nonlinear, strain rate dependent deformation response. The computed results compare well to experimentally obtained values.

Goldberg, Robert K.↗

High-Fidelity Micromechanics Model Enhanced for Multiphase Particulate Materials

This 3-year effort involves the development of a comprehensive micromechanics model and a related computer code, capable of accurately estimating both the average response and the local stress and strain fields in the individual phases, assuming both elastic and inelastic behavior. During the first year (fiscal year 2001) of the investigation, a version of the model called the High-Fidelity Generalized Method of Cells (HFGMC) was successfully completed for the thermo-inelastic response of continuously reinforced multiphased materials with arbitrary periodic microstructures (refs. 1 and 2). The model s excellent predictive capability for both the macroscopic response and the microlevel stress and strain fields was demonstrated through comparison with exact analytical and finite element solutions. This year, HFGMC was further extended in two technologically significant ways. The first enhancement entailed the incorporation of fiber/matrix debonding capability into the two-dimensional version of HFGMC for modeling the response of unidirectionally reinforced composites such as titanium matrix composites, which exhibit poor fiber/matrix bond. Comparison with experimental data validated the model s predictive capability. The second enhancement entailed further generalization of HFGMC to three dimensions to enable modeling the response of particulate-reinforced (discontinuous) composites in the elastic material behavior domain. Next year, the three-dimensional version will be generalized to encompass inelastic effects due to plasticity, viscoplasticity, and damage, as well as coupled electromagnetothermomechanical (including piezoelectric) effects.

Pindera, Marek-Jerzy↗

A simple theoretical model for calculating and parameterizing the ionospheric photoelectron flux

A method for calculating the ionospheric photoelectron flux is developed which uses the concept of average electron energy loss to simplify the calculation of the degraded electron spectrum. This method requires only a knowledge of the total inelastic electron impact cross sections and can be used for the calculation of all secondary ion and excited state production rates. The simple calculation reduces the computing time by a factor of 10 and considerably reduces storage requirements. It is found that the ionospheric photoelectron flux in the local equilibrium region is directly proportional to the attenuated solar EUV flux, which is a function only of the total neutral column density, and is independent of the neutral density composition. Thus, it is shown that electron impact cross section can be chosen so that this method may also be used to parameterize the measured ionospheric photoelectron fluxes.

Richards, P. G.↗

GAMBIT program

GAMBIT computer program provides multigroup flux and current averaged neutron cross sections for input to transport programs. The fast neutron compilation includes P/0/ and P/1/ elastic scattering transfer matrices, inelastic scattering matrices, n,2n cross section matrices, absorption and fission cross sections for all isotopes.

Collier, G.↗

Analytical and experimental studies of graphite-epoxy and boron-epoxy angle ply laminates in compression

The applicability and adequacy of several computer techniques in predicting satisfactorily the nonlinear/inelastic response of angle ply laminates were evaluated. The analytical predictions were correlated with the results of a test program on the inelastic response under axial compression of a large variety of graphite-epoxy and boron-epoxy angle ply laminates. These comparison studies indicate that neither of the abovementioned analyses can satisfactorily predict either the mode of response or the ultimate stress value corresponding to a particular angle ply laminate configuration. Consequently, also the simple failure mechanisms assumed in the analytical models were not verified.

Weller, T.↗

The MHOST finite element program: 3-D inelastic analysis methods for hot section components. Volume 1: Theoretical manual

Formulations and algorithms implemented in the MHOST finite element program are discussed. The code uses a novel concept of the mixed iterative solution technique for the efficient 3-D computations of turbine engine hot section components. The general framework of variational formulation and solution algorithms are discussed which were derived from the mixed three field Hu-Washizu principle. This formulation enables the use of nodal interpolation for coordinates, displacements, strains, and stresses. Algorithmic description of the mixed iterative method includes variations for the quasi static, transient dynamic and buckling analyses. The global-local analysis procedure referred to as the subelement refinement is developed in the framework of the mixed iterative solution, of which the detail is presented. The numerically integrated isoparametric elements implemented in the framework is discussed. Methods to filter certain parts of strain and project the element discontinuous quantities to the nodes are developed for a family of linear elements. Integration algorithms are described for linear and nonlinear equations included in MHOST program.

Nakazawa, Shohei↗

Slow and fast motion of cracks in inelastic solids. Part 1: Slow growth of cracks in a rate sensitive tresca solid. Part 2: Dynamic crack represented by the Dugdale model

An extension is proposed of the classical theory of fracture to viscoelastic and elastic-plastic materials in which the plasticity effects are confined to a narrow band encompassing the crack front. It is suggested that the Griffith-Irwin criterion of fracture, which requires that the energy release rate computed for a given boundary value problem equals the critical threshold, ought to be replaced by a differential equation governing the slow growth of a crack prior to the onset of rapid propagation. A new term which enters the equation of motion in the dissipative media is proportional to the energy lost within the end sections of the crack, and thus reflects the extent of inelastic behavior of a solid. A concept of apparent surface energy is introduced to account for the geometry dependent and the rate dependent phenomena which influence toughness of an inelastic solid. Three hypotheses regarding the condition for fracture in the subcritical range of load are compared. These are: (1) constant fracture energy (Cherepanov), (2) constant opening displacement at instability (Morozov) and (3) final stretch criterion (Wnuk).

Wnuk, M. P.↗

Computer simulation of the Farley-Buneman instability and anomalous electron heating in the auroral ionosphere

The nonlinear saturation of the Farley-Buneman instability in a collisional plasma is studied by a 2 1/2 dimensional electrostatic particle simulation which includes inelastic and elastic collisions of electrons and elastic collision of ions with neutrals. In the simulation, a uniform convection electric field is applied externally so that the relative velocity between the electrons and ions is greater than the ion sound speed and destabilizes the instability. A nonlinear frequency shift from higher to lower frequencies and diffusion of the wave spectrum in two dimensional wave number space are found. It is found that the dominant mechanism for electron heating is due to an enhanced effective electron collision frequency and hence enhanced resistive heating as suggested by Primdahl (1986) and not due to the heating of electrons by the electric field of the waves parallel to the magnetic field. For the ionospheric conditions discussed by Schlegel and St.-Maurice (1981)an anomalous heating rate of about 4 x 10 to the 7th W/cu cm is found.

Machida, S.↗

Comment on the accuracy of Rabitz' effective potential approximation for rotational excitation by collisions

Cross sections for rotational excitation of HCN by low-energy collisions with He have been computed with the effective-potential approximation of Rabitz (1972) and compared with accurate quantum close-coupling results. Elastic cross sections are found to agree to about 20%; inelastic cross sections agree in general magnitude, but not in detailed values for specific quantum transitions.

Green, S.↗

Postseismic rebound due to creep of the lower lithosphere and asthenosphere

Postseismic surface deformations are attributed to the inelastic flow of the subcrustal regions of the Earth following an earthquake. A multilayer representation of the Earth's rheological properties is used in conjunction with a finite element computational scheme to calculate time dependent displacements and strains subsequent to a strike slip earthquake. The deviatoric stress strain relations for the uppermost layer is assumed elastic. Lower layers are assumed to be, in order of increasing depth, a standard, linear, three element, viscoelastic solid; a linear viscoelastic fluid; and another elastic solid. Physically these layers correspond to the upper lithosphere, lower lithosphere, asthenosphere, and lower mantle, respectively. Elastic dilatational properties are assumed throughout. Appreciable postseismic displacements, possibly approaching meters, for large earthquakes, arise from viscoelastic relaxation following sudden coseismic slip. Furthermore, compared to the simpler case of an elastic lithosphere over a viscoelastic asthenosphere and the near fault postseismic shear, strain is increased, by a factor of two or more in some cases, by the presence of a viscoelastic lower lithosphere. Also, the duration of postseismic straining is increased by the relatively slow relaxation of this layer.

Cohen, S. C.↗

Simplified cyclic structural analyses of SSME turbine blades

Anisotropic high-temperature alloys are used to meet the safety and durability requirements of turbine blades for high-pressure turbopumps in reusable space propulsion systems. The applicability to anisotropic components of a simplified inelastic structural analysis procedure developed at the NASA Lewis Research Center is assessed. The procedure uses as input the history of the total strain at the critical crack initiation location computed from elastic finite-element analyses. Cyclic heat transfer and structural analyses are performed for the first stage high-pressure fuel turbopump blade of the space shuttle main engine. The blade alloy is directionally solidified MAR-M 246 (nickel base). The analyses are based on a typical test stand engine cycle. Stress-strain histories for the airfoil critical location are computed using both the MARC nonlinear finite-element computer code and the simplified procedure. Additional cases are analyzed in which the material yield strength is arbitrarily reduced to increase the plastic strains and, therefore, the severity of the problem. Good agreement is shown between the predicted stress-strain solutions from the two methods. The simplified analysis uses about 0.02 percent (5 percent with the required elastic finite-element analyses) of the CPU time used by the nonlinear finite element analysis.

Kaufman, A.↗