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.
Equations for thin elastic conical shells and digital program for analysis of stress and deformation on fixed edge segmental conical shells - solution by finite difference technique
The component element method was used to develop a transient dynamic analysis computer program which is essentially based on modal synthesis combined with a central, finite difference, numerical integration scheme. The methodology leads to a modular or building-block technique that is amenable to computer programming. To verify the analytical method, turbine engine transient response analysis (TETRA), was applied to two blade-out test vehicles that had been previously instrumented and tested. Comparison of the time dependent test data with those predicted by TETRA led to recommendations for refinement or extension of the analytical method to improve its accuracy and overcome its shortcomings. The development of working equations, their discretization, numerical solution scheme, the modular concept of engine modelling, the program logical structure and some illustrated results are discussed. The blade-loss test vehicles (rig full engine), the type of measured data, and the engine structural model are described.
Flow/damage surfaces can be defined in terms of stress, inelastic strain rate, and internal variables using a thermodynamics framework. A macroscale definition relevant to thermodynamics and usable in an experimental program is employed to map out surfaces of constant inelastic power in various stress planes. The inelastic flow of a model silicon carbide/ titanium composite system having rectangular, hexagonal, and square diagonal fiber packing arrays subjected to biaxial stresses is quantified by flow/damage surfaces that are determined numerically from micromechanics, using both finite element analysis and the generalized method of cells. Residual stresses from processing are explicitly included and damage in the form of fiber-matrix debonding under transverse tensile and/or shear loading is represented by a simple interface model. The influence of microstructural architecture is largest whenever fiber-matrix debonding is not an issue; for example in the presence of transverse compressive stresses. Additionally, as the fiber volume fraction increases, so does the effect of microstructural architecture. With regard to the micromechanics analysis, the overall inelastic flow predicted by the generalized method of cells is in excellent agreement with that predicted using a large number of displacement-based finite elements.
Flow/damage surfaces can be defined in terms of stress, inelastic strain rate, and internal variables using a thermodynamics framework. A macroscale definition relevant to thermodynamics and usable in an experimental program is employed to map out surfaces of constant inelastic power in various stress planes. The inelastic flow of a model silicon carbide/ titanium composite system having rectangular, hexagonal, and square diagonal fiber packing, arrays subjected to biaxial stresses is quantified by flow/damage surfaces that are determined numerically from micromechanics. using both finite element analysis and the generalized method of cells. Residual stresses from processing are explicitly included and damage in the form of fiber-matrix debonding under transverse tensile and/or shear loading is represented by a simple interface model. The influence of microstructural architecture is largest whenever fiber-matrix debonding is not an issue, for example in the presence of transverse compressive stresses. Additionally, as the fiber volume fraction increases, so does the effect of microstructural architecture. With regard to the micromechanics analysis, the overall inelastic flow predicted by the generalized method of cells is in excellent agreement with that predicted using a large number of displacement-based finite elements.
The purpose of the Probabilistic Structural Analysis Method (PSAM) project is to develop structural analysis capabilities for the design analysis of advanced space propulsion system hardware. The boundary element method (BEM) is used as the basis of the Probabilistic Advanced Analysis Methods (PADAM) which is discussed. The probabilistic BEM code (PBEM) is used to obtain the structural response and sensitivity results to a set of random variables. As such, PBEM performs analogous to other structural analysis codes such as finite elements in the PSAM system. For linear problems, unlike the finite element method (FEM), the BEM governing equations are written at the boundary of the body only, thus, the method eliminates the need to model the volume of the body. However, for general body force problems, a direct condensation of the governing equations to the boundary of the body is not possible and therefore volume modeling is generally required.
A great deal of research has been performed characterizing the in-plane fiber-dominated properties, under both static and fatigue loading, of advanced composite materials. To the author's knowledge, no study has been performed to date investigating fatigue characteristics in the transverse direction. This information is important in the design of bonded composite airframe structure where repeated, cyclic out-of-plane bending may occur. Recent tests characterizing skin/stringer debond failures in reinforced composite panels where the dominant loading in the skin is flexure along the edge of the frame indicate failure initiated either in the skin or else the flange, near the flange tip. When failure initiated in the skin, transverse matrix cracks formed in the surface skin ply closest to the flange and either initiated delaminations or created matrix cracks in the next lower ply, which in turn initiated delaminations. When failure initiated in the flanges, transverse cracks formed in the flange angle ply closest to the skin and initiated delamination. In no configuration did failure propagate through the adhesive bond layer. For the examined skin/flange configurations, the maximum transverse tension stress at failure correlates very well with the transverse tension strength of the composites. Transverse tension strength (static) data of graphite epoxy composites have been shown to vary with the volume of material stressed. As the volume of material stressed increased, the strength decreased. A volumetric scaling law based on Weibull statistics can be used to predict the transverse strength measurements. The volume dependence reflects the presence of inherent flaws in the microstructure of the lamina. A similar approach may be taken to determine a volume scale effect on the transverse tension fatigue behavior of graphite/epoxy composites. The objective of this work is to generate transverse tension strength and fatigue S-N characteristics for composite materials using 3-point flexure tests of 90 deg graphite/epoxy specimens. Investigations will include the volume scale effect as well as frequency and span-to-thickness ratio effects. Prior to the start of the experimental study, an analytical study using finite element modeling will be performed to investigate the span-to-thickness effect. The ratio of transverse flexure stress to shear stress will be monitored and its values predicted by the FEM analysis compared with the value obtained using a 'strength of materials' based approach.
The cluster expansion technique is used in combination with first-principles calculations of the total energy of ordered compounds to study the energetics of the disordered phase for a number of alloy systems. The effect of short-range order, as seen in the energy differences between an alloy with the configuration corresponding to finite temperature and the perfectly random one is studied. The role of relaxation due to large size mismatch is discussed in terms of an effective cluster volume approximation. Very good agreement with measured energies of formation is obtained when short-range order and relaxation effects are taken into account.
A description of a computer program entitled VORTEX that may be used to determine the aerodynamic performance of horizontal axis wind turbines is given. The computer code implements a vortex method from finite span wind theory and determines the induced velocity at the rotor disk by integrating the Biot-Savart law. It is assumed that the trailing helical vortex filaments form a wake of constant diameter (the rigid wake assumption) and travel downstream at the free stream velocity. The program can handle rotors having any number of blades which may be arbitrarily shaped and twisted. Many numerical details associated with the program are presented. A complete listing of the program is provided and all program variables are defined. An example problem illustrating input and output characteristics is solved.
A transient, semi-coupled, multi-dimensional thermal and flow analysis methodology was developed to predict the thermal/gas dynamic conditions in the field joint region of the Redesigned Solid Rocket Motor (RSRM). Transient temperature response, pressure history, and combustion gas flow rate (within the field joint region), were of principle interest, in the course of this study. The thermal environment in the field joint was modeled using SINDA, a finite difference based thermal network analyzer. The combustion gas flow boundary condition was generated using the FLAP code; this code performs a transient, lumped-parameter, control volume analysis to solve the mass, momentum, and energy conservation equations. The FLAP computer code was modified to account for erosion of the NBR insulation material, following ignition. An independent grid sensitivity study was conducted to determine an appropriate grid distribution near the wall. The predicted results, obtained using an optimum grid distribution and computer generated flow boundary condition, were compared with subscale test data.
Described are the theoretical development and computer implementation of reliable and efficient methods for the analysis of coupled mechanical problems that involve the interaction of mechanical, thermal, phase-change and electromagnetic subproblems. The focus application has been the modeling of superconductivity and associated quantum-state phase-change phenomena. In support of this objective the work has addressed the following issues: (1) development of variational principles for finite elements, (2) finite element modeling of the electromagnetic problem, (3) coupling of thermal and mechanical effects, and (4) computer implementation and solution of the superconductivity transition problem. The main accomplishments have been: (1) the development of the theory of parametrized and gauged variational principles, (2) the application of those principled to the construction of electromagnetic, thermal and mechanical finite elements, and (3) the coupling of electromagnetic finite elements with thermal and superconducting effects, and (4) the first detailed finite element simulations of bulk superconductors, in particular the Meissner effect and the nature of the normal conducting boundary layer. The theoretical development is described in two volumes. This volume, Volume 1, describes mostly formulations for specific problems. Volume 2 describes generalization of those formulations.
A new plasma accelerator concept that employs electrodeless plasma preionization and pulsed inductive acceleration is presented. Preionization is achieved through an electron cyclotron resonance discharge that produces a weakly-ionized plasma at the face of a conical theta pinch-shaped inductive coil. The presence of the preionized plasma allows for current sheet formation at lower discharge voltages than those found in other pulsed inductive accelerators. The location of an electron cyclotron resonance discharge can be controlled through the design of the applied magnetic field in the thruster. A finite-element model of the magnetic field was used as a design tool, allowing for the implementation of an arrangement of permanent magnets that yields a small volume of preionized propellant at the coil face. This allows for current sheet formation at the face of the inductive coil, minimizing the initial inductance of the pulse circuit and maximizing the potential efficiency of the new accelerator.
An artificial convolutional neural network was created to efficiently mimic a micromechanics model, the High Fidelity Generalized Method of Cells, for use in multi-scale structural finite element analysis. The network was found to quickly and accurately replicate the stiffness predicted by the micromechanics model using a 2D image of an idealized representative volume element of a fiber/matrix microstructure. The long-term goal of this work is to efficiently apply multi-scale methods for predicting the damage progression of a composite structure.
A venting equation commonly used to describe the transient evolution of pressure within a volume containing an outgassing or offgassing source and a restrictive vent conductance under conditions of molecular flow has been solved analytically. Solutions are found for sources of finite thickness characterized by classical diffusion-limited behavior (proportional to inverse square root of time), as well as responses for thick material sources often observed in testing that are characterized by a more general form of power-law decay, up to inverse time behavior associated with surface desorption. Solutions involve evaluating integrals where both numerators and denominators of the integrands diverge with time, making wide-ranging transient solutions difficult to directly compute numerically. Usually, one can avoid evaluating these integrals by assuming quasistatic conditions at long time scales. A novel approach is used in this work to analytically produce solutions by generating bespoke mathematical functions, some of which solve integrals that have apparently had no previous analytical solution.
"A venting equation commonly used to describe the transient evolution of pressure within a volume containing an outgassing or offgassing source and a restrictive vent conductance under conditions of molecular flow has been solved analytically. Solutions are found for sources of finite thickness characterized by classical diffusion-limited behavior (proportional to inverse square root of time), as well as responses for thick material sources often observed in testing that are characterized by a more general form of power-law decay, up to inverse time behavior associated with surface desorption. Solutions involve evaluating integrals where both numerators and denominators of the integrands diverge with time, making wide-ranging transient solutions difficult to directly compute numerically. Usually, one can avoid evaluating these integrals by assuming quasistatic conditions at long time scales. A novel approach is used in this work to analytically produce solutions by generating bespoke mathematical functions, some of which solve integrals that have apparently had no previous analytical solution."
The method used to formulate the fluid/structure equations is described in the updates to the NASTRAN Theoretical Manual. The basis for defining the fluid are three dimensional finite elements connected to fluid grid points defining the Eulerian pressure at a point fixed in space. The use of a pressure single degree of freedom at each point rather than three displacements allows a finer mesh of elements with a reasonable matrix order. In the formulation of the fluid/structure system the interior fluid degrees of freedom are transformed and removed from the solution matrices. The eigenvalues of the combination are extracted from small, fully dense, symmetric mass and stiffness matrices, efficiently processed with the Givens method. The solution matrices are defined by the free surface displacements and the reduced structure coordinates.
The analysis performed on the High Pressure Oxidizer Turbopump (HPOTP) preburner pump bearing assembly located on the Space Shuttle Main Engine (SSME) is summarized. An ANSYS finite element model for the inlet assembly was built and executed. Thermal and static analyses were performed.
Parametric studies were performed to assess the effects of various parameters on the free-vibration behavior (natural frequencies) of (+/- theta)(sub 2) angle-ply, fiber composite, thin shell structures in a hygrothermal environment. Knowledge of the natural frequencies of structures is important in considering their response to various kinds of excitation, especially when structures and force systems are complex and when excitations are not periodic. The three dimensional, finite element structural analysis computer code CSTEM was used in the Cray YMP computer environment. The fiber composite shell was assumed to be cylindrical and made from T300 graphite fibers embedded in an intermediate-modulus, high-strength matrix. The following parameters were investigated: the length and the laminate thickness of the shell, the fiber orientation, the fiber volume fraction, the temperature profile through the thickness of the laminate, and laminates with different ply thicknesses. The results indicate that the fiber orientation and the length of the laminated shell had significant effects on the natural frequencies. The fiber volume fraction, the laminate thickness, and the temperature profile through the shell thickness had weak effects on the natural frequencies. Finally, the laminates with different ply thicknesses had an insignificant influence on the behavior of the vibrated laminated shell. Also, a single through-the-thickness, eight-node, three dimensional composite finite element analysis appears to be sufficient for investigating the free-vibration behavior of thin, composite, angle-ply shell structures.
A three dimensional elastoplastic finite element micromechanical model was developed to study the state of stress around a broken fibers in a unidirectional composite. A boron/aluminum composite consisting of 50 percent by volume of fibers in a square array and subjected to an axial loading is taken as a specific example. This loading in the fiber direction is applied in small increments, by prescribing increments of boundary displacement, until the first failure occurs. The effect of reduced material properties of the aluminum matrix material at elevated temperature is also studied. The results are presented in the form of stress contours and stress-strain plots.