Search

This monograph is organized to highlight the successful application of light alloys on aircraft and space launch vehicles, the role of NASA in enabling these applications for each different class of flight vehicles, and a discussion of the major advancements made in discipline areas of research. In each section, key personnel and selected references are included. These references are intended to provide additional information for technical specialists and others who desire a more in-depth discussion of the contributions. Also in each section, lessons learned and future challenges are highlighted to help guide technical personnel either in the conduct or management of current and future research projects related to light-weighting advanced air and space vehicles.

The development of dynamic mesh capability for turbulent flow simulations using the Streamlined Upwind Petrov-Galerkin (SUPG) discretization is described. The current work extends previous research to include high-order spatial accuracy, including the satisfaction of the discrete geometric conservation law (GCL) on curved elements. Two closely-related schemes are described and the ability of these schemes to satisfy the GCL, while also maintaining temporal accuracy and conservation is assessed. Studies indicate that although one scheme discretizes the time derivative in conservative form, both schemes exhibit temporal conservation errors that decrease according to the expected design order of accuracy. The source of the temporal conservation errors is examined, and it is demonstrated that many finite-volume and finite-element schemes can also be expected to have difficulty strictly satisfying conservation in time. The effects on conservation are examined and, while present in the simulations, are seen to be negligible for the problems considered.

A Streamlined Upwind Petrov-Galerkin (SUPG) stabilized finite-element discretization has been implemented as a library into the FUN3D unstructured-grid flow solver. Motivation for the selection of this methodology is given, details of the implementation are provided, and the discretization for the interior scheme is verified for linear and quadratic elements by using the method of manufactured solutions. A methodology is also described for capturing shocks, and simulation results are compared to the finite-volume formulation that is currently the primary method employed for routine engineering applications. The finite-element methodology is demonstrated to be more accurate than the finite-volume technology, particularly on tetrahedral meshes where the solutions obtained using the finite-volume scheme can suffer from adverse effects caused by bias in the grid. Although no effort has been made to date to optimize computational efficiency, the finite-element scheme is competitive with the finite-volume scheme in terms of computer time to reach convergence.

The growth of a fatigue crack through a material is the result of a complex interaction between the applied loading, component geometry, three-dimensional constraint, load history, environment, material microstructure and several other factors. Previous studies have developed experimental and computational methods to relate the fatigue crack growth rate to many of the above conditions, with the intent of discovering some fundamental material response, i.e. crack growth rate as a function of something. Currently, the technical community uses the stress intensity factor solution as a simplistic means to relate fatigue crack growth rate to loading, geometry and all other variables. The stress intensity factor solution is a very simple linear-elastic representation of the continuum mechanics portion of crack growth. In this paper, the authors present fatigue crack growth rate data for two different high strength steel alloys generated using standard methods. The steels exhibit behaviour that appears unexplainable, compared to an aluminium alloy presented as a baseline for comparison, using the stress intensity factor solution.

This paper presents numerical simulations for a high-speed transonic compressor rotor, NASA Rotor 35, using an unstructured grid flow solver U2NCLE. It solves Reynolds averaged Navier-Stokes equations using an arbitrary Mach number solution algorithm. Experiments were conducted at the NASA Glen Research Center to investigate compressor stability enhancement using steady and discrete tip injections. The mass injection upstream of the tip of a high- speed axial compressor rotor has shown effectiveness in suppressing stall in tip-critical rotors. This work presents numerical simulations of a high-speed compressor rotor with and without tip injections. Methods to model the flow injection are described. Numerical simulations are conducted at 80 and 100 percent of the design speed. Results show that numerical solutions are consistent with the experimental trends.

The criteria used to prevent failure of wind-tunnel models and support hardware were revised as part of a project to enhance the capabilities of cryogenic wind tunnel testing at NASA Langley Research Center. Specifically, damage-tolerance fatigue life prediction methods are now required for critical components, and material selection criteria are more general and based on laboratory test data. The suitability of two candidate model alloys (AerMet 100 and C-250 steel) was investigated by obtaining the fatigue crack growth and fracture data required for a damage-tolerance fatigue life analysis. Finally, an example is presented to illustrate the newly implemented damage tolerance analyses required of wind-tunnel model system components.

