ALBERT

Description: Control system experiments are described that were carried out on flexible structure facilities at the Jet Propulsion Laboratory (JPL) and the Wright-Patterson Air Force Base (WPAFB). The feedback controllers were designed using the technique of weighted gap optimization to increase structural damping and achieve robust control.

Description: Tests were performed measuring the locations and geometries of delaminations in Fiberite T300/976 graphite/epoxy, Fiberite IM7/977-2 graphite-toughened epoxy, and ICI APC-2 graphite/PEEK plates subjected to transverse impact loads. The data provide specific information on the effects of impactor velocity, impactor mass, material, thickness of back ply group, difference in fiber orientation between adjacent ply groups, plate thickness, and impactor nose radius. The data were compared to the results of the Finn-Springer model. The model was found to describe the data with reasonable accuracy.

Description: The results of a study of the effects of impact damage on compression-loaded trapezoidal-corrugation sandwich and semi-sandwich graphite-thermoplastic panels are presented. Sandwich panels with two identical face sheets and a trapezoidal corrugated core between them, and semi-sandwich panels with a corrugation attached to a single skin are considered in this study. Panels were designed, fabricated and tested. The panels were made using the manufacturing process of thermoforming, a less-commonly used technique for fabricating composite parts. Experimental results for unimpacted control panels and panels subjected to impact damage prior to loading are presented. Little work can be found in the literature about these configurations of thermoformed panels.

Description: The conditions under which global and local singularities may arise in structural optimization are examined. Examples of these singularities are presented, and a framework is given within which the singularities can be recognized. It is shown, in particular, that singularities can be identified through the analysis of stress-displacement relations together with compatibility conditions or the displacement-stress relations derived by the integrated force method of structural analysis. Methods of eliminating the effects of singularities are suggested and illustrated numerically.

Description: The vibration isolation problem is formulated as a 1D kinematic problem. The geometry of the stochastic wall trajectories arising from the stroke constraint is defined in terms of their significant extrema. An optimal control solution for the minimum acceleration return path determines a lower bound on platform mean square acceleration. This bound is expressed in terms of the probability density function on the significant maxima and the conditional fourth moment of the first passage time inverse. The first of these is found analytically while the second is found using a Monte Carlo simulation. The rms acceleration lower bound as a function of available space is then determined through numerical quadrature.

Description: An account is given of an inspection method that has been successfully used to assess the postimpact damage and residual strength of syntactic (glass microspheres in epoxy matrix) foam-core sandwich panels with hybrid (carbon and glass fiber-reinforced) composite skins, which inherently possess high damage tolerance. SEM establishes that the crushing of the microspheres is responsible for the absorption of most of the impact energy. Damage tolerance is a function of the localization of damage by that high impact energy absorption.

Description: A new mathematical approach for quantifying the mechanical properties of elastomeric materials under biaxial loads is presented. Specific equations relating principal strains and principal true stresses for a homogeneous, isotropic, and elastic material are proposed that resemble the conventional Hooke's law of linear elasticity. The predicted stresses are compared to those from three different sets of experimental data and to stresses from three different theories based on Rivlin's work-of-deformation approach. The proposed approach is considered to be of benefit to design engineers involved in a broad range of rubber products.

Description: A method is presented for calculating the locations, shapes, and sizes of delaminations which occur in a fiber reinforced composite plate subjected to transverse static or dynamic (impact) loads. The plate may be simply supported, clamped, or free along its edges. A model of the delamination formation was developed. This model was then coupled with a finite element analysis. The model and the finite element analysis were implemented by a computer code which can be used to estimate the load at which damage is initiated as well as the locations, shapes, and sizes of the delaminations.

Description: Nonlinear finite-element structural analyses were performed on the first stage high-pressure fuel turbopump blade of the Space Shuttle Main Engine. The analyses examined the structural response and the dynamic characteristics at typical operating conditions. Single crystal material PWA-1480 was considered for the analyses. Structural response and the blade natural frequencies with respect to the crystal orientation were investigated. The analyses were conducted based on typical test stand engine cycle. Influence of combined thermal, aerodynamic, and centrifugal loadings was considered. Results obtained showed that the single crystal secondary orientation effects on the maximum principal stresses are not highly significant.

Description: A method for deriving constrained or fixed-base modes and frequencies from free-free modes of a structure with mass-loaded boundaries is developed. Problems associated with design and development of test fixtures can be avoided with such an approach. The analytical methodology presented is used to assess applicability of the mass-additive method for three types of structures and to determine the accuracy of derived constrained modes and frequencies. Results show that mass loading of the boundaries enables local interface modes to be measured within a desired frequency bandwidth, thus allowing constrained modes to be derived with considerably fewer free-free modes than for unloaded boundaries. Good convergence was obtained for a simple beam and a truss-like Shuttle payload, both of which had well-spaced modes and stiff interface support structures. Slow convergence was obtained for a space station module prototype, a shell-like structure having high modal density.

Description: An algorithm for a general, multilevel structural optimization by substructuring is derived, based on the linear decomposition concept that is rooted in the Bellman's Optimality Criterion enhanced with the optimum sensitivity derivatives used as a means to account for coupling among the subproblems, each of which is limited to optimization of a substructure. The algorithm applies also to those multidisciplinary problems whose subproblems form a hierarchy similar to that of substructures. In systems where the subproblems communicate with each other at the same level, the decomposition becomes non-hierarchic and the system may be optimized as a whole based on the derivatives of the system behavior with respect to the design variables computed by a method that bypasses finite differencing on the system analysis. When a multidisciplinary system includes a structure as its part, a hybrid, hierarchic/non-hierarchic decomposition applies. Numerical examples and references to computational experience accumulated to date illustrate the discussion.

