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: 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: 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: 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.

Description: This research focuses on the computational dynamics of flexible constrained multibody systems. At first a recursive mapping formulation of the kinematical expressions in a minimum dimension as well as the matrix representation of the equations of motion are presented. The method employs Kane's equation, FEM, and concepts of continuum mechanics. The generalized active forces are extended to include the effects of high temperature conditions, such as creep, thermal stress, and elastic-plastic deformation. The time variant constraint relations for rolling/contact conditions between two flexible bodies are also studied. The constraints for validation of MBS simulation of gear meshing contact using a modified Timoshenko beam theory are also presented. The last part deals with minimization of vibration/deformation of the elastic beam in multibody systems making use of time variant boundary conditions. The above methodologies and computational procedures developed are being implemented in a program called DYAMUS.

Description: Since the late 1980's, research in recursive formulations of multibody dynamics has flourished. Historically, much of this research can be traced to applications of low dimensionality in mechanism and vehicle dynamics. Indeed, there is little doubt that recursive order N methods are the method of choice for this class of systems. This approach has the advantage that a minimal number of coordinates are utilized, parallelism can be induced for certain system topologies, and the method is of order N computational cost for systems of N rigid bodies. Despite the fact that many authors have dismissed redundant coordinate formulations as being of order N(exp 3), and hence less attractive than recursive formulations, we present recent research that demonstrates that at least three distinct classes of redundant, nonrecursive multibody formulations consistently achieve order N computational cost for systems of rigid and/or flexible bodies. These formulations are as follows: (1) the preconditioned range space formulation; (2) penalty methods; and (3) augmented Lagrangian methods for nonlinear multibody dynamics. The first method can be traced to its foundation in equality constrained quadratic optimization, while the last two methods have been studied extensively in the context of coercive variational boundary value problems in computational mechanics. Until recently, however, they have not been investigated in the context of multibody simulation, and present theoretical questions unique to nonlinear dynamics. All of these nonrecursive methods have additional advantages with respect to recursive order N methods: (1) the formalisms retain the highly desirable order N computational cost; (2) the techniques are amenable to concurrent simulation strategies; (3) the approaches do not depend upon system topology to induce concurrency; and (4) the methods can be derived to balance the computational load automatically on concurrent multiprocessors. In addition to the presentation of the fundamental formulations, this paper presents new theoretical results regarding the rate of convergence of order N constraint stabilization schemes associated with the newly introduced class of methods.

Description: In the analysis of vibrations of continuous elastic systems, one often encounters complicated transcendental equations with roots directly related to the system's natural frequencies. Typically, these equations contain system parameters whose values must be specified before a numerical solution can be obtained. The present paper presents a method whereby the fundamental frequency can be obtained in analytical form to any desired degree of accuracy. The method is based upon truncation of rapidly converging series involving inverse powers of the system natural frequencies. A straightforward method to developing these series and summing them in closed form is presented. It is demonstrated how Computer Algebra can be exploited to perform the intricate analytical procedures which otherwise would render the technique difficult to apply in practice. We illustrate the method by developing two analytical approximations to the fundamental frequency of a vibrating cantilever carrying a rigid tip body. The results are compared to the numerical solution of the exact (transcendental) frequency equation over a range of system parameters.

Description: A two-stage model reduction methodology, combining the classical Component Mode Synthesis (CMS) method and the newly developed Enhanced Projection and Assembly (EP&A) method, is proposed in this research. The first stage of this methodology, called the COmponent Modes Projection and Assembly model REduction (COMPARE) method, involves the generation of CMS mode sets, such as the MacNeal-Rubin mode sets. These mode sets are then used to reduce the order of each component model in the Rayleigh-Ritz sense. The resultant component models are then combined to generate reduced-order system models at various system configurations. A composite mode set which retains important system modes at all system configurations is then selected from these reduced-order system models. In the second stage, the EP&A model reduction method is employed to reduce further the order of the system model generated in the first stage. The effectiveness of the COMPARE methodology has been successfully demonstrated on a high-order, finite-element model of the cruise-configured Galileo spacecraft.

Description: A useful model for open chains of flexible bodies undergoing large rigid body motions, but small elastic deformations, is one in which the equations of motion are linearized in the small elastic deformations and deformation rates. For slow rigid body motions, the correctly linearized, or consistent, set of equations can be compared to prematurely linearized, or inconsistent, equations and to 'oversimplified,' or ruthless, equations through the use of open loop dynamic simulations. It has been shown that the inconsistent model should never be used, while the ruthless model should be used whenever possible. The consistent and inconsistent models differ by stress stiffening terms. These are due to zeroth-order stresses effecting virtual work via nonlinear strain-displacement terms. In this paper we examine in detail the nature of these stress stiffening terms and conclude that they are significant only when the associated zeroth-order stresses approach 'buckling' stresses. Finally it is emphasized that when the stress stiffening terms are negligible the ruthlessly linearized equations should be used.

