ALBERT

Description: This paper develops a simplified continuum (continuous with respect to time, stress, etc.) fatigue damage model for use in critical design, Life Extending Control and fault prognosis. The work is based on the local strain cyclic damage modeling method. New nonlinear explicit equation forms of cyclic damage in terms of stress amplitude are derived to facilitate the continuum modelling. Stress based continuum models are derived. Extension to plastic strain-strain rate models is also presented. Progress toward a non-zero mean stress based is presented. Also new nonlinear explicit equation forms in terms of stress amplitude are derived for this case. Application of the various models to design, control, and fault prognosis is considered.

Description: The Photogrammetric Appendage Structural Dynamics Experiment (PASDE) is a Hitchhiker payload scheduled to fly as part of the International Space Station (ISS) Phase-1 flight program to the Russian Space Station Mir. The objective of the first flight of PASDE on STS-74 is to obtain video images of the Mir Kvant-2 solar array response to various structural dynamic excitation events. This experiment will demonstrate the use of photogrammetric techniques for on-orbit structural dynamics measurements. Photogrammetric measurements will provide a low cost alternative to appendage mounted accelerometers to the ISS program. The PASDE experiment hardware consists of three instruments each containing two video cameras, two video tape recorders, a modified video signal time inserter, and associated avionics boxes. The instruments were designed and built at the NASA Langley Research Center, and are integrated into standard Hitchhiker canisters at the NASA Goddard Space Flight Center. The Hitchhiker canisters are then installed into the Space Shuttle cargo bay in locations selected to achieve good video coverage and photogrammetric geometry. The measurement resolution of the instruments is expected to be on the order of 0.25 cm (0.1 in.).

Description: The Boeing Company, under contract to the Electric Power Research Institute (EPRI), has completed a test program on the Mod-2 wind turbines at Goodnoe Hills, Washington. The objectives were to update fatigue load spectra, discern site and machine differences, measure vortex generator effects, and to evaluate rotational sampling techniques. This paper shows the test setup and loads instrumentation, loads data comparisons and test/analysis correlations. Test data are correlated with DYLOSAT predictions using both the NASA interim turbulence model and rotationally sampled winds as inputs. The latter is demonstrated to have the potential to improve the test/analysis correlations. The paper concludes with an assessment of the importance of vortex generators, site dependence, and machine differences on fatigue loads. The adequacy of prediction techniques used are evaluated and recommendations are made for improvements to the methodology.

Description: Empirical equations are presented with which to model rotationally-sampled (R-S) turbulence for input to structural-dynamic computer codes and the calculation of wind turbine fatigue loads. These equations are derived from R-S turbulence data which were measured at the vertical-plane array in Clayton, New Mexico. For validation, the equations are applied to the calculation of cyclic flapwise blade loads for the NASA/DOE Mod-2 2.5-MW experimental HAWT's (horizontal-axis wind turbines), and the results compared to measured cyclic loads. Good correlation is achieved, indicating that the R-S turbulence model developed in this study contains the characteristics of the wind which produce many of the fatigue loads sustained by wind turbines. Empirical factors are included which permit the prediction of load levels at specified percentiles of occurrence, which is required for the generation of fatigue load spectra and the prediction of the fatigue lifetime of structures.

Description: A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic axisymmetric response of viscoplastic shells of revolution. The analytical formulation is based on Reissner's large deformation shell theory with the effects of transverse shear deformation, rotatory inertia and moments turning around the normal to the middle surface included. The material model is chosen to be isothermal viscoplasticity, and an associated flow rule is used with a von Mises effective stress. A mixed formulation is used with the fundamental unknowns consisting of six stress resultants, three generalized displacements and three velocity components. Spatial discretization is performed using finite elements, with discontinuous stress resultants across element interfaces. The temporal integration is performed by using an explicit central difference scheme (leap-frog method) with an implicit constitutive update. The sensitivity coefficients are evaluated using a direct differentiation approach. Numerical results are presented for a spherical cap subjected to step loading, and a circular plate subjected to impulsive loading. The sensitivity coefficients are generated by evaluating the derivatives of the response quantities with respect to thickness, mass density, Young's modulus, and two of the material parameters characterizing the viscoplastic response. Time histories of the response and sensitivity coefficients are presented, along with spatial distributions of these quantities at selected times.

