ALBERT

Description: A thin-spherical shell loaded as a cantilever beam-stresses, deflections, and influence coefficients

Description: Modes of nonuniform thin circular rings

Description: Design charts of continuous film, self-acting journal bearings of finite length

Description: Dynamic response of rectangular plate to series of moving loads

Description: Impact performance of composite spacecraft structures

Description: Momentum transfer and cratering from aluminum sphere hypervelocity impact

Description: The analysis of balloon envelopes by the finite element (FE) method is plagued by convergence problems. A pratical FE analysis approach is based on the fact that in thin shells with non-zero Gaussian curvature the membrane solution component is essentially decoupled from the bending solution component. A proxy-problem is solved by using a small artificial bending stiffness that assures convergence without significantly affecting the membrane solution component. This approach has been previously validated on slightly overpressurized balloon envelopes. Extensions of this approach to more difficult problems in the structural analysis of balloon envelopes are presented. The convergence forcing modelling measures are discussed. Implications of the findings of the analysis results to future balloon designs are also discussed.

Description: The Composite Plate Buckling Analysis Program (COMPPAP) was written to help engineers determine buckling loads of orthotropic (or isotropic) irregularly shaped plates without requiring hand calculations from design curves or extensive finite element modeling. COMPPAP is a one element finite element program that utilizes high-order displacement functions. The high order of the displacement functions enables the user to produce results more accurate than traditional h-finite elements. This program uses these high-order displacement functions to perform a plane stress analysis of a general plate followed by a buckling calculation based on the stresses found in the plane stress solution. The current version assumes a flat plate (constant thickness) subject to a constant edge load (normal or shear) on one or more edges. COMPPAP uses the power method to find the eigenvalues of the buckling problem. The power method provides an efficient solution when only one eigenvalue is desired. Once the eigenvalue is found, the eigenvector, which corresponds to the plate buckling mode shape, results as a by-product. A positive feature of the power method is that the dominant eigenvalue is the first found, which is this case is the plate buckling load. The reported eigenvalue expresses a load factor to induce plate buckling. COMPPAP is written in ANSI FORTRAN 77. Two machine versions are available from COSMIC: a PC version (MSC-22428), which is for IBM PC 386 series and higher computers and compatibles running MS-DOS; and a UNIX version (MSC-22286). The distribution medium for both machine versions includes source code for both single and double precision versions of COMPPAP. The PC version includes source code which has been optimized for implementation within DOS memory constraints as well as sample executables for both the single and double precision versions of COMPPAP. The double precision versions of COMPPAP have been successfully implemented on an IBM PC 386 compatible running MS-DOS, a Sun4 series computer running SunOS, an HP-9000 series computer running HP-UX, and a CRAY X-MP series computer running UNICOS. COMPPAP requires 1Mb of RAM and the BLAS and LINPACK math libraries, which are included on the distribution medium. The COMPPAP documentation provides instructions for using the commercial post-processing package PATRAN for graphical interpretation of COMPPAP output. The UNIX version includes two electronic versions of the documentation: one in LaTex format and one in PostScript format. The standard distribution medium for the PC version (MSC-22428) is a 5.25 inch 1.2Mb MS-DOS format diskette. The standard distribution medium for the UNIX version (MSC-22286) is a .25 inch streaming magnetic tape cartridge (Sun QIC-24) in UNIX tar format. For the UNIX version, alternate distribution media and formats are available upon request. COMPPAP was developed in 1992.

Description: The Composite Plate Buckling Analysis Program (COMPPAP) was written to help engineers determine buckling loads of orthotropic (or isotropic) irregularly shaped plates without requiring hand calculations from design curves or extensive finite element modeling. COMPPAP is a one element finite element program that utilizes high-order displacement functions. The high order of the displacement functions enables the user to produce results more accurate than traditional h-finite elements. This program uses these high-order displacement functions to perform a plane stress analysis of a general plate followed by a buckling calculation based on the stresses found in the plane stress solution. The current version assumes a flat plate (constant thickness) subject to a constant edge load (normal or shear) on one or more edges. COMPPAP uses the power method to find the eigenvalues of the buckling problem. The power method provides an efficient solution when only one eigenvalue is desired. Once the eigenvalue is found, the eigenvector, which corresponds to the plate buckling mode shape, results as a by-product. A positive feature of the power method is that the dominant eigenvalue is the first found, which is this case is the plate buckling load. The reported eigenvalue expresses a load factor to induce plate buckling. COMPPAP is written in ANSI FORTRAN 77. Two machine versions are available from COSMIC: a PC version (MSC-22428), which is for IBM PC 386 series and higher computers and compatibles running MS-DOS; and a UNIX version (MSC-22286). The distribution medium for both machine versions includes source code for both single and double precision versions of COMPPAP. The PC version includes source code which has been optimized for implementation within DOS memory constraints as well as sample executables for both the single and double precision versions of COMPPAP. The double precision versions of COMPPAP have been successfully implemented on an IBM PC 386 compatible running MS-DOS, a Sun4 series computer running SunOS, an HP-9000 series computer running HP-UX, and a CRAY X-MP series computer running UNICOS. COMPPAP requires 1Mb of RAM and the BLAS and LINPACK math libraries, which are included on the distribution medium. The COMPPAP documentation provides instructions for using the commercial post-processing package PATRAN for graphical interpretation of COMPPAP output. The UNIX version includes two electronic versions of the documentation: one in LaTex format and one in PostScript format. The standard distribution medium for the PC version (MSC-22428) is a 5.25 inch 1.2Mb MS-DOS format diskette. The standard distribution medium for the UNIX version (MSC-22286) is a .25 inch streaming magnetic tape cartridge (Sun QIC-24) in UNIX tar format. For the UNIX version, alternate distribution media and formats are available upon request. COMPPAP was developed in 1992.

Description: Structural flaws and cracks may grow under fatigue inducing loads and, upon reaching a critical size, cause structural failure to occur. The growth of these flaws and cracks may occur at load levels well below the ultimate load bearing capability of the structure. The Fatigue Crack Growth Computer Program, NASA/FLAGRO, was developed as an aid in predicting the growth of pre-existing flaws and cracks in structural components of space systems. The earlier version of the program, FLAGRO4, was the primary analysis tool used by Rockwell International and the Shuttle subcontractors for fracture control analysis on the Space Shuttle. NASA/FLAGRO is an enhanced version of the program and incorporates state-of-the-art improvements in both fracture mechanics and computer technology. NASA/FLAGRO provides the fracture mechanics analyst with a computerized method of evaluating the "safe crack growth life" capabilities of structural components. NASA/FLAGRO could also be used to evaluate the damage tolerance aspects of a given structural design. The propagation of an existing crack is governed by the stress field in the vicinity of the crack tip. The stress intensity factor is defined in terms of the relationship between the stress field magnitude and the crack size. The propagation of the crack becomes catastrophic when the local stress intensity factor reaches the fracture toughness of the material. NASA/FLAGRO predicts crack growth using a two-dimensional model which predicts growth independently in two directions based on the calculation of stress intensity factors. The analyst can choose to use either a crack growth rate equation or a nonlinear interpolation routine based on tabular data. The growth rate equation is a modified Forman equation which can be converted to a Paris or Walker equation by substituting different values into the exponent. This equation provides accuracy and versatility and can be fit to data using standard least squares methods. Stress-intensity factor numerical values can be computed for making comparisons or checks of solutions. NASA/FLAGRO can check for failure of a part-through crack in the mode of a through crack when net ligament yielding occurs. NASA/FLAGRO has a number of special subroutines and files which provide enhanced capabilities and easy entry of data. These include crack case solutions, cyclic load spectrums, nondestructive examination initial flaw sizes, table interpolation, and material properties. The materials properties files are divided into two types, a user defined file and a fixed file. Data is entered and stored in the user defined file during program execution, while the fixed file contains already coded-in property value data for many different materials. Prompted input from CRT terminals consists of initial crack definition (which can be defined automatically), rate solution type, flaw type and geometry, material properties (if they are not in the built-in tables of material data), load spectrum data (if not included in the loads spectrum file), and design limit stress levels. NASA/FLAGRO output includes an echo of the input with any error or warning messages, the final crack size, whether or not critical crack size has been reached for the specified stress level, and a life history profile of the crack propagation. NASA/FLAGRO is modularly designed to facilitate revisions and operation on minicomputers. The program was implemented on a DEC VAX 11/780 with the VMS operating system. NASA/FLAGRO is written in FORTRAN77 and has a memory requirement of 1.4 MB. The program was developed in 1986.