Description: This paper discusses two schemes for doing finite element K calibration in the frequency domain. The baseline scheme uses the definition of K as a limit toward the crack tip. The limiting process requires using a very fine mesh around the crack tip making the scheme computationally very expensive. It is shown that the behavior of K as a function of frequency is very similar to a modal response. Taking advantage of this, a more efficient scheme involves a modal analysis of the cracked sheet and scaling the response to that of the static stress intensity. In this way, only a static K calibration need be performed. All the examples shown are for a frequency range spanning multiple resonances and with two levels of damping.

Description: An analytical calibration of the Stable Poisson Loaded (SPL) specimen is presented. The specimen configuration is similar to the ASTM E-561 compact-tension specimen with displacement controlled wedge loading used for R-curve determination. The crack mouth opening displacements (CMODs) are produced by the diametral expansion of an axially compressed cylindrical pin located in the wake of a machined notch. Due to the unusual loading configuration, a three-dimensional finite element analysis was performed with gap elements simulating the contact between the pin and specimen. In this report, stress intensity factors, CMODs, and crack displacement profiles, are reported for different crack lengths and different contacting conditions. It was concluded that the computed stress intensity factor decreases sharply with increasing crack length thus making the SPL specimen configuration attractive for fracture testing of brittle, high modulus materials.

Description: This paper presents a design methodology for a laminated composite stiffened panel, subjected to multiple in-plane loads and bending moments. Design variables include the skin and stiffener ply orientation angles and stiffener geometry variables. Optimum designs are sought which minimize structural weight and satisfy mechanical performance requirements. Two types of mechanical performance requirements are placed on the panel, maximum strain and minimum strength. Minimum weight designs are presented which document that the choice of mechanical performance requirements cause changes in the optimum design. The effects of lay-up constraints which limit the ply angles to user specified values, such as symmetric or quasi-isotropic laminates, are also investigated.

Description: Recent advances in structural design for control are reviewed. Attention is given to adaptive structures, passive damping, ground testing, and system identification. Directions for future research are suggested.

Description: The benefits of structural dynamics modeling methods in aerospace structures are reviewed. Four major issues in structural dynamics modeling are discussed which encompass most of its subdisciplines: reduced order modeling, constraints in problems with large motions, computational strategies, and fundamental methods. Directions for future research in these areas are addressed.

Description: To simulate the dynamical motion of articulated, multiflexible body structures, one can use multibody simulation packages such as DISCOS. To this end, one must supply appropriate reduced-order models for all of the flexible components involved. The component modes projection and assembly model reduction (COMPARE) methodology is one way to construct these reduced-order component models, which when reassembled capture important system input-to-output mapping of the full-order model at multiple system configurations of interest. In conjunction, we must also supply component damping matrices which when reassembled generate a system damping matrix that has certain desirable properties. The problem of determining the damping factors of components' modes to achieve a given system damping matrix is addressed here. To this end, we must establish from first principles a matrix-algebraic relation between the system's modal damping matrix and the components' modal damping matrices. An unconstrained/constrained optimization problem can then be formulated to determine the component modes' damping factors that best satisfy that matrix-algebraic relation. The effectiveness of the developed methodology, called ModeDamp, has been successfully demonstrated on a high-order, finite element model of the Galileo spacecraft.

Description: This paper presents an alternative to the correlation of individual components of a mode shape vectors by directly examining the sensitivity of the cross-orthogonality between test and analytical mode shapes. If the test and analysis mode shapes are identical, the diagonal elements of the cross-orthogonality will be identical to the test orthogonality matrix, so the cross-orthogonality matrix provides a concise measure of the 'closeness' between test and analysis mode shapes. There are two major advantages to the cross-orthogonality correlation approach. The first is that a direct correlation of this matrix will more directly meet the goal of the correlation effort (measured by cross-orthogonality). Secondly, and more importantly, the correlation of cross-orthogonality greatly reduces the amount of data that needs to be handled when compared to the correlation of mode shapes.

Description: Random vibration response data for orbiting geophysical observatory - comparative analysis of flight, acoustic, and vibration tests

Description: Development and verification of vibration test requirements for Apollo command and service modules

Description: Acoustic test conditions to verify structural and component integrity of Apollo lunar module

Description: Acoustic qualification testing of aerodynamic fins used on S-1C stage of Saturn 5 vehicle

Description: Saturn S-2, S-4B, and instrument unit subassembly and assembly vibration and acoustic evaluation programs

Description: A finite element solution for the structural behavior of a scientific balloon has been obtained using a non-linear finite element code. The pneumatic skin is modelled by shell elements that are given a small artificial bending stiffness to overcome numerical problems yet the membrane solution remains unaffected. Validation of the analysis approach is provided through strain measurements on a small scale balloon that exhibits all essential features of a full scale balloon.

Description: Amplitude frequency relations of nonlinear vibrations in uniform beams with various boundary conditions using perturbation method

Description: Polymide resin/glass fiber sandwich panels studied for high temperature effects on mechanical properties for application to SST

Description: Simulated micrometeoroid bombardment of white thermal control coating causing changes in optical properties

Description: Thin plate ballistic limit for perforation by cylindrical projectiles, using computational procedures and graphs

Description: Edge restraint coupling effect on buckling mode of ring stiffened shells of revolution

Description: Vibrational characteristics of ring stiffened cylindrical shells including inter-ring motion, comparing results with energy methods and orthotropic approach

Description: Elastic buckling and postbuckling behavior of thin doubly curved shell panels subject to nonuniform heating, noting load vs deflection

Description: Steady state frequency response of conical and cylindrical shells under lateral excitation

Description: Elastic displacements for edge cracked plate specimens used as crack extension indicator in plane strain fracture toughness measurements

Description: Kinetic theory of fracture initiation, discussing changes in relative orientation distribution of network chains

