Search

This program began as a grant from NASA Headquarters, NGT-10025, which was in effect from 10/l/93 until 10/31/96. The remaining funding for this effort was transferred from NASA Headquarters to NASA Langley and a new grant NGT-1-52155 was issued covering the period II/l/96 to 5/15/99. This report serves as the final report of NGT-1-52155. For a number of years, Virginia Tech had been on the forefront of research in the area of multidisciplinary analysis and design. In June of 1994, faculty members from aerospace and ocean engineering, engineering science and mechanics, mechanical engineering, industrial engineering, mathematics and computer sciences, at Virginia Tech joined together to form the Multidisciplinary Analysis and Design (MAD) Center for Advanced Vehicles. The center was established with the single goal: to perform research that is relevant to the needs of the US industry and to foster collaboration between the university, government and industry. In October of 1994, the center was chosen by NASA headquarters as one of the five university centers to establish a fellowship program to develop a graduate program in multidisciplinary analysis and design. The fellowship program provides full stipend and tuition support for seven U. S. students per year during their graduate studies. The grant is currently being administered by the NMO Branch of NASA Langley. To advise us regarding the problems faced by the industry, an industrial advisory board has been formed consisting of representatives from industry as well as government laboratories. The present membership includes major aerospace companies: Aurora Flight Sciences, Boeing: Philadelphia, Boeing: Long Beach, Boeing: Seattle, Boeing: St. Louis, Cessna, Ford, General Electric, Hughes, Lockheed-Martin: Palo Alto, Northrop-Grumman, Sikorsky, smaller, aerospace software companies: Aerosoft, Phoenix Integration and Proteus Engineering, along with representatives from government agencies, including: NASA Ames, Langley and Lewis. The function of the advisory board is to channel information from its member companies to faculty members concerning problems that need research attention in the general area of multidisciplinary design optimization (MDO). The faculty and their graduate students make presentations to the board on their research. The board makes recommendations on the research and suggests new areas and problems which need attention. Many students participating in the program spend 3-6 months in industry working on their research projects. We are completing the fifth year of the fellowship program and have had four advisory board meetings in Blacksburg. Ten students have spent the three month periods in industry. In addition to the research element of the MAD Center efforts we also have an academic component. We have developed a menu of design-related graduate courses and two new courses: one in Aerospace Manufacturing and another in MDO. Some of the MAD Center activities are described on the world-wide web at http://www.aoe.vt.edu/mads.html The MAD Center represents an innovative approach for joint Industry-Government-University cooperation in the development of a comprehensive program in engineering education which addresses the design needs of industry. The following charts list detail of the grant: mission of the MAD center, faculty members, purpose of the advisory board, board members, summary of the graduate and undergraduate program, history of the fellowship program, mission of the fellowship program, requirements of MAD fellows, course requirements, students supported, advisory board participants, and MAD center research papers

The present minimum-mass designs for a two-spar rectangular box beam were derived on the basis of linear-buckling FEM analysis constraints. In order to ascertain the effects of any geometric nonlinearities on these designs, each was subjected to a geometrically nonlinear FEM analysis. In all cases, the structure collapses below the design load, and does so in a mode which differs from that of linear theory. This discrepancy is attributable to such nonlinear panel-interaction mechanisms as rib-crusing loads. The optimized design is highly sensitive to crushing loads, relative to the nonoptimal design.

Genetic algorithms (GAs) readily handle discrete problems, and can be made to generate many optima, as is presently illustrated for the case of design for minimum-weight stiffened panels with buckling constraints. The GA discrete design procedure proved superior to extant alternatives for both stiffened panels with cutouts and without cutouts. High computational costs are, however, associated with this discrete design approach at the current level of its development.

A continuous beam model is developed for approximate analysis of a large, slender, beam-like truss. The model is incorporated in a multi-level optimization scheme for the weight minimization of such trusses. This scheme is tested against traditional optimization procedures for savings in computational cost. Results from both optimization methods are presented for comparison.

An integer programming formulation for the design of symmetric and balanced rectangular composite laminates with simply supported boundary conditions subject to buckling and strain constraints is presented. The design variables that define the stacking sequence of the laminate are ply-identity zero-one integers. The buckling constraint is linear in terms of the ply-identity design variables, but strains are nonlinear functions of these variables. A linear approximation is developed for the strain constraints so that the problem can be solved by sequential linearization using the branch and bound algorithm. Examples of graphite-epoxy plates under biaxial compression are presented. Optimum stacking sequences obtained using the linear approximation are compared with global optimum designs obtained using a genetic search procedure.

The new concept of using a continuum model in the multilevel decomposition of structural optimization problems is presented. The practicality of this concept is demonstrated by application in a scheme for the optimization of beam-like space structures. This scheme is tested against traditional optimization procedures for savings in computational cost. Results from both optimization methods are presented for comparison.