Description: Delaminations near the outer surface of a laminate are susceptible to local buckling and buckling-induced delamination propagation when the laminate is subjected to transverse impact loading. This results in a loss of stiffness and strength. TRBUCKL is an unique dynamic delamination buckling and delamination propagation analysis capability that can be incorporated into the structural analysis program, NASTRAN. This capability will aid engineers in the design of structures incorporating composite laminates. The capability consists of: (1) a modification of the direct time integration solution sequence which provides a new analysis algorithm that can be used to predict delamination buckling in a laminate subjected to dynamic loading; and (2) a new method of modeling the composite laminate using plate bending elements and multipoint constraints. The capability now exists to predict the time at which the onset of dynamic delamination buckling occurs, the dynamic buckling mode shape, and the dynamic delamination strain energy release rate. A procedure file for NASTRAN, TRBUCKL predicts both impact induced buckling in composite laminates with initial delaminations and the strain energy release rate due to extension of the delamination. In addition, the file is useful in calculating the dynamic delamination strain energy release rate for a composite laminate under impact loading. This procedure simplifies the simulation of progressive crack extension. TRBUCKL has been incorporated into COSMIC NASTRAN. TRBUCKL is a DMAP Alter for NASTRAN. It is intended for use only with the COSMIC NASTRAN Direct Transient Analysis (RF 9) solution sequence. The program is available as a listing only. TRBUCKL was developed in 1987.

Description: Predictions of fatigue crack growth behavior can be made with the Fatigue Crack Growth Structural Analysis (FASTRAN II) computer program. As cyclic loads are applied to a selected crack configuration with an initial crack size, FASTRAN II predicts crack growth as a function of cyclic load history until either a desired crack size is reached or failure occurs. FASTRAN II is based on plasticity-induced crack-closure behavior of cracks in metallic materials and accounts for load-interaction effects, such as retardation and acceleration, under variable-amplitude loading. The closure model is based on the Dugdale model with modifications to allow plastically deformed material to be left along the crack surfaces as the crack grows. Plane stress and plane strain conditions, as well as conditions between these two, can be simulated in FASTRAN II by using a constraint factor on tensile yielding at the crack front to approximately account for three-dimensional stress states. FASTRAN II contains seventeen predefined crack configurations (standard laboratory fatigue crack growth rate specimens and many common crack configurations found in structures); and the user can define one additional crack configuration. The baseline crack growth rate properties (effective stress-intensity factor against crack growth rate) may be given in either equation or tabular form. For three-dimensional crack configurations, such as surface cracks or corner cracks at holes or notches, the fatigue crack growth rate properties may be different in the crack depth and crack length directions. Final failure of the cracked structure can be modelled with fracture toughness properties using either linear-elastic fracture mechanics (brittle materials), a two-parameter fracture criterion (brittle to ductile materials), or plastic collapse (extremely ductile materials). The crack configurations in FASTRAN II can be subjected to either constant-amplitude, variable-amplitude or spectrum loading. The applied loads may be either tensile or compressive. Several standardized aircraft flight-load histories, such as TWIST, Mini-TWIST, FALSTAFF, Inverted FALSTAFF, Felix and Gaussian, are included as options. FASTRAN II also includes two other methods that will help the user input spectrum load histories. The two methods are: (1) a list of stress points, and (2) a flight-by-flight history of stress points. Examples are provided in the user manual. Developed as a research program, FASTRAN II has successfully predicted crack growth in many metallic materials under various aircraft spectrum loading. A computer program DKEFF which is a part of the FASTRAN II package was also developed to analyze crack growth rate data from laboratory specimens to obtain the effective stress-intensity factor against crack growth rate relations used in FASTRAN II. FASTRAN II is written in standard FORTRAN 77. It has been successfully compiled and implemented on Sun4 series computers running SunOS and on IBM PC compatibles running MS-DOS using the Lahey F77L FORTRAN compiler. Sample input and output data are included with the FASTRAN II package. The UNIX version requires 660K of RAM for execution. The standard distribution medium for the UNIX version (LAR-14865) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. It is also available on a 3.5 inch diskette in UNIX tar format. The standard distribution medium for the MS-DOS version (LAR-14944) is a 5.25 inch 360K MS-DOS format diskette. The contents of the diskette are compressed using the PKWARE archiving tools. The utility to unarchive the files, PKUNZIP.EXE, is included. The program was developed in 1984 and revised in 1992. Sun4 and SunOS are trademarks of Sun Microsystems, Inc. IBM PC is a trademark of International Business Machines Corp. MS-DOS is a trademark of Microsoft, Inc. F77L is a trademark of the Lahey Computer Systems, Inc. UNIX is a registered trademark of AT&T Bell Laboratories. PKWARE and PKUNZIP are trademarks of PKWare, Inc.

Description: The Panel Analysis and Sizing Code (PASCO) was developed for the buckling and vibration analysis and sizing of prismatic structures having an arbitrary cross section. PASCO is primarily intended for analyzing and sizing stiffened panels made of laminated orthotropic materials and is of particular value in analyzing and sizing filamentary composite structures. When used in the analysis mode, PASCO calculates laminate stiffnesses, lamina stress and strains (including the effects of temperature and panel bending), buckling loads, vibration frequencies, and overall panel stiffness. When used in the sizing mode, PASCO adjusts sizing variables to provide a low-mass panel design that carries a set of specified loadings without exceeding buckling or material strength allowables and that meets other design requirements such as upper and lower bounds on sizing variables, upper and lower bounds on overall bending, extensional and shear stiffnesses, and lower bounds on vibration frequencies. Although emphasis in PASCO is placed on flat panels having several identical bays, the only restriction on configuration modeling is that the structure is assumed to be prismatic. In addition, it is assumed that loads and temperatures do not vary along the length of a panel. Because of their wide application in aerospace structures, stiffened panels are readily handled by PASCO. The panel cross section may be composed of an arbitrary assemblage of thin, flat, rectangular plate elements that are connected together along their longitudinal edges. Each plate element consists of a balanced symmetric laminate of any number of layers of orthotropic material. Any group of element widths, layer thicknesses, and layer orientation angles may be selected as sizing variables. Substructuring is available to increase the efficiency of the analysis and to simplify the modeling of complicated structures. The Macintosh version of PASCO includes an interactive, graphic preprocessor called MacPASCO. The main objective of MacPASCO is to make the use of PASCO faster, simpler, and less error-prone. By using a graphical user interface (GUI), MacPASCO simplifies the specification of panel geometry and reduces user input errors, thus making the modeling and analysis of panel designs more efficient. The user draws the initial structural geometry on the computer screen, then uses a combination of graphic and text inputs to: refine the structural geometry, specify information required for analysis such as panel load conditions, and define design variables and constraints for minimum-mass optimization. Composite panel design is an ideal application because the graphical user interface can: serve as a visual aid, eliminate the tedious aspects of text-based input, and eliminate many sources of input errors. The current version of MacPASCO does not implement all the modeling features of PASCO, but has been found to be sufficient for many users. Many difficulties common to text-based inputs are avoided because MacPASCO uses a GUI. First, the graphic displays eliminate syntax errors, like misplaced commas and incorrect command names, because there is no text-based syntax. Second, graphic displays allow the user to see the geometry as it is created and immediately detect and correct any errors. Third, MacPASCO's drawing tools have been designed to avoid modeling errors. Fourth, the graphic displays make revisions to existing structural designs much easier and less error-prone by eliminating the need for the user to conceptualize the text input as geometry. The user can work directly with the geometry displayed on the screen. Finally, MacPASCO automatically generates the correct PASCO input from the geometry displayed on the screen. This input file can be used with any machine version of PASCO to actually perform the analysis and sizing and to output results. The DEC VAX version of PASCO is written in FORTRAN IV for batch execution and has been implemented on a DEC VAX series computer. The Macintosh version of PASCO was developed for Macintosh II series computers with at least 2Mb of RAM running MPW Pascal 3.0 and Language Systems FORTRAN 2.0 under the MPW programming environment. It includes MPW compatible makefiles for compiling the source code. The Macintosh version uses input files compatible with versions of PASCO running on different platforms. MacPASCO is written in Macintosh Programmers Workbench 3.0, MPW Pascal 3.0, and MacAPP 2.0. The Pascal source code is included on the distribution diskette. MacAPP is a development library which is not included. MacPASCO requires a Mac Plus, SE/30, or MacII, IIx, IIcx, IIci, or IIfx running System 6.0 or greater. MacPASCO is System 7.0 compatible. A minimum of 2Mb of RAM is required for execution. The Macintosh version of PASCO is distributed on four 3.5 inch 800K Macintosh format diskettes. The DEC VAX version is distributed on a 9-track 1600 BPI magnetic tape. The PASCO program was developed in 1981, adapted to the DEC VAX in 1983 and to the Macintosh in 1991. MacPASCO was released in 1992.

Description: PLAN2D is a FORTRAN computer program for the plastic analysis of planar rigid frame structures. Given a structure and loading pattern as input, PLAN2D calculates the ultimate load that the structure can sustain before collapse. Element moments and plastic hinge rotations are calculated for the ultimate load. The location of hinges required for a collapse mechanism to form are also determined. The program proceeds in an iterative series of linear elastic analyses. After each iteration the resulting elastic moments in each member are compared to the reserve plastic moment capacity of that member. The member or members that have moments closest to their reserve capacity will determine the minimum load factor and the site where the next hinge is to be inserted. Next, hinges are inserted and the structural stiffness matrix is reformulated. This cycle is repeated until the structure becomes unstable. At this point the ultimate collapse load is calculated by accumulating the minimum load factor from each previous iteration and multiplying them by the original input loads. PLAN2D is based on the program STAN, originally written by Dr. E.L. Wilson at U.C. Berkeley. PLAN2D has several limitations: 1) Although PLAN2D will detect unloading of hinges it does not contain the capability to remove hinges; 2) PLAN2D does not allow the user to input different positive and negative moment capacities and 3) PLAN2D does not consider the interaction between axial and plastic moment capacity. Axial yielding and buckling is ignored as is the reduction in moment capacity due to axial load. PLAN2D is written in FORTRAN and is machine independent. It has been tested on an IBM PC and a DEC MicroVAX. The program was developed in 1988.