Description: The Long Duration Exposure Facility (LDEF) micrometeoroid/space debris impact data has been reduced in terms that are convenient for evaluating the overall quantitative effect on material properties. Impact crater flux has been evaluated as a function of angle from velocity vector and as a function of crater size. This data is combined with spall data from flight and ground testing to calculate effective solar absorption and emittance values versus time. Results indicate that the surface damage from micrometeoroid/space debris does not significantly affect the overall surface optical thermal physical properties. Of course the local damage around impact craters radically alter optical properties. Damage to composites and solar cells on an overall basis was minimal.

Description: Prediction of whether a pressurized cylinder will fail catastrophically when impacted by a projectile has important applications ranging from perforation of an airplane's skin by a failed turbine blade to meteorite impact of a space station habitation module. This report summarizes the accomplishment of one task for a project whose aim is to simulate numerically the outcome of a high velocity impact of pressure vessels. A finite element patch covering the vicinity of a growing crack has been constructed to estimate the J-integral (crack driving force) during the impact. Explicit expressions for the J-integral through the nodal values of displacement, strain, and stress have been written. The patch is to be used repeatedly to estimate the amount of crack growth during the time of the impact. The resulting crack size is to be compared to an estimated critical crack size for the pressurized cylinder.

Description: A brief overview of the transient dynamics capabilities at Sandia National Laboratories, with an emphasis on recent new developments and current research is presented. In addition, the Sandia National Laboratories (SNL) Engineering Analysis Code Access System (SEACAS), which is a collection of structural and thermal codes and utilities used by analysts at SNL, is described. The SEACAS system includes pre- and post-processing codes, analysis codes, database translation codes, support libraries, Unix shell scripts for execution, and an installation system. SEACAS is used at SNL on a daily basis as a production, research, and development system for the engineering analysts and code developers. Over the past year, approximately 190 days of CPU time were used by SEACAS codes on jobs running from a few seconds up to two and one-half days of CPU time. SEACAS is running on several different systems at SNL including Cray Unicos, Hewlett Packard PH-UX, Digital Equipment Ultrix, and Sun SunOS. An overview of SEACAS, including a short description of the codes in the system, are presented. Abstracts and references for the codes are listed at the end of the report.

Description: Numerical simulation of vehicle crashworthiness and occupant protection are addressed. The vehicle crashworthiness design objectives are to design the vehicle structure for optimum impact energy absorption, and to design the restraint system (seatbelts, airbags, bolsters, etc.) for optimum occupant protection. The following approaches are taken; a major part of the impact energy is to be absorbed by the vehicle structure; the restraint components will provide protection against the remaining crash energy; certain vehicle components are designed to deform under specific types and speeds of impact in a desired mode for sound energy management; structural components such as front side rails, rear rails, door structure and pillars undergo large amounts of deformation; and with properly designed geometry and material these components assist in mitigating the effects of impact.

Description: The impact analysis of composite aircraft structures is discussed. Topics discussed include: background remarks on aircraft crashworthiness; comments on modeling strategies for crashworthiness simulation; initial study of simulation of progressive failure of an aircraft component constructed of composite material; and research direction in composite characterization for impact analysis.

Description: Due to the unavailability and, later, prohibitive cost of the computational power required, many phenomena in nonlinear dynamic systems have in the past been addressed in terms of linear systems. Linear systems respond to periodic inputs with periodic outputs, and may be characterized in the time domain or in the frequency domain as convenient. Reduction to the frequency domain is frequently desireable to reduce the amount of computation required for solution. Nonlinear systems are only soluble in the time domain, and may exhibit a time history which is extremely sensitive to initial conditions. Such systems are termed chaotic. Dynamic buckling, aeroelasticity, fatigue analysis, control systems and electromechanical actuators are among the areas where chaotic vibrations have been observed. Direct transient analysis over a long time period presents a ready means of simulating the behavior of self-excited or externally excited nonlinear systems for a range of experimental parameters, either to characterize chaotic behavior for development of load spectra, or to define its envelope and preclude its occurrence.

Description: The University of Virginia examined the design of actuators for both single-degree-of-freedom (SDOF) and multiple-degree-of-freedom (MDOF) active microgravity isolation systems. For SDOF systems, two actuators were considered: a special large gap magnetic actuator and a large stroke Lorentz actuator. The magnetic actuator was viewed to be of greater difficulty than the Lorentz actuator with little compelling technical advantage and was dropped from consideration. A Lorentz actuator was designed and built for the SDOF test rig using magnetic circuit and finite element analysis. The design and some experimental results are discussed. The University also examined the design of actuators for MDOF isolation systems. This includes design of an integrated 1 cm gap 6-DOF noncontacting magnetic suspension system and of a 'coarse' follower which permits the practical extension of magnetic suspension to large strokes. The proposed 'coarse' actuator was a closed kinematic chain manipulator known as a Stewart Platform. The integration of the two isolation systems together, the isolation tasks assigned to each, and possible control architectures were also explored. The results of this research are examined.

Description: Composite structures have the potential to be cost-effective, structurally efficient primary aircraft structures. The Advanced Composites Technology (ACT) Program has the goal to develop the technology to exploit this potential for heavily loaded aircraft structures. As part of the ACT Program, Lockheed Aeronautical Systems Company completed the design and fabrication of the Technology Integration Box Beam (TIBB). The TIBB is an advanced composite prototype structure for the center wing section of the C-130 aircraft. Lockheed subjected the TIBB to downbending, upbending, torsion and combined upbending and torsion load conditions to verify the design. The TIBB failed at 83 percent of design ultimate load for the combined upbending and torsion load condition. The objective of this paper is to describe the mechanisms that led to the failure of the TIBB. The results of a comprehensive analytical and experimental study are presented. Analytical results include strain and deflection results from both a global analysis of the TIBB and a local analysis of the failure region. These analytical results are validated by experimental results from the TIBB tests. The analytical and experimental results from the TIBB tests are used to determine a sequence of events that resulted in failure of the TIBB. A potential cause of failure is high stresses in a stiffener runout region. Analytical and experimental results are also presented for a stiffener runout specimen that was used to simulate the TIBB failure mechanisms.