An optimization system based on general-purpose finite element code CSM Testbed and optimization program ADS is described. The system can be used to obtain minimum-mass designs of composite shell structures with complex stiffening arrangements. Ply thicknesses, ply orientations, and stiffener heights can be used as design variables. Buckling, displacement, and material failure constraints can be imposed on the design. The system is used to conduct a preliminary design study of geodesically stiffened shells. For comparison purposes, optimal designs of unstiffened shells, and ring and longitudinal stringer stiffened shells are also studied. Trends in the design of geodesically stiffened shells are identified. Features that enhance the capabilities and efficiency of the design system are described.

A procedure is developed analytically to design and produce an optimized stiffened panel with a centrally located hole. The PASCO panel-buckling analysis and sizing code is employed with the Engineering Analysis Language (EAL) code for finite-element analysis to optimize the design variables. Important design considerations are stiffener height and ply thicknesses in the skin and the stiffener laminates. The EAL code calculates the maximum strains around the hole boundary, and the resulting constraints are combined with the buckling constraint developed by the PASCO code. Cooptimization of the two codes is effected until convergence on the optimal solution is achieved. Results of the analytical procedure are experimentally verified by applying a compressive end-load to graphite-fiber laminated plates with center holes. The experimental results support the analytical procedure, and it is found that for the same weight, the soft-skin designs tested are found to have the best tolerance characteristics and resistance to failure.

Experiments were performed on center-notched specimens of unidirectional AS4/3501-6 graphite-epoxy subjected to varied far-field loading conditions. Specimens tested include on-axis and off-axis tensile coupons and Iosipescu shear specimens, allowing the study of normal, shear, and combined far-field loading conditions. From the tests, the direction of initial crack extension and the far-field stresses at crack initiation and fracture are obtained. The normal stress ratio theory is subsequently applied within a macroscopic-level homogeneous anisotropic stress solution to predict the direction of crack extension, critical stresses, and other crack growth behavior for the experiments. Comparison between experimental results and theoretical predictions indicates that application of the normal stress ratio theory within a macroscopic-level stress analysis has significant potential as an approach to modeling composite crack growth.

An analytical study is presented showing the effects of the notch tip geometry on the location and direction of crack growth from an existing notch in a unidirectional fibrous composite modelled as a homogeneous, anisotropic, elastic material. Anisotropic elasticity and the normal stress ratio theory are used to study crack growth from elliptical notches in unidirectional composites. Sharp cracks, circular holes, and ellipses are studied under far-field tension and shear loading. Limited comparisons are made showing good correlation with experimental results.

Using lifting-line theory and beam analysis, the geometry (planiform and twist) and composite material structural sizes (skin thickness, spar cap, and web thickness) were designed for a sailplane wing, subject to both structural and aerodynamic constraints. For all elements, the integrated design (simultaneously designing the aerodynamics and the structure) was superior in terms of performance and weight to the sequential design (where the aerodynamic geometry is designed to maximize the performance, following which a structural/aeroelastic design minimizes the weight). Integrated designs produced less rigid, higher aspect ratio wings with favorable aerodynamic/structural interactions.

The paper introduces a failure model for laminated composite plates with a cutout under combined compressive and shear loads. The model is based on kinking failure of the load-carrying fibers around a cutout, and includes the effect of local shearing and compressive stresses. Comparison of predictions of the model with available experimental results for quasi-isotropic and orthotropic plates with a circular hole indicated a good agreement. Predictions for orthotropic plates under combined loading are compared with the predictions of a point-stress model. The present model indicates significant reductions in axial load-carrying capacity due to shearing loads for plates with principal axis of orthotropy oriented along the axial load direction. A gain in strength is achieved by rotating the axis of orthotropy to counteract the shearing stress, or by eliminating the compressive-shear deformation coupling.

An experimental study was conducted to evaluate the effect of lateral slots on the buckling response, the postbuckling response, and the failure characteristics of flat rectangular graphite-epoxy plates loaded in compression. The slots did not significantly affect the prebuckling and buckling behavior of the plates. The slots caused local changes in the strain distribution and out-of-plane deformations near the slot. Failure loads and modes were strongly affected by slot location. Centrally located slots did not affect the failure loads or failure modes of the plates. Slots located close to an unloaded edge of a plate reduced the postbuckling strength of the plate and affected the failure mode. For one slightly offset slot location, interaction of out-of-plane deformations with the slot tip was observed to affect the failure mode significantly without affecting the failure load.

An automated procedure for designing minimum-weight composite panels subject to a local damage constraint under tensile loading was developed. A finite element program based on linear elastic fracture mechanics for calculating stress intensity factors (SIF) was incorporated in the design cycle. Panel fracture toughness was obtained by using a strain based criterion. A general purpose mathematical optimization algorithm was used for the weight minimization. Analytical sensitivity derivatives of the SIF employing the adjoint variable technique was used to enhance the computational efficiency of the procedure. Design results for both unstiffened and stiffened plates are presented.

Attention is given to the application of the adjoint variable method under the special combination of displacement and stress boundary conditions which occur when static test conditions are simulated. In the two-step problem-solving method employed, calculation of internal stresses, displacements and total reaction force at the boundary is followed by the calculation of the applied load and the total reaction force due to unit displacements.