Description: Predictions of fatigue crack growth behavior can be made with the Fatigue Crack Growth Structural Analysis (FASTRAN II) computer program. As cyclic loads are applied to a selected crack configuration with an initial crack size, FASTRAN II predicts crack growth as a function of cyclic load history until either a desired crack size is reached or failure occurs. FASTRAN II is based on plasticity-induced crack-closure behavior of cracks in metallic materials and accounts for load-interaction effects, such as retardation and acceleration, under variable-amplitude loading. The closure model is based on the Dugdale model with modifications to allow plastically deformed material to be left along the crack surfaces as the crack grows. Plane stress and plane strain conditions, as well as conditions between these two, can be simulated in FASTRAN II by using a constraint factor on tensile yielding at the crack front to approximately account for three-dimensional stress states. FASTRAN II contains seventeen predefined crack configurations (standard laboratory fatigue crack growth rate specimens and many common crack configurations found in structures); and the user can define one additional crack configuration. The baseline crack growth rate properties (effective stress-intensity factor against crack growth rate) may be given in either equation or tabular form. For three-dimensional crack configurations, such as surface cracks or corner cracks at holes or notches, the fatigue crack growth rate properties may be different in the crack depth and crack length directions. Final failure of the cracked structure can be modelled with fracture toughness properties using either linear-elastic fracture mechanics (brittle materials), a two-parameter fracture criterion (brittle to ductile materials), or plastic collapse (extremely ductile materials). The crack configurations in FASTRAN II can be subjected to either constant-amplitude, variable-amplitude or spectrum loading. The applied loads may be either tensile or compressive. Several standardized aircraft flight-load histories, such as TWIST, Mini-TWIST, FALSTAFF, Inverted FALSTAFF, Felix and Gaussian, are included as options. FASTRAN II also includes two other methods that will help the user input spectrum load histories. The two methods are: (1) a list of stress points, and (2) a flight-by-flight history of stress points. Examples are provided in the user manual. Developed as a research program, FASTRAN II has successfully predicted crack growth in many metallic materials under various aircraft spectrum loading. A computer program DKEFF which is a part of the FASTRAN II package was also developed to analyze crack growth rate data from laboratory specimens to obtain the effective stress-intensity factor against crack growth rate relations used in FASTRAN II. FASTRAN II is written in standard FORTRAN 77. It has been successfully compiled and implemented on Sun4 series computers running SunOS and on IBM PC compatibles running MS-DOS using the Lahey F77L FORTRAN compiler. Sample input and output data are included with the FASTRAN II package. The UNIX version requires 660K of RAM for execution. The standard distribution medium for the UNIX version (LAR-14865) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. It is also available on a 3.5 inch diskette in UNIX tar format. The standard distribution medium for the MS-DOS version (LAR-14944) is a 5.25 inch 360K MS-DOS format diskette. The contents of the diskette are compressed using the PKWARE archiving tools. The utility to unarchive the files, PKUNZIP.EXE, is included. The program was developed in 1984 and revised in 1992. Sun4 and SunOS are trademarks of Sun Microsystems, Inc. IBM PC is a trademark of International Business Machines Corp. MS-DOS is a trademark of Microsoft, Inc. F77L is a trademark of the Lahey Computer Systems, Inc. UNIX is a registered trademark of AT&T Bell Laboratories. PKWARE and PKUNZIP are trademarks of PKWare, Inc.

Description: A three-bay, space, cantilever truss is probabilistically evaluated for adaptive/smart/intelligent behavior. For each behavior, the scatter (ranges) in buckling loads, vibration frequencies, and member axial forces are probabilistically determined. Sensitivities associated with uncertainties in the structure, material and load variables that describe the truss are determined for different probabilities. The relative magnitude for these sensitivities are used to identify significant truss variables that control/classify its behavior to respond as an adaptive/smart/intelligent structure. Results show that the probabilistic buckling loads and vibration frequencies increase for each truss classification, with a substantial increase for intelligent trusses. Similarly, the probabilistic member axial forces reduce for adaptive and intelligent trusses and increase for smart trusses.

Description: The effects of different fabrication procedures to increase the damage tolerance of sandwich panels were studied. Baseline panels consisted of a 25.4 mm premolded core, surfaced with 177 C cure film adhesive and carbon-bismaleimide prepreg which were subsequently cocured onto the core. It was found that panels with a prefabricated skin, which was subsequently bonded onto the core with room temperature cure adhesive, showed greatly increased damage tolerance.

Description: The SNIP program is a FORTRAN computer code that generates NASTRAN structural model thermal loads when given SINDA (or similar thermal model) temperature results. SNIP correlates thermal nodes to structural elements to interface SINDA finite difference thermal models with NASTRAN finite element structural models. Node-to-element correlation includes determining which SINDA nodes should be related to each NASTRAN element and calculating a weighing factor for temperatures associated with each element-related thermal node. SNIP provides structural model thermal loads that accurately reflect thermal model results while reducing the time required to interface thermal and structural models as compared to other methods. SNIP uses thermal model geometry to search the three-dimensional space around each structural element for the nearest thermal nodes. Thermal model geometry is the combination of standard thermal model temperature results from SINDA and structural model geometry from NASTRAN. Thermal and structural models must both be defined in the same, single Cartesian coordinate system. The thermal nodes located nearest each element are used to determine element temperature for thermal distortion and stress analysis. The program shapes the three-dimensional search region while the user controls the size. With these region specifications, the numerical coding of thermal nodes, and the structural element numbers; the code can provide for the separation of substructures during correlation. The input to SNIP contains a file of thermal model temperature results and a physical location of each thermal node in three-dimensional space, combined in a SNIP-unique format. The input also contains a standard NASTRAN input deck for a model made up of plate, shell, beam, and bar elements. SNIP supports the CTRIA, CQUAD, CBAR, and CBEAM elements of NASTRAN. The user adjusts the input parameters in the source code which control the node-to-element correlation. The program outputs NASTRAN element temperature load cards for each element and NASTRAN case control cards for each temperature load set. SNIP also outputs a list of elements that contains the numbers of the SINDA nodes related to each NASTRAN element and the weight that is given to each node in temperature calculations. SNIP is written in ANSI standard FORTRAN 77. The PC version requires a PC FORTRAN compiler and has compiled successfully using Lahey FORTRAN v. 3.0. A core memory of 300k is recommended. The program was developed in 1987.

Description: The following constitutes a summary of this paper: on-orbit identification methodology starts with nonparametric techniques for a priori system identification; development of the nonparametric identification and model determination experiment software has been completed; the validation experiments to be performed on the JPL Control and Identification Technology Validation Laboratory have been designed.

Description: NPLOT is an interactive computer graphics program for plotting undeformed and deformed NASTRAN finite element models (FEMs). Although there are many commercial codes already available for plotting FEMs, these have limited use due to their cost, speed, and lack of features to view BAR elements. NPLOT was specifically developed to overcome these limitations. On a vector type graphics device the two best ways to show depth are by hidden line plotting or haloed line plotting. A hidden line algorithm generates views of models with all hidden lines removed, and a haloed line algorithm displays views with aft lines broken in order to show depth while keeping the entire model visible. A haloed line algorithm is especially useful for plotting models composed of many line elements and few surface elements. The most important feature of NPLOT is its ability to create both hidden line and haloed line views accurately and much more quickly than with any other existing hidden or haloed line algorithms. NPLOT is also capable of plotting a normal wire frame view to display all lines of a model. NPLOT is able to aid in viewing all elements, but it has special features not generally available for plotting BAR elements. These features include plotting of TRUE LENGTH and NORMALIZED offset vectors and orientation vectors. Standard display operations such as rotation and perspective are possible, but different view planes such as X-Y, Y-Z, and X-Z may also be selected. Another display option is the Z-axis cut which allows a portion of the fore part of the model to be cut away to reveal details of the inside of the model. A zoom function is available to terminals with a locator (graphics cursor, joystick, etc.). The user interface of NPLOT is designed to make the program quick and easy to use. A combination of menus and commands with help menus for detailed information about each command allows experienced users greater speed and efficiency. Once a plot is on the screen the interface becomes command driven, enabling the user to manipulate the display or execute a command without having to return to the menu. NPLOT is also able to plot deformed shapes allowing it to perform post-processing. The program can read displacements, either static displacements or eigenvectors, from a MSC/NASTRAN F06 file or a UAI/NASTRAN PRT file. The displacements are written into a unformatted scratch file where they are available for rapid access when the user wishes to display a deformed shape. All subcases or mode shapes can be read in at once. Then it is easy to enable the deformed shape, to change subcases or mode shapes and to change the scale factor for subsequent plots. NPLOT is written in VAX FORTRAN for DEC VAX series computers running VMS. As distributed, the NPLOT source code makes calls to the DI3000 graphics package from Precision Visuals; however, a set of interface routines is provided to translate the DI3000 calls into Tektronix PLOT10/TCS graphics library calls so that NPLOT can use the standard Tektronix 4010 which many PC terminal emulation software programs support. NPLOT is available in VAX BACKUP format on a 9-track 1600 BPI DEC VAX BACKUP format magnetic tape (standard media) or a TK50 tape cartridge. This program was developed in 1991. DEC, VAX, VMS, and TK50 are trademarks of Digital Equipment Corporation. Tektronix, PLOT10, and TCS are trademarks of Tektronix, Inc. DI3000 is a registered trademark of Precision Visuals, Inc. NASTRAN is a registered trademark of the National Aeronautics and Space Administration. MSC/ is a trademark of MacNeal-Schwendler Corporation. UAI is a trademark of Universal Analytics, Inc.

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