Description: Approximately 20.4 sq m of Teflon thermal blankets on the nonspinning Long Duration Exposure Facility (LDEF) were exposed to the orbital debris and micrometeoroid environment in low-Earth orbit (LEO) for approximately 5.7 years. Each blanket consisted of an outer layer (approximately 125 micron thick) of FEP Teflon that was backed by a vapor-deposited metal mirror (Inconel; less than 1 micron thick). The inner surface consisted of organic binders and Chemglaze thermal protective paint (approximately 50 micron thick) resulting in a somewhat variable, total blanket thickness of approximately 180 to 200 microns. There was at least one of these blankets, each exposing approximately 1.2 sq m of surface area, on nine of LDEF's 12 principal pointing directions, the exceptions being Rows 3, 9, and 12. As a consequence, these blankets represent a significant opportunity for micrometeoroid and debris studies, in general, and specifically they provide an opportunity to address those issues that require information about pointing direction (i.e., spatial density of impact events as a function of instrument orientation). During deintegration of the LDEF spacecraft at KSC, all penetration holes greater than or equal to 300 micron in diameter were documented and were recently synthesized in terms of spatial density as a function of LDEF viewing direction by. The present report describes ongoing cratering and penetration experiments in pure Teflon targets, which are intended to establish the relationships between crater or penetration-hole diameters and the associated projectile dimensions at laboratory velocities (i.e., 6 km/s). The ultimate objective of these efforts is to extract reliable mass-frequencies and associated fluxes of hypervelocity particles in LEO.

Description: For practical reasons, experimental studies of collisional fragmentation must at times rely on explosives to fragment a target body. For example, Housen et al., described experiments in which spheres were fragmented in a pressurized atmosphere. Explosives were used because impacts could not be performed in the pressure chamber. Explosives can also be used to study targets much larger than those which can be disrupted by conventional light-gas guns, thereby allowing size- and rate-effects to be investigated. The purpose of this study is to determine the charge burial depth required to simulate various aspects of collisions.

Description: Dust-sized olivine particles were fired at a copper plate using the Space Power Institute hypervelocity facility, simulating micrometeoroid damage from natural debris to spacecraft in low-Earth orbit (LEO). Techniques were developed for measuring crater volume, particle volume, and particle velocity, with the particle velocities ranging from 5.6 to 8.7 km/s. A roughly linear correlation was found between crater volume and particle energy which suggested that micrometeoroids follow standard hypervelocity relationships. The residual debris analysis showed that for olivine impacts of up to 8.7 km/s, particle residue is found in the crater. By using the Space Power Institute hypervelocity facility, micrometeoroid damage to satellites can be accurately modeled.

Description: The floor of the Core Module Simulator (CMS) is required to support various combinations of dead load and live load during the testing process. Even though there is published data on the structural capability of the grating, it is not always evident if the combined loadings with joint loads will cause structural failure.

Description: NASA Langley developments in response calculations needed for failure and life predictions are discussed. Topics covered include: structural failure analysis in concurrent engineering; accuracy of independent regional modeling demonstrated on classical example; functional interface method accurately joins incompatible finite element models; interface method for insertion of local detail modeling extended to curve pressurized fuselage window panel; interface concept for joining structural regions; motivation for coupled 2D-3D analysis; compression panel with discontinuous stiffener coupled 2D-3D model and axial surface strains at the middle of the hat stiffener; use of adaptive refinement with multiple methods; adaptive mesh refinement; and studies on quantity effect of bow-type initial imperfections on reliability of stiffened panels.

Description: A recently developed high-temperature fatigue life prediction computer code is presented and an example of its usage given. The code discussed is based on the Total Strain version of Strainrange Partitioning (TS-SRP). Included in this code are procedures for characterizing the creep-fatigue durability behavior of an alloy according to TS-SRP guidelines and predicting cyclic life for complex cycle types for both isothermal and thermomechanical conditions. A reasonably extensive materials properties database is included with the code.

Description: Three parallel computational simulation methods are being developed at the LeRC Structural Mechanics Branch (SMB) for composite structures failure and life analysis: progressive fracture CODSTRAN; hierarchical methods for high-temperature composites; and probabilistic evaluation. Results to date demonstrate that these methods are effective in simulating composite structures failure/life/reliability.

Description: Titanium metal matrix composites are being evaluated for structural applications on advanced hypersonic vehicles. These composites are reinforced with ceramic fibers such as silicon carbide, SCS-6. This combination of matrix and fiber results in a high stiffness, high strength composite that has good retention of properties even at elevated temperatures. However, significant thermal stresses are developed within the composite between the fiber and the matrix due to the difference in their respective coefficients of thermal expansion. In addition to the internal stresses that are generated due to thermal cycling, the overall laminate will be subjected to considerable mechanical loads during the thermal cycling. In order to develop life prediction methodology, one must be able to predict the stresses and strains that occur in the composite's constituents during the complex loading. Thus the purpose is to describe such an analytical tool, VISCOPLY.

Description: With the advent of advanced materials in rotating gas turbine engine components, the methodologies for life prediction of these parts must also increase in sophistication and capability. Pratt & Whitney's view of generic requirements for composite component life prediction systems are presented, efforts underway to develop these systems are discussed, and industry participation in key areas requiring development is solicited.

Description: An analytical model was developed for predicting the response of laminated composites with or without a cutout and subjected to in-plane tensile and shear loads. Material damage resulting from the loads in terms of matrix cracking, fiber-matrix shearing, and fiber breakage was considered in the model. Delamination, an out-of-plane failure mode, was excluded from the model.