Description: In this paper we discuss the issue of geometric stiffening as it arises in the context of multibody dynamics. This topic has been treated in a number of previous publications in this journal and appears to be a debated subject. The controversy revolves primarily around the 'correct' methodology for incorporating the stiffening effect into dynamics formulations. The main goal of this work is to present the different approaches that have been developed for this problem through an in-depth review of several publications dealing with this subject. This is done with the goal of contributing to a precise understanding of the existing methodologies for modelling the stiffening effects in multibody systems. Thus, in presenting the material we attempt to illuminate the key characteristics of the various methods as well as show how they relate to each other. In addition, we offer a number of novel insights and clarifying interpretations of these schemes. The paper is completed with a general classification and comparison of the different approaches.

Description: The low-cycle fatigue behavior of Haynes 188 and B1900+Hf under isothermal, bithermal, and thermomechanical loading conditions has been analyzed on the basis of the total hysteresis energy expended per cycle. It has been observed that in the case of isothermal fatigue the total hysteresis energy correlates well with the fatigue life. In the case of bithermal 'high rate' fatigue, for a given total hysteresis energy per cycle, the fatigue life is equal to or greater than the isothermal fatigue life at the maximum bithermal temperature. This observation could be used to establish a lower bound on life for design purposes. In one case of bithermal creep-fatigue and in thermomechanical fatigue, the life is shorter than that corresponding to the isothermal life at the maximum temperature. The energy supplied, per se, may not always give a systematic correlation with the fatigue life in the cases where time-dependent creep and environmental effects are encountered. Thus, in bithermal creep-fatigue and thermomechanical fatigue, the role of creep and environment and their dependence on the energy supplied have to be properly accounted for before the energy term can be used for life prediction.

Description: This report reviews aeroelastic analyses for propulsion components (propfans, compressors and turbines) being developed and used at NASA LeRC. These aeroelastic analyses include both structural and aerodynamic models. The structural models include a typical section, a beam (with and without disk flexibility), and a finite-element blade model (with plate bending elements). The aerodynamic models are based on the solution of equations ranging from the two-dimensional linear potential equation to the three-dimensional Euler equations for multibladed configurations. Typical calculated results are presented for each aeroelastic model. Suggestions for further research are made. Many of the currently available aeroelastic models and analysis methods are being incorporated in a unified computer program, APPLE (Aeroelasticity Program for Propulsion at LEwis).

Description: Tensile and pure bending tests were conducted on specimens having welded joints made from 2219-T87 aluminum alloy and 2319 filler. Data were collected using photoelastic coatings and strain gages. Stress-strain relationships and contraction ratios were determined at several points in a grid covering the weld material and heat affected zone. Material behavior was nonlinear and nonuniform at all points in the grid and contraction ratios did not conform to those predicted by Chakrabarty's plasticity theory. Yielding in joints made using four new welding procedures was examined. None of the new procedures produced more uniform yielding in the joint.

Description: The merits and limitations of the optimality criteria (OC) method for the minimum weight design of structures subjected to multiple load conditions under stress, displacement, and frequency constraints were investigated by examining several numerical examples. The examples were solved utilizing the Optimality Criteria Design Code that was developed for this purpose at NASA Lewis Research Center. This OC code incorporates OC methods available in the literature with generalizations for stress constraints, fully utilized design concepts, and hybrid methods that combine both techniques. Salient features of the code include multiple choices for Lagrange multiplier and design variable update methods, design strategies for several constraint types, variable linking, displacement and integrated force method analyzers, and analytical and numerical sensitivities. The performance of the OC method, on the basis of the examples solved, was found to be satisfactory for problems with few active constraints or with small numbers of design variables. For problems with large numbers of behavior constraints and design variables, the OC method appears to follow a subset of active constraints that can result in a heavier design. The computational efficiency of OC methods appears to be similar to some mathematical programming techniques.

Description: Two different methods are used in finite element structural analysis to define element mass distribution: consistent and lumped mass approximations. The lumped mass method is currently favored in multibody dynamics software packages because of its computational efficiency. This paper, however, presents a multibody dynamics formulation based on the consistent mass approach which should prove to be as computationally efficient as the lumped mass approach. It is shown that when component flexibility is modelled using modal coordinates, the equations of motion, as formulated by these two methods, are very similar in appearance, thereby facilitating their comparison. It is also shown, through simulation of a simple multibody system, that a system formulated with the lumped mass approximation will often require more finite elements and smaller integration step sizes than the corresponding consistent mass system.