Description: The rotational stiffness of hinge joints, and the gap of the joints applied in large deployable trusses, have been experimentally shown to have a significant role in determining such structures' dynamic behavior; an analytical validation of these results is presented for the case where linear rotation springs are used to model the hinge joints employed in a simple beam in trusses. The results obtained indicate that the natural frequencies of these structures depend not only on joint stiffness but also on joint location. Such gap parameters as gap size, stiffness, position, and excitation-force levels, are discussed with a view to a deeper understanding of their effects on a space interferometry system's simulated dynamic responses.

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

Description: Microfracture (fiber/matrix fracture, interphase debonding and interply delamination) in high temperature metal matrix composites (HTMMC), subjected to thermal loading, is computationally simulated. Both unidirectional and crossply SiC/Ti15 composites are evaluated for microfracture driven by thermal loads, using multicell finite element models. Results indicate that under thermal loads alone, microfracture propagation is not as sensitive as it is under mechanical loads.

Description: This paper presents results from a finite element micromechanics analysis of thermally induced stresses in composites at cryogenic temperatures typical of spacecraft operating environments. The influence of microstructural geometry, constituent and interphase properties, and laminate orientation were investigated. Stress field results indicated that significant matrix stresses occur in composites exposed to typical spacecraft thermal excursions; these stresses varied with laminate orientation and circumferential position around the fiber. The major difference in the predicted response of unidirectional and multidirectional laminates was the presence of tensile radial stresses, at the fiber/matrix interface, in multidirectional laminates with off-axis ply angles greater than 15 deg. The predicted damage initiation temperatures and modes were in good agreement with experimental data for both low (207 GPa) and high (517 GPa) modulus carbon fiber/epoxy composites.

Description: The structural integrity of proposed high speed aircraft can be seriously affected by the extremely high surface temperatures and large temperature gradients throughout the vehicle's structure. Variations in the structure's elastic characteristics as a result of thermal effects can be observed by changes in vibration frequency, damping, and mode shape. Analysis codes that predict these changes must be correlated and verified with experimental data. The experimental modal test techniques and procedures used to conduct uniform, nonuniform, and transient thermoelastic vibration tests are presented. Experimental setup and elevated temperature instrumentation considerations are also discussed. Modal data for a 12 by 50 inch aluminum plate heated to a temperature of 475 F are presented. These data show the effect of heat on the plate's modal characteristics. The results indicated that frequency decreased, damping increased, and mode shape remained unchanged as the temperature of the plate was increased.

Description: This paper presents two approaches that will improve the accuracy of explicit approximations of frequency responses. They are important components of structural design optimization and system identification. The first is a product form approximation that improves approximated frequency responses near resonant conditions. This approach is applicable to both design optimization and system identification problems. The second is aimed at improving the computed sensitivity coefficients that are used to solve system identification problems.

Description: LATDYN is a computer code for modeling the Large Angle Transient DYNamics of flexible articulating structures and mechanisms involving joints about which members rotate through large angles. LATDYN extends and brings together some of the aspects of Finite Element Structural Analysis, Multi-Body Dynamics, and Control System Analysis; three disciplines that have been historically separate. It combines significant portions of their distinct capabilities into one single analysis tool. The finite element formulation for flexible bodies in LATDYN extends the conventional finite element formulation by using a convected coordinate system for constructing the equation of motion. LATDYN's formulation allows for large displacements and rotations of finite elements subject to the restriction that deformations within each are small. Also, the finite element approach implemented in LATDYN provides a convergent path for checking solutions simply by increasing mesh density. For rigid bodies and joints LATDYN borrows extensively from methodology used in multi-body dynamics where rigid bodies may be defined and connected together through joints (hinges, ball, universal, sliders, etc.). Joints may be modeled either by constraints or by adding joint degrees of freedom. To eliminate error brought about by the separation of structural analysis and control analysis, LATDYN provides symbolic capabilities for modeling control systems which are integrated with the structural dynamic analysis itself. Its command language contains syntactical structures which perform symbolic operations which are also interfaced directly with the finite element structural model, bypassing the modal approximation. Thus, when the dynamic equations representing the structural model are integrated, the equations representing the control system are integrated along with them as a coupled system. This procedure also has the side benefit of enabling a dramatic simplification of the user interface for modeling control systems. Three FORTRAN computer programs, the LATDYN Program, the Preprocessor, and the Postprocessor, make up the collective LATDYN System. The Preprocessor translates user commands into a form which can be used while the LATDYN program provides the computational core. The Postprocessor allows the user to interactively plot and manage a database of LATDYN transient analysis results. It also includes special facilities for modeling control systems and for programming changes to the model which take place during analysis sequence. The documentation includes a Demonstration Problem Manual for the evaluation and verification of results and a Postprocessor guide. Because the program should be viewed as a byproduct of research on technology development, LATDYN's scope is limited. It does not have a wide library of finite elements, and 3-D Graphics are not available. Nevertheless, it does have a measure of "user friendliness". The LATDYN program was developed over a period of several years and was implemented on a CDC NOS/VE & Convex Unix computer. It is written in FORTRAN 77 and has a virtual memory requirement of 1.46 MB. The program was validated on a DEC MICROVAX operating under VMS 5.2.

Description: An integrated approach to dynamic testing and mathematical model analysis is described. The overall approach addresses four key tasks, namely, pretest planning and analysis, test data acquisition, data reduction and analysis, and test/analysis correlation and mathematical model updates. Several key software programs are employed to accomplish this task. They are a leading finite element code, a sophisticated data analysis processor and a graphical pre- and post-processor along with an advanced interface utility. Several practical structures are used to illustrate tools and concepts employed in the integrated test analysis process.

Description: The paper covers two distinct parts: theory and application. The goal of this work was the reduction of model size with an increase in eigenvalue/vector accuracy. This method is ideal for the condensation of large truss- or beam-type structures. The theoretical approach involves the conversion of a continuum transfer matrix beam element into an 'Exact' dynamic stiffness element. This formulation is implemented in a finite element environment. This results in the need to solve a transcendental eigenvalue problem. Once the eigenvalue is determined the eigenvectors can be reconstructed with any desired spatial precision. No discretization limitations are imposed on the reconstruction. The results of such a combined finite element and transfer matrix formulation is a much smaller FEM eigenvalue problem. This formulation has the ability to extract higher eigenvalues as easily and as accurately as lower eigenvalues. Moreover, one can extract many more eigenvalues/vectors from the model than the number of degrees of freedom in the FEM formulation. Typically, the number of eigenvalues accurately extractable via the 'Exact' element method are at least 8 times the number of degrees of freedom. In contrast, the FEM usually extracts one accurate (within 5 percent) eigenvalue for each 3-4 degrees of freedom. The 'Exact' element results in a 20-30 improvement in the number of accurately extractable eigenvalues and eigenvectors.

Description: The increasing number of applications of fiber-reinforced composites in industry demands a detailed understanding of their material properties and behavior. A three-dimensional finite-element computer program called PAFAC (Plastic and Failure Analysis of Composites) has been developed for the elastic-plastic analysis of fiber-reinforced composite materials and structures. The evaluation of stresses and deformations at edges, cut-outs, and joints is essential in understanding the strength and failure for metal-matrix composites since the onset of plastic yielding starts very early in the loading process as compared to the composite's ultimate strength. Such comprehensive analysis can only be achieved by a finite-element program like PAFAC. PAFAC is particularly suited for the analysis of laminated metal-matrix composites. It can model the elastic-plastic behavior of the matrix phase while the fibers remain elastic. Since the PAFAC program uses a three-dimensional element, the program can also model the individual layers of the laminate to account for thickness effects. In PAFAC, the composite is modeled as a continuum reinforced by cylindrical fibers of vanishingly small diameter which occupy a finite volume fraction of the composite. In this way, the essential axial constraint of the phases is retained. Furthermore, the local stress and strain fields are uniform. The PAFAC finite-element solution is obtained using the displacement method. Solution of the nonlinear equilibrium equations is obtained with a Newton-Raphson iteration technique. The elastic-plastic behavior of composites consisting of aligned, continuous elastic filaments and an elastic-plastic matrix is described in terms of the constituent properties, their volume fractions, and mutual constraints between phases indicated by the geometry of the microstructure. The program uses an iterative procedure to determine the overall response of the laminate, then from the overall response determines the stress state in each phase of the composite material. Failure of the fibers or matrix within an element can also be modeled by PAFAC. PAFAC is written in FORTRAN IV for batch execution and has been implemented on a CDC CYBER 170 series computer with a segmented memory requirement of approximately 66K (octal) of 60 bit words. PAFAC was developed in 1982.