Description: The need for enhanced and improved performance of structural components subject to severe cyclic thermal/mechanical loadings, such as in the aerospace industry, requires development of appropriate solution technologies involving time-dependent inelastic analyses. Such analyses are mandatory to predict local stress-strain response and to assess more accurately the cyclic life time of structural components. The NASA-Lewis Research Center is cognizant of this need. As a result of concerted efforts at Lewis during the last few years, several such finite element solution technologies (in conjunction with the finite element program MARC) were developed and successfully applied to numerous uniaxial and multiaxial problems. These solution technologies, although developed for use with MARC program, are general in nature and can easily be extended for adaptation with other finite element programs such as ABAQUS, ANSYS, etc. The description and results obtained from two such inelastic finite element solution technologies are presented. The first employs a classical (non-unified) creep-plasticity model. An application of this technology is presented for a hypersonic inlet cowl-lip problem. The second of these technologies uses a unified creep-plasticity model put forth by Freed. The structural component for which this finite element solution technology is illustrated, is a cylindrical rocket engine thrust chamber. The advantages of employing a viscoplastic model for nonlinear time-dependent structural analyses are demonstrated. The life analyses for cowl-lip and cylindrical thrust chambers are presented. These analyses are conducted by using the stress-strain response of these components obtained from the corresponding finite element analyses.

Description: The results of the so-called energetic approach to fracture with particular attention to the issue of energy dissipation due to crack propagation are applied to the case of a crack with cohesive zone. The thermodynamic admissibility of subcritical crack growth (SCG) is discussed together with some hypotheses that lead to the derivation of SCG laws. A two-phase cohesive zone model for discontinuous crack growth is presented and its thermodynamics analyzed, followed by an example of its possible application.

Description: The evolution of high-temperature, creep-fatigue, life-prediction methods used for cyclic crack initiation is traced from inception in the late 1940's. The methods reviewed are material models as opposed to structural life prediction models. Material life models are used by both structural durability analysts and by material scientists. The latter use micromechanistic models as guidance to improve a material's crack initiation resistance. Nearly one hundred approaches and their variations have been proposed to date. This proliferation poses a problem in deciding which method is most appropriate for a given application. Approaches were identified as being combinations of thirteen different classifications. This review is intended to aid both developers and users of high-temperature fatigue life prediction methods by providing a background from which choices can be made. The need for high-temperature, fatigue-life prediction methods followed immediately on the heels of the development of large, costly, high-technology industrial and aerospace equipment immediately following the second world war. Major advances were made in the design and manufacture of high-temperature, high-pressure boilers and steam turbines, nuclear reactors, high-temperature forming dies, high-performance poppet valves, aeronautical gas turbine engines, reusable rocket engines, etc. These advances could no longer be accomplished simply by trial and error using the 'build-em and bust-em' approach. Development lead times were too great and costs too prohibitive to retain such an approach. Analytic assessments of anticipated performance, cost, and durability were introduced to cut costs and shorten lead times. The analytic tools were quite primitive at first and out of necessity evolved in parallel with hardware development. After forty years more descriptive, more accurate, and more efficient analytic tools are being developed. These include thermal-structural finite element and boundary element analyses, advanced constitutive stress-strain-temperature-time relations, and creep-fatigue-environmental models for crack initiation and propagation. The high-temperature durability methods that have evolved for calculating high-temperature fatigue crack initiation lives of structural engineering materials are addressed. Only a few of the methods were refined to the point of being directly useable in design. Recently, two of the methods were transcribed into computer software for use with personal computers.

Description: The ideas associated with the subsequent viewgraphs are summarized. The primary motivation behind this presentation is to observe that certain macroscopic, microscopic, and submicroscopic phenomena are being understood that have basic influence on understanding the durability and high temperature sensitivity of polymers and polymer-based composites. This understanding covers important issues of long term stability with respect to residual stresses and deformations which can have very deleterious effects on structures used for long periods of time as a result of the heat-involving manufacturing process. Beyond this, important progress is being made in understanding the nonlinear material response of polymers in the fracture context, because the nonlinear mechanics of the material at the tip of a crack, either propagating or ready to do so, is being understood with increasing precision.

Description: The purpose is to provide an end-user's perspective on the state of the art in life prediction and failure analysis by focusing on subsonic transport fuselage issues being addressed in the NASA/Boeing Advanced Technology Composite Aircraft Structure (ATCAS) contract and a related task-order contract. First, some discrepancies between the ATCAS tension-fracture test database and classical prediction methods is discussed, followed by an overview of material modeling work aimed at explaining some of these discrepancies. Finally, analysis efforts associated with a pressure-box test fixture are addressed, as an illustration of modeling complexities required to model and interpret tests.

Description: The explicit transient dynamics technology in use today for simulating the impact and subsequent transient dynamic response of a structure has its origins in the 'hydrocodes' dating back to the late 1940's. The growth in capability in explicit transient dynamics technology parallels the growth in speed and size of digital computers. Computer software for simulating the explicit transient dynamic response of a structure is characterized by algorithms that use a large number of small steps. In explicit transient dynamics software there is a significant emphasis on speed and simplicity. The finite element technology used to generate the spatial discretization of a structure is based on a compromise between completeness of the representation for the physical processes modelled and speed in execution. That is, since it is expected in every calculation that the deformation will be finite and the material will be strained beyond the elastic range, the geometry and the associated gradient operators must be reconstructed, as well as complex stress-strain models evaluated at every time step. As a result, finite elements derived for explicit transient dynamics software use the simplest and barest constructions possible for computational efficiency while retaining an essential representation of the physical behavior. The best example of this technology is the four-node bending quadrilateral derived by Belytschko, Lin and Tsay. Today, the speed, memory capacity and availability of computer hardware allows a number of the previously used algorithms to be 'improved.' That is, it is possible with today's computing hardware to modify many of the standard algorithms to improve their representation of the physical process at the expense of added complexity and computational effort. The purpose is to review a number of these algorithms and identify the improvements possible. In many instances, both the older, faster version of the algorithm and the improved and somewhat slower version of the algorithm are found implemented together in software. Specifically, the following seven algorithmic items are examined: the invariant time derivatives of stress used in material models expressed in rate form; incremental objectivity and strain used in the numerical integration of the material models; the use of one-point element integration versus mean quadrature; shell elements used to represent the behavior of thin structural components; beam elements based on stress-resultant plasticity versus cross-section integration; the fidelity of elastic-plastic material models in their representation of ductile metals; and the use of Courant subcycling to reduce computational effort.