A series of fracture tests were conducted on Middle-crack tension M(T) and compact tension C(T) specimens to determine the effects of specimen type, specimen width, notch tip sharpness and buckling on the fracture behavior of cracked thin sheet (0.04 inch thick) 2024-T3 aluminum alloy material. A series of M(T) specimens were tested with three notch tip configurations: (1) a fatigue pre-cracked notch, (2) a 0.010-inch-diameter wire electrical discharge machined (EDM) notch, and (3) a EDM notch sharpened with a razor blade. The test procedures are discussed and the experimental results for failure stress, load vs. crack extension and the material stress-strain response are reported.

The fatigue crack growth threshold, defining crack growth as either very slow or nonexistent, has been traditionally determined with standardized load reduction methodologies. Some experimental procedures tend to induce load history effects that result in remote crack closure from plasticity. This history can affect the crack driving force, i.e. during the unloading process the crack will close first at some point along the wake, reducing the effective load at the crack tip. One way to reduce the effects of load history is to propagate a crack under constant amplitude loading. As a crack propagates under constant amplitude loading, the stress intensity factor, K, will increase, as will the crack growth rate, da/dN. A fatigue crack growth threshold test procedure is developed and experimentally validated that does not produce load history effects and can be conducted at a specified stress ratio, R.

A summary of experimental data is presented for a series of fracture tests and tensile tests that were conducted on a 2324-T39 aluminum alloy material. Static tensile tests were performed to determine the stress-strain response for the material. For the fracture tests three thicknesses (B = 0.5, 0.7, and 0.9-inches) and two specimen types (middle crack tension and compact tension specimens) were used. The middle crack tension specimens were 12-inches wide, and the compact tension specimens were 2, 4, and 6-inches wide. The fracture tests were conducted in the L-T orientation (crack perpendicular to the rolling direction and load parallel to the rolling direction).

A series of fracture test were conducted to determine the effects of specimen type specimen width and buckling on the fracture behavior of cracked thin sheet (0.063 inch thick) 2024-T3 aluminum alloy. A summary of the experimental measurements is presented for fracture tests conducted on two specimen types and various widths. Middle-crack tension M(T) and compact tension C(T) specimens were tested in the L-T and T-L orientation with duplicate tests for each condition. Four widths (W= 3, 12, 24, and 40 inch) were tested for the middle-crack tension specimens, and three widths (W=2, 4, and 6 inch) were tested for the compact tension specimens. The M(T) specimens were tested in either a constrained (out-of-plane displacements restrained with antibuckling guides) or unconstrained conditions were the specimen was free to buckle out of plane Measurements were made of load against crack extension for all specimens.

A multidisciplinary sensitivity analysis technique that has been shown to be independent of step-size selection is examined further. The accuracy of this step-size independent technique, which uses complex variables for determining sensitivity derivatives, has been previously established. The primary focus of this work is to validate the aero-structural analysis procedure currently being used. This validation consists of comparing computed and experimental data obtained for an Aeroelastic Research Wing (ARW-2). Since the aero-structural analysis procedure has the complex variable modifications already included into the software, sensitivity derivatives can automatically be computed. Other than for design purposes, sensitivity derivatives can be used for predicting the solution at nearby conditions. The use of sensitivity derivatives for predicting the aero-structural characteristics of this configuration is demonstrated.

An extensive experimental database has been assembled from very detailed teardown examinations of fatigue cracks found in rivet holes of fuselage structural components. Based on this experimental database, a comprehensive analysis methodology was developed to predict the onset of widespread fatigue damage in lap joints of fuselage structure. Several computer codes were developed with specialized capabilities to conduct the various analyses that make up the comprehensive methodology. Over the past several years, the authors have interrogated various aspects of the analysis methods to determine the degree of computational rigor required to produce numerical predictions with acceptable engineering accuracy. This study led to the formulation of a practical engineering approach to predicting fatigue crack growth in riveted lap joints. This paper describes the practical engineering approach and compares predictions with the results from several experimental studies.