Description: The Panel Analysis and Sizing Code (PASCO) was developed for the buckling and vibration analysis and sizing of prismatic structures having an arbitrary cross section. PASCO is primarily intended for analyzing and sizing stiffened panels made of laminated orthotropic materials and is of particular value in analyzing and sizing filamentary composite structures. When used in the analysis mode, PASCO calculates laminate stiffnesses, lamina stress and strains (including the effects of temperature and panel bending), buckling loads, vibration frequencies, and overall panel stiffness. When used in the sizing mode, PASCO adjusts sizing variables to provide a low-mass panel design that carries a set of specified loadings without exceeding buckling or material strength allowables and that meets other design requirements such as upper and lower bounds on sizing variables, upper and lower bounds on overall bending, extensional and shear stiffnesses, and lower bounds on vibration frequencies. Although emphasis in PASCO is placed on flat panels having several identical bays, the only restriction on configuration modeling is that the structure is assumed to be prismatic. In addition, it is assumed that loads and temperatures do not vary along the length of a panel. Because of their wide application in aerospace structures, stiffened panels are readily handled by PASCO. The panel cross section may be composed of an arbitrary assemblage of thin, flat, rectangular plate elements that are connected together along their longitudinal edges. Each plate element consists of a balanced symmetric laminate of any number of layers of orthotropic material. Any group of element widths, layer thicknesses, and layer orientation angles may be selected as sizing variables. Substructuring is available to increase the efficiency of the analysis and to simplify the modeling of complicated structures. The Macintosh version of PASCO includes an interactive, graphic preprocessor called MacPASCO. The main objective of MacPASCO is to make the use of PASCO faster, simpler, and less error-prone. By using a graphical user interface (GUI), MacPASCO simplifies the specification of panel geometry and reduces user input errors, thus making the modeling and analysis of panel designs more efficient. The user draws the initial structural geometry on the computer screen, then uses a combination of graphic and text inputs to: refine the structural geometry, specify information required for analysis such as panel load conditions, and define design variables and constraints for minimum-mass optimization. Composite panel design is an ideal application because the graphical user interface can: serve as a visual aid, eliminate the tedious aspects of text-based input, and eliminate many sources of input errors. The current version of MacPASCO does not implement all the modeling features of PASCO, but has been found to be sufficient for many users. Many difficulties common to text-based inputs are avoided because MacPASCO uses a GUI. First, the graphic displays eliminate syntax errors, like misplaced commas and incorrect command names, because there is no text-based syntax. Second, graphic displays allow the user to see the geometry as it is created and immediately detect and correct any errors. Third, MacPASCO's drawing tools have been designed to avoid modeling errors. Fourth, the graphic displays make revisions to existing structural designs much easier and less error-prone by eliminating the need for the user to conceptualize the text input as geometry. The user can work directly with the geometry displayed on the screen. Finally, MacPASCO automatically generates the correct PASCO input from the geometry displayed on the screen. This input file can be used with any machine version of PASCO to actually perform the analysis and sizing and to output results. The DEC VAX version of PASCO is written in FORTRAN IV for batch execution and has been implemented on a DEC VAX series computer. The Macintosh version of PASCO was developed for Macintosh II series computers with at least 2Mb of RAM running MPW Pascal 3.0 and Language Systems FORTRAN 2.0 under the MPW programming environment. It includes MPW compatible makefiles for compiling the source code. The Macintosh version uses input files compatible with versions of PASCO running on different platforms. MacPASCO is written in Macintosh Programmers Workbench 3.0, MPW Pascal 3.0, and MacAPP 2.0. The Pascal source code is included on the distribution diskette. MacAPP is a development library which is not included. MacPASCO requires a Mac Plus, SE/30, or MacII, IIx, IIcx, IIci, or IIfx running System 6.0 or greater. MacPASCO is System 7.0 compatible. A minimum of 2Mb of RAM is required for execution. The Macintosh version of PASCO is distributed on four 3.5 inch 800K Macintosh format diskettes. The DEC VAX version is distributed on a 9-track 1600 BPI magnetic tape. The PASCO program was developed in 1981, adapted to the DEC VAX in 1983 and to the Macintosh in 1991. MacPASCO was released in 1992.

Description: A solution to the correlation between structural dynamic test results and finite element analyses of the same components is presented in this paper. Basically, the method can be categorized as a Levenberg-Marquardt type Gauss-Newton method which requires only the differences between FE modal analyses and test results and their first derivatives with respect to preassigned design variables. With proper variable normalization and equation scaling, the method has been made numerically better-conditioned and the inclusion of the Levenberg-Marquardt technique overcomes any remaining difficulty encountered in inverting singular or near-singular matrices. An important feature is that each iteration requires only one function evaluation along with the associated design sensitivity analysis and so the procedure is computationally efficient.

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

Description: The strain gage blind hole-drilling technique may be used to determine residual stresses at and below the surface of components. In this paper, the hole-drilling analysis methodology for thick plates is reviewed, and experimental data are used to evaluate the methodology and to assess its applicability to thin plates. Data on the effects of gage pattern, surface preparation, hole spacing, hole eccentricity, and stress level are also presented.

Description: The energy absorption response and crushing characteristics of geometrically scaled graphite-Kevlar epoxy composite plates were investigated. Three different trigger mechanisms including chamfer, notch, and steeple geometries were incorporated into the plate specimens to initiate crushing. Sustained crushing was achieved with a simple test fixture which provided lateral support to prevent global buckling. Values of specific sustained crushing stress (SSCS) were obtained which were comparable to values reported for tube specimens from previously published data. Two sizes of hybrid plates were fabricated; a baseline or model plate, and a full-scale plate with in-plane dimensions scaled by a factor of two. The thickness dimension of the full-scale plates was increased using two different techniques; the ply-level method in which each ply orientation in the baseline laminate stacking sequence is doubled, and the sublaminate technique in which the baseline laminate stacking sequence is repeated as a group. Results indicated that the SSCS is independent of trigger mechanism geometry. However, a reduction in the SSCS of 10-25 percent was observed for the full-scale plates as compared with the baseline specimens, indicating a scaling effect in the crushing response.

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

Description: The natural vibration frequencies of a structure can be affected by inducing stress in the structure. The success of this kind of control of the resonant frequencies of a truss structure depends on the geometry of the structure. It is shown that in adaptive truss structures the method is effective for vibrations in less stiff directions, such as the normal direction of the plane containing all of the bars of a node, suggesting its applicability for cable, membrane, and thin plate and shell structures.

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

Description: An experimental technique was used to measure structural intensity through an aircraft fuselage with an excitation load applied near one of the wing attachment locations. The fuselage was relatively large, requiring several measurement locations to analyze the intensity flow through the whole of the structure. For the measurement of structural intensity, the use of a transducer array was necessary at every location of interest. A trade-off was therefore required between the number of measurement transducers, the mounting of these transducers, and the accuracy of the measurements. Using four accelerometers mounted on a bakelite platform, structural intensity vectors were measured at locations distributed throughout the fuselage. The results of these measurements, together with a discussion on the suitability of the approach for measuring structural intensity on a real structure, are presented in this paper.

Description: The feasibility of using scale model testing for predicting the full-scale behavior of flat composite coupons loaded in tension and beam-columns loaded in flexure is examined. Classical laws of similitude are applied to fabricate and test replica model specimens to identify scaling effects in the load response, strength, and mode of failure. Experiments were performed on graphite-epoxy composite specimens having different laminate stacking sequences and a range of scaled sizes. From the experiments it was deduced that the elastic response of scaled composite specimens was independent of size. However, a significant scale effect in strength was observed. In addition, a transition in failure mode was observed among scaled specimens of certain laminate stacking sequences. A Weibull statistical model and a fracture mechanics based model were applied to predict the strength scale effect since standard failure criteria cannot account for the influence of absolute specimen size on strength.

Description: BUCKO is a computer program developed to predict the buckling load of a rectangular compression-loaded orthotropic plate with a centrally located cutout. The plate is assumed to be a balanced, symmetric laminate of uniform thickness. The cutout shape can be elliptical, circular, rectangular, or square. The BUCKO package includes sample data that demonstrates the essence of the program and its ease of usage. BUCKO uses an approximate one-dimensional formulation of the classical two-dimensional buckling problem following the Kantorovich method. The boundary conditions are considered to be simply supported unloaded edges and either clamped or simply supported loaded edges. The plate is loaded in uniaxial compression by either uniformly displacing or uniformly stressing two opposite edges of the plate. The BUCKO analysis consists of two parts: calculation of the inplane stress distribution prior to buckling, and calculation of the plate axial load and displacement at buckling. User input includes plate planform and cutout geometry, plate membrane and bending stiffnesses, finite difference parameters, boundary condition data, and loading data. Results generated by BUCKO are the prebuckling strain energy, inplane stress resultants, buckling mode shape, critical end shortening, and average axial and transverse strains at buckling. BUCKO is written in FORTRAN V for batch execution and has been implemented on a CDC CYBER 170 series computer operating under NOS with a central memory requirement of approximately 343K of 60 bit words. This program was developed in 1984 and was last updated in 1990.