Description: The objective is to describe three research thrusts in crashworthiness analysis: adaptivity; mixed time integration, or subcycling, in which different timesteps are used for different parts of the mesh in explicit methods; and methods for contact-impact which are highly vectorizable. The techniques are being developed to improve the accuracy of calculations, ease-of-use of crashworthiness programs, and the speed of calculations. The latter is still of importance because crashworthiness calculations are often made with models of 20,000 to 50,000 elements using explicit time integration and require on the order of 20 to 100 hours on current supercomputers. The methodologies are briefly reviewed and then some example calculations employing these methods are described. The methods are also of value to other nonlinear transient computations.

Description: The automotive industry has used computational methods for crashworthiness since the early 1970's. These methods have ranged from simple lumped parameter models to full finite element models. The emergence of the full finite element models in the mid 1980's has significantly altered the research direction. However, there remains a need for both simple, rapid modeling methods and complex detailed methods. Some directions for continuing research are discussed.

Description: Two components of a structure which are located side by side, will come in contact by certain force and will transfer the compressive force along the contact area. If the force acts in the opposite direction, the elements will separate and no force will be transferred. If this contact is modeled, the load path will be correctly represented, and the load redistribution results in more realistic stresses in the structure. This is accomplished by using different sets of rigid elements for different loading conditions, or by creating multipoint constraint sets. Comparison of these two procedures is presented for a 4 panel unit (PU) stowage drawer installed in an experiment rack in the Spacelab Life Sciences (SLS-2) payload.

Description: The purpose of this study is to create, test and document a procedure to integrate mathematical optimization algorithms with COSMIC NASTRAN. This procedure is very important to structural design engineers who wish to capitalize on optimization methods to ensure that their design is optimized for its intended application. The OPTNAST computer program was created to link NASTRAN and design optimization codes into one package. This implementation was tested using two truss structure models and optimizing their designs for minimum weight, subject to multiple loading conditions and displacement and stress constraints. However, the process is generalized so that an engineer could design other types of elements by adding to or modifying some parts of the code.

Description: In determining the natural modes and frequencies of a linear elastic structure, Guyan reduction is often used to reduce the size of the mass and stiffness matrices and the solution of the reduced system is obtained first. The reduced system modes are then expanded to the size of the original system by using a static transformation linking the retained degrees of freedom to the omitted degrees of freedom. In the present paper, the transformation matrix of Guyan reduction is modified to include additional terms from a series accounting for the inertial effects. However, the inertial terms are dependent on the unknown frequencies. A practical approximation is employed to compute the inertial terms without any iteration. This new transformation is implemented in NASTRAN using a DMAP sequence alter. Numerical examples using a cantilever beam illustrate the necessary condition for allowing a large number of additional terms in the proposed series correction of Guyan reduction. A practical example of a large model of the Plasma Motor Generator module to be flown on a Delta launch vehicle is also presented.

Description: The QUAD4 and TRIA3 elements are the primary plate/shell elements in NASTRAN. These elements enable the user to analyze thin plate/shell structures for membrane, bending and shear phenomena. They are also very new elements in the NASTRAN library. These elements are extremely versatile and constitute a substantially enhanced analysis capability in NASTRAN. However, with the versatility comes the burden of understanding a myriad of modeling implications and their effect on accuracy and analysis quality. The validity of many aspects of these elements were established through a series of benchmark problem results and comparison with those available in the literature and obtained from other programs like MSC/NASTRAN and CSAR/NASTRAN. Never-the-less such a comparison is never complete because of the new and creative use of these elements in complex modeling situations. One of the important features of QUAD4 and TRIA3 elements is the offset capability which allows the midsurface of the plate to be noncoincident with the surface of the grid points. None of the previous elements, with the exception of bar (beam), has this capability. The offset capability played a crucial role in the design of QUAD4 and TRIA3 elements. It allowed modeling layered composites, laminated plates and sandwich plates with the metal and composite face sheets. Even though the basic implementation of the offset capability is found to be sound in the previous applications, there is some uncertainty in relatively simple applications. The main purpose of this paper is to test the integrity of the offset capability and provide guidelines for its effective use. For the purpose of simplicity, references in this paper to the QUAD4 element will also include the TRIA3 element.

Description: This paper explains how NASTRAN can be utilized advantageously in the preliminary design cycle. The initial portion of the preliminary design process lends itself to programs that can produce multiple configurations or variations on a particular design with minimal cost or effort. The latter portion of the process encompasses refining the design and adding more detailed analyses (particularly for other disciplines). A method for quickly generating balanced spacecraft loading conditions for use in preliminary design and analysis also is explained. The following additional sections are included: Background, Symbols, Analytical Process, Aerodynamic Load Distributions, NASTRAN Applications, Conclusion and References.

Description: The following are presented: the experimental procedure and the results that have been performed to determine the axial rigidity of the strut-node joint; and the method for modifying a simple testing machine to make it capable of performing more accurate tests over a specific load range and able to accept larger test assemblies.

Description: The objective of this experiment is to determine the level of tensile strength of uncemented, dry, granular materials. The experimental apparatus does not lend itself to a direct measurement of the material's tensile strength, but must be analyzed as a stress field problem in order to arrive at a tensile strength value. The experiment, and subsequent analysis, serve to instruct the student on the influence of gravitationally induced stresses in frictional granular materials, the importance and difficulty of accurately describing the entire failure envelope for granular materials in the low mean stress range, and the fundamental principles of material modeling.

