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 application of NASTRAN/COSMIC in predicting the transient motion of ship structures to underwater, non-contact explosions is discussed. Examples illustrate the finite element models, mathematical formulations of loading functions and, where available, comparisons between analytical and experimental results.
A coupled finite element/boundary element capability is described for calculating the sound pressure field scattered by an arbitrary submerged 3-D elastic structure. Structural and fluid impedances are calculated with no approximation other than discretization. The surface fluid pressures and normal velocities are first calculated by coupling a NASTRAN finite element model of the structure with a discretized form of the Helmholtz surface integral equation for the exterior field. Far field pressures are then evaluated from the surface solution using the Helmholtz exterior integral equation. The overall approach is illustrated and validated using a known analytic solution for scattering from submerged spherical shells.
The linear-quadratic Gaussian/loop-transfer-recovery method is used to synthesize a fine-pointing control system for a large space antenna. A finite-element model for the 122-m hoop/column antenna is employed, and a compensator, utilizing attitudes sensors and torque actuators, is designed which achieves pointing performance while maintaining stability robustness to unmodeled dynamics. Inclusion of the rigid-body modes plus the first three elastic modes is found to be necessary to achieve a 0.1-rad/s bandwidth. Results are obtained by employing a modification of the standard robustness recovery procedure, which reduces the conservative nature of the design methodology. Performance degradation is encountered due to the presence of unavoidable invariant zeros within the design bandwidth.
Solution techniques for handling large scale engineering optimization problems are reviewed. Potentials for practical applications as well as their limited capabilities are discussed. A new solution algorithm for design sensitivity is proposed. The algorithm is based upon the multilevel substructuring concept to be coupled with the adjoint method of sensitivity analysis. There are no approximations involved in the present algorithm except the usual approximations introduced due to the discretization of the finite element model. Results from the six- and thirty-bar planar truss problems show that the proposed multilevel scheme for sensitivity analysis is more effective (in terms of computer incore memory and the total CPU time) than a conventional (one level) scheme even on small problems. The new algorithm is expected to perform better for larger problems and its applications on the new generation of computer hardwares with 'parallel processing' capability is very promising.
An extension of the program STARS (a general-purpose structural analysis program) has been developed; this extension implements a complete aeroservoelastic analysis capability. Previous capabilities included finite-element modeling as well as statics, buckling, vibration, dynamic response, and flutter analyses. This paper presents a description and the formulation of STARS in its current state along with example dynamic, aeroelastic, and aeroservoelastic analyses pertaining to the X-29A aircraft. These examples include vibration analysis results as well as flutter analysis results obtained by the conventional k method and the velocity root-contour solution. Finally, selected open- and closed-loop aeroservoelastic analysis results based on a hybrid formulation are compared to illustrate, using the calculated frequency responses, the interactions of structures, aerodynamics, and flight controls.
This paper describes a new methodology for the formulation of the aeroelastic stability and response problem for helicopter rotor blades. The mathematical expressions for the aerodynamic loads need not be explicit functions of the blade displacement quantities. This methodology is combined with a finite element model of the blade, and a quasilinearization solution technique. The resulting computer program is used to study the behavior of blades with noncoincident elastic axis, aerodynamic centers, and centers of mass.
The structural dynamic characteristics of the International Space Station (ISS), the interim reference configuration established for NASA's Space Station developmental program, are discussed, and a finite element model is described. Modes and frequencies of the station below 2.0 Hz are derived, and the dynamic response of the station is simulated for an external impulse load corresponding to a failed shuttle-docking maneuver. A three-axis attitude control system regulates the ISS orientation, with control moment gyros responding to attitude and attitude rate signals. No instabilities were found in the attitude control system.
Meniscus-defined techniques for the growth of thin silicon sheets fall into two categories: vertical and horizontal growth. The interactions of the temperature field and the crystal shape are analyzed for both methods using two-dimensional finite-element models which include heat transfer and capillarity. Heat transfer in vertical growth systems is dominated by conduction in the melt and the crystal, with almost flat melt/crystal interfaces that are perpendicular to the direction of growth. The high axial temperature gradients characteristic of vertical growth lead to high thermal stresses. The maximum growth rate is also limited by capillarity which can restrict the conduction of heat from the melt into the crystal. In horizontal growth the melt/crystal interface stretches across the surface of the melt pool many times the crystal thickness, and low growth rates are achievable with careful temperature control. With a moderate axial temperature gradient in the sheet a substantial portion of the latent heat conducts along the sheet and the surface of the melt pool becomes supercooled, leading to dendritic growth. The thermal supercooling is surpressed by lowering the axial gradient in the crystal; this configuration is the most desirable for the growth of high quality crystals. An expression derived from scaling analysis relating the growth rate and the crucible temperature is shown to be reliable for horizontal growth.
Thermal stresses developed during the growth of silicon crystal ribbon have been shown to be negligible, thus eliminating residual stresses and dislocations, if the temperature profile satisfies a second-order partial differential equation inside the ribbon. This has been numerically verified through a finite element model, an outline of which is presented here. This model shows that, for homogeneous isotropic material with temperature independent thermal expansion coefficients, thermal stresses will vanish if the temperature profile satisfies the Laplacian. A comparison of stresses due to uniform and nonuniform temperature gradients in the plane of the ribbon is also presented. The strategies employed to control the round-off error and to validate the computer model are discussed.
The analysis code BUMPER executes a numerical solution to the problem of calculating the probability of no penetration (PNP) of a spacecraft subject to man-made orbital debris or meteoroid impact. The codes were developed on a DEC VAX 11/780 computer that uses the Virtual Memory System (VMS) operating system, which is written in FORTRAN 77 with no VAX extensions. To help illustrate the steps involved, a single sample analysis is performed. The example used is the space station reference configuration. The finite element model (FEM) of this configuration is relatively complex but demonstrates many BUMPER features. The computer tools and guidelines are described for constructing a FEM for the space station under consideration. The methods used to analyze the sensitivity of PNP to variations in design, are described. Ways are suggested for developing contour plots of the sensitivity study data. Additional BUMPER analysis examples are provided, including FEMs, command inputs, and data outputs. The mathematical theory used as the basis for the code is described, and illustrates the data flow within the analysis.
Temperature, thermal stresses, and residual creep stresses were studied by comparing laboratory values measured on a built-up titanium structure with values calculated from finite-element models. Several such models were used to examine the relationship between computational thermal stresses and thermal stresses measured on a built-up structure. Element suitability, element density, and computational temperature discrepancies were studied to determine their impact on measured and calculated thermal stress. The optimum number of elements is established from a balance between element density and suitable safety margins, such that the answer is acceptably safe yet is economical from a computational viewpoint. It is noted that situations exist where relatively small excursions of calculated temperatures from measured values result in far more than proportional increases in thermal stress values. Measured residual stresses due to creep significantly exceeded the values computed by the piecewise linear elastic strain analogy approach. The most important element in the computation is the correct definition of the creep law. Computational methodology advances in predicting residual stresses due to creep require significantly more viscoelastic material characterization.
The ability to accurately characterize propellant in a finite element model is a concern of engineers tasked with studying the dynamic response of the Space Shuttle Solid Rocket Motor (SRM). THe uncertainties arising from propellant characterization through specimem testing led to the decision to perform a model survey and model correlation of a single segment of the Shuttle SRM. Multiple input methods were used to excite and define case/propellant modes of both an inert segment and, later, a live propellant segment. These tests were successful at defining highly damped, flexible modes, several pairs of which occured with frequency spacing of less than two percent.
Shape optimization of unconstrained and constrained damping layers is completed. The specific problem analyzed is a cantilever beam loaded at its tip by a harmonic force. Finite element modeling and mathematical programming techniques are used to obtain the solution. Performance measures are taken to be reduction of maximum diplacement and increase in fatigue lifetime. Results include the improvement, over the uniform treatment case, of these measures when the profile of the damping layer is optimized.
This report covers the design analysis of a very thin, highly swept, propeller blade to be used in the Large-Scale Advanced Prop-Fan (LAP) test program. The report includes: design requirements and goals, a description of the blade configuration which meets requirements, a description of the analytical methods utilized/developed to demonstrate compliance with the requirements, and the results of these analyses. The methods described include: finite element modeling, predicted aerodynamic loads and their application to the blade, steady state and vibratory response analyses, blade resonant frequencies and mode shapes, bird impact analysis, and predictions of stalled and unstalled flutter phenomena. Summarized results include deflections, retention loads, stress/strength comparisons, foreign object damage resistance, resonant frequencies and critical speed margins, resonant vibratory mode shapes, calculated boundaries of stalled and unstalled flutter, and aerodynamic and acoustic performance calculations.
The requirement for low vibrations has achieved the status of a critical design consideration in modern helicopters. There is now a recognized need to account for vibrations during both the analytical and experimental phases of design. Research activities in this area were both broad and varied and notable advances were made in recent years in the critical elements of the technology base needed to achieve the goal of a jet smooth ride. The purpose is to present an overview of accomplishments and current activities of govern and government-sponsored research in the area of rotorcraft vibrations and structural dynamics, focusing on NASA and Army contributions over the last decade or so. Specific topics addressed include: airframe finite-element modeling for static and dynamic analyses, analysis of coupled rotor-airframe vibrations, optimization of airframes subject to vibration constraints, active and passive control of vibrations in both the rotating and fixed systems, and integration of testing and analysis in such guises as modal analysis, system identification, structural modification, and vibratory loads measurement.
To evaluate the potential for spacecraft damage caused by activation of pyrotechnic mechanisms, pyrotechnic shock tests were conducted on three configurations: pin-pullers on an orthogonal double Hopkinson bar arrangement, pin-pullers on a mock-up of the HALOE (Halogen Occultation Experiment) structure, and a section of separation joint on a single Hopkinson bar. Strains and accelerations were measured. The strains were converted to output stresses, forces and moments. Acceleration shock response spectra are presented for both acceleration- and force-time signals. The devices were attached to the Hopkinson bars with adaptors, which are typical of attachments used in practice. In order to predict the effect of the adaptors, finite element models of the adaptor-bar combinations were analyzed and results are presented with regard to grid size, time interval, and frequency range required to predict the response to pulses of 10 and 100 microsecond duration.
The Centaur G Prime modal test resulted in sets of modes (frequencies, mode shapes and damping) with an accuracy similiar to or better than that normally obtained from the modal testing of linear structures with no backlash and small damping. In other words, performing the test at high level greatly minimized the backlash effect and provided a valid, simple linearization of the trunnion friction problem for the Centaur in the Shuttle Cargo Bay. All the most important modes (target modes) were measured and provided the data base for updating the finite element model for the pre-flight verification loads analysis.
A 3-D finite element model was developed to analyze the fiber-resin behavior near the delamination front in a graphite-epoxy double cantilever beam (DCB) specimen. The specimen interior was analyzed using a typical one-fiber slice, represented by a local 3-D fiber-resin model. The resin stresses were computed for the resin-rich layer at the ply interface as well as for the regions between the fibers close to the delamination front. However, the computed strain energy release rate G sub I along the delamination front varied by less than two percent, and was within about four percent of the plane-strain value. The von Mises yield criterion was used to estimate the extent of yielding near the delamination front. The yielding extended ahead of the delamination and also developed between the fibers. Although the fibers had only a negligible effect on G sub I, they caused yielding within the ply and therefore could influence delamination fracture toughness. The normal and shear stresses at the fiber-resin interface were computed near the delamination front. These results suggest that multi-axial stress criteria may be required to analyze fiber-resin interfaces.