Description: A model identification methodology for structural dynamics has been applied to simulated vibrational data as a first step in evaluating its accuracy. The evaluation has taken into account a wide variety of factors which affect the accuracy of the procedure. The effects of each of these factors were observed in both the response time histories and the estimates of the parameters of the model by comparing them with the exact values of the system. Each factor was varied independently but combinations of these have also been considered in an effort to simulate real situations. The results of the tests have shown that for the chain model, the procedure yields robust estimates of the stiffness parameters under the conditions studied whenever uniqueness is ensured. When inaccuracies occur in the results, they are intimately related to non-uniqueness conditions inherent in the inverse problem and not to shortcomings in the methodology.

Description: An elastic-plastic finite-element analysis with a critical crack-tip-opening displacement criterion was used to simulate fracture of various size compact and bend specimens made of HY-130 steel. From the calculated load-crack-extension and load-displacement curves, J-resistance (J-R) curves were determined by several methods. The simulated 3-R curves were insensitive to specimen size up to maximum load but were sensitive to specimen configuration for crack extensions greater than 10 percent of the initial uncracked ligament length.

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

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

Description: The concept of pressure proof testing of fuselage structures with fatigue cracks to insure structural integrity was evaluated from a fracture mechanics viewpoint. A generic analytical and experimental investigation was conducted on uniaxially loaded flat panels with crack configurations and stress levels typical of longitudinal lap splice joints in commercial transport aircraft fuselages. The results revealed that the remaining fatigue life after a proof cycle was longer than that without the proof cycle because of crack growth retardation due to increased crack closure. However, based on a crack length that is slightly less than the critical value at the maximum proof stress, the minimum assured life or proof test interval must be no more than 550 pressure cycles for a 1.33 proof factor and 1530 pressure cycles for a 1.5 proof factor to prevent in-flight failures.

Description: The experiment is aimed at controlling the boundary layer transition location and the plate vibration when excited by a flow and an upstream sound source. Sound has been found to affect the flow at the leading edge and the response of a flexible plate in a boundary layer. Because the sound induces early transition, the panel vibration is acoustically coupled to the turbulent boundary layer by the upstream radiation. Localized surface heating at the leading edge delays the transition location downstream of the flexible plate. The response of the plate excited by a turbulent boundary layer (without sound) shows that the plate is forced to vibrate at different frequencies and with different amplitudes as the flow velocity changes indicating that the plate is driven by the convective waves of the boundary layer. The acoustic disturbances induced by the upstream sound dominate the response of the plate when the boundary layer is either turbulent or laminar. Active vibration control was used to reduce the sound induced displacement amplitude of the plate.

Description: A spectral analysis method is presented for the extraction of frequency dependent strain energy release rates and stress intensity factors. The approach is based on the crack closure technique and is formulated directly in the frequency domain. Because of this, it is computationally more efficient than an alternative time domain approach for systems exhibiting peak behavior. The method is demonstrated for a center cracked panel subject to static in-plane and random acoustic loading.

Description: A finite element approach is developed for beam type pipes undergoing large deflections subjected to random loadings. The influence of fluid velocity on the random response is investigated. The root-mean-square (rms) deflections and frequencies for different sound spectrum level values are determined for pipes with both ends either simply supported or clamped. The required number of modes to achieve accurate rms deflections is studied. The prediction of fatigue life is then based on the maximum rms stress. This analytical investigation will help to broaden the basic understanding of the role of fluid flow within structures subjected to random excitations.

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