Description: A method for automated data collection has been developed for a Charpy impact tester. A potentiometer is connected to the pivot point of the hammer and measures the angular displacement of the hammer. This data is collected with a computer and, through appropriate software, accurately records the energy absorbed by the specimen. The device can be easily calibrated with minimal effort.

Description: Strength of sheet metal parts containing cracks studied by simulating fail safe design

Description: HF plate vibration damping due to gas pumping in riveted joints, considering pressure dependence of viscosity

Description: Nonlinear boundary layer for shells of revolution, deriving nonlinear differential equations for infinitesimal strains and small rotation

Description: Plastic strain pulse propagation in cylindrical lead bars, discussing stress-strain-time deformation

Description: Time dependent theory of mechanical strength of homogeneous and partially oriented and stressed media based on microscopic molecular structure

Description: Structural analysis modification by force method accounting for elastic property changes of some elements

Description: Explicit expressions for stiffness matrix of triangular torus elements associated with linear displacement fields and generalized Hookian material

Description: Four pole parameters of truncated thin conical shell under axial excitations obtained by applying membrane shell theory

Description: Approximate buckling load for thin ogive shell of revolution under uniform external pressure using differential equations

Description: Doubly curved shallow thin shells buckling under edge loading, discussing curvatures role and shell equation

Description: Metal alloys high temperature low cycle strain fatigue resistance in creep range estimated from tensile and stress rupture data

Description: Modern military aircraft transparency systems, windshields and canopies, are complex systems which must meet a large and rapidly growing number of requirements. Many of these transparency system requirements are conflicting, presenting difficult balances which must be achieved. One example of a challenging requirements balance or trade is shaping for stealth versus aircrew vision. The large number of requirements involved may be grouped in a variety of areas including man-machine interface; structural integration with the airframe; combat hazards; environmental exposures; and supportability. Some individual requirements by themselves pose very difficult, severely nonlinear analysis problems. One such complex problem is that associated with the dynamic structural response resulting from high energy bird impact. An improved analytical capability for soft-body impact simulation was developed.

Description: During the 1960's over 30 full-scale aircraft crash tests were conducted by the Flight Safety Foundation under contract to the Aviation Applied Technology Directorate (AATD) of the U.S. Army Aviation Systems Command (AVSCOM). The purpose of these tests were to conduct crash injury investigations that would provide a basis for the formulation of sound crash resistance design criteria for light fixed-wing and rotary wing aircraft. This resulted in the Crash Survival Design Criteria Designer's Guide which was first published in 1967 and has been revised numerous times, the last being in 1989. Full-scale aircraft crash testing is an expensive way to investigate structural deformations of occupied spaces and to determine the decelerative loadings experienced by occupants in a crash. This gave initial impetus to the U.S. Army to develop analytical methods to predict the dynamic response of aircraft structures in a crash. It was believed that such analytical tools could be very useful in the preliminary design stage of a new helicopter system which is required to demonstrate a level of crash resistance and had to be more cost effective than full-scale crash tests or numerous component design support tests. From an economic point of view, it is more efficient to optimize for the incorporation of crash resistance features early in the design stage. However, during preliminary design it is doubtful if sufficient design details, which influence the exact plastic deformation shape of structural elements, will be available. The availability of simple procedures to predict energy absorption and load-deformation characteristics will allow the designer to initiate valuable cost, weight, and geometry tradeoff studies. The development of these procedures will require some testing of typical specimens. This testing should, as a minimum, verify the validity of proposed procedures for providing pertinent nonlinear load-deformation data. It was hoped that through the use of these analytical models, the designer could optimize aircraft design for crash resistance from both a weight and cost increment standpoint, thus enhancing the acceptance of the design criteria for crash resistance.

Description: The Landing and Impact Dynamics Branch of NASA Langley Research Center has been involved in impact dynamics research since the early 1970's. For the first ten years, the emphasis of the research was on metal aircraft structures in both the General Aviation Crash Dynamics Program and the Controlled Impact Demonstration (CID) Program, a transport aircraft program culminating in the controlled crash test of a Boeing 720 aircraft in 1984. Subsequent to the transport work, the emphasis has been on composite structures with efforts directed at understanding the behavior, responses, failure mechanisms, and general loads associated with the composite material systems under crash type loadings. Considerable work has been conducted to address the energy absorption characteristics and it indicates that composites can absorb as much if not considerably more energy than comparable aluminum structures. However, due to their brittle nature, attention must be given to proper geometry and designs to take advantage of the good energy absorbing properties while providing desired structural integrity. Achieving the desired new designs often requires an understanding of how more conventional designs behave under crash type loadings. The purpose is to present a review of the composite impact dynamics research being conducted at NASA Langley Research Center. Examples are presented of experimental and analytical data to illustrate the activities in the four program elements of the composite research.

Description: The Magnetic Gimbal Proof-of-Concept Hardware activities, accomplishments, and test results are discussed. The Magnetic Gimbal Fabrication and Test (MGFT) program addressed the feasibility of using a magnetic gimbal to isolate an Electro-Optical (EO) sensor from the severe angular vibrations induced during the firing of divert and attitude control system (ACS) thrusters during space flight. The MGFT effort was performed in parallel with the fabrication and testing of a mechanically gimballed, flex pivot based isolation system by the Hughes Aircraft Missile Systems Group. Both servo systems supported identical EO sensor assembly mockups to facilitate direct comparison of performance. The results obtained from the MGFT effort indicate that the magnetic gimbal exhibits the ability to provide significant performance advantages over alternative mechanically gimballed techniques.