This paper presents a brief overview of some of the more recent advances in steady aerodynamic shape-design sensitivity analysis and optimization, based on advanced computational fluid dynamics (CFD). The focus here is on those methods particularly well-suited to the study of geometrically complex configurations and their potentially complex associated flow physics. When nonlinear state equations are considered in the optimization process, difficulties are found in the application of sensitivity analysis. Some techniques for circumventing such difficulties are currently being explored and are included here. Attention is directed to methods that utilize automatic differentiation to obtain aerodynamic sensitivity derivatives for both complex configurations and complex flow physics. Various examples of shape-design sensitivity analysis for unstructured-grid CFD algorithms are demonstrated for different formulations of the sensitivity equations. Finally, the use of advanced, unstructured-grid CFDs in multidisciplinary analyses and multidisciplinary sensitivity analyses within future optimization processes is recommended and encouraged.

An assessment of the plasticity-induced crack-closure concept is made, in light of some of the questions that have been raised on the validity of the concept, and the assumptions that have been made concerning crack-dp damage below the crack-opening stress. The impact of using other crack-tip parameters, such as the cyclic crack-tip displacement, to model crack-growth rate behavior was studied. Crack-growth simulations, using a crack-closure model, showed a close relation between traditional Delta K eff, and the cyclic crack-tip displacement (Delta eff) for an aluminum alloy and a steel. Evaluations of the cyclic hysteresis energy demonstrated that the cyclic plastic damage below the crack-opening stress was negligible in the Paris crack-growth regime. Some of the standard and newly proposed remote measurement methods to determine the 'effective' crack-tip driving parameter were evaluated on middle-crack tension specimens. A potential source of the Kmax effect on crack-growth rates was studied on an aluminum alloy. Results showed that the ratio of Kmax to Kc had a strong effect on crack-growth rates at high stress ratios and at low stress ratios for very high stress levels. The crack-closure concept and the traditional crack-growth rate equations were able to correlate and predict crack-growth rates under these extreme conditions.

This paper presents a brief overview of some of the more recent advances in steady aerodynamic shape-design sensitivity analysis and optimization, based on advanced computational fluid dynamics. The focus here is on those methods particularly well- suited to the study of geometrically complex configurations and their potentially complex associated flow physics. When nonlinear state equations are considered in the optimization process, difficulties are found in the application of sensitivity analysis. Some techniques for circumventing such difficulties are currently being explored and are included here. Attention is directed to methods that utilize automatic differentiation to obtain aerodynamic sensitivity derivatives for both complex configurations and complex flow physics. Various examples of shape-design sensitivity analysis for unstructured-grid computational fluid dynamics algorithms are demonstrated for different formulations of the sensitivity equations. Finally, the use of advanced, unstructured-grid computational fluid dynamics in multidisciplinary analyses and multidisciplinary sensitivity analyses within future optimization processes is recommended and encouraged.

A three-dimensional unstructured grid approach to aerodynamic shape sensitivity analysis and design optimization has been developed and is extended to model geometrically complex configurations. The advantage of unstructured grids (when compared with a structured-grid approach) is their inherent ability to discretize irregularly shaped domains with greater efficiency and less effort. Hence, this approach is ideally suited for geometrically complex configurations of practical interest. In this work the nonlinear Euler equations are solved using an upwind, cell-centered, finite-volume scheme. The discrete, linearized systems which result from this scheme are solved iteratively by a preconditioned conjugate-gradient-like algorithm known as GMRES for the two-dimensional geometry and a Gauss-Seidel algorithm for the three-dimensional; similar procedures are used to solve the accompanying linear aerodynamic sensitivity equations in incremental iterative form. As shown, this particular form of the sensitivity equation makes large-scale gradient-based aerodynamic optimization possible by taking advantage of memory efficient methods to construct exact Jacobian matrix-vector products. Simple parameterization techniques are utilized for demonstrative purposes. Once the surface has been deformed, the unstructured grid is adapted by considering the mesh as a system of interconnected springs. Grid sensitivities are obtained by differentiating the surface parameterization and the grid adaptation algorithms with ADIFOR (which is an advanced automatic-differentiation software tool). To demonstrate the ability of this procedure to analyze and design complex configurations of practical interest, the sensitivity analysis and shape optimization has been performed for a two-dimensional high-lift multielement airfoil and for a three-dimensional Boeing 747-200 aircraft.