Description: The equilibrium equations and the compatibility conditions are fundamental to the analyses of structures. However, anyone who undertakes even a cursory generic study of the compatibility conditions can discover, with little effort, that historically this facet of structural mechanics had not been adequately researched by the profession. Now the compatibility conditions (CC's) have been researched and are understood to a great extent. For finite element discretizations, the CC's are banded and can be divided into three distinct categories: (1) the interface CC's, (2) the cluster or field CC's, and (3) the external CC's. The generation of CC's requires the separating of a local region, then writing the deformation displacement relation (ddr) for the region, and finally, the eliminating of the displacements from the ddr. The procedure to generate all three types of CC's is presented and illustrated through examples of finite element models. The uniqueness of the CC's thus generated is shown. The utilization of CC's has resulted in the novel integrated force method (IFM). The solution that is obtained by the IFM converges with a significantly fewer number of elements, compared to the stiffness and the hybrid methods.

Description: The efficient solution of the eigenvalue problem that results from inserting passive dampers with variable stiffness and damping coefficients into a structure is addressed. Eigenanalysis of reduced models obtained by retaining a number of normal modes augmented with Ritz vectors corresponding to the static solutions resulting from the load patterns introduced by the dampers has been empirically shown to yield excellent approximations to the full eigenvalue problem. An analysis of this technique in the case of a single damper is presented. A priori and a posteriori error estimates are generated and tested on numerical examples. Comparison theorems with modally truncated models and a Markov parameter matching reduced-order model are derived. These theorems corroborate the heuristic that residual flexibility methods improve low-frequency approximation of the system. The analysis leads to other techniques for eigenvalue approximation. Approximate closed-form solutions are derived that include a refinement to eigenvalue derivative methods for approximation. An efficient Newton scheme is also developed. A numerical example is presented demonstrating the effectiveness of each of these methods.

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

Description: Identification of large space structures' distributed mass, stiffness, and energy dissipation characteristics poses formidable analytical, numerical, and implementation difficulties. Development of reliable on-orbit structural identification methods is important for implementing active vibration suppression concepts which are under widespread study in the large space structures community. Near the heart of the identification problem lies the necessity of making a large number of spatially distributed measurements of the structure's vibratory response and the associated force/moment inputs with sufficient spatial and frequency resolution. In the present paper, we discuss a method whereby tens of active or passive (retro-reflecting) targets on the structure are tracked simultaneously by the focal planes of two or more video cameras mounted on an adjacent platform. Triangulation (optical ray intersection) of the conjugate image centroids yield inertial trajectories of each target on the structure. Given the triangulated motion of the targets, we apply and extend methodology developed by Creamer, Junkins, and Juang to identify the frequencies, mode shapes, and updated estimates for the mass/stiffness/damping parameterization of the structure. The methodology is semi-automated, for example, the post experiment analysis of the video imagery to determine the inertial trajectories of the targets typically requires less than thirty minutes of real time. Using methodology discussed herein, the frequency response of a large number of points on the structure (where reflective targets are mounted) on the structure can be determined from optical measurements alone. For comparison purposes, we also utilize measurements from accelerometers and a calibrated impulse hammer. While our experimental work remains in a research stage of development, we have successfully tracked and stereo triangulated 20 targets (on a vibrating cantilevered grid structure) at a sample frequency of 200 HZ, and have established conclusively the feasibility and desirability of this approach. We discuss, in summary, recent advances in analog and digital video processing methodology, actuation methods, and bring them to bear on the structural identification problem. We include a brief discussion of our experimental hardware and some recent experimental results which support the practical feasibility of this structural vibration sensing approach.

Description: Limitations of the frequency domain methods in analyzing structura1 vibrations has created an awareness of the comparative merits of the time domain methods. Although time domain methods would be ideal for modeling large precisions space systems, the popular methods based on fitting theoretical response to actual data by least squares are too sensitive to noise and require too much data to be suitable for orbiting space crafts. This paper briefly reviews the theory and illustrative applications of a time domain methodology called Data Dependent Systems (DDS) that eliminates these limitations. Simulation results are presented to demonstrate a better than 4-place accuracy in the identifications of all system parameters, both modal (frequencies, damping ratios, and mode shapes) and physical (mass, stiffness, and damping matrices).

Description: The objective of the 12-meter truss modal test is to experimentally determine the frequencies, damping, and mode shapes for the first 6 modes in both principle axes and to use this information to update the FEM. These objectives will lead us to our goal of actively controlling the flexible modes of the truss. A secondary objective is to evaluate our capabilities to ground test this class of structures.

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

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

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

Description: Model equations are presented for approximate methods simulating the long-term behavior of composite materials and structures in hot/humid service environments. These equations allow laminate property upgradings with time, and can account for the effects of service environments on creep response. These methodologies are illustrated for various individual and coupled temperature/moisture, longitudinal/transverse, and composite material type cases. Creep deformation is noted to rise dramatically for cases of matrix-borne, but not of fiber-borne, loading in hot, humid environments; the coupled influence of temperature and moisture is greater than a mere combination of their individual influences.

Description: The models employed in the present computational methods for evaluating severe service-environment effects on adhesively bonded joints in composites are based on composite analyses and structural mechanics, encompassing nonlinear environmental degradation. The methods are demonstrated for the case of a butt joint with a single doubler, subjected to the environmental effects as well as static and cyclic loads. The highest joint strength is noted to be required in the case of cyclic loads and hygrothermal service environments; margins of safety for adhesive material stresses decline rapidly in such cases.

Description: Simplified predictive methods and models to computationally simulate durability and damage in polymer matrix composite materials/structures are described. The models include (1) progressive fracture, (2) progressively damaged structural behavior, (3) progressive fracture in aggressive environments, (4) stress concentrations, and (5) impact resistance. Several examples are included to illustrate applications of the models and to identify significant parameters and sensitivities. Comparisons with limited experimental data are made.

Description: The present paper describes the application of the computer code IPACS (Integrated Probabilistic Assessment of Composite Structures) to air craft wing type structures. The code performs a complete probabilistic structural analysis for composites taking into account the uncertainties in geometry, boundary conditions, material properties, laminate lay-ups and loads. Results of the analysis are presented in terms of cumulative distribution functions (CDF) and probability density function (PDF) of life of a wing type composite structure under different hygrothermal environments subjected to random pressure. The sensitivity of fatigue life to a number of critical structural/material variables is also computed from the analysis.

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

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

Description: After briefly describing the principles of frozen stress photoelastic and moire interferometric analyses, and the corresponding algorithms for converting optical data from each method into stress intensity factors (SIF), the methods are applied to the determination of crack shapes, SIF determination, crack closure displacement fields, and pre-crack damage mechanisms in typical aircraft component configurations.

Description: This paper presents a review of some common small-crack test specimens, the underlying causes of the small-crack effect, and the fracture-mechanics parameters that have been used to correlate or predict their growth behavior. This review concentrates on continuum mechanics concepts and on the nonlinear behavior of small cracks. The paper reviews some stress-intensity factor solutions for small-crack test specimens and develops some simple elastic-plastic J integral and cyclic J integral expressions that include the influence of crack-closure. These parameters were applied to small-crack growth data on two aluminum alloys, and a fatigue life prediction methodology is demonstrated. For these materials, the crack-closure transient from the plastic wake was found to be the major factor in causing the small-crack effect.

Description: The bounds on the equivalent elastic material properties of a composite are presently addressed by a unified energy approach which is valid for both unidirectional and 2D and 3D woven composites. The unit cell considered is assumed to consist, first, of the actual composite arrangement of the fibers and matrix material, and then, of an equivalent pseudohomogeneous material. Equating the strain energies due to the two arrangements yields an estimate of the upper bound for the material equivalent properties; successive increases in the order of displacement field that is assumed in the composite arrangement will successively produce improved upper bound estimates.

Description: An inverse technique was used to calculate through-thickness fatigue crack closure behavior. The through-thickness variation in crack opening stress-intensity factor was calculated by considering the variation in the three-dimensional stress-intensity factor, the variation in crack growth rate along the crack front, and a relationship between the crack growth rate and effective stress-intensity factor range (da/dN-Delta-K(eff)). The three-dimensional stress-intensity factor variation was obtained from an elastic finite element analysis of specific crack-front profiles observed experimentally. The variation in crack growth rate along the crack front was obtained experimentally from comparison of observed crack front changes. The da/dN-Delta-K(eff) relationship was estimated from high stress ratio, constant load amplitude, and fatigue crack growth tests. The through-thickness crack opening stress-intensity factor results agreed with crack opening measurements obtained from fatigue striations, near-tip strain gages, and remote strain and displacement gages.

Description: The evolution of high-temperature, creep-fatigue, life-prediction methods used for cyclic crack initiation is traced from inception in the late 1940s. The methods reviewed are material models as opposed to structural life prediction models. Material life models are used by both structural durability analysts and by material scientists. The latter use micromechanistic models as guidance to improve a material's crack initiation resistance. Nearly one hundred approaches and their variations have been proposed to date. This proliferation poses a problem in deciding which method is most appropriate for a given application. Approaches have been identified as being combinations of fourteen different classifications. This review is intended to aid both developers and users of high-temperature fatigue life prediction methods by providing a background from which choices can be made.

Description: State-of-the-art nonlinear finite element analysis techniques are evaluated by applying them to a realistic aircraft structural component. A wing panel from the V-22 tiltrotor aircraft is chosen because it is a typical modern aircraft structural component for which there is experimental data for comparison of results. From blueprints and drawings, a very detailed finite element model containing 2284 9-node Assumed Natural-Coordinate Strain elements was generated. A novel solution strategy which accounts for geometric nonlinearity through the use of corotating element reference frames and nonlinear strain-displacement relations is used to analyze this detailed model. Results from linear analyses using the same finite element model are presented in order to illustrate the advantages and costs of the nonlinear analysis as compared with the more traditional linear analysis.

Description: A thermoviscoplastic finite element method employing the Bodner-Partom constitutve model is used to investigate the response of simplified thermal-structural models to intense local heating. The computational method formulates the problem in rate and advances the solution in time by numerical integration. The thermoviscoplastic response of simplified structures with prescribed temperatures is investigated. With rapid rises of temperature, the nickel alloy structures display initially higher yield stresses due to strain rate effects. As temperatures approach elevated values, yield stress and stiffness degrade rapidly and pronounced plastic deformation occurs.

Description: A three-bay, space, cantilever truss is probabilistically evaluated to describe progressive buckling and truss collapse in view of the numerous uncertainties associated with the structural, material, and load variables that describe the truss. Initially, the truss is deterministically analyzed for member forces, and members in which the axial force exceeds the Euler buckling load are identified. These members are then discretized with several intermediate nodes, and a probabilistic buckling analysis is performed on the truss to obtain its probabilistic buckling loads and the respective mode shapes. Furthermore, sensitivities associated with the uncertainties in the primitive variables are investigated, margin of safety values for the truss are determined, and truss end node displacements are noted. These steps are repeated by sequentially removing buckled members until onset of truss collapse is reached. Results show that this procedure yields an optimum truss configuration for a given loading and for a specified reliability.

Description: In this paper, the post impact compressive behavior of polymeric composites is studied both analytically and experimentally. In the analytical study, a closed-form solution is obtained for postbuckling of composites with a circular delamination. Using this solution, the reduced stiffness of the impact damaged region is calculated and the residual compressive strengths of quasi-isotropic laminates are predicted as a function of damage size. In the experimental study, in-plane displacements near the damaged region are determined by employing a micro moire interferometry technique. The strain concentrations as well as the compression-after-impact strength have been measured and compared with analytical predictions. Good agreement between the predicted results and experimental data was observed.

Description: A study was conducted to identify one of the mechanisms that contributes to the reduced compression strength of composite materials with through-the-thickness (TTT) reinforcements. In this study a series of thick (0/90) laminates with stitched and integrally woven TTT reinforcements were fabricated and statically tested. In both the stitching and weaving process a surface loop of TTT reinforcement yarn is created between successive TTT penetrations. It was shown that the surface loop of the TTT reinforcement 'kinked' the in-plane fibers in such a manner that they were made ineffective in carrying compressive load. The improvement in strength by removal of the surface loop and 'kinked' in-plane fibers was between 7 and 35 percent.

Description: The nonlinear behavior of a high-temperature metal-matrix composite (HT-MMC) was simulated by using the metal matrix composite analyzer (METCAN) computer code. The simulation started with the fabrication process, proceeded to thermomechanical cyclic loading, and ended with the application of a monotonic load. Classical laminate theory and composite micromechanics and macromechanics are used in METCAN, along with a multifactor interaction model for the constituents behavior. The simulation of the stress-strain behavior from the macromechanical and the micromechanical points of view, as well as the initiation and final failure of the constituents and the plies in the composite, were examined in detail. It was shown that, when the fibers and the matrix were perfectly bonded, the fracture started in the matrix and then propagated with increasing load to the fibers. After the fibers fractured, the composite lost its capacity to carry additional load and fractured.

Description: Impacters of various masses were dropped from various heights onto thick graphite/epoxy filament-wound cylinders. The cylinders represented filament-wound cases made for the booster motors of the Space Shuttle. Tups of various shapes were affixed to the impacters. Some of the cylinders were filled with inert propellant, and some were empty. The cylinders were impacted numerous times around the circumference and then cut into tension coupons, each containing an impact site. The size of the damage and the residual tension strength were measured. For hemispherical tups, strength was reduced as much as 30 percent by nonvisible damage. The damage consisted of matrix cracking and broken fibers. Analytical methods were used to predict the damage and residual tension strength. A factor of safety to account for nonvisible damage was determined. For corner and rod shaped tups, any damage that resulted in strength loss was readily visible.

Description: The paper first presents the details of the development of a new six-noded plane triangular finite dynamic element. A block Lanczos algorithm is developed next for the accurate and efficient solution of the quadratic matrix eigenvalue problem associated with the finite dynamic element formulation. The resulting computer program fully exploits matrix sparsity inherent in such a discretization and proves to be most efficient for the extraction of the usually required first few roots and vectors, including repeated ones. Most importantly, the present eigenproblem solution is shown to be comparable to that of the corresponding finite element analysis, thereby rendering the associated dynamic element method rather attractive owing to superior convergence characteristics of such elements, presented herein.

Description: As a means of determining a stress intensity factor solution, the compliance properties of an ARALL-2 laminated-sheet composite were investigated. Fatigue crack growth rate (FCGR) tests were conducted on middle crack tension (MT) specimens fabricated from a layup consisting of three sheets of 2024-T3 aluminum bonded together with unidirectional aramid fibers embedded in epoxy. Excellent fatigue crack growth properties are obtained by the presence of unbroken aramid fibers in the wake of the crack tip. These unbroken fibers act as a bridging mechanism to inhibit further crack growth. To quantify the effect of maximum fatigue load on compliance, a series of FCGR tests were performed. Effective crack lengths were determined to be at least 10 mm shorter than surface measured crack lengths for a 76-mm-wide specimen. The bridging zone was estimated to be at least 5 mm. Compliance and stress intensity factor as functions of effective crack length were determined.

Description: Tapered (0 deg) laminates of S2/CE9000 and S2/SP250 glass/epoxies, and IM6/1827I graphite/epoxy were tested in cyclic tension. The specimens usually showed some initial stable delaminations in the tapered region, but these did not affect the stiffness of the specimens, and loading was continued until the specimens either delaminated unstably, or reached 10(exp 6) to 2 x 10 (exp 7) million cycles with no unstable delamination. The final unstable delamination originated at the junction of the thin and tapered regions. A finite-element model was developed for the tapered laminate with and without the initial stable delaminations observed in the tests. The analysis showed that for both cases the most likely place for an opening (Mode I) delamination to originate is at the junction of the taper and thin regions. For each material type, the models were used to calculate the strain energy release rate, G, associated with delaminations originating at that junction and growing either into the thin region or tapered region. For the materials tested, cyclic G(sub Imax) values from DCB tests were used with the maximum strain energy release rates calculated from the finite-element analysis to predict the onset of unstable delamination at the junction as a function of fatigue cycles. The predictions were compared to experimental values of maximum cyclic load as a function of cycles to unstable delamination from fatigue tests in tapered laminates. For the IM6/1827I and S2/SP250 laminates, the predictions agreed very well with the test data. Predicted values for the S2/CE9000 were conservative compared to the test data.

Description: The present study proposes a mathematical model utilizing the internal state variable concept for predicting the upper bounds of the reduced axial and shear stiffnesses in cross-ply laminates with matrix cracks. The displacement components at the matrix crack surfaces are explicitly expressed in terms of the observable axial and shear strains and the undamaged material properties. The reduced axial and shear stiffnesses are predicted for glass/epoxy and graphite/epoxy laminates. Comparison of the model with other theoretical and experimental studies is also presented to confirm direct applicability of the model to angle-ply laminates with matrix cracks subjected to general in-plane loading.

Description: A simple strength of materials analysis was developed for a double-cantilever beam (DCB) specimen to account for geometric nonlinearity effects due to large deflections and T-tabs. A new DCB data analysis procedure was developed to include the effects of these nonlinearities. The results of the analysis were evaluated by DCB tests performed for materials having a wide range of toughnesses. The materials used in the present study were T300/5208, IM7/8551-7, and AS4/PEEK. Based on the present analysis, for a typical deflection/crack length ratio of 0.3 (for AS4/PEEK), T-tabs and large deflections cause a 15 and 3 percent error, respectively, in the computer Mode I strain energy release rate. Design guidelines for DCB specimen thickness and T-tab height were also developed in order to keep errors due to these nonlinearities within 2 percent. Based on the test results, for both hinged and tabbed specimens, the effects of large deflection on the Mode I fracture toughness (G sub Ic) were almost negligible (less than 1 percent) in the case of T300/5208 and IM7/8551-7; however, AS4/PEEK showed a 2 to 3 percent effect. The effects of T-tabs G sub Ic were more significant for all the materials with T300/5208 showing a 5 percent error, IM7/8551-7 a 15 percent error, and, AS4/PEEK a 20 percent error.

Description: A test rig for testing a thick split cantilever beam for scissoring delamination (mode 3) fracture toughness was developed. A 3-D finite element analysis was conducted on the test specimen to determine the strain energy release rate, G, distribution along the delamination front. The virtual crack closure technique was used to calculate the G components resulting from interlaminar tension, GI, interlaminar sliding shear, GII, and interlaminar tearing shear, GIII. The finite element analysis showed that at the delamination front no GI component existed, but a GII component was present in addition to a GIII component. Furthermore, near the free edges, the GII component was significantly higher than the GIII component. The GII/GIII ratio was found to increase with delamination length but was insensitive to the beam depth. The presence of GII at the delamination front was verified experimentally by examination of the failure surfaces. At the center of the beam, where the failure was in mode 3, there was significant fiber bridging. However, at the edges of the beam where the failure was in mode 3, there was no fiber bridging and mode 2 shear hackles were observed. Therefore, it was concluded that the split cantilever beam configuration does not represent a pure mode 3 test. The experimental work showed that the mode 2 fracture toughness, GIIc, must be less than the mode 3 fracture toughness, GIIIc. Therefore, a conservative approach to characterizing mode 3 delamination is to equate GIIIc to GIIc.

Description: A new solution method has recently been developed to analyze the thermal and structural response of polymeric composites during chemical decomposition (Sullivan and Salamon, 1992a). This method couples the equations of energy, mass transport and momentum and solves these equations simultaneously at each time step. In this paper, the couples solution method is applied to analyze the response of a glass phenolic composite during its chemical decomposition. The numerical solution to this particular problem was originally performed by Henderson, et al. (1985) and the actual response was measured by Henderson and Hagen (1985) and Ramamurthy (1988). In addition to further verification of the coupled solution method, the present solution offers additional insight into certain physical events which occur during chemical decomposition.

Description: Thermal deformations and stresses were studied in a silicon-carbide/aluminum filamentary composite at temperatures up to 370 C (700 F). Longitudinal and transverse thermal strains were measured with strain gages and a dilatometer. An elastoplastic micromechanical analysis based on a one-dimensional rule-of-mixtures model and an axisymmetric two-material composite cylinder model was performed. It was established that beyond a critical temperature thermal strains become nonlinear with decreasing longitudinal and increasing transverse thermal-expansion coefficients. This behavior was attributed to the plastic stresses in the aluminum matrix above the critical temperature. An elastoplastic analysis of both micromechanical models was performed to determine the stress distributions and thermal deformation in the fiber and matrix of the composite. While only axial stresses can be determined by the rule-of-mixtures model, the complete triaxial state of stress is established by the composite cylinder model. Theoretical predictions for the two thermal-expansion coefficients were in satisfactory agreement with experimental results.

Description: The problem of optimal placement of active members which are used for vibration control in adaptive truss structures is investigated. The control scheme is based on the method of eigenvalue assignment as a means of shaping the transient response of the controlled adaptive structures, and the minimization of required control action is considered as the optimization criterion. To this end, a performance index which measures the control strokes of active members is formulated in an efficient way. In order to reduce the computation burden, particularly for the case where the locations of active members have to be selected from a large set of available sites, several heuristic searching schemes are proposed for obtaining the near-optimal locations. The proposed schemes significantly reduce the computational complexity of placing multiple active members to the order of that when a single active member is placed.

Description: The mixed-mode bending (MMB) test uses a lever to apply simultaneously mode I and mode II loading to a split-beam specimen. An iterative analysis that accounts for the geometric nonlinearity of the MMB test was developed. The analysis accurately predicted the measured load-displacement response and the strain energy release rate, G, of an MMB test specimen made of AS4/PEEK. The errors in G when calculated using linear theory were found to be as large as 30 percent in some cases. Because it would be inconvenient to use a nonlinear analysis to analyze MMB data, the MMB apparatus was redesigned to minimize the nonlinearity. With the improved apparatus, loads are applied just above the midplane of the test specimen through a roller attached to the lever. This apparatus was demonstrated by measuring the mixed-mode delamination fracture toughhess of the test specimen. The nonlinearity errors associated with testing this tough composite material were less than +/- 3 percent. The data from the improved MMB apparatus analyzed with a linear analysis were similar to those found with the original apparatus and the nonlinear analysis.

Description: Mixed-mode matrix fracture in central notched off-axis unidirectional composite laminates was investigated. A limited number of unidirectional tensile type specimens with a central, horizontal, notch were tested. Crack initiation and propagation were examined under various local stress fields that were controlled by fiber orientations. The tested specimens were simulated using a two dimensional finite element method with constant strain loading. The strain energy release rates along the crack were evaluated via crack closure technique. The variation of critical strain energy rates with off-axis angle was studied. The results from single (one-sided) and double (two-sided) crack simulations were presented and compared.

Description: Aluminum plates of various length/width ratios loaded in compression are examined analytically in terms of buckling and postbuckling for comparison with classical theory. The plates are considered to be subjected to longitudinal compressive displacement with no edge stress and are assumed to be simply supported at the side edges. The average longitudinal direct stresses are computed to plot characteristic curves by means of the following theories: 3D flexibility, higher-order shear deformation, first-order shear deformation, and classical von Karman. The 3D flexibility approach yields the lowest results with more unknowns than the higher-order shear-deformation method. The 3D flexibility approach is considered to be the most accurate, and it is shown that the calculated resultants and displacements tend to vary when consideration is given to the effects of transverse shear.

Description: An attempt is made to establish the validity of the antiplane solution and to calculate all components of the stress intensity factor (SIF) at the crack front for thin plates subjected to mode III loading. The stress analysis was performed for a range of thickness to crack length (b/a) ratio, which covered long cracks in a thin aircraft fuselage structure (b/a = 0.05).

Description: Some of the basic issues of ratchetting behavior that are being addressed by the viscoplastic modeling community are discussed. Some of the shortcomings of existing viscoplastic models are examined in the light of the difficulty involved in using established viscoplastic modeling techniques to predict ratchetting accurately.