Description: CASES (Controls, Astrophysics and Structures Experiment in Space) is a proposed space experiment to collect x-ray images of the galactic center and solar disk with unprecedented resolution. This requires precision pointing and suppression of vibrations in the long flexible structure that comprises the 32-m x-ray telescope optical bench. Two separate electro-optical sensor systems are provided for the ground test facility (GTF). The Boom Motion Tracker (BMT) measures eigenvector data for post-mission use in system identification. The Tip Displacement Sensor (TDS) measures boom tip position and is used as feedback for the closed-loop control system that stabilizes the boom. Both the BMT and the TDS have met acceptance specifications and were delivered to MSFC in February 1992. This paper describes the sensor concept, the sensor configuration as implemented in the GTF, and the results of characterization and performance testing.

Description: Objectives of the research are: (1) to develop design requirements for damped struts to stabilize control system in the high frequency cross-over and spill-over range; (2) to design, fabricate and test viscously damped strut and viscoelastically damped strut; (3) to verify accuracy of design and analysis methodology of damped struts; and (4) to design and build test apparatus, and develop data reduction algorithm to measure strut complex stiffness. In order to meet the stringent performance requirements of the SPICE experiment, the active control system is used to suppress the dynamic responses of the low order structural modes. However, the control system also inadvertently drives some of the higher order modes unstable in the cross-over and spill-over frequency range. Passive damping is a reliable and effective way to provide damping to stabilize the control system. It also improves the robustness of the control system. Damping is designed into the SPICE testbed as an integral part of the control-structure technology.

Description: The experimental study of a piezoceramic active truss is presented. This active strut is unique in that the piezoceramic configurations allow the stroke length of the strut not to be dependent on the piezoceramic material's expansion range but on the deflection range of the piezoceramic bender segment. A finite element model of a piezoceramic strut segment was constructed. Piezoceramic actuation was simulated using thermally induced strains. This model yielded information on the stiffness and force range of a bender element. The static and dynamic properties of the strut were identified experimentally. Feedback control was used to vary the stiffness of the strut. The experimentally verified model was used to explore implementation possibilities of the strut.

Description: Thermal stress analyses are an important aspect in the development of aerospace vehicles such as the National Aero-Space Plane (NASP) and the High-Speed Civil Transport (HSCT) at NASA-LaRC. These analyses require knowledge of the temperature within the structures which consequently necessitates the need for thermal property data. The initial goal of this research effort was to develop a methodology for the estimation of thermal properties of aerospace structural materials at room temperature and to develop a procedure to optimize the estimation process. The estimation procedure was implemented utilizing a general purpose finite element code. In addition, an optimization procedure was developed and implemented to determine critical experimental parameters to optimize the estimation procedure. Finally, preliminary experiments were conducted at the Aircraft Structures Branch (ASB) laboratory.

Description: A status report of Remote Manipulator System (RMS) active damping augmentation is presented. Topics covered include: active damping augmentation; benefits of RMS ADA; simulated payload definition; sensor and actuator definition; ADA control law design; Shuttle Engineering Simulator (SES) real-time simulation; and astronaut evaluation.

Description: The NASA STRuctural ANalysis (NASTRAN) program is one of the most extensively used engineering applications software in the world. It contains a wealth of matrix operations and numerical solution techniques, and they were used to construct efficient eigenvalue routines. The purpose of this paper is to examine the current eigenvalue routines in NASTRAN and to make efficiency comparisons with a more recent implementation of the Block Lanczos algorithm by Boeing Computer Services (BCS). This eigenvalue routine is now available in the BCS mathematics library as well as in several commercial versions of NASTRAN. In addition, CRAY maintains a modified version of this routine on their network. Several example problems, with a varying number of degrees of freedom, were selected primarily for efficiency bench-marking. Accuracy is not an issue, because they all gave comparable results. The Block Lanczos algorithm was found to be extremely efficient, in particular, for very large size problems.

Description: Two related methods for improving the dependent (OMIT set) displacements after performing a Guyan reduction are presented. The theoretical bases for the methods are derived. The NASTRAN DMAP ALTERs used to implement the methods in a NASTRAN execution are described. Data are presented that verify the methods and the NASTRAN DMAP ALTERs.

Description: This paper describes the interface/integration between FEM/SINDA, a general purpose geometry driven thermal analysis code, and the FEM software: I-DEAS, PATRAN, and NASTRAN. FEM/SINDA brings together the advantages of the finite element method to model arbitrary geometry and anisotropic materials and SINDA's finite difference capability to model thermal properties, loads, and boundary conditions that vary with time or temperature. I-DEAS and PATRAN thermal entities are directly supported since FEM/SINDA uses the nodes of the FEM model as the point at which the temperature is determined. Output from FEM/SINDA (as well as the FEM/SINDA input deck) can be used directly by NASTRAN for structural analysis.

Description: A method for selecting optimum NASTRAN analysis set degrees of freedom for the dynamic eigenvalue problem is described. Theoretical development of the Guyan reduction procedure on which the method is based is first summarized. The algorithm used to select the analysis set degrees of freedom is then developed. Two example problems are provided to demonstrate the accuracy of the algorithm.

Description: The purpose of this paper is to describe the implementation and use of a consistent family of two and three dimensional elements in NASTRAN. The elements which are based on a mixed formulation include a replacement of the original NASTRAN shear element and the addition of triangular quadrilateral shell elements and tetrahedral, pentahedral and hexahedral solid elements. These elements support all static loads including temperature gradient and pressure load. The mass matrix is also generated to support all dynamic rigid formats.

Description: The effectiveness of viscous elements in introducing damping in a structure is a function of several variables including their number, their location in the structure, and their physical properties. In this paper, the optimal damper placement and tuning problem is posed to optimize these variables. Both discrete and continuous optimization problems are formulated and solved corresponding, respectively, to the problems of placement of passive elements and to the tuning of their parameters. The paper particularly emphasizes the critical computational issues resulting from the optimization formulations. Numerical results involving a lightly damped testbed structure are presented.