Description: Hypervelocity impact experiments were performed to further test the survivability of carbonaceous impactors and to determine potential products that may have been synthesized during impact. Diamonds were launched by the Ames two-stage light gas gun into Al plate at velocities of 2.75 and 3.1 km sec(exp -1). FESEM imagery confirms that diamond fragments survived in both experiments. Earlier experiments found that diamonds were destroyed on impact above 4.3 km sec(exp -1). Thus, the upper stability limit for diamond on impact into Al, as determined from our experimental conditions, is between 3.1 and 4.3 km sec(exp -1). Particles of the carbonaceous chondrite Nogoya were also launched into Al at a velocity of 6.2 km sec (exp -1). Laser desorption (L (exp 2) MS) analyses of the impactor residues indicate that the lowest and highest mass polycyclic aromatic hydrocarbons (PAH's) were largely destroyed on impact; those of intermediate mass (202-220 amu) remained at the same level or increased in abundance. In addition, alkyl-substituted homologs of the most abundant pre-impacted PAH's were synthesized during impact. These results suggest that an unknown fraction of some organic compounds can survive low to moderate impact velocities and that synthesized products can be expected to form up to velocities of, at least, 6.5 km sec(exp -1). We also present examples of craters formed by a unique microparticle accelerator that could launch micron-sized particles of almost any coherent material at velocities up to approximately 15 km sec(exp -1). Many of the experiments have a direct bearing on the interpretation of LDEF craters.

Description: Interpretation of the wealth of impact data available from the Long Duration Exposure Facility, in terms of the absolute and relative populations of space debris and natural micrometeoroids, requires three dimensional models of the distribution of impact directions, velocities and masses of such particles, as well as understanding of the impact processes. Although the stabilized orbit of LDEF provides limited directional information, it is possible to determine more accurate impact directions from detailed crater morphology. The applicability of this technique has already been demonstrated but the relationship between crater shape and impactor direction and velocity has not been derived in detail. We present the results of impact experiments and simulations: (1) impacts at micron dimensions using the Unit's 2MV Van de Graaff accelerator; (2) impacts at mm dimensions using a Light Gas Gun; and (3) computer simulations using AUTODYN-3D from which an empirical relationship between crater shape and impactor velocity, direction and particle properties we aim to derive. Such a relationship can be applied to any surface exposed to space debris or micrometeoroid particles for which a detailed pointing history is available.

Description: A 5.2 mm crater in Al-metal represents the largest found on LDEF. We have examined this crater by field emission scanning electron microscopy (FESEM), energy dispersive spectroscopy (EDS) and time-of-flight/secondary ion mass spectroscopy (TOF-SIMS) in order to determine if there is any evidence of impactor residue. Droplet and dome-shaped columns, along with flow features, are evidence of melting. EDS from the crater cavity and rim show Mg, C, O and variable amounts of Si, in addition to Al. No evidence for a chondritic impactor was found, and it hypothesized that the crater may be the result of impact with space debris.

Description: The finite deformation field of a plane strain Mode 1 crack in a hyperelastic and incompressible material was examined under the assumptions of small scale nonlinearity. Finite element analyses were performed for two different material laws, a Neo-Hookean material and a third order invariant of a Rivlin material. The numerical results for both materials were compared to the appropriate theoretical asymptotic solution. A local cavitation locus surrounding the crack tip was identified for the Neo-Hookean material. For the third order invariant Rivlin material, maximum values of the dominant stress component were found close to the surface of the crack, above and below the deformed crack tip.

Description: Approximately 20 sq m of protective thermal blankets, largely composed of Teflon, were retrieved from the Long Duration Exposure Facility (LDEF) after the spacecraft had spent approximately 5.7 years in space. Examination of these blankets revealed that they contained thousands of hypervelocity impact features ranging from micron-sized craters to penetration holes several millimeters in diameter. We conducted impact experiments in an effort to reproduce such features and to -- hopefully -- understand the relationships between projectile size and the resulting crater or penetration-hole diameter over a wide range of impact velocity. Such relationships are needed to derive the size- and mass-frequency distribution and flux of natural and man-made particles in low-Earth orbit. Powder propellant and light-gas guns were used to launch soda-lime glass spheres of 3.175 mm (1/8 inch) nominal diameter (Dp) into pure Teflon FEP targets at velocities ranging from 1 to 7 km/s. Target thickness (T) was varied over more than three orders of magnitude from infinite halfspace targets (Dp/T less than 0.1) to very thin films (Dp/T greater than 100). Cratering and penetration of massive Teflon targets is dominated by brittle failure and the development of extensive spall zones at the target's front and, if penetrated, the target's rear side. Mass removal by spallation at the back side of Teflon targets may be so severe that the absolute penetration-hole diameter (Dh) can become larger than that of a standard crater (Dc) at relative target thicknesses of Dp/T = 0.6-0.9. The crater diameter is infinite halfspace Teflon targets increases -- at otherwise constant impact conditions -- with encounter velocity by a factor of V0.44. In contrast, the penetration-hole size is very thin foils (Dp/T greater than 50) is essentially unaffected by impact velocity. Penetrations at target thicknesses intermediate to these extremes will scale with variable exponents of V. Our experimental matrix is sufficiently systematic and complete, up to 7 km/s, to make reasonable recommendations for the velocity-scaling of Teflon craters and penetrations. We specifically suggest that cratering behavior and associated equations dominate all impacts in which the shock-pulse duration of the projectile (tp) is shorter than that of the target (tt). We also demonstrate that each penetration hole from space-retrieved surfaces may be assigned a unique projectile size, provided an impact velocity is known or assumed. This calibration seems superior to the traditional ballistic-limit approach.

Description: A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic frictional contact/impact response of axisymmetric composite structures. The structures are assumed to consist of an arbitrary number of perfectly bonded homogeneous anisotropic layers. The material of each layer is assumed to be hyperelastic, and the effect of geometric nonlinearity is included. The sensitivity coefficients measure the sensitivity of the response to variations in different material, lamination and geometric parameters of the structure. A displacement finite element model is used for the discretization. The normal contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the nodal displacements, and the Lagrange multipliers associated with the contact conditions. The Lagrange multipliers are allowed to be discontinuous at interelement boundaries. Tangential contact conditions are incorporated by using a penalty method in conjunction with the classical Coulomb's friction model. Temporal integration is performed by using Newmark method. The Newton-Raphson iterative scheme is used for the solution of the resulting nonlinear algebraic equations, and for the determination of the contact region, contact conditions (sliding or sticking), and the contact pressures. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented from the frictional contact/impact response of a composite spherical cap impacting on a rigid plate.

Description: A computational strategy is presented for postbuckling and nonlinear static analyses of large complex structures on distributed-memory parallel computers. The strategy is designed for message-passing parallel computer systems. The key elements of the proposed strategy are: (1) a multiple-parameter reduced basis technique; (2) a nested dissection (or multilevel substructuring) ordering scheme; (3) parallel assembly of global matrices; and (4) a parallel sparse equation solver. The effectiveness of the strategy is assessed by performing thermomechanical postbuckling analyses of stiffened composite panels with cutouts, and nonlinear large-deflection analyses of High Speed Civil Transport models on three distributed-memory computers. The numerical studies presented demonstrate the advantages of nested dissection-based solvers over traditional skyline-based solvers on distributed-memory machines.

Description: A computational strategy is presented for the nonlinear static and postbuckling analyses of large complex structures on massively parallel computers. The strategy is designed for distributed-memory, message-passing parallel computer systems. The key elements of the proposed strategy are: (1) a multiple-parameter reduced basis technique; (2) a nested dissection (or multilevel substructuring) ordering scheme; (3) parallel assembly of global matrices; and (4) a parallel sparse equation solver. The effectiveness of the strategy is assessed by applying it to thermo-mechanical postbuckling analyses of stiffened composite panels with cutouts, and nonlinear large-deflection analyses of HSCT models on Intel Paragon XP/S computers. The numerical studies presented demonstrate the advantages of nested dissection-based solvers over traditional skyline-based solvers on distributed memory machines.

Description: A reduced basis technique and a computational procedure are presented for generating the nonlinear vibrational response, and evaluating the first-order sensitivity coefficients of thin-walled composite frames. The sensitivity coefficients are the derivatives of the nonlinear frequency with respect to the material and lamination parameters of the frame. A mixed formulation is used with the fundamental unknowns consisting of both the generalized displacements and stress resultants in the frame. The flanges and webs of the frames are modeled by using geometrically nonlinear two-dimensional shell and plate finite elements. The computational procedure can be conveniently divided into three distinct steps. The first step involves the generation of various-order perturbation vectors, and their derivatives with respect to the material and lamination parameters of the frame, using the Linstedt-Poincare perturbation technique. The second step consists of using the perturbation vectors as basis vectors, computing the amplitudes of these vectors and the nonlinear frequency of vibration, via a direct variational procedure. The third step consists of using the perturbation vectors, and their derivatives, as basis vectors and computing the sensitivity coefficients of the nonlinear frequency via a second application of the direct variational procedure. Numerical results are presented for semicircular thin-walled frames with I and J section, showing the convergence of the nonlinear frequency and the sensitivity coefficients obtained by both the reduced-basis and perturbation techniques.

Description: A computational strategy is presented for calculating sensitivity coefficients for the nonlinear large-deflection and postbuckling responses of laminated composite structures on distributed-memory parallel computers. The strategy is applicable to any message-passing distributed computational environment. The key elements of the proposed strategy are: (1) a multiple-parameter reduced basis technique; (2) a parallel sparse equation solver based on a nested dissection (or multilevel substructuring) node ordering scheme; and (3) a multilevel parallel procedure for evaluating hierarchical sensitivity coefficients. The hierarchical sensitivity coefficients measure the sensitivity of the composite structure response to variations in three sets of interrelated parameters; namely, laminate, layer and micromechanical (fiber, matrix, and interface/interphase) parameters. The effectiveness of the strategy is assessed by performing hierarchical sensitivity analysis for the large-deflection and postbuckling responses of stiffened composite panels with cutouts on three distributed-memory computers. The panels are subjected to combined mechanical and thermal loads. The numerical studies presented demonstrate the advantages of the reduced basis technique for hierarchical sensitivity analysis on distributed-memory machines.

Description: A reduction technique and a computational procedure are presented for predicting the tire contact response and evaluating the sensitivity coefficients of the different response quantities. The sensitivity coefficients measure the sensitivity of the contact response to variations in the geometric and material parameters of the tire. The tire is modeled using a two-dimensional laminated anisotropic shell theory with the effects of variation in geometric and material parameters, transverse shear deformation, and geometric nonlinearities included. The contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the stress resultants, the generalized displacements, and the Lagrange multipliers associated with the contact conditions. The elemental arrays are obtained by using a modified two-field, mixed variational principle. For the application of the reduction technique, the tire finite element model is partitioned into two regions. The first region consists of the nodes that are likely to come in contact with the pavement, and the second region includes all the remaining nodes. The reduction technique is used to significantly reduce the degrees of freedom in the second region. The effectiveness of the computational procedure is demonstrated by a numerical example of the frictionless contact response of the space shuttle nose-gear tire, inflated and pressed against a rigid flat surface.

Description: Analytic three-dimensional solutions are presented for the coupled thermoelectroelastic response of multilayered hybrid composite plates. The plates consist of a combination of fiber-reinforced cross-ply and piezothermoelastic layers. Both the thermoelectroelastic static response and its sensitivity coefficients are computed. The sensitivity coefficients measure the sensitivity of the response to variations in different mechanical, thermal and piezoelectric material properties of the plate. A linear constitutive model is used, and the material properties are assumed to be independent of the temperature and the electric field. The plates are assumed to have rectangular geometry and special material symmetries. A mixed formulation is used with the fundamental unknowns consisting of the three transverse stress components; three displacement components; transverse component of the electric displacement field; electric potential; transverse heat flux component, and temperature change. Each of the fundamental unknowns is expressed in terms of a double Fourier series in the Cartesian surface coordinates. A state space approach is used to generate the static response and to evaluate the sensitivity coefficients. Extensive numerical results are presented showing the effects of variation in the geometric parameters of the plate on the different response quantities and their sensitivity coefficients.

Description: A computational procedure is presented for evaluating the sensitivity coefficients of the static frictional contact response of axisymmetric composite structures. The structures are assumed to consist of an arbitrary number of perfectly bonded homogeneous anisotropic layers. The material of each layer is assumed to be hyperelastic, and the effect of geometric nonlinearity is included. The sensitivity coefficients measure the sensitivity of the response variations in different material, lamination and geometric parameters of the structure. A displacement finite element model is used for the discretization. The normal contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of nodal displacements, and Lagrange multipliers associated with the contact conditions. The Lagrange multipliers are allowed to be discontinuous at interelement boundaries. Tangential contact conditions are incorporated by using a penalty method in conjunction with the classical Coulomb's friction model. The Newton-Raphson iterative scheme is used for the solution of the resulting nonlinear algebraic equations, and for the determination of the contact region, contact conditions (sliding or sticking), and the contact pressures. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented for the frictional contact of a composite spherical cap pressed against a rigid plate.

Description: The results of a study of the detailed thermomechanical postbuckling response characteristics of flat unstiffened composite panels with central circular cutouts are presented. The panels are subjected to combined temperature changes and applied edge loading (or edge displacements). The analysis is based on a first-order shear deformation plate theory. A mixed formulation is used with the fundamental unknowns consisting of the generalized displacements and the stress resultants of the plate. The postbuckling displacements, transverse shear stresses, transverse shear strain energy density, and their sensitivity coefficients are evaluated. The sensitivity coefficients measure the sensitivity of the post-buckling response to variations in the different lamination and material parameters of the panel. Numerical results are presented showing the effects of the variations in the hole diameter, laminate stacking sequence, fiber orientation, and aspect ratio of the panel on the thermomechanical postbuckling response and its sensitivity to changes in panel parameters.

Description: A study is made of the effect of mesh distortion on the accuracy of transverse shear stresses and their first-order and second-order sensitivity coefficients in multilayered composite panels subjected to mechanical and thermal loads. The panels are discretized by using a two-field degenerate solid element, with the fundamental unknowns consisting of both displacement and strain components, and the displacement components having a linear variation throughout the thickness of the laminate. A two-step computational procedure is used for evaluating the transverse shear stresses. In the first step, the in-plane stresses in the different layers are calculated at the numerical quadrature points for each element. In the second step, the transverse shear stresses are evaluated by using piecewise integration, in the thickness direction, of the three-dimensional equilibrium equations. The same procedure is used for evaluating the sensitivity coefficients of transverse shear stresses. Numerical results are presented showing no noticeable degradation in the accuracy of the in-plane stresses and their sensitivity coefficients with mesh distortion. However, such degradation is observed for the transverse shear stresses and their sensitivity coefficients. The standard of comparison is taken to be the exact solution of the three-dimensional thermoelasticity equations of the panel.

Description: A study is made of the buckling and postbuckling responses of flat, unstiffened composite panels subjected to various combinations of mechanical and thermal loads. The analysis is based on a first-order shear deformation von Karman-type plate theory. A mixed formulation is used with the fundamental unknowns consisting of the strain components, stress resultants and the generalized displacements of the plate. The stability boundary, postbuckling response and the sensitivity coefficients are evaluated. The sensitivity coefficients measure the sensitivity of the buckling and postbuckling responses to variations in the different lamination and material parameters of the panel. Numerical results are presented for both solid panels and panels with central circular cutouts. The results show the effects of the variations in the fiber orientation angels, aspect ratio of the panel, and the hole diameter (for panels with cutouts) on the stability boundary, postbuckling response and sensitivity coefficients.

Description: A study is made of the accuracy of the steady-state (static) thermoelectroelastic response of multilayered hybrid composite plates predicted by five modeling approaches, based on two-dimensional plate theories. The plates consist of a combination of fiber-reinforced and piezothermoelastic layers. The standard of comparison is taken to be the exact three-dimensional thermoelectroelastic solutions, and the quantities compared include gross response characteristics (e.g., strain energy components, and average through-the-thickness displacements); detailed, through-the-thickness distributions of displacements and stresses; and sensitivity coefficients of the response quantities (derivatives of the response quantities with respect to material parameters of the plate). The modeling approaches considered include first-order theory; third-order theory; discrete-layer theory (with piecewise linear variation of the in-plane displacements, temperature and electric potential, in the thickness direction); and two predictor-corrector procedures.

Description: A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic axisymmetric, fully-coupled, thermoviscoplastic response of shells of revolution. The analytical formulation is based on Reissner's large deformation shell theory with the effects of large-strain, transverse shear deformation, rotatory inertia and moments turning around the normal to the middle surface included. The material model is chosen to be viscoplasticity with strain hardening and thermal hardening, and an associated flow rule is used with a von Mises effective stress. A mixed formulation is used for the shell equations with the fundamental unknowns consisting of six stress resultants, three generalized displacements and three velocity components. The energy-balance equation is solved using a Galerkin procedure, with the temperature as the fundamental unknown. Spatial discretization is performed in one dimension (meridional direction) for the momentum and constitutive equations of the shell, and in two dimensions (meridional and thickness directions) for the energy-balance equation. The temporal integration is performed by using an explicit central difference scheme (leap-frog method) for the momentum equation; a predictor-corrector version of the trapezoidal rule is used for the energy-balance equation; and an explicit scheme consistent with the central difference method is used to integrate the constitutive equations. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented for a spherical cap subjected to step loading.

Description: The results of a detailed study of the buckling and postbuckling responses of composite panels with central circular cutouts are presented. The panels are subjected to combined edge shear and temperature change. The panels are discretized by using a two-field degenerate solid element with each of the displacement components having a linear variation throughout the thickness of the panel. The fundamental unknowns consist of the average mechanical strains through the thickness and the displacement components. The effects of geometric nonlinearities and laminated anisotropic material behavior are included. The stability boundary, postbuckling response and the hierarchical sensitivity coefficients are evaluated. The hierarchical sensitivity coefficients measure the sensitivity of the buckling and postbuckling responses to variations in the panel stiffnesses, and the material properties of both the individual layers and the constituents (fibers and matrix). Numerical results are presented for composite panels with central circular cutouts subjected to combined edge shear and temperature change, showing the effects of variations in the hole diameter, laminate stacking sequence and fiber orientation, on the stability boundary and postbuckling response and their sensitivity to changes in the various panel parameters.

Description: A hybrid neurocomputing/numerical strategy is presented for geometrically nonlinear analysis of structures. The strategy combines model-free data processing capabilities of computational neural networks with a Pade approximants-based perturbation technique to predict partial information about the nonlinear response of structures. In the hybrid strategy, multilayer feedforward neural networks are used to extend the validity of solutions by using training samples produced by Pade approximations to the Taylor series expansion of the response function. The range of validity of the training samples is taken to be the radius of convergence of Pade approximants and is estimated by setting a tolerance on the diverging approximants. The norm of residual vector of unbalanced forces in a given element is used as a measure to assess the quality of network predictions. To further increase the accuracy and the range of network predictions, additional training data are generated by either applying linear regression to weight matrices or expanding the training data by using predicted coefficients in a Taylor series. The effectiveness of the hybrid strategy is assessed by performing large-deflection analysis of a doubly-curved composite panel with a circular cutout, and postbuckling analyses of stiffened composite panels subjected to an in-plane edge shear load. In all the problems considered, the hybrid strategy is used to predict selective information about the structural response, namely the total strain energy and the maximum displacement components only.

Description: Three recent developments in the sensitivity analysis for thermomechanical postbuckling response of composite panels are reviewed. The three developments are: (1) effective computational procedure for evaluating hierarchical sensitivity coefficients of the various response quantities with respect to the different laminate, layer, and micromechanical characteristics; (2) application of reduction methods to the sensitivity analysis of the postbuckling response; and (3) accurate evaluation of the sensitivity coefficients to transverse shear stresses. Sample numerical results are presented to demonstrate the effectiveness of the computational procedures presented. Some of the future directions for research on sensitivity analysis for the thermomechanical postbuckling response of composite and smart structures are outlined.

Description: A computational procedure is presented for evaluating the sensitivity coefficients of the viscoplastic response of structures subjected to dynamic loading. A state of plane stress is assumed to exist in the structure, a velocity strain-Cauchy stress formulation is used, and the geometric non-linearities arising from large strains are incorporated. The Jaumann rate is used as a frame indifferent stress rate. The material model is chosen to be isothermal viscoplasticity, and an associated flow rule is used with a von Mises effective stress. The equations of motion emanating from a finite element semi-discretization are integrated using an explicit central difference scheme with an implicit stress update. The sensitivity coefficients are evaluated using a direct differentiation approach. Since the domain of integration is the current configuration, the sensitivity coefficients of the spatial derivatives of the shape functions must be included. Numerical results are presented for a thin plate with a central cutout subjected to an in-plane compressive loading. The sensitivity coefficients are generated by evaluating the derivatives of the response quantities with respect to Young's modulus, and two of the material parameters characterizing the viscoplastic response. Time histories of the response and sensitivity coefficients, and spatial distributions at selected times are presented.

Description: Analytic three-dimensional thermoelasticity solutions are presented for static problems of simply supported sandwich panels and cylindrical shells subjected to mechanical and thermal loads. The panels and shells have laminated composite face sheets of arbitrary thickness separated by a core. Each of the individual layers of the face sheets and the core is modeled as a three-dimensional continuum. Analytic first-order sensitivity coefficients are evaluated to assess the sensitivity of the responses to variations in material parameters of the face sheets and the core, as well as to variations in the curvatures and thicknesses of the sandwich and face sheets. Also, the strain energy associated with various stress components in the face sheets and core are calculated and compared. The information obtained in the present study can aid the development and assessment of two-dimensional models for sandwich structures and illuminate the role of particular material parameters in an equivalent model for the core.

Description: Three stiffened panel concepts are evaluated to find optimized designs for integral stiffeners in the barrels of Reusable Launch Vehicle fuel tanks. The three panel concepts considered are a T-stiffened panel, a panel with one blade stiffener centered between each pair of T-stiffeners, and a panel with two blade stiffeners equally spaced between each pair of T-stiffeners. The panels are optimized using PASCO for a range of compressive loads, and the computed areal weight for each panel is used to compare the concepts and predict tank weights. The areal weight of the T-stiffened panel with one blade is up to seven-percent lower than the other panel concepts. Two tank construction methods are compared for a representative tank design with three barrels. In the first method, 45-degree circumferential sections of a barrel are each designed to carry the same maximum load in the barrel. In the second method, each barrel section is designed for the maximum load in that section. Representative tanks designed with the first method are over 250 lb heavier than tanks designed using the second method. Optimized panel designs and areal weights are also computed for variation of the nominal panel length and skin thickness.

Description: A parametric study is presented of the buckling behavior of infinitely long symmetrically laminated anisotropic plates subjected to combined loads. The study focuses on the interaction of a subcritical (stable) secondary loading state of constant magnitude and a primary destabilizing load that is increased in magnitude until buckling occurs. The loads, considered in this report are uniform axial compression, pure in-plane bending, transverse tension and compression, and shear. Results are presented that were obtained by using a special purpose nondimensional analysis that is well suited for parametric studies of clamped and simply supported plates. In particular, results are presented for a +/- 45(sub S) graphite-epoxy laminate that is highly anisotropic and representative of a laminate used for spacecraft applications. In addition, generic buckling-design charts are presented for a wide range of nondimensional parameters that are applicable to a broad class of laminate constructions. These results show the general behavioral trends of specially orthotropic plates and the effects of flexural anisotropy on plates subjected to various combined loading conditions. An important finding of the present study is that the effects of flexural anisotropy on the buckling resistance of a plate can be significantly more important for plates subjected to combined loads than for plates subjected to single-component loads.

Description: This document contains presentations and discussions from the joint UVA/NASA Workshop on Computational Modeling of Tires. The workshop attendees represented NASA, the Army and Air force, tire companies, commercial software developers, and academia. The workshop objectives were to assess the state of technology in the computational modeling of tires and to provide guidelines for future research.

Description: Three critical compression splice joint locations in a stitched graphite-epoxy transport wing stub box have been analyzed to determine their expected structural performance. The wing box is representative of a section of a commercial transport wing box and was designed and constructed by McDonnell Douglas Aerospace Company as part of the NASA Advanced Composites Technology (ACT) program. The results of the finite element analyses of the splice joints are presented. The analysis results indicate that failure will not occur in the splice joint regions for loads less than the Design Ultimate Load of the wing box.

Description: Thermal buckling characteristics of hypersonic aircraft sandwich panels of various aspect ratios were investigated. The panel is fastened at its four edges to the substructures under four different edge conditions and is subjected to uniform temperature loading. Minimum potential energy theory and finite element methods were used to calculate the panel buckling temperatures. The two methods gave fairly close buckling temperatures. However, the finite element method gave slightly lower buckling temperatures than those given by the minimum potential energy theory. The reasons for this slight discrepancy in eigensolutions are discussed in detail. In addition, the effect of eigenshifting on the eigenvalue convergence rate is discussed.

Description: An investigation was performed to ascertain the feasibility of using IPACS (Integrated Probabilistic Assessment of Composite Structures) for probabilistic analysis of a composite fan blade, the development of which is being pursued by various industries for the next generation of aircraft engines. A model representative of the class of fan blades used in the GE90 engine has been chosen as the structural component to be analyzed with IPACS. In this study, typical uncertainties are assumed in the level, and structural responses for ply stresses and frequencies are evaluated in the form of cumulative probability density functions. Because of the geometric complexity of the blade, the number of plies varies from several hundred at the root to about a hundred at the tip. This represents a extremely complex composites application for the IPACS code. A sensitivity study with respect to various random variables is also performed.

Description: The article begins by examining the fundamentals of traditional deterministic design philosophy. The initial section outlines the concepts of failure criteria and limit state functions two traditional notions that are embedded in deterministic design philosophy. This is followed by a discussion regarding safety factors (a possible limit state function) and the common utilization of statistical concepts in deterministic engineering design approaches. Next the fundamental aspects of a probabilistic failure analysis are explored and it is shown that deterministic design concepts mentioned in the initial portion of the article are embedded in probabilistic design methods. For components fabricated from ceramic materials (and other similarly brittle materials) the probabilistic design approach yields the widely used Weibull analysis after suitable assumptions are incorporated. The authors point out that Weibull analysis provides the rare instance where closed form solutions are available for a probabilistic failure analysis. Since numerical methods are usually required to evaluate component reliabilities, a section on Monte Carlo methods is included to introduce the concept. The article concludes with a presentation of the technical aspects that support the numerical method known as fast probability integration (FPI). This includes a discussion of the Hasofer-Lind and Rackwitz-Fiessler approximations.

Description: An effort is currently underway to develop an experiment titled joint Damping E_periment (JDX) to fly on the Space Shuttle as Get Away Special Payload G-726. This project is funded by NASA's IN-Space Technology Experiments Program and is scheduled to fly in July 1995 on STS-69. JDX will measure the influence of gravity on the structural damping of a three bay truss having clearance fit pinned joints. Structural damping is an important parameter in the dynamics of space structures. Future space structures will require more precise knowledge of structural damping than is currently available. The mission objectives are to develop a small-scale shuttle flight experiment that allows researchers to: (1) characterize the influence of gravity and joint gaps on structural damping and dynamic behavior of a small-scale truss model, and (2) evaluate the applicability of low-g aircraft test results for predicting on-orbit behavior. Completing the above objectives will allow a better understanding and/or prediction of structural damping occurring in a pin jointed truss. Predicting damping in joints is quite difficult. One of the important variables influencing joint damping is gravity. Previous work has shown that gravity loads can influence damping in a pin jointed truss structure. Flying this experiment as a GAS payload will allow testing in a microgravity environment. The on-orbit data (in micro-gravity) will be compared with ground test results. These data will be used to help develop improved models to predict damping due to pinned joints. Ground and low-g aircraft testing of this experiment has been completed. This paper describes the experiment and presents results of both ground and low-g aircraft tests which demonstrate that damping of the truss is dramatically influenced by gravity.

Description: NASA is developing a 'tool box' that includes a number of advanced structural analysis computer codes which, taken together, represent the comprehensive fracture mechanics capability required to predict the onset of widespread fatigue damage. These structural analysis tools have complementary and specialized capabilities ranging from a finite-element-based stress-analysis code for two- and three-dimensional built-up structures with cracks to a fatigue and fracture analysis code that uses stress-intensity factors and material-property data found in 'look-up' tables or from equations. NASA is conducting critical experiments necessary to verify the predictive capabilities of the codes, and these tests represent a first step in the technology-validation and industry-acceptance processes. NASA has established cooperative programs with aircraft manufacturers to facilitate the comprehensive transfer of this technology by making these advanced structural analysis codes available to industry.

Description: Bending strength, fracture toughness, and elastic modulus data were acquired for YBa2Cu3O(x), Bi2Sr2CaCu2O(x) (Bi,Pb)2Sr2Ca2Cu3O(x), and Tl2Ba2Ca2Cu3O(x) bars. These data and thermal expansion coefficients strongly suggest that the maximum possible tensile strain without fracture of bulk tapes or wires is approximately equals 0.2%. In Ag-clad conductors, residual stresses will be of limited benefit, but fractures produced by larger strains can be accommodated by shunting current through the Ag.

Description: Aircraft tire wear results from operating conditions that are quite different from those encountered in land vehicles. One of the most important of these is touchdown, where the tire suddenly spins up from zero to a large angular velocity. This phenomenon is studied from both the analytical and experimental points of view. The analysis is basic, using elementary properties of the tire and wheel. It results in a new dimensionless description of the process. The experimental study consists primarily of small scale laboratory data, although limited full scale tire data is also presented. The data show increasing weight loss during touchdown as the dimensionless severity increases.

Description: The present invention related generally to severing materials, and more particularly to severing or weakening materials through explosively induced, augmented shock waves. Explosive cords are placed in grooves on the upper surface of the material to be severed or weakened. The explosive cords are initiated simultaneously to introduce explosive shock waves into the material. These shock waves progress toward the centerline between the explosive cords and the lower surface of the material. Intersecting and reflected waves produce a rarefaction zone on the centerline to fail the material in tension. A groove may also be cut in the lower surface of the material to aid in severing or weakening the material.

Description: A FORTRAN program for calculating the stresses of n cracks embedded in an isotropic plate is presented. Formulas are given for calculating the stresses of one crack, two cracks, and n cracks in an isotropic plate. Then the program code that accomplishes this is provided.

Description: The objects of the first, exploratory, stage of the project were listed as: (1) to make a detailed and critical review of the Boundary Element method as already published and with regard to elastic-plastic fracture mechanics, to assess its potential for handling present concepts in two-dimensional and three-dimensional cases. To this was subsequently added the Finite Volume method and certain aspects of the Finite Element method for comparative purposes; (2) to assess the further steps needed to apply the methods so far developed to the general field, covering a practical range of geometries, work hardening materials, and composites: to consider their application under higher temperature conditions; (3) to re-assess the present stage of development of the energy dissipation rate, crack tip opening angle and J-integral models in relation to the possibilities of producing a unified technology with the previous two items; and (4) to report on the feasibility and promise of this combined approach and, if appropriate, make recommendations for the second stage aimed at developing a generalized crack growth technology for its application to real-life problems.

Description: An experimental program was conducted to study the damaging effects of tensile and compressive prestrains on the fatigue life of nickel-base, Inconel 718 superalloy at room temperature. To establish baseline fatigue behavior, virgin specimens with a solid uniform gage section were fatigued to failure under fully-reversed strain-control. Additional specimens were prestrained to 2 percent, 5 percent, and 10 percent (engineering strains) in the tensile direction and to 2 percent (engineering strain) in the compressive direction under stroke-control, and were subsequently fatigued to failure under fully-reversed strain-control. Experimental results are compared with estimates of remaining fatigue lives (after prestraining) using three life prediction approaches: (1) the Linear Damage Rule; (2) the Linear Strain and Life Fraction Rule; and (3) the nonlinear Damage Curve Approach. The Smith-Watson-Topper parameter was used to estimate fatigue lives in the presence of mean stresses. Among the cumulative damage rules investigated, best remaining fatigue life predictions were obtained with the nonlinear Damage Curve Approach.

Description: Design loads are presented for the General Electric MOD-SA wind turbine. The MOD-SA system consists of a 400 ft. diameter, upwind, two-bladed, teetered rotor connected to a 7.3 mW variable-speed generator. Fatigue loads are specified in the form of histograms for the 30 year life of the machine, while limit (or maximum) loads have been derived from transient dynamic analysis at critical operating conditions. Loads prediction was accomplished using state of the art aeroelastic analyses developed at General Electric. Features of the primary predictive tool - the Transient Rotor Analysis Code (TRAC) are described in the paper. Key to the load predictions are the following wind models: (1) yearly mean wind distribution; (2) mean wind variations during operation; (3) number of start/shutdown cycles; (4) spatially large gusts; and (5) spatially small gusts (local turbulence). The methods used to develop statistical distributions from load calculations represent an extension of procedures used in past wind programs and are believed to be a significant contribution to Wind Turbine Generator analysis. Test/theory correlations are presented to demonstrate code load predictive capability and to support the wind models used in the analysis. In addition MOD-5A loads are compared with those of existing machines. The MOD-5A design was performed by the General Electric Company, Advanced Energy Program Department, under Contract DEN3-153 with NASA Lewis Research Center and sponsored by the Department of Energy.

Description: The implementation of an interface element capability within the COMET-AR software system is described. The report is intended for use by both users of currently implemented interface elements and developers of new interface element formulations. Guidance on the use of COMET-AR is given. A glossary is provided as an Appendix to this report for readers unfamiliar with the jargon of COMET-AR. A summary of the currently implemented interface element formulation is presented in Section 7.3 of this report.

Description: Impact dampers belong to the category of passive vibration devices used to attenuate the vibration of discrete and continuous systems. An impact damper generally consists of a mass which is allowed to travel freely between two defined stops. Under the right conditions, the vibration of the structure to which the impact damper is attached will cause the mass of the impact damper to strike the structure. Previous analytical and experimental research work on the effect of impact dampers in attenuating the vibration of discrete and continuous systems have demonstrated their effectiveness. It has been shown in this study that impact dampers can increase the intrinsic damping of a lightly-damped flexible structure. The test structure consists of a slender flexible beam supported by a pin-type support at one end and supported by a linear helical flexible spring at another location. Sinusoidal excitation spanning the first three natural frequencies was applied in the horizontal plane. The orientation of the excitation and the test structure in the horizontal plane minimizes the effect of gravity on the behavior of the test structure. The excitation was applied using a linear sine sweep technique. The span of the test structure, the mass of the impact damper, the distance of travel, and the location of the impact damper along the span of the test structure were varied. The damping ratio are estimated for sixty test configurations. The results show that the impact damper significantly increases the damping ratio of the test structure. Statistical analysis of the results using the method of multiple linear regression indicates that a reasonable fit has been accomplished. It is concluded that additional experimental analysis of flexible structures in microgravity environment is needed in order to achieve a better understanding of the behavior of impact damper under conditions of microgravity. Numerical solution of the behavior of flexible structures equipped with impact dampers is also needed to predict stresses and deformations under operating conditions of microgravity in space applications.

Description: The FPCAS3D computer code has been developed for aeroelastic stability analysis of bladed disks such as those in fans, compressors, turbines, propellers, or propfans. The aerodynamic analysis used in this code is based on the unsteady three-dimensional full potential equation which is solved for a blade row. The structural analysis is based on a finite-element model for each blade. Detailed explanations of the aerodynamic analysis, the numerical algorithms, and the aeroelastic analysis are not given in this report. This guide can be used to assist in the preparation of the input data required by the FPCAS3D code. A complete description of the input data is provided in this report. In addition, six examples, including inputs and outputs, are provided.

Description: A probabilistic design method is described and demonstrated using a smart composite wing. Probabilistic structural design incorporates naturally occurring uncertainties including those in constituent (fiber/matrix) material properties, fabrication variables, structure geometry and control-related parameters. Probabilistic sensitivity factors are computed to identify those parameters that have a great influence on a specific structural reliability. Two performance criteria are used to demonstrate this design methodology. The first criterion requires that the actuated angle at the wing tip be bounded by upper and lower limits at a specified reliability. The second criterion requires that the probability of ply damage due to random impact load be smaller than an assigned value. When the relationship between reliability improvement and the sensitivity factors is assessed, the results show that a reduction in the scatter of the random variable with the largest sensitivity factor (absolute value) provides the lowest failure probability. An increase in the mean of the random variable with a negative sensitivity factor will reduce the failure probability. Therefore, the design can be improved by controlling or selecting distribution parameters associated with random variables. This can be implemented during the manufacturing process to obtain maximum benefit with minimum alterations.

Description: The Generalized Method of Cells (GMC), a micromechanics based constitutive model, is implemented into the finite element code MARC using the user subroutine HYPELA. Comparisons in terms of transverse deformation response, micro stress and strain distributions, and required CPU time are presented for GMC and finite element models of fiber/matrix unit cell. GMC is shown to provide comparable predictions of the composite behavior and requires significantly less CPU time as compared to a finite element analysis of the unit cell. Details as to the organization of the HYPELA code are provided with the actual HYPELA code included in the appendix.

Description: Because the 2195 aluminum-lithium material of the super lightweight external tank (SLWT ET) has a lower toughness than the 2219 aluminum used in previous ET's, careful attention must be paid to stress concentrations. This report details the analysis performed on some of the stress concentrations in the orthogrid panels of the liquid hydrogen tank.

Description: Passive isolator, active vibration absorber, and an integrated passive/active (hybrid) control are studied for their effectiveness in reducing structural vibration under seismic excitations. For the passive isolator, a laminated rubber bearing base isolator which has been studied and used extensively by researchers and seismic designers is considered. An active vibration absorber concept, which can provide guaranteed closed-loop stability with minimum knowledge of the controlled system, is used to reduce the passive isolator displacement and to suppress the top floor vibration. A three-story building model is used for the numerical simulation. The performance of an active vibration absorber and a hybrid vibration controller in reducing peak structural responses is compared with the passively isolated structural response and with absence of vibration control systems under the N00W component of El Centro 1940 and N90W component of the Mexico City earthquake excitation records. The results show that the integrated passive/active vibration control system is most effective in suppressing the peak structural acceleration for the El Centro 1940 earthquake when compared with the passive or active vibration absorber alone. The active vibration absorber, however, is the only system that suppresses the peak acceleration of the structure for the Mexico City 1985 earthquake.

Description: A high-temperature optical window is essential to the optical diagnostics of high-temperature combustion rigs. Laser Doppler velocimetry, schlieren photography, light sheet visualization, and laser-induced fluorescence spectroscopy are a few of the tests that require optically clear access to the combustor flow stream. A design was developed for a high-temperature window that could withstand the severe environment of the NASA Lewis 3200 F Lean Premixed Prevaporized (LPP) Flame Tube Test Rig. The development of this design was both time consuming and costly. This report documents the design process and the lessons learned, in an effort to reduce the cost of developing future designs for high-temperature optical windows.

Description: The Space Power Institute (SPI) at Auburn University has conducted experiments on the effects of impact angle on crater morphology and impactor residue retention for hypervelocity impacts. Copper target plates were set at angles of 30 deg, 45 deg, 60 deg, and 75 deg from the particle flight path. For the 30 deg and 45 deg impacts, in the velocity regime greater than 8 km s(exp -1) the resultant craters are almost identical to normal incidence impacts. The only difference found was in the apparent distribution of particle residue within the crater, and further research is needed to verify this. The 60 deg and 75 deg impacts showed marked differences in crater symmetry, crater lip shape, and particle residue distribution in the same velocity regime. Impactor residue shock fractionation effects have been quantified in first-order. It is concluded that a combination of analysis techniques can yield further information on impact velocity, direction, and angle of incidence.

Description: The Center for Computational Structures Technology (CST) is intended to serve as a focal point for the diverse CST research activities. The CST activities include the use of numerical simulation and artificial intelligence methods in modeling, analysis, sensitivity studies, and optimization of flight-vehicle structures. The Center is located at NASA Langley and is an integral part of the School of Engineering and Applied Science of the University of Virginia. The key elements of the Center are: (1) conducting innovative research on advanced topics of CST; (2) acting as pathfinder by demonstrating to the research community what can be done (high-potential, high-risk research); (3) strong collaboration with NASA scientists and researchers from universities and other government laboratories; and (4) rapid dissemination of CST to industry, through integration of industrial personnel into the ongoing research efforts.

Description: This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present first part of the report, we focus on the theoretical developments, and discussions of the results of numerical-performance studies using the integration schemes for GVIPS and NAV models.

Description: This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present second part of the report, we focus on the specific details of the numerical schemes, and associated computer algorithms, for the finite-element implementation of GVIPS and NAV models.

Description: A new method for commanding machines to move with increased dynamic performance was developed. This method is an enhanced version of input shaping, a patented vibration suppression algorithm. This technique intercepts a command input to a system command that moves the mechanical system with increased performance and reduced residual vibration. This document describes many advanced methods for generating highly optimized shaping sequences which are tuned to particular systems. The shaping sequence is important because it determines the trade off between move/settle time of the system and the insensitivity of the input shaping algorithm to variations or uncertainties in the machine which can be controlled. For example, a system with a 5 Hz resonance that takes 1 second to settle can be improved to settle instantaneously using a 0.2 shaping sequence (thus improving settle time by a factor of 5). This system could vary by plus or minus 15% in its natural frequency and still have no apparent vibration. However, the same system shaped with a 0.3 second shaping sequence could tolerate plus or minus 40% or more variation in natural frequency. This document describes how to generate sequences that maximize performance, sequences that maximize insensitivity, and sequences that trade off between the two. Several software tools are documented and included.

Description: Solutions to increasingly larger structural optimization problems are desired. However, computational resources are strained to meet this need. New methods will be required to solve increasingly larger problems. The present approaches to solving large-scale problems involve approximations for the constraints of structural optimization problems and/or decomposition of the problem into multiple subproblems that can be solved in parallel. An area of game theory, equilibrium programming (also known as noncooperative game theory), can be used to unify these existing approaches from a theoretical point of view (considering the existence and optimality of solutions), and be used as a framework for the development of new methods for solving large-scale optimization problems. Equilibrium programming theory is described, and existing design techniques such as fully stressed design and constraint approximations are shown to fit within its framework. Two new structural design formulations are also derived. The first new formulation is another approximation technique which is a general updating scheme for the sensitivity derivatives of design constraints. The second new formulation uses a substructure-based decomposition of the structure for analysis and sensitivity calculations. Significant computational benefits of the new formulations compared with a conventional method are demonstrated.

Description: Because the 2195 aluminum-lithium of the super lightweight external tank (SLWT ET) has a lower toughness than the 2219 aluminum used in previous ET's, careful attention must be paid to stress concentration in the SLWT ET. This report details the initial analysis performed by NASA to determine the material properties required to ensure structural integrity in these critical areas.

Description: This guide describes the input data required for using ECAP2D (Euler Cascade Aeroelastic Program-Two Dimensional). ECAP2D can be used for steady or unsteady aerodynamic and aeroelastic analysis of two dimensional cascades. Euler equations are used to obtain aerodynamic forces. The structural dynamic equations are written for a rigid typical section undergoing pitching (torsion) and plunging (bending) motion. The solution methods include harmonic oscillation method, influence coefficient method, pulse response method, and time integration method. For harmonic oscillation method, example inputs and outputs are provided for pitching motion and plunging motion. For the rest of the methods, input and output for pitching motion only are given.

Description: A recently developed micromechanical theory for the thermoelastic response of functionally graded composites with nonuniform fiber spacing in the through-thickness direction is further extended to enable analysis of material architectures characterized by arbitrarily nonuniform fiber spacing in two directions. In contrast to currently employed micromechanical approaches applied to functionally graded materials, which decouple the local and global effects by assuming the existence of a representative volume element at every point within the composite, the new theory explicitly couples the local and global effects. The analytical development is based on volumetric averaging of the various field quantities, together with imposition of boundary and interfacial conditions in an average sense. Results are presented that illustrate the capability of the derived theory to capture local stress gradients at the free edge of a laminated composite plate due to the application of a uniform temperature change. It is further shown that it is possible to reduce the magnitude of these stress concentrations by a proper management of the microstructure of the composite plies near the free edge. Thus by an appropriate tailoring of the microstructure it is possible to reduce or prevent the likelihood of delamination at free edges of standard composite laminates.

Description: This paper describes the work that is done at the Lehrstuhl fur Raumfahrttechnik (lrt) at the Technische Universitat Munchen to examine particle impacts into germanium surfaces which were flown on board the LDEF satellite. Besides the description of the processing of the samples, a brief overview of the particle launchers at our institute is given together with descriptions of impact morphology of high- and hypervelocity particles into germanium. Since germanium is a brittle, almost glass-like material, the impact morphology may also be interesting for anyone dealing with materials such as optics and solar cells. The main focus of our investigations is to learn about the impacting particle's properties, for example mass, velocity and direction. This is done by examining the morphology, various geometry parameters, crater obliqueness and crater volume.

Description: The investigation of the vibrational disturbances of the Hubble Space Telescope that were discovered soon after deployment in orbit is described in detail. It was found that the disturbances were particularly evident during orbital day-night crossings, and that the magnitudes of the disturbances were considerably larger than the design jitter requirements. This paper describes the process by which the vibrations were characterized and isolated to a particular mechanism. The analysis of the flight data and comparisons with computer simulation results showed that the source of the disturbances was the thermally driven deformation of the solar arrays in conjunction with frictional effects in the array mechanisms. The control system was successfully modified to attenuate the disturbances to tolerable levels pending mechanical and thermal redesign of the solar arrays. The new arrays were installed during the first space telescope servicing mission and, in combination with the enhanced control system algorithm, reduced the disturbances to satisfactory levels.

Description: This paper describes a computer program developed for structural dynamic analysis of horizontal axis wind turbines (HAWT's). It is based on the finite element method through its reliance on NASTRAN for the development of mass, stiffness, and damping matrices of the tower end rotor, which are treated in NASTRAN as separate structures. The tower is modeled in a stationary frame and the rotor in one rotating at a constant angular velocity. The two structures are subsequently joined together (external to NASTRAN) using a time-dependent transformation consistent with the hub configuration. Aerodynamic loads are computed with an established flow model based on strip theory. Aeroelastic effects are included by incorporating the local velocity and twisting deformation of the blade in the load computation. The turbulent nature of the wind, both in space and time, is modeled by adding in stochastic wind increments. The resulting equations of motion are solved in the time domain using the implicit Newmark-Beta integrator. Preliminary comparisons with data from the Boeing/NASA MOD2 HAWT indicate that the code is capable of accurately and efficiently predicting the response of HAWT's driven by turbulent winds.

Description: Underintegrated methods are investigated with respect to their stability and convergence properties. The focus was on identifying regions where they work and regions where techniques such as hourglass viscosity and hourglass control can be used. Results obtained show that underintegrated methods typically lead to finite element stiffness with spurious modes in the solution. However, problems exist (scalar elliptic boundary value problems) where underintegrated with hourglass control yield convergent solutions. Also, stress averaging in underintegrated stiffness calculations does not necessarily lead to stable or convergent stress states.

Description: A statistically generated weighting function for a second-order polynomial curve fit of residual functions has been developed. The residual flexibility test method, from which a residual function is generated, is a procedure for modal testing large structures in an external constraint-free environment to measure the effects of higher order modes and interface stiffness. This test method is applicable to structures with distinct degree-of-freedom interfaces to other system components. A theoretical residual function in the displacement/force domain has the characteristics of a relatively flat line in the lower frequencies and a slight upward curvature in the higher frequency range. In the test residual function, the above-mentioned characteristics can be seen in the data, but due to the present limitations in the modal parameter evaluation (natural frequencies and mode shapes) of test data, the residual function has regions of ragged data. A second order polynomial curve fit is required to obtain the residual flexibility term. A weighting function of the data is generated by examining the variances between neighboring data points. From a weighted second-order polynomial curve fit, an accurate residual flexibility value can be obtained. The residual flexibility value and free-free modes from testing are used to improve a mathematical model of the structure. The residual flexibility modal test method is applied to a straight beam with a trunnion appendage and a space shuttle payload pallet simulator.

Description: In elasticity, the method of forces, wherein stress parameters are considered as the primary unknowns, is known as the Beltrami-Michell formulation (BMF). The existing BMF can only solve stress boundary value problems; it cannot handle the more prevalent displacement of mixed boundary value problems of elasticity. Therefore, this formulation, which has restricted application, could not become a true alternative to the Navier's displacement method, which can solve all three types of boundary value problems. The restrictions in the BMF have been alleviated by augmenting the classical formulation with a novel set of conditions identified as the boundary compatibility conditions. This new method, which completes the classical force formulation, has been termed the completed Beltrami-Michell formulation (CBMF). The CBMF can solve general elasticity problems with stress, displacement, and mixed boundary conditions in terms of stresses as the primary unknowns. The CBMF is derived from the stationary condition of the variational functional of the integrated force method. In the CBMF, stresses for kinematically stable structures can be obtained without any reference to the displacements either in the field or on the boundary. This paper presents the CBMF and its derivation from the variational functional of the integrated force method. Several examples are presented to demonstrate the applicability of the completed formulation for analyzing mixed boundary value problems under thermomechanical loads. Selected example problems include a cylindrical shell wherein membrane and bending responses are coupled, and a composite circular plate.

Description: Structural analyses are developed to determine the linear elastic and the geometrically nonlinear elastic response of an internally pressurized, orthogonally stiffened, composite material cylindrical shell. The configuration is a long circular cylindrical shell stiffened on the inside by a regular arrangement of identical stringers and identical rings. Periodicity permits the analysis of a unit cell model consisting of a portion of the shell wall centered over one stringer-ring joint. The stringer-ring-shell joint is modeled in an idealized manner; the stiffeners are mathematically permitted to pass through one another without contact, but do interact indirectly through their mutual contact with the shell at the joint. Discrete beams models of the stiffeners include a stringer with a symmetrical cross section and a ring with either a symmetrical or an asymmetrical open section. Mathematical formulations presented for the linear response include the effect of transverse shear deformations and the effect of warping of the ring's cross section due to torsion. These effects are important when the ring has an asymmetrical cross section because the loss of symmetry in the problem results in torsion and out-of-plane bending of the ring, and a concomitant rotation of the joint at the stiffener intersection about the circumferential axis. Data from a composite material crown panel typical of a large transport fuselage structure are used for two numerical examples. Although the inclusion of geometric nonlinearity reduces the 'pillowing' of the shell, it is found that bending is localized to a narrow region near the stiffener. Including warping deformation of the ring into the analysis changes the sense of the joint rotation. Transverse shear deformation models result in increased joint flexibility.

Description: A methodology to compute cumulative probability distribution functions (CDF) of fatigue life for different ratios, r of applied stress to the laminate strength based on first ply failure criteria has been developed and demonstrated. Degradation effects due to long term environmental exposure and mechanical cyclic loads are considered in the simulation process. A unified time-stress dependent multi-factor interaction equation model developed at NASA Lewis Research Center has been used to account for the degradation/aging of material properties due to cyclic loads. Fast probability integration method is used to perform probabilistic simulation of uncertainties. Sensitivity of fatigue life reliability to uncertainties in the primitive random variables are computed and their significance in the reliability based design for maximum life is discussed. The results show that the graphite/epoxy (0/+45/90) deg laminate with ply thickness 0.125 in. has 500,000 cycles life for applied stress to laminate strength ratio of 0.6 and a reliability of 0.999. Also, the fatigue life reliability has been found to be most sensitive to the ply thickness and matrix tensile strength. Tighter quality controls must therefore be enforced on ply thickness and matrix strength in order to achieve high reliability of the structure.

Description: Specific forms for both the Gibb's and complementary dissipation potentials are chosen such that a complete (i.e., fully associative) potential base multiaxial, nonisothermal unified viscoplastic model is obtained. This model possesses one tensorial internal state variable (that is, associated with dislocation substructure) and an evolutionary law that has nonlinear kinematic hardening and both thermal and strain induced recovery mechanisms. A unique aspect of the present model is the inclusion of nonlinear hardening through the use of a compliance operator, derived from the Gibb's potential, in the evolution law for the back stress. This nonlinear tensorial operator is significant in that it allows both the flow and evolutionary laws to be fully associative (and therefore easily integrated), greatly influences the multiaxial response under non-proportional loading paths, and in the case of nonisothermal histories, introduces an instantaneous thermal softening mechanism proportional to the rate of change in temperature. In addition to this nonlinear compliance operator, a new consistent, potential preserving, internal strain unloading criterion has been introduced to prevent abnormalities in the predicted stress-strain curves, which are present with nonlinear hardening formulations, during unloading and reversed loading of the external variables. The specific model proposed is characterized for a representative titanium alloy commonly used as the matrix material in SiC fiber reinforced composites, i.e., TIMETAL 21S. Verification of the proposed model is shown using 'specialized' non-standard isothermal and thermomechanical deformation tests.

Description: Solving for the displacements of free-free coupled systems acted upon by static loads is commonly performed throughout the aerospace industry. Many times, these problems are solved using static analysis with inertia relief. This solution technique allows for a free-free static analysis by balancing the applied loads with inertia loads generated by the applied loads. For some engineering applications, the displacements of the free-free coupled system induce additional static loads. Hence, the applied loads are equal to the original loads plus displacement-dependent loads. Solving for the final displacements of such systems is commonly performed using iterative solution techniques. Unfortunately, these techniques can be time-consuming and labor-intensive. Since the coupled system equations for free-free systems with displacement-dependent loads can be written in closed-form, it is advantageous to solve for the displacements in this manner. Implementing closed-form equations in static analysis with inertia relief is analogous to implementing transfer functions in dynamic analysis. Using a MSC/NASTRAN DMAP Alter, displacement-dependent loads have been included in static analysis with inertia relief. Such an Alter has been used successfully to solve efficiently a common aerospace problem typically solved using an iterative technique.

Description: The fast fracture strength distribution of uniaxially ground, alpha silicon carbide was investigated as a function of grinding angle relative to the principal stress direction in flexure. Both as-ground and ground/annealed surfaces were investigated. The resulting flexural strength distributions were used to verify reliability models and predict the strength distribution of larger plate specimens tested in biaxial flexure. Complete fractography was done on the specimens. Failures occurred from agglomerates, machining cracks, or hybrid flaws that consisted of a machining crack located at a processing agglomerate. Annealing eliminated failures due to machining damage. Reliability analyses were performed using two and three parameter Weibull and Batdorf methodologies. The Weibull size effect was demonstrated for machining flaws. Mixed mode reliability models reasonably predicted the strength distributions of uniaxial flexure and biaxial plate specimens.

Description: Probabilistic predictions based on the IPACS code are presented for the material and structural response of unnotched and notched, IM6/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 judged poor because IPACS did not have a progressive failure capability at the time this work was performed. The report 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 cyclic deformation behavior of a wrought cobalt-base superalloy, Haynes 188, has been investigated over a range of temperatures between 25 and 1000 C under isothermal and in-phase thermomechanical fatigue (TMF) conditions. Constant mechanical strain rates (epsilon-dot) of 10(exp -3)/s and 10(exp -4)/s were examined with a fully reversed strain range of 0.8%. Particular attention was given to the effects of dynamic strain aging (DSA) on the stress-strain response and low cycle fatigue life. A correlation between cyclic deformation behavior and microstructural substructure was made through detailed transmission electron microscopy. Although DSA was found to occur over a wide temperature range between approximately 300 and 750 C the microstructural characteristics and the deformation mechanisms responsible for DSA varied considerably and were dependent upon temperature. In general, the operation of DSA processes led to a maximum of the cyclic stress amplitude at 650 C and was accompanied by pronounced planar slip, relatively high dislocation density, and the generation of stacking faults. DSA was evidenced through a combination of phenomena, including serrated yielding, an inverse dependence of the maximum cyclic hardening with epsilon-dot, and an instantaneous inverse epsilon-dot sensitivity verified by specialized epsilon-dot -change tests. The TMF cyclic hardening behavior of the alloy appeared to be dictated by the substructural changes occuring at the maximum temperature in the TMF cycle.

Description: In spite of great improvements in accuracy through the use of computers, design methods, which can be equally critical in establishing the commercial success of a material, have been treated as afterthoughts. Early investment in design and development technologies can easily reduce manufacturing costs later in the product cycle. To avoid lengthy product development times for ceramic composites, funding agencies for materials research must commit resources to support design and development technologies early in the material life cycle. These technologies need not focus on designing the material, rather, the technology must focus on designing with the material, i. e., developing methods to design components fabricated from the new material. Thus a basic tenet that motivated this research effort is that a persistent need exists for improvements in the analysis of components fabricated from CMC material systems. From an aerospace design engineer's perspective the new generation of ceramic composites offers a significant potential for raising the thrust/weight ratio and reducing NOx emissions of gas turbine engines. Continuous ceramic fiber composites exhibit an increase in work of fracture, which allows for 'graceful' rather than catastrophic failure. When loaded in the fiber direction, these composites retain substantial strength capacity beyond the initiation of transverse matrix cracking despite the fact that neither of its constituents would exhibit such behavior if tested alone. As additional load is applied beyond first matrix cracking, the matrix tends to break in a series of cracks bridged by the ceramic fibers. Thus any additional load is born increasingly by the fibers until the ultimate strength of the composite is reached. Establishing design protocols that enable the engineer to analyze and predict this type of behavior in ceramic composites was the general goal of this project.

Description: A time-domain study of the random response of a laminated plate subjected to combined acoustic and thermal loads is carried out. The features of this problem also include given uniform static inplane forces. The formulation takes into consideration a possible initial imperfection in the flatness of the plate. High decibel sound pressure levels along with high thermal gradients across thickness drive the plate response into nonlinear regimes. This calls for the analysis to use von Karman large deflection strain-displacement relationships. A finite element model that combines the von Karman strains with the first-order shear deformation plate theory is developed. The development of the analytical model can accommodate an anisotropic composite laminate built up of uniformly thick layers of orthotropic, linearly elastic laminae. The global system of finite element equations is then reduced to a modal system of equations. Numerical simulation using a single-step algorithm in the time-domain is then carried out to solve for the modal coordinates. Nonlinear algebraic equations within each time-step are solved by the Newton-Raphson method. The random gaussian filtered white noise load is generated using Monte Carlo simulation. The acoustic pressure distribution over the plate is capable of accounting for a grazing incidence wavefront. Numerical results are presented to study a variety of cases.

Description: Launch of ionosphere explorer-a /ie-a/ satellite structure, orbit, scout launch vehicle, tracking and data aquisition, and scientific experiments

Description: Summary of mission of ranger vii lunar probe and nature of lunar surface

Description: Ranger vii lunar probe - flight chronology, launch vehicle, spacecraft, and photography

Description: Agreement for esro i and esro ii satellites between nasa and european space research organization /esro/

Description: Nasa launching of satellite to gather information on meteoroid hazards in space flight

Description: Description and scientific experiments of imp-b interplanetary explorer satellite

Description: Echo ii communications satellite launching

Description: Relay satellite destruction failure, anticipated lifespan, and history of transmissions

Description: Radiosonde for measuring intensity of cosmic rays in stratosphere

Description: Stratosphere-mesosphere meteorological rocket

Description: Model for determining mean cumulative flux of meteoroids in cislunar space

Description: Transient response of multistage space vehicle to external load

Description: Wall pressure fluctuations and skin vibrations - free shear layers and oscillating shocks

Description: Fluctuating pressure requirements for calculation of flight vehicle vibration environments - structural response equations

Description: Nonlinear partial differential equations of two variables for vibratory motions of beam

Description: Bending moment of vehicle in terms of control parameters - vehicle elastic and sloshing modes

Description: Relay communications satellite launching aboard delta rocket

Description: Mission and preparation of ranger b lunar probe

Description: Probability that stress is less than strength at prescribed confidence levels for normally distributed data

Description: Earth-landing systems for mercury, gemini and apollo spacecraft, including parachutes and paragliders

Description: Palmgren-Miner, Manson, Fuller, and Valluri methods of predicting fatigue life in spectrum loading

Description: Plastic strains and energy density in thin cracked plates, discussing ductile fracture

Description: Spin stabilization and ion thrustor engines for space electric rocket test /sert/ spacecraft

Description: Nimbus satellite - spacecraft, subsystems, and launch vehicle