Description: This paper describes the refinements of a previously published method for estimating a full modal damping matrix from complex test modes. It also documents application of the refined method to a structure where complex test modes were derived by the ERA method from multi-input random vibration test data. A numerical example based on simulated test data is presented to demonstrate the validity of the method. The application using real data was not successful, presumably because of noise in the small phase angles of the measured complex modes. Alternative test and data reduction procedures are suggested as possible remedies to the problem. A careful analysis of measurement and data processing errors should be made to examine basic feasibility before implementing the alternative procedures. The ability to estimate a full modal damping matrix is considered important for the preflight estimation of on-orbit damping, and for the synthesis of structural damping from substructure tests.

Description: This paper explores the possibility of performing system identification at the substructure level and then synthesizing the results to obtain a mathematical model for the assembled structure. The study here shows that in order to enforce interface compatibility and equilibrium conditions to the substructure test data, it is necessary to place collocated actuator/sensor pair at every interface degree-of-freedom. Procedures for assembling substructure transfer function data, substructure state-space models, and substructure Markov parameters are presented. Testing difficulties and possible solutions are also discussed. A numerical simulation example is included to illustrate the proposed substructure synthesis methods.

Description: A sensitivity analysis of flutter response of a two-degree of freedom airfoil with plunging and pitching degrees of freedom in transonic flow was performed using a state-space representation of the unsteady aerodynamic behavior. The structural equations of motion of the airfoil with bending and torsional degrees of freedom are coupled to the unsteady airloads, and the aeroelastic system so modeled is solved as an eigenvalue problem to determine the stability. The results of eigenanalysis showed good agreement with flutter calculations performed using a time-integration of the aeroelastic equations. The sensitivities of the flutter speed with respect to the mass and stiffness parameters wee computed by both the analytical and finite difference methods, showing excellent agreement.

Description: A new approach, which allows the physical parameters to be adjusted and preserves the structural connectivity condition, is developed in this paper. This approach is based on the sensitivities of the mass and stiffness matrices with respect to physical parameters in an attempt to use test data to correct the analytical model. The analytical mass and stiffness matrices can be corrected simultaneously. The modal data of the updated model dramatically converges to that of the test model. A numerical example is presented to illustrate the feasibility of this approach. Comparisons with the other methods are also given.

Description: The debond configurations presently subjected to 3D FEM fracture mechanics analyses are respectively of the flange-skin strip and skin-stiffener configuration type. Two methods employing the virtual crack closure technique were used to evaluate the strain energy release rate, or 'G-value' distributions across the debond front. Both methods yielded nearly identical G-value distributions for the debond configurations studied; they were compared with plane strain and shell analyses results from the literature for the flange skin strip configuration, and found to be in good agreement. Mode II is dominant for the skin-stiffener debond configuration.

Description: Verification of analytical models through correlation with ground test results of a complex space truss structure is demonstrated. A multi-component, dynamically scaled space station model configuration is the focus structure for this work. Previously established test/analysis correlation procedures are used to develop improved component analytical models. Integrated system analytical models, consisting of updated component analytical models, are compared with modal test results to establish the accuracy of system-level dynamic predictions. Design sensitivity model updating methods are shown to be effective for providing improved component analytical models. Also, the effects of component model accuracy and interface modeling fidelity on the accuracy of integrated model predictions is examined.

Description: Previous analytical studies have shown the feasibility of using a limited number of active members to adaptively alter the structural shape or behavior to other desired states. As a demonstration of this concept, this paper describes the results of a sequence of validation tests, their correlation with the analytical results, and the problems encountered. An existing full scale, space-erectable, high precision truss structure was used. For the most part, the results showed good agreement with the analysis. However, micron level nonlinearities in the structural behavior were discovered. The significance of the presence of these nonlinearities in precision structures and their impact on the basic premise of adaptivity are discussed.

Description: A composite wing with spars, bulkheads and built-in sensor/control devices is probabilistically configured with a methodology for the probabilistic assessment of smart composite structures. Structural responses such as changes in angle of attack, vertical displacements and stresses in the non-control and control plies are probabilistically assessed to quantify their respective scatter ranges. Sensitivity factors are evaluated to identify those parameters that have the greatest influence on a specific structural response. Results show that smart composite structures can be configured to control distortions and ply stresses, and to have specified scatter ranges in the frequencies, buckling loads in the presence of defects to satisfy specified design requirements.