ZIP2DL is a two-dimensional, elastic-plastic finte element program for stress analysis and crack growth simulations, developed for the NASA Langley Research Center. It has many of the salient features of the ZIP2D program. For example, ZIP2DL contains five material models (linearly elastic, elastic-perfectly plastic, power-law hardening, linear hardening, and multi-linear hardening models), and it can simulate mixed-mode crack growth for prescribed crack growth paths under plane stress, plane strain and mixed state of stress conditions. Further, as an extension of ZIP2D, it also includes a number of new capabilities. The large-deformation kinematics in ZIP2DL will allow it to handle elastic problems with large strains and large rotations, and elastic-plastic problems with small strains and large rotations. Loading conditions in terms of surface traction, concentrated load, and nodal displacement can be applied with a default linear time dependence or they can be programmed according to a user-defined time dependence through a user subroutine. The restart capability of ZIP2DL will make it possible to stop the execution of the program at any time, analyze the results and/or modify execution options and resume and continue the execution of the program. This report includes three sectons: a theoretical manual section, a user manual section, and an example manual secton. In the theoretical secton, the mathematics behind the various aspects of the program are concisely outlined. In the user manual section, a line-by-line explanation of the input data is given. In the example manual secton, three types of examples are presented to demonstrate the accuracy and illustrate the use of this program.

The seventh Jerry L. Swedlow Memorial Lecture presents a review of some of the technical developments, that have occurred during the past 40 years, which have led to the merger of fatigue and fracture mechanics concepts. This review is made from the viewpoint of 'crack propagation.' As methods to observe the 'fatigue' process have improved, the formation of fatigue micro-cracks have been observed earlier in life and the measured crack sizes have become smaller. These observations suggest that fatigue damage can now be characterized by 'crack size.' In parallel, the crack-growth analysis methods, using stress-intensity factors, have also improved. But the effects of material inhomogeneities, crack-fracture mechanisms, and nonlinear behavior must now be included in these analyses. The discovery of crack-closure mechanisms, such as plasticity, roughness, and oxide/corrosion/fretting product debris, and the use of the effective stress-intensity factor range, has provided an engineering tool to predict small- and large-crack-growth rate behavior under service loading, conditions. These mechanisms have also provided a rationale for developing, new, damage-tolerant materials. This review suggests that small-crack growth behavior should be viewed as typical behavior, whereas large-crack threshold behavior should be viewed as the anomaly. Small-crack theory has unified 'fatigue' and 'fracture mechanics' concepts; and has bridged the cap between safe-life and durability/damage-tolerance design concepts.

NASA has developed a comprehensive analytical methodology for predicting the onset of widespread fatigue damage in fuselage structure. The determination of the number of flights and operational hours of aircraft service life that are related to the onset of widespread fatigue damage includes analyses for crack initiation, fatigue crack growth, and residual strength. Therefore, the computational capability required to predict analytically the onset of widespread fatigue damage must be able to represent a wide range of crack sizes from the material (microscale) level to the global structural-scale level. NASA studies indicate that the fatigue crack behavior in aircraft structure can be represented conveniently by the following three analysis scales: small three-dimensional cracks at the microscale level, through-the-thickness two-dimensional cracks at the local structural level, and long cracks at the global structural level. The computational requirements for each of these three analysis scales are described in this paper.