Description: A sensitivity analysis scheme of the static aeroelastic response of a wing is developed, by incorporating a piecewise panel-based pressure representation into an existing wing aeroelastic model to improve the model's fidelity, including the sensitivity of the wing static aeroelastic response with respect to various shape parameters. The new formulation is quite general and accepts any aerodynamics and structural analysis capability. A program is developed which combines the local sensitivities, such as the sensitivity of the stiffness matrix or the aerodynamic kernel matrix, into global sensitivity derivatives.

Description: The effects of damage on a model space truss are investigated using experimentally measured modes along with a procedure that uses the vibratory residual forces and a weighted sensitivity analysis to estimate the extent of mass and/or stiffness variations, where damage is characterized as a reduction in structural stiffness. It is shown that this approach yields excellent results for a four-bay space truss under consideration and was able to accurately predict not only the location of the damage region but also estimate the relative severity for three different configurations. In addition, the method was able to accurately detect and differentiate a shift in structural mass.

Description: The Timoshenko beam theory was used to model cracked beams and to calculate the total strain energy release rate. The root rotation of the beam segments at the crack tip were estimated based on an approximate 2D elasticity solution. By including the strain energy released due to the root rotations of the beams during crack extension, the strain energy release rate obtained using beam theory agrees very well with the 2D finite element solution. Numerical examples were given for various beam geometries and loading conditions. Comparisons with existing beam models were also given.

Description: Probabilistic predictions based on the Integrated Probabilistic Assessment of Composite Structures (IPACS) code are presented for the material and structural response of unnotched and notched, 1M6/3501-6 Gr/Ep laminates. Comparisons of predicted and measured modulus and strength distributions are given for unnotched unidirectional, cross-ply, and quasi-isotropic laminates. The predicted modulus distributions were found to correlate well with the test results for all three unnotched laminates. Correlations of strength distributions for the unnotched laminates are judged good for the unidirectional laminate and fair for the cross-ply laminate, whereas the strength correlation for the quasi-isotropic laminate is deficient because IPACS did not yet have a progressive failure capability. The paper also presents probabilistic and structural reliability analysis predictions for the strain concentration factor (SCF) for an open-hole, quasi-isotropic laminate subjected to longitudinal tension. A special procedure was developed to adapt IPACS for the structural reliability analysis. The reliability results show the importance of identifying the most significant random variables upon which the SCF depends, and of having accurate scatter values for these variables.

Description: The concept of stiffness tailoring for improved shear buckling resistance of rectangular composite plates is investigated analytically. The tailoring involves only the redistribution of the given material with given orientations to create beneficial stiffening patterns across the planform of the plate. The resulting local nonuniformities in thickness and membrane and bending stiffness combine to change the buckling response of the plate. The weight and average membrane shear stiffness are essentially unaffected by the tailoring. Practical limitations on the degree to which the tailoring may be carried out are shown to govern most designs. Improvements in the shear buckling load on the order of 50 percent are shown possible with monolithic tailoring. Tailored sandwich concepts, in which a light-weight core material is added to keep both plate surfaces flat, can produce improvements well over 100 percent in specific buckling loads compared with uniform composite plates.

Description: Mortar, porous alumina, and sand targets, which were spherical in shape and from 11 to 15 cm in diameter, were impacted normally by plastic (polycarbonate) projectiles of nearly 1 g in mass at velocities about 6 km/s. Fragment velocity at the antipole of impact site (antipodal velocity, V(sub a)), for each experiment, was obtained from two Flash X-ray images recorded prior to and at predetermined delayed time after impact event. It has been revealed that the velocities for the same E/M(sub t) (impact energy divided by target mass) depend strongly on target material, and differ about an order of magnitude between the sand and basalt.

Description: Realizing and understanding the effects of the near-Earth space environment on a spacecraft during its mission lifetime is becoming more important with the regeneration of America's space program. Included among these potential effects are the following: erosion and surface degradation due to atomic oxygen impingement; ultraviolet exposure embrittlement; and delamination, pitting, cratering, and ring formation due to micrometeoroid and debris impacts. These effects may occur synergistically and may alter the spacecraft materials enough to modify the resultant crater, star crack, and/or perforation. This study concentrates on modelling the effects of micrometeoroid and debris hypervelocity impacts into aluminum materials (6061-T6). Space debris exists in all sizes, and has the possibility of growing into a potentially catastrophic problem, particularly since self-collisions between particles can rapidly escalate the number of small impactors. We have examined the morphologies of the Long Duration Exposure Facility (LDEF) impact craters and the relationship between the observed impact damage on LDEF versus the existing models for both the natural (micrometeoroid) and manmade (debris) environments in order to better define these environments.